Search

We study the problem of minimising the connection time between non-equilibrium steady states of the Brownian Gyrator. This is a paradigmatic model in non-equilibrium statistical mechanics, an overdamped Brownian particle trapped in a two-dimensional elliptical potential, with the two degrees of freedom $(x,y)$ coupled to two, in principle different, thermal baths with temperatures $T_x$ and $T_y$, respectively. Application of Pontryagin's Maximum Principle reveals that shortest protocols belong to the boundaries of the control set defined by the limiting values of the parameters $(k,u)$ characterising the elliptical potential. We identify two classes of optimal minimum time protocols, i.e. brachistochrones: (i) regular bang-bang protocols, for which $(k,u)$ alternatively take their minimum and maximum values allowed, and (ii) infinitely degenerate singular protocols. We thoroughly investigate the minimum connection time over the brachistochrones in the limit of having infinite capacity for compression. A plethora of striking phenomena emerge: sets of states attained at null connection times, discontinuities in the connection time along adjacent target states, and the fact that, starting from a state in which the oscillators are coupled, uncoupled states are impossible to reach in a finite time., Comment: 31 pages, 8 figures

In the absence of external forcing, all trajectories on the phase plane of the van der Pol oscillator tend to a closed, periodic, trajectory -- the limit cycle -- after infinite time. Here, we drive the van der Pol oscillator with an external time-dependent force to reach the limit cycle in a given finite time. Specifically, we are interested in minimising the non-conservative contribution to the work when driving the system from a given initial point on the phase plane to any final point belonging to the limit cycle. There appears a speed limit inequality, which expresses a trade-off between the connection time and cost -- in terms of the non-conservative work. We show how the above results can be { generalized to the broader family of non-linear oscillators given by} the Li\'enard equation. Finally, we also look into the problem of minimising the total work done by the external force., Comment: 13 pages, 6 figures; minor revision, new perspectives added; accepted for publication in Chaos

In this perspective paper, we look into memory effects in out-of-equilibrium systems. To be concrete, we exemplify memory effects with the paradigmatic case of granular fluids, although extensions to other contexts such as molecular fluids with non-linear drag are also considered. The focus is put on two archetypal memory effects: the Kovacs and Mpemba effects. In brief, the first is related to imperfectly reaching a steady state -- either equilibrium or non-equilibrium, whereas the second is related to reaching a steady state faster despite starting further. Connections to optimal control theory thus naturally emerge and are briefly discussed, Comment: Perspective paper for EPL, 7 pages, 6 figures; accepted version, minor revision

Scanning tunneling microscopy experiments have revealed an 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 analyse 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. Minimisation of this free energy gives rise to the emergence of different mechanical phases, whose thermodynamic stability is thoroughly analysed, both analytically and numerically. All phases are characterised 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., Comment: 18 pages, 8 figures

Obtaining adiabatic processes that connect equilibrium states in a given time represents a challenge for mesoscopic systems. In this paper, we explicitly show how to build these finite-time adiabatic processes for an overdamped Brownian particle in an arbitrary potential, a system that is relevant both at the conceptual and the practical level. This is achieved by jointly engineering the time evolutions of the binding potential and the fluid temperature. Moreover, we prove that the second principle imposes a speed limit for such adiabatic transformations: there appears a minimum time to connect the initial and final states. This minimum time can be explicitly calculated for a general compression/decompression situation., Comment: Main text: 5 pages; 18 pages with appendices and references; major revision with results for a general non-linear potential and study of fluctuations added; Physical Review E in press

In kinetic theory, a system is usually described by its one-particle distribution function $f(\mathbf{r},\mathbf{v},t)$, such that $f(\mathbf{r},\mathbf{v},t)d\mathbf{r} d\mathbf{v}$ is the fraction of particles with positions and velocities in the intervals $(\mathbf{r}, \mathbf{r}+d\mathbf{r})$ and $(\mathbf{v}, \mathbf{v}+d\mathbf{v})$, respectively. Therein, global stability and the possible existence of an associated Lyapunov function or $H$-theorem are open problems when non-conservative interactions are present, as in granular fluids. Here, we address this issue in the framework of a lattice model for granular-like velocity fields. For a quite general driving mechanism, including both boundary and bulk driving, we show that the steady state reached by the system in the long time limit is globally stable. This is done by proving analytically that a certain $H$-functional is non-increasing in the long time limit. Moreover, for two specific energy injection mechanisms, we are able to demonstrate that the proposed $H$-functional is non-increasing for all times. Also, we put forward a proof that clearly illustrates why the "classical" Boltzmann functional $H_{B}[f]=\int\! d\mathbf{r} \, d\mathbf{v} f(\mathbf{r},\mathbf{v},t) \ln f(\mathbf{r},\mathbf{v},t)$ is inadequate for systems with non-conservative interactions. Not only is this done for the simplified kinetic description that holds in the lattice models analysed here but also for a general kinetic equation, like Boltzmann's or Enskog's., Comment: 13 pages, 5 figures; minor revision, accepted for publication in PRE

We comment on the paper "Critique and correction of the currently accepted solution of the infinite spherical well in quantum mechanics" by Huang Young-Sea and Thomann Hans-Rudolph, EPL 115, 60001 (2016) ., Comment: 2 pages; Submitted to the Comments Section of EPL

Long-range spatial correlations in the velocity and energy fields of a granular fluid are discussed in the framework of a 1d lattice model. The dynamics of the velocity field occurs through nearest-neighbour inelastic collisions that conserve momentum but dissipate energy. A set of equations for the fluctuating hydrodynamics of the velocity and energy mesoscopic fields give a first approximation for (i) the velocity structure factor and (ii) the finite-size correction to the Haff law, both in the homogeneous cooling regime. At a more refined level, we have derived the equations for the two-site velocity correlations and the total energy fluctuations. First, we seek a perturbative solution thereof, in powers of the inverse of system size. On the one hand, when scaled with the granular temperature, the velocity correlations tend to a stationary value in the long time limit. On the other hand, the scaled standard deviation of the total energy diverges, that is, the system shows multiscaling. Second, we find an exact solution for the velocity correlations in terms of the spectrum of eigenvalues of a certain matrix. The results of numerical simulations of the microscopic model confirm our theoretical results, including the above described multiscaling phenomenon.

The dependence of the unfolding pathway of proteins on the pulling speed is investigated. This is done by introducing a simple one-dimensional chain comprising $N$ units, with different characteristic bistable free energies. These units represent either each of the modules in a modular protein or each of the intermediate "unfoldons" in a protein domain, which can be either folded or unfolded. The system is pulled by applying a force to the last unit of the chain, and the units unravel following a preferred sequence. We show that the unfolding sequence strongly depends on the pulling velocity $v_{p}$. In the simplest situation, there appears a critical pulling speed $v_{c}$: for pulling speeds $v_{p}v_{c}$ it is the pulled unit that unfolds first. By means of a perturbative expansion, we find quite an accurate expression for this critical velocity., Comment: accepted for publication in JSTAT