Efficient Four-Component Dirac-Coulomb-Gaunt Hartree-Fock in the Pauli Spinor Representation.

Four-component Dirac-Hartree-Fock is an accurate mean-field method for treating molecular systems where relativistic effects are important. However, the computational cost and complexity of the two-electron interaction make this method less common, even though we can consider the Dirac-Hartree-Fock Hamiltonian the "ground truth" of the electronic structure, barring explicit quantum electrodynamical effects. Being able to calculate these effects is then vital to the design of lower scaling methods for accurate predictions in computational spectroscopy and properties of heavy element complexes that must include relativistic effects for even qualitative accuracy. In this work, we present a Pauli quaternion formalism of maximal component and spin separation for computing the Dirac-Coulomb-Gaunt Hartree-Fock ground state, with a minimal floating point operation count algorithm. This approach also allows one to explicitly separate different spin physics from the two-body interactions, such as spin-free, spin-orbit, and spin-spin contributions. Additionally, we use this formalism to examine relativistic trends in the periodic table and analyze the basis set dependence of atomic gold and gold dimer systems.

Generalization of Block-Localized Wave Function for Constrained Optimization of Excited Determinants.

The block-localized wave function method is useful to provide insights on chemical bonding and intermolecular interactions through energy decomposition analysis. The method relies on block localization of molecular orbitals (MOs) by constraining the orbitals to basis functions within given blocks. Here, a generalized block-localized orbital (GBLO) method is described to allow both physically localized and delocalized MOs to be constrained in orbital-block definitions. Consequently, GBLO optimization can be conveniently tailored by imposing specific constraints. The GBLO method is illustrated by three examples: (1) constrained polarization response orbitals through dipole and quadrupole perturbation in a water dimer complex, (2) the ground and first excited-state potential energy curves of ethene about its C-C bond rotation, and (3) excitation energies of double electron excited states. Multistate density functional theory is used to determine the energies of the adiabatic ground and excited states using a minimal active space (MAS) comprising specifically charge-constrained and excited determinant configurations that are variationally optimized by the GBLO method. We find that the GBLO expansion that includes delocalized MOs in configurational blocks significantly reduces computational errors in comparison with physical block localization, and the computed ground- and excited-state energies are in good accordance with experiments and results obtained from multireference configuration interaction calculations.

Modeling L2,3-edge X-ray absorption spectroscopy with linear response exact two-component relativistic time-dependent density functional theory.

X-ray absorption spectroscopy (XAS) is a powerful tool that can provide physical insights into element-specific chemical processes and reactivities. Although relativistic time-dependent density functional theory (TDDFT) has been previously applied to model the L-edge region in XAS, there has not been a more comprehensive study of the choices of basis sets and density functional kernels available for variational relativistic excited state methods. In this work, we introduce the implementation of the generalized preconditioned locally harmonic residual algorithm to solve the complex-valued relativistic TDDFT for modeling the L-edge X-ray absorption spectra. We investigate the L2,3-edge spectra of a series of molecular complexes using relativistic linear response TDDFT with a hybrid iterative diagonalization algorithm. A systematic error analysis was carried out with a focus on the energetics, intensities, and magnitude of L2-L3 splitting compared to experiments. Additionally, the results from relativistic TDDFT calculations are compared to those computed using other theoretical methods, and the multideterminantal effects on the L-edge XAS were investigated.

The Role of Excited-State Proton Relays in the Photochemical Dynamics of Water Nanodroplets.

In this work, we applied nonadiabatic excited-state molecular dynamics in tandem with ab initio electronic structure theory to illustrate a complete mechanistic landscape underpinning the ultraviolet absorption-initiated photochemical dynamics in water nanodroplets. The goal is to understand the nonequilibrium excited-state molecular dynamics initiated by the relaxation of a solvated photoelectron and consequential photochemical processes. The lowest-lying excited state shows the proton dissociation for a single water molecule forming intermediate hydronium complexes through a proton relay. At approximately 100 fs, the proton relay process gives rise to the relaxation of the excited state accompanied by a rapid increase in the nonadiabatic coupling strength with the ground state, and the nanodroplet nonradiatively decays. The nonadiabatic transition to the ground state produces excited vibrational states that facilitate the recombination of the dissociated proton and hydroxyl group, eventually leading to the desorption of water molecules from the nanodroplet. Additionally, lifetimes of transient photochemical events are also resolved for the relaxation of a solvated electron, excited-state proton relay, and nonradiative transition.

Generalized Hartree-Fock with Nonperturbative Treatment of Strong Magnetic Fields: Application to Molecular Spin Phase Transitions.

In this work, we present a framework of an ab initio variational approach to effectively explore electronic spin phase transitions in molecular systems inside of a homogeneous magnetic field. In order to capture this phenomenon, the complex generalized Hartree-Fock ([Formula: see text]) method is used in the spinor formalism with London orbitals. Recursive algorithms for computing the one- and two-electron integrals of London orbitals are also provided. A Pauli matrix representation of the [Formula: see text] method is introduced to separate spin contributions from the scalar part of the Fock matrix. Next, spin phase transitions in two different molecular systems are investigated in the presence of a strong magnetic field. Noncollinear spin configurations are observed during the spin phase transitions in H2 and a dichromium complex, [(H3N)4Cr(OH)2Cr(NH3)4]4+, with an increase in magnetic field strength. The competing driving forces of exchange coupling and the spin Zeeman effect have been shown to govern the spin phase transition and its transition rate. Additionally, the energetic contributions of the spin Zeeman, orbital Zeeman, and diamagnetic terms to the potential energy surface are also analyzed.

Modeling L2,3-Edge X-ray Absorption Spectroscopy with Real-Time Exact Two-Component Relativistic Time-Dependent Density Functional Theory.

X-ray absorption spectroscopy is a powerful technique to probe local electronic and nuclear structure. There has been extensive theoretical work modeling K-edge spectra from first principles. However, modeling L-edge spectra directly with density functional theory poses a unique challenge requiring further study. Spin-orbit coupling must be included in the model, and a noncollinear density functional theory is required. Using the real-time exact two-component method, we are able to variationally include one-electron spin-orbit coupling terms when calculating the absorption spectrum. The abilities of different basis sets and density functionals to model spectra for both closed- and open-shell systems are investigated using SiCl4 and three transition metal complexes, TiCl4, CrO2Cl2, and [FeCl6]3-. Although we are working in the real-time framework, individual molecular orbital transitions can still be recovered by projecting the density onto the ground state molecular orbital space and separating contributions to the time evolving dipole moment.