Charge Transfer with a Spin. II: A Framework for Diabatization which Localizes Charge and Spin

This paper presents a rapidly converging, two-step diabatization framework using complex-unitary rotations to localize both charge and spin in open-shell systems with spin-orbit coupling, thereby generating smooth diabatic potential energy surfaces that preserve dipole and spin properties while maintaining time-reversal symmetry.

The Gist

The Big Picture: Untangling a Messy Dance

Imagine you are trying to watch a complex dance performance. The dancers (electrons) are moving incredibly fast, and they are holding hands with invisible partners (their "spin" and their "charge").

In the world of chemistry, scientists usually try to describe these dancers by watching their exact, instantaneous positions. This is called the Adiabatic view. It's accurate, but it's chaotic. As the dancers move, they suddenly swap partners, spin in different directions, and their energy levels jump around wildly. It's like trying to film a dance where the camera shakes so much you can't tell who is leading or who is following.

This paper introduces a new way to film the dance. The authors want to create a Diabatic view. Think of this as a "steady cam" shot. They want to smooth out the camera work so that:

1. Charge Localization: You can clearly see which dancer is on the "Left" side of the stage and which is on the "Right" side, even as they move.
2. Spin Localization: You can clearly see which dancer is spinning "Up" and which is spinning "Down," without the camera spinning wildly around them.

The Problem: The "Magic Spin" (Spin-Orbit Coupling)

In many molecules, especially those with an odd number of electrons (like a single unpaired dancer), there is a tricky force called Spin-Orbit Coupling (SOC).

Imagine the dancers are wearing magnetic gloves. As they move across the stage (nuclear motion), the magnetic field changes. This causes their "spin" (the direction they are pointing) to twist and turn in a way that doesn't make intuitive sense.

• The Old Way: Previous methods tried to track these dancers, but because the "spin" was twisting so fast, the data looked like static on a TV screen. You couldn't tell if the spin changed because of the physics or just because the camera angle was weird.
• The Challenge: The authors needed a way to untangle the "charge" (where the electron is) from the "spin" (which way it's pointing) without breaking the laws of physics (specifically, Time-Reversal Symmetry, which is like a rule that says the dance looks the same whether played forward or backward).

The Solution: The "Two-Step" Smoothie Maker

The authors developed a mathematical recipe (a framework) to fix this. They call it a two-step optimization. Imagine you are making a smoothie, and you want to separate the fruit chunks from the juice perfectly.

Step 1: The "Charge" Filter (Dipole Localization)
First, they apply a filter that forces the electrons to stay on their respective sides.

• Analogy: Imagine two buckets, one on the left and one on the right. The math forces the electrons to stay in the left bucket or the right bucket, maximizing the difference between them. This creates a smooth path where the electron clearly moves from Left to Right without getting confused.

Step 2: The "Spin" Filter (Spin Localization)
Once the electrons are in the right buckets, the authors apply a second filter to the spin.

• Analogy: Imagine the dancers are wearing hats. Sometimes the hats spin wildly as the dancers move. The authors' math gently "tunes" the hats so they point in a smooth, consistent direction as the dancer moves across the stage. They don't stop the spin from changing (because physics demands it), but they stop the random, chaotic spinning caused by the camera angle. They make the spin axis "drift" slowly and smoothly, like a compass needle slowly turning as you walk north, rather than spinning like a top.

The "Kramers Pair" Mystery

The paper deals with a specific rule in quantum mechanics called Kramers Degeneracy.

• The Metaphor: Imagine every dancer has a "shadow twin." If you have one dancer, you must have a twin. They are identical in energy but opposite in spin.
• The Problem: In the old methods, you could swap the dancer and their twin at any time, and the math would still work, but the "spin direction" would jump around randomly.
• The Fix: The authors' method treats the dancer and the twin as a single unit. They rotate them together in a way that keeps the "shadow twin" relationship intact while still smoothing out the spin direction. This ensures the physics remains honest (Time-Reversal Symmetry is preserved).

Why Does This Matter? (The "Chiral" Connection)

Why do we care about smoothing out these spins?
The authors hint at a phenomenon called Chiral Induced Spin Selectivity (CISS).

• The Analogy: Imagine a spiral staircase (a chiral molecule). If you walk up the stairs, you naturally twist your body to the right. In the molecular world, electrons moving through a "spiral" molecule might naturally prefer to spin one way (Up) over the other (Down).
• The Impact: By using this new "steady cam" method, scientists can finally see exactly how the electron's spin twists as it moves through these spiral molecules. This could help us build better solar cells, faster computer chips, and new medical treatments that rely on controlling electron spin.

Summary

This paper is about building a better camera for the quantum world.

1. The Problem: Electrons in certain molecules are messy; their location and spin get tangled up, making it hard to predict how they move.
2. The Tool: A new mathematical algorithm that acts like a "steady cam."
3. The Result: It produces a smooth, easy-to-understand map of how electrons move and spin. It separates the "where" (charge) from the "how it spins" (spin) while respecting the fundamental rules of the universe.
4. The Future: This tool will help scientists understand how nature uses spin to control chemical reactions, potentially leading to revolutionary new technologies.

Technical Summary

1. Problem Statement

The paper addresses the challenge of constructing diabatic states for open-shell molecular systems undergoing charge transfer (CT) in the presence of Spin-Orbit Coupling (SOC).

• The Adiabatic Limitation: Standard adiabatic electronic states (eigenstates of the Hamiltonian) are convenient for calculation but ill-suited for modeling CT and nonadiabatic dynamics. Near avoided crossings, their charge distributions and spin expectation values change abruptly, making it difficult to track electron transfer or interpret spin dynamics.
• The SOC Challenge: Traditional diabatization methods (e.g., Boys localization, Edmiston–Ruedenberg, fragment charge diabatization) are designed for spin-free systems. They fail in the presence of strong SOC because:
  1. SOC breaks spin conservation and entangles spatial and spin degrees of freedom.
  2. Open-shell systems with an odd number of electrons exhibit Kramers degeneracy (time-reversal symmetry), resulting in complex-valued spinors rather than real orbitals.
  3. Existing methods often cannot handle the degenerate, complex-valued nature of these states or fail to distinguish between different spin states within a manifold.
• The Specific Goal: The authors aim to develop a unified framework that simultaneously localizes charge (in real space) and localizes spin (in spin space) while preserving Time-Reversal Symmetry (TRS). This ensures that the spin quantization axis varies smoothly along the reaction coordinate rather than twisting arbitrarily due to gauge freedom.

2. Methodology

The authors propose a two-step diabatization procedure based on complex-unitary Jacobi sweeps. The method transforms a set of adiabatic states (specifically two Kramers doublets, forming a 4-state crossing) into a diabatic basis.

A. Theoretical Framework

• System: The method targets systems with an odd number of electrons (Kramers doublets). The transformation matrix $U$ must be a complex-valued unitary matrix (specifically within the $SU(2)$ subgroup for Kramers pairs) to preserve the Kramers structure and TRS.
• Two-Step Optimization: The procedure sequentially maximizes two objective functions:
  1. Charge Localization (Dipole Maximization): Maximizes the magnitude of the dipole moment difference between states. This is achieved by approximately joint diagonalizing the dipole moment matrices ($\mu_x, \mu_y, \mu_z$).
  2. Spin Localization: Maximizes the magnitude of the spin expectation values ($\langle S_x \rangle, \langle S_y \rangle, \langle S_z \rangle$) within each Kramers pair. This aligns the local spin quantization axis smoothly.

B. Algorithm: Jacobi Sweeps with TRS Constraints

The optimization is performed via iterative Jacobi sweeps (successive $2\times2$ or $4\times4$ rotations) rather than optimizing the full matrix at once.

• Type I Rotations (Spin Diabatization): Rotations occur within a single Kramers pair $(i, \bar{i})$. These are $2\times2$ $SU(2)$ rotations that mix the time-reversed partners to maximize spin localization without breaking the Kramers relation ($\hat{T}|\xi_i\rangle = |\bar{\xi}_i\rangle$).
• Type II Rotations (Charge Diabatization): Rotations occur between two different Kramers pairs $(i, \bar{i})$ and $(j, \bar{j})$. To preserve TRS, a rotation mixing $|i\rangle$ and $|j\rangle$ must be accompanied by a conjugate rotation mixing their time-reversed partners $|\bar{i}\rangle$ and $|\bar{j}\rangle$. This requires a $4\times4$ block-diagonal complex unitary matrix.
• Solution: The optimal rotation angles ($\gamma, \phi$) are found by solving a cubic equation derived from maximizing the localization functional, extending previous real-valued solutions to the complex domain.

C. Implementation

• Electronic Structure: The adiabatic states are generated using the electron/hole-transfer Dynamically-weighted State-Averaged Constrained CASSCF (eDSC/hDSC) method, which handles odd-electron systems and complex orbitals.
• Software: Implemented in a development version of Q-Chem.

3. Key Contributions

1. Unified Charge-Spin Diabatization: The first framework to simultaneously localize charge and spin for open-shell systems with SOC, generalizing the Boys localization concept to include spin operators.
2. Complex-Unitary Formalism: Development of specific update rules for complex-unitary rotations that strictly preserve Time-Reversal Symmetry and Kramers degeneracy, addressing a gap in existing diabatization literature.
3. Smooth Pseudospin Texture: The method yields diabatic states where the spin expectation vector (pseudospin) varies smoothly along the reaction coordinate, removing arbitrary gauge fluctuations inherent in degenerate manifolds.
4. Analytical Solution for Complex Rotations: Derivation of a cubic equation to solve for optimal rotation parameters in the complex domain, enabling rapid convergence.

4. Results

The method was applied to a phenoxy-phenol (ph-ph) system, a proton-coupled electron transfer (PCET) model where a hydrogen atom bridges two phenyl fragments.

• Potential Energy Surfaces (PES): The resulting diabatic PESs are smooth and exhibit well-defined crossings, even across a wide range of SOC strengths (scaling factors $\eta = 1$ to $200$). The charge-localized states (Left vs. Right) maintain distinct electronic characters.
• Spin Expectation Values:
  • In the diabatic basis, the spin expectation vectors $\langle \mathbf{S} \rangle$ vary smoothly with nuclear geometry.
  • The method successfully fixes the "gauge" of the spin quantization axis, preventing the spurious, rapid fluctuations seen in adiabatic representations.
  • Time-reversal symmetry is preserved: time-reversed partners exhibit equal magnitude but opposite sign spin expectation values.
• Spin Contamination ($\langle S^2 \rangle$):
  • The study analyzes the deviation of $\langle S^2 \rangle$ from the ideal doublet value ($0.75$).
  • Results show that the diabatization procedure introduces negligible spin contamination (deviation $< 0.05\%$ for weak SOC).
  • Crucially, the diabatic and adiabatic states have nearly identical $\langle S^2 \rangle$ values, confirming that the transformation does not artificially break spin symmetry.
• Pseudospin Texture: The authors quantify the rotation of the pseudospin axis along the reaction coordinate using a signed angle $\theta(R)$. This reveals a non-trivial "pseudospin texture" generated by the interplay of nuclear motion and SOC.

5. Significance and Future Implications

• Modeling Nonadiabatic Dynamics: By providing smooth diabatic states with well-defined charge and spin characters, this framework enables accurate modeling of nonadiabatic population transfer in radical systems, which is essential for understanding photochemistry and electron transfer.
• Chiral Induced Spin Selectivity (CISS): The ability to track the smooth rotation of the spin axis is critical for investigating CISS, where electron transfer in chiral environments exhibits spin selectivity. The authors suggest this method can correlate orbital helicity with spin-axis rotation.
• Phase Space Electronic Structure: The smooth diabatic basis serves as an ideal starting point for advanced theories involving nuclear-electronic correlation and internal magnetic fields (phase space electronic structure theory).
• Generalizability: While tested on a 4-state crossing, the framework is designed to be extensible to larger systems with multiple charge centers and complex spin manifolds.

In summary, this paper provides a rigorous mathematical and computational solution to the "spin localization" problem in open-shell charge transfer, bridging the gap between accurate electronic structure calculations and physically intuitive diabatic models for spin-active systems.