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. more...

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. more...

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. more...

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. more...

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). more...

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. more...

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. more...

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. more...

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. more...