Modeling interactions between rubidium atom and magnetometer cell wall molecules.

Magnetometer cell wall coat molecules play an important role in preserving the lifetime of pumped alkali metal atoms for use in magnetometers that are capable of measuring very small magnetic fields. The goal of this study is to help rationalize the design of the cell coat molecules. Rubidium-87 is studied in terms of its interaction with three template cell coat molecules: ethane, ethene, and methyltrichlorosilane (MeTS). Ab initio electronic structure methods are applied to investigate the effect that the coat molecules have on the 2S ground state and 2P excited state of 87Rb. We find that, from the ab initio results, the three template molecules have differing effects, with MeTS having the largest effect on the ground state and ethane or ethene having the largest effect on the non-degenerate excited states.

Symmetry and reactivity of π-systems in electric and magnetic fields: a perspective from conceptual DFT.

The extension of conceptual density-functional theory (conceptual DFT) to include external electromagnetic fields in chemical systems is utilised to investigate the effects of strong magnetic fields on the electronic charge distribution and its consequences on the reactivity of π-systems. Formaldehyde, H2CO, is considered as a prototypical example and current-density-functional theory (current-DFT) calculations are used to evaluate the electric dipole moment together with two principal local conceptual DFT descriptors, the electron density and the Fukui functions, which provide insight into how H2CO behaves chemically in a magnetic field. In particular, the symmetry properties of these quantities are analysed on the basis of group, representation, and corepresentation theories using a recently developed automatic program for symbolic symmetry analysis, QSYM2. This allows us to leverage the simple symmetry constraints on the macroscopic electric dipole moment components to make profound predictions on the more nuanced symmetry transformation properties of the microscopic frontier molecular orbitals (MOs), electron densities, and Fukui functions. This is especially useful for complex-valued MOs in magnetic fields whose detailed symmetry analyses lead us to define the new concepts of modular and phasal symmetry breaking. Through these concepts, the deep connection between the vanishing constraints on the electric dipole moment components and the symmetry of electron densities and Fukui functions can be formalised, and the inability of the magnetic field in all three principal orientations considered to induce asymmetry with respect to the molecular plane of H2CO can be understood from a molecular perspective. Furthermore, the detailed forms of the Fukui functions reveal a remarkable reversal in the direction of the dipole moment along the CO bond in the presence of a parallel or perpendicular magnetic field, the origin of which can be attributed to the mixing between the frontier MOs due to their subduced symmetries in magnetic fields. The findings in this work are also discussed in the wider context of a long-standing debate on the possibility to create enantioselectivity by external fields.

Capturing the electron-electron cusp with the coupling-constant averaged exchange-correlation hole: A case study for Hooke's atoms.

In density-functional theory, the exchange-correlation (XC) energy can be defined exactly through the coupling-constant (λ) averaged XC hole n̄xc(r,r'), representing the probability depletion of finding an electron at r' due to an electron at r. Accurate knowledge of n̄xc(r,r') has been crucial for developing XC energy density-functional approximations and understanding their performance for molecules and materials. However, there are very few systems for which accurate XC holes have been calculated since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of λ while the electron density remains fixed at the physical (λ = 1) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron-electron cusp using finite basis sets. Focusing on Hooke's atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the λ-integrated XC hole. Our results indicate that, for Hooke's atoms, the error introduced by the description of the electron-electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle-, and coupling-constant-averaged XC holes are also calculated and provide a benchmark against which the Perdew-Burke-Ernzerhof and local density approximation XC hole models are assessed.

QSym2: A Quantum Symbolic Symmetry Analysis Program for Electronic Structure.

Symmetry provides a powerful machinery to classify, interpret, and understand quantum-mechanical theories and results. However, most contemporary quantum chemistry packages lack the ability to handle degeneracy and symmetry breaking effects, especially in non-Abelian groups, and they are not able to characterize symmetry in the presence of external magnetic or electric fields. In this article, a program written in Rust entitled QSym2 that makes use of group and representation theories to provide symmetry analysis for a wide range of quantum-chemical calculations is introduced. With its ability to generate character tables symbolically on-the-fly and by making use of a generic symmetry-orbit-based representation analysis method formulated in this work, QSym2 is able to address all of these shortcomings. To illustrate these capabilities of QSym2, four sets of case studies are examined in detail in this article: (i) high-symmetry C84H64, C60, and B9- to demonstrate the analysis of degenerate molecular orbitals (MOs); (ii) octahedral Fe(CN)63- to demonstrate the analysis of symmetry-broken determinants and MOs; (iii) linear hydrogen fluoride in a magnetic field to demonstrate the analysis of magnetic symmetry; and (iv) equilateral H3+ to demonstrate the analysis of density symmetries.

Magnetic optical rotation from real-time simulations in finite magnetic fields.

We present a numerical approach to magnetic optical rotation based on real-time time-dependent electronic-structure theory. Not relying on perturbation expansions in the magnetic field strength, the formulation allows us to test the range of validity of the linear relation between the rotation angle per unit path length and the magnetic field strength that was established empirically by Verdet 160 years ago. Results obtained from time-dependent coupled-cluster and time-dependent current density-functional theory are presented for the closed-shell molecules H2, HF, and CO in magnetic fields up to 55 kT at standard temperature and pressure conditions. We find that Verdet's linearity remains valid up to roughly 10-20 kT, above which significant deviations from linearity are observed. Among the three current density-functional approximations tested in this work, the current-dependent Tao-Perdew-Staroverov-Scuseria hybrid functional performs the best in comparison with time-dependent coupled-cluster singles and doubles results for the magnetic optical rotation.

Semiempirical Methods for Molecular Systems in Strong Magnetic Fields.

A general scheme is presented to extend semiempirical methods to include the effects of arbitrary strength magnetic fields, while maintaining computational efficiency. The approach utilizes three main modifications; a London atomic orbital (LAO) basis set is introduced, field-dependent kinetic energy corrections are added to the model Hamiltonian, and spin-Zeeman interaction energy terms are included. The approach is applied to the widely available density-functional tight-binding method GFN1-xTB. Considering the basis set requirements for the kinetic energy corrections in a magnetic field leads to two variants: a single-basis approach GFN1-xTB-M0 and a dual-basis approach GFN1-xTB-M1. The LAO basis in the latter includes the appropriate nodal structure for an accurate representation of the kinetic energy corrections. The variants are assessed by benchmarking magnetizabilities and nuclear magnetic resonance shielding constants calculated using weak magnetic fields. Remarkably, the GFN1-xTB-M1 approach also exhibits excellent performance for strong fields, |B| ≤ 0.2B0 (B0 = 2.3505 × 105 T), recovering exotic features such as the para- to dia-magnetic transition in the BH molecule and the preferred electronic configuration, molecular conformation, and orientation of benzene. At stronger field strengths, |B| > 0.2B0, a degradation in the quality of the results is observed. The utility of GFN1-xTB-M1 is demonstrated by performing conformer searches in a range of field strengths for the cyclooctatetraene molecule, with GFN1-xTB-M1 capturing the transition from tub to planar conformations at high field, consistent with much more computationally demanding current-density functional theory calculations. Magnetically induced currents are also shown to be well described for the benzene and infinitene molecules, the latter demonstrating the flexibility and computational efficiency of the approach. The GFN1-xTB-M1 approach is a useful tool for the study of structure, conformation, and dynamics of large systems in magnetic fields at the semiempirical level as well as for preoptimization of molecular structure in ab initio calculations, enabling more efficient exploration of complex potential energy surfaces and reactivity in the presence of external fields.

Real-time time-dependent self-consistent field methods with dynamic magnetic fields.

The first finite basis set implementation of the real-time time-dependent self-consistent field method in a dynamic (time-dependent) magnetic field using London atomic orbitals (LAOs) is presented. The accuracy of the finite basis approach using LAOs is benchmarked against numerical results from the literature for the hydrogen atom and H2 in the presence of rapidly oscillating magnetic fields. This comparison is used to inform the choice of appropriate basis sets for studies under such conditions. Remarkably, relatively modest compact LAO basis sets are sufficient to obtain accurate results. Analysis of electron dynamics in the hydrogen atom shows that LAO calculations correctly capture the time evolution of orbital occupations. The Fourier transformation of the autocorrelation function yields a power spectrum exhibiting harmonics associated with coherent emission, which closely matches the literature and further confirms the accuracy of this approach. The dynamical response of the electron density in H2 for a magnetic field parallel to the internuclear axis shows similar behavior to benchmark studies. The flexibility of this implementation is then demonstrated by considering how the dynamical response changes as a function of the orientation of the molecule relative to the applied field. At non-parallel orientations, the symmetry of the system is lowered and numerical benchmark data, which exploit cylindrical symmetry, are no-longer readily available. The present study demonstrates the utility of LAO-based calculations for extreme dynamic magnetic fields, providing a stress-test on the choice of basis. Future applications of this approach for less extreme dynamic magnetic fields are briefly discussed.

Exact constraints and appropriate norms in machine-learned exchange-correlation functionals.

Machine learning techniques have received growing attention as an alternative strategy for developing general-purpose density functional approximations, augmenting the historically successful approach of human-designed functionals derived to obey mathematical constraints known for the exact exchange-correlation functional. More recently, efforts have been made to reconcile the two techniques, integrating machine learning and exact-constraint satisfaction. We continue this integrated approach, designing a deep neural network that exploits the exact constraint and appropriate norm philosophy to de-orbitalize the strongly constrained and appropriately normed (SCAN) functional. The deep neural network is trained to replicate the SCAN functional from only electron density and local derivative information, avoiding the use of the orbital-dependent kinetic energy density. The performance and transferability of the machine-learned functional are demonstrated for molecular and periodic systems.

An Embedded Fragment Method for Molecules in Strong Magnetic Fields.

An extension of the embedded fragment method for calculations on molecular clusters is presented, which includes strong external magnetic fields. The approach is flexible, allowing for calculations at the Hartree-Fock, current-density-functional theory, Møller-Plesset perturbation theory, and coupled-cluster levels using London atomic orbitals. For systems consisting of discrete molecular subunits, calculations using London atomic orbitals can be performed in a computationally tractable manner for systems beyond the reach of conventional calculations, even those accelerated by resolution-of-the-identity or Cholesky decomposition methods. To assess the applicability of the approach, applications to water clusters are presented, showing how strong magnetic fields enhance binding within the clusters. However, our calculations suggest that, contrary to previous suggestions in the literature, this enhanced binding may not be directly attributable to strengthening of hydrogen bonding. Instead, these results suggest that this arises for larger field strengths as a response of the system to the presence of the external field, which induces a charge density build up between the monomer units. The approach is embarrassingly parallel and its computational tractability is demonstrated for clusters of up to 103 water molecules in triple-ζ basis sets, which would correspond to conventional calculations with more than 12 000 basis functions.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Structural and Electronic Studies of Substituted m-Terphenyl Group 12 Complexes.

The effects of para-substitution on the structural and electronic properties of four series of two-coordinate m-terphenyl Group 12 complexes (R-Ar#)2M (M = Zn, Cd, Hg; R = t-Bu 1-3, SiMe3 4-6, Cl 7-9, CF3 10-12, where R-Ar# = 2,6-{2,6-Xyl}2-4-R-C6H2 and 2,6-Xyl = 2,6-Me2C6H3) have been investigated. X-ray crystallography shows little structural variation across the series, with no significant change in the C-M-C bond distances and angles. However, considerable electronic differences are revealed by heteronuclear nuclear magnetic resonance (NMR) spectroscopy; a linear correlation is observed between the 113Cd, 199Hg, and 1H (2,6-Xyl methyl protons) NMR chemical shifts of the para-substituted complexes and the Hammett constants for the R-substituents. Specifically, an upfield shift in the NMR signal is observed with increasingly electron-withdrawing R-substituents. Density functional theory (DFT) calculations are employed to attempt to rationalize these trends.

Topological Analysis of Functions on Arbitrary Grids: Applications to Quantum Chemistry.

Algorithms are presented for performing a topological analysis of an arbitrary function, evaluated on an arbitrary grid of points. These algorithms work strictly by post-processing the data and require no additional function evaluations. This is achieved by connecting the grid points with a neighborhood graph, allowing the topological analysis to be recast as a problem in the graph theory. The flexibility of the approach is demonstrated for various applications involving analysis of the charge and magnetically induced current densities in molecules, where features of the neighborhood graph are found to correspond to chemically relevant topographical properties, such as Bader charges. These properties converge using orders of magnitude fewer grid points than uniform-grid approaches while exhibiting an appealing O[N log(N)] scaling of the computational cost. The issue of grid bias is discussed in the context of graph-based algorithms and strategies for avoiding this bias are presented. Python implementations of the algorithms are provided.

Revealing the exotic structure of molecules in strong magnetic fields.

A novel implementation for the calculation of molecular gradients under strong magnetic fields is employed at the current-density functional theory level to optimize the geometries of molecular structures, which change significantly under these conditions. An analog of the ab initio random structure search is utilized to determine the ground-state equilibrium geometries for Hen and CHn systems at high magnetic field strengths, revealing the most stable structures to be those in high-spin states with a planar geometry aligned perpendicular to the field. The electron and current densities for these systems have also been investigated to develop an explanation of chemical bonding in the strong field regime, providing an insight into the exotic chemistry present in these extreme environments.

Extending conceptual DFT to include external variables: the influence of magnetic fields.

An extension of conceptual DFT to include the influence of an external magnetic field is proposed in the context of a program set up to cope with the ever increasing variability of reaction conditions and concomitant reactivity. The two simplest global reactivity descriptors, the electronic chemical potential (µ) and the hardness (η), are considered for the main group atoms H-Kr using current density-functional theory. The magnetic field strength, |B|, is varied between 0.0 and 1.0 B 0 = he -1 a 0 -2 ≈ 2.3505 × 105 T, encompassing the Coulomb and intermediate regimes. The carbon atom is studied as an exemplar system to gain insight into the behaviour of the neutral, cationic and anionic species under these conditions. Their electronic configurations change with increasing |B|, leading to a piecewise behaviour of the ionization energy (I) and electron affinity (A) values as a function of |B|. This results in complex behaviour of properties such as the electronegativity χ = -1/2(I + A) = -µ and hardness η = 1/2(I - A). This raises an interesting question: to what extent are atomic properties periodic in the presence of a magnetic field? In the Coulomb regime, close to |B| = 0, we find the familiar periodicity of the atomic properties, and make the connections to response functions central to conceptual DFT. However, as the field increases in the intermediate regime configurational changes of the atomic species lead to discontinuous changes in their properties; fundamentally changing their behaviour, which is illustrated by constructing a periodic table of χ and η values at |B| = 0.5 B 0. These values tend to increase for groups 1-2 and decrease for groups 16-18, leading to a narrower range overall and suggesting substantial changes in the chemistry of the main group elements. Changes within each group are also examined as a function of |B|. These are more complex to interpret due to the larger number of configurations accessible to heavier elements at high field. This is illustrated for group 17 where Cl and Br have qualitatively different configurations to their lighter cogener at |B| = 0.5 B 0. The insight into periodic trends in strong magnetic fields may provide a crucial starting point for predicting chemical reactivity under these exotic conditions.

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: a Comparison between Energy- and Variance-Based Approaches.

Self-consistent field methods for excited states offer an attractive low-cost route to study not only excitation energies but also properties of excited states. Here, we present the generalization of two self-consistent field methods, the maximum overlap method (MOM) and the σ-SCF method, to calculate excited states in strong magnetic fields and investigate their stability and accuracy in this context. These methods use different strategies to overcome the well-known variational collapse of energy-based optimizations to the lowest solution of a given symmetry. The MOM tackles this problem in the definition of the orbital occupations to constrain the self-consistent field procedure to converge on excited states, while the σ-SCF method is based on the minimization of the variance instead of the energy. To overcome the high computational cost of the variance minimization, we present a new implementation of the σ-SCF method with the resolution of identity approximation, allowing the use of large basis sets, which is an important requirement for calculations in strong magnetic fields. The accuracy of these methods is assessed by comparison with the benchmark literature data for He, H2, and CH+. The results reveal severe limitations of the variance-based scheme, which become more acute in large basis sets. In particular, many states are not accessible using variance optimization. Detailed analysis shows that this is a general feature of variance optimization approaches due to the masking of local minima in the optimization. In contrast, the MOM shows promising performance for computing excited states under these conditions, yielding results consistent with available benchmark data for a diverse range of electronic states.

Modeling Ultrafast Electron Dynamics in Strong Magnetic Fields Using Real-Time Time-Dependent Electronic Structure Methods.

An implementation of real-time time-dependent Hartree-Fock (RT-TDHF) and current density functional theory (RT-TDCDFT) for molecules in strong uniform magnetic fields is presented. In contrast to earlier implementations, the present work enables the use of the RT-TDCDFT formalism, which explicitly includes field-dependent terms in the exchange-correlation functional. A range of current-dependent exchange-correlation functionals based on the TPSS functional are considered, including a range-separated variant, which is particularly suitable for application to excited state calculations. The performance of a wide range of propagator algorithms for real-time methods is investigated in this context. A recently proposed molecular orbital pair decomposition analysis allows for assignment of electronic transitions, providing detailed information about which molecular orbitals are involved in each excitation. The application of these methods is demonstrated for the electronic absorption spectra of N2 and H2O both in the absence and in the presence of a magnetic field. The dependence of electronic spectra on the magnetic field strength and its orientation relative to the molecule is studied. The complex evolution of the absorption spectra with magnetic field is rationalized using the molecular orbital pair decomposition analysis, which provides crucial insight in strong fields where the spectra are radically different from their zero-field counterparts.

Optimizing Molecular Geometries in Strong Magnetic Fields.

An efficient implementation of geometrical derivatives at the Hartree-Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie-Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Structural and electronic studies of substituted m-terphenyl lithium complexes.

The effect of para-substitution upon the structural and electronic properties of a series of m-terphenyl lithium complexes [R-Ar#-Li]2 (R = t-Bu 1, SiMe32, H 3, Cl 4, CF35; where R-Ar# = 2,6-{2,6-Xyl}2-4-R-C6H2 and 2,6-Xyl = 2,6-Me2C6H3) has been investigated. X-ray crystallography reveals the complexes to be structurally similar, with little variation in C-M-C bond lengths and angles across the series. However, in-depth NMR spectroscopic studies reveal notable electronic differences, showing a linear correlation between the 7Li{1H} NMR chemical shifts of the para-substituted complexes and their Hammett constants. The flanking methyl protons exhibit a similar electronic shift in the 1H NMR spectra, which has been rationalised by the presence of through-space LiH interactions, as evidenced by two-dimensional 7Li-1H heteronuclear Overhauser spectroscopy (HOESY). In both cases, electron-withdrawing substituents are found to cause an upfield peak shift. A computational analysis is employed to account for these trends.

Robust All-Electron Optimization in Orbital-Free Density-Functional Theory Using the Trust-Region Image Method.

We present a Gaussian-basis implementation of orbital-free density-functional theory (OF-DFT) in which the trust-region image method (TRIM) is used for optimization. This second-order optimization scheme has been constructed to provide benchmark all-electron results with very tight convergence of the particle-number constraint, associated chemical potential, and electron density. It is demonstrated that, by preserving the saddle-point nature of the optimization and simultaneously optimizing the density and chemical potential, an order of magnitude reduction in the number of iterations required for convergence is obtained. The approach is compared and contrasted with a new implementation of the nested optimization scheme put forward by Chan, Cohen, and Handy. Our implementation allows for semilocal kinetic-energy (and exchange-correlation) functionals to be handled self-consistently in all-electron calculations. The all-electron Gaussian-basis setting for these calculations will enable direct comparison with a wide range of standard high-accuracy quantum-chemical methods as well as with Kohn-Sham density-functional theory. We expect that the present implementation will provide a useful tool for analyzing the performance of approximate kinetic-energy functionals in finite systems.