Many-body quantum electrodynamics networks: Non-equilibrium condensed matter physics with light

We review recent developments regarding the quantum dynamics and many-body physics with light, in superconducting circuits and Josephson analogues, by analogy with atomic physics. We start with quantum impurity models addressing dissipative and driven systems. Both theorists and experimentalists are making efforts towards the characterization of these non-equilibrium quantum systems. We show how Josephson junction systems can implement the equivalent of the Kondo effect with microwave photons. The Kondo effect can be characterized by a renormalized light frequency and a peak in the Rayleigh elastic transmission of a photon. We also address the physics of hybrid systems comprising mesoscopic quantum dot devices coupled with an electromagnetic resonator. Then, we discuss extensions to Quantum Electrodynamics (QED) Networks allowing one to engineer the Jaynes–Cummings lattice and Rabi lattice models through the presence of superconducting qubits in the cavities. This opens the door to novel many-body physics with light out of equilibrium, in relation with the Mott–superfluid transition observed with ultra-cold atoms in optical lattices. Then, we summarize recent theoretical predictions for realizing topological phases with light. Synthetic gauge fields and spin–orbit couplings have been successfully implemented in quantum materials and with ultra-cold atoms in optical lattices — using time-dependent Floquet perturbations periodic in time, for example — as well as in photonic lattice systems. Finally, we discuss the Josephson effect related to Bose–Hubbard models in ladder and two-dimensional geometries, producing phase coherence and Meissner currents. The Bose–Hubbard model is related to the Jaynes–Cummings lattice model in the large detuning limit between light and matter (the superconducting qubits). In the presence of synthetic gauge fields, we show that Meissner currents subsist in an insulating Mott phase.     

Non—equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method

In this paper, we propose a new approach to implementing boundary conditions in the lattice Boltzmann method (LBM). The basic idea is to decompose the distribution function at the boundary node into its equilibrium and non-equilibrium parts, and then to approximate the non-equilibrium part with a first-order extrapolation of the non-equilibrium part of the distribution at the neighbouring fluid node. Schemes for velocity and pressure boundary conditions are constructed based on this method. The resulting schemes are of second-order accuracy. Numerical tests show that the numerical solutions of the LBM together with the present boundary schemes are in excellent agreement with the analytical solutions. Second-order convergence is also verified from the results. It is also found that the numerical stability of the present schemes is much better than that of the original extrapolation schemes proposed by Chen et al. (1996 Phys. Fluids 8 2527).     

Accuracy analysis of high-order lattice Boltzmann models for rarefied gas flows

In this work, we have theoretically analyzed and numerically evaluated the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows. In the incompressible limit, the LB equation is shown to be able to reduce to the linearized Bhatnagar–Gross–Krook (BGK) equation. Therefore, when the same Gauss–Hermite quadrature is used, LB method closely resembles the discrete velocity method (DVM). In addition, the order of Hermite expansion for the equilibrium distribution function is found not to be directly correlated with the approximation order in terms of the Knudsen number to the BGK equation for incompressible flows. Meanwhile, we have numerically evaluated the LB models for a standing-shear-wave problem, which is designed specifically for assessing model accuracy by excluding the influence of gas molecule/surface interactions at wall boundaries. The numerical simulation results confirm that the high-order terms in the discrete equilibrium distribution function play a negligible role in capturing non-equilibrium effect for low-speed flows. By contrast, appropriate Gauss–Hermite quadrature has the most significant effect on whether LB models can describe the essential flow physics of rarefied gas accurately. Our simulation results, where the effect of wall/gas interactions is excluded, can lead to conclusion on the LB modeling capability that the models with higher-order quadratures provide more accurate results. For the same order Gauss–Hermite quadrature, the exact abscissae will also modestly influence numerical accuracy. Using the same Gauss–Hermite quadrature, the numerical results of both LB and DVM methods are in excellent agreement for flows across a broad range of the Knudsen numbers, which confirms that the LB simulation is similar to the DVM process. Therefore, LB method can offer flexible models suitable for simulating continuum flows at the Navier–Stokes level and rarefied gas flows at the linearized Boltzmann model equation level.     

Glassy dynamics of kinetically constrained models

We review the use of kinetically constrained models (KCMs) for the study of dynamics in glassy systems. The characteristic feature of KCMs is that they have trivial, often non-interacting, equilibrium behaviour but interesting slow dynamics due to restrictions on the allowed transitions between configurations. The basic question which KCMs ask is therefore how much glassy physics can be understood without an underlying 'equilibrium glass transition'. After a brief review of glassy phenomenology, we describe the main model classes, which include spin-facilitated (Ising) models, constrained lattice gases, models inspired by cellular structures such as soap froths, models obtained via mappings from interacting systems without constraints, and finally related models such as urn, oscillator, tiling and needle models. We then describe the broad range of techniques that have been applied to KCMs, including exact solutions, adiabatic approximations, projection and mode-coupling techniques, diagrammatic approaches and mappings to quantum systems or effective models. Finally, we give a survey of the known results for the dynamics of KCMs both in and out of equilibrium, including topics such as relaxation time divergences and dynamical transitions, nonlinear relaxation, ageing and effective temperatures, cooperativity and dynamical heterogeneities, and finally non-equilibrium stationary states generated by external driving. We conclude with a discussion of open questions and possibilities for future work.     

Asymmetric simple exclusion processes with complex lattice geometries: A review of models and phenomena

We summarize the findings of a large number of studies concerning the totally asymmetric simple exclusion process (TASEP) with complex lattice geometries. The TASEP has been recognized as a paradigm in modeling and analyzing non-equilibrium traffic systems. The paper surveys both the observed physical phenomena and several popular meanfield approaches used to analyze the extended TASEP models. Several interesting physical phenomena, such as phase separation, spontaneous symmetry breaking, and the finite-size effect, have been identified and explained. The future investigations of the extended TASEP with complex lattice geometries are also introduced. This paper may help to obtain a better understanding of non-equilibrium systems.     

Non-equilibrium critical phenomena and phase transitions into absorbing states

This review addresses recent developments in non-equilibrium statistical physics. Focusing on phase transitions from fluctuating phases into absorbing states, the universality class of directed percolation is investigated in detail. The survey gives a general introduction to various lattice models of directed percolation and studies their scaling properties, field-theoretic aspects, numerical techniques, as well as possible experimental realizations. In addition, several examples of absorbing-state transitions which do not belong to the directed percolation universality class will be discussed. As a closely related technique, we investigate the concept of damage spreading. It is shown that this technique is ambiguous to some extent, making it impossible to define chaotic and regular phases in stochastic non-equilibrium systems. Finally, we discuss various classes of depinning transitions in models for interface growth which are related to phase transitions into absorbing states.     

About the Definition of the Local Equilibrium Lattice Temperature in Suspended Monolayer Graphene

The definition of temperature in non-equilibrium situations is among the most controversial questions in thermodynamics and statistical physics. In this paper, by considering two numerical experiments simulating charge and phonon transport in graphene, two different definitions of local lattice temperature are investigated: one based on the properties of the phonon–phonon collision operator, and the other based on energy Lagrange multipliers. The results indicate that the first one can be interpreted as a measure of how fast the system is trying to approach the local equilibrium, while the second one as the local equilibrium lattice temperature. We also provide the explicit expression of the macroscopic entropy density for the system of phonons, by which we theoretically explain the approach of the system toward equilibrium and characterize the nature of the equilibria, in the spatially homogeneous case.     

Deterministic derivation of non-equilibrium free energy theorems for natural isothermal isobaric systems

The non-equilibrium free energy theorems show how distributions of work along non-equilibrium paths are related to free energy differences between the equilibrium states at the end points of these paths. In this paper we develop a natural way of barostatting a system and give the first deterministic derivation of the Crooks and Jarzynski relations for these isothermal isobaric systems. We illustrate these relations by applying them to molecular dynamics simulations of a model polymer undergoing stretching.     

A non-equilibrium Ising model of turbulence

We introduce a model of interacting lattices at different resolutions driven by the two-dimensional Ising dynamics with a nearest-neighbor interaction. We study this model both with tools borrowed from equilibrium statistical mechanics as well as non-equilibrium thermodynamics. Our findings show that this model keeps the signature of the equilibrium phase transition. The critical temperature of the equilibrium models corresponds to the state maximizing the entropy and delimits two out-of-equilibrium regimes, one satisfying the Onsager relations for systems close to equilibrium and one resembling convective turbulent states. Since the model preserves the entropy and energy fluxes in the scale space, it seems a good candidate for parametric studies of out-of-equilibrium turbulent systems.     

Mixtures of anisotropic and spherical colloids: Phase behavior, confinement, percolation phenomena and kinetics

Purely entropic systems such as suspensions of hard rods, platelets and spheres show rich phase behavior. Rods and platelets have successfully been used as models to predict the equilibrium properties of liquid crystals for several decades. Over the past years hard particle models have also been studied in the context of non-equilibrium statistical mechanics, in particular regarding the glass transition, jamming, sedimentation and crystallization. Recently suspensions of hard anisotropic particles also moved into the focus of materials scientists who work on conducting soft matter composites. An insulating polymer resin that is mixed with conductive filler particles becomes conductive when the filler percolates. In this context the mathematical topic of connectivity percolation finds an application in modern nano-technology. In this article, we briefly review recent work on the phase behavior, confinement effects, percolation transition and phase transition kinetics in hard particle models. In the first part, we discuss the effects that particle anisotropy and depletion have on the percolation transition. In the second part, we present results on the kinetics of the liquid-to-crystal transition in suspensions of spheres and of ellipsoids.     

Criticality in a non-equilibrium,driven system: Charged colloidal rods (fd-viruses) in electric fields

Experiments on suspensions of charged colloidal rods (fd-virus particles) in external electric fields are performed, which show that a non-equilibrium critical point can be identified. Several transition lines of field-induced phases and states meet at this point and it is shown that there is a length- and time-scale which diverge at the non-equilibrium critical point. The off-critical and critical behavior is characterized, with both power law and logarithmic divergencies. These experiments show that analogous features of the classical, critical divergence of correlation lengths and relaxation times in equilibrium systems are also exhibited by driven systems that are far out of equilibrium, related to phases/states that do not exist in the absence of the external field.     

Macroscopic Fluctuation Theory for Stationary Non-Equilibrium States

We formulate a dynamical fluctuation theory for stationary non-equilibrium states (SNS) which is tested explicitly in stochastic models of interacting particles. In our theory a crucial role is played by the time reversed dynamics. Within this theory we derive the following results: the modification of the Onsager–Machlup theory in the SNS; a general Hamilton–Jacobi equation for the macroscopic entropy; a non-equilibrium, nonlinear fluctuation dissipation relation valid for a wide class of systems; an H theorem for the entropy. We discuss in detail two models of stochastic boundary driven lattice gases: the zero range and the simple exclusion processes. In the first model the invariant measure is explicitly known and we verify the predictions of the general theory. For the one dimensional simple exclusion process, as recently shown by Derrida, Lebowitz, and Speer, it is possible to express the macroscopic entropy in terms of the solution of a nonlinear ordinary differential equation; by using the Hamilton–Jacobi equation, we obtain a logically independent derivation of this result.     

Non-Equilibrium Entropy and Irreversibility in Generalized Stochastic Loewner Evolution from an Information-Theoretic Perspective

In this study, we theoretically investigated a generalized stochastic Loewner evolution (SLE) driven by reversible Langevin dynamics in the context of non-equilibrium statistical mechanics. Using the ability of Loewner evolution, which enables encoding of non-equilibrium systems into equilibrium systems, we formulated the encoding mechanism of the SLE by Gibbs entropy-based information-theoretic approaches to discuss its advantages as a means to better describe non-equilibrium systems. After deriving entropy production and flux for the 2D trajectories of the generalized SLE curves, we reformulated the system’s entropic properties in terms of the Kullback–Leibler (KL) divergence. We demonstrate that this operation leads to alternative expressions of the Jarzynski equality and the second law of thermodynamics, which are consistent with the previously suggested theory of information thermodynamics. The irreversibility of the 2D trajectories is similarly discussed by decomposing the entropy into additive and non-additive parts. We numerically verified the non-equilibrium property of our model by simulating the long-time behavior of the entropic measure suggested by our formulation, referred to as the relative Loewner entropy.     

Flux expressions and NEMD perturbations for models of semi-flexible molecules

We derive energy and momentum flux expressions, for systems composed of a general class of semi-flexible molecules, in the Ciccotti-Ferrario-Ryckaert linear constraint formalism. According to this formalism, the whole set of Cartesian coordinates is divided into basic (independent) and secondary (dependent) subsets. It is found that energy and momentum flux vectors have a simple and general expression using both basic and secondary coordinates. In the case of non-equilibrium molecular dynamics, we give general and simple heat and shear flow algorithms, deriving the dissipative fluxes in the space of all Cartesian coordinates. In comparison with previous derivations for some models of flexible molecules, the present approach has more general applicability. Moreover, it leads to more efficient equilibrium and non-equilibrium molecular dynamics calculations of transport coefficients, and requires minimal programming effort.     

Response Theory for Equilibrium and Non-Equilibrium Statistical Mechanics: Causality and Generalized Kramers-Kronig Relations

We consider the general response theory recently proposed by Ruelle for describing the impact of small perturbations to the non-equilibrium steady states resulting from Axiom A dynamical systems. We show that the causality of the response functions entails the possibility of writing a set of Kramers-Kronig (K-K) relations for the corresponding susceptibilities at all orders of nonlinearity. Nonetheless, only a special class of directly observable susceptibilities obey K-K relations. Specific results are provided for the case of arbitrary order harmonic response, which allows for a very comprehensive K-K analysis and the establishment of sum rules connecting the asymptotic behavior of the harmonic generation susceptibility to the short-time response of the perturbed system. These results set in a more general theoretical framework previous findings obtained for optical systems and simple mechanical models, and shed light on the very general impact of considering the principle of causality for testing self-consistency: the described dispersion relations constitute unavoidable benchmarks that any experimental and model generated dataset must obey. The theory exposed in the present paper is dual to the time-dependent theory of perturbations to equilibrium states and to non-equilibrium steady states, and has in principle similar range of applicability and limitations. In order to connect the equilibrium and the non equilibrium steady state case, we show how to rewrite the classical response theory by Kubo so that response functions formally identical to those proposed by Ruelle, apart from the measure involved in the phase space integration, are obtained. These results, taking into account the chaotic hypothesis by Gallavotti and Cohen, might be relevant in several fields, including climate research. In particular, whereas the fluctuation-dissipation theorem does not work for non-equilibrium systems, because of the non-equivalence between internal and external fluctuations, K-K relations might be robust tools for the definition of a self-consistent theory of climate change.     

Reconfigurable computing for Monte Carlo simulations: Results and prospects of the Janus project

We describe Janus, a massively parallel FPGA-based computer optimized for the simulation of spin glasses, theoretical models for the behavior of glassy materials. FPGAs (as compared to GPUs or many-core processors) provide a complementary approach to massively parallel computing. In particular, our model problem is formulated in terms of binary variables, and floating-point operations can be (almost) completely avoided. The FPGA architecture allows us to run many independent threads with almost no latencies in memory access, thus updating up to 1024 spins per cycle. We describe Janus in detail and we summarize the physics results obtained in four years of operation of this machine; we discuss two types of physics applications: long simulations on very large systems (which try to mimic and provide understanding about the experimental non-equilibrium dynamics), and low-temperature equilibrium simulations using an artificial parallel tempering dynamics. The time scale of our non-equilibrium simulations spans eleven orders of magnitude (from picoseconds to a tenth of a second). On the other hand, our equilibrium simulations are unprecedented both because of the low temperatures reached and for the large systems that we have brought to equilibrium. A finite-time scaling ansatz emerges from the detailed comparison of the two sets of simulations. Janus has made it possible to perform spin-glass simulations that would take several decades on more conventional architectures. The paper ends with an assessment of the potential of possible future versions of the Janus architecture, based on state-of-the-art technology.     

Lattice Boltzmann model for combustion and detonation

In this paper we present a lattice Boltzmann model for combustion and detonation. In this model the fluid behavior is described by a finite-difference lattice Boltzmann model by Gan et al. [Physica A, 2008, 387: 1721]. The chemical reaction is described by the Lee-Tarver model [Phys. Fluids, 1980, 23: 2362]. The reaction heat is naturally coupled with the flow behavior. Due to the separation of time scales in the chemical and thermodynamic processes, a key technique for a successful simulation is to use the operator-splitting scheme. The new model is verified and validated by well-known benchmark tests. As a specific application of the new model, we studied the simple steady detonation phenomenon. To show the merit of LB model over the traditional ones, we focus on the reaction zone to study the non-equilibrium effects. It is interesting to find that, at the von Neumann peak, the system is nearly in its thermodynamic equilibrium. At the two sides of the von Neumann peak, the system deviates from its equilibrium in opposite directions. In the front of von Neumann peak, due to the strong compression from the reaction product behind the von Neumann peak, the system experiences a sudden deviation from thermodynamic equilibrium. Behind the von Neumann peak, the release of chemical energy results in thermal expansion of the matter within the reaction zone, which drives the system to deviate the thermodynamic equilibrium in the opposite direction. From the deviation from thermodynamic equilibrium, Δ _m *, defined in this paper, one can understand more on the macroscopic effects of the system due to the deviation from its thermodynamic equilibrium.