DFT and TD-DFT studies to elucidate the configurational isomers of ferric aerobactin, ferric petrobactin, and their ferric photoproducts.

Iron-chelating siderophores such as aerobactin and petrobactin are produced by marine bacteria to sequester iron under low iron stress. Those that contain a citrate moiety undergo light-catalyzed ligand-to-metal charge transfer, inducing decarboxylation and formation of photoproducts. In this work, we employed density functional theory to obtain the optimized geometries and determine the relative energies and geometric parameters of different configurations of Fe(III)-coordinated aerobactin, petrobactin, and their photoproducts. Time-dependent density functional theory was then used to compute the UV-Vis absorption spectra of these species, and the comparison against experimental spectra further elucidated the structural configurations most likely to be adopted by these compounds. Frequency calculations provided Fe-O force constants on the same order as other siderophores. The relative energies and predicted spectra support the cis-cis C-fac configuration for ferric aerobactin and the cis-trans C-mer configuration for its photoproduct, while only mild support is found for specific configurations of the ferric petrobactin structures (meta-mer and meta-fac for the precursor, cis-cis para-fac for the photoproduct). The predicted ferric petrobactin spectra are found to be fairly insensitive to the configuration of the ferric complex.

Environment- and Conformation-Induced Frequency Shifts of C-D Vibrational Stark Probes in NAD(P)H Cofactors.

NAD(P)H cofactors are found in all forms of life and are essential for electron and hydrogen atom transfer. The linear response of a carbon-deuterium (C-D) vibration based on the vibrational Stark effect can facilitate measurements of electric fields for critical biological reactions including cofactor-mediated hydride transfer. We find both inter- and intramolecular electric fields influence the C-D frequency in NAD(P)H and nicotinamide-like models where the reactive C4-hydrogen has been deuterated. Hence, the C-D frequency can report both environmental electrostatics and conformational changes of the nicotinamide ring. Conformation-dependent effects are mediated through space as electrostatic effects, rather than through-bond. A Stark tuning rate of ∼0.57 cm-1/(MV/cm) was determined using both experimental and computational approaches, including vibrational solvatochromism, molecular dynamics simulations, and in silico Stark calculations. The vibrational probe's Stark tuning rate is shown to be robust and suitable for measuring fields along hydride transfer reaction coordinates in enzymes.

Environmentally Driven Symmetry Breaking Quenches Dual Fluorescence in Proflavine.

Nonadiabatic couplings between several electronic excited states are ubiquitous in many organic chromophores and can significantly influence optical properties. A recent experimental study demonstrated that the proflavine molecule exhibits surprising dual fluorescence in the gas phase, which is suppressed in polar solvent environments. Here, we uncover the origin of this phenomenon by parametrizing a linear-vibronic coupling Hamiltonian from spectral densities of system-bath coupling constructed along molecular dynamics trajectories, fully accounting for interactions with the condensed-phase environment. The finite-temperature absorption, steady-state emission, and time-resolved emission spectra are then computed using powerful, numerically exact tensor network approaches. We find that the dual fluorescence in vacuum is driven by a single well-defined coupling mode but is quenched in solution due to dynamic solvent-driven symmetry breaking that mixes the two low-lying electronic states. We expect the computational framework developed here to be widely applicable to the study of non-Condon effects in complex condensed-phase environments.

Training machine learning potentials for reactive systems: A Colab tutorial on basic models.

In the last several years, there has been a surge in the development of machine learning potential (MLP) models for describing molecular systems. We are interested in a particular area of this field - the training of system-specific MLPs for reactive systems - with the goal of using these MLPs to accelerate free energy simulations of chemical and enzyme reactions. To help new members in our labs become familiar with the basic techniques, we have put together a self-guided Colab tutorial (https://cc-ats.github.io/mlp_tutorial/), which we expect to be also useful to other young researchers in the community. Our tutorial begins with the introduction of simple feedforward neural network (FNN) and kernel-based (using Gaussian process regression, GPR) models by fitting the two-dimensional Müller-Brown potential. Subsequently, two simple descriptors are presented for extracting features of molecular systems: symmetry functions (including the ANI variant) and embedding neural networks (such as DeepPot-SE). Lastly, these features will be fed into FNN and GPR models to reproduce the energies and forces for the molecular configurations in a Claisen rearrangement reaction.

Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions.

In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.

Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure.

Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. We show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions.

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Machine learning based implicit solvent model for aqueous-solution alanine dipeptide molecular dynamics simulations.

Inspired by the recent work from Noé and coworkers on the development of machine learning based implicit solvent model for the simulation of solvated peptides [Chen et al., J. Chem. Phys., 2021, 155, 084101], here we report another investigation of the possibility of using machine learning (ML) techniques to "derive" an implicit solvent model directly from explicit solvent molecular dynamics (MD) simulations. For alanine dipeptide, a machine learning potential (MLP) based on the DeepPot-SE representation of the molecule was trained to capture its interactions with its average solvent environment configuration (ASEC). The predicted forces on the solute deviated only by an RMSD of 0.4 kcal mol-1 Å-1 from the reference values, and the MLP-based free energy surface differed from that obtained from explicit solvent MD simulations by an RMSD of less than 0.9 kcal mol-1. Our MLP training protocol could also accurately reproduce combined quantum mechanical molecular mechanical (QM/MM) forces on the quantum mechanical (QM) solute in ASEC environment, thus enabling the development of accurate ML-based implicit solvent models for ab initio-QM MD simulations. Such ML-based implicit solvent models for QM calculations are cost-effective in both the training stage, where the use of ASEC reduces the number of data points to be labelled, and the inference stage, where the MLP can be evaluated at a relatively small additional cost on top of the QM calculation of the solute.

SPICE, A Dataset of Drug-like Molecules and Peptides for Training Machine Learning Potentials.

Machine learning potentials are an important tool for molecular simulation, but their development is held back by a shortage of high quality datasets to train them on. We describe the SPICE dataset, a new quantum chemistry dataset for training potentials relevant to simulating drug-like small molecules interacting with proteins. It contains over 1.1 million conformations for a diverse set of small molecules, dimers, dipeptides, and solvated amino acids. It includes 15 elements, charged and uncharged molecules, and a wide range of covalent and non-covalent interactions. It provides both forces and energies calculated at the ωB97M-D3(BJ)/def2-TZVPPD level of theory, along with other useful quantities such as multipole moments and bond orders. We train a set of machine learning potentials on it and demonstrate that they can achieve chemical accuracy across a broad region of chemical space. It can serve as a valuable resource for the creation of transferable, ready to use potential functions for use in molecular simulations.

Analytic high-order energy derivatives for metal nanoparticle-mediated infrared and Raman scattering spectra within the framework of quantum mechanics/molecular mechanics model with induced charges and dipoles.

This work is devoted to deriving and implementing analytic second- and third-order energy derivatives with respect to the nuclear coordinates and external electric field within the framework of the hybrid quantum mechanics/molecular mechanics method with induced charges and dipoles (QM/DIM). Using these analytic energy derivatives, one can efficiently compute the harmonic vibrational frequencies, infrared (IR) and Raman scattering (RS) spectra of the molecule in the proximity of noble metal clusters/nanoparticles. The validity and accuracy of these analytic implementations are demonstrated by the comparison of results obtained by the finite-difference method and the analytic approaches and by the full QM and QM/DIM calculations. The complexes formed by pyridine and two sizes of gold clusters (Au18 and Au32) at varying intersystem distances of 3, 4, and 5 Å are used as the test systems, and Raman spectra of 4,4'-bipyridine in the proximity of Au2057 and Ag2057 metal nanoparticles (MNP) are calculated by the QM/DIM method and compared with experimental results as well. We find that the QM/DIM model can well reproduce the IR spectra obtained from full QM calculations for all the configurations, while although it properly enhances some of the vibrational modes, it artificially overestimates RS spectral intensities of several modes for the systems with very short intersystem distance. We show that this could be improved, however, by incorporating the hyperpolarizability of the gold metal cluster in the evaluation of RS intensities. Additionally, we address the potential impact of charge migration between the adsorbate and MNPs.

Computational Investigation of Substituent Effects on the Fluorescence Wavelengths of Oxyluciferin Analogs.

Oxyluciferin, which is the light emitter for firefly bioluminescence, has been subjected to extensive chemical modifications to tune its emission wavelength and quantum yield. However, the exact mechanisms for various electron-donating and withdrawing groups to perturb the photophysical properties of oxyluciferin analogs are still not fully understood. To elucidate the substituent effects on the fluorescence wavelength of oxyluciferin analogs, we applied the absolutely localized molecular orbitals (ALMO)-based frontier orbital analysis to assess various types of interactions (i.e. permanent electrostatics/exchange repulsion, polarization, occupied-occupied orbital mixing, virtual-virtual orbital mixing, and charge-transfer) between the oxyluciferin and substituent orbitals. We suggested two distinct mechanisms that can lead to red-shifted oxyluciferin emission wavelength, a design objective that can help increase the tissue penetration of bioluminescence emission. Within the first mechanism, an electron-donating group (such as an amino or dimethylamino group) can contribute its highest occupied molecular orbital (HOMO) to an out-of-phase combination with oxyluciferin's HOMO, thus raising the HOMO energy of the substituted analog and narrowing its HOMO-LUMO gap. Alternatively, an electron-withdrawing group (such as a nitro or cyano group) can participate in an in-phase virtual-virtual orbital mixing of fragment LUMOs, thus lowering the LUMO energy of the substituted analog. Such an ALMO-based frontier orbital analysis is expected to lead to intuitive principles for designing analogs of not only the oxyluciferin molecule, but also many other functional dyes.

Solvent Organization and Electrostatics Tuned by Solute Electronic Structure: Amide versus Non-Amide Carbonyls.

The ability to exploit carbonyl groups to measure electric fields in enzymes and other complex reactive environments by using the vibrational Stark effect has inspired growing interest in how these fields can be measured, tuned, and ultimately designed. Previous studies have concentrated on the role of the solvent in tuning the fields exerted on the solute. Here, we explore instead the role of the solute electronic structure in modifying the local solvent organization and electric field exerted on the solute. By measuring the infrared absorption spectra of amide-containing molecules, as prototypical peptides, and contrasting them with non-amide carbonyls in a wide range of solvents, we show that these solutes experience notable differences in their frequency shifts in polar solvents. Using vibrational Stark spectroscopy and molecular dynamics simulations, we demonstrate that while some of these differences can be rationalized by using the distinct intrinsic Stark tuning rates of the solutes, the larger frequency shifts for amides and dimethylurea primarily result from the larger solvent electric fields experienced by their carbonyl groups. These larger fields arise due to their stronger p-π conjugation, which results in larger C═O bond dipole moments that further induce substantial solvent organization. Using electronic structure calculations, we decompose the electric fields into contributions from solvent molecules that are in the first solvation shell and those from the bulk and show that both of these contributions are significant and become larger with enhanced conjugation in solutes. These results show that structural modifications of a solute can be used to tune both the solvent organization and electrostatic environment, indicating the importance of a solute-centric paradigm in modulating and designing the electrostatic environment in condensed-phase chemical processes.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges.

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

A two-directional vibrational probe reveals different electric field orientations in solution and an enzyme active site.

The catalytic power of an electric field depends on its magnitude and orientation with respect to the reactive chemical species. Understanding and designing new catalysts for electrostatic catalysis thus requires methods to measure the electric field orientation and magnitude at the molecular scale. We demonstrate that electric field orientations can be extracted using a two-directional vibrational probe by exploiting the vibrational Stark effect of both the C=O and C-D stretches of a deuterated aldehyde. Combining spectroscopy with molecular dynamics and electronic structure partitioning methods, we demonstrate that, despite distinct polarities, solvents act similarly in their preference for electrostatically stabilizing large bond dipoles at the expense of destabilizing small ones. In contrast, we find that for an active-site aldehyde inhibitor of liver alcohol dehydrogenase, the electric field orientation deviates markedly from that found in solvents, which provides direct evidence for the fundamental difference between the electrostatic environment of solvents and that of a preorganized enzyme active site.

Cavity quantum-electrodynamical time-dependent density functional theory within Gaussian atomic basis. II. Analytic energy gradient.

Following the formulation of cavity quantum-electrodynamical time-dependent density functional theory (cQED-TDDFT) models [Flick et al., ACS Photonics 6, 2757-2778 (2019) and Yang et al., J. Chem. Phys. 155, 064107 (2021)], here, we report the derivation and implementation of the analytic energy gradient for polaritonic states of a single photochrome within the cQED-TDDFT models. Such gradient evaluation is also applicable to a complex of explicitly specified photochromes or, with proper scaling, a set of parallel-oriented, identical-geometry, and non-interacting molecules in the microcavity.

Relativistic Orbital-Optimized Density Functional Theory for Accurate Core-Level Spectroscopy.

Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital-optimized density functional theory (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding a ∼0.5 eV root-mean-square error versus experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the >50 eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Correction: Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations.

Correction for 'Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations' by Junjie Yang et al., Phys. Chem. Chem. Phys., 2020, 22, 26838-26851, DOI: .

Elucidating the Electronic Structure of a Delayed Fluorescence Emitter via Orbital Interactions, Excitation Energy Components, Charge-Transfer Numbers, and Vibrational Reorganization Energies.

Recently, Wang and co-workers carried out frontier molecule orbital engineering in the design of m-Cz-BNCz, a thermally activated delayed fluorescence (TADF) molecule that emits pure green light at an external quantum efficiency of 27%. To further understand the underlying molecular design principles, we employed four advanced electronic structure analysis tools. First, an absolutely localized molecular orbitals (ALMO-) based analysis indicates an antibonding combination between the highest occupied molecular orbitals (HOMOs) of the donor 3,6-di-tert-butylcarbazole fragment and the acceptor BNCz fragment, which raises the HOMO energy and red-shifts the fluorescence emission wavelength. Second, excitation energy component analysis reveals that the S1-T1 gap is dominated by two-electron components of the excitation energies. Third, charge transfer number analysis, which is extended to use fragment-based Hirshfeld weights, indicates that the S1 and T1 excited states of m-Cz-BNCz (within time-dependent density functional theory) have notable charge transfer characters (27% for S1 and 12% for T1). This provides a balance between a small single-triplet gap and a substantial fluorescence intensity. Last, a vibrational reorganization energy analysis pinpoints the torsional motion between the BNCz and Cz moieties of m-Cz-BNCz as the source for its wider emission peak than that of p-Cz-BNCz. These four types of analyses are expected to be very valuable in the study and design of other TADF and functional dye molecules.

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis.

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N2, and BF bound to a [Ru(II)(NH3)5]2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.