Measurement of Prandtl number and thermal conductivity Summary report

Prandtl numbers and thermal conductivity of air, argon, and hydrocarbon fuel combustion product

On the Theoretical Determination of Photolysis Properties for Atmospheric Volatile Organic Compounds

Volatile organic compounds (VOC) are ubiquitous atmospheric molecules that generate a complex network of chemical reactions in the troposphere, often triggered by the absorption of sunlight. Understanding the VOC composition of the atmosphere relies on our ability to characterize all of their possible reaction pathways. When considering reactions of (transient) VOCs with sunlight, the availability of photolysis rate constants, utilized in general atmospheric models, is often out of experimental reach due to the unstable nature of these molecules. Here, we show how recent advances in computational photochemistry allow us to calculate \textit{in silico} the different ingredients of a photolysis rate constant, namely the photoabsorption cross-section and wavelength-dependent quantum yields. The rich photochemistry of tert-butyl hydroperoxide, for which experimental data is available, is employed to test our protocol and highlight the strengths and weaknesses of different levels of electronic structure and nonadiabatic molecular dynamics to study the photochemistry of (transient) VOCs

Tracking the ultraviolet-induced photochemistry of thiophenone during and after ultrafast ring opening

Photoinduced isomerization reactions lie at the heart of many chemical processes in nature. The mechanisms of such reactions are determined by a delicate interplay of coupled electronic and nuclear dynamics occurring on the femtosecond scale, followed by the slower redistribution of energy into different vibrational degrees of freedom. Here we apply time-resolved photoelectron spectroscopy with a seeded extreme ultraviolet free-electron laser to trace the ultrafast ring opening of gas-phase thiophenone molecules following ultraviolet photoexcitation. When combined with ab initio electronic structure and molecular dynamics calculations of the excited- and ground-state molecules, the results provide insights into both the electronic and nuclear dynamics of this fundamental class of reactions. The initial ring opening and non-adiabatic coupling to the electronic ground state are shown to be driven by ballistic S–C bond extension and to be complete within 350 fs. Theory and experiment also enable visualization of the rich ground-state dynamics that involve the formation of, and interconversion between, ring-opened isomers and the cyclic structure, as well as fragmentation over much longer timescales

Direct Measurements of Colloidal Solvophoresis under Imposed Solvent and Solute Gradients

We describe a microfluidic system that enables direct visualization and measurement of diffusiophoretic migration of colloids in response to imposed solution gradients. Such measurements have proven difficult or impossible in macroscopic systems due to difficulties in establishing solution gradients that are sufficiently strong yet hydrodynamically stable. We validate the system with measurements of the concentration-dependent diffusiophoretic mobility of polystyrene colloids in NaCl gradients, confirming that diffusiophoretic migration velocities are proportional to gradients in the logarithm of electrolyte concentration. We then perform the first direct measurement of the concentration-dependent "solvophoretic" mobility of colloids in ethanol-water gradients, whose dependence on concentration and gradient strength was not known either theoretically or experimentally, but which our measurements reveal to be proportional to the gradient in the logarithm of ethanol mole fraction. Finally, we examine solvophoretic migration under a variety of qualitatively distinct chemical gradients, including solvents that are miscible or have finite solubility with water, an electrolyte for which diffusiophoresis proceeds down concentration gradients (unlike for most electrolytes), and a nonelectrolyte (sugar). Our technique enables the direct characterization of diffusiophoretic mobilities of various colloids under various solvent and solute gradients, analogous to the electrophoretic ζ-potential measurements that are routinely used to characterize suspensions. We anticipate that such measurements will provide the feedback required to test and develop theories for solvophoretic and diffusiophoretic migration and ultimately to the conceptual design and engineering of particles that respond in a desired way to their chemical environments

On the Nature of Geometric and Topological Phases in the Presence of Conical Intersections

The observable nature of topological phases related to conical intersections in molecules is studied. Topological phases should be ubiquitous in molecular processes, but their elusive character has often made them a topic of discussion. To shed some light on this issue, we simulate the dynamics governed by a Jahn–Teller Hamiltonian and analyze it employing two theoretical representations of the molecular wave function: the adiabatic and the exact factorization. We find fundamental differences between effects related to topological phases arising exclusively in the adiabatic representation, and thus not related to any physical observable, and geometric phases within the exact factorization that can be connected to an observable quantity. We stress that while the topological phase of the adiabatic representation is an intrinsic property of the Hamiltonian, the geometric phase of the exact factorization depends on the dynamics that the system undergoes and is connected to the circulation of the nuclear momentum field

Dynamics near a conical intersection - a diabolical compromise for the approximations of ab initio multiple spawning

Full multiple spawning (FMS) offers an exciting framework for the development of strategies to simulate the excited-state dynamics of molecular systems. FMS proposes to depict the dynamics of nuclear wavepackets by using a growing set of traveling multidimensional Gaussian functions called trajectory basis functions (TBFs). Perhaps the most recognized method emanating from FMS is the so-called ab initio multiple spawning (AIMS). In AIMS, the couplings between TBFs—in principle exact in FMS—are approximated to allow for the on-the-fly evaluation of required electronic-structure quantities. In addition, AIMS proposes to neglect the so-called second-order nonadiabatic couplings and the diagonal Born–Oppenheimer corrections. While AIMS has been applied successfully to simulate the nonadiabatic dynamics of numerous complex molecules, the direct influence of these missing or approximated terms on the nonadiabatic dynamics when approaching and crossing a conical intersection remains unknown to date. It is also unclear how AIMS could incorporate geometric-phase effects in the vicinity of a conical intersection. In this work, we assess the performance of AIMS in describing the nonadiabatic dynamics through a conical intersection for three two-dimensional, two-state systems that mimic the excited-state dynamics of bis(methylene)adamantyl, butatriene cation, and pyrazine. The population traces and nuclear density dynamics are compared with numerically exact quantum dynamics and trajectory surface hopping results. We find that AIMS offers a qualitatively correct description of the dynamics through a conical intersection for the three model systems. However, any attempt at improving the AIMS results by accounting for the originally neglected second-order nonadiabatic contributions appears to be stymied by the hermiticity requirement of the AIMS Hamiltonian and the independent first-generation approximation

A Molecular Perspective on Tully Models for Nonadiabatic Dynamics

Over the past decades, an important number of methods have been developed to simulate the nonadiabatic dynamics of molecules, that is, the dynamics of molecules beyond the Born–Oppenheimer approximation. These nonadiabatic methods differ in the way they approximate the dynamics emanating from the time-dependent molecular Schrödinger equation. In 1990, Tully devised a series of three one-dimensional model systems to test the approximations of the method called trajectory surface hopping. The Tully models were designed to probe different scenarios of nonadiabatic processes, such as single and multiple nonadiabatic (re)crossings. These one-dimensional models rapidly became the testbed for any new nonadiabatic dynamics strategy. In this work, we present a molecular perspective to the Tully models by highlighting a correspondence between these simple one-dimensional models and processes happening during the excited-state dynamics of molecules. More importantly, each of these nonadiabatic processes can be connected to a given exemplary molecular system, and we propose here three molecules that could serve as molecular Tully models, reproducing some of the key features of the original models but this time in a high-dimensional space. We compare trajectory surface hopping with the ab initio multiple spawning for the three molecular Tully models and highlight particular features and differences between these methods resulting from their distinct approximations. We also provide all the necessary information – initial conditions and all required parameters for the dynamics as well as the electronic structure – employed in our simulations such that the molecular Tully models can become in the future a unified and standardized test for ab initio nonadiabatic molecular dynamics methods. The molecular Tully models also offer an exciting link between the world of low-dimensional model systems for nonadiabatic dynamics and the excited-state dynamics of molecular systems in their full dimensionality