Transport properties of stochastic fluids

Basic Information

• Paper ID: 2510.12557
• Title: Transport properties of stochastic fluids
• Authors: Chandrodoy Chattopadhyay, Josh Ott, Thomas Schäfer, Vladimir V. Skokov
• Classification: nucl-th cond-mat.quant-gas hep-ph
• Publication Date: October 14, 2025
• Paper Link: https://arxiv.org/abs/2510.12557

Abstract

This paper investigates thermal conductivity and momentum transport properties in stochastic hydrodynamics. The authors study fluids described by Model H in the Hohenberg-Halperin classification, examining transport properties of critical and non-critical fluids in two and three dimensions. The research is based on numerical simulations of Model H using the Metropolis algorithm, with transport coefficients extracted via Kubo relations. The expected logarithmic divergence of shear viscosity in two-dimensional non-critical fluids is observed. At the critical point, transport coefficients exhibit power-law scaling with system size L. The thermal conductivity κ shows the strongest divergence in two dimensions, κ ∼ L^{x_κ}, where x_κ = 1.6 ± 0.1. Weaker divergence is observed in three dimensions, x_κ = 1.25 ± 0.3, with the shear viscosity scaling exponent x_η significantly smaller than x_κ in both dimensions.

Research Background and Motivation

Problem Definition

The core problems addressed in this research are understanding transport properties in stochastic hydrodynamics, particularly behavior near the critical point. Specifically:

1. Behavior of thermal conductivity and momentum transport under stochastic fluctuations
2. Scaling behavior of transport coefficients in critical fluids
3. Differences in transport properties between two- and three-dimensional systems

Significance

The importance of stochastic effects in hydrodynamics is manifested in several aspects:

1. Effective Field Theory Foundation: Fluctuations and fluctuation-dissipation relations are core elements for formulating hydrodynamics as an effective field theory
2. Discovery of New Phenomena: Fluctuations lead to new phenomena, including non-analytic long-time tail effects that are missed in naive gradient expansion
3. Phase Transition Description: Fluctuation dynamics are crucial for describing fluids near phase transitions, particularly critical slowing down, dynamic scaling, and critical transport behavior
4. QCD Phase Diagram Applications: Model H is expected to describe dynamics near a possible endpoint in the QCD phase diagram

Limitations of Existing Methods

Traditional research approaches have the following limitations:

1. ε-expansion Methods: Historically, critical dynamics have been studied primarily using ε-expansion, but this involves assumptions and approximations
2. Mode Coupling Theory: Based on non-systematic approximations
3. Functional Renormalization Group: While successful, difficult to generalize to non-trivial backgrounds
4. Numerical Simulation Constraints: Previous studies were largely limited to linearized theory

Core Contributions

1. Verification of Theoretical Predictions: Successfully verified logarithmic divergence of shear viscosity in two-dimensional non-critical fluids
2. Critical Scaling Measurement: First precise measurement through numerical simulation of power-law scaling exponents of transport coefficients near the critical point
3. Kubo Relation Application: Successful application of Kubo relations within the stochastic hydrodynamics framework to extract transport coefficients
4. Multi-dimensional Comparison: Systematic comparison of transport properties differences between two- and three-dimensional systems
5. Model Validation: Provided numerical validation and benchmarking for Model H and Model H0

Methodology Details

Task Definition

The research task is to determine transport coefficients in stochastic fluids through numerical simulation, including:

• Input: Stochastic partial differential equations of Model H
• Output: Numerical values of shear viscosity η and thermal conductivity κ and their scaling behavior
• Constraints: Satisfaction of fluctuation-dissipation relations and conservation laws

Model Architecture

Model H Equations

Model H describes the coupling of a scalar order parameter field φ with fluid momentum density π⃗ in conservative form:

∂_t φ = -∇⃗ · j⃗                    (1)
∂_t π^T_i = -P^T_ij ∇_k Π_jk       (2)

where P^T_ij = δ_ij - ∇_i∇_j/∇² is the transverse projection operator.

Flux Expressions

j_i = (1/ρ)φπ^T_i - κ∇_i(δH/δφ) + θ_i                    (3)
Π_ij = (1/ρ)π^T_i π^T_j + ∇_iφ∇_jφ - (η/ρ)(∇_iπ^T_j + ∇_jπ^T_i) + Λ_ij    (4)

Effective Hamiltonian

H = ∫d^d x [1/(2ρ)(π^T_i)² + 1/2(∇φ)² + 1/2 m²φ² + 1/4 λφ⁴]    (5)

Noise Correlations

⟨θ_i(t,x⃗)θ_j(t',x⃗')⟩ = 2κT δ_ij δ(x⃗-x⃗')δ(t-t')    (6)
⟨Λ_ij(t,x⃗)Λ_kl(t',x⃗')⟩ = 2ηT(δ_ikδ_jl + δ_ilδ_jk - 2/3 δ_ijδ_kl)δ(x⃗-x⃗')δ(t-t')    (7)

Technical Innovations

Kubo Relation Application

Transport coefficients are extracted using Kubo relations:

Shear Viscosity:

η = 1/(TV) ∫₀^∞ dt ⟨T_xy(0)T_xy(t)⟩    (20)

Thermal Conductivity:

κ = 1/(3TV) ∫₀^∞ dt ⟨J_i(0)J^i(t)⟩    (23)

Metropolis Algorithm

The Metropolis algorithm is employed to ensure:

1. Automatic satisfaction of fluctuation-dissipation relations
2. Combination of diffusion and stochastic updates into a single Metropolis step
3. Updates of approximately conserved deterministic advection terms

Experimental Setup

Numerical Parameters

• Lattice Spacing: a
• System Sizes: L = 16, 24, 36, 40, 48, 64
• Time Step: Δt = 0.04 t_lat
• Temperature Units: T₀
• Typical Parameters: η₀ = 1, ρ = 1 (lattice units)

Evaluation Metrics

1. Scaling Exponents: x_κ, x_η (power-law exponents of critical transport coefficients)
2. Logarithmic Divergence Coefficients: Verification of theoretical prediction proportionality constants
3. Scaling Relations: Verification of x_κ + x_η = 4 - d - η* and z = 4 - x_κ - η*

Comparison Methods

1. Theoretical Predictions: ε-expansion, mode coupling theory
2. Different Models: Model H vs Model H₀ (without self-advection terms)
3. Different Dimensions: Two-dimensional vs three-dimensional systems

Experimental Results

Main Results

Two-Dimensional Non-Critical Fluids

• Logarithmic Divergence Verification: Expected δη ∝ log(L) behavior observed
• Fitting Coefficient: 0.0158 (theoretical prediction: 1/(16π) ≈ 0.0198)
• Relative Error: ~20%, consistent with lattice effects

Critical Scaling Behavior

Two-Dimensional Systems:

• Model H: x_κ = 1.603 ± 0.1, x_η = 0.247 ± 0.1
• Model H₀: x_κ = 1.70 ± 0.1, x_η = 0.404 ± 0.1

Three-Dimensional Systems:

• Model H: x_κ = 1.266 ± 0.3
• Model H₀: x_κ = 1.545 ± 0.3

Scaling Relation Verification

Two-Dimensional Results

Model H: x_κ + x_η = 1.85, theoretical prediction: 2 - η* = 1.75 Dynamic scaling: z = 2.11, theoretical prediction: 4 - x_κ - η* = 2.15

Three-Dimensional Results

Model H: z = 3.02, theoretical prediction: 4 - x_κ - η* ≈ 2.7 Some discrepancy exists, possibly related to systematic errors

Key Findings

1. Stronger Thermal Conductivity Divergence: Critical divergence of κ is significantly stronger than η, consistent with ε-expansion predictions
2. Dimensional Effects: Divergence in two dimensions is stronger than in three dimensions
3. Model Consistency: Models H and H₀ yield consistent critical exponents within error margins
4. Kubo Relation Validity: Successfully applied to extract IR-sensitive transport properties

Related Work

Theoretical Foundation

1. Hohenberg-Halperin Classification: Standard framework for dynamic critical phenomena
2. ε-Expansion Theory: Renormalization group treatment by Siggia, Halperin, and Hohenberg
3. Mode Coupling Theory: Kawasaki approximation and related extensions

Development of Numerical Methods

1. Early Work: Numerical methods for stochastic Navier-Stokes equations by Bell et al.
2. Lattice Methods: Spectral function lattice calculations by Berges et al.
3. Metropolis Algorithm: Series of works by the Schäfer-Skokov group

Application Background

1. QCD Phase Diagram: Search for critical points in heavy-ion collisions
2. Cold Atom Systems: Transport phenomena near critical points
3. Statistical Physics: Universality classes of critical dynamics

Conclusions and Discussion

Main Conclusions

1. Successful Numerical Verification: Kubo relations effectively extract transport coefficients of stochastic fluids
2. Confirmation of Critical Scaling: Clear power-law scaling behavior observed, κ ∼ L^{x_κ}, η ∼ L^{x_η}
3. Theoretical Consistency: Results are basically consistent with ε-expansion and scaling relations
4. Method Reliability: Provides benchmarks for studying critical transport in complex backgrounds

Limitations

1. UV Sensitivity: Kubo formula cannot accurately extract viscosity corrections related to ultraviolet cutoff
2. Systematic Errors: Uncertainties from finite-size effects and lattice discretization
3. Computational Complexity: Large dynamic exponent z ≳ 3 in three-dimensional systems leads to high computational demands
4. Precision Limitations: Difficulty in extracting x_η in three dimensions

Future Directions

1. Discretization Improvements: Seek discretization schemes that automatically satisfy Kubo relations
2. Precision Enhancement: Improve numerical precision of x_κ and x_η, particularly measurement of x_η in three dimensions
3. Finite-Size Corrections: Better understanding of finite-size scaling corrections to Kubo integrals
4. Practical Applications: Study the role of critical transport in real systems (e.g., quark-gluon plasma)

In-Depth Evaluation

Strengths

1. Methodological Innovation: First systematic application of Kubo relations to stochastic hydrodynamics
2. Theoretical Verification: Successfully verified multiple important theoretical predictions (logarithmic divergence, power-law scaling)
3. Systematic Study: Comprehensive comparison across multiple dimensions and models
4. Numerical Accuracy: Achieves relatively high numerical precision within technical constraints
5. Physical Insight: Reveals stronger critical divergence of thermal conductivity compared to viscosity

Weaknesses

1. UV Issues: Systematic problems in lattice discretization when handling UV-sensitive quantities
2. Error Control: Relatively large systematic errors in some results, particularly in three dimensions
3. Computational Cost: High dynamic exponent limits system sizes due to computational complexity
4. Theoretical Discrepancies: Notable deviations in verification of certain scaling relations

Impact

1. Academic Value: Provides important numerical benchmarks for critical dynamics research
2. Methodological Contribution: Application of Kubo relations in stochastic fluids has methodological significance
3. Cross-Disciplinary Applications: Results have value for QCD, cold atoms, statistical physics, and other fields
4. Technical Advancement: Promotes development of numerical solution techniques for stochastic partial differential equations

Applicable Scenarios

1. Basic Research: Theoretical verification and prediction of critical point dynamics
2. Heavy-Ion Physics: Theoretical support for QCD critical point search
3. Cold Atom Experiments: Theoretical guidance for critical behavior of ultracold Fermi gases
4. Computational Physics: Methodological reference for stochastic fluid numerical simulations

References

The paper cites 62 important references, spanning from classical Martin-Siggia-Rose theory to the latest functional renormalization group methods, providing a solid theoretical foundation and comprehensive literature background for the research.

This paper makes important contributions to the fields of stochastic hydrodynamics and critical phenomena research, particularly in numerical methods and theoretical verification. Despite some technical limitations, it establishes an important foundation for future related research.