Theoretical investigation of Jahn–Teller and pseudo-Jahn–Teller interactions in the ammonia cation

The ultrafast vibronic dynamics of ammonia's D̃ state

Using vacuum-ultraviolet time-resolved velocity map imaging of photoelectrons, we study ultrafast coupled electronic and nuclear dynamics in low-lying Rydberg states of ammonia. Vibrationally-resolved internal vibrational relaxation (IVR) is observed in a progression of the e' bending modes. This vibrational progression is only observed in the D̃ state, and is lost upon ultrafast internal conversion to the C̃ and B̃ electronic states. Due to the ultrashort time scale of the internal conversion (ca. 64 fs), and the vibronic resolution, the non-adiabatic coupling vectors are identified and verified with ab initio calculations. The time-scale of this IVR process is highly surprising and significant because IVR is usually treated as an incoherent process that proceeds statistically, according to a "Fermi's Golden Rule"-like model, where the process scales with the available degrees of freedom. Here, we show that it can be highly non-statistical, restricted to only a very small subset of vibrational motions.

Laser-induced electron diffraction of the ultrafast umbrella motion in ammonia

Visualizing molecular transformations in real-time requires a structural retrieval method with Ångström spatial and femtosecond temporal atomic resolution. Imaging of hydrogen-containing molecules additionally requires an imaging method sensitive to the atomic positions of hydrogen nuclei, with most methods possessing relatively low sensitivity to hydrogen scattering. Laser-induced electron diffraction (LIED) is a table-top technique that can image ultrafast structural changes of gas-phase polyatomic molecules with sub-Ångström and femtosecond spatiotemporal resolution together with relatively high sensitivity to hydrogen scattering. Here, we image the umbrella motion of an isolated ammonia molecule (NH₃) following its strong-field ionization. Upon ionization of a neutral ammonia molecule, the ammonia cation (NH₃ ⁺) undergoes an ultrafast geometrical transformation from a pyramidal ( Φ HNH = 107 ° ) to planar ( Φ HNH = 120 ° ) structure in approximately 8 femtoseconds. Using LIED, we retrieve a near-planar ( Φ HNH = 117   ±   5 ° ) field-dressed NH₃ ⁺ molecular structure 7.8 - 9.8 femtoseconds after ionization. Our measured field-dressed NH₃ ⁺ structure is in excellent agreement with our calculated equilibrium field-dressed structure using quantum chemical ab initio calculations.

Simulation of Vibronic Spectra of Flexible Systems: Hybrid DVR-Harmonic Approaches

Our general framework for the simulation of vibrational signatures in electronic spectra has been extended to treat one large-amplitude motion (LAM) at the anharmonic level, coupled to the other small-amplitude motions (SAM) treated as harmonic. The coupling between LAM and SAM is minimized thanks to the use of delocalized internal coordinates, which are built automatically from the molecular topology. General LAMs can be employed, ranging from intrinsic reaction coordinates to rigid or flexible paths based on the distinguished coordinate approach. The anharmonic model is based on a fully numerical method based on the discrete variable representation (DVR) theory, supporting different types of boundary conditions. The inclusion of this model in a general-purpose electronic structure code makes available to the user a large panel of quantum chemistry models, for both isolated and condensed phases. The flexibility and reliability of the new framework are illustrated by some case studies, covering various types of LAMs, ranging from a small test case, the photoelectron spectrum of ammonia, to larger systems, such as phenylanthracene and cyclobutanone.

Revisiting the (E + A) ⊗ (e + a) problems of polyatomic systems with trigonal symmetry: general expansions of their vibronic Hamiltonians

In this work, we derive general expansions in vibrational coordinates for the (E + A) ⊗ (e + a) vibronic Hamiltonians of molecules with one and only one C₃ axis.

Nonadiabatic nuclear dynamics of the ammonia cation studied by surface hopping classical trajectory calculations

The Landau-Zener (LZ) type classical-trajectory surface-hopping algorithm is applied to the nonadiabatic nuclear dynamics of the ammonia cation after photoionization of the ground-state neutral molecule to the excited states of the cation. The algorithm employs a recently proposed formula for nonadiabatic LZ transition probabilities derived from the adiabatic potential energy surfaces. The evolution of the populations of the ground state and the two lowest excited adiabatic states is calculated up to 200 fs. The results agree well with quantum simulations available for the first 100 fs based on the same potential energy surfaces. Three different time scales are detected for the nuclear dynamics: Ultrafast Jahn-Teller dynamics between the excited states on a 5 fs time scale; fast transitions between the excited state and the ground state within a time scale of 20 fs; and relatively slow partial conversion of a first-excited-state population to the ground state within a time scale of 100 fs. Beyond 100 fs, the adiabatic electronic populations are nearly constant due to a dynamic equilibrium between the three states. The ultrafast nonradiative decay of the excited-state populations provides a qualitative explanation of the experimental evidence that the ammonia cation is nonfluorescent.

The Jahn-Teller effect in the electron momentum spectroscopy of ammonia

The 1e and 3a(1) bands of the ammonia molecule have been studied using the high-resolution electron momentum spectroscopy at impact energies of 1200 and 600 eV. Several slices of 1e and 3a(1) bands in the different binding energy ranges were selected, and their electron-momentum distributions were carefully compared. The discernable difference among the distributions of the selected slices of the 1e band shows that the Jahn-Teller effect indeed influences the electron momentum distribution of the 1e orbital of ammonia.

Efficient computation of adiabatic populations in multi-mode Jahn-Teller systems through the use of effective vibrational modes

A highly efficient scheme for computing adiabatic electronic populations in multi-mode Jahn-Teller systems is presented. It relies on the transformation to an effective-mode vibrational basis in which the relevant quantities depend on the coordinates of a single mode only. In this way, the generally tedious numerical evaluation of high-dimensional integrals is avoided and replaced by one-dimensional integrations. The effective-mode scheme is applied to a variety of two-mode and three-mode Jahn-Teller systems and gives a typical speedup of about two to three orders of magnitude as compared to the direct evaluation of the adiabatic populations. The gain grows rapidly with the numbers of modes.

Imaging studies of the photodissociation of NH3+ and ND3+ cations

Velocity map ion imaging methods have been used to study the photofragmentation dynamics of state-selected NH3+ and ND3+ cations. The cations were prepared in selected nu2+ bending vibrational levels of the ground (x2A'') electronic state by two-photon resonant, three-photon ionization of NH3(ND3), via several different nu2' levels of the and ' Rydberg states. Subsequent excitation to the A2E state by absorption of a 207.6 nm photon resulted in N-H(D) bond fission and NH2+(ND2+) fragment ion formation. These fragments exhibit isotropic recoil velocity distributions, which peak at low kinetic energy but extend to the maximum allowed by energy conservation. Such findings accord with conclusions from earlier electron induced photoionization and photoelectron-photoion coincidence studies of NH3 at similar total energies (defined relative to the ground-state neutral) and, as previously, can be rationalized in terms of excitation to the Jahn-Teller distorted state, rapid radiationless transfer via one or more conical intersections linking the and state potential energy surfaces (PESs) and subsequent unimolecular decay on the latter PES. Weak NH2+ and NH+ fragment ion signals are also observed when exciting with the ionization laser only; imaging these fragment ions provides some insights into their likely formation mechanisms.

Photoionization-induced dynamics of ammonia: ab initio potential energy surfaces and time-dependent wave packet calculations for the ammonia cation

An analytical anharmonic six-dimensional three-sheeted potential energy surface for the ground and first excited states of the ammonia cation has been developed which is tailored to model the ultrafast photoinduced dynamics. Selected ab initio cuts, obtained by multireference configuration interaction calculations, have been used to determine the parameters of a diabatic representation for this Jahn-Teller and pseudo-Jahn-Teller system. The model includes higher-order coupling terms both for the Jahn-Teller and for the pseudo-Jahn-Teller matrix elements. The relaxation to the ground state is possible via dynamical pseudo-Jahn-Teller couplings involving the asymmetric bending and stretching coordinates. The photoelectron spectrum of NH3 and the internal conversion dynamics of NH3+ have been determined by wave packet propagation calculations employing the multiconfigurational time-dependent Hartree method. Three different time scales are found in the dynamics calculations for the second absorption band. The ultrafast Jahn-Teller dynamics of the two excited states occurs on a 5 fs time scale. The major part of the internal conversion to the ground state takes place within a short time scale of 20 fs. This fast internal conversion is, however, incomplete and the remaining excited state population does not decay completely even within 100 fs.

Inversion vibration of PH3+(X2A2") studied by zero kinetic energy photoelectron spectroscopy

We report the first rotationally resolved spectroscopic studies on PH3+(X2A2") using zero kinetic energy photoelectron spectroscopy and coherent VUV radiation. The spectra about 8000 cm(-1) above the ground vibrational state of PH3+(X2A2") have been recorded. We observed the vibrational energy level splittings of PH3+(X2A2") due to the tunneling effect in the inversion (symmetric bending) vibration (nu2+). The energy splitting for the first inversion vibrational state (0+/0-) is 5.8 cm(-1). The inversion vibrational energy levels, rotational constants, and adiabatic ionization energies (IEs) for nu2+ = 0-16 have been determined. The bond angles between the neighboring P-H bonds and the P-H bond lengths are also obtained using the experimentally determined rotational constants. With the increasing of the inversion vibrational excitations (nu2+), the bond lengths (P-H) increase a little and the bond angles (H-P-H) decrease a lot. The inversion vibrational energy levels have also been calculated by using one dimensional potential model and the results are in good agreement with the experimental data for the first several vibrational levels. In addition to inversion vibration, we also observed firstly the other two vibrational modes: the symmetric P-H stretching vibration (nu1+) and the degenerate bending vibration (nu4+). The fundamental frequencies for nu1+ and nu4+ are 2461.6 (+/-2) and 1043.9 (+/-2) cm(-1), respectively. The first IE for PH3 was determined as 79670.9 (+/-1) cm(-1).

Multistate multimode vibronic dynamics: Entanglement of electronic and vibrational degrees of freedom in the benzene radical cation

An earlier theoretical treatment of multimode and multistate vibronic coupling in the benzene radical cation [Koppel et al., J. Chem. Phys. 117, 2657 (2002)] is extended to investigate also the behavior of the nuclear degrees of freedom and to include additional electronic states. The five lowest doublet electronic states are considered which have been shown earlier to be all interconnected through a series of conical intersections of their potential-energy surfaces. In the most extensive calculations, they are all included simultaneously in the quantum dynamical calculations performed, which represent a system of unprecedented complexity treated in this way. The results are compared with various reduced-dimensionality treatments (i.e., employing reduced vibrational and electronic function spaces). The different temporal behavior of the various electronic populations is emphasized and traced to the different locations of the various seams of conical intersections: due to the coherent oscillations of the time-dependent wave packet this leads to an oscillatory behavior in some cases and to monotonous behavior in others. A seemingly irreversible behavior of the system dynamics in this strictly microscopic treatment is confirmed. The importance of this benchmark system to highlight complex, entangled multimode, and multistate vibronic dynamics is pointed out.

Exploring the Jahn-Teller and pseudo-Jahn-Teller conical intersections in the ethane radical cation

We report a theoretical account on the static and dynamic aspects of the Jahn-Teller (JT) and pseudo-Jahn-Teller (PJT) interactions in the ground and first excited electronic states of the ethane radical cation. The findings are compared with the experimental photoionization spectrum of ethane. The present theoretical approach is based on a model diabatic Hamiltonian and with the parameters derived from ab initio calculations. The optimized geometry of ethane in its electronic ground state (1A1g) revealed an equilibrium staggered conformation belonging to the D3d symmetry point group. At the vertical configuration, the ethane radical cation belongs to this symmetry point group. The ground and low-lying electronic states of this radical cation are of 2Eg, 2A1g, 2Eu, and 2A2u symmetries. Elementary symmetry selection rule suggests that the degenerate electronic states of the radical cation are prone to the JT distortion when perturbed along the degenerate vibrational modes of eg symmetry. The 2A1g state is estimated to be approximately 0.345 eV above the 2Eg state and approximately 2.405 eV below the 2Eu state at the vertical configuration. The symmetry selection rule also suggests PJT crossings of the 2A1g and the 2Eg electronic states of the radical cation along the vibrational modes of eg symmetry and such crossings appear to be energetically favorable also. The irregular vibrational progressions, with numerous shoulders and small peaks, observed below 12.55 eV in the experimental recording are manifestations of the dynamic (E x e)-JT effect. Our findings revealed that the PJT activity of the degenerate vibrational modes is particularly strong in the 2Eg-2A1g electronic manifold which leads to a broad and diffuse structure of the observed photoelectron band.