Simulations of Condensed Phase Photochemistry: Cage Effect and Internal Conversion in Azoalkanes and Nitrosamines

Investigation of nonadiabatic dynamics in the photolysis of methyl nitrate (CH3ONO2) by on-the-fly surface hopping simulation

The photolysis mechanism of methyl nitrate (CH₃ONO₂) was studied using the on-the-fly surface hopping dynamics at the XMS-CASPT2 level. Several critical geometries, including electronic state minima and conical intersections, were obtained, which play essential roles in the nonadiabatic dynamics of CH₃ONO₂. The ultrafast nonadiabatic decay dynamics to the ground state were simulated, which gives a proper explanation on the broad and structureless absorption spectra of CH₃ONO₂. The photodissociation channels, including CH₃O + NO₂, CH₃O + NO + O, and others, as well as their branching ratios, were identified. When the dynamics starts from the lowest two electronic states (S₁ and S₂), the CH₃O + NO₂ channel is the dominant photolysis pathway, although we observed the minor contributions of other channels. In contrast, when the trajectories start from the third excited state S₃, both CH₃O + NO₂ and CH₃O + NO + O channels become important. Here the CH₃O-NO₂ bond dissociation takes place first, and then for some trajectories, the N-O bond of the NO₂ part breaks successively. The quasi-degeneracy of electronic states may exist in the dissociation limits of both CH₃O + NO₂ and CH₃O + NO + O channels. The current work provides valuable information in the understanding of experimental findings of the wavelength-dependent photolysis mechanism of CH₃ONO₂.

Electronic relaxation of aqueous aminoazobenzenes studied by time-resolved photoelectron spectroscopy and surface hopping TDDFT dynamics calculations

Studies of ultrafast relaxation of molecular chromophores are complicated by the fact that most chromophores of biological and technological importance are rather large molecules and are strongly affected by their environment, either solvent or a protein cage. Here we present an approach which allows us to follow transient electronic structure of complex photoexcited molecules. We use the method of time-resolved photoelectron spectroscopy in solution to follow relaxation of two prototypical aqueous chromophores, Methyl Orange and Metanil Yellow, both of which are aminoazobenzene derivatives. Using excitation by 400 nm laser pulses and ionization by wavelength-selected 46.7 nm XUV pulses from high-order harmonic generation we follow relaxation of both molecules via the dark S1 state. The photoelectron spectra yield binding energies of both ground and excited states. We combine the experimental results with surface hopping time-dependent density functional theory (TDDFT) calculations employing B3LYP+D3 and ωB97X-D functionals. The results demonstrate that the method is generally suitable for description of ultrafast dynamics in these molecules and can recover absolute binding energies observed in the experiment. The B3LYP+D3 functional appears to be better suited for these systems, especially in the case of Metanil Yellow, where it indicates the importance of an intramolecular charge transfer state. Our results pave the way towards quantitative understanding of evolving electronic structure in photo-induced relaxation processes.

Performance of trajectory surface hopping method in the treatment of ultrafast intersystem crossing dynamics

We carried out extensive studies to examine the performance of the fewest-switches surface hopping method in the description of the ultrafast intersystem crossing dynamic of various singlet-triplet (S-T) models by comparison with the results of the exact full quantum dynamics. Different implementation details and some derivative approaches were examined. As expected, it is better to perform the trajectory surface hopping calculations in the spin-adiabatic representation or by the local diabatization approach, instead of in the spin-diabatic representation. The surface hopping method provides reasonable results for the short-time dynamics in the S-T model with weak spin-orbital coupling (diabatic coupling), although it does not perform well in the models with strong spin-orbital coupling (diabatic coupling). When the system accesses the S-T potential energy crossing with rather high kinetic energy, the trajectory surface hopping method tends to produce a good description of the nonadiabatic intersystem crossing dynamics. The impact of the decoherence correction on the performance of the trajectory surface hopping is system dependent. It improves the result accuracy in many cases, while its influence may also be minor for other cases.

On the importance of initial conditions for excited-state dynamics

The vast majority of ab initio excited-state simulations are performed within semiclassical, trajectory-based approaches. Apart from the underlying electronic-structure theory, the reliability of the simulations is controlled by a selection of initial conditions for the classical trajectories. We discuss appropriate choices of initial conditions for simulations of different experimental arrangements: dynamics initiated by continuum-wave (CW) laser fields or triggered by ultrashort laser pulses.

An On-the-Fly Surface-Hopping Program JADE for Nonadiabatic Molecular Dynamics of Polyatomic Systems: Implementation and Applications

Nonadiabatic dynamics simulations have rapidly become an indispensable tool for understanding ultrafast photochemical processes in complex systems. Here, we present our recently developed on-the-fly nonadiabatic dynamics package, JADE, which allows researchers to perform nonadiabatic excited-state dynamics simulations of polyatomic systems at an all-atomic level. The nonadiabatic dynamics is based on Tully's surface-hopping approach. Currently, several electronic structure methods (CIS, TDHF, TDDFT(RPA/TDA), and ADC(2)) are supported, especially TDDFT, aiming at performing nonadiabatic dynamics on medium- to large-sized molecules. The JADE package has been interfaced with several quantum chemistry codes, including Turbomole, Gaussian, and Gamess (US). To consider environmental effects, the Langevin dynamics was introduced as an easy-to-use scheme into the standard surface-hopping dynamics. The JADE package is mainly written in Fortran for greater numerical performance and Python for flexible interface construction, with the intent of providing open-source, easy-to-use, well-modularized, and intuitive software in the field of simulations of photochemical and photophysical processes. To illustrate the possible applications of the JADE package, we present a few applications of excited-state dynamics for various polyatomic systems, such as the methaniminium cation, fullerene (C20), p-dimethylaminobenzonitrile (DMABN) and its primary amino derivative aminobenzonitrile (ABN), and 10-hydroxybenzo[h]quinoline (10-HBQ).

Non-adiabatic dynamics close to conical intersections and the surface hopping perspective

Conical intersections play a major role in the current understanding of electronic de-excitation in polyatomic molecules, and thus in the description of photochemistry and photophysics of molecular systems. This article reviews aspects of the basic theory underlying the description of non-adiabatic transitions at conical intersections, with particular emphasis on the important case when the dynamics of the nuclei are treated classically. Within this classical nuclear motion framework, the main aspects of the surface hopping methodology in the conical intersection context are presented. The emerging picture from this treatment is that of electronic transitions around conical intersections dominated by the interplay of the nuclear velocity and the derivative non-adiabatic coupling vector field.

Vibrational dynamics of carboxylic acid dimers in gas and dilute solution

Ultrafast mid-IR transient absorption spectroscopy has been used to study the vibrational dynamics of hydrogen-bonded cyclic dimers of trifluoroacetic acid and formic acid in both the gas and solution phases (0.05 M in CCl(4)). Ultrafast excitation of the broad O-H cyclic dimer band leads, in the gas phase, to large-scale structural changes of the dimer creating a species with a distinct free O-H stretching band on 20 ps and 200 ps timescales. These timescales are assigned to ring-opening and dissociation of the dimer, respectively. In the solution phase, no such structural rearrangement occurs and our results are consistent with previous studies. The gas phase dynamics are insensitive to both the specific excitation energy (over a span of 550 cm(-1)) and the chemical identity of the dimer.

Femtochemistry of trans-azomethane: a combined experimental and theoretical study

The dissociation dynamics of trans-azomethane upon excitation to the S,(n, pi) state with a total energy of 93 kcalmol-1 is investigated using femtosecond-resolved mass spectrometry in a molecular beam. The transient signal shows an opposite pump-probe excitation feature for the UV (307 nm) and the visible (615 nm) pulses at the perpendicular polarization in comparison with the signal obtained at the parallel polarization: The one-photon symmetry-forbidden process excited by the UV pulse is dominant at the perpendicular polarization, whereas the two-photon symmetry-allowed process initiated by the visible pulse prevails at the parallel polarization. At the perpendicular polarization, we found that the two C-N bonds of the molecule break in a stepwise manner, that is, the first C-N bond breaks in approximately 70 fs followed by the second one in approximately 100 fs, with the intermediate characterized. At the parallel polarization, the first C-N bond cleavage was found to occur in 100 fs with the intensity of the symmetry-allowed transition being one order of magnitude greater than the intensity of the symmetry-forbidden transition at the perpendicular polarization. Theoretical calculations using time-dependent density functional theory (TDDFT) and the complete active space self-consistent field (CASSCF) method have been carried out to characterize the potential energy surface for the ground state, the low-lying excited states, and the cationic ground state at various levels of theory. Combining the experimental and theoretical results, we identified the elementary steps in the mechanism. The initial driving force of the ultrafast bond-breaking process of trans-azomethane (at the perpendicular polarization) is due to the CNNC torsional motion initiated by the vibronic coupling through an intensity-borrowing mechanism for the symmetry-forbidden n-pi transition. Following this torsional motion and the associated molecular symmetry breaking, an S0/S1 conical intersection (CI) can be reached at a torsional angle of 93.1 degrees (predicted at the CASSCF(8.7)/cc-pVDZ level of theory). Funneling through the S0/S1, CI could activate the asymmetric C-N stretching motion, which is the key motion for the consecutive C-N bond breakages on the femtosecond time scale.

On the dissociation of aromatic radical anions in solution

A new theoretical formulation is given for the reaction rate and path for the important reaction class of aromatic radical anion dissociation in solution [Ar-X]-.-->Ar. + X-, and is illustrated for the case of the cyanochlorobenzene radical anion [CN-phi-CI]-. in dimethylformamide. Among the theory's novel features is the inclusion of the conical intersection aspect of this ground electronic state problem, which is key in allowing the reaction to occur and which has a significant impact on the reaction barrier height. Reasonable agreement with the experimental rate is found.

Semiclassical simulations of azomethane photochemistry in the gas phase and in solution

We present semiclassical simulations of the dynamical events which follow the n --> pi excitation of an azomethane molecule: decay to the ground state, isomerization, and dissociation. Ab initio potential energy surfaces and couplings, slightly modified according to available experimental data, are employed in conjunction with a trajectory-surface-hopping method. For an isolated molecule, we show that dissociation takes place in the ground state, because the radiationless decay is very fast. The dissociation is mainly sequential: the methyldiazenyl radical first produced fragmentates very rapidly to N(2) + CH(3)*. The results of previous experiments by Lee's and Zewail's groups, from which partially conflicting conclusions had been drawn, are reproduced and interpreted. Molecular dynamics, coupled to the surface-hopping algorithm, allows us to simulate the photochemistry in acqueous solution. Here the fragmentation is suppressed because of vibrational energy loss to the solvent: only the isomerization takes place.