103 Search Results

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn–Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids canmore » be approximated with individual CSFs, despite the presence of multireference character in the ground state. Here, we also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O₃. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.« less

Controlling excited-state reactivity is a long-standing challenge in photochemistry, as a desired pathway may be inaccessible or compete with other unwanted channels. An important example is internal conversion of the anionic green fluorescent protein (GFP) chromophore where non-selective progress along two competing torsional modes (P: phenolate and I: imidazolinone) impairs and enables Z-to-E photoisomerization, respectively. Developing strategies to promote photoisomerization could drive new areas of applications of GFP-like proteins. Motivated by the charge-transfer dichotomy of the torsional modes, we explore chemical substitution on the P-ring of the chromophore as a way to control excited-state pathways and improve photoisomerization. As demonstratedmore » by methoxylation, selective P-twisting appears difficult to achieve because the electron-donating potential effects of the substituents are counteracted by inertial effects that directly retard the motion. Conversely, these effects act in concert to promote I-twisting when introducing electron-withdrawing groups. Specifically, 2,3,5-trifluorination leads to both pathway selectivity and a more direct approach to the I-twisted intersection which, in turn, doubles the photoisomerization quantum yield. Our results suggest P-ring engineering as an effective approach to boost photoisomerization of the anionic GFP chromophore.« less

Scalable numerical solutions to the time dependent Schrodinger equation remain an outstanding goal in theoretical chemistry. Here we present a method which utilises recent breakthroughs in signal processing to consistently adapt a dictionary of basis functions to the dynamics of the system. Further, we show that for two low-dimensional model problems the size of the basis set does not grow quickly with time and appears only weakly dependent on dimensionality. The generality of this finding remains to be seen. The method primarily uses energies and gradients of the potential, opening the possibility for its use in on-the-fly ab initio quantummore » wavepacket dynamics.« less

It has been well-established that the topography around conical intersections between excited electronic states is incorrectly described by coupled cluster and many other single reference theories (the intersections are “defective”). Despite this, we show both analytically and numerically that the geometric phase effect (GPE) is correctly reproduced upon traversing a path around a defective excited-state conical intersection (CI) in coupled cluster theory. The theoretical analysis is carried out by using a non-Hermitian generalization of the linear vibronic coupling approach. Interestingly, the approach qualitatively explains the characteristic (incorrect) shape of the defective CIs and CI seams. Moreover, the validity of themore » approach and the presence of the GPE indicate that defective CIs are local (and not global) artifacts. This implies that a sufficiently accurate coupled cluster method could predict nuclear dynamics, including geometric phase effects, as long as the nuclear wavepacket never gets too close to the conical intersections.« less

Electrocyclic reactions are characterized by the concerted formation and cleavage of both σ and π bonds through a cyclic structure. This structure is known as a pericyclic transition state for thermal reactions and a pericyclic minimum in the excited state for photochemical reactions. However, the structure of the pericyclic geometry has yet to be observed experimentally. We use a combination of ultrafast electron diffraction and excited state wavepacket simulations to image structural dynamics through the pericyclic minimum of a photochemical electrocyclic ring-opening reaction in the molecule α-terpinene. The structural motion into the pericyclic minimum is dominated by rehybridization of twomore » carbon atoms, which is required for the transformation from two to three conjugated π bonds. The σ bond dissociation largely happens after internal conversion from the pericyclic minimum to the electronic ground state. These findings may be transferrable to electrocyclic reactions in general.« less

Machine learning (ML) offers an attractive method for making predictions about molecular systems while circumventing the need to run expensive electronic structure calculations. Once trained on ab initio data, the promise of ML is to deliver accurate predictions of molecular properties that were previously computationally infeasible. In this work, we develop and train a graph neural network model to correct the basis set incompleteness error (BSIE) between a small and large basis set at the RHF and B3LYP levels of theory. Our results show that, when compared to fitting to the total potential, an ML model fitted to correct themore » BSIE is better at generalizing to systems not seen during training. We test this ability by training on single molecules while evaluating on molecular complexes. We also show that ensemble models yield better behaved potentials in situations where the training data is insufficient. However, even when only fitting to the BSIE, acceptable performance is only achieved when the training data sufficiently resemble the systems one wants to make predictions on. As a result, the test error of the final model trained to predict the difference between the cc-pVDZ and cc-pV5Z potential is 0.184 kcal/mol for the B3LYP density functional, and the ensemble model accurately reproduces the large basis set interaction energy curves on the S66x8 dataset.« less

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theorymore » based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.« less

Light-driven molecular motors have generated considerable interest due to their potential applications in material and biological systems. Recently, Greb and Lehn reported a new class of molecular motors, chiral N-alkyl imines, which undergo unidirectional rotation induced by light and heat. The mechanism of unidirectional motion in molecular motors containing a C═N group has been assumed to consist of photoinduced torsion about the double bond. In this work, we present a computational study of the photoisomerization dynamics of a chiral N-alkyl imine motor. We find that the location and energetics of minimal energy conical intersections (MECIs) alone are insufficient to understandmore » the mechanism of the motor. Furthermore, a key part of the mechanism consists of out-of-plane distortions of the N atom (followed by isomerization about the double bond). Dynamic effects and out-of-plane distortions are critical to understand the observed (rather low) quantum yield for photoisomerization. Finally, our results provide hints as to how the photoisomerization quantum yield might be increased, improving the efficiency of this class of molecular motors.« less

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to themore » GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.« less

The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. Here, we describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google’s protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs withmore » electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ~250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI^-)—a model chromophore in green fluorescent proteins—on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.« less