Unifying renormalized and bare viscosity in two-dimensional molecular dynamics simulations

This paper utilizes two-dimensional molecular dynamics simulations to introduce a wavenumber-dependent viscosity that unifies renormalized and bare viscosity by demonstrating how equilibrium stress correlations diverge at small wavenumbers while determining bare viscosity at large wavenumbers, thereby bridging mesoscopic fluctuations and macroscopic transport through microscopic dynamics.

The Gist

Imagine you are trying to understand how thick a fluid is—like honey versus water. In physics, this "thickness" is called viscosity.

Usually, we think of viscosity as a fixed property of a material, like a label on a bottle of syrup. But in the microscopic world of atoms and molecules, things get weird, especially in flat, two-dimensional systems (like a soap film or a very thin layer of gas).

Here is the story of what this paper discovered, told through simple analogies.

1. The Problem: The "Zoom-In" Paradox

Imagine you are looking at a crowd of people dancing in a square room.

• The Macro View (The Big Picture): If you stand far away and watch the whole room, the crowd moves smoothly. You can calculate how "thick" or resistant the crowd is to moving as a whole. This is the Renormalized Viscosity.
• The Micro View (The Close-Up): If you zoom in and watch just two people bumping into each other, you see chaotic, jerky movements. This is the Bare Viscosity.

The Catch: In 2D systems, the "Big Picture" viscosity doesn't stay constant. The bigger the room (the system size), the thicker the fluid appears to be. It's as if the crowd gets more resistant the more people you add, eventually becoming infinitely thick if the room were infinite. This is a mathematical nightmare for scientists because it means the standard rules of fluid dynamics break down.

Scientists knew there must be a fundamental, underlying "Bare Viscosity" (the true thickness of the fluid at the atomic level), but they couldn't find it. It was like trying to find the weight of a single grain of sand while standing on a beach that keeps growing larger and larger.

2. The Solution: The "Frequency Tuner"

The authors (Yokota, Itami, and Sasa) came up with a clever trick. Instead of looking at the whole room or just one grain of sand, they decided to look at the crowd through a tunable lens that changes based on the "wavelength" (or frequency) of the movement.

Think of it like an equalizer on a music system:

• Low Frequencies (Bass): These represent slow, large-scale movements of the whole crowd. In this range, the viscosity looks huge and divergent (the "Renormalized" value).
• High Frequencies (Treble): These represent fast, tiny, jittery movements of individual particles.

The team defined a new concept: Wavenumber-Dependent Viscosity. Imagine this as a dial that lets you measure the fluid's thickness at different "zoom levels" simultaneously.

3. The Discovery: Connecting the Dots

By running massive computer simulations (molecular dynamics) of particles bouncing around, they measured this viscosity dial. They found two amazing things:

1. The Low-End Match: When they tuned the dial to look at slow, large movements (low frequency), their new measurement perfectly matched the old, divergent "Renormalized Viscosity." This confirmed their math was right.
2. The High-End Goldmine: When they turned the dial all the way up to the fastest, smallest movements (high frequency), the viscosity stopped changing. It settled on a flat, stable number.

This stable number is the "Bare Viscosity."

It's like finding the true weight of the grain of sand. Once you look at the movement fast enough (high frequency), the chaotic crowd noise disappears, and you see the fundamental, unchanging property of the material.

4. Why This Matters

Before this paper, scientists had to guess where to "cut off" their calculations to find the bare viscosity. It was like trying to guess the price of a house by looking at a blurry photo; you had to guess how much detail to ignore.

This paper provides a systematic map. It shows that:

• The "messy" large-scale behavior (Renormalized) and the "clean" small-scale behavior (Bare) are actually two sides of the same coin.
• You can find the fundamental constants of a material (the Bare Viscosity and the "cutoff" scale) just by looking at how the fluid behaves at different speeds and sizes in a simulation.

The Takeaway Analogy

Imagine you are trying to measure the roughness of a mountain.

• From space (Macro), the mountain looks smooth.
• From a helicopter (Micro), it looks jagged.
• If you keep zooming in, the jaggedness seems to get worse and worse (divergence).

This paper invented a special camera that can take a picture of the mountain at every zoom level at once. They discovered that if you zoom in far enough to see the individual rocks (high frequency), the jaggedness stops increasing and settles on a specific, measurable roughness. That specific roughness is the "true" nature of the mountain, and now we have a reliable way to measure it without getting lost in the math.

This is a huge step forward because it allows scientists to predict how 2D fluids (like those in advanced electronics or biological membranes) will behave, bridging the gap between the chaotic world of atoms and the smooth world of everyday fluids.

Technical Summary

Here is a detailed technical summary of the paper "Unifying renormalized and bare viscosity in two-dimensional molecular dynamics simulations" by Yokota, Itami, and Sasa.

1. Problem Statement

In low-dimensional systems (specifically two dimensions), standard hydrodynamics fails to provide a complete description because transport coefficients, such as viscosity, exhibit anomalous behavior.

• Renormalized Viscosity ($\eta_R$): In 2D, the macroscopic viscosity measured in experiments or simulations diverges logarithmically with system size ($L$) due to long-time tails in stress correlations caused by thermal fluctuations. This makes $\eta_R$ a non-material parameter dependent on the system size.
• Bare Viscosity ($\eta_0$): Fluctuating hydrodynamics (FH) introduces a "bare" viscosity $\eta_0$ as a fundamental parameter in the stochastic equations of motion. However, $\eta_0$ is distinct from the observed $\eta_R$.
• The Gap: While theoretical frameworks (like dynamical renormalization group analysis) predict the relationship between $\eta_R$ and $\eta_0$, there has been no systematic, microscopic method to extract the intrinsic bare viscosity $\eta_0$ and the ultraviolet (UV) cutoff scale ($a_{uv}$) directly from molecular dynamics (MD) simulations without ad hoc assumptions. Previous attempts to define bare coefficients via projection operators were numerically intractable.

2. Methodology

The authors propose a novel framework to bridge the gap between microscopic dynamics and macroscopic transport by introducing a wavenumber-dependent viscosity, $\eta^*(k)$.

• System Setup:

  • 2D periodic system of $N$ particles interacting via a soft-core potential (chosen to match previous studies on viscosity renormalization).
  • Simulations performed using LAMMPS with the velocity-Verlet algorithm.
  • System sizes ranging up to $N=12,288$ ($L=128$) to capture the divergence regime.

• Key Definitions:

  • Fine-grained Shear Stress: Defined via Irving-Kirkwood theory, $\Pi_{xy}^{IK}(r, t)$.
  • Fourier Components: The authors analyze the time-averaged Fourier components of the shear stress, $\tilde{\Pi}_{xy}^{IK}(k, t)$.
  • Shear Stress Correlation Function: They define a finite-wavenvector correlation function:
    $$S(k, L, \tau) = \frac{1}{2\tau k_B T L^2} \int_0^\tau dt \int_0^\tau dt' \langle \tilde{\Pi}_{xy}^{IK}(k, t) \tilde{\Pi}_{xy}^{IK}(-k, t') \rangle_{eq}$$
  • Limiting Functions:
    • Renormalized Viscosity: $\eta_R(L) = \lim_{\tau \to \infty} S(k=0, L, \tau)$.
    • Thermodynamic Limit Correlation: $S^\infty(k) = \lim_{L \to \infty} S(k, L)$ for $k \neq 0$.
    • Wavenumber-Dependent Viscosity: $\eta^*(k) = \max_\theta S^\infty(k)$, where the maximum is taken over the angle $\theta$ of the wavevector $k$. This step is crucial because $S^\infty(k)$ is anisotropic (vanishes along axes due to momentum conservation) but has a well-defined maximum at $\theta = \pi/4$.

• Theoretical Connection:

  • Small-$k$ Limit: The behavior of $\eta^*(k)$ as $k \to 0$ is shown to correspond to the system-size dependence of $\eta_R(L)$ via the relation $L \approx 2\pi\sqrt{2}/k$.
  • Large-$k$ Limit: In the high-wavenumber regime (near the UV cutoff $k \approx 2\pi/a_{uv}$), nonlinear fluctuating hydrodynamic terms are negligible. The theory predicts that $S^\infty(k)$ approaches the bare viscosity $\eta_0$.

3. Key Contributions

1. Definition of $\eta^*(k)$: The authors define a robust, $L$-independent viscosity function $\eta^*(k)$ derived directly from microscopic stress correlations, which serves as a bridge between the divergent macroscopic viscosity and the intrinsic bare viscosity.
2. Resolution of Anisotropy: They address the singularity of shear stress correlations at $k=0$ and along coordinate axes by identifying the directional maximum ($\theta = \pi/4$) as the physically relevant quantity for viscosity.
3. Unified Framework: They establish a continuous mapping:
   • $\eta_R(L) \leftrightarrow \eta^*(k)$ (Small $k$ / Large $L$)
   • $\eta^*(k) \leftrightarrow \eta_0$ (Large $k$ / UV cutoff)
4. Extraction of Intrinsic Parameters: The method allows for the simultaneous determination of the bare viscosity $\eta_0$ and the UV cutoff length scale $a_{uv}$ from a single set of MD simulations, without needing separate non-equilibrium boundary measurements.

4. Results

• Divergence of $\eta_R(L)$: Simulations confirmed that $\eta_R(L)$ diverges logarithmically with system size $L$, consistent with theoretical predictions for 2D fluids.
• Convergence of $\eta^*(k)$:
  • For small $k$, $\eta^*(k)$ matches $\eta_R(L)$ when $k = 2\pi\sqrt{2}/L$, validating the equivalence between the thermodynamic limit of the correlation function and the system-size dependent viscosity.
  • For intermediate $k$ ($0 < k \lesssim \pi$),$\eta^*(k)$ exhibits a plateau, corresponding to the hydrodynamic regime.
  • For large $k$ ($k \gtrsim \pi$), $\eta^*(k)$ approaches a constant value.
• Quantitative Values:
  • The plateau at large $k$ yields the bare viscosity: $\eta_0 \approx 0.28$.
  • The transition to this plateau occurs at a wavenumber corresponding to a UV cutoff length scale: $a_{uv} \approx 2$ (in units of particle diameter $\sigma$).
• Consistency: These values are consistent with previous estimates obtained via non-equilibrium measurements for a similar system, validating the new equilibrium-based approach.

5. Significance

• Fundamental Physics: The paper resolves a long-standing issue in statistical mechanics by providing a microscopic definition and calculation method for "bare" transport coefficients in fluctuating hydrodynamics. It demonstrates that while macroscopic transport is anomalous in 2D, the underlying microscopic parameters are well-defined and intrinsic to the material.
• Methodological Advancement: The proposed method is computationally efficient (equilibrium MD) and avoids the numerical intractability of projection operator methods or the complexity of non-equilibrium boundary setups.
• Generalizability: The authors argue that this framework can be extended to other transport coefficients (e.g., thermal conductivity) and other systems, including active matter and electronic fluids, providing a general tool for connecting microscopic dynamics to macroscopic fluctuating hydrodynamics.
• Experimental Relevance: The findings offer a theoretical basis for interpreting experiments in 2D systems (e.g., soap films, dusty plasmas, 2D Fermi gases) where viscosity renormalization is significant, allowing researchers to distinguish between system-size effects and intrinsic material properties.