Source author record

Juan P. Garrahan

Juan P. Garrahan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

77works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

77 published item(s)

preprint2022arXiv

Finite time large deviations via matrix product states

Recent work has shown the effectiveness of tensor network methods for computing large deviation functions in constrained stochastic models in the infinite time limit. Here we show that these methods can also be used to study the statistics of dynamical observables at arbitrary finite time. This is a harder problem because, in contrast to the infinite time case where only the extremal eigenstate of a tilted Markov generator is relevant, for finite time the whole spectrum plays a role. We show that finite time dynamical partition sums can be computed efficiently and accurately in one dimension using matrix product states, and describe how to use such results to generate rare event trajectories on demand. We apply our methods to the Fredrickson-Andersen (FA) and East kinetically constrained models, and to the symmetric simple exclusion process (SSEP), unveiling dynamical phase diagrams in terms of counting field and trajectory time. We also discuss extensions of this method to higher dimensions.

preprint2022arXiv

Topological phases in the dynamics of the simple exclusion process

We study the dynamical large deviations of the classical stochastic symmetric simple exclusion process (SSEP) by means of numerical matrix product states. We show that for half-filling, long-time trajectories with a large enough imbalance between the number hops in even and odd bonds of the lattice belong to distinct symmetry protected topological (SPT) phases. Using tensor network techniques, we obtain the large deviation (LD) phase diagram in terms of counting fields conjugate to the dynamical activity and the total hop imbalance. We show the existence of high activity trivial and non-trivial SPT phases (classified according to string-order parameters) separated by either a critical phase or a critical point. Using the leading eigenstate of the tilted generator, obtained from infinite-system density matrix renormalisation group (DMRG) simulations, we construct a near-optimal dynamics for sampling the LDs, and show that the SPT phases manifest at the level of rare stochastic trajectories. We also show how to extend these results to other filling fractions, and discuss generalizations to asymmetric SEPs.

preprint2021arXiv

Exact solution of the Rule 150 reversible cellular automaton

We study the dynamics of the Rule 150 reversible cellular automaton (RCA). This is a one-dimensional lattice system of binary variables with synchronous (Floquet) dynamics, corresponding to a bulk deterministic and reversible discrete version of the kinetically constrained XOR-Fredrickson-Andersen model, whereby the local dynamics is restricted: a site flips if and only if the states of its neighbouring sites are different from each other. Like other RCA which have been studied recently, such as Rule 54 and Rule 201, Rule 150 is integrable, however, in contrast is noninteracting. In particular, the emergent quasiparticles - the domain walls - behave as free fermions. This then allows us to solve the model by means of matrix product ansätze. We find the exact equilibrium and nonequilibrium stationary states for systems with closed (periodic) and open (stochastic) boundaries, respectively, resolve the full spectrum of the time evolution operator and, therefore, gain access to the relaxation dynamics, and obtain the exact large deviation statistics of dynamical observables in the long time limit.

preprint2021arXiv

Reinforcement learning of rare diffusive dynamics

We present a method to probe rare molecular dynamics trajectories directly using reinforcement learning. We consider trajectories that are conditioned to transition between regions of configuration space in finite time, like those relevant in the study of reactive events, as well as trajectories exhibiting rare fluctuations of time-integrated quantities in the long time limit, like those relevant in the calculation of large deviation functions. In both cases, reinforcement learning techniques are used to optimize an added force that minimizes the Kullback-Leibler divergence between the conditioned trajectory ensemble and a driven one. Under the optimized added force, the system evolves the rare fluctuation as a typical one, affording a variational estimate of its likelihood in the original trajectory ensemble. Low variance gradients employing value functions are proposed to increase the convergence of the optimal force. The method we develop employing these gradients leads to efficient and accurate estimates of both the optimal force and the likelihood of the rare event for a variety of model systems.

preprint2020arXiv

A Deep Learning Functional Estimator of Optimal Dynamics for Sampling Large Deviations

In stochastic systems, numerically sampling the relevant trajectories for the estimation of the large deviation statistics of time-extensive observables requires overcoming their exponential (in space and time) scarcity. The optimal way to access these rare events is by means of an auxiliary dynamics obtained from the original one through the so-called ``generalised Doob transformation''. While this optimal dynamics is guaranteed to exist its use is often impractical, as to define it requires the often impossible task of diagonalising a (tilted) dynamical generator. While approximate schemes have been devised to overcome this issue they are difficult to automate as they tend to require knowledge of the systems under study. Here we address this problem from the perspective of deep learning. We devise an iterative semi-supervised learning scheme which converges to the optimal or Doob dynamics with the clear advantage of requiring no prior knowledge of the system. We test our method in a paradigmatic statistical mechanics model with non-trivial dynamical fluctuations, the fully packed classical dimer model on the square lattice, showing that it compares favourably with more traditional approaches. We discuss broader implications of our results for the study of rare dynamical trajectories.

preprint2020arXiv

A reinforcement learning approach to rare trajectory sampling

Very often when studying non-equilibrium systems one is interested in analysing dynamical behaviour that occurs with very low probability, so called rare events. In practice, since rare events are by definition atypical, they are often difficult to access in a statistically significant way. What are required are strategies to "make rare events typical" so that they can be generated on demand. Here we present such a general approach to adaptively construct a dynamics that efficiently samples atypical events. We do so by exploiting the methods of reinforcement learning (RL), which refers to the set of machine learning techniques aimed at finding the optimal behaviour to maximise a reward associated with the dynamics. We consider the general perspective of dynamical trajectory ensembles, whereby rare events are described in terms of ensemble reweighting. By minimising the distance between a reweighted ensemble and that of a suitably parametrised controlled dynamics we arrive at a set of methods similar to those of RL to numerically approximate the optimal dynamics that realises the rare behaviour of interest. As simple illustrations we consider in detail the problem of excursions of a random walker, for the case of rare events with a finite time horizon; and the problem of a studying current statistics of a particle hopping in a ring geometry, for the case of an infinite time horizon. We discuss natural extensions of the ideas presented here, including to continuous-time Markov systems, first passage time problems and non-Markovian dynamics.

preprint2020arXiv

A Tensor Network Approach to Finite Markov Decision Processes

Tensor network (TN) techniques - often used in the context of quantum many-body physics - have shown promise as a tool for tackling machine learning (ML) problems. The application of TNs to ML, however, has mostly focused on supervised and unsupervised learning. Yet, with their direct connection to hidden Markov chains, TNs are also naturally suited to Markov decision processes (MDPs) which provide the foundation for reinforcement learning (RL). Here we introduce a general TN formulation of finite, episodic and discrete MDPs. We show how this formulation allows us to exploit algorithms developed for TNs for policy optimisation, the key aim of RL. As an application we consider the issue - formulated as an RL problem - of finding a stochastic evolution that satisfies specific dynamical conditions, using the simple example of random walk excursions as an illustration.

preprint2020arXiv

Matrix product state of multi-time correlations

For an interacting spatio-temporal lattice system we introduce a formal way of expressing multi-time correlation functions of local observables located at the same spatial point with a time state, i.e. a statistical distribution of configurations observed along a time lattice. Such a time state is defined with respect to a particular equilibrium state that is invariant under space and time translations. The concept is developed within the Rule 54 reversible cellular automaton, for which we explicitly construct a matrix product form of the time state, with matrices that act on the 3-dimensional auxiliary space. We use the matrix-product state to express equal-space time-dependent density-density correlation function, which, for special maximum-entropy values of equilibrium parameters, agrees with the previous results. Additionally, we obtain an explicit expression for the probabilities of observing all multi-time configurations, which enables us to study distributions of times between consecutive excitations and prove the absence of decoupling of timescales in the Rule 54 model.

preprint2020arXiv

Quantum East model: localization, non-thermal eigenstates and slow dynamics

We study in detail the properties of the quantum East model, an interacting quantum spin chain inspired by simple kinetically-constrained models of classical glasses. Through a combination of analytics, exact diagonalization and tensor-network methods we show the existence of a transition, from a fast to a slow thermalization regime, which manifests itself throughout the spectrum. On the slow side, by exploiting the localization of the ground state and the form of the Hamiltonian, we explicitly construct a large (exponential in size) number of non-thermal states which become exact finite-energy-density eigenstates in the large-size limit, as expected for a true phase transition. A ``super-spin'' generalization allows us to find a further large class of area-law states proved to display very slow relaxation. These states retain memory of their initial conditions for extremely long times. Our numerical analysis reveals that the localization properties are not limited to the ground state and that} many eigenstates have large overlap with product states and can be approximated well by matrix product states at arbitrary energy densities. The mechanism that induces localization to the ground state, and hence the non-thermal behavior of the system, can be extended to a wide range of models including a number of simple spin chains. We discuss implications of our results for slow thermalization and non-ergodicity more generally in disorder-free systems with constraints and we give numerical evidence that these results may be extended to two dimensional systems.

preprint2020arXiv

Symmetry-induced fluctuation relations for dynamical observables irrespective of their behaviour under time-reversal

We extend previous work to describe a class of fluctuation relations (FRs) that emerge as a consequence of symmetries at the level of stochastic trajectories in Markov chains. We prove that given such a symmetry, and for a suitable dynamical observable, it is always possible to obtain a FR under a biased dynamics corresponding to the so-called generalized Doob transform. The general transformations of the dynamics that we consider go beyond time-reversal or spatial isometries, and an implication is the existence of FRs for observables irrespective of their behaviour under time-reversal, for example for time-symmetric observables rather than currents. We further show how to deduce in the long-time limit these FRs from the symmetry properties of the generator of the dynamics. We illustrate our results with four examples that highlight the novel features of our work.

preprint2019arXiv

A non-equilibrium quantum many-body Rydberg atom engine

The standard approach to quantum engines is based on equilibrium systems and on thermodynamic transformations between Gibbs states. However, non-equilibrium quantum systems offer enhanced experimental flexibility in the control of their parameters and, if used as engines, a more direct interpretation of the type of work they deliver. Here we introduce an out-of-equilibrium quantum engine inspired by recent experiments with cold atoms. Our system is connected to a single environment and produces mechanical work from many-body interparticle interactions arising between atoms in highly excited Rydberg states. As such, it is not a heat engine but an isothermal one. We perform many-body simulations to show that this system can produce work. The setup we introduce and investigate represents a promising platform for devising new types of microscopic machines and for exploring quantum effects in thermodynamic processes.

preprint2019arXiv

Trajectory phase transitions in non-interacting spin systems

We show that a collection of independent Ising spins evolving stochastically can display surprisingly large fluctuations towards ordered behaviour, as quantified by certain types of time-integrated plaquette observables, despite the underlying dynamics being non-interacting. In the large deviation (LD) regime of long times and large system size, this can give rise to a phase transition in trajectory space. As a non-interacting system we consider a collection of spins undergoing single spin-flip dynamics at infinite-temperature. For the dynamical observables we study, the associated tilted generators have an exact and explicit spin-plaquette duality. Such setup suggests the existence of a transition (in the large size limit) at the self-dual point of the tilted generator. The nature of the LD transition depends on the observable. We consider explicitly two situations: (i) for a pairwise bond observable the LD transition is continuous, and equivalent to that of the transverse field Ising model; (ii) for a higher order plaquette observable, in contrast, the LD transition is first order. Case (i) is easy to prove analytically, while we confirm case (ii) numerically via an efficient trajectory sampling scheme that exploits the non-interacting nature of the original dynamics.

preprint2016arXiv

Classical stochastic dynamics and continuous matrix product states: gauge transformations, conditioned and driven processes, and equivalence of trajectory ensembles

Borrowing ideas from open quantum systems, we describe a formalism to encode ensembles of trajectories of classical stochastic dynamics in terms of continuous matrix product states (cMPSs). We show how to define in this approach "biased" or "conditioned" ensembles where the probability of trajectories is biased from that of the natural dynamics by some condition on trajectory observables. In particular, we show that the generalised Doob transform which maps a conditioned process to an equivalent "auxiliary" or "driven" process (one where the same conditioned set of trajectories is generated by a proper stochastic dynamics) is just a gauge transformation of the corresponding cMPS. We also discuss how within this framework one can easily prove properties of the dynamics such as trajectory ensemble equivalence and fluctuation theorems.

preprint2016arXiv

Emergence of cooperative dynamics in fully packed classical dimers

We study the behavior of classical dimer coverings of the square lattice - a paradigmatic model for systems subject to constraints - evolving under local stochastic dynamics, by means of Monte Carlo simulations and theoretical arguments. We observe clear signatures of correlated dynamics in both global and local observables and over a broad range of time scales, indicating a breakdown of the simple continuum description that approximates well the statics. We show that this collective dynamics can be understood in terms of one-dimensional "strings" of high mobility, which govern both local and long-wavelength dynamical properties. We introduce a coarse-grained description of the strings, based on the Edwards-Wilkinson model, which leads to exact results in the limit of low string density and provides a detailed qualitative understanding of the dynamics in all flux sectors. We discuss the implications of our results for the dynamics of constrained systems more generally.

preprint2016arXiv

Front propagation versus bulk relaxation in the annealing dynamics of a kinetically constrained model of ultrastable glasses

Glasses prepared by physical vapour deposition have been shown to be remarkably more stable than those prepared by standard cooling protocols, with properties that appear to be similar to systems aged for extremely long times. When subjected to a rapid rise in temperature, ultrastable glasses anneal towards the liquid in a qualitatively different manner than ordinary glasses, with the seeming competition of different timescales and lengthscales. We numerically reproduce the phenomenology of ultrastable glass annealing with a kinetically constrained model, a three dimensional East model with soft constraints, in a setting where the bulk is in an ultrastable configuration and a free surface is permanently excited. Annealing towards the liquid state is given by the competition between the ballistic propagation of a front from the free surface and a much slower nucleation-like relaxation in the bulk. The crossover between these mechanisms also explains the change in behaviour with film thickness seen experimentally.

preprint2016arXiv

Metastability in an open quantum Ising model

We apply a recently developed theory for metastability in open quantum systems to a one-dimensional dissipative quantum Ising model. Earlier results suggest this model features either a non-equilibrium phase transition or a smooth but sharp crossover, where the stationary state changes from paramagnetic to ferromagnetic, accompanied by strongly intermittent emission dynamics characteristic of first-order coexistence between dynamical phases. We show that for a range of parameters close to this transition/crossover point the dynamics of the finite system displays pronounced metastability, i.e., the system relaxes first to long-lived metastable states, before eventual relaxation to the true stationary state. From the spectral properties of the quantum master operator we characterise the low-dimensional manifold of metastable states, which are shown to be probability mixtures of two, paramagnetic and ferromagnetic, metastable phases. We also show that for long times the dynamics can be approximated by a classical stochastic dynamics between the metastable phases that is directly related to the intermittent dynamics observed in quantum trajectories and thus the dynamical phases.

preprint2016arXiv

Non-equilibrium fluctuations and metastability arising from non-additive interactions in dissipative multi-component Rydberg gases

We study the out-of-equilibrium dynamics of dissipative gases of atoms excited to two or more high-lying Rydberg states. This situation bears interesting similarities to classical binary (in general $p$-ary) mixtures of particles. The effective forces between the components are determined by the inter-level and intra-level interactions of Rydberg atoms. These systems permit to explore new parameter regimes which are physically inaccessible in a classical setting, for example one in which the mixtures exhibit non-additive interactions. In this situation the out-of-equilibrium evolution is characterized by the formation of metastable domains that reach partial equilibration long before the attainment of stationarity. In experimental settings with mesoscopic sizes, this collective behavior may in fact take the appearance of dynamic symmetry breaking.

preprint2016arXiv

Non-equilibrium universality in the dynamics of dissipative cold atomic gases

The theory of continuous phase transitions predicts the universal collective properties of a physical system near a critical point, which for instance manifest in characteristic power-law behaviours of physical observables. The well-established concept at or near equilibrium, universality, can also characterize the physics of systems out of equilibrium. The most fundamental instance of a genuine non-equilibrium phase transition is the directed percolation universality class, where a system switches from an absorbing inactive to a fluctuating active phase. Despite being known for several decades it has been challenging to find experimental systems that manifest this transition. Here we show theoretically that signatures of the directed percolation universality class can be observed in an atomic system with long range interactions. Moreover, we demonstrate that even mesoscopic ensembles --- which are currently studied experimentally --- are sufficient to observe traces of this non-equilibrium phase transition in one, two and three dimensions.

preprint2016arXiv

Robustness of many-body localization in the presence of dissipation

Many-body localization (MBL) has emerged as a novel paradigm for robust ergodicity breaking in closed quantum many-body systems. However, it is not yet clear to which extent MBL survives in the presence of dissipative processes induced by the coupling to an environment. Here we study heating and ergodicity for a paradigmatic MBL system---an interacting fermionic chain subject to quenched disorder---in the presence of dephasing. We find that, even though the system is eventually driven into an infinite-temperature state, heating as monitored by the von Neumann entropy can progress logarithmically slowly, implying exponentially large time scales for relaxation. This slow loss of memory of initial conditions make signatures of non-ergodicity visible over a long, but transient, time regime. We point out a potential controlled realization of the considered setup with cold atomic gases held in optical lattices.

preprint2016arXiv

Self-similar non-equilibrium dynamics of a many-body system with power-law interactions

The influence of power-law interactions on the dynamics of many-body systems far from equilibrium is much less explored than their effect on static and thermodynamic properties. To gain insight into this problem we introduce and analyze here an out-of-equilibrium deposition process in which the deposition rate of a given particle depends as a power-law on the distance to previously deposited particles. This model draws its relevance from recent experimental progress in the domain of cold atomic gases which are studied in a setting where atoms that are excited to high-lying Rydberg states interact through power-law potentials that translate into power-law excitation rates. The out-of-equilibrium dynamics of this system turns out to be surprisingly rich. It features a self-similar evolution which leads to a characteristic power-law time dependence of observables such as the particle concentration, and results in a scale invariance of the structure factor. Our findings show that in dissipative Rydberg gases out of equilibrium the characteristic distance among excitations --- often referred to as the blockade radius --- is not a static but rather a dynamic quantity.

preprint2016arXiv

Signatures of many-body localisation in a system without disorder and the relation to a glass transition

We study a quantum spin system with local bilinear interactions and without quenched disorder which seems to display characteristic signatures of a many-body localisation (MBL) transition. From direct diagonalisation of small systems, we find a change in certain dynamical and spectral properties at a critical value of a coupling, from those characteristic of a thermalising phase to those characteristic of a MBL phase. The system we consider is known to have a quantum phase transition in its ground-state in the limit of large size, related to a first-order active-to-inactive phase transition in the stochastic trajectories of an associated classical model of glasses. Our results here suggest that this transition is present throughout the spectrum of the system in the large size limit. These findings may help understand the connection between MBL and structural glass transitions.

preprint2016arXiv

Towards a theory of metastability in open quantum dynamics

By generalising concepts from classical stochastic dynamics, we establish the basis for a theory of metastability in Markovian open quantum systems. Partial relaxation into long-lived metastable states - distinct from the asymptotic stationary state - is a manifestation of a separation of timescales due to a splitting in the spectrum of the generator of the dynamics. We show here how to exploit this spectral structure to obtain a low dimensional approximation to the dynamics in terms of motion in a manifold of metastable states constructed from the low-lying eigenmatrices of the generator. We argue that the metastable manifold is in general composed of disjoint states, noiseless subsystems and decoherence-free subspaces.

preprint2015arXiv

A generic map from non-Lindblad to Lindblad master equations

Many current problems of interest in quantum non-equilibrium are described by time-local master equations (TLMEs) for the density matrix that are not of the Lindblad form, that is, that are not strictly probability conserving and/or Markovian. Here we describe an generic approach by which the system of interest that obeys the TLME is coupled to an ancilla, such that the dynamics of the combined system-plus-ancilla is Markovian and thus described by a Lindblad equation. This in turn allows us to recover the properties of the original TLME dynamics from a physical unravelling of this associated Lindblad dynamics. We discuss applications of this generic mapping in two areas of current interest. The first one is that of "thermodynamics of trajectories", where non-Lindblad master equations encode the large-deviation properties of the dynamics, and we show that the relevant large-deviation functions (i.e. dynamical free-energies) can be recovered from appropriate observables of the ancilla. The second one is that of quantum filters, where we show tracking a quantum system undergoing a continuous homodyne measurement with another quantum system of the same size will inherently be inefficient in our framework.

preprint2015arXiv

Crystalline structures and frustration in a two-component Rydberg gas

We study the static behavior of a gas of atoms held in a one-dimensional lattice where two distinct electronically high-lying Rydberg states are simultaneously excited by laser light. We focus on a situation where interactions of van-der-Waals type take place only among atoms that are in the same Rydberg state. We analytically investigate at first the so-called classical limit of vanishing laser driving strength. We show that the system exhibits a surprisingly complex ground state structure with a sequence of compatible to incompatible transitions. The incompatibility between the species leads to mutual frustration, a feature which pertains also in the quantum regime. We perform an analytical and numerical investigation of these features and present an approximative description of the system in terms of a Rokhsar-Kivelson Hamiltonian which permits the analytical understanding of the frustration effects even beyond the classical limit.

preprint2015arXiv

Dynamical symmetries and crossovers in a three-spin system with collective dissipation

We consider the non-equilibrium dynamics of a simple system consisting of interacting spin-$1/2$ particles subjected to a collective damping. The model is close to situations that can be engineered in hybrid electro/opto-mechanical settings. Making use of large-deviation theory, we find a Gallavotti-Cohen symmetry in the dynamics of the system as well as evidence for the coexistence of two dynamical phases with different activity levels. We show that additional damping processes smoothen out this behavior. Our analytical results are backed up by Monte Carlo simulations that reveal the nature of the trajectories contributing to the different dynamical phases.

preprint2015arXiv

Dynamics of many-body localisation in a translation invariant quantum glass model

We study the real-time dynamics of a translationally invariant quantum spin chain, based on the East kinetically constrained glass model, in search for evidence of many-body localisation in the absence of disorder. Numerical simulations indicate a change, controlled by a coupling parameter, from a regime of fast relaxation---corresponding to thermalisation---to a regime of very slow relaxation. This slowly relaxing regime is characterised by dynamical features usually associated with non-ergodicity and many-body localisation (MBL): memory of initial conditions, logarithmic growth of entanglement entropy, and non-exponential decay of time-correlators. We show that slow relaxation is a consequence of sensitivity to spatial fluctuations in the initial state. While numerics indicate that certain relaxation timescales grow markedly with size, our finite size results are consistent both with an MBL transition, expected to only occur in disordered systems, or with a pronounced quasi-MBL crossover.

preprint2015arXiv

Equivalence of matrix product ensembles of trajectories in open quantum systems

The equivalence of thermodynamic ensembles is at the heart of statistical mechanics and central to our understanding of equilibrium states of matter. Recently, it has been shown that there is a formal connection between the dynamics of open quantum systems and the statistical mechanics in an extra dimension. This is established through the fact that an open system dynamics generates a Matrix Product state (MPS) which encodes the set of all possible quantum jump trajectories and permits the construction of generating functions in the spirit of thermodynamic partition functions. In the case of continuous-time Markovian evolution, such as that generated by a Lindblad master equation, the corresponding MPS is a so-called continuous MPS which encodes the set of continuous measurement records terminated at some fixed total observation time. Here we show that if one instead terminates trajectories after a fixed total number of quantum jumps, e.g. emission events into the environment, the associated MPS is discrete. This establishes an interesting analogy: The continuous and discrete MPS correspond to different ensembles of quantum trajectories, one characterised by total time the other by total number of quantum jumps. Hence they give rise to quantum versions of different thermodynamic ensembles, akin to "grand-canonical" and "isobaric", but for trajectories. Here we prove that these trajectory ensembles are equivalent in a suitable limit of long time or large number of jumps. This is in direct analogy to equilibrium statistical mechanics where equivalence between ensembles is only strictly established in the thermodynamic limit. An intrinsic quantum feature is that the equivalence holds only for all observables that commute with the number of quantum jumps.

preprint2015arXiv

Intermittency and dynamical Lee-Yang zeros of open quantum systems

We use high-order cumulants to investigate the Lee-Yang zeros of generating functions of dynamical observables in open quantum systems. At long times the generating functions take on a large deviation form with singularities of the associated cumulant generating functions - or dynamical free energies - signifying phase transitions in the ensemble of dynamical trajectories. We consider a driven three-level system as well as the dissipative Ising model. Both systems exhibit dynamical intermittency in the statistics of quantum jumps. From the short-time behavior of the dynamical Lee-Yang zeros we identify critical values of the counting field which we attribute to the observed intermittency and dynamical phase co-existence. Furthermore, for the dissipative Ising model we construct a trajectory phase diagram and estimate the value of the transverse field where the stationary state changes from being ferromagnetic (inactive) to being paramagnetic (active).

preprint2015arXiv

Overlap and activity glass transitions in plaquette spin models with hierarchical dynamics

We consider thermodynamic and dynamic phase transitions in plaquette spin models of glasses. The thermodynamic transitions involve coupled (annealed) replicas of the model. We map these coupled-replica systems to a single replica in a magnetic field, which allows us to analyse the resulting phase transitions in detail. For the triangular plaquette model (TPM), we find for the coupled-replica system a phase transition between high- and low-overlap phases, occuring at a coupling eps*(T), which vanishes in the low-temperature limit. Using computational path sampling techniques, we show that a single TPM also displays space-time transitions between active and inactive dynamical phases. These first-order dynamical transitions occur at a critical counting field s_c(T)>=0 that appears to vanish at zero temperature, in a manner reminiscent of the thermodynamic overlap transition. In order to extend the ideas to three dimensions we introduce the square pyramid model which also displays both overlap and activity transitions. We discuss a possible common origin of these various phase transitions, based on long-lived (metastable) glassy states.

preprint2015arXiv

Phase transition for quenched coupled replicas in a plaquette spin model of glasses

We study a three-dimensional plaquette spin model whose low temperature dynamics is glassy, due to localised defects and effective kinetic constraints. While the thermodynamics of this system is smooth at all temperatures, we show that coupling it to a second system with a fixed (quenched) configuration can lead to a phase transition, at finite coupling. The order parameter is the overlap between the copies, and the transition is between phases of low and high overlap. We find critical points whose properties are consistent with random-field Ising universality. We analyse the interfacial free energy cost between the high- and low-overlap states that coexist at (and below) the critical point, and we use this cost as the basis for a finite-size scaling analysis. We discuss these results in the context of mean-field and dynamical facilitation theories of the glass transition.

preprint2014arXiv

Comment on "Abrupt transition in the structural formation of interconnected networks", F. Radicchi and A. Arenas, Nature Phys. 9, 717 (2013)

A recent paper [F. Radicchi and A. Arenas, Nature Phys. 9, 717 (2013)] presented the finding of an abrupt transition in the structure of interconnected networks. This transition was said to be generic and to occur even in networks of finite size. Furthermore, it was remarked that this singular behaviour could be understood in the spirit of a first-order phase transition. We show here that the generic singularity found in that paper is a trivial consequence of the reducibility of the "supra-Laplacian" operator studied. The singular changes observed are therefore not related to any collective abrupt structural transformation in the interconnected networks.

preprint2014arXiv

Common physical framework explains phase behavior and dynamics of atomic, molecular and polymeric network-formers

We show that the self-assembly of a diverse collection of building blocks can be understood within a common physical framework. These building blocks, which form periodic honeycomb networks and nonperiodic variants thereof, range in size from atoms to micron-scale polymers, and interact through mechanisms as different as hydrogen bonds and covalent forces. A combination of statistical mechanics and quantum mechanics shows that one can capture the physics that governs assembly of these networks by resolving only the geometry and strength of building block interactions. The resulting framework reproduces a broad range of phenomena seen experimentally, including periodic and nonperiodic networks in thermal equilibrium, and nonperiodic supercooled and glassy networks away from equilibrium. Our results show how simple `design criteria' control assembly of a wide variety of networks, and suggest that kinetic trapping can be a useful way of making functional assemblies.

preprint2014arXiv

Cumulants of time-integrated observables of closed quantum systems and PT-symmetry, with an application to the quantum Ising chain

We study the connection between the cumulants of a time-integrated observable of a quantum system and the PT-symmetry properties of the non-Hermitian deformation of the Hamiltonian from which the generating function of these cumulants is obtained. This non-Hermitian Hamiltonian can display regimes of broken and of unbroken PT-symmetry, depending on the parameters of the problem and on the counting field that sets the strength of the non-Hermitian perturbation. This in turn determines the analytic structure of the long-time cumulant generating function (CGF) for the time-integrated observable. We consider in particular the case of the time-integrated (longitudinal) magnetisation in the one-dimensional Ising model in a transverse field. We show that its long-time CGF is singular on a curve in the magnetic field/counting field plane that delimits a regime where PT-symmetry is spontaneously broken (which includes the static ferromagnetic phase), from one where it is preserved (which includes the static paramagnetic phase). In the paramagnetic phase, conservation of PT -symmetry implies that all cumulants are sub-linear in time, a behaviour usually associated to the absence of decorrelation.

preprint2014arXiv

Dynamical phase transitions as a resource for quantum enhanced metrology

We consider the general problem of estimating an unknown control parameter of an open quantum system. We establish a direct relation between the evolution of both system and environment and the precision with which the parameter can be estimated. We show that when the open quantum system undergoes a first-order dynamical phase transition the quantum Fisher information (QFI), which gives the upper bound on the achievable precision of any measurement of the system and environment, becomes quadratic in observation time (cf. "Heisenberg scaling"). In fact, the QFI is identical to the variance of the dynamical observable that characterises the phases that coexist at the transition, and enhanced scaling is a consequence of the divergence of the variance of this observable at the transition point. This identification allows to establish the finite time scaling of the QFI. Near the transition the QFI is quadratic in time for times shorter than the correlation time of the dynamics. In the regime of enhanced scaling the optimal measurement whose precision is given by the QFI involves measuring both system and output. As a particular realisation of these ideas, we describe a theoretical scheme for quantum enhanced phase estimation using the photons being emitted from a quantum system near the coexistence of dynamical phases with distinct photon emission rates.

preprint2014arXiv

Emergent rhombus tilings from molecular interactions with $M$-fold rotational symmetry

We show that model molecules with particular rotational symmetries can self-assemble into network structures equivalent to rhombus tilings. This assembly happens in an emergent way, in the sense that molecules spontaneously select irregular 4-fold local coordination from a larger set of possible local binding geometries. The existence of such networks can be rationalized by simple geometrical arguments, but the same arguments do not guarantee networks' spontaneous self-assembly. This class of structures must in certain regimes of parameter space be able to reconfigure into networks equivalent to triangular tilings.

preprint2014arXiv

Excitations are localized and relaxation is hierarchical in glass-forming liquids

For several atomistic models of glass formers, at conditions below their glassy dynamics onset temperatures, ${T_\mathrm{o}}$, we use importance sampling of trajectory space to study the structure, statistics and dynamics of excitations responsible for structural relaxation. Excitations are detected in terms of persistent particle displacements of length $a$. At supercooled conditions, for $a$ of the order of or smaller than a particle diameter, we find that excitations are associated with correlated particle motions that are sparse and localized, occupying a volume with an average radius that is temperature independent and no larger than a few particle diameters. We show that the statistics and dynamics of these excitations are facilitated and hierarchical. Excitation energy scales grow logarithmically with $a$. Excitations at one point in space facilitate the birth and death of excitations at neighboring locations, and space-time excitation structures are microcosms of heterogeneous dynamics at larger scales. This nature of dynamics becomes increasingly dominant as temperature $T$ is lowered. We show that slowing of dynamics upon decreasing temperature below $T_\mathrm{o}$ is the result of a decreasing concentration of excitations and concomitant growing hierarchical length scales, and further that the structural relaxation time $τ$ follows the parabolic law, $\log(τ/ τ_\mathrm{o}) = J^2(1/T - 1/T_\mathrm{o})^2$, for $T<T_\mathrm{o}$, where $J$, $τ_\mathrm{o}$ and $T_\mathrm{o}$ can be predicted quantitatively from dynamics at short time scales. Particle motion is facilitated and directional, and we show this becomes more apparent with decreasing $T$. We show that stringlike motion is a natural consequence of facilitated, hierarchical dynamics.

preprint2014arXiv

Localisation in space and time in disordered-lattice open quantum dynamics

We study a two-dimensional tight-binding lattice for excitons with on-site disorder, coupled to a thermal environment at infinite temperature. The disorder acts to localise an exciton spatially, while the environment generates dynamics which enable exploration of the lattice. Although the steady state of the system is trivially uniform, we observe a rich dynamics and uncover a dynamical phase transition in the space of temporal trajectories. This transition is identified as a localisation in the dynamics generated by the bath. We explore spatial features in the dynamics and employ a generalisation of the inverse participation ratio to deduce an ergodic timescale for the lattice.

preprint2014arXiv

Meta-work and the analogous Jarzynski relation in ensembles of dynamical trajectories

Recently there has been growing interest in extending the thermodynamic method from static configurations to dynamical trajectories. In this approach, ensembles of trajectories are treated in an analogous manner to ensembles of configurations in equilibrium statistical mechanics: generating functions of dynamical observables are interpreted as partition sums, and the statistical properties of trajectory ensembles are encoded in free-energy functions that can be obtained through large-deviation methods in a suitable large time limit. This establishes what one can call a 'thermodynamics of trajectories'. In this paper we go a step further, and make a first connection to fluctuation theorems by generalising them to this dynamical context. We show that an effective 'meta-dynamics' in the space of trajectories gives rise to the celebrated Jarzynski relation connecting an appropriately defined 'meta-work' with changes in dynamical generating functions. We demonstrate the potential applicability of this result to computer simulations for two open quantum systems, a two-level system and the micromaser. We finally discuss the behavior of the Jarzynski relation across a first-order trajectory phase transition.

preprint2014arXiv

Open quantum reaction-diffusion dynamics: absorbing states and relaxation

We consider an extension of classical stochastic reaction-diffusion (RD) dynamics to open quantum systems. We study a class of models of hard core particles on a one-dimensional lattice whose dynamics is generated by a quantum master operator and where particle hopping is coherent while reactions, such as pair annihilation or pair coalescence, are dissipative. These are quantum open generalisations of the $A+A \to \varnothing$ and $A+A \to A$ classical RD models. We characterise the relaxation of the state towards the stationary regime via a decomposition of the system Hilbert space into transient and recurrent subspaces. We provide a complete classification of the structure of the recurrent subspace (and the non-equilibrium steady states) in terms of the dark states associated to the quantum master operator and its general spectral properties. We also show that, in one dimension, relaxation towards these absorbing dark states is slower than that predicted by a mean-field analysis due to fluctuation effects, in analogy with what occurs in classical RD systems. Numerical simulations of small systems suggest that the decay of the density in one dimension, in both the open quantum $A+A \to \varnothing$ and $A+A \to A$ cases, may go asymptotically as $t^{-b}$ with $1/2 < b < 1$.

preprint2014arXiv

Out-of-equilibrium evolution of kinetically constrained many-body quantum systems under purely dissipative dynamics

We explore the relaxation dynamics of quantum many-body systems that undergo purely dissipative dynamics through non-classical jump operators that can establish quantum coherence. Our goal is to shed light on the differences in the relaxation dynamics that arise in comparison to systems evolving via classical rate equations. In particular, we focus on a scenario where both quantum and classical dissipative evolution lead to a stationary state with the same values of diagonal or "classical" observables. As a basis for illustrating our ideas we use spin systems whose dynamics becomes correlated and complex due to dynamical constraints, inspired by kinetically constrained models (KCMs) of classical glasses. We show that in the quantum case the relaxation can be orders of magnitude slower than the classical one due to the presence of quantum coherences. Aspects of these idealized quantum KCMs become manifest in a strongly interacting Rydberg gas under electromagnetically induced transparency (EIT) conditions in an appropriate limit. Beyond revealing a link between this Rydberg gas and the rather abstract dissipative KCMs of quantum glassy systems, our study sheds light on the limitations of the use of classical rate equations for capturing the non-equilibrium behavior of this many-body system.

preprint2014arXiv

Out-of-equilibrium structures in strongly interacting Rydberg gases with dissipation

The non-equilibrium dynamics of a gas of cold atoms in which Rydberg states are off-resonantly excited is studied in the presence of noise. The interplay between interaction and off-resonant excitation leads to an initial dynamics where aggregates of excited Rydberg atoms slowly nucleate and grow, eventually reaching long-lived meta-stable arrangements which then relax further on much longer timescales. This growth dynamics is governed by an effective Master equation which permits a transparent and largely analytical understanding of the underlying physics. By means of extensive numerical simulations we study the many-body dynamics and the correlations of the resulting non-equilibrium states in various dimensions. Our results provide insight into the dynamical richness of strongly interacting Rydberg gases in noisy environments, and highlight the usefulness of these kind of systems for the exploration of soft-matter-type collective behaviour.

preprint2014arXiv

Using the s-Ensemble to Probe Glasses Formed by Cooling and Aging

From length scale distributions characterizing frozen amorphous domains, we relate the s-ensemble method with standard cooling and aging protocols for forming glass. We show that in a general class of models, where space-time scaling is in harmony with that of experiment, the domain size distributions obtained with the s-ensemble are identical to those obtained through cooling or aging, but the computational effort for applying the s-ensemble is generally many orders of magnitude smaller than that of straightforward numerical simulation of cooling or aging.

preprint2013arXiv

Calorimetric glass transition explained by hierarchical dynamic facilitation

The glass transition refers to the non-equilibrium process by which an equilibrium liquid is transformed to a non-equilibrium disordered solid, or vice versa. Associated response functions, such as heat capacities, are markedly different on cooling than on heating, and the response to melting a glass depends markedly on the cooling protocol by which the glass was formed. This paper shows how this irreversible behavior can be interpreted quantitatively in terms of an East-model picture of localized excitations (or soft spots) in which molecules can move with a specific direction, and from which excitations with the same directionality of motion can appear or disappear in adjacent regions. As a result of this facilitated dynamics, excitations become correlated in a hierarchical fashion. These correlations are manifested in the dynamic heterogeneity of the supercooled liquid phase. While equilibrium thermodynamics is virtually featureless, a non-equilibrium glass phase emerges when the model is driven out of equilibrium with a finite cooling rate. The correlation length of this emergent phase is large and increases with decreasing cooling rate. A spatially and temporally resolved fictive temperature encodes memory of its preparation. Parameters characterizing the model can be determined from reversible transport data, and with these parameters, predictions of the model agree well with irreversible differential scanning calorimetry.

preprint2013arXiv

Comment on "Is there a breakdown of the Stokes-Einstein relation in Kinetically Constrained Models at low temperature?" by O. Blondel and C. Toninelli, arXiv:1307.1651

We consider the diffusion constant, D, of a probe particle coupled to the East model, extending previous numerical results for this model to encompass a total of twelve orders of magnitude in relaxation time, τ. Our considerations thus cover the full range of relaxation times relevant to experimental studies of reversible glass formers. Throughout this range, we show that the results of the model are consistent with a fractional Stokes-Einstein law, D \propto τ^{-ξ} with ξ \approx 0.77. The same data can also be fit by relating Dτ to the East model excitation concentration q, specifically, Dτ \sim 1/q^α with α \approx 1.6. In either case, breakdown of the Stokes-Einstein relation in this kinetically constrained model is unquestionable and large.

preprint2013arXiv

Dynamical phase transitions, time-integrated observables and geometry of states

We show that there exist dynamical phase transitions (DPTs), as defined in [Phys. Rev. Lett. 110 135704 (2013)], in the transverse-field Ising model (TFIM) away from the static quantum critical points. We study a class of special states associated with singularities in the generating functions of time-integrated observables found in [Phys. Rev. B 88 184303 (2013)]. Studying the dynamics of these special states under the evolution of the TFIM Hamiltonian, we find temporal non-analtyicities in the initial-state return probability associated with dynamical phase transitions. By calculating the Berry phase and Chern number we show the set of special states have interesting geometric features similar to those associated with static quantum critical points.

preprint2013arXiv

Fluctuating observation time ensembles in the thermodynamics of trajectories

The dynamics of stochastic systems, both classical and quantum, can be studied by analysing the statistical properties of dynamical trajectories. The properties of ensembles of such trajectories for long, but fixed, times are described by large-deviation (LD) rate functions. These LD functions play the role of dynamical free-energies: they are cumulant generating functions for time-integrated observables, and their analytic structure encodes dynamical phase behaviour. This "thermodynamics of trajectories" approach is to trajectories and dynamics what the equilibrium ensemble method of statistical mechanics is to configurations and statics. Here we show that, just like in the static case, there is a variety of alternative ensembles of trajectories, each defined by their global constraints, with that of trajectories of fixed total time being just one of these. We show that an ensemble of trajectories where some time-extensive quantity is constant (and large) but where total observation time fluctuates, is equivalent to the fixed-time ensemble, and the LD functions that describe one ensemble can be obtained from those that describe the other. We discuss how the equivalence between generalised ensembles can be exploited in path sampling schemes for generating rare dynamical trajectories.

preprint2013arXiv

Kinetic constraints, hierarchical relaxation and onset of glassiness in strongly interacting and dissipative Rydberg gases

We show that the dynamics of a laser driven Rydberg gas in the limit of strong dephasing is described by a master equation with manifest kinetic constraints. The equilibrium state of the system is uncorrelated but the constraints in the dynamics lead to spatially correlated collective relaxation reminiscent of glasses. We study and quantify the evolution towards equilibrium in one and two dimensions, and analyze how the degree of glassiness and the relaxation time are controlled by the interaction strength between Rydberg atoms. We also find that spontaneous decay of Rydberg excitations leads to an interruption of glassy relaxation that takes the system to a highly correlated non-equilibrium stationary state. The results presented here, which are in principle also applicable other systems such as polar molecules and atoms with large magnetic dipole moments, show that the collective behavior of cold atomic and molecular ensembles can be similar to that found in soft condensed-matter systems.

preprint2013arXiv

Perspective: The Glass Transition

We provide here a brief perspective on the glass transition field. It is an assessment, written from the point of view of theory, of where the field is and where it seems to be heading. We first give an overview of the main phenomenological characteristics, or "stylised facts", of the glass transition problem, i.e. the central observations that a theory of the physics of glass formation should aim to explain in a unified manner. We describe recent developments, with a particular focus on real space properties, including dynamical heterogeneity and facilitation, the search for underlying spatial or structural correlations, and the relation between the thermal glass transition and athermal jamming. We then discuss briefly how competing theories of the glass transition have adapted and evolved to account for such real space issues. We consider in detail two conceptual and methodological approaches put forward recently, that aim to access the fundamental critical phenomenon underlying the glass transition, be it thermodynamic or dynamic in origin, by means of biasing of ensembles, of configurations in the thermodynamic case, or of trajectories in the dynamic case. We end with a short outlook.

preprint2013arXiv

Rare-event trajectory ensemble analysis reveals metastable dynamical phases in lattice proteins

We explore the dynamical large-deviations of a lattice heteropolymer model of a protein by means of path sampling of trajectories. We uncover the existence of non-equilibrium dynamical phase-transitions in ensembles of trajectories between active and inactive dynamical phases, whose nature depends on properties of the interaction potential. When the full heterogeneity of interactions due to the amino-acid sequence is preserved, as in a fully interacting model or in a heterogeneous version of the Gō model where only native interactions are considered, the transition is between the equilibrium native state and a highly native but kinetically trapped state. In contrast, for the homogeneous Gō model, where there is a single native energy and the sequence plays no role, the dynamical transition is a direct consequence of the static bi-stability between unfolded and native states. In the heterogeneous case the native-active and native-inactive states, despite their static similarity, have widely varying dynamical properties, and the transition between them occurs even in lattice proteins whose sequences are designed to make them optimal folders.

preprint2013arXiv

Trajectory phase transitions and dynamical Lee-Yang zeros of the Glauber-Ising chain

We examine the generating function of the time-integrated energy for the one-dimensional Glauber-Ising model. At long times, the generating function takes on a large-deviation form and the associated cumulant generating function has singularities corresponding to continuous trajec- tory (or "space-time") phase transitions between paramagnetic trajectories and ferromagnetically or anti-ferromagnetically ordered trajectories. In the thermodynamic limit, the singularities make up a whole curve of critical points in the complex plane of the counting field. We evaluate analyti- cally the generating function by mapping the generator of the biased dynamics to a non-Hermitian Hamiltonian of an associated quantum spin chain. We relate the trajectory phase transitions to the high-order cumulants of the time-integrated energy which we use to extract the dynamical Lee-Yang zeros of the generating function. This approach offers the possibility to detect continuous trajectory phase transitions from the finite-time behaviour of measurable quantities.

preprint2013arXiv

Transition in coupled replicas may not imply a finite temperature ideal glass transition in glass forming systems

A key open question in the glass transition field is whether a finite temperature thermodynamic transition to the glass state exists or not. Recent simulations of coupled replicas in atomistic models have found signatures of a static transition as a function of replica coupling. This can be viewed as evidence of an associated thermodynamic glass transition in the uncoupled system. We demonstrate here that a different interpretation is possible. We consider the triangular plaquette model, an interacting spin system which displays (East model-like) glassy dynamics in the absence of any static transition. We show that when two replicas are coupled there is a curve of equilibrium phase transitions, between phases of small and large overlap, in the temperature-coupling plane (located on the self-dual line of an exact temperature/coupling duality of the system) which ends at a critical point. Crucially, in the limit of vanishing coupling the finite temperature transition disappears, and the uncoupled system is in the disordered phase at all temperatures. We discuss interpretation of atomistic simulations in light of this result.

preprint2012arXiv

Broken symmetry and the variation of critical properties in the phase behaviour of supramolecular rhombus tilings

The degree of randomness, or partial order, present in two-dimensional supramolecular arrays of isophthalate tetracarboxylic acids is shown to vary due to subtle chemical changes such as the choice of solvent or small differences in molecular dimensions. This variation may be quantified using an order parameter and reveals a novel phase behaviour including random tiling with varying critical properties as well as ordered phases dominated by either parallel or non-parallel alignment of neighbouring molecules, consistent with long-standing theoretical studies. The balance between order and randomness is driven by small differences in the intermolecular interaction energies, which we show, using numerical simulations, can be related to the measured order parameter. Significant variations occur even when the energy difference is much less than the thermal energy highlighting the delicate balance between entropic and energetic effects in complex self-assembly processes.

preprint2012arXiv

Characterization of dynamical phase transitions in quantum jump trajectories beyond the properties of the stationary state

We describe how to characterize dynamical phase transitions in open quantum systems from a purely dynamical perspective, namely, through the statistical behavior of quantum jump trajectories. This approach goes beyond considering only properties of the steady state. While in small quantum systems dynamical transitions can only occur trivially at limiting values of the controlling parameters, in many-body systems they arise as collective phenomena and within this perspective they are reminiscent of thermodynamic phase transitions. We illustrate this in open models of increasing complexity: a three-level system, a dissipative version of the quantum Ising model, and the micromaser. In these examples dynamical transitions are accompanied by clear changes in static behavior. This is however not always the case, and in general dynamical phase behavior needs to be uncovered by observables which are strictly dynamical, e.g. dynamical counting fields. We demonstrate this via the example of a class of models of dissipative quantum glasses, whose dynamics can vary widely despite having identical (and trivial) stationary states.

preprint2012arXiv

Dynamical phases and intermittency of the dissipative quantum Ising model

We employ the concept of a dynamical, activity order parameter to study the Ising model in a transverse magnetic field coupled to a Markovian bath. For a certain range of values of the spin-spin coupling, magnetic field and dissipation rate, we identify a first order dynamical phase transition between active and inactive {\em dynamical phases}. We demonstrate that dynamical phase-coexistence becomes manifest in an intermittent behavior of the bath quanta emission. Moreover, we establish the connection between the dynamical order parameter that quantifies the activity, and the longitudinal magnetization that serves as static order parameter. The system we consider can be implemented in current experiments with Rydberg atoms and trapped ions.

preprint2012arXiv

Dynamical Phases in the Full Counting Statistics of the Resonant-Level Model

We present a thermodynamic formalism to study the full counting statistics (FCS) of charge transport through a quantum dot coupled to two leads in the resonant-level model. We show that a close analogue of equilibrium phase transitions exists for the statistics of transferred charge; by tuning an appropriate `counting field', crossovers to different dynamical phases are possible. Our description reveals a mapping between the FCS of a given device and current measurements over a range of devices with different dot-lead coupling strengths. Further insight into features in the FCS is found by studying the occupation of the dot conditioned on the transported charge between the leads.

preprint2012arXiv

Facilitated spin models of dissipative quantum glasses

We introduce a class of dissipative quantum spin models with local interactions and without quenched disorder that show glassy behaviour. These models are the quantum analogs of the classical facilitated spin models. Just like their classical counterparts, quantum facilitated models display complex glassy dynamics despite the fact that their stationary state is essentially trivial. In these systems, dynamical arrest is a consequence of kinetic constraints and not of static ordering. These models display a quantum version of dynamic heterogeneity: the dynamics towards relaxation is spatially correlated despite the absence of static correlations. Associated dynamical fluctuation phenomena such as decoupling of timescales is also observed. Moreover, we find that close to the classical limit quantum fluctuations can enhance glassiness, as recently reported for quantum liquids.

preprint2012arXiv

Inactive dynamical phase of a symmetric exclusion process on a ring

We investigate the nature of the dynamically inactive phase of a simple symmetric exclusion process on a ring. We find that as the system's activity is tuned to a lower-than-average value the particles progressively lump into a single cluster, thereby forming a kink in the density profile. All dynamical regimes, and their finite size range of validity, are explicitly determined.

preprint2012arXiv

Phase Transitions in Trajectories of a Superconducting Single-Electron Transistor Coupled to a Resonator

Recent progress in the study of dynamical phase transitions has been made with a large-deviation approach to study trajectories of stochastic jumps using a thermodynamic formalism. We study this method applied to an open quantum system consisting of a superconducting single-electron transistor, near the Josephson quasiparticle resonance, coupled to a resonator. We find that the dynamical behavior shown in rare trajectories can be rich even when the mean dynamical activity is small and thus the formalism gives insights into the form of fluctuations. The structure of the dynamical phase diagram found from the quantum-jump trajectories of the resonator is studied and we see that sharp transitions in the dynamical activity may be related to the appearance and disappearance of bistabilities in the state of the resonator as system parameters are changed. We also demonstrate that for a fast resonator, the trajectories of quasiparticles are similar to the resonator trajectories.

preprint2012arXiv

Random and ordered phases of off-lattice rhombus tiles

We study the covering of the plane by non-overlapping rhombus tiles, a problem well-studied only in the limiting case of dimer coverings of regular lattices. We go beyond this limit by allowing tiles to take any position and orientation on the plane, to be of irregular shape, and to possess different types of attractive interactions. Using extensive numerical simulations we show that at large tile densities there is a phase transition from a fluid of rhombus tiles to a solid packing with broken rotational symmetry. We observe self-assembly of broken-symmetry phases, even at low densities, in the presence of attractive tile-tile interactions. Depending on tile shape and interactions the solid phase can be random, possessing critical orientational fluctuations, or crystalline. Our results suggest strategies for controlling tiling order in experiments involving `molecular rhombi'.

preprint2012arXiv

Thermodynamics of Quadrature Trajectories in Open Quantum Systems

We apply a large-deviation method to study the diffusive trajectories of the quadrature operators of light within a reservoir connected to dissipative quantum systems. We formulate the study of quadrature trajectories in terms of characteristic operators and show that in the long time limit the statistics of such trajectories obey a large-deviation principle. We take our motivation from homodyne detection schemes which allow the statistics of quadrature operator of the light field to be measured. We illustrate our approach with four examples of increasing complexity: a driven two-level system, a `blinking' three-level system, a pair of weakly-coupled two-level driven systems, and the micromaser. We discuss how quadrature operators can serve as alternative order parameters for the classification of dynamical phases, which is particularly useful in cases where the statistics of quantum jumps cannot distinguish between such phases. The formalism we introduce also allows us to analyse the properties of the light emitted by quantum jump trajectories which fluctuate far from the typical dynamics.

preprint2012arXiv

Time-integrated observables as order parameters for dynamical phase transitions in closed quantum systems

The dynamical behaviour of many-body systems is often richer than what can be anticipated from their static properties. Here we show that in closed quantum systems this becomes evident by considering time-integrated observables as order parameters. In particular, the analytic properties of their generating functions, as estimated by full-counting statistics, allow to identify dynamical phases, i.e. phases with specific fluctuation properties of time-integrated observables, and to locate the transitions between these phases. We discuss in detail the case of the quantum Ising chain in a transverse field. We show that this model displays a continuum of quantum dynamical transitions, of which the static transition is just an end point. These singularities are not a consequence of particular choices of initial conditions or other external non-equilibrium protocols such as quenches in coupling constants. They can be probed generically through quantum jump statistics of an associated open problem, and for the case of the quantum Ising chain we outline a possible experimental realisation of this scheme by digital quantum simulation with cold ions.

preprint2012arXiv

Trajectory phase transitions, Lee-Yang zeros, and high-order cumulants in full counting statistics

We investigate Lee-Yang zeros of generating functions of dynamical observables and establish a general relation between phase transitions in ensembles of trajectories of stochastic many-body systems and the time evolution of high-order cumulants of such observables. This connects dynamical free-energies for full counting statistics in the long-time limit, which can be obtained via large-deviation methods and whose singularities indicate dynamical phase transitions, to observables that are directly accessible in simulation and experiment. As an illustration we consider facilitated spin models of glasses and show that from the short-time behavior of high-order cumulants it is possible to infer the existence and location of dynamical or "space-time" transitions in these systems.

preprint2011arXiv

Preparation and relaxation of very stable glassy states of a simulated liquid

We prepare metastable glassy states in a model glass-former made of Lennard-Jones particles by sampling biased ensembles of trajectories with low dynamical activity. These trajectories form an inactive dynamical phase whose `fast' vibrational degrees of freedom are maintained at thermal equilibrium by contact with a heat bath, while the `slow' structural degrees of freedom are located in deep valleys of the energy landscape. We examine the relaxation to equilibrium and the vibrational properties of these metastable states. The glassy states we prepare by our trajectory sampling method are very stable to thermal fluctuations and also more mechanically rigid than low-temperature equilibrated configurations.

preprint2011arXiv

Quantum trajectory phase transitions in the micromaser

We study the dynamics of the single atom maser, or micromaser, by means of the recently introduced method of thermodynamics of quantum jump trajectories. We find that the dynamics of the micromaser displays multiple space-time phase transitions, i.e., phase transitions in ensembles of quantum jump trajectories. This rich dynamical phase structure becomes apparent when trajectories are classified by dynamical observables that quantify dynamical activity, such as the number of atoms that have changed state while traversing the cavity. The space-time transitions can be either first-order or continuous, and are controlled not just by standard parameters of the micromaser but also by non-equilibrium "counting" fields. We discuss how the dynamical phase behavior relates to the better known stationary state properties of the micromaser.

preprint2010arXiv

Corresponding States of Structural Glass Formers. II

The earlier paper of this same title demonstrated a collapse of relaxation data of fragile supercooled glass forming liquids [\textit{J. Phys. Chem. B} \textbf{113}, 5563-5567 (2009)]. For temperature $T$ below that of the onset to supercooled behavior, $T_{\mathrm{o}}$, the logarithm of structural relaxation time, $\log τ$, is given by the parabolic form $\log (τ/τ_{\mathrm{o}}) \, =\, J^2(1/T - 1/T_{\mathrm{o}})^{2}$, where $J$ and $τ_{\mathrm{o}}$ are temperature independent. This paper presents further applications of this formula. In particular, it is shown that the effects of attractive forces in numerical simulation of glass forming liquids can be logically organized in terms of $J$ and $T_{\mathrm{o}}$. Further, analysis of experimental data for several systems suggests that $J$ and $T_{\mathrm{o}}$ are material properties. In contrast, values of similar parameters for other fitting formulas are shown to depend not only upon the material but also upon the range of data used in fitting these formulas. Expressions demonstrated to fail in this way include the Vogel-Fulcher-Tammann formula, a double-exponential formula, and a fractional exponential formula.

preprint2010arXiv

Entropically-stabilised growth of a two-dimensional random tiling

The assembly of molecular networks into structures such as random tilings and glasses has recently been demonstrated for a number of two-dimensional systems. These structures are dynamically-arrested on experimental timescales so the critical regime in their formation is that of initial growth. Here we identify a transition from energetic to entropic stabilisation in the nucleation and growth of a molecular rhombus tiling. Calculations based on a lattice gas model show that clustering of topological defects and the formation of faceted boundaries followed by a slow relaxation to equilibrium occurs under conditions of energetic stabilisation. We also identify an entropically-stabilised regime in which the system grows directly into an equilibrium configuration without the need for further relaxation. Our results provide a methodology for identifying equilibrium and non-equilibrium randomness in the growth of molecular tilings, and we demonstrate that equilibrium spatial statistics are compatible with exponentially slow dynamical behaviour.

preprint2010arXiv

Finite-temperature critical point of a glass transition

We generalize the simplest kinetically constrained model of a glass-forming liquid by softening kinetic constraints, allowing them to be violated with a small finite rate. We demonstrate that this model supports a first-order dynamical (space-time) phase transition, similar to those observed with hard constraints. In addition, we find that the first-order phase boundary in this softened model ends in a finite-temperature dynamical critical point, which we expect to be present in natural systems. We discuss links between this critical point and quantum phase transitions, showing that dynamical phase transitions in $d$ dimensions map to quantum transitions in the same dimension, and hence to classical thermodynamic phase transitions in $d+1$ dimensions. We make these links explicit through exact mappings between master operators, transfer matrices, and Hamiltonians for quantum spin chains.

preprint2010arXiv

Kinetically Constrained Models

In this chapter we summarize recent developments in the study of kinetically constrained models (KCMs) as models for glass formers. After recalling the definition of the KCMs which we cover we study the possible occurrence of ergodicity breaking transitions and discuss in some detail how, before any such transition occurs, relaxation timescales depend on the relevant control parameter (density or temperature). Then we turn to the main issue: the prediction of KCMs for dynamical heterogeneities. We focus in particular on multipoint correlation functions and susceptibilities, and decoupling in the transport coefficients. Finally we discuss the recent view of KCMs as being at first order coexistence between an active and an inactive space-time phase.

preprint2010arXiv

Space-time Phase Transitions in Driven Kinetically Constrained Lattice Models

Kinetically constrained models (KCMs) have been used to study and understand the origin of glassy dynamics. Despite having trivial thermodynamic properties, their dynamics slows down dramatically at low temperatures while displaying dynamical heterogeneity as seen in glass forming supercooled liquids. This dynamics has its origin in an ergodic-nonergodic first-order phase transition between phases of distinct dynamical "activity". This is a "space-time" transition as it corresponds to a singular change in ensembles of trajectories of the dynamics rather than ensembles of configurations. Here we extend these ideas to driven glassy systems by considering KCMs driven into non-equilibrium steady states through non-conservative forces. By classifying trajectories through their entropy production we prove that driven KCMs also display an analogous first-order space-time transition between dynamical phases of finite and vanishing entropy production. We also discuss how trajectories with rare values of entropy production can be realized as typical trajectories of a mapped system with modified forces.

preprint2010arXiv

Thermalization in a coherently driven ensemble of two-level systems

We study the coherent quantum evolution of a closed and driven mesoscopic chain of two-level systems that interact via the van-der-Waals interaction in their excited state. The Hamiltonian consists of a part corresponding to a classical lattice gas and an off-diagonal driving term without classical counterpart. We show that in a certain parameter range the latter leads to a thermalization of the system with respect to observables of the classical lattice gas such as the interaction energy and particle number distribution. We investigate the evolution of the system into this thermal state and discuss how to determine the corresponding temperature. Our findings can be applied to understand thermalization in strongly interacting systems of laser-driven Rydberg atoms, ions or polar molecules.

preprint2009arXiv

Dynamics on the Way to Forming Glass: Bubbles in Space-time

We review a theoretical perspective of the dynamics of glass forming liquids and the glass transition. It is a perspective we have developed with our collaborators during this decade. It is based upon the structure of trajectory space. This structure emerges from spatial correlations of dynamics that appear in disordered systems as they approach non-ergodic or jammed states. It is characterized in terms of dynamical heterogeneity, facilitation and excitation lines. These features are associated with a newly discovered class of non-equilibrium phase transitions. Equilibrium properties have little if anything to do with it. The broken symmetries of these transitions are obscure or absent in spatial structures, but they are vivid in space-time (i.e., trajectory space). In our view, the glass transition is an example of this class of transitions. The basic ideas and principles we review were originally developed through the analysis of idealized and abstract models. Nevertheless, the central ideas are easily illustrated with reference to molecular dynamics of more realistic atomistic models, and we use that illustrative approach here.

preprint2009arXiv

Metastable states and space-time phase transitions in a spin-glass model

We study large deviations of the dynamical activity in the random orthogonal model (ROM). This is a fully connected spin-glass model with one-step replica symmetry breaking behaviour, consistent with the random first-order transition scenario for structural glasses. We show that this model displays dynamical (space-time) phase-transitions between active and inactive phases, as demonstrated by singularities in large deviation functions. We argue that such transitions are generic in systems with long-lived metastable states.

preprint2009arXiv

Thermodynamics of Quantum Jump Trajectories

We apply the large-deviation method to study trajectories in dissipative quantum systems. We show that in the long time limit the statistics of quantum jumps can be understood from thermodynamic arguments by exploiting the analogy between large-deviation and free-energy functions. This approach is particularly useful for uncovering properties of rare dissipative trajectories. We also prove, via an explicit quantum mapping, that rare trajectories of one system can be realized as typical trajectories of an alternative system.

preprint2008arXiv

Corresponding States of Structural Glass Formers

The variation with respect to temperature T of transport properties of 58 fragile structural glass forming liquids (68 data sets in total) are analyzed and shown to exhibit a remarkable degree of universality. In particular, super-Arrhenius behaviors of all super-cooled liquids appear to collapse to one parabola for which there is no singular behavior at any finite temperature. This behavior is bounded by an onset temperature To above which liquid transport has a much weaker temperature dependence. A similar collapse is also demonstrated, over the smaller available range, for existing numerical simulation data.

preprint2008arXiv

Relationship between vibrations and dynamical heterogeneity in a model glass former: extended soft modes but local relaxation

We study the relation between short-time vibrational modes and long-time relaxational dynamics in a kinetically constrained lattice gas with harmonic interactions between neighbouring particles. We find a correlation between the location of the low (high) frequency vibrational modes and regions of high (low) propensity for motion. This is similar to what was observed in continuous force systems, but our interpretation is different: in our case relaxation is due to localised excitations which propagate through the system; these localised excitations act as background disorder for the elastic network, giving rise to anomalous vibrational modes. Our results show that a correlation between spatially extended low frequency modes and high propensity regions does not imply that relaxational dynamics originates in extended soft modes. We consider other measures of elastic heterogeneity, such as non-affine displacement fields and mode localisation lengths, and discuss implications of our results to interpretations of dynamic heterogeneity more generally.

preprint2007arXiv

De-coupling of Exchange and Persistence Times in Atomistic Models of Glass Formers

With molecular dynamics simulations of a fluid mixture of classical particles interacting with pair-wise additive Weeks-Chandler-Andersen potentials, we consider the time series of particle displacements and thereby determine distributions for local persistence times and local exchange times. These basic characterizations of glassy dynamics are studied over a range of super-cooled conditions and shown to have behaviors, most notably de-coupling, similar to those found in kinetically constrained lattice models of structural glasses. Implications are noted.

preprint1998arXiv

Instantons and the QCD Vacuum Wavefunctional

We analyze the instanton transitions in the framework of the gauge invariant variational calculation in the pure Yang-Mills theory. Instantons are identified with the saddle points in the integration over the gauge group which projects the Gaussian wave functional onto the gauge invariant physical Hilbert space. We show that the dynamical mass present in the best variational state provides an infrared cutoff for the instanton sizes. The instantons of the size $ρ<1/M$ are suppressed and the large size instanton problem arising in the standard WKB calculation is completely avoided in the present variational framework.