ABSTRACT
Fitting coupled potential energy surfaces is a critical step in simulating electronically nonadiabatic chemical reactions and energy transfer processes. Analytic representation of coupled potential energy surfaces enables one to perform detailed dynamics calculations. Traditionally, fitting is performed in a diabatic representation to avoid fitting the cuspidal ridges of coupled adiabatic potential energy surfaces at conical intersection seams. In this work, we provide an alternative approach by carrying out fitting in the adiabatic representation using a modified version of the Frobenius companion matrices, whose usage was first proposed by Opalka and Domcke. Their work involved minimizing the errors in fits of the characteristic polynomial coefficients (CPCs) and diagonalizing the resulting companion matrix, whose eigenvalues are adiabatic potential energies. We show, however, that this may lead to complex eigenvalues and spurious discontinuities. To alleviate this problem, we provide a new procedure for the automatic discovery of CPCs and the diagonalization of a companion matrix by using a special neural network architecture. The method effectively allows analytic representation of global coupled adiabatic potential energy surfaces and their gradients with only adiabatic energy input and without experience-based selection of a diabatization scheme. We demonstrate that the new procedure, called the companion matrix neural network (CMNN), is successful by showing applications to LiH, H3, phenol, and thiophenol.
ABSTRACT
Broad and diverse sets of accurate data provide useful metrics for assessing the performance of new theoretical methods. However, assessing methods against large databases can be an arduous task. Here, we present 17 representative energetic databases, defined as small databases whose errors and error spreads are representative of larger databases and which therefore can serve as efficient benchmarks for developing and testing electronic structure methods and density functionals. In 15 cases, the representative databases have 6 entries while being representative of larger databases with 14-107 entries, and in the other two cases, they have 14 entries while being representative of larger databases with 418-455 entries. The mean unsigned error (MUE) of 100 electronic structure methods on a given representative database is typically within about 8% of the MUE on its parent database, and the root-mean-square error (RMSE) is typically within about 11% of the RMSE on the parent database. Thus, the representative databases are quite successful in indicating accuracy while maintaining good diversity. The databases include both main-group and transition-metal compounds and reactions, and they include bond energies, reaction energies, barrier heights, noncovalent interactions, ionization potentials, and absolute energies.
ABSTRACT
We combine in situ laser spectroscopy, quantum chemistry, and kinetic calculations to study the reaction of a singlet oxygen atom with dimethyl ether. Infrared laser absorption spectroscopy and Faraday rotation spectroscopy are used for the detection and quantification of the reaction products OH, H2O, HO2, and CH2O on submillisecond time scales. Fitting temporal profiles of products with simulations using an in-house reaction mechanism allows product branching to be quantified at 30, 60, and 150 Torr. The experimentally determined product branching agrees well with master equation calculations based on electronic structure data and transition state theory. The calculations demonstrate that the dimethyl peroxide (CH3OOCH3) generated via O-insertion into the C-O bond undergoes subsequent dissociation to CH3O + CH3O through energetically favored reactions without an intrinsic barrier. This O-insertion mechanism can be important for understanding the fate of biofuels leaking into the atmosphere and for plasma-based biofuel processing technologies.
ABSTRACT
Based on quantum mechanically guided experiments that observed elusive intermediates in the domain of inception that lies between large molecules and soot particles, we provide a new mechanism for the formation of carbonaceous particles from gas-phase molecular precursors. We investigated the clustering behavior of resonantly stabilized radicals (RSRs) and their interactions with unsaturated hydrocarbons through a combination of gas-phase reaction experiments and theoretical calculations. Our research directly observed a sequence of covalently bound clusters (CBCs) as key intermediates in the evolution from small RSRs, such as benzyl (C7H7), indenyl (C9H7), 1-methylnaphthyl (1-C11H9), and 2-methylnaphthyl (2-C11H9), to large polycyclic aromatic hydrocarbons (PAHs) consisting of 28 to 55 carbons. We found that hydrogen abstraction and RSR addition drive the formation and growth of CBCs, leading to progressive H-losses, the generation of large PAHs and PAH radicals, and the formation of white smoke (incipient carbonaceous particles). This mechanism of progressive H-losses from CBCs (PHLCBC) elucidates the crucial relationship among RSRs, CBCs, and PAHs, and this study provides an unprecedentedly seamless path of observed assembly from small RSRs to large nanoparticles. Understanding the PHLCBC mechanism over a wide temperature range may enhance the accuracy of multiscale models of soot formation, guide the synthesis of carbonaceous nanomaterials, and deepen our understanding of the origin and evolution of carbon within our galaxy.
ABSTRACT
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
ABSTRACT
Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.
ABSTRACT
We present an improvement of the local-pair zero-point-energy (LP-ZPE) scheme of Mukherjee and Barbatti. The new approximation is called the improved LP-ZPE scheme or iLP-ZPE. This scheme can produce trajectories that do not have unphysical leaking of zero-point energy from high-frequency spectator modes into low-frequency modes. We illustrate the method with a successful direct dynamics application to the Ne···HF van der Waals molecule. The method is well suited for direct dynamics calculations because it does not require costly evaluations of local Hessians or instantaneous normal modes along the trajectories.
ABSTRACT
Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.
ABSTRACT
We report new potential energy surfaces for six coupled 5A' states and 14 coupled 3A' states of O3. The new surfaces are created by parametrically managed diabatization by deep neural network (PM-DDNN). The PM-DDNN method uses calculated adiabatic potential energy surfaces to discover and fit an underlying adiabatic-equivalent set of diabatic surfaces and their couplings and obtains the fit to the adiabatic surfaces by diagonalization of the diabatic potential energy matrix (DPEM). The procedure yields the adiabatic surfaces and their gradients, as well as the DPEM and its gradient. If desired one can also compute the nonadiabatic coupling due to the transformation. The present work improves on previous work by using a new coordinate to guide the decay of the neural network contribution to the many-body fit to the whole DPEM. The main objective was to obtain smoother potentials than the previous ones with better suitability for dynamics calculations, and this was achieved. Furthermore, we obtained suitably small deviations from the input reference data. For the six coupled 5A' surfaces, the 60,366 data below 10 eV are fit with a mean unsigned error (MUE) of 49 meV, and for the 14 coupled 3A' surfaces, the 76,733 data below 10 eV are fit with an MUE of 28 meV. The data below 5 eV fit even more accurately with MUEs of 37 meV (5A') and 20 meV (3A').
ABSTRACT
Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree-Fock exchange. These results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.
ABSTRACT
Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. In this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S1/S0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.
ABSTRACT
The accuracy of Kohn-Sham density functional theory depends strongly on the approximation to the exchange-correlation functional. In this work, we present a new exchange-correlation functional called M11pz (M11 plus rung-3.5 terms with zero Hartree-Fock exchange) that is built on the M11plus functional with the goal of using its rung-3.5 terms without a Hartree-Fock exchange term, especially to improve the accuracy for strongly correlated systems. The M11pz functional is optimized with the same local and rung-3.5 ingredients that are used in M11plus but without any percentage of Hartree-Fock exchange. The performance of M11pz is compared with eight local functionals, and M11pz is found to be in top three when the errors or ranks are averaged over eight grouped and partially overlapping databases: AME418/22, atomic and molecular energies; MGBE172, main-group bond energies; TMBE40, transition-metal bond energies; SR309, single-reference systems; MR54, multireference systems; BH192, barrier heights; NC579, noncovalent interaction energies; and MS20, molecular structures. For calculations of band gaps of solids, M11pz is the second best of the nine tested functionals that have zero Hartree-Fock exchange.
ABSTRACT
Constructing analytic representations of global and semiglobal potential energy surfaces is difficult and can be laborious, and it is even harder when one needs coupled potential energy surfaces and their electronically nonadiabatic couplings. When accomplished, however, the resulting potential functions are a valuable resource. To facilitate the convenient use of potentials that have been developed, we provide a collection of existing surfaces in a library with consistent units and formats. A potential energy surface library of this type, namely PotLib, was built more than 20 years ago. However, that library only provided pristine Fortran subroutines for each potential energy surface, and therefore, it is not as user-friendly as would be desirable. Here, we report the creation of ChemPotPy, a CHEMical library of POTential energy surfaces in PYthon. ChemPotPy is a user-friendly library for analytic representation of single-state and multistate potential energy surfaces and couplings. A given entry in the library contains an analytic potential energy function or analytic functions for a set of coupled potential energy surfaces, and depending on the case, it may also include analytic or numerical gradients, nonadiabatic coupling vectors, and/or diabatic potential energy matrices and their gradients. Only three inputs, namely, the chemical formula of the system, the name of the potential energy surface or surface set, and the Cartesian geometry, are required. ChemPotPy uses the same units for input and output quantities of all surfaces and surface sets to facilitate general interfaces with the dynamics programs. The initial version of the library contains 338 entries, and we anticipate that more will be added in the future.
ABSTRACT
The predictive ability of density functional theory is fundamental to its usefulness in chemical applications. Recent work has compared solution-phase enthalpies of activation for metal-ligand bond dissociation to enthalpies of reaction for bond dissociation, and the present work continues those comparisons for 43 density functional methods. The results for ligand dissociation enthalpies of 30 metal-ligand complexes tested in this work reveal significant inadequacies of some functionals as well as challenges from the dispersion corrections to some functionals. The analysis presented here demonstrates the excellent performance of a recent density functional, M11plus, which contains nonlocal rung-3.5 correlation. We also find a good agreement between theory and experiment for some functionals without empirical dispersion corrections such as M06, r2SCAN, M06-L, and revM11, as well as good performance for some functionals with added dispersion corrections such as ωB97X-D (which always has a correction) and BLYP, B3LYP, CAM-B3LYP, and PBE0 when the optional dispersion corrections are added.
ABSTRACT
The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.
ABSTRACT
Activating the C-H bonds of alkanes without further oxidation to more thermodynamically stable products, CO and CO2, is a long-sought goal of catalytic chemistry. Inspired by the monocopper active site of methane monooxygenase, we synthesized a Cu-doped ZIF-8 metal-organic framework with 25% Cu and 75% Zn in the nodes and activated it by heating to 200 °C and dosing in a stepwise fashion with O2, methane, and steam. We found that it does oxidize methane to methanol and formaldehyde. The catalysis persists through at least five cycles, and beyond the third cycle, the selectivity improves to the extent that no CO2 can be detected. Experimental characterization and analysis were carried out by PXRD, DRUV-vis, SEM, and XAS (XANES and EXAFS). The reaction is postulated to proceed at open-coordination copper sites generated by defects, and the mechanism of methanol production was explicated by density functional calculations with the revMO6-L exchange-correlation functional. The calculations reveal a catalytic cycle of oxygen-activated CuI involving the conversion of two molecules of CH4 to two molecules of CH3OH by a sequence of hydrogen atom transfer reactions and rebound steps. For most steps in the cycle, the reaction is more favored by singlet species than by triplets.
ABSTRACT
Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.
ABSTRACT
We present hybrid multiconfiguration density coherence functional theory (HMC-DCFT), and we optimize a density coherence functional by parametrization against a diverse data set of 59 bond energies and 60 barrier heights. We compare the results to calculations on the same data set by CASSCF, CASPT2, six Kohn-Sham and hybrid Kohn-Sham exchange-correlation functionals, and three on-top functionals for pair-density functional theory (PDFT) and hybrid PDFT. The new functional has better accuracy than all compared methods.
ABSTRACT
Although experimental methods can be used to obtain the quantitative kinetics of atmospheric reactions, experimental data are often limited to a narrow temperature range. The reaction of SO3 with water vapor is important for elucidating the formation of sulfuric acid in the atmosphere; however, the kinetics is uncertain at low temperatures. Here, we calculate rate constants for reactions of sulfur trioxide with two water molecules. We consider two mechanisms: the SO3···H2O + H2O reaction and the SO3 + (H2O)2 reaction. We find that beyond-CCSD(T) contributions to the barrier heights are very large, and multidimensional tunneling, unusually large anharmonicity of high-frequency modes, and torsional anharmonicity are important for obtaining quantitative kinetics. We find that at lower temperatures, the formation of the termolecular precursor complexes, which is often neglected, is rate-limiting compared to passage through the tight transition states. Our calculations show that the SO3···H2O + H2O mechanism is more important than the SO3 + (H2O)2 mechanism at 5-50 km altitudes. We find that the rate ratio between SO3···H2O + H2O and SO3 + (H2O)2 is greater than 20 at altitudes between 10 and 35 km, where the concentration of SO3 is very high.
ABSTRACT
The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH2OH)CHOO], for its reactions with H2O and (H2O)2, and for the reaction of the E-(CH2OH)CHOO H2O complex with H2O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH2OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH2OH)CHOO H2O + H2O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H2O)2. In addition, we find that the atmospheric lifetime of E-(CH2OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10-6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH2OH)CHOO enhances its reactivity.