Cavity-modified molecular dipole switching dynamics.

Polaritonic states, which are formed by resonances between a molecular excitation and the photonic mode of a cavity, have a number of useful properties that offer new routes to control molecular photochemistry using electric fields. To provide a theoretical description of how polaritonic states affect the real-time electron dynamics in molecules, a new method is described where the effects of strong light-molecule coupling are implemented using real-time electronic structure theory. The coupling between the molecular electronic states and the cavity is described by the Pauli-Fierz Hamiltonian, and transitions between polaritonic states are induced via an external time-dependent electric field using time-dependent configuration interaction (TDCI) theory, producing quantum electrodynamics TDCI (QED-TDCI). This method is used to study laser-induced ultrafast charge transfer and dipole-switching dynamics of the LiCN molecule inside a cavity. The increase in cavity coupling strength is found to have a significant impact on the energies and transition dipole moments of the molecule-cavity system. The convergence of the polaritonic state energies as a function of the number of included electronic and photonic basis states is discussed.

Multi-configuration electron-nuclear dynamics: An open-shell approach.

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.

Many-photon excitation of organic molecules in a cavity-Superradiance as a measure of coherence.

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.

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium-Ligand Bonds.

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.

Dynamics and spectroscopy of molecular ensembles in a lossy microcavity.

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.

The first microsolvation step for furans: New experiments and benchmarking strategies.

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.

Modifying the Nonradiative Decay Dynamics through Conical Intersections via Collective Coupling to a Cavity Mode.

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.

Spin trapping and flipping in FeCO through relativistic electron dynamics.

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 furan microsolvation blind challenge for quantum chemical methods: First steps.

Herein we present the results of a blind challenge to quantum chemical methods in the calculation of dimerization preferences in the low temperature gas phase. The target of study was the first step of the microsolvation of furan, 2-methylfuran and 2,5-dimethylfuran with methanol. The dimers were investigated through IR spectroscopy of a supersonic jet expansion. From the measured bands, it was possible to identify a persistent hydrogen bonding OH-O motif in the predominant species. From the presence of another band, which can be attributed to an OH-π interaction, we were able to assert that the energy gap between the two types of dimers should be less than or close to 1 kJ/mol across the series. These values served as a first evaluation ruler for the 12 entries featured in the challenge. A tentative stricter evaluation of the challenge results is also carried out, combining theoretical and experimental results in order to define a smaller error bar. The process was carried out in a double-blind fashion, with both theory and experimental groups unaware of the results on the other side, with the exception of the 2,5-dimethylfuran system which was featured in an earlier publication.

The role of the CI expansion length in time-dependent studies.

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.

Quantum and quasi-classical collisional dynamics of O2-Ar at high temperatures.

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.

Effects of roaming trajectories on the transition state theory rates of a reduced-dimensional model of ketene isomerization.

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.

Remarks on the validity of the fixed nuclei approximation in quantum electron dynamics.

In this paper we report quantum dynamical simulations that test the fixed nuclei approximation, which is usually invoked in ab initio correlated electron dynamics. We do so by employing a recently developed method, called multiconfiguration electron nuclear dynamics. Additionally, we discuss the influence of the multiconfiguration expansion length on the results.

The multi-configuration electron-nuclear dynamics method applied to LiH.

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.

Correlated electron dynamics: how aromaticity can be controlled.

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.