Ripples in a string coupled to Glauber spins

Buckling in a rotationally invariant spin-elastic model

Scanning tunneling microscopy experiments have revealed a spontaneous rippled-to-buckled transition in heated graphene sheets, in absence of any mechanical load. Several models relying on a simplified picture of the interaction between elastic and internal, electronic, degrees of freedom have been proposed to understand this phenomenon. Nevertheless, these models are not fully consistent with the classical theory of elasticity, since they do not preserve rotational invariance. Herein, we develop and analyze an alternative classical spin-elastic model that preserves rotational invariance while giving a qualitative account of the rippled-to-buckled transition. By integrating over the internal degrees of freedom, an effective free energy for the elastic modes is derived, which only depends on the curvature. Minimization of this free energy gives rise to the emergence of different mechanical phases, whose thermodynamic stability is thoroughly analyzed, both analytically and numerically. All phases are characterized by a spatially homogeneous curvature, which plays the role of the order parameter for the rippled-to-buckled transition, in both the one- and two-dimensional cases. In the latter, our focus is put on the honeycomb lattice, which is representative of actual graphene.

Martingale-induced local invariance in progressive quenching

Progressive quenching (PQ) is a stochastic process during which one fixes, one after another, the degrees of freedom of a globally coupled Ising spin system while letting it thermalize through a heat bath. It has previously been shown that during PQ, the mean equilibrium spin value follows a martingale process and this process can characterize the memory of the system. In the present study, we find that the aforementioned martingale implies a local invariance of the path weight for the total quenched magnetization, the Markovian process whose increment is the spin that is fixed last. Consequently, PQ lets the probability distribution for the total quenched magnetization evolve while keeping the Boltzmann-like factor, or a canonical structure, under constraint, which consists of a path-independent potential and a path-counting entropy. Moreover, when the PQ starts from full equilibrium, the probability distribution at each stage of PQ is found to be the limit distribution of what we call recycled quenching, the process in which a randomly chosen quenched spin is unquenched after a single step of PQ. The local invariance is directly derived from the martingale property, and not from other known theorems on martingale processes.

Fluctuation-induced current from freestanding graphene

At room temperature, micron-sized sheets of freestanding graphene are in constant motion, even in the presence of an applied bias voltage. We quantify the out-of-plane movement by collecting the displacement current using a nearby small-area metal electrode and present an Ito-Langevin model for the motion coupled to a circuit containing diodes. Numerical simulations show that the system reaches thermal equilibrium and the average rates of heat and work provided by stochastic thermodynamics tend quickly to zero. However, there is power dissipated by the load resistor, and its time average is exactly equal to the power supplied by the thermal bath. The exact power formula is similar to Nyquist's noise power formula, except that the rate of change of diode resistance significantly boosts the output power, and the movement of the graphene shifts the power spectrum to lower frequencies. We have calculated the equilibrium average of the power by asymptotic and numerical methods. Excellent agreement is found between experiment and theory.

Bifurcation analysis and phase diagram of a spin-string model with buckled states

We analyze a one-dimensional spin-string model, in which string oscillators are linearly coupled to their two nearest neighbors and to Ising spins representing internal degrees of freedom. String-spin coupling induces a long-range ferromagnetic interaction among spins that competes with a spin-spin antiferromagnetic coupling. As a consequence, the complex phase diagram of the system exhibits different flat rippled and buckled states, with first or second order transition lines between states. This complexity translates to the two-dimensional version of the model, whose numerical solution has been recently used to explain qualitatively the rippled to buckled transition observed in scanning tunneling microscopy experiments with suspended graphene sheets. Here we describe in detail the phase diagram of the simpler one-dimensional model and phase stability using bifurcation theory. This gives additional insight into the physical mechanisms underlying the different phases and the behavior observed in experiments.

Statistical mechanics of the Huxley-Simmons model

The chemomechanical model of Huxley and Simmons (HS) [A. F. Huxley and R. M. Simmons, Nature 233, 533 (1971)NATUAS0028-083610.1038/233533a0] provides a paradigmatic description of mechanically induced collective conformational changes relevant in a variety of biological contexts, from muscles power stroke and hair cell gating to integrin binding and hairpin unzipping. We develop a statistical mechanical perspective on the HS model by exploiting a formal analogy with a paramagnetic Ising model. We first study the equilibrium HS model with a finite number of elements and compute explicitly its mechanical and thermal properties. To model kinetics, we derive a master equation and solve it for several loading protocols. The developed formalism is applicable to a broad range of allosteric systems with mean-field interactions.

Spin-oscillator model for the unzipping of biomolecules by mechanical force

A spin-oscillator system models unzipping of biomolecules (such as DNA, RNA, or proteins) subject to an external force. The system comprises a macroscopic degree of freedom, represented by a one-dimensional oscillator, and internal degrees of freedom, represented by Glauber spins with nearest-neighbor interaction and a coupling constant proportional to the oscillator position. At a critical value F(c) of an applied external force F, the oscillator rest position (order parameter) changes abruptly and the system undergoes a first-order phase transition. When the external force is cycled at different rates, the extension given by the oscillator position exhibits a hysteresis cycle at high loading rates, whereas it moves reversibly over the equilibrium force-extension curve at very low loading rates. Under constant force, the logarithm of the residence time at the stable and metastable oscillator rest position is proportional to F-F(c) as in an Arrhenius law.