Breakthrough in Quantum Materials Modeling Incorporates Relativistic Effects

A significant advancement in computational materials science has enabled more accurate modeling of materials containing heavier elements. Researchers have successfully extended the phaseless auxiliary-field quantum Monte Carlo (pw-AFQMC) method to accurately incorporate spin-orbit coupling (SOC), a relativistic effect crucial for understanding the behavior of these materials.

The Challenge of Relativistic Effects in Materials Science

Relativistic effects, particularly spin-orbit coupling, become increasingly important in materials with heavier elements, influencing their magnetic, electronic, and topological characteristics. However, incorporating these effects into advanced simulations has historically been computationally challenging. This novel method bridges a longstanding gap in materials science, allowing for the concurrent treatment of both electronic correlation and SOC effects.

How the New Method Works

The research team, including Zheng Liu and Shiwei Zhang from The Center for Advanced Quantum Studies and School of Physics and Astronomy at Beijing Normal University, and the Center for Computational Quantum Physics at the Flatiron Institute, along with Fengjie Ma, integrated SOC using fully-relativistic pseudopotentials derived from Dirac-like equations. These pseudopotentials, optimized for multiple-projector norm-conservation, allow for accurate calculations even with the increased computational demands of including SOC.

The modified pw-AFQMC method utilizes a two-component Hamiltonian in the spinor basis, effectively doubling the size of the computational space to account for electron spin. The phaseless approximation, a key methodological innovation, avoids the ‘sign problem’ inherent in many quantum Monte Carlo simulations by restricting the trial wave function to be non-negative, enhancing the efficiency and stability of the calculations.

Validation and Applications

To demonstrate the accuracy of this approach, researchers computed the dissociation energy of iodine molecules (I₂) and the cohesive energy of lead (Pb). The calculations revealed a substantial influence of SOC on both properties. Specifically, the calculated dissociation energy of I₂ was 7.637 eV, an improvement over scalar-relativistic pseudopotential calculations which yielded 7.513 eV. The cohesive energy of lead was calculated as 2.341 eV, closely aligning with both density functional theory (DFT) calculations (2.378 eV) and experimental data (2.365 eV).

The team also applied their method to predict the transition pressure of indium phosphide (InP) as it changes from the zinc-blende to rock-salt crystal structure, accurately determining the pressure at which this phase transition occurs (14.2 GPa).

Implications for Future Technologies

This advancement promises to accelerate the discovery and design of novel materials for spintronics, topological quantum computing, and advanced magnetic storage, where precise understanding and modeling of SOC are paramount. The ability to accurately model materials with relativistic effects opens new avenues for exploring and engineering materials with tailored properties.

Research Team and Affiliations

The research was conducted by scientists affiliated with:

• The Center for Advanced Quantum Studies and School of Physics and Astronomy, Beijing Normal University
• The Center for Computational Quantum Physics, Flatiron Institute
• Department of Physics, Harvard University
• Department of Physics and Astronomy, Johns Hopkins University
• Institute for Theoretical Physics, University of Innsbruck

The work was contributed by Maine Christos, Zhu-Xi Luo, Henry Shackleton, Ya-Hui Zhang, Mathias S Scheurer, and Subir Sachdev.