Predictive Simulations of Ionization Energies of Solvated Halide Ions with Relativistic Embedded Equation of Motion Coupled Cluster Theory

Completely Multipolar Model for Many-Body Water-Ion and Ion-Ion Interactions

This work constructs an advanced force field, the Completely Multipolar Model (CMM), to quantitatively reproduce each term of an energy decomposition analysis (EDA) for aqueous solvated alkali metal cations and halide anions and their ion pairings. We find that all individual EDA terms remain well-approximated in the CMM for ion-water and ion-ion interactions, except for polarization, which shows errors due to the partial covalency of ion interactions near their equilibrium. We quantify the onset of the dative bonding regime by examining the change in molecular polarizability and Mayer bond indices as a function of distance, showing that partial covalency manifests by breaking the symmetry of atomic polarizabilities while strongly damping them at short-range. This motivates an environment-dependent atomic polarizability parameter that depends on the strength of the local electric field experienced by the ions to account for strong damping, with anisotropy introduced by atomic multipoles. The resulting CMM model for ions provides accurate dimer surfaces and three-body polarization and charge transfer compared to EDA, and shows excellent performance on various ion benchmarks including vibrational frequencies and cluster geometries.

Probing the Local Environment in Potassium Salts and Potassium-Promoted Catalysts by Potassium Valence-to-Core X-ray Emission Spectroscopy

Potassium plays an important role in biology as well as a promoter in heterogeneous catalysis. There are, however, limited characterization techniques for potassium available in the literature. This study elucidates the potential of element-selective X-ray emission spectroscopy (XES) for characterizing the coordination environment and the electronic properties of potassium. A series of XES measurements were conducted, primarily focusing on the VtC transition (Kβ2,5) of potassium halides (KCl, KBr, and KI) and oxide-bound potassium salts, including potassium nitrate (KNO3) and potassium carbonate (K2CO3). Across the series of potassium halides, the VtC transition energy is observed to increase, as accurately reproduced by TDDFT calculations. Molecular orbital analysis suggests that the Kβ2,5 transition is primarily derived from halide np contributions, with the primary factor influencing the energy shift being the metal-ligand distances. For oxide ligands, an additional Kβ″ transition appears alongside the Kβ2,5, which is attributed to a low-energy ligand ns, as elucidated by theoretical calculations. Finally, the XES spectra of two potassium-promoted catalysts for ammonia decomposition/synthesis were measured. These spectra show that potassium within the catalyst is distinct from other K salts in the VtC region, which could be promising for understanding the role of potassium as an electronic promoter.

Coupled-cluster treatment of complex open-shell systems: the case of single-molecule magnets

We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

Interoperable workflows by exchanging grid-based data between quantum-chemical program packages

Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability.

Static embedding with pair coupled cluster doubles based methods

Quantum embedding methods have recently been significantly developed to model large molecular structures. This work proposes a novel wave function theory in a density functional theory (WTF-in-DFT) embedding scheme based on pair-coupled cluster doubles (pCCD)-type methods. While pCCD can reliably describe strongly-correlated systems with mean-field-like computational cost, the large extent of the dynamic correlation can be accounted for by (linearized) coupled-cluster corrections on top of the pCCD wave function. Here we focus on the linearized coupled-cluster singles and doubles (LCCSD) ansatz for electronic ground states and its extension to excited states within the equation of motion (EOM) formalism. We test our EOM-pCCD-LCCSD-in-DFT approach for the vertical excitation energies of the hydrogen-bonded water-ammonia complex, micro-solvated thymine, and uranyl tetrahalides (UO2X42-, X = F, Cl, Br). Furthermore, we assess the quality of the embedding potential using an orbital entanglement and correlation analysis. The approximate embedding models successfully capture changes in the excitation energies going from bare fragments to supramolecular structures and represent a promising computational method for excited states in large molecular systems.

Core Excitations of Uranyl in Cs2UO2Cl4 from Relativistic Embedded Damped Response Time-Dependent Density Functional Theory Calculations

X-ray spectroscopies, by their high selectivity and sensitivity to the chemical environment around the atoms probed, provide significant insights into the electronic structures of molecules and materials. Interpreting experimental results requires reliable theoretical models, accounting for environmental, relativistic, electron correlation, and orbital relaxation effects in a balanced manner. In this work, we present a protocol for the simulation of core excited spectra with damped response time-dependent density functional theory based on the Dirac-Coulomb Hamiltonian (4c-DR-TD-DFT), in which environmental effects are accounted for through the frozen density embedding (FDE) method. We showcase this approach for the uranium M4- and L3-edges and oxygen K-edge of the uranyl tetrachloride (UO2Cl42-) unit as found in a host Cs2UO2Cl4 crystal. We have found that the 4c-DR-TD-DFT simulations yield excitation spectra that very closely match the experiment for the uranium M4-edge and the oxygen K-edge, with good agreement for the broad experimental spectra for the L3-edge. By decomposing the complex polarizability in terms of its components, we have been able to correlate our results with angle-resolved spectra. We have observed that for all edges, but in particular the uranium M4-edge, an embedded model in which the chloride ligands are replaced by an embedding potential reproduces rather well the spectral profile obtained for UO2Cl42-. Our results underscore the importance of the equatorial ligands to simulating core spectra at both uranium and oxygen edges.

Solvation, Surface Propensity, and Chemical Reactions of Solutes at Atmospheric Liquid-Vapor Interfaces

surface is covered by oceans, a large number of liquid aerosol particles fill the air, and clouds hold a tiny but critical fraction of Earth's water in the air to influence our climate and hydrology, enabling the lives of humans and ecosystems. The surfaces of these liquids provide the interface for the transfer of gases, for nucleation processes, and for catalyzing important chemical reactions. Coupling a range of spectroscopic tools to liquid microjets has become an important approach to better understanding dynamics, structure, and chemistry at liquid interfaces. Liquid microjets offer stability in vacuum and ambient pressure environments, thus also allowing X-ray photoelectron spectroscopy (XPS) with manageable efforts in terms of differential pumping. Liquid microjets are operated at speeds sufficient to allow for a locally equilibrated surface in terms of water dynamics and solute surface partitioning. XPS is based on the emission of core-level electrons, the binding energy of which is selective for the element and its chemical environment. Inelastic scattering of electrons establishes the probing depth of XPS in the nanometer range and thus its surface sensitivity.In this Account, we focus on aqueous solutions relevant to the surface of oceans, aqueous aerosols, or cloudwater. We are interested in understanding solvation and acid dissociation at the interface, interfacial aspects of reactions with gas-phase reactants, and the interplay of ions with organic molecules at the interface. The strategy is to obtain a link between the molecular-level picture and macroscopic properties and reactivity in the atmospheric context.We show consistency between surface tension and XPS for a range of surface-active organic species as an important proof for interrogating an equilibrated liquid surface. Measurements with organic acids and amines offer important insight into the question of apparent acidity or basicity at the interface. Liquid microjet XPS has settled the debate of the surface enhancement of halide ions, shown using the example of bromide and its oxidation products. Despite the absence of a strong enhancement for the bromide ion, its rate of oxidation by ozone is surface catalyzed through the stabilization of the bromide ozonide intermediate at the interface. In another reaction system, the one between Fe2+ and H2O2, a similar intermediate in the form of highly valent iron species could not be detected by XPS under the experimental conditions employed, shedding light on the abundance of this intermediate in the environment but also on the constraints within which surface species can be detected. Emphasizing the importance of electrostatic effects, we show how a cationic surfactant attracts charged bromide anions to the interface, accompanied by enhanced oxidation rates by ozone, overriding the role of surfactants as a barrier for the access of gas-phase reactants. The reactivity and structure at interfaces thus result from a subtle balance between hygroscopic and hydrophobic interactions, electrostatic effects, and the structural properties of both liquids and solutes.

Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques

Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac-Kohn-Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA's python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3-H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.

Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution via relativistic quantum embedding calculations

In this work, we investigate the effects of the environment on the X-ray photoelectron spectra of hydrogen chloride and chloride ions adsorbed on ice surfaces, as well as of chloride ions in water droplets. In our approach, we combine a density functional theory (DFT) description of the ice surface with that of halogen species using the recently developed relativistic core-valence separation equation of motion coupled cluster (CVS-EOM-IP-CCSD) via the frozen density embedding formalism (FDE), to determine the K and L_1,2,3 edges of chlorine. Our calculations, which incorporate temperature effects through snapshots from classical molecular dynamics simulations, are shown to reproduce experimental trends in the change of the core binding energies of Cl^- upon moving from a liquid (water droplets) to an interfacial (ice quasi-liquid layer) environment. Our simulations yield water valence band binding energies in good agreement with experiment, which vary little between the droplets and the ice surface. For halide core binding energies there is an overall trend for overestimating experimental values, though good agreement between theory and experiment is found for Cl^- in water droplets and on ice. For HCl on the other hand there are significant discrepancies between experimental and calculated core binding energies when we consider structural models that maintain the H-Cl bond more or less intact. An analysis of models that allow for pre-dissociated and dissociated structures suggests that experimentally observed chemical shifts in binding energies between Cl^- and HCl would require that H⁺ (in the form of H₃O⁺) and Cl^- are separated by roughly 4-6 Å.

Environment Effects on X-Ray Absorption Spectra With Quantum Embedded Real-Time Time-Dependent Density Functional Theory Approaches

In this work we implement the real-time time-dependent block-orthogonalized Manby-Miller embedding (rt-BOMME) approach alongside our previously developed real-time frozen density embedding time-dependent density functional theory (rt-TDDFT-in-DFT FDE) code, and investigate these methods' performance in reproducing X-ray absorption spectra (XAS) obtained with standard rt-TDDFT simulations, for model systems comprised of solvated fluoride and chloride ions ([X@( H 2 O ) 8 - , X = F, Cl). We observe that for ground-state quantities such as core orbital energies, the BOMME approach shows significantly better agreement with supermolecular results than FDE for the strongly interacting fluoride system, while for chloride the two embedding approaches show more similar results. For the excited states, we see that while FDE (constrained not to have the environment densities relaxed in the ground state) is in good agreement with the reference calculations for the region around the K and L1 edges, and is capable of reproducing the splitting of the 1s1 (n + 1)p1 final states (n + 1 being the lowest virtual p orbital of the halides), it by and large fails to properly reproduce the 1s1 (n + 2)p1 states and misses the electronic states arising from excitation to orbitals with important contributions from the solvent. The BOMME results, on the other hand, provide a faithful qualitative representation of the spectra in all energy regions considered, though its intrinsic approximation of employing a lower-accuracy exchange-correlation functional for the environment induces non-negligible shifts in peak positions for the excitations from the halide to the environment. Our results thus confirm that QM/QM embedding approaches are viable alternatives to standard real-time simulations of X-ray absorption spectra of species in complex or confined environments.

Core and Valence Level Photoelectron Spectroscopy of Nanosolvated KCl

The solvation of alkali and halide ions in the aqueous environment has been a subject of intense experimental and theoretical research with multidisciplinary interests; yet, a comprehensive molecular-level understanding has still not been obtained. In recent years, electron spectroscopy has been increasingly applied to study the electronic and structural properties of aqueous ions with implications, especially in atmospheric chemistry. In this work, we report core and valence level (Cl 2p, Cl 3p, and K 3p) photoelectron spectra of the common alkali halide, KCl, doped in gas-phase water clusters in the size range of a few hundred water molecules. The results indicate that the electronic structure of these nanosolutions shows a distinct character from that observed at the liquid-vapor interface in liquid microjets and ambient pressure setups. Insights are provided into the unique solvation properties of ions in a nanoaqueous environment, emerging properties of bulk electrolyte solutions with growing cluster size, and sensitivity of the electronic structure to varying solvation configurations.

Magnetic Anisotropy in a Cubane-like Ni4 Complex: An Ab Initio Perspective

Magnetic anisotropy, in the absence of an external magnetic field, relates to the degeneracy lift of energy levels. In the standard case of transition metal complexes, this property is usually modeled by an anisotropic spin Hamiltonian and one speaks of "zero-field splitting" (ZFS) of spin states. While the case of mononuclear complexes has been extensively described by means of ab initio quantum mechanical calculations, the literature on polynuclear complexes studied with these methodologies is rather scarce. In this work, advanced multiconfigurational wave function theory methods are applied to compute the ZFS of the ground S = 4 state of an actual tetranickel(II) complex, displaying a magnet behavior below 0.5 K. First, the isotropic couplings are computed in the absence of the spin-orbit coupling operator, in the full complex and also in clusters with only two active nickel(II) centers, confirming the occurrence of weak ferromagnetic couplings in this system. Second, the single-site magnetic anisotropies are computed on a cluster bearing only one active nickel(II) site, showing that the single-site anisotropy axes are not oriented in an optimal fashion for generating a large uniaxial molecular anisotropy. Furthermore, the possibility for involving only a few local orbital excited states in the calculation is assessed, actually opening the way for a consistent and manageable treatment of the ZFS of the ground S = 4 state. Third, multiconfigurational calculations are performed on the full complex, confirming the weak uniaxial anisotropy occurring for this state and also, interestingly, revealing a significant contribution of the lowest-lying orbitally excited S = 3 states. Overall, by comparison with the experiment, the reported results question the common habit of using only one structure, in particular derived from a crystallography experiment, to compute magnetic anisotropy parameters.

Generalization of Intrinsic Orbitals to Kramers-Paired Quaternion Spinors, Molecular Fragments, and Valence Virtual Spinors

Localization of molecular orbitals finds its importance in the representation of chemical bonding (and antibonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed on this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic acid, and ferrocene molecules, and then demonstrate the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium, and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.

Environmental Effects with Frozen-Density Embedding in Real-Time Time-Dependent Density Functional Theory Using Localized Basis Functions

Frozen-density embedding (FDE) represents a versatile embedding scheme to describe the environmental effect on electron dynamics in molecular systems. The extension of the general theory of FDE to the real-time time-dependent Kohn-Sham method has previously been presented and implemented in plane waves and periodic boundary conditions [Pavanello, M.; J. Chem. Phys. 2015, 142, 154116]. In the current paper, we extend our recent formulation of the real-time time-dependent Kohn-Sham method based on localized basis set functions and developed within the Psi4NumPy framework to the FDE scheme. The latter has been implemented in its "uncoupled" flavor (in which the time evolution is only carried out for the active subsystem, while the environment subsystems remain at their ground state), using and adapting the FDE implementation already available in the PyEmbed module of the scripting framework PyADF. The implementation was facilitated by the fact that both Psi4NumPy and PyADF, being native Python API, provided an ideal framework of development using the Python advantages in terms of code readability and reusability. We employed this new implementation to investigate the stability of the time-propagation procedure, which is based on an efficient predictor/corrector second-order midpoint Magnus propagator employing an exact diagonalization, in combination with the FDE scheme. We demonstrate that the inclusion of the FDE potential does not introduce any numerical instability in time propagation of the density matrix of the active subsystem, and in the limit of the weak external field, the numerical results for low-lying transition energies are consistent with those obtained using the reference FDE calculations based on the linear-response TDDFT. The method is found to give stable numerical results also in the presence of a strong external field inducing nonlinear effects. Preliminary results are reported for high harmonic generation (HHG) of a water molecule embedded in a small water cluster. The effect of the embedding potential is evident in the HHG spectrum reducing the number of the well-resolved high harmonics at high energy with respect to the free water. This is consistent with a shift toward lower ionization energy passing from an isolated water molecule to a small water cluster. The computational burden for the propagation step increases approximately linearly with the size of the surrounding frozen environment. Furthermore, we have also shown that the updating frequency of the embedding potential may be significantly reduced, much less than one per time step, without jeopardizing the accuracy of the transition energies.

Hydration of NH_{4}^{+} in Water: Bifurcated Hydrogen Bonding Structures and Fast Rotational Dynamics

Understanding the hydration and diffusion of ions in water at the molecular level is a topic of widespread importance. The ammonium ion (NH_{4}^{+}) is an exemplar system that has received attention for decades because of its complex hydration structure and relevance in industry. Here we report a study of the hydration and the rotational diffusion of NH_{4}^{+} in water using ab initio molecular dynamics simulations and quantum Monte Carlo calculations. We find that the hydration structure of NH_{4}^{+} features bifurcated hydrogen bonds, which leads to a rotational mechanism involving the simultaneous switching of a pair of bifurcated hydrogen bonds. The proposed hydration structure and rotational mechanism are supported by existing experimental measurements, and they also help to rationalize the measured fast rotation of NH_{4}^{+} in water. This study highlights how subtle changes in the electronic structure of hydrogen bonds impacts the hydration structure, which consequently affects the dynamics of ions and molecules in hydrogen bonded systems.

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.