Do Molecular Geometries Change Under Vibrational Strong Coupling?

As pioneering experiments have shown, strong coupling between molecular vibrations and light modes in an optical cavity can significantly alter molecular properties and even affect chemical reactivity. However, the current theoretical description is limited and far from complete. To explore the origin of this exciting observation, we investigate how the molecular structure changes under strong light-matter coupling using an ab initio method based on the cavity Born-Oppenheimer Hartree-Fock ansatz. By optimizing H2O and H2O2 resonantly coupled to cavity modes, we study the importance of reorientation and geometric relaxation. In addition, we show that the inclusion of one or two cavity modes can change the observed results. On the basis of our findings, we derive a simple concept to estimate the effect of the cavity interaction on the molecular geometry using the molecular polarizability and the dipole moments.

Unraveling a Cavity-Induced Molecular Polarization Mechanism from Collective Vibrational Strong Coupling.

We demonstrate that collective vibrational strong coupling of molecules in thermal equilibrium can give rise to significant local electronic polarizations in the thermodynamic limit. We do so by first showing that the full nonrelativistic Pauli-Fierz problem of an ensemble of strongly coupled molecules in the dilute-gas limit reduces in the cavity Born-Oppenheimer approximation to a cavity-Hartree equation for the electronic structure. Consequently, each individual molecule experiences a self-consistent coupling to the dipoles of all other molecules, which amount to non-negligible values in the thermodynamic limit (large ensembles). Thus, collective vibrational strong coupling can alter individual molecules strongly for localized "hotspots" within the ensemble. Moreover, the discovered cavity-induced polarization pattern possesses a zero net polarization, which resembles a continuous form of a spin glass (or better polarization glass). Our findings suggest that the thorough understanding of polaritonic chemistry, requires a self-consistent treatment of dressed electronic structure, which can give rise to numerous, so far overlooked, physical mechanisms.

Electronic Fingerprint of the Protonated Imidazole Dimer Probed by X-ray Absorption Spectroscopy.

Protons in low-barrier superstrong hydrogen bonds are typically delocalized between two electronegative atoms. Conventional methods to characterize such superstrong hydrogen bonds are vibrational spectroscopy and diffraction techniques. We introduce soft X-ray spectroscopy to uncover the electronic fingerprints for proton sharing in the protonated imidazole dimer, a prototypical building block enabling effective proton transport in biology and high-temperature fuel cells. Using nitrogen core excitations as a sensitive probe for the protonation status, we identify the X-ray signature of a shared proton in the solvated imidazole dimer in a combined experimental and theoretical approach. The degree of proton sharing is examined as a function of structural variations that modify the shape of the low-barrier potential in the superstrong hydrogen bond. We conclude by showing how the sensitivity to the quantum distribution of proton motion in the double-well potential is reflected in the spectral signature of the shared proton.

Ab Initio Vibro-Polaritonic Spectra in Strongly Coupled Cavity-Molecule Systems.

Recent experiments have revealed the profound effect of strong light-matter interactions in optical cavities on the electronic ground state of molecular systems. This phenomenon, known as vibrational strong coupling, can modify reaction rates and induce the formation of molecular vibrational polaritons, hybrid states involving both photon modes, and vibrational modes of molecules. We present an ab initio methodology based on the cavity Born-Oppenheimer Hartree-Fock ansatz, which is specifically powerful for ensembles of molecules, to calculate vibro-polaritonic IR spectra. This method allows for a comprehensive analysis of these hybrid states. Our semiclassical approach, validated against full quantum simulations, reproduces key features of the vibro-polaritonic spectra. The underlying analytic gradients also allow for optimization of cavity-coupled molecular systems and performing semiclassical dynamics simulations.

Following the Nonadiabatic Ultrafast Dynamics of Uracil via Simulated X-ray Absorption Spectra.

The nucleobase uracil exhibits high photostability due to ultrafast relaxation processes mediated by conical intersections (CoIns), where the interplay between nuclear and electron dynamics becomes crucial. In our previous study, we observed seemingly long-lived traces of electronic coherence for the relaxation process through the S2/S1 CoIn by applying our ansatz for coupled nuclear and electron dynamics in molecules (NEMol). In this work, we theoretically investigate how time-dependent transient X-ray absorption spectroscopy can be used to observe this ultrafast dynamics. Therefore, we calculated X-ray absorption spectra (XAS) for the oxygen K-edge, using a multireference protocol in combination with NEMol dynamics. Thus, we have access to both the transient XAS based on the nuclear wavepacket dynamics and the modulation of the signals caused by the electronic coherence induced by the excitation process and the presence of a CoIn seam. In both cases, we were able to qualitatively predict its influence on the resulting XAS.

Cavity Born-Oppenheimer Hartree-Fock Ansatz: Light-Matter Properties of Strongly Coupled Molecular Ensembles.

Experimental studies indicate that optical cavities can affect chemical reactions through either vibrational or electronic strong coupling and the quantized cavity modes. However, the current understanding of the interplay between molecules and confined light modes is incomplete. Accurate theoretical models that take into account intermolecular interactions to describe ensembles are therefore essential to understand the mechanisms governing polaritonic chemistry. We present an ab initio Hartree-Fock ansatz in the framework of the cavity Born-Oppenheimer approximation and study molecules strongly interacting with an optical cavity. This ansatz provides a nonperturbative, self-consistent description of strongly coupled molecular ensembles, taking into account the cavity-mediated dipole self-energy contributions. To demonstrate the capability of the cavity Born-Oppenheimer Hartree-Fock ansatz, we study the collective effects in ensembles of strongly coupled diatomic hydrogen fluoride molecules. Our results highlight the importance of the cavity-mediated intermolecular dipole-dipole interactions, which lead to energetic changes of individual molecules in the coupled ensemble.

Nonadiabatic Wave Packet Dynamics with Ab Initio Cavity-Born-Oppenheimer Potential Energy Surfaces.

Strong coupling of molecules with quantized electromagnetic fields can reshape their potential energy surfaces by forming dressed states. In such a scenario, it is possible to manipulate the dynamics of the molecule and open new photochemical reaction pathways. A theoretical approach to describe such coupled molecular-photon systems is the Cavity-Born-Oppenheimer (CBO) approximation. Similarly to the standard Born-Oppenheimer (BO) approximation, the system is partitioned and the electronic part of the system is treated quantum mechanically. This separation leads to CBO surfaces that depend on both nuclear and photonic coordinates. In this work, we demonstrated, for two molecular examples, how the concept of the CBO approximation can be used to perform nonadiabatic wave packet dynamics of a coupled molecular-cavity system. The light-matter interaction is incorporated in the CBO surfaces and the associated nonadiabatic coupling elements. We show that molecular and cavity contributions can be treated on the same numerical footing. This approach gives a new perspective on the description of light-matter coupling in molecular systems.

Time-resolved X-ray and XUV based spectroscopic methods for nonadiabatic processes in photochemistry.

The photochemistry of numerous molecular systems is influenced by conical intersections (CIs). These omnipresent nonadiabatic phenomena provide ultra-fast radiationless relaxation channels by creating degeneracies between electronic states and decide over the final photoproducts. In their presence, the Born-Oppenheimer approximation breaks down, and the timescales of the electron and nuclear dynamics become comparable. Due to the ultra-fast dynamics and the complex interplay between nuclear and electronic degrees of freedom, the direct experimental observation of nonadiabatic processes close to CIs remains challenging. In this article, we give a theoretical perspective on novel spectroscopic techniques capable of observing clear signatures of CIs. We discuss methods that are based on ultra-short laser pulses in the extreme ultraviolet and X-ray regime, as their spectral and temporal resolution allow for resolving the ultra-fast dynamics near CIs.

Ultrafast strong-field dissociation of vinyl bromide: An attosecond transient absorption spectroscopy and non-adiabatic molecular dynamics study.

Attosecond extreme ultraviolet (XUV) and soft x-ray sources provide powerful new tools for studying ultrafast molecular dynamics with atomic, state, and charge specificity. In this report, we employ attosecond transient absorption spectroscopy (ATAS) to follow strong-field-initiated dynamics in vinyl bromide. Probing the Br M edge allows one to assess the competing processes in neutral and ionized molecular species. Using ab initio non-adiabatic molecular dynamics, we simulate the neutral and cationic dynamics resulting from the interaction of the molecule with the strong field. Based on the dynamics results, the corresponding time-dependent XUV transient absorption spectra are calculated by applying high-level multi-reference methods. The state-resolved analysis obtained through the simulated dynamics and related spectral contributions enables a detailed and quantitative comparison with the experimental data. The main outcome of the interaction with the strong field is unambiguously the population of the first three cationic states, D 1, D 2, and D 3. The first two show exclusively vibrational dynamics while the D 3 state is characterized by an ultrafast dissociation of the molecule via C-Br bond rupture within 100 fs in 50% of the analyzed trajectories. The combination of the three simulated ionic transient absorption spectra is in excellent agreement with the experimental results. This work establishes ATAS in combination with high-level multi-reference simulations as a spectroscopic technique capable of resolving coupled non-adiabatic electronic-nuclear dynamics in photoexcited molecules with sub-femtosecond resolution.

Coupled nuclear and electron dynamics in the vicinity of a conical intersection.

Ultrafast optical techniques allow us to study ultrafast molecular dynamics involving both nuclear and electronic motion. To support interpretation, theoretical approaches are needed that can describe both the nuclear and electron dynamics. Hence, we revisit and expand our ansatz for the coupled description of the nuclear and electron dynamics in molecular systems (NEMol). In this purely quantum mechanical ansatz, the quantum-dynamical description of the nuclear motion is combined with the calculation of the electron dynamics in the eigenfunction basis. The NEMol ansatz is applied to simulate the coupled dynamics of the molecule NO2 in the vicinity of a conical intersection (CoIn) with a special focus on the coherent electron dynamics induced by the non-adiabatic coupling. Furthermore, we aim to control the dynamics of the system when passing the CoIn. The control scheme relies on the carrier envelope phase of a few-cycle IR pulse. The laser pulse influences both the movement of the nuclei and the electrons during the population transfer through the CoIn.

Waveform control of molecular dynamics close to a conical intersection.

Conical intersections are ubiquitous in chemical systems but, nevertheless, extraordinary points on the molecular potential energy landscape. They provide ultra-fast radiationless relaxation channels, their topography influences the product branching, and they equalize the timescales of the electron and nuclear dynamics. These properties reveal optical control possibilities in the few femtosecond regime. In this theoretical study, we aim to explore control options that rely on the carrier envelope phase of a few-cycle IR pulse. The laser interaction creates an electronic superposition just before the wave packet reaches the conical intersection. The imprinted phase information is varied by the carrier envelope phase to influence the branching ratio after the conical intersection. We test and analyze this scenario in detail for a model system and show to what extent it is possible to transfer this type of control to a realistic system like uracil.

Visualizing conical intersection passages via vibronic coherence maps generated by stimulated ultrafast X-ray Raman signals.

The rates and outcomes of virtually all photophysical and photochemical processes are determined by conical intersections. These are regions of degeneracy between electronic states on the nuclear landscape of molecules where electrons and nuclei evolve on comparable timescales and thus become strongly coupled, enabling radiationless relaxation channels upon optical excitation. Due to their ultrafast nature and vast complexity, monitoring conical intersections experimentally is an open challenge. We present a simulation study on the ultrafast photorelaxation of uracil, based on a quantum description of the nuclei. We demonstrate an additional window into conical intersections obtained by recording the transient wavepacket coherence during this passage with an X-ray free-electron laser pulse. Two major findings are reported. First, we find that the vibronic coherence at the conical intersection lives for several hundred femtoseconds and can be measured during this entire time. Second, the time-dependent energy-splitting landscape of the participating vibrational and electronic states is directly extracted from Wigner spectrograms of the signal. These offer a physical picture of the quantum conical intersection pathways through visualizing their transient vibronic coherence distributions. The path of a nuclear wavepacket in the vicinity of the conical intersection is directly mapped by the proposed experiment.

Photoprotecting Uracil by Coupling with Lossy Nanocavities.

We analyze how the photorelaxation dynamics of a molecule can be controlled by modifying its electromagnetic environment using a nanocavity mode. In particular, we consider the photorelaxation of the RNA nucleobase uracil, which is the natural mechanism to prevent photodamage. In our theoretical work, we identify the operative conditions in which strong coupling with the cavity mode can open an efficient photoprotective channel, resulting in a relaxation dynamics twice as fast as the natural one. We rely on a state-of-the-art chemically detailed molecular model and a non-Hermitian Hamiltonian propagation approach to perform full-quantum simulations of the system dissipative dynamics. By focusing on the photon decay, our analysis unveils the active role played by cavity-induced dissipative processes in modifying chemical reaction rates, in the context of molecular polaritonics. Remarkably, we find that the photorelaxation efficiency is maximized when an optimal trade-off between light-matter coupling strength and photon decay rate is satisfied. This result is in contrast with the common intuition that increasing the quality factor of nanocavities and plasmonic devices improves their performance. Finally, we use a detailed model of a metal nanoparticle to show that the speedup of the uracil relaxation could be observed via coupling with a nanosphere pseudomode, without requiring the implementation of complex nanophotonic structures.

Intersystem Crossing as a Key Component of the Nonadiabatic Relaxation Dynamics of Bithiophene and Terthiophene.

We present a nonadiabatic dynamics study concerning the subpicosecond relaxation of excited states in dimeric and trimeric thiophene chains. The influence of the triplet states in the overall process is, for the first time, taken into account by explicitly including spin-orbit couplings and hence allowing intersystem crossing phenomena. We observe the fundamental role of the triplet state manifold in driving the full relaxation process. In particular we evidence the effect of both, inter-ring rotation and ring-opening, in the process, as compared to the monomer, where the ring-opening process appears as the dominant one. In addition, the evolution of the open structures allows for trans to cis isomerization in the dimer and trimer. The overall relaxation process slows down with chain elongation. The complex decay mechanism characterized by the presence of different competing channels, due to the presence of a quasi degenerate manifold, is explained allowing the rationalization of oligothiophenes photophysics.

Ab initio molecular dynamics of thiophene: the interplay of internal conversion and intersystem crossing.

The fast and slow components of the relaxation of photoexcited thiophene have been investigated by means of SHARC (surface hopping including arbitrary couplings) molecular dynamics based on multiconfiguration electronic structure calculations. Triplet states are included to ascertain their role in the relaxation process. After thiophene is excited to the S1 state, ultrafast dynamics (τfast = 96 fs) initiates a ring opening due to cleavage of a carbon sulfur bond and simultaneous ring puckering. This time constant is in agreement with previous experimental and theoretical studies. The subsequent dynamics of the open-ring structures is characterized by the interplay of internal conversion and intersystem crossing. For the open-ring structures, the S0, S1, T1 and T2 states are nearly degenerate and the spin-orbit couplings are large. The underlying potential energy surface is flat and long-lived open-ring structures in the singlet as well as in the triplet states are formed. Both the participation of triplet states and the shape of the energy surface explain the experimentally observed slow ring closure in the ground state.

Record Broken: A Copper Peroxide Complex with Enhanced Stability and Faster Hydroxylation Catalysis.

Tyrosinase model systems pinpoint pathways to translating Nature's synthetic abilities for useful synthetic catalysts. Mostly, they use N-donor ligands which mimic the histidine residues coordinating the two copper centres. Copper complexes with bis(pyrazolyl)methanes with pyridinyl or imidazolyl moieties are already reported as excellent tyrosinase models. Substitution of the pyridinyl donor results in the new ligand HC(3-tBuPz)2 (4-CO2 MePy) which stabilises a room-temperature stable µ-η2 :η2 -peroxide dicopper(II) species upon oxygenation. It reveals highly efficient catalytic activity as it hydroxylates 8-hydroxyquinoline in high yields (TONs of up to 20) and much faster than all other model systems (max. conversion within 7.5 min). Stoichiometric reactions with para-substituted sodium phenolates show saturation kinetics which are nearly linear for electron-rich substrates. The resulting Hammett correlation proves the electrophilic aromatic substitution mechanism. Furthermore, density functional theory (DFT) calculations elucidate the influence of the substituent at the pyridinyl donor: the carboxymethyl group adjusts the basicity and nucleophilicity without additional steric demand. This substitution opens up new pathways in reactivity tuning.

Deactivation pathways of thiophene and oligothiophenes: internal conversion versus intersystem crossing.

Oligothiophenes and polythiophenes are building blocks of organic-based energy conversion materials. Therefore the lifetime of the excited states plays a central role. As a first step to understand the factors influencing the performance, we investigated the deactivation processes from the first excited state S1 of thiophene and small oligothiophenes containing up to four rings using quantum chemical calculations. For thiophene a low-lying S1/S0 conical intersection seam is easily accessible and drives the fast internal conversion. In oligothiophenes barriers inhibit this passage while deactivation pathways via intersystem crossing channels open. The first one is responsible for the high triplet quantum yields and takes place shortly after the Franck-Condon region. The second one occurs in the vicinity of a local S1 minimum. The calculated spin-orbit coupling strength together with the singlet-triplet energy gaps can explain the decreasing triplet and increasing fluorescence quantum yields for growing chain length. From the triplets the ground state is reachable by inter-ring torsions and T1/S0 intersections. The present results allow a deeper understanding of the deactivation pathways of thiophene and small oligothiophenes and are of potential interest for the photophysics of longer oligothiophenes and polythiophenes used in optical devices.