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. The method can deal with both massless and massive particles, thereby encompassing ultra- and mildly relativistic regimes alike. The computational performance of the method for the simulation of relativistic flows across the aforementioned regimes is discussed in detail, along with prospects of future applications.

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. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer.

Probing bulk viscosity in relativistic flows.

We derive an analytical connection between kinetic relaxation rate and bulk viscosity of a relativistic fluid in d 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. Finally, we present numerical simulations of transport processes in quark-gluon plasmas, with special focus on the effects of bulk viscosity which might prove amenable to future experimental verification. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

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. These results provide a solid basis for accurately calibrated computational studies of relativistic dissipative flows.

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. Second, we show how to implement an importance sampling protocol to significantly enhance, by orders of magnitudes, the probability to sample extreme and rare events, making it possible to estimate moments of field variables of extremely high order (up to 30 and more). By employing reweighting techniques, we map the biased configurations back to the original probability measure in order to probe their statistical importance. Finally, we show that by biasing the system towards very intense negative gradients, the HMC algorithm is able to explore the statistical fluctuations around instanton configurations. Our results will also be interesting and relevant in lattice gauge theory since they provide unique insights into reweighting techniques.

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. Furthermore, theoretical insight is gained by studying a crossover between the T=T_{c} and T=0 fixed points.

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. Our numerical results provide neat evidence in favor of the Chapman-Enskog [The Mathematical Theory of Non-Uniform Gases, 3rd ed. (Cambridge University Press, Cambridge, U.K., 1970)] procedure as suggested by recent theoretical analyses along with qualitative hints at the basic reasons why the Chapman-Enskog expansion might be better suited than Grad's method [Commun. Pure Appl. Math. 2, 331 (1949)0010-364010.1002/cpa.3160020403] to capture the emergence of dissipative effects in relativistic fluids.

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. This work provides a significant step towards the formulation of a unified relativistic lattice kinetic scheme, covering both massive and near-massless particles regimes.

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. Orbach, Phys. Rev. Lett. 118, 157203 (2017)].PRLTAO0031-900710.1103/PhysRevLett.118.157203 The value of the coherence length estimated through the analysis of microscopic correlation functions turns out to be quantitatively consistent with its measurement through macroscopic response functions. Further, nonlinear susceptibilities, recently measured in glass-forming liquids, scale as powers of the same microscopic length.

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. We also provide details on the optimized lattice-Boltzmann code that we have run on a cluster of GPUs.

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.01 s in experiments. We have studied the system relaxation both for high and for low temperatures, clearly identifying a dynamical transition point. This dynamical temperature is strictly positive and depends on the external applied magnetic field. We discuss different possibilities for the underlying physics, which include a thermodynamical spin-glass transition, a mode-coupling crossover, or an interpretation reminiscent of the random first-order picture of structural glasses.

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). Data are taken from a Rayleigh-Taylor system confined by the existence of an adiabatic gradient.

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. We compute the critical exponents and the transition point in the equilibrium setting, and use them to show that aging in dynamic heterogeneities can be described by a finite-time scaling ansatz, with potential implications for experimental work.

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). We introduce novel analysis techniques to compute the coherence length in a model-independent way. We present strong evidence for a replicon correlator and for overlap equivalence. The emerging picture is compatible with noncoarsening behavior.