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.

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.

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.

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.

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, ||, is varied between 0.

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

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.

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.

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.

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.

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] (R = t-Bu 1, SiMe2, H 3, Cl 4, CF5; where R-Ar = 2,6-{2,6-Xyl}-4-R-CH and 2,6-Xyl = 2,6-MeCH) 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 Li{H} NMR chemical shifts of the para-substituted complexes and their Hammett constants.

Analyzing Magnetically Induced Currents in Molecular Systems Using Current-Density-Functional Theory.

A suite of tools for the analysis of magnetically induced currents is introduced. These are applicable to both the weak-field regime, well described by linear response perturbation theory, and to the strong-field regime, which is inaccessible to such methods. A disc-based quadrature scheme is proposed for the analysis of magnetically induced current susceptibilities, providing quadratures that are consistently defined between different molecular systems and applicable to both planar 2D and general 3D molecular systems in a black-box manner.

GW quasiparticle energies of atoms in strong magnetic fields.

Quasiparticle energies of the atoms H-Ne have been computed in the GW approximation in the presence of strong magnetic fields with field strengths varying from 0 to 0.25 atomic units (0.25 B=0.

Uniform magnetic fields in density-functional theory.

We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional density functional theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term linear vector potential-DFT (LDFT), the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre-Fenchel transformations are constructed.

Connections between variation principles at the interface of wave-function and density-functional theories.

A recently proposed variation principle [N. I. Gidopoulos, Phys.

Magnetic-Field Density-Functional Theory (BDFT): Lessons from the Adiabatic Connection.

We study the effects of magnetic fields in the context of magnetic field density-functional theory (BDFT), where the energy is a functional of the electron density ρ and the magnetic field B. We show that this approach is a worthwhile alternative to current-density functional theory (CDFT) and may provide a viable route to the study of many magnetic phenomena using density-functional theory (DFT). The relationship between BDFT and CDFT is developed and clarified within the framework of the four-way correspondence of saddle functions and their convex and concave parents in convex analysis.

Efficient Calculation of Molecular Integrals over London Atomic Orbitals.

The use of London atomic orbitals (LAOs) in a nonperturbative manner enables the determination of gauge-origin invariant energies and properties for molecular species in arbitrarily strong magnetic fields. Central to the efficient implementation of such calculations for molecular systems is the evaluation of molecular integrals, particularly the electron repulsion integrals (ERIs). We present an implementation of several different algorithms for the evaluation of ERIs over Gaussian-type LAOs at arbitrary magnetic field strengths.

Interpolated energy densities, correlation indicators and lower bounds from approximations to the strong coupling limit of DFT.

We investigate the construction of approximated exchange-correlation functionals by interpolating locally along the adiabatic connection between the weak- and the strong-coupling regimes, focussing on the effect of using approximate functionals for the strong-coupling energy densities. The gauge problem is avoided by dealing with quantities that are all locally defined in the same way. Using exact ingredients at weak coupling we are able to isolate the error coming from the approximations at strong coupling only.

Exchange-Correlation Functionals via Local Interpolation along the Adiabatic Connection.

The construction of density-functional approximations is explored by modeling the adiabatic connection locally, using energy densities defined in terms of the electrostatic potential of the exchange-correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to assess the accuracy of a range of local interpolation models against accurate exchange-correlation energy densities.

Current Density Functional Theory Using Meta-Generalized Gradient Exchange-Correlation Functionals.

We present the self-consistent implementation of current-dependent (hybrid) meta-generalized gradient approximation (mGGA) density functionals using London atomic orbitals. A previously proposed generalized kinetic energy density is utilized to implement mGGAs in the framework of Kohn-Sham current density functional theory (KS-CDFT). A unique feature of the nonperturbative implementation of these functionals is the ability to seamlessly explore a wide range of magnetic fields up to 1 au (∼235 kT) in strength.

Fractional Electron Loss in Approximate DFT and Hartree-Fock Theory.

Plots of electronic energy vs electron number, determined using approximate density functional theory (DFT) and Hartree-Fock theory, are typically piecewise convex and piecewise concave, respectively. The curves also commonly exhibit a minimum and maximum, respectively, in the neutral → anion segment, which lead to positive DFT anion HOMO energies and positive Hartree-Fock neutral LUMO energies. These minima/maxima are a consequence of using basis sets that are local to the system, preventing fractional electron loss.