Trajectory surface-hopping study of methane photodissociation dynamics

Unraveling the Photodissociation Branching and Pathways of Methane at 118 nm by Imaging the CH3, CH2, and CH Fragments

Under irradiation of a vacuum ultraviolet (VUV) photon, methane dissociates and yields multiple fragments. This photochemical behavior is not only of fundamental importance but also with wide-ranging implications in several branches of science. Despite that and numerous previous investigations, the product channel branching is still under debate, and the underlying dissociation mechanisms remain elusive. In this study, the photofragment imaging technique was exploited for the first time to map out the momentum and anisotropy parameter distributions of the CH3, CH2, and CH fragments at the 118 nm photolysis wavelength (10.48 eV photon energy). In conjunction with previously reported results of the H atom fragment at 121.6 nm (10.2 eV), a complete set of product channel branching in both two-body and three-body fragmentations is accurately determined. We concluded that extensive nonadiabatic transitions partake in the processes with two-body fragmentations accounting for more than 90% of overall photodissociation, for which the channel branching values for CH2 + H2 and CH3 + H are about 0.66 and 0.25, respectively. Careful kinematic analysis enables us to untangle the intertwined triple fragmentations into the CH2(X̃ 3B1 and ã 1A1) + H + H and CH(X2Π) + H + H2 channels and to evidence their underlying sequential (or stepwise) mechanisms. With the aid of electronic correlation and prior theoretical calculations of the potential energy surfaces, the most probable or dominant dissociation pathways are elucidated. Comparisons with fragmentary reports in the literature on various photochemical aspects are also documented, and discrepancies are clarified. This comprehensive study benchmarks the VUV photochemistry of methane and advances our understanding of this important process.

Density functionals for core excitations

The core excitation energies and related principal ionization energies are obtained for selected molecules using several density functionals and compared with benchmark equation-of-motion coupled cluster (EOM-CC) results. Both time-dependent and time-independent formulations of excitation spectra in the time-dependent density functional theory and the EOM-CC are employed to obtain excited states that are not always easily accessible with the time-independent method. Among those functionals, we find that the QTP(00) functional, which is only parameterized to reproduce the five IPs of water, provides excellent core IPs and core excitation energies, consistently yielding better excitation and ionization energies. We show that orbital eigenvalues of KS density functional theory play an important role in determining the accuracy of the excitation and photoelectron spectra.

Low-temperature gas-phase formation of indene in the interstellar medium

Polycyclic aromatic hydrocarbons (PAHs) are fundamental molecular building blocks of fullerenes and carbonaceous nanostructures in the interstellar medium and in combustion systems. However, an understanding of the formation of aromatic molecules carrying five-membered rings-the essential building block of nonplanar PAHs-is still in its infancy. Exploiting crossed molecular beam experiments augmented by electronic structure calculations and astrochemical modeling, we reveal an unusual pathway leading to the formation of indene (C₉H₈)-the prototype aromatic molecule with a five-membered ring-via a barrierless bimolecular reaction involving the simplest organic radical-methylidyne (CH)-and styrene (C₆H₅C₂H₃) through the hitherto elusive methylidyne addition-cyclization-aromatization (MACA) mechanism. Through extensive structural reorganization of the carbon backbone, the incorporation of a five-membered ring may eventually lead to three-dimensional PAHs such as corannulene (C₂₀H₁₀) along with fullerenes (C₆₀, C₇₀), thus offering a new concept on the low-temperature chemistry of carbon in our galaxy.

Ultraviolet photochemistry of ethane: implications for the atmospheric chemistry of the gas giants

Chemical processing in the stratospheres of the gas giants is driven by incident vacuum ultraviolet (VUV) light. Ethane is an important constituent in the atmospheres of the gas giants in our solar system. The present work describes translational spectroscopy studies of the VUV photochemistry of ethane using tuneable radiation in the wavelength range 112 ≤ λ ≤ 126 nm from a free electron laser and event-triggered, fast-framing, multi-mass imaging detection methods. Contributions from at least five primary photofragmentation pathways yielding CH2, CH3 and/or H atom products are demonstrated and interpreted in terms of unimolecular decay following rapid non-adiabatic coupling to the ground state potential energy surface. These data serve to highlight parallels with methane photochemistry and limitations in contemporary models of the photoinduced stratospheric chemistry of the gas giants. The work identifies additional photochemical reactions that require incorporation into next generation extraterrestrial atmospheric chemistry models which should help rationalise hitherto unexplained aspects of the atmospheric ethane/acetylene ratios revealed by the Cassini-Huygens fly-by of Jupiter.

Photoinduced C–H bond fission in prototypical organic molecules and radicals

We survey and assess current knowledge regarding the primary photochemistry of hydrocarbon molecules and radicals.

Quantum transition probabilities during a perturbing pulse: Differences between the nonadiabatic results and Fermi's golden rule forms

For a perturbed quantum system initially in the ground state, the coefficient c_k(t) of excited state k in the time-dependent wave function separates into adiabatic and nonadiabatic terms. The adiabatic term a_k(t) accounts for the adjustment of the original ground state to form the new ground state of the instantaneous Hamiltonian H(t), by incorporating excited states of the unperturbed Hamiltonian H₀ without transitions; a_k(t) follows the adiabatic theorem of Born and Fock. The nonadiabatic term b_k(t) describes excitation into another quantum state k; b_k(t) is obtained as an integral containing the time derivative of the perturbation. The true transition probability is given by b_k(t) ², as first stated by Landau and Lifshitz. In this work, we contrast b_k(t) ² and c_k(t) ². The latter is the norm-square of the entire excited-state coefficient which is used for the transition probability within Fermi's golden rule. Calculations are performed for a perturbing pulse consisting of a cosine or sine wave in a Gaussian envelope. When the transition frequency ω_k0 is on resonance with the frequency ω of the cosine wave, b_k(t) ² and c_k(t) ² rise almost monotonically to the same final value; the two are intertwined, but they are out of phase with each other. Off resonance (when ω_k0 ≠ ω), b_k(t) ² and c_k(t) ² differ significantly during the pulse. They oscillate out of phase and reach different maxima but then fall off to equal final values after the pulse has ended, when a_k(t) ≡ 0. If ω_k0 < ω, b_k(t) ² generally exceeds c_k(t) ², while the opposite is true when ω_k0 > ω. While the transition probability is rising, the midpoints between successive maxima and minima fit Gaussian functions of the form a exp[-b(t - d)²]. To our knowledge, this is the first analysis of nonadiabatic transition probabilities during a perturbing pulse.

Coherence in nonradiative transitions: internal conversion in Rydberg-excited N-methyl and N-ethyl morpholine

The survival of coherent wavepacket motion during internal conversions is observed in relatively large molecules, N-methyl morpholine and N-ethyl morpholine, where standard models imply fast decoherence in a statistical limit.

Methane ice photochemistry and kinetic study using laser desorption time-of-flight mass spectrometry at 20 K

Laser Desorption Post-Ionization Time-Of-Flight Mass Spectrometry is used to perform a systematic kinetic study on the pure methane photolysis in the condensed phase at 20 K and provides for the first time effective rate constants and branching ratios for primary processes leading to CH₃, CH₂, and CH radicals upon irradiation by VUV light in the 120–170 nm domain.

Improved catalytic activity of rhodium monolayer modified nickel (110) surface for the methane dehydrogenation reaction: a first-principles study

The catalytic activity of pure Ni (110) and single Rh layer deposited Ni (110) surface for the complete dehydrogenation of methane is theoretically investigated by means of gradient-corrected periodic density functional theory. A detailed kinetic study, based on the analysis of the optimal reaction pathway for the transformation of CH4 to C and H through four elementary steps (CH4 → CH3 + H; CH3 → CH2 + H; CH2 → CH + H; CH → C + H) is presented for pure Ni (110) and Rh/Ni (110) surfaces and compared with pure Rh (110) surface. Through systematic examination of adsorbed geometries and transition states, we show that single layer deposition of Rh on Ni (110) surface has a striking influence on lowering the activation energy barrier of the dehydrogenation reaction. Moreover, it is found that a pure Ni (110) surface has a tendency for carbon deposition on the catalytic surface during the methane dissociation reaction which decreases the stability of the catalyst. However, the deposition of carbon is largely suppressed by the addition of a Rh overlayer on the pure Ni (110) surface. The physical origin of stronger chemisorption of carbon on Ni (110) relative to Rh/Ni (110) has been elucidated by getting insight into the electronic structures and d-band model of the catalytic surfaces. Considering the balance in both the catalytic activity as well as the catalyst stability, we propose that the Rh/Ni (110) surface possesses much improved catalytic property compared to pure Ni (110) and pure Rh (110) surfaces.

Photodissociation dynamics of propargylene, HCCCH

We report a joint theoretical and experimental study on the photodissociation of the C₃H₂ isomer propargylene, HCCCH, combining velocity map imaging with nonadiabatic surface hopping calculations.

Photoinduced ultrafast Wolff rearrangement: a non-adiabatic dynamics perspective

One reaction, two routes: full-dimensional non-adiabatic dynamics simulations shed light on the ultrafast photoinduced Wolff rearrangement in an α-diazocarbonyl compound. The trajectories show both concerted asynchronous and stepwise processes leading to the corresponding ketene.

Adiabatic and nonadiabatic contributions to the energy of a system subject to a time-dependent perturbation: complete separation and physical interpretation

When a time-dependent perturbation acts on a quantum system that is initially in the nondegenerate ground state ∣0> of an unperturbed Hamiltonian H(0), the wave function acquires excited-state components ∣k> with coefficients c(k)(t) exp(-iE(k)t/ℏ), where E(k) denotes the energy of the unperturbed state ∣k>. It is well known that each coefficient c(k)(t) separates into an adiabatic term a(k)(t) that reflects the adjustment of the ground state to the perturbation--without actual transitions--and a nonadiabatic term b(k)(t) that yields the probability amplitude for a transition to the excited state. In this work, we prove that the energy at any time t also separates completely into adiabatic and nonadiabatic components, after accounting for the secular and normalization terms that appear in the solution of the time-dependent Schrödinger equation via Dirac's method of variation of constants. This result is derived explicitly through third order in the perturbation. We prove that the cross-terms between the adiabatic and nonadiabatic parts of c(k)(t) vanish, when the energy at time t is determined as an expectation value. The adiabatic term in the energy is identical to the total energy obtained from static perturbation theory, for a system exposed to the instantaneous perturbation λH'(t). The nonadiabatic term is a sum over excited states ∣k> of the transition probability multiplied by the transition energy. By evaluating the probabilities of transition to the excited eigenstates ∣k'(t)> of the instantaneous Hamiltonian H(t), we provide a physically transparent explanation of the result for E(t). To lowest order in the perturbation parameter λ, the probability of finding the system in state ∣k'(t)> is given by λ(2) ∣b(k)(t)∣(2). At third order, the transition probability depends on a second-order transition coefficient, derived in this work. We indicate expected differences between the results for transition probabilities obtained from this work and from Fermi's golden rule.

Nitrogen incorporation in CH(4)-N(2) photochemical aerosol produced by far ultraviolet irradiation

Nitrile incorporation into Titan aerosol accompanying hydrocarbon chemistry is thought to be driven by extreme UV wavelengths (λ<120 nm) or magnetospheric electrons in the outer reaches of the atmosphere. Far UV radiation (120-200 nm), which is transmitted down to the stratosphere of Titan, is expected to affect hydrocarbon chemistry only and not initiate the formation of nitrogenated species. We examined the chemical properties of photochemical aerosol produced at far UV wavelengths, using a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS), which allows for elemental analysis of particle-phase products. Our results show that aerosol formed from CH(4)/N(2) photochemistry contains a surprising amount of nitrogen, up to 16% by mass, a result of photolysis in the far UV. The proportion of nitrogenated organics to hydrocarbon species is shown to be correlated with that of N(2) in the irradiated gas. The aerosol mass greatly decreases when N(2) is removed, which indicates that N(2) plays a major role in aerosol production. Because direct dissociation of N(2) is highly improbable given the immeasurably low cross section at the wavelengths studied, the chemical activation of N(2) must occur via another pathway. Any chemical activation of N(2) at wavelengths >120 nm is presently unaccounted for in atmospheric photochemical models. We suggest that reaction with CH radicals produced from CH(4) photolysis may provide a mechanism for incorporating N into the molecular structure of the aerosol. Further work is needed to understand the chemistry involved, as these processes may have significant implications for how we view prebiotic chemistry on early Earth and similar planets.