Nonequilibrium Dynamics of Dirac Quantum Criticality in Imaginary Time

Through large-scale quantum Monte Carlo simulations and a generalized nonequilibrium scaling theory, this study reveals unconventional negative critical exponents in the imaginary-time relaxation of chiral Heisenberg Dirac quantum criticality, establishing an efficient short-time dynamics framework for investigating fermionic quantum criticality.

The Gist

Imagine you are watching a pot of water on a stove. If you heat it slowly, it eventually boils. That boiling point is a "phase transition." In the world of quantum physics, particles can also undergo sudden, dramatic changes in their behavior, called Quantum Phase Transitions.

This paper is about studying a very specific, exotic type of quantum transition involving Dirac fermions (particles that act like massless light, found in materials like graphene). The researchers wanted to understand how these particles behave right at the moment of this transition, but with a twist: they looked at it in Imaginary Time.

Here is the breakdown of their discovery using simple analogies:

1. The Setting: A Quantum Dance Floor

Think of the material (a honeycomb lattice of atoms) as a giant dance floor.

• The Dancers: Electrons (fermions) are the dancers.
• The Music: The "interaction" between them is the music.
• The Transition: At a specific volume of music (interaction strength), the dancers suddenly switch from a chaotic, free-flowing jam session (a "Dirac Semimetal") to a rigid, synchronized marching band (an "Antiferromagnet").

2. The Trick: "Imaginary Time"

In physics, we usually study how things change in real time (seconds ticking by). But to simulate these quantum systems on a computer, scientists use a mathematical trick called Imaginary Time.

• The Analogy: Imagine you are trying to find the lowest point in a foggy valley (the ground state). In real time, you might wander around aimlessly. In "Imaginary Time," it's like gravity suddenly gets super strong, and you slide down the hill instantly to the bottom.
• The Problem: Usually, to get to the bottom, you have to slide for a very long time. If the hill is very flat near the bottom (which happens at a critical point), it takes forever. This is called "critical slowing down."

3. The Discovery: The "Initial Slip"

The researchers realized they didn't need to wait for the system to settle down. They could learn everything they needed by watching the first few seconds of the slide.

They started the simulation with three different "starting positions" for the dancers:

1. The Marching Band: Everyone is already in perfect order.
2. The Jam Session: Everyone is dancing freely.
3. The Random Crowd: Everyone is dancing randomly, with no pattern.

The Big Surprise:
When they started with the Random Crowd, they observed a weird phenomenon called the "Initial Slip."

• In Normal Physics (Classical): If you start with a random crowd, they usually start to organize themselves a little bit before the chaos takes over. It's like a "slip" where order briefly increases. This is described by a positive number.
• In This Quantum World: The researchers found the opposite! The "slip" went backward. The order didn't just fail to increase; it actively decreased faster than expected. They measured a negative number (specifically -0.84) to describe this.

Why?
The paper explains that in these Dirac systems, the "fermion dancers" (the electrons) are incredibly fast and sensitive. They react to the chaos almost instantly, preventing the "order" (the magnetic alignment) from ever getting a foothold. It's like trying to build a sandcastle while a tsunami is hitting it; the water (fermion fluctuations) destroys the castle before you can even lay the first brick.

4. The New Toolkit: A Shortcut for Scientists

The most important part of this paper isn't just the weird negative number; it's the new method they created.

• Old Way: To find the rules of a quantum system, you usually have to simulate it until it reaches a perfect, calm state. This takes massive computing power and time, and often fails because of a "sign problem" (a mathematical glitch that makes calculations impossible for certain materials).
• New Way (This Paper): They proved you can figure out the rules of the system by looking only at the very beginning of the process (the short-time dynamics).
  • Analogy: Instead of waiting for a cake to bake for an hour to see if it's good, you can tell if the recipe is right just by smelling the batter in the first 10 seconds.

Why Does This Matter?

1. Speed: It allows scientists to study complex quantum materials much faster.
2. Solving the "Sign Problem": Because the simulation stops early (before the math gets too messy), they can study materials that were previously impossible to calculate.
3. Future Tech: This helps us understand materials like graphene and topological insulators, which are the building blocks for future super-fast computers and quantum devices.

In a Nutshell:
The authors found that when quantum particles undergo a dramatic change, they behave in a counter-intuitive way right at the start, actively resisting order due to their unique speed. By harnessing this "initial slip," they created a fast, efficient shortcut to understand the deepest secrets of quantum matter without waiting for the system to settle down.

Technical Summary

1. Problem Statement

The paper addresses the gap in understanding nonequilibrium dynamics within Dirac quantum criticality, specifically focusing on the chiral Heisenberg universality class.

• Context: While equilibrium properties of Dirac fermions (e.g., in graphene and topological materials) and their quantum phase transitions (QPTs) are well-studied, the short-time relaxation dynamics in imaginary time (ITCD) for systems with coupled fermionic and bosonic critical modes remains largely unexplored.
• Challenge: Traditional equilibrium methods (like standard Quantum Monte Carlo) suffer from critical slowing down and the fermion sign problem as systems approach the ground state, making the study of critical properties computationally expensive and often intractable for large systems.
• Specific Question: How do gapless Dirac fermion fluctuations influence the short-time relaxation dynamics in imaginary time, and can a unified scaling theory be established to extract critical exponents efficiently without reaching equilibrium?

2. Methodology

The authors employ large-scale, sign-problem-free Determinant Quantum Monte Carlo (DQMC) simulations to study the imaginary-time relaxation dynamics of the Honeycomb Hubbard Model.

• Model: The system is defined by the Hamiltonian:
  $$H = -t \sum_{\langle ij \rangle, \sigma} c^\dagger_{i\sigma} c_{j\sigma} + U \sum_i \left(n_{i\uparrow} - \frac{1}{2}\right)\left(n_{i\downarrow} - \frac{1}{2}\right)$$
  The system undergoes a phase transition from a Dirac Semimetal (DSM) to an Antiferromagnetic (AFM) insulator at a critical interaction $U_c/t \approx 3.9$.
• Dynamical Protocol:
  • The system evolves according to the imaginary-time Schrödinger equation: $|\psi(\tau)\rangle = e^{-\tau H} |\psi(0)\rangle / Z(\tau)$.
  • Three distinct initial states are prepared to probe different fixed points under renormalization group (RG) transformations:
    1. Saturated AFM state: Ordered magnetic state.
    2. Non-interacting DSM state: Gapless Dirac semimetal.
    3. Random Spin (RS) state: Uncorrelated state with random spin orientations.
• Observables:
  • Order Parameter ($m^2$): Square of the staggered magnetization (AFM structure factor).
  • Fermion Correlation ($G$): Correlation function at the Dirac point momentum, related to the quasi-particle weight $Z$.
  • Auto-correlation ($A$): Used to characterize the initial slip behavior.
• Scaling Analysis: The authors generalize the short-time scaling theory to include both bosonic (order parameter) and fermionic (correlation) modes, analyzing data in the short-time window $1 \ll \tau \ll L^z$.

3. Key Contributions

1. Generalized Nonequilibrium Scaling Theory: The authors extend the standard short-time scaling theory (previously applied mainly to bosonic systems) to Dirac systems. They derive specific scaling forms for observables ($m^2$ and $G$) that account for the distinct relaxation behaviors of fermionic and bosonic critical modes.
2. Discovery of Negative Critical Initial Slip: They identify a negative critical initial slip exponent ($\theta \approx -0.84$) in the relaxation of the auto-correlation function starting from a random spin state. This contrasts sharply with classical and bosonic quantum systems where $\theta$ is typically positive.
3. Efficient Critical Property Determination: They demonstrate that critical exponents ($\nu, \beta, \eta_\psi, z$) and the critical point ($U_c$) can be accurately determined using only short-time imaginary-time data, bypassing the need for long-time equilibrium evolution.
4. Sign Problem Mitigation Strategy: The work proposes a framework where the severity of the fermion sign problem (which grows exponentially with imaginary time) is circumvented by extracting critical information from the short-time regime, potentially enabling the study of sign-problematic fermionic models.

4. Key Results

• Critical Point and Exponents:

  • Determined the critical interaction ratio: $U_c/t = 3.91(3)$.
  • Correlation length exponent: $\nu = 1.17(7)$.
  • Order parameter scaling dimension: $\beta/\nu = 0.80(3)$.
  • Fermion anomalous dimension: $\eta_\psi = 0.15(4)$.
  • Dynamic exponent: $z = 1$ (consistent with Lorentz symmetry).
  • Note: These values align with previous equilibrium studies but were obtained with significantly less computational effort.

• Negative Initial Slip Exponent ($\theta$):

  • From the Random Spin (RS) initial state, the auto-correlation function $A(\tau)$ scales as $A(\tau) \sim \tau^{-d/z + \theta}$.
  • The extracted exponent is $\theta = -0.84(4)$.
  • Physical Mechanism: In classical systems, domains grow initially ($\theta > 0$). In Dirac criticality, gapless fermions relax much faster than bosonic modes. These rapid fermionic fluctuations suppress the formation of magnetic domains in the early stage, leading to a decay of correlations and a negative $\theta$. This confirms that fermionic critical fluctuations dominate the early-time dynamics.

• Scaling Collapse:

  • Data for $m^2$ and $G$ from different system sizes ($L$) and initial states collapse onto universal curves when rescaled according to the derived scaling forms (e.g., $m^2 \sim \tau^{-2\beta/\nu z}$ for AFM initial state; $m^2 \sim L^{-d}\tau^{d/z - 2\beta/\nu z}$ for DSM/RS states).
  • The scaling forms successfully distinguish between the relaxation behaviors of the DSM and AFM initial states.

5. Significance

• Theoretical Advancement: The paper provides the first unified theoretical framework for nonequilibrium imaginary-time dynamics in Dirac systems, explicitly handling the interplay between fermionic and bosonic critical modes.
• Computational Efficiency: By utilizing short-time dynamics, the method offers a highly efficient route to determine critical exponents, avoiding the computational bottleneck of critical slowing down.
• Overcoming the Sign Problem: The approach suggests a viable path to study fermionic quantum critical points in models where the sign problem is severe, as the simulation can be terminated before the sign problem becomes prohibitive (short $\tau$).
• Experimental Relevance: With the advent of quantum simulators (e.g., Rydberg atoms, quantum circuits) capable of simulating fermionic systems and realizing imaginary-time evolution, this work provides a blueprint for experimentally detecting quantum criticality and measuring critical exponents with high efficiency and scalability.

In summary, this work establishes a powerful new paradigm for investigating quantum criticality in strongly correlated fermionic systems, revealing unique non-stationary dynamics driven by gapless fermions and offering a practical solution to long-standing computational challenges.