Engineering Papers⌕ Search

Engineering topics

Schäfer, Thomas

Publications and source records attributed to Schäfer, Thomas.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations of Stochastic Fluid Dynamics near a Critical Point in the Phase Diagram

Here, we present simulations of stochastic fluid dynamics in the vicinity of a critical endpoint belonging to the universality class of the Ising model. This study is motivated by the challenge of modeling the dynamics of critical fluctuations near a conjectured critical endpoint in the phase diagram of quantum chromodynamics (QCD). We focus on the interaction of shear modes with a conserved scalar density, which is known as model H. We show that the observed dynamical scaling behavior depends on the correlation length and the shear viscosity of the fluid. As the correlation length is increased or the viscosity is decreased we observe a crossover from the dynamical exponent of critical diffusion, z≃4, to the expected scaling exponent of model H, z≃3. We use our method to investigate the time-dependent correlation function of non-Gaussian moments M n (t) of the order parameter. We find that the relaxation time depends in a nontrivial manner on the power n.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Magnetism and metallicity in moiré transition metal dichalcogenides

The ability to control the properties of twisted bilayer transition metal dichalcogenides in situ makes them an ideal platform for investigating the interplay of strong correlations and geometric frustration. Of particular interest are the low energy scales, which make it possible to experimentally access both temperature and magnetic fields that are of the order of the bandwidth or the correlation scale. In this manuscript, we analyze the moiré Hubbard model, believed to describe the low energy physics of an important subclass of the twisted bilayer compounds. We establish its magnetic and the metal–insulator phase diagram for the full range of magnetic fields up to the fully spin-polarized state. We find a rich phase diagram including fully and partially polarized insulating and metallic phases of which we determine the interplay of magnetic order, Zeeman-field, and metallicity, and make connection to recent experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamics of non-Gaussian fluctuations in model A

Motivated by the experimental search for the QCD critical point, we perform simulations of a stochastic field theory with purely relaxational dynamics (model A). We verify the expected dynamic scaling of correlation functions. Using a finite size scaling analysis, we obtain the dynamic critical exponent z = 2.026(56). We investigate time dependent correlation functions of higher moments M n (t) of the order parameter M(t) for n = 1, 2, 3, 4. We obtain dynamic scaling with the same critical exponent z for all n, but the relaxation constant depends on n. We also study the relaxation of M n (t) after a quench, where the simulation is initialized in the high temperature phase, and the dynamics is studied at the critical temperature T c . Finally, we find that the evolution does not follow simple scaling with the dynamic exponent z, and that it involves an early time rise followed by late stage relaxation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Lastly, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Dissipative superfluid hydrodynamics for the unitary Fermi gas

In this work we establish constraints on the temperature dependence of the shear viscosity η in the superfluid phase of a dilute Fermi gas in the unitary limit. Our results are based on analyzing experiments that measure the aspect ratio of a deformed cloud after release from an optical trap. Here, we discuss how to apply the two-fluid formalism to the unitary gas and provide a suitable parametrization of the equation of state. We show that in expansion experiments the difference between the normal and superfluid velocities remains small and can be treated as a perturbation. We find that expansion experiments favor a shear viscosity that decreases significantly in the superfluid regime. Using an exponential parametrization, we find η(T c /2T F ) ≲ 0.37[η(T c /T F )], where T c is the critical temperature and T F is the local Fermi temperature of the gas.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multiplicative noise and the diffusion of conserved densities

Stochastic fluid dynamics governs the long time tails of hydrodynamic correlation functions, and the critical slowing down of relaxation phenomena in the vicinity of a critical point in the phase diagram. In this work we study the role of multiplicative noise in stochastic fluid dynamics. Multiplicative noise arises from the dependence of transport coefficients, such as the diffusion constants for charge and momentum, on fluctuating hydrodynamic variables. We study long time tails and relaxation in the diffusion of a conserved density (model B), and a conserved density coupled to the transverse momentum density (model H). Careful attention is paid to fluctuation-dissipation relations. We observe that multiplicative noise contributes at the same order as non-linear interactions in model B, but is a higher order correction to the relaxation of a scalar density and the tail of the stress tensor correlation function in model H.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Diffusive dynamics of critical fluctuations near the QCD critical point

A quantitatively reliable theoretical description of the dynamics of fluctuations in nonequilibrium is indispensable in the experimental search for the QCD critical point by means of ultrarelativistic heavy-ion collisions. Herein we consider the fluctuations of the net-baryon density which becomes the slow, critical mode near the critical point. Due to net-baryon number conservation the dynamics is described by the fluid dynamical diffusion equation, which we extend to contain a white noise stochastic current. Including nonlinear couplings from the 3d Ising model universality class in the free energy functional, we solve the fully interacting theory in a finite size system. We observe that purely Gaussian white noise generates non-Gaussian fluctuations, but finite size effects and exact net-baryon number conservation lead to significant deviations from the expected behavior in equilibrated systems. In particular the skewness shows a qualitative deviation from infinite volume expectations.With this benchmark established we study the real-time dynamics of the fluctuations. We recover the expected dynamical scaling behavior and observe retardation effects and the impact of critical slowing down near the pseudocritical temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Critical points at infinity, non-Gaussian saddles, and bions

It has been argued that many non-perturbative phenomena in quantum mechanics (QM) and quantum field theory (QFT) are determined by complex field configurations, and that these contributions should be understood in terms of Picard-Lefschetz theory. In this work we compute the contribution from non-BPS multi-instanton configurations, such as instanton-anti-instanton [$I\overline{I}$] pairs, and argue that these contributions should be interpreted as exact critical points at infinity. The Lefschetz thimbles associated with such critical points have a specific structure arising from the presence of non-Gaussian, quasi-zero mode (QZM), directions. When fermion degrees of freedom are present, as in supersymmetric theories, the effective bosonic potential can be written as the sum of a classical and a quantum potential. We show that in this case the semi-classical contribution of the critical point at infinity vanishes, but there is a non-trivial contribution that arises from its associated non-Gaussian QZM-thimble. This approach resolves several puzzles in the literature concerning the semi-classical contribution of correlated [$I\overline{I}$] pairs. It has the surprising consequence that the configurations dominating the expansion of observables, and the critical points defining the Lefschetz thimble decomposition need not be the same, a feature not present in the traditional Picard-Lefschetz approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Hydrodynamic tails and a fluctuation bound on the bulk viscosity

Here, we study the small frequency behavior of the bulk viscosity spectral function using stochastic fluid dynamics. We obtain a number of model independent results, including the long-time tail of the bulk stress correlation function and the leading nonanalyticity of the spectral function at small frequency. We also establish a lower bound on the bulk viscosity which is weakly dependent on assumptions regarding the range of applicability of fluid dynamics. The bound on the bulk viscosity $ζ$ scales as $ζ_{min} ~( P - \frac{2}{3}\mathscr{E}) ^2Σ_iD_i^{-2}$, where $D_i$ are the diffusion constants for energy and momentum and $P - \frac{2}{3}\mathscr{E}$, where $P$ is the pressure and $\mathscr{E}$ is the energy density, is a measure of scale breaking. Applied to the cold Fermi gas near unitarity, $|λ/a_s| ≳ 1$, where $λ$ is the thermal de Broglie wavelength and $a_s$ is the s-wave scattering length, this bound implies that the ratio of bulk viscosity to entropy density satisfies $ζ/s ≳ 0.1ℏ/k_B$. Here, $ℏ$ is Planck's constant and $k_B$ is Boltzmann's constant.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Properties of hot and dense matter from relativistic heavy ion collisions

Here, we review the progress achieved in extracting the properties of hot and dense matter from relativistic heavy ion collisions at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory and the large hadron collider (LHC) at CERN. We focus on bulk properties of the medium, in particular the evidence for thermalization, aspects of the equation of state, transport properties, as well as fluctuations and correlations. We also discuss the in-medium properties of hadrons with light and heavy quarks, and measurements of dileptons and quarkonia. Lastly, this review is dedicated to the memory of Gerald E. Brown.

43 PARTICLE ACCELERATORS↗