 Award ID(s):
 2107923
 NSFPAR ID:
 10414820
 Author(s) / Creator(s):
 ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
 Date Published:
 Journal Name:
 Journal of chemical theory and computation
 ISSN:
 15499618
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spincrossover transitions. First, we discuss different representations of electronic states employed in the gridbased and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spinorbit couplings. Second, we describe three NAMD methods that have been used to simulate spincrossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration timedependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spindependent processes. These questions should be possible to address with future methodological developments in NAMD.more » « less

null (Ed.)Since the postulation of carbenes by Buchner (1903) and Staudinger (1912) as electrondeficient transient species carrying a divalent carbon atom, carbenes have emerged as key reactive intermediates in organic synthesis and in molecular mass growth processes leading eventually to carbonaceous nanostructures in the interstellar medium and in combustion systems. Contemplating the short lifetimes of these transient molecules and their tendency for dimerization, free carbenes represent one of the foremost obscured classes of organic reactive intermediates. Here, we afford an exceptional glance into the fundamentally unknown gasphase chemistry of preparing two prototype carbenes with distinct multiplicities—triplet pentadiynylidene (HCCCCCH) and singlet ethynylcyclopropenylidene (cC 5 H 2 ) carbene—via the elementary reaction of the simplest organic radical—methylidyne (CH)—with diacetylene (HCCCCH) under singlecollision conditions. Our combination of crossed molecular beam data with electronic structure calculations and quasiclassical trajectory simulations reveals fundamental reaction mechanisms and facilitates an intimate understanding of bondbreaking processes and isomerization processes of highly reactive hydrocarbon intermediates. The agreement between experimental chemical dynamics studies under singlecollision conditions and the outcome of trajectory simulations discloses that molecular beam studies merged with dynamics simulations have advanced to such a level that polyatomic reactions with relevance to extreme astrochemical and combustion chemistry conditions can be elucidated at the molecular level and expanded to higherorder homolog carbenes such as butadiynylcyclopropenylidene and triplet heptatriynylidene, thus offering a versatile strategy to explore the exotic chemistry of novel higherorder carbenes in the gas phase.more » « less

Among various techniques designed for studying openshell species, electron paramagnetic resonance (EPR) spectroscopy plays an important role. The key quantity measured by EPR is the gtensor, describing the coupling between an external magnetic field and molecular electronic spin. One theoretical framework for quantum chemistry calculations of gtensors is based on response theory, which involves substantial developments that are specific to the underlying electronic structure models. A simplified and easiertoimplement approach is based on the stateinteraction scheme, in which perturbation is included by considering a small number of states. We describe and benchmark the stateinteraction approach using equationofmotion coupledcluster and restrictedactivespace configuration interaction wave functions. The analysis confirms that this approach can deliver accurate results and highlights caveats of applying it, such as a choice of the reference state, convergence with respect to the number of states used in calculations, etc. The analysis also contributes toward a better understanding of challenges in calculations of higherorder properties using approximate wave functions.more » « less

Abstract Brownian dynamics (BD) is a computational method to simulate molecular diffusion processes. Although the BD method has been developed over several decades and is well established, new methodological developments are improving its accuracy, widening its scope, and increasing its application. In biological applications, BD is used to investigate the diffusive behavior of molecules subject to forces due to intermolecular interactions or interactions with material surfaces. BD can be used to compute rate constants for diffusional association, generate structures of encounter complexes for molecular binding partners, and examine the transport properties of geometrically complex molecules. Often, a series of simulations is performed, for example, for different protein mutants or environmental conditions, so that the effects of the changes on diffusional properties can be estimated. While biomolecules are commonly described at atomic resolution and internal molecular motions are typically neglected, coarse‐graining and the treatment of conformational flexibility are increasingly employed. Software packages for BD simulations of biomolecules are growing in capabilities, with several new packages providing novel features that expand the range of questions that can be addressed. These advances, when used in concert with experiment or other simulation methods, such as molecular dynamics, open new opportunities for application to biochemical and biological systems. Here, we review some of the latest developments in the theory, methods, software, and applications of BD simulations to study biomolecular diffusional association processes and provide a perspective on their future use and application to outstanding challenges in biology, bioengineering, and biomedicine.
This article is categorized under:
Structure and Mechanism > Computational Biochemistry and Biophysics
Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods
Software > Simulation Methods

ABSTRACT In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT‐PBC). Calculations are compared with experimental measurements and show semi‐quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X‐ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.