Taming Disulfide Bonds with Laser Fields. Nonadiabatic Surface-Hopping Simulations in a Ruthenium Complex

Laser control of chemical reactions is a challenging field of research. In particular, the theoretical description of coupled electronic and nuclear motion in the presence of laser fields is not a trivial task and simulations are mostly restricted to small systems or molecules treated within reduced dimensionality. Here, we demonstrate how the excited state dynamics of [Ru­(S–Sbpy)­(bpy)2]2+ can be controlled using explicit laser fields in the context of fewest-switches surface hopping. In particular, the transient properties along the excited state dynamics leading to population of the T1 minimum energy structure are exploited to define simple laser fields capable of slowing and even completely stopping the onset of S–S bond dissociation. The use of a linear vibronic coupling model to parametrize the potential energy surfaces showcases the strength of the surface-hopping methodology to study systems including explicit laser fields using many nuclear degrees of freedom and a large amount of close-lying electronic excited states

Vibrational Sampling and Solvent Effects on the Electronic Structure of the Absorption Spectrum of 2‑Nitronaphthalene

The influence of vibrational motion on electronic excited state properties is investigated for the organic chromophore 2-nitronaphtalene in methanol. Specifically, the performance of two vibrational sampling techniques – Wigner sampling and sampling from an ab initio molecular dynamics trajectory– is assessed, in combination with implicit and explicit solvent models. The effects of the different sampling/solvent combinations on the energy and electronic character of the absorption bands are analyzed in terms of charge transfer and exciton size, computed from the electronic transition density. The absorption spectra obtained using sampling techniques and its underlying properties are compared to those of the electronic excited states calculated at the Franck–Condon equilibrium geometry. It is found that the absorption bands of the vibrational ensembles are red-shifted compared to the Franck–Condon bright states, and this red-shift scales with the displacement from the equilibrium geometry. Such displacements are found larger and better described when using ensembles from the harmonic Wigner distribution than snapshots from the molecular dynamics trajectory. Particularly relevant is the torsional motion of the nitro group that quenches the charge transfer character of some of the absorption bands. This motion, however, is better described in the molecular dynamics trajectory. Thus, none of the vibrational sampling approaches can satisfactorily capture all important aspects of the nuclear motion. The inclusion of solvent also red-shifts the absorption bands with respect to the gas phase. This red-shift scales with the charge-transfer character of the bands and is found larger for the implicit than for the explicit solvent model. The advantages and drawbacks of the different sampling and solvent models are discussed to guide future research on the calculation of UV–vis spectra of nitroaromatic compounds

Surface hopping dynamics on vibronic coupling models

CONSPECTUS: The simulation of photoinduced non-adiabatic dynamics is of great relevance in many scientific disciplines, ranging from physics and materials science to chemistry and biology. Upon light irradiation, different relaxation processes take place in which electronic and nuclear motion are intimately coupled. These are best described by the time-dependent molecular Schrödinger equation, but its solution poses fundamental practical challenges to contemporary theoretical chemistry. Two widely used and complementary approaches to this problem are multiconfigurational time-dependent Hartree (MCTDH) and trajectory surface hopping (SH). MCTDH is an accurate fully quantum-mechanical technique but often is feasible only in reduced dimensionality, in combination with approximate vibronic coupling (VC) Hamiltonians, or both (i.e., reduced-dimensional VC potentials). In contrast, SH is a quantum−classical technique that neglects most nuclear quantum effects but allows nuclear dynamics in full dimensionality by calculating potential energy surfaces on the fly. If nuclear quantum effects do not play a central role and a linear VC (LVC) Hamiltonian is appropriatee.g., for stiff molecules that generally keep their conformation in the excited statethen it seems advantageous to combine the efficient LVC and SH techniques. In this Account, we describe how surface hopping based on an LVC Hamiltonian (SH/LVC)as recently implemented in the SHARC surface hopping packagecan provide an economical and automated approach to simulate nonadiabatic dynamics. First, we illustrate the potential of SH/LVC in a number of showcases, including intersystem crossing in SO2, intra-Rydberg dynamics in acetone, and several photophysical studies on large transition-metal complexes, which would be much more demanding or impossible to perform with other methods. While all of the applications provide very useful insights into lightinduced phenomena, they also hint at difficulties faced by the SH/LVC methodology that need to be addressed in the future. Second, we contend that the SH/LVC approach can be useful to benchmark SH itself. By the use of the same (LVC) potentials as MCTDH calculations have employed for decades and by relying on the efficiency of SH/LVC, it is possible to directly compare multiple SH test calculations with a MCTDH reference and ponder the accuracy of various correction algorithms behind the SH methodology, such as decoherence corrections or momentum rescaling schemes. Third, we demonstrate how the efficiency of SH/ LVC can also be exploited to identify essential nuclear and electronic degrees of freedom to be employed in more accurate MCTDH calculations. Lastly, we show that SH/LVC is able to advance the development of SH protocols that can describe nuclear dynamics including explicit laser fieldsa very challenging endeavor for trajectory-based schemes. To end, this Account compiles the typical costs of contemporary SH simulations, evidencing the great advantages of using parametrized potentials. The LVC model is a sleeping beauty that, kissed by SH, is fueling the field of excited-state molecular dynamics. We hope that this Account will stimulate future research in this direction, leveraging the advantages of the SH/VC schemes to larger extents and extending their applicability to uncharted territories

Surface hopping dynamics on vibronic coupling models

CONSPECTUS: The simulation of photoinduced non-adiabatic dynamics is of great relevance in many scientific disciplines, ranging from physics and materials science to chemistry and biology. Upon light irradiation, different relaxation processes take place in which electronic and nuclear motion are intimately coupled. These are best described by the time-dependent molecular Schrödinger equation, but its solution poses fundamental practical challenges to contemporary theoretical chemistry. Two widely used and complementary approaches to this problem are multiconfigurational time-dependent Hartree (MCTDH) and trajectory surface hopping (SH). MCTDH is an accurate fully quantum-mechanical technique but often is feasible only in reduced dimensionality, in combination with approximate vibronic coupling (VC) Hamiltonians, or both (i.e., reduced-dimensional VC potentials). In contrast, SH is a quantum−classical technique that neglects most nuclear quantum effects but allows nuclear dynamics in full dimensionality by calculating potential energy surfaces on the fly. If nuclear quantum effects do not play a central role and a linear VC (LVC) Hamiltonian is appropriatee.g., for stiff molecules that generally keep their conformation in the excited statethen it seems advantageous to combine the efficient LVC and SH techniques. In this Account, we describe how surface hopping based on an LVC Hamiltonian (SH/LVC)as recently implemented in the SHARC surface hopping packagecan provide an economical and automated approach to simulate nonadiabatic dynamics. First, we illustrate the potential of SH/LVC in a number of showcases, including intersystem crossing in SO2, intra-Rydberg dynamics in acetone, and several photophysical studies on large transition-metal complexes, which would be much more demanding or impossible to perform with other methods. While all of the applications provide very useful insights into lightinduced phenomena, they also hint at difficulties faced by the SH/LVC methodology that need to be addressed in the future. Second, we contend that the SH/LVC approach can be useful to benchmark SH itself. By the use of the same (LVC) potentials as MCTDH calculations have employed for decades and by relying on the efficiency of SH/LVC, it is possible to directly compare multiple SH test calculations with a MCTDH reference and ponder the accuracy of various correction algorithms behind the SH methodology, such as decoherence corrections or momentum rescaling schemes. Third, we demonstrate how the efficiency of SH/ LVC can also be exploited to identify essential nuclear and electronic degrees of freedom to be employed in more accurate MCTDH calculations. Lastly, we show that SH/LVC is able to advance the development of SH protocols that can describe nuclear dynamics including explicit laser fieldsa very challenging endeavor for trajectory-based schemes. To end, this Account compiles the typical costs of contemporary SH simulations, evidencing the great advantages of using parametrized potentials. The LVC model is a sleeping beauty that, kissed by SH, is fueling the field of excited-state molecular dynamics. We hope that this Account will stimulate future research in this direction, leveraging the advantages of the SH/VC schemes to larger extents and extending their applicability to uncharted territories