The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion-water interactions, typified by F(HO) with ≳ 15.
View Article and Find Full Text PDFJ Chem Theory Comput
November 2024
Phys Chem Chem Phys
September 2024
Correction for 'Visualizing and characterizing excited states from time-dependent density functional theory' by John M. Herbert , 2024, , 3755-3794, https://doi.org/10.
View Article and Find Full Text PDFThe titular domain-based local pair natural orbital (DLPNO) approximation is the most widely used method for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies in large supramolecular complexes has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density and involve nonlocal dispersion that is discarded or approximated if the screening of pair natural orbitals (PNOs) is too aggressive. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases.
View Article and Find Full Text PDFTime-dependent density functional theory (TD-DFT) within a restricted excitation space is an efficient means to compute core-level excitation energies using only a small subset of the occupied orbitals. However, core-to-valence excitation energies are significantly underestimated when standard exchange-correlation functionals are used, which is partly traceable to systemic issues with TD-DFT's description of Rydberg and charge-transfer excited states. To mitigate this, we have implemented an empirically modified combination of configuration interaction with single substitutions (CIS) based on Kohn-Sham orbitals, which is known as "DFT/CIS.
View Article and Find Full Text PDFCharge-transfer (CT) excited states are crucial to organic light-emitting diodes (OLEDs), particularly to those based on thermally activated delayed fluorescence (TADF). However, accurately modeling CT states remains challenging, even with modern implementations of (time-dependent) density functional theory [(TD-)DFT], especially in a dielectric environment. To identify shortcomings and improve the methodology, we previously established the STGABS27 benchmark set with highly accurate experimental references for the adiabatic energy gap between the lowest singlet and triplet excited states (Δ).
View Article and Find Full Text PDFIn constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only.
View Article and Find Full Text PDFThe effects of temperature and chemical environment on a pentanuclear cyanide-bridged, trigonal bipyramidal molecular paramagnet have been investigated. Using element- and oxidation state-specific near-ambient pressure X-ray photoemission spectroscopy (NAP-XPS) to probe charge transfer and second order, nonlinear vibrational spectroscopy, which is sensitive to symmetry changes based on charge (de)localization coupled with DFT, a detailed picture of environmental effects on charge-transfer-induced spin transitions is presented. The molecular cluster, CoFe(tmphen)(μ-CN)(-CN), abbrev.
View Article and Find Full Text PDFPhys Chem Chem Phys
January 2024
Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization.
View Article and Find Full Text PDFThe hydrated electron, e, has attracted much attention as a central species in radiation chemistry. However, much less is known about e at the water/air surface, despite its fundamental role in electron transfer processes at interfaces. Using time-resolved electronic sum-frequency generation spectroscopy, the electronic spectrum of e at the water/air interface and its dynamics are measured here, following photo-oxidation of the phenoxide anion.
View Article and Find Full Text PDFJ Phys Chem Lett
December 2023
Donor-acceptor (D-A) materials can exhibit a wide range of unique photophysical properties with applications in next-generation optoelectronics. Electronic structure calculations of D-A dimers are often employed to predict the properties of D-A materials. One of the most important D-A dimer quantities is the degree of charge transfer (DCT) in the S state, which correlates with properties such as fluorescence lifetimes and intersystem crossing rates in D-A materials.
View Article and Find Full Text PDFThe many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening.
View Article and Find Full Text PDFWe simulate X-ray absorption spectra at elemental K-edges using time-dependent density functional theory (TDDFT) in both its conventional linear-response implementation and its explicitly time-dependent or "real-time" formulation. Real-time TDDFT simulations enable broadband spectra calculations without the need to invoke frozen occupied orbitals ("core/valence separation"), but we find that these spectra are often contaminated by transitions to the continuum that originate from lower-energy core and semicore orbitals. This problem becomes acute in triple-ζ basis sets, although it is sometimes sidestepped in double-ζ basis sets.
View Article and Find Full Text PDFMethods for computing X-ray absorption spectra based on a constrained core hole (possibly containing a fractional electron) are examined. These methods are based on Slater's transition concept and its generalizations, wherein core-to-valence excitation energies are determined using Kohn-Sham orbital energies. Methods examined here avoid promoting electrons beyond the lowest unoccupied molecular orbital, facilitating robust convergence.
View Article and Find Full Text PDFElectronic structure calculations on enzymes require hundreds of atoms to obtain converged results, but fragment-based approximations offer a cost-effective solution. We present calculations on enzyme models containing 500-600 atoms using the many-body expansion, comparing to benchmarks in which the entire enzyme-substrate complex is described at the same level of density functional theory. When the amino acid fragments contain ionic side chains, the many-body expansion oscillates under vacuum boundary conditions but rapid convergence is restored using low-dielectric boundary conditions.
View Article and Find Full Text PDFJ Chem Theory Comput
November 2022
A widespread belief persists that the Boys-Bernardi function counterpoise (CP) procedure "overcorrects" supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic.
View Article and Find Full Text PDFWith the large numbers of man-made chemicals produced and released in the environment, there is a need to provide assessments on their potential effects on environmental safety and human health. Current regulatory frameworks rely on a mix of both hazard and risk-based approaches to make safety decisions, but the large number of chemicals in commerce combined with an increased need to conduct assessments in the absence of animal testing makes this increasingly challenging. This challenge is catalysing the use of more mechanistic knowledge in safety assessment from both in silico and in vitro approaches in the hope that this will increase confidence in being able to identify modes of action (MoA) for the chemicals in question.
View Article and Find Full Text PDFModeling polaron defects is an important aspect of computational materials science, but the description of unpaired spins in density functional theory (DFT) often suffers from delocalization error. To diagnose and correct the overdelocalization of spin defects, we report an implementation of density-corrected (DC-)DFT and its analytic energy gradient. In DC-DFT, an exchange-correlation functional is evaluated using a Hartree-Fock density, thus incorporating electron correlation while avoiding self-interaction error.
View Article and Find Full Text PDF