State-selective predissociation dynamics of methylamines: The vibronic and H∕D effects on the conical intersection dynamics

N-H and N-C Bond Dissociation Pathways in Ultraviolet Photodissociation of Dimethylamine

We investigated the interlinked N-H and N-C photochemistry of primary and secondary amines via the state-resolved detection of vibrationally excited CH3 product and H atom product by 200-235 nm dimethylamine photodissociation using resonance-enhanced multiphoton ionization (REMPI) and velocity map imaging (VMI) techniques. The out-of-plane bending (ν2) vibrationally excited CH3 showed a bimodal translational energy distribution that became unimodal with a near-zero product yield at longer photolysis wavelengths (λphotolysis). In contrast, a unimodal distribution was observed for the C-H stretching (νCH) vibrationally excited CH3 products with an almost constant product yield in the examined λphotolysis region. We ascribed the state-specific energy releases of the CH3 products to two reaction pathways based on calculations of the potential energy surface (PES): the direct N-CH3 dissociation pathway and the indirect N-CH3 dissociation pathway via the N-H bond conical intersection. Meanwhile, the H atom product showed a bimodal energy distribution similar to the ammonia photodissociation model, with an excited-state counterproduct channel that became accessible at a shorter λphotolysis. These results suggest that the N-H and N-C bond dissociations are connected, and these dissociations cause different photochemistry between primary/secondary amines and tertiary amines.

Hybrid Density Functional Valence Bond Method with Multistate Treatment

Recently, a hybrid density functional valence bond (VB) method, λ-DFVB(U), has been proposed and shown to give accuracy that is comparable to that of CASPT2 in calculations of atomization energies, atomic excitation energies, and reaction barriers, while its computational cost is approximately the same as the valence bond self-consistent-field (VBSCF) method. However, the interaction between electronic states is not included in λ-DFVB(U) since the last step of λ-DFVB(U) is not a diagonalization of the Hamiltonian matrix on the electronic state basis. Therefore, λ-DFVB(U) gives the wrong topology of the potential energy surfaces (PESs) near the conical intersection region. In the present paper, we propose a novel hybrid density functional VB method with multistate treatment, named λ-DFVB(MS), in which an effective Hamiltonian matrix is constructed on the basis of the diabatic states obtained by the valence-bond-based compression approach for the diabatization scheme, and the interaction between electronic states can be included through the diagonalization of the effective Hamiltonian matrix. Test calculations show that λ-DFVB(MS) gives the correct topology of the PESs near the conical intersection region. We also show that the VBSCF wave function with selected VB structures can be applied as a reference in λ-DFVB(MS).

Photodissociation dynamics of methylamine in the blue edge of the A-band. II. The NH2 + CH3 channel

The photodissociation dynamics leading to the C-N bond cleavage in methylamine (CH3NH2) are investigated upon photoexcitation in the blue edge of the first absorption A-band, in the 198-204 nm range. Velocity map images of the generated methyl (CH3) fragment detected in specific vibrational modes, i.e., ν = 0, ν1 = 1, and ν2 = 1, through resonance enhanced multiphoton ionization, are presented along with the corresponding translational energy distributions and the angular analysis. The experimental results are complemented by high-level ab initio calculations of potential energy curves as a function of the C-N bond distance. While a similar single Boltzmann-type contribution is observed in all the translational energy distributions measured, the speed-dependent anisotropy parameter obtained through the angular analysis reveals the presence of two different mechanisms. Prompt dissociation through the conical intersection between the Ã1A' first excited state and the ground state located in the exit channel is, indeed, revealed as a minor channel. In contrast, slow dissociation on the ground state, presumably from frustrated N-H bond cleavage trajectories, constitutes the major reaction pathway leading to the methyl formation.

Photodissociation dynamics of methylamine in the blue edge of the A-band. I. The H-atom elimination channel

The photodissociation dynamics of methylamine (CH3NH2) upon excitation in the blue edge of the first absorption A-band, in the 198-203 nm range, are investigated by means of nanosecond pump-probe laser pulses and velocity map imaging combined with H(2S)-atom detection through resonance enhanced multiphoton ionization. The images and corresponding translational energy distributions for the H-atoms produced show three different contributions associated with three reaction pathways. The experimental results are complemented by high-level ab initio calculations. The potential energy curves computed as a function of the N-H and C-H bond distances allow us to draw a picture of the different mechanisms. Major dissociation occurs through N-H bond cleavage and it is triggered by an initial geometrical change, i.e., from a pyramidal configuration of the C-NH2 with respect to the N atom to a planar geometry. The molecule is then driven into a conical intersection (CI) seam where three outcomes can take place: first, threshold dissociation into the second dissociation limit, associated with the formation of CH3NH(Ã), is observed; second, direct dissociation after passage through the CI leading to the formation of ground state products; and third, internal conversion into the ground state well in advance to dissociation. While the two last pathways were previously reported at a variety of wavelengths in the 203-240 nm range, the former had not been observed before to the best of our knowledge. The role of the CI and the presence of an exit barrier in the excited state, which modify the dynamics leading the two last mechanisms, are discussed considering the different excitation energies used.

Linearized Pair-Density Functional Theory

Multiconfiguration pair-density functional theory (MC-PDFT) is a post-SCF multireference method that has been successful at computing ground- and excited-state energies. However, MC-PDFT is a single-state method in which the final MC-PDFT energies do not come from diagonalization of a model-space Hamiltonian matrix, and this can lead to inaccurate topologies of potential energy surfaces near locally avoided crossings and conical intersections. Therefore, in order to perform physically correct ab initio molecular dynamics with electronically excited states or to treat Jahn-Teller instabilities, it is necessary to develop a PDFT method that recovers the correct topology throughout the entire nuclear configuration space. Here we construct an effective Hamiltonian operator, called the linearized PDFT (L-PDFT) Hamiltonian, by expanding the MC-PDFT energy expression to first order in a Taylor series of the wave function density. Diagonalization of the L-PDFT Hamiltonian gives the correct potential energy surface topology near conical intersections and locally avoided crossings for a variety of challenging cases including phenol, methylamine, and the spiro cation. Furthermore, L-PDFT outperforms MC-PDFT and previous multistate PDFT methods for predicting vertical excitations from a variety of representative organic chromophores.

Photochemical carbon-sulfur bond cleavage in thioethers mediated via excited state Rydberg-to-valence evolution

Time-resolved photoelectron imaging and supporting ab initio quantum chemistry calculations were used to investigate non-adiabatic excess energy redistribution dynamics operating in the saturated thioethers diethylsulfide, tetrahydrothiophene and thietane. In all cases, 200 nm excitation leads to molecular fragmentation on an ultrafast (<100 fs) timescale, driven by the evolution of Rydberg-to-valence orbital character along the S-C stretching coordinate. The C-S-C bending angle was also found to be a key coordinate driving initial internal conversion through the excited state Rydberg manifold, although only small angular displacements away from the ground state equilibrium geometry are required. Conformational constraints imposed by the cyclic ring structures of tetrahydrothiophene and thietane do not therefore influence dynamical timescales to any significant extent. Through use of a high-intensity 267 nm probe, we were also able to detect the presence of some transient (bi)radical species. These are extremely short lived, but they appear to confirm the presence of two competing excited state fragmentation channels - one proceeding directly from the initially prepared 4p manifold, and one involving non-adiabatic population of the 4s state. This is in addition to a decay pathway leading back to the S₀ electronic ground state, which shows an enhanced propensity in the 5-membered ring system tetrahydrothiophene over the other two species investigated.

Compressed-State Multistate Pair-Density Functional Theory

Multiconfiguration pair-density functional theory (MC-PDFT) is a multireference method that can be used to calculate excited states. However, MC-PDFT potential energy surfaces have the wrong topology at conical intersections because the last step of MC-PDFT is not a diagonalization of a model-space Hamiltonian matrix, as done in, for example, multistate second-order perturbation theory (MS-CASPT2). We have previously proposed methods that solve this problem by diagonalizing a model-space effective Hamiltonian matrix, where the diagonal elements are MC-PDFT energies for intermediate states, and the off-diagonal elements are evaluated by wave function theory. One previous method is called variational multistate PDFT (VMS-PDFT), whose intermediate states maximize the trace of the effective Hamiltonian, namely, the sum of the MC-PDFT energies of the model-space states; the VMS-PDFT is very robust but is more computationally expensive than another method, extended multistate PDFT (XMS-PDFT), in which the transformation to intermediate states is accomplished without needing any density functional evaluations. However, although VMS-PDFT was accurate in all cases tested, XMS-PDFT was accurate in only some of them. In the present paper, we propose a new method, called compressed-state multistate PDFT (CMS-PDFT), that is as efficient as XMS-PDFT and as accurate as VMS-PDFT. The new method maximizes the trace of the classical Coulomb energy of the intermediate states such that the electron densities of the intermediate states are compressed. We show that CMS-PDFT performs robustly even where XMS-PDFT fails.

Multi-state pair-density functional theory

Multi-configuration pair-density functional theory (MC-PDFT) has previously been applied successfully to carry out ground-state and excited-state calculations. However, because they include no interaction between electronic states, MC-PDFT calculations in which each state's PDFT energy is calculated separately can give an unphysical double crossing of potential energy surfaces (PESs) in a region near a conical intersection. We have recently proposed state-interaction pair-density functional theory (SI-PDFT) to treat nearly degenerate states by creating a set of intermediate states with state interaction; although this method is successful, it is inconvenient because two SCF calculations and two sets of orbitals are required and because it puts the ground state on an unequal footing with the excited states. Here we propose two new methods, called extended-multi-state-PDFT (XMS-PDFT) and variational-multi-state-PDFT (VMS-PDFT), that generate the intermediate states in a balanced way with a single set of orbitals. The former uses the intermediate states proposed by Granovsky for extended multi-configuration quasi-degenerate perturbation theory (XMC-QDPT); the latter obtains the intermediate states by maximizing the sum of the MC-PDFT energies for the intermediate states. We also propose a Fourier series expansion to make the variational optimizations of the VMS-PDFT method convenient, and we implement this method (FMS-PDFT) both for conventional configuration-interaction solvers and for density-matrix-renormalization-group solvers. The new methods are tested for eight systems, exhibiting avoided crossings among two to six states. The FMS-PDFT method is successful for all cases for which it has been tested (all cases in this paper except O3 for which it was not tested), and XMS-PDFT is successful for all eight cases except the mixed-valence case. Since both XMS-PDFT and VMS-PDFT are less expensive than XMS-CASPT2, they will allow well-correlated calculations on much larger systems for which perturbation theory is unaffordable.

State-interaction pair density functional theory for locally avoided crossings of potential energy surfaces in methylamine

State-interaction pair-density functional theory agrees well with extended MS-CASPT2 in regions of strong state coupling near conical intersections.

Control of Nonadiabatic Passage through a Conical Intersection by a Dynamic Resonance

Nonadiabatic processes, dominated by dynamic passage of reactive fluxes through conical intersections (CIs), are considered to be appealing means for manipulating reaction paths, particularly via initial vibrational preparation. Nevertheless, obtaining direct experimental evidence of whether specific-mode excitation affects the passage at the CI is challenging, requiring well-resolved time- or frequency-domain experiments. Here promotion of methylamine-d2 (CH3ND2) molecules to spectral-resolved rovibronic states on the excited S1 potential energy surface, coupled to sensitive D photofragment probing, allowed us to follow the N-D bond fission dynamics. The branching ratios between slow and fast D photofragments and the internal energies of the CH3ND(X̃) photofragments confirm correlated anomalies for predissociation initiated from specific rovibronic states. These anomalies reflect the existence of a dynamic resonance that strongly depends on the energy of the initially excited rovibronic states, the evolving vibrational mode on the repulsive S1 part during N-D bond elongation, and the manipulated passage through the CI that leads to CH3ND radicals excited with C-N-D bending. This resonance plays an important role in the bifurcation dynamics at the CI and can be foreseen to exist in other photoinitiated processes and to control their outcome.

Evidence for quantum effects in the predissociation of methylamine isotopologues

The H product distributions obtained from the predissociation of methylamine isotopologues are extremely sensitive to the energy difference between the initially prepared vibrational states and the conical intersections and not only to the nature of the pre-excited nuclear motions.

Time-resolved photoionization spectroscopy of mixed Rydberg-valence states: indole case study

Time-resolved photoelectron imaging reveals subtle new mechanistic insight into the ultraviolet relaxation dynamics of gas-phase indole.

Vibronic structures and dynamics of the predissociating dimethyl sulfide and its isotopomers (CH₃SCH₃, CD₃SCD₃, CH₃SCD₃) at the conical intersection

Conical intersection seam comprised of crossing surfaces of two lowest excited states of dimethyl sulfide (DMS) has been directly accessed by the one-photon excitation from the ground equilibrium state. Since the S-C bond rupture takes place promptly, the molecular structure on the excited state effectively belongs to C(S) symmetry. Namely, excited states of 1(1)B1 and 1(1)A2 in C(2)V become 1(1)A'' and 2(1)A'' states in C(S), respectively, and the optical transition from the ground equilibrium state to the dissociating molecule at the conical intersection seam is symmetry-allowed to facilitate the nonadiabatic transition on the 2(1)A'' state, leading eventually to the CH3S + CH3 products. The dynamic study of DMS, in this sense, gives the great opportunity to unravel the vibronic structure of the conical intersection seam by the conventional one-photon excitation method. In this work, utilizing the photofragment excitation (PHOFEX) spectroscopic method, the vibronic structures of DMS and its isotope analogs (CD3SCD3, CH3SCD3) at the conical intersection seam have been revealed, providing accurate lifetimes and detailed dynamics associated with individual vibronic transitions. The lifetime of the excited DMS is estimated to be ~100 fs, indicating that the dissociation is complete within one single oscillation in the conical intersection region. It is also found that the symmetric CSC stretching mode is strongly coupled to the reaction coordinate, as manifested by our experimental finding that the fragmentation yield of the S-CD3 bond is enhanced compared to that of the S-CH3 bond in the CH3SCD3 dissociation reaction only when the CSC symmetric stretching vibrational mode is excited at the conical intersection region. This work demonstrates that the better understanding of the excited state could make the bond-selective chemistry into reality.

Ultrafast dynamics of aniline following 269-238 nm excitation and the role of the S2(pi3s/pi sigma) state

Femtosecond time-resolved photoelectron imaging is employed to investigate ultrafast electronic relaxation in aniline, a prototypical aromatic amine. The molecule is excited at wavelengths between 269 and 238 nm. We observe that the S2(pi3s/pi sigma*) state is populated directly during the excitation process at all wavelengths and that the population bifurcates to two decay pathways. One of these involves ultrafast relaxation from the Rydberg component of S2(pi3s/pi sigma*) to the S1(Pi Pi)* state, from which it relaxes back to the electronic ground state on a much longer timescale. The other appears to involve motion along the pi sigma* dissociative potential energy surface. At higher excitation energies, the dominant excitation is to the S3(pi pi*) state, which undergoes extremely efficient electronic relaxation back to the ground state. Our study supports some conclusions reached from H-atom photofragment translational spectroscopy measurements and pump-probe photoionization measurements and contradicts some others.

Shedding new light on the role of the Rydberg state in the photochemistry of aniline

Efficient electronic relaxation following the absorption of ultraviolet light is crucial for the photostability of biological chromophores, so understanding the microscopic details of the decay pathways is of considerable interest. Here, we employ femtosecond time-resolved photoelectron imaging to investigate the ultrafast intramolecular dynamics of aniline, a prototypical aromatic amine, following excitation just below the second absorption maximum. We find that both the second ππ* state and the Rydberg state are populated during the excitation process. Surprisingly, the dominant non-radiative decay pathway is an ultrafast relaxation mechanism that transfers population straight back to the electronic ground-state. The vibrational energy resolution and photoelectron angular distributions obtained in our experiments reveal an interesting bifurcation of the Rydberg population to two non-radiative decay channels. The existence of these competing non-radiative relaxation channels in aniline illustrates how its photostability arises from a subtle balance between dynamics on different electronically excited states and importantly between Rydberg and valence states.

Observation of NH X(3)Σ(-) as a Primary Product of Methylamine Photodissociation: Evidence of Roaming-Mediated Intersystem Crossing?

3+1 Resonance-enhanced multiphoton ionization and photofragment excitation spectroscopy have been used to identify NH X(3)Σ(-) as a primary product of methylamine photodissociation after state-specific excitation to the S1 state. On the basis of standard thermochemical data, NH X(3)Σ(-) can be formed only in conjunction with closed-shell CH4 coproducts, indicating that dissociation must occur on the T1 surface. It is proposed that the mechanism for the formation of triplet NH and CH4 involves intramolecular abstraction between frustrated radical products and is an example of roaming-mediated intersystem crossing.