Geometric Phase Effect in Thermodynamic Properties and in the Imaginary-Time Multi-Electronic-State Path Integral Formulation

This paper demonstrates that the previously developed imaginary-time multi-electronic-state path integral (MES-PI) formulation naturally captures the geometric phase effect arising from conical intersections, and quantifies its impact on low-temperature thermodynamic properties using an ad hoc GP-excluded construction as a comparison baseline.

The Gist

The Big Idea: The "Ghost" in the Machine

Imagine you are walking through a forest. Usually, if you walk in a circle and return to your starting point, you feel exactly the same as when you left. You haven't changed.

But in the quantum world of molecules, there are special spots called Conical Intersections (CIs). Think of these as invisible "portals" or "hills" in the landscape of a molecule's energy. If a molecule walks around one of these portals in a circle, something strange happens: it comes back with a "ghostly" twist.

This twist is called the Geometric Phase (GP). It's like the molecule returns home wearing a different hat, or perhaps its shadow has flipped upside down, even though it took the exact same path. This isn't just a visual trick; it changes how the molecule vibrates, reacts, and stores energy.

The Problem:
For a long time, scientists had a standard way of simulating molecules (called the Born-Oppenheimer approximation). It was like using a map that ignored these "portals." If you used this old map, you would calculate the molecule's energy and heat properties correctly at high temperatures, but at low temperatures, the map would be wrong because it missed the "ghostly twist."

The Solution:
This paper does not introduce a new simulation method from scratch. Instead, it shines a light on a powerful tool that was already developed by the authors in 2018: the Multi-Electronic-State Path Integral (MES-PI).

Think of MES-PI not as a new invention, but as a sophisticated engine that was built to handle complex molecules. The breakthrough in this paper is the discovery that this existing engine naturally captures the ghostly twist without needing to know exactly where the portals are beforehand. It happens automatically because of how the "beads" in the chain talk to each other. The authors used this pre-existing method to prove that the "ghost" fundamentally changes the molecule's thermodynamics, a property that was previously unquantified in this context.

---

Key Concepts Explained with Analogies

1. The Ring Polymer (The Temperature Chain)

In quantum mechanics, to understand how a molecule behaves at a certain temperature, we don't look at it moving through real time. Instead, we use a mathematical concept called Imaginary Time, which is directly linked to temperature.

• The Beads: Imagine the molecule as a ring of beads (like a necklace). These beads are not snapshots of the molecule at different moments in real time; they are copies of the molecule coupled together.
• The Springs: The beads are connected by springs.
• The Temperature Connection: This ring represents the molecule's quantum "fuzziness" at a specific temperature.
  • Hot: The ring is tight and small (the molecule acts more like a classical particle).
  • Cold: The ring stretches out and becomes long and floppy (the molecule acts more like a wave).
• The Twist: If this necklace loops around a "Conical Intersection" (the portal), the whole necklace gets a twist. If you don't account for this twist, your calculation of the molecule's "heat" (thermodynamics) will be wrong. This is a static picture of quantum statistics, not a movie of the molecule moving.

2. The "Ad Hoc" Baseline (The Fake Reality)

To prove that the "ghostly twist" (Geometric Phase) actually matters, the authors created a counter-factual simulation.

• They took their advanced simulation (MES-PI) and manually forced the necklace to ignore the twist. They added a rule: "If the necklace loops around the portal, pretend it didn't."
• The Result: When they compared the "Real" simulation (with the twist) to the "Fake" simulation (without the twist), they saw a huge difference in the molecule's Heat Capacity (how much energy it takes to warm it up) at low temperatures.
• The Lesson: The "ghost" isn't just a curiosity; it fundamentally changes how the molecule stores heat.

3. The "Cusp" Problem (Why it's hard to calculate)

Here is a crucial detail about the math: Including the Geometric Phase actually makes the energy landscape smooth and easy to navigate.

• The Analogy: Imagine trying to walk up a hill. If you include the "ghostly twist" correctly, the hill is smooth. However, if you try to force the simulation to ignore the twist (or use older, simplified methods that try to add the twist back in artificially), the landscape develops a sharp, jagged spike (a "cusp") right at the center.
• The Problem: Standard walking algorithms (mathematical approximations) stumble on these sharp spikes. They take a very long time to figure out the right path because the math becomes unstable.
• The Fix: The authors found that the standard MES-PI method (which includes the twist naturally) handles the landscape smoothly without any issues. However, for the "Fake" simulations where the twist is removed, or for specific simplified methods that use a "winding number" to force the twist, the math gets jagged. To fix those specific cases, the authors developed a special trick called GPA-SP. This trick smooths out the path for the artificial methods, making them faster and more accurate. It is not needed for the standard, natural MES-PI method, which is already smooth.

4. The "Single-State" vs. "Multi-State" View

• The Simple View (SES): Sometimes, a molecule is so cold that it only cares about its lowest energy state. In this case, you can use a simpler version of the simulation. The paper shows that even in this simple case, you must account for the "twist" (Geometric Phase) to get the right answer.
• The Complex View (MES): In real, complex systems (like large molecules or materials), the molecule might jump between different energy states. The authors' main method (MES-PI) is the "Swiss Army Knife" that works for both simple and complex cases. It automatically handles the twists and turns without needing a manual map of where the portals are.

Why Does This Matter?

1. Accuracy at Low Temperatures: If you are studying super-cold molecules (like in quantum computing or ultracold chemistry), ignoring this "ghostly twist" gives you the wrong answer. This paper provides the proof and the tool to get it right.
2. No Need for a Map: In complex molecules, we often don't know exactly where these "portals" (Conical Intersections) are located. The existing MES-PI method finds the answer without needing to know the map in advance. It just lets the simulation run, and the physics does the rest.
3. Thermodynamics: Most previous studies focused on how molecules move (dynamics). This paper proves that the Geometric Phase also changes how molecules store heat (thermodynamics), which is crucial for understanding chemical reactions and material properties.

The Takeaway

Think of the Geometric Phase as a molecular signature. Just as a signature proves a document is authentic, this phase proves a molecule has interacted with the deep topology of its own energy landscape.

This paper says: "Don't ignore the signature." If you want to accurately predict how molecules behave, especially when they are cold and quiet, you must use a simulation method that respects this signature. The authors have demonstrated that their pre-existing engine (MES-PI) naturally captures this signature, ensuring our digital models of the molecular world are as real as the molecules themselves.

Technical Summary

1. Problem Statement

The Geometric Phase (GP), also known as the Berry phase, is a fundamental quantum effect arising when a molecular system traverses a closed loop encircling a Conical Intersection (CI) between potential energy surfaces (PESs). While the GP is well-known to influence dynamic and spectroscopic properties (e.g., vibronic levels, reaction rates), its impact on thermodynamic properties (such as heat capacity and internal energy) at low temperatures has been less rigorously quantified.

Standard Imaginary-Time Path Integral Molecular Dynamics (PIMD) based on the Born-Oppenheimer (BO) approximation often fails to capture the GP. This is because conventional implementations typically rely on real-valued, single-valued adiabatic electronic wavefunctions. However, these wavefunctions are inherently discontinuous (due to branch cuts) around a CI, leading to a sign change ($\pi$ phase) that standard numerical differentiation misses. Consequently, standard PIMD simulations may yield significant errors in low-temperature thermodynamic properties for systems with CIs.

2. Methodology

The multi-electronic-state path integral (MES-PI) formulation was established previously by Liu and Liu (J. Chem. Phys. 2018, 148, 102319). The present work does not introduce MES-PI; rather, it demonstrates that the previously developed imaginary-time MES-PI framework naturally captures the geometric phase, and quantifies the resulting thermodynamic consequences. The study employs the archetypal Jahn-Teller $E \otimes e$ model (a "Mexican hat" potential) as a benchmark system.

• Exact Benchmarks: The authors first solved the time-independent Schrödinger equation using a basis-set-expansion method for both:
  1. GP-Included Hamiltonian: Correctly accounts for the topological phase (using phase-factor-modified eigenvectors).
  2. GP-Excluded Hamiltonian: Uses brute-force differentiation of real-valued eigenvectors, artificially omitting the GP.
• MES-PI Formulation: They applied the MES-PI approach (Liu and Liu, 2018), which represents the partition function as a trace over nuclear and electronic degrees of freedom. Crucially, this formulation uses overlap matrices ($C_{j,j+1}$) between successive imaginary-time slices (beads) of the ring polymer.
  • Clarification: The ring polymer beads represent a discretisation of imaginary time $\tau \in [0, \beta\hbar]$, not a real-time trajectory; the formalism is one of quantum statistical mechanics rather than quantum dynamics.
• Artificial GP-Exclusion: To isolate the GP effect, the authors constructed an ad hoc GP-excluded MES-PI method. This involves calculating the winding number ($W$) of the ring polymer around the CI and multiplying the partition function by a topological phase factor, $( -1 )^W$, effectively canceling the natural GP.
• Single-Electronic-State (SES) Limits: The study also examined simplified schemes (BO-PI, BHA-PI, and SES-overlap-PI) to determine when the full multi-state treatment is necessary versus when single-state approximations suffice.

3. Key Contributions

• Natural Capture of GP in MES-PI: The paper demonstrates that the previously established MES-PI formulation (Liu and Liu, 2018) naturally captures the GP effect without requiring explicit knowledge of the CI location or complex gauge fields. The GP emerges intrinsically from the electronic trace of the product of statistically weighted overlap matrices along the ring polymer. Even if the input electronic wavefunctions are discontinuous (as provided by standard quantum chemistry software), the product of overlaps along the closed loop automatically accumulates the correct topological phase.
• Quantification of Thermodynamic Impact: The authors explicitly show that the GP significantly alters low-temperature thermodynamic properties. Specifically, the GP changes the degeneracy of the ground state (from non-degenerate to doubly degenerate), which manifests as distinct differences in heat capacity ($C_V$) at low temperatures.
• Convergence Analysis: A critical finding is the impact of the GP on numerical convergence.
  • GP-Included: The system behaves smoothly at the origin, allowing for standard $O(1/N_{bead}^2)$ convergence in Trotter decomposition.
  • GP-Excluded (Artificial): Removing the GP creates a "cusp" or singularity in the action at the CI. This destroys the second-order error cancellation, leading to a significantly slower convergence rate of $O(1/N_{bead}^{0.5})$ (half an order).
• Algorithmic Improvement: The authors propose and validate the GPA-SP (Geometric Phase Absorbed into Spring Potential) algorithm. This method is introduced primarily to accelerate the artificial GP-excluded MES-PI scheme and certain SES schemes where the GP is incorporated via winding-number phase factors. By absorbing the GP into the spring potential, GPA-SP restores the favorable $O(1/N_{bead}^2)$ convergence rate for these specific cases, making them computationally efficient. (Note: The standard GP-included MES-PI already converges well without this modification).

4. Key Results

• Heat Capacity Anomalies: At low temperatures (e.g., $\beta = 8$ r.u., approx. 70 K for the benchmark system), the GP-included system exhibits a lower heat capacity compared to the GP-excluded system. The GP-excluded case shows an artificial peak in heat capacity due to the incorrect ground-state degeneracy.
• Winding Number Distribution: MES-PIMD simulations reveal that ring polymers with non-zero winding numbers (encircling the CI) have a non-negligible probability (approx. 35% for $W=1$). This confirms that topological effects are statistically significant in thermal equilibrium.
• Validity of Approximations:
  • In the Multi-Electronic-State (MES) regime, the full MES-PI is required for accuracy.
  • In the Single-Electronic-State (SES) limit (large energy gaps, low temperature), simplified schemes like SES-overlap-PI (which includes ground-state overlaps) can capture the GP.
  • However, standard BO-PI and BHA-PI (without overlap terms) fail to capture the GP unless an explicit winding-number phase factor is manually added.
• Truncation Error: The study confirms that truncating the electronic Hilbert space (keeping only the lowest few states) yields reliable results provided the thermal energy ($k_B T$) is significantly smaller than the energy gap to the first discarded state ($\Delta E$).

5. Significance

This work establishes a rigorous framework for simulating nonadiabatic thermodynamic properties in complex molecular systems where Conical Intersections are present but their topology may not be known a priori.

• General Applicability: Unlike methods requiring explicit gauge potentials or continuous complex basis sets, the MES-PI approach works with standard, discontinuous real-valued electronic wavefunctions from quantum chemistry codes.
• Thermodynamic Correction: It corrects a long-standing oversight in low-temperature thermodynamics, showing that ignoring the GP leads to qualitatively wrong predictions for heat capacity and other observables.
• Foundation for Dynamics: By providing accurate equilibrium properties and initial conditions, this framework lays the groundwork for future real-time nonadiabatic dynamics simulations (e.g., surface hopping, exact factorization) in complex systems like single-molecule magnets and ultracold molecules.
• Algorithmic Efficiency: The introduction of the GPA-SP algorithm resolves the severe convergence bottleneck associated with GP-excluded simulations (and specific SES variants), which otherwise suffer from a slow $O(1/N_{bead}^{0.5})$ convergence rate, making these rigorous calculations feasible for larger systems.

In summary, the paper proves that the topological nature of the electronic wavefunction is not just a dynamic curiosity but a critical factor in the statistical mechanics of molecules, and the MES-PI formulation (established in 2018) is the most robust tool currently available to capture these effects naturally.