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

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

We developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants

We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S₁ → S₀ transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 10⁶-10¹⁶ s^-1 compared with an experimental rate of the order 10⁷ s^-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.

Communication: Impact of the current density on paramagnetic NMR properties

Meta-generalized gradient approximations (meta-GGAs) and local hybrid functionals generally depend on the kinetic energy density τ. For magnetic properties, this necessitates generalizations to ensure gauge invariance. In most implementations, τ is generalized by incorporating the external magnetic field. However, this introduces artifacts in the response of the density matrix and does not satisfy the iso-orbital constraint. Here, we extend previous approaches based on the current density to paramagnetic NMR shieldings and EPR g-tensors. The impact is assessed for main-group compounds and transition-metal complexes considering 25 density functional approximations. It is shown that the current density leads to substantial improvements-especially for the popular Minnesota and SCAN functional families. Thus, we strongly recommend to use the current density generalized τ in paramagnetic NMR and EPR calculations with meta-GGAs.

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.

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.

Assessment of the performance of DFT functionals in the fulfillment of off-diagonal hypervirial relationships

Off-diagonal hypervirial relationships, combined with quantum mechanical sum rules of charge-current conservation, offer a way to test electronic excited-state transition energies and moments, which does not need any external reference. A number of fundamental relationships were recast into absolute deviations from zero, which have been used to assess the performance of some popular DFT functionals. Extended TD-DFT calculations have been carried out for a pool of molecules chosen for this purpose, adopting a large basis set to ensure high quality results. A partial agreement with previous benchmarks is observed.