Core-electron contributions to the magnetic response of molecules with heavy elements and their significance in aromaticity assessments

This study delves into the magnetic response of core electrons and their influence on the global magnetic response of planar and three-dimensional systems containing heavy elements, employing the removing valence electron (RVE) approximation. We also explore electronic aromaticity indices to understand the potential role of core electrons on electron delocalization in the absence of an external perturbation. The study reveals that core electrons significantly contribute to the overall magnetic response, especially to the magnetic shielding, affecting the interpretation of aromaticity. In contrast, the calculation of the electronic aromaticity indices suggests a negligible participation of the core electrons on electron delocalization. Despite their widespread use, the study emphasizes caution in labeling systems as strongly aromatic based solely on shielding function computations. It is noteworthy to emphasize the limitations associated with each aromaticity criterion; particularly in the context of magnetic shielding function calculations, the core-electron effect contamination is undeniable. Hence, the integration of various criteria becomes imperative for attaining a comprehensive understanding of magnetic responses within complex systems.

Time evolution as an optimization problem: The hydrogen atom in strong laser fields in a basis of time-dependent Gaussian wave packets

Recent advances in attosecond science have made it increasingly important to develop stable, reliable, and accurate algorithms and methods to model the time evolution of atoms and molecules in intense laser fields. A key process in attosecond science is high-harmonic generation, which is challenging to model with fixed Gaussian basis sets, as it produces high-energy electrons, with a resulting rapidly varying and highly oscillatory wave function that extends over dozens of ångström. Recently, Rothe's method, where time evolution is rephrased as an optimization problem, has been applied to the one-dimensional Schrödinger equation. Here, we apply Rothe's method to the hydrogen wave function and demonstrate that thawed, complex-valued Gaussian wave packets with time-dependent width, center, and momentum parameters are able to reproduce spectra obtained from essentially exact grid calculations for high-harmonic generation with only 50-181 Gaussians for field strengths up to 5 × 1014 W/cm2. This paves the way for the inclusion of continuum contributions into real-time, time-dependent electronic-structure theory with Gaussian basis sets for strong fields and eventually accurate simulations of the time evolution of molecules without the Born-Oppenheimer approximation.

Direct Probe of Conical Intersection Photochemistry by Time-Resolved X-ray Magnetic Circular Dichroism

The direct probing of photochemical dynamics by detecting the electronic coherence generated during passage through conical intersections is an intriguing challenge. The weak coherence signal and the difficulty in preparing purely excited wave packets that exclude coherence from other sources make it experimentally challenging. We propose to use time-resolved X-ray magnetic circular dichroism to probe the wave packet dynamics around the conical intersection. The magnetic field amplifies the relative strength of the electronic coherence signal compared to populations through the magnetic field response anisotropy. More importantly, since the excited state relaxation through conical intersections involves a change of parity, the magnetic coupling matches the symmetry of the response function with the electronic coherence, making the coherence signal only sensitive to the conical intersection induced coherence and excludes the pump pulse induced coherence between the ground state and excited state. In this theoretical study, we apply this technique to the photodissociation dynamics of a pyrrole molecule and demonstrate its capability of probing electronic coherence at a conical intersection as well as population transfer. We demonstrate that a magnetic field can be effectively used to extract novel information about electron and nuclear molecular dynamics.

Efficient approximate screening techniques for integrals over London atomic orbitals

Efficient integral screening techniques are essential for the investigation of extended molecular structures. This work presents a critical assessment of well-established approximate screening techniques and extends them for integrals over London atomic orbitals, which are required in the presence of strong, external magnetic fields. Through the examination of helium clusters in such extreme environments, we demonstrate that seemingly straightforward extensions of field-free screening techniques as proposed in the recent literature can lead to significant errors. To rectify this, we propose two alternative screening techniques that lead to the desired speedups while still maintaining strict error control.

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.

The approximate coupled-cluster methods CC2 and CC3 in a finite magnetic field

In this paper, we report on the implementation of CC2 and CC3 in the context of molecules in finite magnetic fields. The methods are applied to the investigation of atoms and molecules through spectroscopic predictions and geometry optimizations for the study of the atmosphere of highly magnetized White Dwarf stars. We show that ground-state finite-field (ff) CC2 is a reasonable alternative to CCSD for energies and, in particular, for geometrical properties. For excited states, ff-CC2 is shown to perform well for states with predominant single-excitation character. Yet, for cases in which the excited state wavefunction has double-excitation character with respect to the reference, ff-CC2 can easily lead to completely unphysical results. Ff-CC3, however, is shown to reproduce the CCSDT behavior very well and enables the treatment of larger systems at a high accuracy.

Insight on Gaussian Basis Set Truncation Errors in Weak to Intermediate Magnetic Fields with an Approximate Hamiltonian

Strong magnetic fields such as those found on white dwarfs have significant effects on the electronic structures of atoms and molecules. However, the vast majority of molecular studies in the literature in such fields are carried out with Gaussian basis sets designed for zero field, leading to large basis set truncation errors [Lehtola et al., Mol. Phys. 2020, 118, e1597989]. In this work, we aim to identify the failures of the Gaussian basis sets in atomic calculations to guide the design of new basis sets for strong magnetic fields. We achieve this by performing fully numerical electronic structure calculations at the complete basis set (CBS) limit for the ground state and low lying excited states of the atoms 1 ≤ Z ≤ 18 in weak to intermediate magnetic fields. We also carry out finite-field calculations for a variety of Gaussian basis sets, introducing a real-orbital approximation for the magnetic-field Hamiltonian. Our primary focus is on the aug-cc-pVTZ basis set, which has been used in many works in the literature. A study of the differences in total energies of the fully numerical CBS limit calculations and the approximate Gaussian basis calculations is carried out to provide insight into basis set truncation errors. Examining a variety of states over the range of magnetic field strengths from B = 0 to B = 0.6B0, we observe significant differences for the aug-cc-pVTZ basis set, while much smaller errors are afforded by the benchmark-quality AHGBSP3-9 basis set [Lehtola, J. Chem. Phys. 2020, 152, 134108]. This suggests that there is considerable room to improve Gaussian basis sets for calculations at finite magnetic fields.

Wandering through quantum-mechanochemistry: from concepts to reactivity and switches

Mechanochemistry has experienced a renaissance in recent years witnessing, at the molecular level, a remarkable interplay between theory and experiment. Molecular mechanochemistry has welcomed a broad spectrum of quantum-chemical methods to evaluate the influence of an external mechanical force on molecular properties. In this contribution, an overview is given on recent work on quantum mechanochemistry in the Brussels Quantum Chemistry group (ALGC). The effect of an external force was scrutinized both in fundamental topics, like reactivity descriptors in Conceptual DFT, and in applied topics, such as designing molecular force probes and tuning the stereoselectivity of certain types of reactions. In the conceptual part, a brief overview of the techniques introducing mechanical forces into a quantum-mechanical description of a molecule is followed by an introduction to conceptual DFT. The evolution of the electronic chemical potential (or electronegativity), chemical hardness and electrophilicity are investigated when a chemical bond in a series of diatomics is put under mechanical stress. Its counterpart, the influence of mechanical stress on bond angles, is analyzed by varying the strain present in alkyne triple bonds by applying a bending force, taking the strain promoted alkyne-azide coupling cycloaddition as an example. The increase of reactivity of the alkyne upon bending is probed by Fukui functions and the local softness. In the applied part, a new molecular force probe is presented based on an intramolecular 6π-electrocyclization in constrained polyenes operating under thermal conditions. A cyclic process is conceived where ring opening and closure are triggered by applying or removing an external pulling force. The efficiency of mechanical activation strongly depends on the magnitude of the applied force and the distance between the pulling points. The idea of pulling point distances as a tool to identify new mechanochemical processes is then tested in [28]hexaphyrins with an intricate equilibrium between Möbius aromatic and Hückel antiaromatic topologies. A mechanical force is shown to trigger the interconversion between the two topologies, using the distance matrix as a guide to select appropriate pulling points. In a final application, the Felkin-Anh model for the addition of nucleophiles to chiral carbonyls under the presence of an external mechanical force is scrutinized. By applying a force for restricting the conformational freedom of the chiral ketone, otherwise inaccessible reaction pathways are promoted on the force-modified potential energy surfaces resulting in a diastereoselectivity different from the force-free reaction.

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.

Unitary coupled-cluster for quantum computation of molecular properties in a strong magnetic field

In truncated coupled-cluster (CC) theories, non-variational and/or generally complex ground-state energies can occur. This is due to the non-Hermitian nature of the similarity transformed Hamiltonian matrix in combination with CC truncation. For chemical problems that deal with real-valued Hamiltonian matrices, complex CC energies rarely occur. However, for complex-valued Hamiltonian matrices, such as those that arise in the presence of strong magnetic fields, complex CC energies can be regularly observed unless certain symmetry conditions are fulfilled. Therefore, in the presence of magnetic fields, it is desirable to pursue CC methods that are guaranteed to give upper-bound, real-valued energies. In this work, we present the first application of unitary CC to chemical systems in a strong magnetic field. This is achieved utilizing the variational quantum eigensolver algorithm applied to the unitary coupled-cluster singles and doubles (UCCSD) method. We benchmark the method on the H2 molecule in a strong magnetic field and then calculate UCCSD energies for the H4 molecule as a function of both geometry and field angle. We show that while standard CCSD can yield generally complex energies that are not an upper-bound to the true energy, UCCSD always results in variational and real-valued energies. We also show that the imaginary components of the CCSD energy are largest in the strongly correlated region. Last, the UCCSD calculations capture a large percentage of the correlation energy.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Eliminating all bonds from the ground state gives rise to ionic bonding in high-spin states of heterodiatomics

Understanding chemical bonding in second-row diatomics has been central to elucidating the basics of bonding itself. Bond strength and the number of bonds are the two factors that decide the reactivity of molecules. While bond strengths have been theoretically computed accurately and experimentally determined, the number of bonds is a more contentious issue, especially for complicated multi-reference systems like C2. We have developed an experimentally verifiable approach to determine bond numbers from excited spin state potential energy surfaces. On applying this to a series of second-row heterodiatomics, we obtain the surprising phenomenon of an inverted charge transfer ionic state after all the ground state bonds are broken via higher spin states. These ionic states are ubiquitous in all heterodiatomics and unexpected in non-metallic systems.

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.

Practical Post-Kohn-Sham Methods for Time-Reversal Symmetry Breaking References

The applicability of reduced scaling algorithms based on auxiliary subspace methods for the correlation energy from the random phase approximation (RPA) as well as the correlation part of the self-energy obtained from the GW method is outlined for time-reversal symmetry breaking Kohn-Sham (KS) references. The updated algorithms allow for an efficient evaluation of RPA energies and GW quasiparticle energies for molecular systems with KS references that break time-reversal symmetry. The latter occur, for example, in magnetic fields. Furthermore, KS references for relativistic open-shell molecules also break time-reversal symmetry due to the single determinant ansatz used. Errors of the updated reduced-scaling algorithms are shown to be negligible compared to reference implementations, while the overall computational scaling is reduced by 2 orders of magnitude. Ionization energies obtained from the GW approximation are shown to be robust even for the electronically complicated group of trivalent lanthanoid ions. Starting from GW quasiparticle energies, it is subsequently shown that light-matter interactions of these systems can be calculated using the Bethe-Salpeter equation (BSE). Using the combined GW-BSE method, the absorption and emission spectra of a molecular europium(III) complex can be obtained including spin-orbit coupling.

Bonding in the high spin lithium clusters: Non-nuclear attractors play a crucial role

The bonding in lithium high-spin clusters contradicts the usual chemical bonding concept since there are no electron pairs between the atoms, and they are bound with parallel spin electrons. Quantum theory of atoms in molecules and interacting quantum atom analysis (IQA) were used to investigate the nature of bonding in the high-spin Li n n + 1 n = 2 - 5 clusters. Our findings demonstrate that the non-nuclear attractors (NNAs) are an essential component of the high-spin lithium clusters and play a key role in keeping them stable. Based on IQA energy terms, an electrostatic destabilizing interaction between the lithium atoms works against the cluster formation. On the other hand, the interactions between lithium atoms and NNA basins are stabilizing and outweigh the lithium-lithium destabilizing effects. In fact, NNAs tend to draw lithium atoms together and stabilize the resulting cluster. The high-spin clusters of lithium can be regarded as electrostatically driven compounds since the electrostatic components are primarily responsible for the stabilizing interactions between NNAs and Li atoms. The only exception is ³ Li₂ , which lacks NNA and has a non-repellent lithium-lithium interaction. Indeed, in the ³ Li₂ , the interatomic electrostatic component is negligibly small, and the exchange-correlation term leads to a weak bonding interaction.

Molecular vibrations in the presence of velocity-dependent forces

A semiclassical theory of small oscillations is developed for nuclei that are subject to velocity-dependent forces in addition to the usual interatomic forces. When the velocity-dependent forces are due to a strong magnetic field, novel effects arise-for example, the coupling of vibrational, rotational, and translational modes. The theory is first developed using Newtonian mechanics and we provide a simple quantification of the coupling between these types of modes. We also discuss the mathematical structure of the problem, which turns out to be a quadratic eigenvalue problem rather than a standard eigenvalue problem. The theory is then re-derived using the Hamiltonian formalism, which brings additional insight, including a close analogy to the quantum-mechanical treatment of the problem. Finally, we provide numerical examples for the H₂, HT, and HCN molecules in a strong magnetic field.

Time-dependent nuclear-electronic orbital Hartree-Fock theory in a strong uniform magnetic field

In an ultrastrong magnetic field, with field strength B ≈ B₀ = 2.35 × 10⁵ T, molecular structure and dynamics differ strongly from that observed on the Earth. Within the Born-Oppenheimer (BO) approximation, for example, frequent (near) crossings of electronic energy surfaces are induced by the field, suggesting that nonadiabatic phenomena and processes may play a more important role in this mixed-field regime than in the weak-field regime on Earth. To understand the chemistry in the mixed regime, it therefore becomes important to explore non-BO methods. In this work, the nuclear-electronic orbital (NEO) method is employed to study protonic vibrational excitation energies in the presence of a strong magnetic field. The NEO generalized Hartree-Fock theory and time-dependent Hartree-Fock (TDHF) theory are derived and implemented, accounting for all terms that result as a consequence of the nonperturbative treatment of molecular systems in a magnetic field. The NEO results for HCN and FHF^- with clamped heavy nuclei are compared against the quadratic eigenvalue problem. Each molecule has three semi-classical modes owing to the hydrogen-two precession modes that are degenerate in the absence of a field and one stretching mode. The NEO-TDHF model is found to perform well; in particular, it automatically captures the screening effects of the electrons on the nuclei, which are quantified through the difference in energy of the precession modes.

External fields in conceptual density functional theory

The necessity of the recent incorporation of new external variables in the context of conceptual DFT (CDFT) is discussed based on the ever-increasing portfolio of experimental reaction conditions in the endeavor of experimentalists to synthesize new molecules with unprecedented properties. Electric and magnetic fields (ε and B), mechanical forces (F), and confinement are proposed as valuable new variables, extending conventional CDFT and its associated response functions. A finite field approach is used to calculate the evolution of both global and local descriptors in a selected series of atomic and molecular applications, and from it derive new response function involving, with one exception, the first derivative to the field considered. The electric field results, displaying, for example, a case of a field-induced enantioselectivity in the Fukui function, may be instrumental in the recent upsurge of chemistry in oriented external electric fields. The study of atomic electronegativity and hardness in magnetic fields displays a piecewise behavior, associated to configurational jumps upon increasing field strength and reveals an overall compression of their ranges for stronger fields, which may be guiding upon investigating chemistry in extremely high fields like in white dwarfs. The evolution of the electronegativity and hardness of diatomics under mechanical force can elegantly be traced back to differences in their equilibrium distance in the neutral, cationic, and anionic state. The well-known reduction of the polarizability under confinement can be seen as a fore-runner of the increasing hardness of atoms under pressure, presently under investigation. Periodicity showing up in a spontaneous way in the variety of properties is a leitmotiv in this study, as well as the interconnections/analogies between the different response functions.

Berry Population Analysis: Atomic Charges from the Berry Curvature in a Magnetic Field

The Berry curvature is essential in Born-Oppenheimer molecular dynamics, describing the screening of the nuclei by the electrons in a magnetic field. Parts of the Berry curvature can be understood as the external magnetic field multiplied by an effective charge so that the resulting Berry force behaves like a Lorentz force during the simulations. Here, we investigate whether these effective charges can provide insight into the electronic structure of a given molecule or, in other words, whether we can perform a population analysis based on the Berry curvature. To develop our approach, we first rewrite the Berry curvature in terms of charges that partially capture the effective charges and their dependencies on the nuclear velocities. With these Berry charges and charge fluctuations, we then construct our population analysis yielding atomic charges and overlap populations. Calculations at the Hartree-Fock level reveal that the atomic charges are similar to those obtained from atomic polar tensors. However, since we additionally obtain an estimate for the fluctuations of the charges and a partitioning of the atomic charges into contributions from all atoms, we conclude that the Berry population analysis is a useful alternative tool to analyze the electronic structures of molecules.

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.

Magnetic-translational sum rule and approximate models of the molecular Berry curvature

The Berry connection and curvature are key components of electronic structure calculations for atoms and molecules in magnetic fields. They ensure the correct translational behavior of the effective nuclear Hamiltonian and the correct center-of-mass motion during molecular dynamics in these environments. In this work, we demonstrate how these properties of the Berry connection and curvature arise from the translational symmetry of the electronic wave function and how they are fully captured by a finite basis set of London orbitals but not by standard Gaussian basis sets. This is illustrated by a series of Hartree-Fock calculations on small molecules in different basis sets. Based on the resulting physical interpretation of the Berry curvature as the shielding of the nuclei by the electrons, we introduce and test a series of approximations using the Mulliken fragmentation scheme of the electron density. These approximations will be particularly useful in ab initio molecular dynamics calculations in a magnetic field since they reduce the computational cost, while recovering the correct physics and up to 95% of the exact Berry curvature.

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.

Molecular Dynamics of Linear Molecules in Strong Magnetic Fields

Molecular rotations and vibrations have been extensively studied by chemists for decades, both experimentally using spectroscopic methods and theoretically with the help of quantum chemistry. However, the theoretical investigation of molecular rotations and vibrations in strong magnetic fields requires computationally more demanding tools. As such, proper calculations of rotational and vibrational spectra were not feasible up until very recently. In this work, we present rotational and vibrational spectra for two small linear molecules, H2 and LiH, in strong magnetic fields. By treating the nuclei as classical particles, trajectories for rotations and vibrations are simulated from ab initio molecular dynamics. Born-Oppenheimer potential energy surfaces are calculated at the Hartree-Fock and MP2 levels of theory, using London atomic orbitals to ensure gauge origin invariance. For the calculation of nuclear trajectories, a highly efficient Tajima propagator is introduced, incorporating the Berry curvature tensor accounting for the screening of nuclear charges.

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 He n and CH n 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.

Analyzing the Behavior of Spin Phases in External Magnetic Fields by Means of Spin-Constrained States

During molecular dissociation in the presence of an external uniform magnetic field, electrons flip their spin antiparallel to the magnetic field because of the stabilizing influence of the spin Zeeman operator. Although generalized Hartree-Fock descriptions furnish the optimal mean-field energetic description of such bond-breaking processes, they are allowed to break Ŝ_z symmetry, leading to intricate and unexpected spin phases and phase transitions. In this work, we show that the behavior of these molecular spin phases can be interpreted in terms of spin phase diagrams constructed by constraining states to target expectation values of projected spin. The underlying constrained states offer a complete electronic characterization of the spin phases and spin phase transitions, as they can be analyzed using standard quantum chemical tools. Because the constrained states effectively span the entire phase space, they could provide an excellent starting point for post-Hartree-Fock methods aimed at gaining more electron correlation or regaining spin symmetry.

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 GW/BSE Method in Magnetic Fields

The GW approximation and the Bethe–Salpeter equation have been implemented into the Turbomole program package for computations of molecular systems in a strong, finite magnetic field. Complex-valued London orbitals are used as basis functions to ensure gauge-invariant computational results. The implementation has been benchmarked against triplet excitation energies of 36 small to medium-sized molecules against reference values obtained at the approximate coupled-cluster level (CC2 approximation). Finally, a spectacular change of colour from orange to green of the tetracene molecule is induced by applying magnetic fields between 0 and 9,000 T perpendicular to the molecular plane.

Molecular point groups and symmetry in external magnetic fields

As quantum-chemical calculations of molecules in static external magnetic fields are becoming increasingly popular, the description of molecular symmetry under such conditions is also becoming more and more relevant. Using group theory, a general scheme of identifying the molecular point group in an external magnetic field is constructed. For both point groups that are non-existent in the absence of a field (C_∞ and C_{∞ h}) and their double groups, the character tables are presented. General properties of all possible point groups are discussed, and it is mathematically proven that they are all Abelian.

Current density and molecular magnetic properties

We give an overview of the molecular response to an external magnetic field perturbing quantum mechanical systems. We present state-of-the-art methods for calculating magnetically-induced current-density susceptibilities. We discuss the essence and properties of current-density susceptibilities and how molecular magnetic properties can be calculated from them. We also review the theory of spin-current densities, how relativity affects current densities and magnetic properties. An overview of the magnetic ring-current criterion for aromaticity is given, which has implications on theoretical and experimental research. The recently reported theory of antiaromaticity and how molecular symmetry affects the magnetic response are discussed and applied to closed-shell paramagnetic molecules. The topology of magnetically induced current densities and its consequences for molecular magnetic properties are also presented with twisted and toroidal molecules as examples.

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, H₂, 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.

Ab Initio molecular dynamics with screened Lorentz forces. II. Efficient propagators and rovibrational spectra in strong magnetic fields

Strong magnetic fields have a large impact on the dynamics of molecules. In addition to the changes in the electronic structure, the nuclei are exposed to the Lorentz force with the magnetic field being screened by the electrons. In this work, we explore these effects using ab initio molecular dynamics simulations based on an effective Hamiltonian calculated at the Hartree-Fock level of theory. To correctly include these non-conservative forces in the dynamics, we have designed a series of novel propagators that show both good efficiency and stability in test cases. As a first application, we analyze simulations of He and H₂ at two field strengths characteristic of magnetic white dwarfs (0.1 B₀ = 2.35 × 10⁴ T and B₀ = 2.35 × 10⁵ T). While the He simulations clearly demonstrate the importance of electron screening of the Lorentz force in the dynamics, the extracted rovibrational spectra of H₂ reveal a number of fascinating features not observed in the field-free case: couplings of rotations/vibrations with the cyclotron rotation, overtones with unusual selection rules, and hindered rotations that transmute into librations with increasing field strength. We conclude that our presented framework is a powerful tool to investigate molecules in these extreme environments.

Ab initio molecular dynamics with screened Lorentz forces. I. Calculation and atomic charge interpretation of Berry curvature

The dynamics of a molecule in a magnetic field is significantly different from its zero-field counterpart. One important difference in the presence of a field is the Lorentz force acting on the nuclei, which can be decomposed as the sum of the bare nuclear Lorentz force and a screening force due to the electrons. This screening force is calculated from the Berry curvature and can change the dynamics qualitatively. It is therefore important to include the contributions from the Berry curvature in molecular dynamics simulations in a magnetic field. In this work, we present a scheme for calculating the Berry curvature numerically using a finite-difference technique, addressing challenges related to the arbitrary global phase of the wave function. The Berry curvature is calculated as a function of bond distance for H₂ at the restricted and unrestricted Hartree-Fock levels of theory and for CH⁺ as a function of the magnetic field strength at the restricted Hartree-Fock level of theory. The calculations are carried out using basis sets of contracted Gaussian functions equipped with London phase factors (London orbitals) to ensure gauge-origin invariance. In this paper, we also interpret the Berry curvature in terms of atomic charges and discuss its convergence in basis sets with and without London phase factors. The calculation of the Berry curvature allows for its inclusion in ab initio molecular dynamics simulations in a magnetic field.

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.

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 N₂ and H₂O 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.

Benchmarking Density Functional Approximations for Diamagnetic and Paramagnetic Molecules in Nonuniform Magnetic Fields

In this article, correlated studies on a test set of 36 small molecules are carried out with both wavefunction (HF, MP2, CCSD) and density functional (LDA, KT3, cTPSS, cM06-L) methods. The effect of correlation on exotic response properties such as molecular electronic anapole susceptibilities is studied and the performance of the various density functional approximations are benchmarked against CCSD and/or MP2. Atoms and molecules are traditionally classified into "diamagnetic" and "paramagnetic" based on their isotropic response to uniform magnetic fields. However, in this article, we propose a more fine-grained classification of molecular systems on the basis of their response to generally nonuniform magnetic fields. The relation of orientation to different qualitative responses is also considered.

Comprehending the quadruple bonding conundrum in C2 from excited state potential energy curves

The question of quadruple bonding in C₂ has emerged as a hot button issue, with opinions sharply divided between the practitioners of Valence Bond (VB) and Molecular Orbital (MO) theory.

The question of quadruple bonding in C₂ has emerged as a hot button issue, with opinions sharply divided between the practitioners of Valence Bond (VB) and Molecular Orbital (MO) theory. Here, we have systematically studied the Potential Energy Curves (PECs) of low lying high spin sigma states of C₂, N₂, Be₂ and HC Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> CH using several MO based techniques such as CASSCF, RASSCF and MRCI. The analyses of the PECs for the ^2S+1Σ_g/u (with 2S + 1 = 1, 3, 5, 7, 9) states of C₂ and comparisons with those of relevant dimers and the respective wavefunctions were conducted. We contend that unlike in the case of N₂ and HC Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> CH, the presence of a deep minimum in the ⁷Σ⁺ state of C₂ and CN⁺ suggests a latent quadruple bonding nature in these two dimers. Our investigations reveal that the number of bonds in the ground state can be determined for 2^nd row dimers by figuring out at what value of spin symmetry a purely dissociative PEC is obtained. For N₂ and HC Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> CH the purely dissociative PEC appears for the septet spin symmetry as compared to that for the nonet in C₂. This is indicative of a higher number of bonds between the two 2^nd row atoms in C₂ as compared to those of N₂ and HC Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> CH. Hence, we have struck a reconciliatory note between the MO and VB approaches. The evidence provided by us can be experimentally verified, thus providing the window so that the narrative can move beyond theoretical conjectures.

Calculation of Diamagnetic Susceptibility Tensors of Organic Crystals: From Coronene to Pharmaceutical Polymorphs

Understanding why crystallization in strong magnetic fields can lead to new polymorphs requires methods to calculate the diamagnetic response of organic molecular crystals. We develop the calculation of the macroscopic diamagnetic susceptibility tensor, χ^cryst, for organic molecular crystals using periodic density functional methods. The crystal magnetic susceptibility tensor, χ^cryst, for all experimentally known polymorphs, and its molecular counterpart, χ^mol, are calculated for flexible pharmaceuticals such as carbamazepine, flufenamic acid, and chalcones, and rigid molecules, such as benzene, pyridine, acridine, anthracene, and coronene, whose molecular magnetic properties have been traditionally studied. A tensor addition method is developed to approximate the crystal diamagnetic susceptibility tensor, χ^cryst, from the molecular one, χ^mol, giving good agreement with those calculated directly using the more costly periodic density functional method for χ^cryst. The response of pharmaceutical molecules and crystals to magnetic fields, as embodied by χ^cryst, is largely determined by the packing in the crystal, as well as the molecular conformation. The anisotropy of χ^cryst can vary considerably between polymorphs though the isotropic terms are fairly constant. The implications for developing a computational method for predicting whether crystallization in a magnetic field could produce a novel or different polymorph are discussed.

Coherent spin control of s-, p-, d- and f-electrons in a silicon quantum dot

Once the periodic properties of elements were unveiled, chemical behaviour could be understood in terms of the valence of atoms. Ideally, this rationale would extend to quantum dots, and quantum computation could be performed by merely controlling the outer-shell electrons of dot-based qubits. Imperfections in semiconductor materials disrupt this analogy, so real devices seldom display a systematic many-electron arrangement. We demonstrate here an electrostatically confined quantum dot that reveals a well defined shell structure. We observe four shells (31 electrons) with multiplicities given by spin and valley degrees of freedom. Various fillings containing a single valence electron-namely 1, 5, 13 and 25 electrons-are found to be potential qubits. An integrated micromagnet allows us to perform electrically-driven spin resonance (EDSR), leading to faster Rabi rotations and higher fidelity single qubit gates at higher shell states. We investigate the impact of orbital excitations on single qubits as a function of the dot deformation and exploit it for faster qubit control.