All Publications

In a recent paper [Phys. Rev. Lett. 129, 120601 (2022)0031-900710.1103/PhysRevLett.129.120601], we have shown that the dynamics of interfaces, in the symmetry-broken phase of the two-dimensional ferromagnetic quantum Ising model, displays a robust form of ergodicity breaking. In this paper, we elaborate more on the issue. First, we discuss two classes of initial states on the square lattice, the dynamics of which is driven by complementary terms in the effective Hamiltonian and may be solved exactly: (a) Strips of consecutive neighboring spins aligned in the opposite direction of the surrounding spins and (b) a large class of initial states, characterized by the presence of a well-defined "smooth"interface separating two infinitely extended regions with oppositely aligned spins. The evolution of the latter states can be mapped onto that of an effective one-dimensional fermionic chain, which is integrable in the infinite-coupling limit. In this case, deep connections with noteworthy results in mathematics emerge, as well as with similar problems in classical statistical physics. We present a detailed analysis of the evolution of these interfaces both on the lattice and in a suitable continuum limit, including the interface fluctuations and the dynamics of entanglement entropy. Second, we provide analytical and numerical evidence supporting the conclusion that the observed nonergodicity - arising from Stark localization of the effective fermionic excitations - persists away from the infinite-Ising-coupling limit, and we highlight the presence of a timescale T∼ecLlnL for the decay of a region of large linear size L. The implications of our work for the classic problem of the decay of a false vacuum are also discussed.

We study the nonequilibrium dynamics of two particles confined in two spatially separated harmonic potentials and linearly coupled to the same thermally fluctuating scalar field, a cartoon for optically trapped colloids in contact with a medium close to a continuous phase transition. When an external periodic driving is applied to one of these particles, a nonequilibrium periodic state is eventually reached in which their motion synchronizes thanks to the field-mediated effective interaction, a phenomenon already observed in experiments. We fully characterize the nonlinear response of the second particle as a function of the driving frequency, in particular far from the adiabatic regime in which the field can be assumed to relax instantaneously. We compare the perturbative, analytic solution to its adiabatic approximation, thus determining the limits of validity of the latter, and we qualitatively test our predictions against numerical simulations.

The effective dynamics of a colloidal particle immersed in a complex medium is often described in terms of an overdamped linear Langevin equation for its velocity with a memory kernel which determines the effective (time-dependent) friction and the correlations of fluctuations. Recently, it has been shown in experiments and numerical simulations that this memory may depend on the possible optical confinement the particle is subject to, suggesting that this description does not capture faithfully the actual dynamics of the colloid, even at equilibrium. Here, we propose a different approach in which we model the medium as a Gaussian field linearly coupled to the colloid. The resulting effective evolution equation of the colloidal particle features a non-linear memory term which extends previous models and which explains qualitatively the experimental and numerical evidence in the presence of confinement. This non-linear term is related to the correlations of the effective noise via a novel fluctuation-dissipation relation which we derive.

We study the statistics of the first-passage time of a single run-and-tumble particle (RTP) in one spatial dimension, with or without resetting, to a fixed target located at L>0. First, we compute the first-passage time distribution of a free RTP, without resetting or in a confining potential, but averaged over the initial position drawn from an arbitrary distribution p(x). Recent experiments used a noninstantaneous resetting protocol that motivated us to study in particular the case where p(x) corresponds to the stationary non-Boltzmann distribution of an RTP in the presence of a harmonic trap. This distribution p(x) is characterized by a parameter ν>0, which depends on the microscopic parameters of the RTP dynamics. We show that the first-passage time distribution of the free RTP, drawn from this initial distribution, develops interesting singular behaviors, depending on the value of ν. We then switch on resetting, mimicked by relaxation of the RTP in the presence of a harmonic trap. Resetting leads to a finite mean first-passage time and we study this as a function of the resetting rate for different values of the parameters ν and b=L/c, where c is the position of the right edge of the initial distribution p(x). In the diffusive limit of the RTP dynamics, we find a rich phase diagram in the (b,ν) plane, with an interesting reentrance phase transition. Away from the diffusive limit, qualitatively similar rich behaviors emerge for the full RTP dynamics.

We study the nonequilibrium evolution of coexisting ferromagnetic domains in the two-dimensional quantum Ising model - a setup relevant in several contexts, from quantum nucleation dynamics and false-vacuum decay scenarios to recent experiments with Rydberg-atom arrays. We demonstrate that the quantum-fluctuating interface delimiting a large bubble can be studied as an effective one-dimensional system through a "holographic"mapping. For the considered model, the emergent interface excitations map to an integrable chain of fermionic particles. We discuss how this integrability is broken by geometric features of the bubbles and by corrections in inverse powers of the ferromagnetic coupling, and provide a lower bound to the timescale after which the bubble is ultimately expected to melt. Remarkably, we demonstrate that a symmetry-breaking longitudinal field gives rise to a robust ergodicity breaking in two dimensions, a phenomenon underpinned by Stark many-body localization of the emergent fermionic excitations of the interface.

Modeling noisy oscillations of active systems is one of the current challenges in physics and biology. Because the physical mechanisms of such processes are often difficult to identify, we propose a linear stochastic model driven by a non-Markovian bistable noise that is capable of generating self-sustained periodic oscillation. We derive analytical predictions for most relevant dynamical and thermodynamic properties of the model. This minimal model turns out to describe accurately bistablelike oscillatory motion of hair bundles in bullfrog sacculus, extracted from experimental data. Based on and in agreement with these data, we estimate the power required to sustain such active oscillations to be of the order of 100 kBT per oscillation cycle.

We discuss how to calculate non-equilibrium universal amplitude ratios in the functional renormalization group approach, extending its applicability. In particular, we focus on the critical relaxation of the Ising model with non-conserved dynamics (model A) and calculate the universal amplitude ratio associated with the fluctuation-dissipation ratio of the order parameter, considering a critical quench from a high-temperature initial condition. Our predictions turn out to be in good agreement with previous perturbative renormalization-group calculations and Monte Carlo simulations.

We study the nonequilibrium relaxational dynamics of a probe particle linearly coupled to a thermally fluctuating scalar field and subject to a harmonic potential, which provides a cartoon for an optically trapped colloid immersed in a fluid close to its bulk critical point. The average position of the particle initially displaced from the position of mechanical equilibrium is shown to feature long-time algebraic tails as the critical point of the field is approached, the universal exponents of which are determined in arbitrary spatial dimensions. As expected, this behavior cannot be captured by adiabatic approaches which assumes fast field relaxation. The predictions of the analytic, perturbative approach are qualitatively confirmed by numerical simulations.

The interplay between cellular growth and cell-cell signaling is essential for the aggregation and proliferation of bacterial colonies, as well as for the self-organization of cell tissues. To investigate this interplay, we focus here on the collective properties of dividing chemotactic cell colonies by studying their long-time and large-scale dynamics through a renormalization group (RG) approach. The RG analysis reveals that a relevant but unconventional chemotactic interaction - corresponding to a polarity-induced mechanism - is generated by fluctuations at macroscopic scales, even when an underlying mechanism is absent at the microscopic level. This emerges from the interplay of the well-known Keller-Segel (KS) chemotactic nonlinearity and cell birth and death processes. At one-loop order, we find no stable fixed point of the RG flow equations. We discuss a connection between the dynamics investigated here and the celebrated Kardar-Parisi-Zhang (KPZ) equation with long-range correlated noise, which points at the existence of a strong-coupling, nonperturbative fixed point. ©

The collective and purely relaxational dynamics of quantum many-body systems after a quench at temperature T = 0, from a disordered state to various phases is studied through the exact solution of the quantum Langevin equation of the spherical and the O(n)-model in the limit n → ∞. The stationary state of the quantum dynamics is shown to be a non-equilibrium state. The quantum spherical and the quantum O(n)-model for n → ∞ are in the same dynamical universality class. The long-time behaviour of single-time and two-time correlation and response functions is analysed and the universal exponents which characterise quantum coarsening and quantum ageing are derived. The importance of the non-Markovian long-time memory of the quantum noise is elucidated by comparing it with an effective Markovian noise having the same scaling behaviour and with the case of non-equilibrium classical dynamics.

The critical Casimir force (CCF) arises from confining fluctuations in a critical fluid and thus it is a fluctuating quantity itself. While the mean CCF is universal, its (static) variance has previously been found to depend on the microscopic details of the system which effectively set a large-momentum cutoff in the underlying field theory, rendering it potentially large. This raises the question how the properties of the force variance are reflected in experimentally observable quantities, such as the thickness of a wetting film or the position of a suspended colloidal particle. Here, based on a rigorous definition of the instantaneous force, we analyze static and dynamic correlations of the CCF for a conserved fluid in film geometry for various boundary conditions within the Gaussian approximation. We find that the dynamic correlation function of the CCF is independent of the momentum cutoff and decays algebraically in time. Within the Gaussian approximation, the associated exponent depends only on the dynamic universality class but not on the boundary conditions. We furthermore consider a fluid film, the thickness of which can fluctuate under the influence of the time-dependent CCF. The latter gives rise to an effective non-Markovian noise in the equation of motion of the film boundary and induces a distinct contribution to the position variance. Within the approximations used here, at short times, this contribution grows algebraically in time whereas, at long times, it saturates and contributes to the steady-state variance of the film thickness.

The critical properties characterizing the formation of the Floquet time crystal in the prethermal phase are investigated analytically in the periodically driven O(N) model. In particular, we focus on the critical line separating the trivial phase with period synchronized dynamics and the absence of long-range spatial order from the nontrivial phase where long-range spatial order is accompanied by period-doubling dynamics. In the vicinity of the critical line, with a combination of dimensional expansion and exact solution for N→∞, we determine the exponent ν that characterizes the divergence of the spatial correlation length of the equal-time correlation functions, the exponent β characterizing the growth of the amplitude of the order parameter, as well as the initial-slip exponent θ of the aging dynamics when a quench is performed from deep in the trivial phase to the critical line. The exponents ν,β,θ are found to be identical to those in the absence of the drive. In addition, the functional form of the aging is found to depend on whether the system is probed at times that are small or large compared to the drive period. The spatial structure of the two-point correlation functions, obtained as a linear response to a perturbing potential in the vicinity of the critical line, is found to show algebraic decays that are longer ranged than in the absence of a drive, and besides being period doubled are also found to oscillate in space at the wave vector ω/(2v), v being the velocity of the quasiparticles, and ω being the drive frequency.

Critical Casimir forces emerge between objects, such as colloidal particles, whenever their surfaces spatially confine the fluctuations of the order parameter of a critical liquid used as a solvent. These forces act at short but microscopically large distances between these objects, reaching often hundreds of nanometers. Keeping colloids at such distances is a major experimental challenge, which can be addressed by the means of optical tweezers. Here, we review how optical tweezers have been successfully used to quantitatively study critical Casimir forces acting on particles in suspensions. As we will see, the use of optical tweezers to experimentally study critical Casimir forces can play a crucial role in developing nano-technologies, representing an innovative way to realize self-assembled devices at the nano- and microscale.

A generically observed mechanism that drives the self-organization of living systems is interaction via chemical signals among the individual elements - which may represent cells, bacteria, or even enzymes. Here we propose an unconventional mechanism for such interactions, in the context of chemotaxis, which originates from the polarity of the particles and which generalizes the well-known Keller-Segel interaction term. We study the resulting large-scale dynamical properties of a system of such chemotactic particles using the exact stochastic formulation of Dean and Kawasaki along with dynamical renormalization group analysis of the critical state of the system. At this critical point, an emergent "Galilean"symmetry is identified, which allows us to obtain the dynamical scaling exponents exactly. These exponents reveal superdiffusive density fluctuations and non-Poissonian number fluctuations. We expect our results to shed light on how molecular regulation of chemotactic circuits can determine large-scale behavior of cell colonies and tissues.

The periodically driven O(N) model is studied near the critical line separating a disordered paramagnetic phase from a period doubled phase, the latter being an example of a Floquet time crystal. The time evolution of one-point and two-point correlation functions are obtained within the Gaussian approximation and perturbatively in the drive amplitude. The correlations are found to show not only period doubling, but also power-law decays at large spatial distances. These features are compared with the undriven O(N) model, within the Gaussian approximation, in the vicinity of the paramagnetic-ferromagnetic critical point. The algebraic decays in space are found to be qualitatively different in the driven and the undriven cases. In particular, the spatiotemporal order of the Floquet time crystal leads to position-momentum and momentum-momentum correlation functions which are more long-ranged in the driven than in the undriven model. The light-cone dynamics associated with the correlation functions is also qualitatively different as the critical line of the Floquet time crystal shows a light cone with two distinct velocities, with the ratio of these two velocities scaling as the square-root of the dimensionless drive amplitude. The Floquet unitary, which describes the time evolution due to a complete cycle of the drive, is constructed for modes with small momenta compared to the drive frequency, but having a generic relationship with the square-root of the drive amplitude. At intermediate momenta, which are large compared to the square-root of the drive amplitude, the Floquet unitary is found to simply rotate the modes. On the other hand, at momenta which are small compared to the square-root of the drive amplitude, the Floquet unitary is found to primarily squeeze the modes, to an extent which increases upon increasing the wavelength of the modes, with a power-law dependence on it.

We investigate the Brownian diffusion of particles in one spatial dimension and in the presence of finite regions within which particles can either evaporate or be reset to a given location. For open boundary conditions, we highlight the appearance of a Brownian yet non-Gaussian diffusion: At long times, the particle distribution is non-Gaussian but its variance grows linearly in time. Moreover, we show that the effective diffusion coefficient of the particles in such systems is bounded from below by 1-2/π times their bare diffusion coefficient. For periodic boundary conditions, i.e., for diffusion on a ring with resetting, we demonstrate a "gauge invariance"of the spatial particle distribution for different choices of the resetting probability currents, in both stationary and nonstationary regimes. Finally, we apply our findings to a stochastic biophysical model for the motion of RNA polymerases during transcriptional pauses, deriving analytically the distribution of the length of cleaved RNA transcripts and the efficiency of RNA cleavage in backtrack recovery.

We study the dynamics of the statistics of the energy transferred across a point along a quantum chain which is prepared in the inhomogeneous initial state obtained by joining two identical semi-infinite parts thermalized at two different temperatures. In particular, we consider the transverse field Ising and harmonic chains as prototypical models of noninteracting fermionic and bosonic excitations, respectively. Within the so-called hydrodynamic limit of large space-time scales we first discuss the mean values of the energy density and current, and then, aiming at the statistics of fluctuations, we calculate exactly the scaled cumulant generating function of the transferred energy. From the latter, the evolution of the associated large deviation function is obtained. A natural interpretation of our results is provided in terms of a semiclassical picture of quasiparticles moving ballistically along classical trajectories. Similarities and differences between the transferred energy scaled cumulant and the large deviation functions in the cases of noninteracting fermions and bosons are discussed.

Confinement of excitations induces quasilocalized dynamics in disorder-free isolated quantum many-body systems in one spatial dimension. This occurrence is signaled by severe suppression of quantum correlation spreading and of entanglement growth, long-time persistence of spatial inhomogeneities, and long-lived coherent oscillations of local observables. In this work, we present a unified understanding of these dramatic effects. The slow dynamical behavior is shown to be related to the Schwinger effect in quantum electrodynamics. We demonstrate that it is quantitatively captured for long-time scales by effective Hamiltonians exhibiting Stark localization of excitations and weak growth of the entanglement entropy for arbitrary coupling strength. This analysis explains the phenomenology of real-time string dynamics investigated in a number of lattice gauge theories, as well as the anomalous dynamics observed in quantum Ising chains after quenches. Our findings establish confinement as a robust mechanism for hindering the approach to equilibrium in translationally invariant quantum statistical systems with local interactions.

The authors regret an error in the grant number for one of the authors in the Acknowledgements section. The Acknowledgements section should read as follows: This work was partially supported by the ERC Starting Grant ComplexSwimmers (grant no. 677511) and by Vetenskapsrådet (grant no. 2016-03523). A. C. acknowledges partial financial support from TUBITAK (grant no. 116F111). J. P. S. acknowledges partial financial support from FONDECYT (grant no. 1171013). The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and reade

Gauge theories are the cornerstone of our understanding of fundamental interactions among elementary particles. Their properties are often probed in dynamical experiments, such as those performed at ion colliders and high-intensity laser facilities. Describing the evolution of these strongly coupled systems is a formidable challenge for classical computers and represents one of the key open quests for quantum simulation approaches to particle physics phenomena. In this work, we show how recent experiments done on Rydberg atom chains naturally realize the real-time dynamics of a lattice gauge theory at system sizes at the boundary of classical computational methods. We prove that the constrained Hamiltonian dynamics induced by strong Rydberg interactions maps exactly onto the one of a U(1) lattice gauge theory. Building on this correspondence, we show that the recently observed anomalously slow dynamics corresponds to a string-inversion mechanism, reminiscent of the string breaking typically observed in gauge theories. This underlies the generality of this slow dynamics, which we illustrate in the context of one-dimensional quantum electrodynamics on the lattice. Within the same platform, we propose a set of experiments that generically show long-lived oscillations, including the evolution of particle-antiparticle pairs, and discuss how a tunable topological angle can be realized, further affecting the dynamics following a quench. Our work shows that the state of the art for quantum simulation of lattice gauge theories is at 51 qubits and connects the recently observed slow dynamics in atomic systems to archetypal phenomena in particle physics.