@article{chattopadhyay_ott_schaefer_v. skokov_2024, title={Simulations of Stochastic Fluid Dynamics near a Critical Point in the Phase Diagram}, volume={133}, ISSN={["1079-7114"]}, url={https://doi.org/10.1103/PhysRevLett.133.032301}, DOI={10.1103/PhysRevLett.133.032301}, abstractNote={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 cross-over 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 time-dependent correlation function of non-Gaussian moments M n (t) of the order parameter.We find that the relaxation time depends in non-trivial manner on the power n.}, number={3}, journal={PHYSICAL REVIEW LETTERS}, author={Chattopadhyay, Chandrodoy and Ott, Josh and Schaefer, Thomas and V. Skokov, Vladimir}, year={2024}, month={Jul} }
 @article{chattopadhyay_ott_schafer_skokov_2023, title={Dynamic scaling of order parameter fluctuations in model B}, volume={108}, ISSN={["2470-0029"]}, url={https://doi.org/10.1103/PhysRevD.108.074004}, DOI={10.1103/PhysRevD.108.074004}, abstractNote={We describe numerical simulations of the stochastic diffusion equation with a conserved charge. 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\simeq 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.}, number={7}, journal={PHYSICAL REVIEW D}, author={Chattopadhyay, Chandrodoy and Ott, Josh and Schafer, Thomas and Skokov, Vladimir}, year={2023}, month={Oct} }