Photochemistry of the water dimer: Time-dependent quantum wave-packet description of the dynamics at the S1-S0 conical intersection

The "Hole" Story in Ionized Water from the Perspective of Ehrenfest Dynamics

The radiolysis of liquid water and the radiation-matter interactions that happen in aqueous environments are important to the fields of chemistry, materials, and environmental sciences, as well as the biological and physiological response to extreme conditions and medical treatments. The initial stage of radiolysis is the ultrafast response, or hole dynamics, that triggers chemical processes within complex energetic landscapes that may include reactivity. A fundamental understanding necessitates the use of theoretical methods that are capable of simulating both ultrafast coherence and non-adiabatic energy transfer pathways. In this work, we carry out an ab initio Ehrenfest dynamics study to provide a more complete description of the ultrafast dynamics and reactive events initiated by photoionization of water. After sudden ionization, a range of processes, including hole trapping and transfer, large OH oscillations, proton transfer and subsequent relay, formation of the metastable Zundel complex, and long-lived coherence, are identified and new insight into their driving forces is elucidated.

The Role of Excited-State Proton Relays in the Photochemical Dynamics of Water Nanodroplets

In this work, we applied nonadiabatic excited-state molecular dynamics in tandem with ab initio electronic structure theory to illustrate a complete mechanistic landscape underpinning the ultraviolet absorption-initiated photochemical dynamics in water nanodroplets. The goal is to understand the nonequilibrium excited-state molecular dynamics initiated by the relaxation of a solvated photoelectron and consequential photochemical processes. The lowest-lying excited state shows the proton dissociation for a single water molecule forming intermediate hydronium complexes through a proton relay. At approximately 100 fs, the proton relay process gives rise to the relaxation of the excited state accompanied by a rapid increase in the nonadiabatic coupling strength with the ground state, and the nanodroplet nonradiatively decays. The nonadiabatic transition to the ground state produces excited vibrational states that facilitate the recombination of the dissociated proton and hydroxyl group, eventually leading to the desorption of water molecules from the nanodroplet. Additionally, lifetimes of transient photochemical events are also resolved for the relaxation of a solvated electron, excited-state proton relay, and nonradiative transition.

Excitation Energy Transfer and Low-Efficiency Photolytic Splitting of Water Ice by Vacuum UV Light

Experimental estimates of photolytic efficiency (yield per photon) for photodissociation and photodesorption from water ice range from about 10(-3) to 10(-1). However, in the case of photodissociation of water in the gas phase, it is close to unity. Exciton dynamics carried out by a quantum mechanical time-dependent propagator shows that in the eight most stable water hexamers, the excitation diffuses away from the initially excited molecule within a few femtoseconds. On the basis of these quantum dynamics simulations, it is hypothesized that the ultrafast exciton energy transfer process, which in general gives rise to a delocalized exciton within these clusters, may contribute to the low efficiency of photolytic processes in water ice. It is proposed that exciton diffusion inherently competes with the nuclear dynamics that drives the photodissociation process in the repulsive S1 state on the sub-10 fs time scale.

Matrix-free application of Hamiltonian operators in Coifman wavelet bases

A means of evaluating the action of Hamiltonian operators on functions expanded in orthogonal compact support wavelet bases is developed, avoiding the direct construction and storage of operator matrices that complicate extension to coupled multidimensional quantum applications. Application of a potential energy operator is accomplished by simple multiplication of the two sets of expansion coefficients without any convolution. The errors of this coefficient product approximation are quantified and lead to use of particular generalized coiflet bases, derived here, that maximize the number of moment conditions satisfied by the scaling function. This is at the expense of the number of vanishing moments of the wavelet function (approximation order), which appears to be a disadvantage but is shown surmountable. In particular, application of the kinetic energy operator, which is accomplished through the use of one-dimensional (1D) [or at most two-dimensional (2D)] differentiation filters, then degrades in accuracy if the standard choice is made. However, it is determined that use of high-order finite-difference filters yields strongly reduced absolute errors. Eigensolvers that ordinarily use only matrix-vector multiplications, such as the Lanczos algorithm, can then be used with this more efficient procedure. Applications are made to anharmonic vibrational problems: a 1D Morse oscillator, a 2D model of proton transfer, and three-dimensional vibrations of nitrosyl chloride on a global potential energy surface.