Fast kinetic simulator for relativistic matter.

Relativistic kinetic theory is ubiquitous to several fields of modern physics, finding application at large scales in systems in astrophysical contexts, all of the way down to subnuclear scales and into the realm of quark-gluon plasmas. This motivates the quest for powerful and efficient computational methods that are able to accurately study fluid dynamics in the relativistic regime as well as the transition to beyond hydrodynamics-in principle all of the way down to ballistic regimes. We present a family of relativistic lattice kinetic schemes for the efficient simulation of relativistic flows in both strongly (fluid) and weakly (rarefied gas) interacting regimes.

Scaling Law Describes the Spin-Glass Response in Theory, Experiments, and Simulations.

The correlation length ξ, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length.

Probing bulk viscosity in relativistic flows.

We derive an analytical connection between kinetic relaxation rate and bulk viscosity of a relativistic fluid in spatial dimensions, all the way from the ultra-relativistic down to the near non-relativistic regime. Our derivation is based on both Chapman-Enskog asymptotic expansion and Grad's method of moments. We validate our theoretical results against a benchmark flow, providing further evidence of the correctness of the Chapman-Enskog approach; we define the range of validity of this approach and provide evidence of mounting departures at increasing Knudsen number.

The Mpemba effect in spin glasses is a persistent memory effect.

The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom-built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a nonequilibrium process, governed by the coherence length ξ of the system. The effect occurs when the bath temperature lies in the glassy phase, but it is not necessary for the thermal protocol to cross the critical temperature.

Hybrid Monte Carlo algorithm for sampling rare events in space-time histories of stochastic fields.

We introduce a variant of the Hybrid Monte Carlo (HMC) algorithm to address large-deviation statistics in stochastic hydrodynamics. Based on the path-integral approach to stochastic (partial) differential equations, our HMC algorithm samples space-time histories of the dynamical degrees of freedom under the influence of random noise. First, we validate and benchmark the HMC algorithm by reproducing multiscale properties of the one-dimensional Burgers equation driven by Gaussian and white-in-time noise.

Relativistic dissipation obeys Chapman-Enskog asymptotics: Analytical and numerical evidence as a basis for accurate kinetic simulations.

We present an analytical derivation of the transport coefficients of a relativistic gas in (2+1) dimensions for both Chapman-Enskog (CE) asymptotics and Grad's expansion methods. We further develop a systematic calibration method, connecting the relaxation time of relativistic kinetic theory to the transport parameters of the associated dissipative hydrodynamic equations. Comparison of our analytical results and numerical simulations shows that the CE method correctly captures dissipative effects, while Grad's method does not, in agreement with previous analyses performed in the (3+1)-dimensional case.

Prospects for the Detection of Electronic Preturbulence in Graphene.

Based on extensive numerical simulations, accounting for electrostatic interactions and dissipative electron-phonon scattering, we propose experimentally realizable geometries capable of sustaining electronic preturbulence in graphene samples. In particular, preturbulence is predicted to occur at experimentally attainable values of the Reynolds number between 10 and 50, over a broad spectrum of frequencies between 10 and 100 GHz.

Aging Rate of Spin Glasses from Simulations Matches Experiments.

Experiments on spin glasses can now make precise measurements of the exponent z(T) governing the growth of glassy domains, while our computational capabilities allow us to make quantitative predictions for experimental scales. However, experimental and numerical values for z(T) have differed. We use new simulations on the Janus II computer to resolve this discrepancy, finding a time-dependent z(T,t_{w}), which leads to the experimental value through mild extrapolations.

Kinetic approach to relativistic dissipation.

Despite a long record of intense effort, the basic mechanisms by which dissipation emerges from the microscopic dynamics of a relativistic fluid still elude complete understanding. In particular, several details must still be finalized in the pathway from kinetic theory to hydrodynamics mainly in the derivation of the values of the transport coefficients. In this paper, we approach the problem by matching data from lattice-kinetic simulations with analytical predictions.

Towards a unified lattice kinetic scheme for relativistic hydrodynamics.

We present a systematic derivation of relativistic lattice kinetic equations for finite-mass particles, reaching close to the zero-mass ultrarelativistic regime treated in the previous literature. Starting from an expansion of the Maxwell-Jüttner distribution on orthogonal polynomials, we perform a Gauss-type quadrature procedure and discretize the relativistic Boltzmann equation on space-filling Cartesian lattices. The model is validated through numerical comparison with standard tests and solvers in relativistic fluid dynamics such as Boltzmann approach multiparton scattering and previous relativistic lattice Boltzmann models.

Matching Microscopic and Macroscopic Responses in Glasses.

We first reproduce on the Janus and Janus II computers a milestone experiment that measures the spin-glass coherence length through the lowering of free-energy barriers induced by the Zeeman effect. Secondly, we determine the scaling behavior that allows a quantitative analysis of a new experiment reported in the companion Letter [S. Guchhait and R.

A statics-dynamics equivalence through the fluctuation-dissipation ratio provides a window into the spin-glass phase from nonequilibrium measurements.

We have performed a very accurate computation of the nonequilibrium fluctuation-dissipation ratio for the 3D Edwards-Anderson Ising spin glass, by means of large-scale simulations on the special-purpose computers Janus and Janus II. This ratio (computed for finite times on very large, effectively infinite, systems) is compared with the equilibrium probability distribution of the spin overlap for finite sizes. Our main result is a quantitative statics-dynamics dictionary, which could allow the experimental exploration of important features of the spin-glass phase without requiring uncontrollable extrapolations to infinite times or system sizes.

Evolution of a double-front Rayleigh-Taylor system using a graphics-processing-unit-based high-resolution thermal lattice-Boltzmann model.

We study the turbulent evolution originated from a system subjected to a Rayleigh-Taylor instability with a double density at high resolution in a two-dimensional geometry using a highly optimized thermal lattice-Boltzmann code for GPUs. Our investigation's initial condition, given by the superposition of three layers with three different densities, leads to the development of two Rayleigh-Taylor fronts that expand upward and downward and collide in the middle of the cell. By using high-resolution numerical data we highlight the effects induced by the collision of the two turbulent fronts in the long-time asymptotic regime.

Dynamical transition in the D=3 Edwards-Anderson spin glass in an external magnetic field.

We study the off-equilibrium dynamics of the three-dimensional Ising spin glass in the presence of an external magnetic field. We have performed simulations both at fixed temperature and with an annealing protocol. Thanks to the Janus special-purpose computer, based on field-programmable gate array (FPGAs), we have been able to reach times equivalent to 0.

Thermodynamic glass transition in a spin glass without time-reversal symmetry.

Spin glasses are a longstanding model for the sluggish dynamics that appear at the glass transition. However, spin glasses differ from structural glasses in a crucial feature: they enjoy a time reversal symmetry. This symmetry can be broken by applying an external magnetic field, but embarrassingly little is known about the critical behavior of a spin glass in a field.

Second-order closure in stratified turbulence: simulations and modeling of bulk and entrainment regions.

The parametrization of small-scale turbulent fluctuations in convective systems and in the presence of strong stratification is a key issue for many applied problems in oceanography, atmospheric science, and planetology. In the presence of stratification, one needs to cope with bulk turbulent fluctuations and with inversion regions, where temperature, density, or both develop highly nonlinear mean profiles due to the interactions between the turbulent boundary layer and the unmixed-stable-flow above or below it. We present a second-order closure able to cope simultaneously with both bulk and boundary layer regions, and we test it against high-resolution state-of-the-art two-dimensional numerical simulations in a convective and stratified belt for values of the Rayleigh number up to Ra∼10(10).

Numerical simulations of Rayleigh-Taylor front evolution in turbulent stratified fluids.

We present state-of-the-art numerical simulations of a two-dimensional Rayleigh-Taylor instability for a compressible stratified fluid. We describe the computational algorithm and its implementation on the QPACE supercomputer. High resolution enables the statistical properties of the evolving interface that we characterize in terms of its fractal dimension to be studied.

Static versus dynamic heterogeneities in the D = 3 Edwards-Anderson-Ising spin glass.

We numerically study the aging properties of the dynamical heterogeneities in the Ising spin glass. We find that a phase transition takes place during the aging process. Statics-dynamics correspondence implies that systems of finite size in equilibrium have static heterogeneities that obey finite-size scaling, thus signaling an analogous phase transition in the thermodynamical limit.

Nonequilibrium spin-glass dynamics from picoseconds to a tenth of a second.

We study numerically the nonequilibrium dynamics of the Ising spin glass, for a time spanning 11 orders of magnitude, thus approaching the experimentally relevant scale (i.e., seconds).

Evidences of Bolgiano-Obhukhov scaling in three-dimensional Rayleigh-Bénard convection.

We present different results from high-resolution high-statistics direct numerical simulations of a three-dimensional convective cell. We test the fundamental physical picture of the presence of both a Bolgiano-Obhukhov-like and a Kolmogorov-like regime. We find that the dimensional predictions for these two distinct regimes (characterized, respectively, by an active and passive role of the temperature field) are consistent with our analysis.