Search

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy De and bond length re are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15-0.19. For the single-reference methods, Kohn-Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the De and re of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core-valence terms; this results in highly accurate dissociation energies and good bond lengths. Averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

The multi-configuration electron–nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.

The Born–Oppenheimer approximation is the fundamental approximation in the quantum-mechanical description of molecules, and holds true in most applications for ground-state properties and to a less...

The equilibrium geometries and spectroscopic properties of two key singlet carbenes, buta-1,3-diynylcarbene (6) and 2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (9), which have not been experimentally observed to date, are investigated using high-level coupled-cluster (CC) methods. The current theoretical study necessitates new experimental data on C5H2 isomers considering the relevance of these molecules to interstellar chemistry. Bent-pentadiynylidene (4) has been missing in the laboratory and the prime focus of our earlier theoretical work. The present theoretical study indicates that isomers 6 and 9 are also viable experimental targets. Apart from ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 3-(didehydrovinylidene)cyclopropene (8), which are identified by Fourier transform microwave spectroscopy, the dipole moments of elusive 4, 6, and 9 are also nonzero (μ ≠ 0). The relative energies of these isomers, calculated at the CCSDT(Q)/CBS level of theory, with respect to linear triplet pentadiynylidene (1) reveal that they all lie within 25.1 kcal mol-1. Therefore, geometric, energetic, aromatic, and spectroscopic parameters are reported here, which may assist the efforts of molecular spectroscopists in the future. Anharmonic vibrational calculations on isomers 6 and 9 indicate that the former is loosely bound and would be challenging to be detected experimentally. Among the undetected carbenes, 9 may be considered as a potential target molecule considering its higher polarity and aromatic nature.

Coherent excitation of a molecular ensemble coupled to a common radiation mode can lead to the collective emission of radiation known as superradiance. This collective emission only occurs if there is an entanglement between the molecules in their ground and excited state and can, therefore, serve as a macroscopic measure of coherence in the ensemble. Reported here are wave packet propagations for various pyrazine models of increasing complexity and molecular ensembles thereof. We show that ensemble coherence upon photoexcitation can prevail up to relatively long time scales although the effect can diminish quickly with increasing ensemble size. Coherence can also build up over time and even reemerge after the molecules have passed through a conical intersection. The effect of the pump pulse characteristics on the collective response of the molecular ensemble is also studied. A broadband pulse imprints a large amount of initial coherence to the system, as compared to a longer pulse with a smaller spread in the frequency domain. However, the differential effects arising from a different pulse duration and coherent bandwidth become less prominent if the emission of light from the ensemble takes place after a non-adiabatic decay process.

Coherent excitation of a molecular ensemble coupled to a common radiation mode can lead to the collective emission of radiation known as superradiance. This collective emission only occurs if there is an entanglement between the molecules in their ground and excited state and can therefore serve as a macroscopic measure of coherence in the ensemble. Reported here are wave packet propagations for various pyrazine models of increasing complexity and molecular ensembles thereof. We show that ensemble coherence upon photoexcitation can prevail up to relatively long time scales, although the effect can diminish quickly with increasing ensemble size. Coherence can also build up over time and even reemerge after the molecules have passed through a conical intersection. The effect of the pump-pulse characteristics on the collective response of the molecular ensemble is also studied. A broad-band pulse imprints a large amount of initial coherence to the system, as compared to a longer pulse with a smaller spread in the frequency domain. However, the differential effects arising from a different pulse duration and coherent bandwidth become less prominent if the emission of light from the ensemble takes place after a non-adiabatic decay process.

The radiative and nonradiative relaxation dynamics of an ensemble of molecules in a microcavity are investigated with emphasis on the impact of the cavity lifetime on reactive and spectroscopic properties. Extending a previous study [I. S. Ulusoy et al., J. Phys. Chem. A 123, 8832-8844 (2019)], it is shown that the dynamics of the ensemble and of single molecules are influenced by the presence of a cavity resonance as long as the polariton splitting can be resolved spectroscopically, which critically depends on the lifetime of the system. Our simulations illustrate how the branching between nonradiative intersystem crossing and radiative decay through the cavity can be tuned by selecting specific cavity photon energies resonant at specific molecular geometries. In the case of cavity-photon energies that are not resonant at the Franck-Condon geometry of the molecules, it is demonstrated numerically and analytically that collective effects are limited to a handful of molecules in the ensemble.

Transition metal compounds are very versatile, and their characteristics can differ profoundly depending on their electronic structure. Compounds in which a spin transition from a low-spin to a high-spin state can be achieved through means of an optical excitation are particularly intriguing, as a controlled spin-flip opens promising avenues in areas such as sensing, information technology, molecular switches and energy technology. The fundamental mechanisms in spin crossover and spin transitions remain unanswered, due to the complexity of electronic structure and interplay of relativistic effects. Presented here is a new approach that allows the first direct study of spin flip dynamics through a mapping of spin-mixed to spin-pure states. The method is applied to FeCO and addresses the spin-flip dynamics during a spin transition. Wave packets that combine different spin states are generated through optical excitation and relevant mechanisms in optically triggered spin transitions are discussed.

The journal of chemical physics 152(16), 164303 (2020). doi:10.1063/5.0004465, The site-specific first microsolvation step of furan and some of its derivatives with methanol is explored to benchmark the ability of quantum-chemical methods to describe the structure, energetics, and vibrational spectrum at low temperature. Infrared and microwave spectra in supersonic jet expansions are used to quantify the docking preference and some relevant quantum states of the model complexes. Microwave spectroscopy strictly rules out in-plane docking of methanol as opposed to the top coordination of the aromatic ring. Contrasting comparison strategies, which emphasize either the experimental or the theoretical input, are explored. Within the harmonic approximation, only a few composite computational approaches are able to achieve a satisfactory performance. Deuteration experiments suggest that the harmonic treatment itself is largely justified for the zero-point energy, likely and by design due to the systematic cancellation of important anharmonic contributions between the docking variants. Therefore, discrepancies between experiment and theory for the isomer abundance are tentatively assigned to electronic structure deficiencies, but uncertainties remain on the nuclear dynamics side. Attempts to include anharmonic contributions indicate that for systems of this size, a uniform treatment of anharmonicity with systematically improved performance is not yet in sight., Published by American Institute of Physics, Melville, NY

The coupling of a molecular ensemble to the confined electromagnetic modes of a microcavity can strongly modify the photophysics and photochemistry of the molecules upon photoexcitation. We investigate here how collective coupling effects lead to modifications of the mechanisms and rates of photochemical processes, in particular, photodissociation and nonradiative decay in NaI and pyrazine, respectively. We show that, after direct excitation into the lower polaritonic states, the lower-energy light-matter hybrid states, the dynamics of the molecular ensemble coupled to light is very similar to the dynamics of the corresponding isolated molecules. Conversely, excitation into the upper polaritonic states results in more complex dynamics that involve as a first step the population transfer toward the manifold of intermediate dark states. These dynamics differ substantially from those of the isolated molecules and may result in measurable time delays for nonradiative decay or excited-state reaction mechanisms. Similarly, we describe how addition of a buffer of nonreactive molecules coupled to the cavity mode can be used to delay the onset of the photochemical processes of the reactive part of the ensemble, where the buffer medium is more effective in inhibiting the reactive process than only reactive molecules in the cavity.

Electron motion in molecules can be initiated by a laser pulse and, through the induced charge migration, can, in turn, trigger a response from the nuclei. The generated electron-nuclear wave packets can quickly decohere, but also lead to recurrence of coherences. Presented here is a quantum dynamics study of ${\mathrm{H}}_{2}$ and LiH in strong fields. No separation of electronic and nuclear motion is assumed, and no precomputed potential energy surface is required. It is shown that electron and nuclear motion are strongly coupled through the laser-pulse excitation, and already moderate nuclear excitation can substantially impact the electron dynamics. Through the full quantum dynamical treatment, coherence properties of electronic and nuclear wave function can be investigated. The nuclear motion leads to decoherences while it can also introduce recurrence of coherences.

Basis functions and molecular integrals are critical components utilized in electronic structure (ab initio and density functional theory) calculations. In this chapter, a mathematical overview is provided about the construction of basis functions and computation of molecular integrals. Included is a discussion about the selection of basis functions for use in computational chemistry calculations, providing insight about factors that should be considered. The behavior of energies and properties with increasing basis set size, as well as extrapolation schemes are presented. Also, a discussion of the computational scaling and more efficient approaches, such as pseudopotential basis sets and the density fitting, are described.

With the recent advances in experimental attosecond science, theoretical predictions of electron dynamics can now be validated against experiment. Time-dependent studies of the electron motion in molecules can be used to obtain information about electronic transitions and the interaction of the electrons with electromagnetic fields. Often, these approaches rely on single-excited wave functions. Presented here is a first attempt to evaluate the accuracy of the time-dependent configuration interaction method so that the optimal representation of the electronic wave function for time-dependent studies can be assessed. A quantifier is determined that can aid in finding this optimal representation. The approach is demonstrated on a variety of molecules that include both localized and intramolecular charge transfer electron excitations. Observables including excitation energies, dipole moments, strengths, and static polarizabilities are obtained from time-independent and time-dependent calculations and are compared to experimental data. In this way, a rigorous routine is developed by which the reliability and accuracy of the CI wave function can be assessed and which represents a first step to a more quantitative description of electron dynamics in molecules.

The journal of chemical physics 148(1), 1-13 (2018). doi:10.1063/1.5009011, Published by American Institute of Physics, Melville, NY

The rates of chemical reactions (or any activated process) are by definition determined by the flux of reactants (or initial states) that end up as products (or final states). The forward flux through any surface that divides reactants from products is a constant as long as only those trajectories that were reactants in the infinite past and products in the infinite future are included in the flux once and only once. Transition state theory (TST) ignores this last clause, thereby overestimating the rate if any of the trajectories recross the dividing surface. However, its advantage is that it replaces a dynamical calculation with a statistical integral over the TST geometry. The recent identification of roaming trajectories-those that persist for a long time as neither reactant nor product without ever visiting near the col on the energy landscape-apparently challenges the dogma that TST's only error lies in the omission of recrossing trajectories. This question is investigated using the isomerization reaction of ketene in which the experimental values are verified to be in reasonable agreement with both the exact and TST values. We have found two trajectories for the ketene isomerization that carry the signature of roaming, but their effect on the calculation of the reaction rate constant using classical transition state theory is small. Indeed, the existence of roaming trajectories is seen to impose a limitation on which dividing surfaces are appropriate for the calculation of either exact or approximate TST rates, but in this case, they do not unseat the existence of dividing surfaces that can be used safely to calculate TST rates.

We present a theoretical framework for the computation of anharmonic vibrational frequencies for large systems, with a particular focus on determining adsorbate frequencies from first principles. We give a detailed account of our local implementation of the vibrational self-consistent field approach and its correlation corrections. We show that our approach is both robust, accurate and can be easily deployed on computational grids in order to provide an efficient computational tool. We also present results on the vibrational spectrum of hydrogen fluoride on pyrene, on the thiophene molecule in the gas phase, and on small neutral gold clusters.

Roaming dynamics have been observed in a three-dimensional model of the ketene isomerization reaction. The roaming trajectories sample the region between the outer potential barriers closest to the respective ketene isomers and involve turning points along the reaction coordinate in a polar representation. These roaming trajectories avoid the intrinsic reaction coordinate and the intermediates to which it is associated. Thus, one-dimensional transition state theory (TST) is generally insufficient as has been confirmed through an analysis of the reactive flux along the dividing surface (DS). A global representation of the DS, however, leads to accurate TST rate constants. The exact and TST microcanonical rates of isomerization have been obtained for the three-dimensional model and compare well to experiment. The global DS is therefore particularly important for obtaining rates in reactions that exhibit roaming. This work thus confirms the findings of our previous two-dimensional treatment of ketene isomerization (Ulusoy et al. in J. Phys. Chem. A 117:7553–7560, 2013).

In this paper, we present quantum dynamical calculations on electron correlation dynamics in atoms and molecules using explicitly time-dependent ab initio configuration interaction theory. The goals are (i) to show that in which cases it is possible to switch off the electronic correlation by ultrashort laser pulses, and (ii) to understand the temporal evolution and the time scale on which it reappears. We characterize the appearance of correlation through electron-electron scattering when starting from an uncorrelated state, and we identify pathways for the preparation of a Hartree-Fock state from the correlated, true ground state. Exemplary results for noble gases, alkaline earth elements, and selected molecules are provided. For Mg we show that the uncorrelated state can be prepared using a shaped ultrashort laser pulse.

A hypersonic vehicle traveling at a high speed disrupts the distribution of internal states in the ambient flow and introduces a nonequilibrium distribution in the post-shock conditions. We investigate the vibrational relaxation in diatom-atom collisions in the range of temperatures between 1000 and 10 000 K by comparing results of extensive fully quantum-mechanical and quasi-classical simulations with available experimental data. The present paper simulates the interaction of molecular oxygen with argon as the first step in developing the aerothermodynamics models based on first principles. We devise a routine to standardize such calculations also for other scattering systems. Our results demonstrate very good agreement of vibrational relaxation time, derived from quantum-mechanical calculations with the experimental measurements conducted in shock tube facilities. At the same time, the quasi-classical simulations fail to accurately predict rates of vibrationally inelastic transitions at temperatures lower than 3000 K. This observation and the computational cost of adopted methods suggest that the next generation of high fidelity thermochemical models should be a combination of quantum and quasi-classical approaches.

In this Article, we show that the aromaticity of a molecule can be turned off by controlling the electron dynamics. We present a controlled switching from the aromatic ground state of benzene to two different nonaromatic states, using a laser pulse. The propagation of the molecular wave function is carried out with the time-dependent configuration interaction method. The laser pulse for switching between the ground and excited states is optimized using optimal control theory. Bond orders and Mulliken charges serve as an aromaticity criterion. The nonaromatic target states exhibit localized bonds and partial charges on the carbon atoms; these localized electrons circulate on an attosecond time scale in the ring system.

Vibrational spectroscopy is a powerful tool to identify molecules and to characterise their chemical state. Inelastic electron tunnelling spectroscopy (IETS) combined with scanning tunnelling microscopy (STM) allows the application of vibrational analysis to a single molecule. Up to now, IETS was restricted to small species due to the complexity of vibration spectra for larger molecules. We extend the horizon of IETS for both experiment and theory by measuring the STM-IETS spectra of mercaptopyridine adsorbed on the (111) surface of gold and comparing it to theoretical spectra. Such complex spectra with more than 20 lines can be reliably determined and computed leading to completely new insights. Experimentally, the vibrational spectra exhibit a dependence on the specific adsorption site of the molecules. Theoretically, this dependence is only accessible if anharmonic contributions to the interaction potentials are included. These joint experimental and theoretical advances open new perspectives for structure determination of organic adlayers.

The multi-configuration electron-nuclear dynamics (MCEND) method is a nonadiabatic quantum dynamics approach to the description of molecular processes. MCEND is a combination of the multi-configuration time-dependent Hartree (MCTDH) method for atoms and its antisymmetrized equivalent MCTDHF for electrons. The purpose of this method is to simultaneously describe nuclear and electronic wave packets in a quantum dynamical way, without the need to calculate potential energy surfaces and diabatic coupling functions. In this paper we present first exemplary calculations of MCEND applied to the LiH molecule, and discuss computational and numerical details of our implementation.