The non-adiabatic nanoreactor: towards the automated discovery of photochemistry

The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S₀/S₁ and S₁/S₂ conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.

The nonadiabatic nanoreactor is a tool for automated photochemical reaction discovery that extensively explores intersection seams and links conical intersections to photoproduct distributions.     

Nonadiabatic relaxation dynamics of water anion cluster and its isotope effects by ring‐polymer molecular dynamics simulation

We have applied a recently developed hybrid quantum ring‐polymer molecular dynamics method to the nonadiabatic p → s relaxation dynamics in water anion clusters to understand the isotope effects observed in previous experiments. The average relaxation times for (H₂O)₅₀^? and (D₂O)₅₀^? were calculated at 120 and 207 fs, respectively, and are comparable to the experimental results. Therefore, we conclude that nuclear quantum effects play an essential role in understanding the observed isotope effects for water anion cluster nonadiabatic dynamics. The nonadiabatic relaxation mechanisms are also discussed in detail. © 2014 Wiley Periodicals, Inc.     

Research on the rupture of the covalent bond of the dichlorine molecule

Abstract  

Photoinduced reaction of gas-phase dichlorine molecules on the short wavelength side of the A-band gives a negative value for the asymmetry parameter of Cl (² P _1/2) fragments, which conflicts with the intrinsic electronic transition mechanism. In this paper, the dissociation process of the dichlorine molecule has been investigated at numerous excitation wavelengths in the range 310–470 nm using the numerical method and the frontier molecular orbital maps are drawn to obtain insight into the character of the relevant molecular orbitals. The possibilities of radial nonadiabatic transition from the C¹Π_u to the third Ω = 1_u (σ _u* ← σ _g) electronic state are also examined, and found to cause large variations for the angular distribution functions. The wavelength dependence of the beta parameter β ₂(Cl*), which is computed from the partial cross-section with the RZD transition mechanism, agrees with the experimental behavior and further justifies the conclusion that the decrease of beta in the Cl + Cl* channel is because of the radial nonadiabatic interaction between the C and 1_u(III) excited states and this interaction is a key mechanism decreasing the beta parameter value. At last, the kinetic energy distributions of fragments are obtained in the asymptotic region.     

Excitation of Electronic States of N<Subscript>2</Subscript> in Radio-Frequency Electric Field by Electron Impact

Excitation of electronic states of the N₂ molecule by electron impact is recognized as an essential process in nitrogen plasmas that strongly impacts their chemical reactivity and other properties. Many surface and coating technologies are based on radio-frequency plasma discharges in nitrogen. In this paper the electron impact excitation rate coefficients for singlet and triplet electronic states of the N₂ molecule have been calculated in non-equilibrium conditions in the presence of a radio-frequency electric field. A Monte Carlo simulation has been performed in order to determine non-equilibrium electron energy distribution functions within one period of the electric field. By using these distribution functions, the excitation rate coefficients have been obtained in the frequency range from 13.56 up to 500 MHz, at reduced electric field values from 200 to 700 Td.     

A numerical simulation of nonadiabatic electron excitation in the strong field regime: Linear polyenes

Time-dependent Hartree-Fock theory has been used to study of the electronic optical response of a series of linear polyenes in strong laser fields. Ethylene, butadiene, and hexatriene have been calculated with 6-31G(d,p) in the presence of a field corresponding to 8.75 x 10(13) W/cm2 and 760 nm. Time evolution of the electron population indicates not only the pi electrons, but also lower lying valence electrons are involved in electronic response. When the field is aligned with the long axis of the molecule, L?wdin population analysis shows large charges at each end of the molecule. For ethylene, the instantaneous dipole moment followed the field adiabatically, but for hexatriene, nonadiabatic effects were very pronounced. For constant intensity, the nonadiabatic effects in the charge distribution, instantaneous dipole, and orbital populations increased nonlinearly with the length of the polyene. These calculations elucidate the mechanism of the strong field nonadiabatic electron excitation of polyatomic molecules leading to their eventual ionization and fragmentation. The described computational methods are a viable tool for studying the complex processes in multielectron atomic and molecular systems in strong laser fields.     

Photodissociation Spectroscopy and Dynamics of Weakly Bonded Metal Ion—Ethylene Complexes

The photodissociation spectroscopy of weakly bonded bimolecular complexes can give important insight into fundamental molecular interactions and dynamics. We have applied these techniques to a study of metal ion‐ethylene interactions in the Mg⁺(3s)‐C₂H₄ and Al⁺(3s²)‐C₂H₄ π‐bonded complexes. Experimental work is supported by ab‐initio electronic structure calculations. These experiments allow us to explore and compare the chemical binding, electronic structure, and nonadiabatic dissociation dynamics of these complexes.     

Isotope effect of hydrogen and lithium hydride molecules. Application of the dynamic extended molecular orbital method and energy component analysis

In order to explore the isotope effect including the nuclear–electronic coupling and nuclear quantum effects under the one-particle approximation, we apply the dynamic extended molecular orbital (DEMO) method and energy component analysis to the hydrogen and lithium hydride isotope molecules. Since the DEMO method determines both electronic and nuclear wave functions simultaneously by variationally optimizing all parameters embedded in the basis sets, the virial theorem is completely satisfied and guarantees the relation of the kinetic and potential energies. We confirm the isotope effect on internuclear distances, nuclear and electronic wave functions, dipole moment, the polarizability, and each energy component. In the case of isotopic species of the hydrogen molecule, the total energy decreases from the H₂ to the T₂ molecule due to the stabilization of the nuclear–electronic potential component, as well as the nuclear kinetic one. In the case of the lithium hydride molecule, the energy lowering by replacing ⁶Li with ⁷Li is calculated to be greater than that by replacing H with D. This is mainly caused by the small destabilization of electron–electron and nuclear–nuclear repulsion in ⁷LiH compared to ⁶LiH, while the change in the repulsive components from ⁶LiH to ⁶LiD increases. Received: 24 March 1999 / Accepted: 5 August 1999 / Published online: 15 December 1999     

Chemical Bonding as a New Avenue for Controlling Excited-State Properties and Excitation Energy-Transfer Processes in Zinc Phthalocyanine–Fullerene Dyads

Whether chemical bonding can regulate the excited-state and optoelectronic properties of donor–acceptor dyads has been largely elusive. In this work, we used electronic structure and nonadiabatic dynamics methods to explore the excited-state properties of covalently bonded zinc phthalocyanine (ZnPc)-fullerene (C₆₀) dyads with a 6–6 (or 5–6) bonding configuration in which ZnPc is bonded to two carbon atoms shared by the two hexagonal rings (or a pentagonal and a hexagonal ring) in C₆₀. In both cases, the locally excited (LE) states on ZnPc are spectroscopically bright. However, their different chemical bonding differentiates the electronic interactions between ZnPc and C₆₀. In the 5–6 bonding configuration, the LE states on ZnPc are much higher in energy than the LE states on C₆₀. Thus, the excitation energy transfer from ZnPc to C₆₀ is thermodynamically favorable. On the other hand, in the 6–6 bonding configuration, such a process is inhibited because the LE states on ZnPc are the lowest ones. More detailed mechanisms are elucidated from nonadiabatic dynamics simulations. In the 6–6 bonding configuration, no excitation energy transfer was observed. In contrast, in the 5–6 bonding configuration, several LE and charge-transfer (CT) excitons were shown to participate in the energy-transfer process. Further analysis reveals that the photoinduced energy transfer is mediated by a CT exciton, such that electron- and hole-transfer processes take place in a concerted but asynchronous manner in the excitation energy transfer. It is also found that high-level electronic structure methods including exciton effects are indispensable to accurately describe photoinduced energy- and electron-transfer processes. Furthermore, this work opens up new avenues for regulating the excited-state properties of molecular donor–acceptor dyads by means of chemical bonding.     

Ab initio surface hopping excited-state molecular dynamics approach on the basis of spin–orbit coupled states: An application to the A-band photodissociation of CH3I

Ab initio molecular dynamics approach has been extended to multi-state dynamics on the basis of the spin–orbit coupled electronic states that are obtained through diagonalization of the spin–orbit coupling matrix with the multi-state second-order multireference perturbation theory energies in diagonal elements and the spin–orbit coupling terms at the state-averaged complete active space self-consistent field level in off-diagonal elements. Nonadiabatic transitions over the spin–orbit coupled states were taken into account explicitly by a surface hopping scheme with utilizing the nonadiabatic coupling terms calculated by numerical differentiation of the spin–orbit coupled wavefunctions and analytical nonadiabatic coupling terms. The present method was applied to the A-band photodissociation of methyl iodide, CH₃I + hv → CH₃ + I (²P_3/2)/I* (²P_1/2), for which a pioneering theoretical work was reported by Amatatsu, Yabushita, and Morokuma. The present results reproduced well the experimental branching ratio and energy distributions in the dissociative products. © 2018 Wiley Periodicals, Inc.     

A multi-dimensional microcanonical Monte Carlo study of S0 → T1 intersystem crossing of isocyanic acid

A general formula for the multi-dimensional Monte Carlo microcanonical nonadiabatic rate constant expressed in configuration space is applied to calculate the rate of intersystem crossing (ISC) between the ground (S₀) and first excited triplet (T₁) states for isocyanic acid. One-, two- and three-dimensional potential energy surfaces are constructed by coupled-cluster single-double CCSD calculations, which are used for Monte Carlo sampling. The calculated S₀→T₁ ISC rate is in good agreement with experimental findings, which gives us a reason to believe that the multi-dimensional Monte Carlo microcanonical nonadiabatic rate theory is a very effective method for calculating nonadiabatic transition rate of a polyatomic molecule.     

Avoided crossings in metal (M)–gas (X) reactions (M = Hg, and X = SiH4, GeH4)

A study of nonadiabatic transitions through avoided crossings between two potential energy curves, associated to the approach of a mercury atom to an organic gas molecule (silane or germane) is presented. We study the Si–H and Ge–H bond breaking in the molecules SiH₄ and GeH₄, which are an important subject in the production of hydrogenated amorphous thin films. We here emphasize the importance of the excited states, the avoided crossings generated during the molecule–metal approach and the nonadiabatic transition probabilities. We have developed a model to extend the Landau–Zener theory utilizing the angle instead of the distance as the main parameter of the reaction, which is particularly adapted for tetrahedral molecules (as silane and germane). The activation process of these molecules requires several stages; first, we solve the Schrödinger equation (within the Born-Oppenheimer approximation) for the metal–molecule system during interaction. We always take into account all those states that can play a role in the reaction, even those that because of their energetic separation from the ground state are forgotten by other groups. The calculations begin at a LCAO-MO approximation and thenceforth variational and perturbative CI including of the order of a million determinants are carried out. Usually, some states of the metal repel the gas molecule and others attract it. This produces a series of avoided crossings among the curves, demanding that the nonadiabatic transition probabilities are obtained. This is the ultimate goal of the present study.     

Theoretical Shaping of Femtosecond Laser Pulses for Molecular Photodissociation with Control Techniques Based on Ehrenfest′s Dynamics and Time‐Dependent Density Functional Theory

The combination of nonadiabatic Ehrenfest‐path molecular dynamics (EMD) based on time‐dependent density functional theory (TDDFT) and quantum optimal control formalism (QOCT) was used to optimize the shape of ultra‐short laser pulses to achieve photodissociation of a hydrogen molecule and the trihydrogen cation H₃⁺. This work completes a previous one [A. Castro, ChemPhysChem, 2013 , 14, 1488–1495], in which the same objective was achieved by demonstrating the combination of QOCT and TDDFT for many‐electron systems on static nuclear potentials. The optimization model, therefore, did not include the nuclear movement and the obtained dissociation mechanism could only be sequential: fast laser‐assisted electronic excitation to nonbonding states (during which the nuclei are considered to be static), followed by field‐free dissociation. Here, in contrast, the optimization was performed with the QOCT constructed on top of the full dynamic model comprised of both electrons and nuclei, as described within EMD based on TDDFT. This is the first numerical demonstration of an optimal control formalism for a hybrid quantum–classical model, that is, a molecular dynamics method.     

FT-IR, FT-Raman, ab initio and DFT structural, vibrational frequency and HOMO-LUMO analysis of 1-naphthaleneacetic acid methyl ester

In this work, the FT-IR and FT-Raman spectra of 1-naphthaleneacetic acid methyl ester (abbreviated as 1-NAAME, C₁₀H₇CH₂CO₂CH₃) have been recorded in the region 3600–10 cm⁻¹. The optimum molecular geometry, normal mode wavenumbers, infrared and Raman intensities, Raman scattering activities, corresponding vibrational assignments, Mullikan atomic charges and other thermo-dynamical parameters were investigated with the help of HF and B3LYP (DFT) method using 6-31G(d,p), 6-311G(d,p) basis sets. Reliable vibrational assignments were made on the basis of total energy distribution (TED) calculated with scaled quantum mechanical (SQM) method. From the calculations, the molecules are predicted to exist predominantly as the C1 conformer. The correlation equations between heat capacity, entropy, enthalpy changes and temperatures were fitted by quadratic formulae. Lower value in the HOMO and LUMO energy gap explains the eventual charge transfer interactions taking place within the molecule. UV–VIS spectral analyses of 1NAAME have been researched by theoretical calculations. In order to understand electronic transitions of the compound, TD-DFT calculations on electronic absorption spectra in gas phase and solvent (DMSO and chloroform) were performed. The calculated frontier orbital energies, absorption wavelengths (λ), oscillator strengths (f) and excitation energies (E) for gas phase and solvent (DMSO and chloroform) are also illustrated.     

An experimental study of collisional transfer of electronic excitation,collisional dissociation of excited molecules and photodissociation in alkali vapour

A single-frequency laser is used to excite Na₂ molecules to the electronic B state. Besides the molecular fluorescence also atomic Na resonance radiation is observed. This is caused by: collisional transfer of electronic excitation from a Na₂(B) molecule to a Na atom, collisional dissociation of Na₂(B) molecules and photodissociation of Na₂ from very high vibration—rotation (v. J) levels of the ground state. We show how the contributions of these processes can be separated experimentally and characterized quantitativily over a wide range of temperatures, using a free-jet expansion. Illustrative results are given for one laser frequency (i.e. one molecular transition). Effective collision cross sections for excitation transfer and for collisional dissociation are given. The probability of photodissociation is compared to the probability of the discrete (B ← X) transition. The relaxation of the number of molecules in high (v. J) levels in a free jet is obtained.     

Numerical simulation of nonadiabatic electron excitation in the strong field regime. 2. Linear polyene cations

Time-dependent Hartree-Fock theory has been used to study the electronic optical response of a series of linear polyene cations (+1 and +2) in strong laser fields. The interaction of ethylene, butadiene, and hexatriene, with pulsed and CW fields corresponding to 8.75 x 10(13) W/cm(2) and 760 nm, have been calculated using the 6-31G(d,p) basis set. Nonadiabatic processes including nonlinear response of the dipole moment to the field and non-resonant energy deposition into excited states were more pronounced for the monocations in comparison with dications. For a given charge state and geometry, the nonadiabatic effects in the charge distribution and instantaneous dipole increased with the length of the polyene. For pulsed fields, the instantaneous dipole continued to oscillate after the field returned to zero and corresponded to a non-resonant electronic excitation involving primarily the lowest electronic transition. For a given molecule and fixed charge state, the degree of nonadiabatic coupling and excitation was greater for geometries with lower excitation energies.     

Molecular electric properties in electronic excited states: multipole moments and polarizabilities of H2O in the lowest1 B 1 and3 B 1 excited states

Summary The dipole and quadrupole moments and the dipole polarizability tensor components are calculated for the¹ B ₁ and³ B ₁ excited states of the water molecule by using the complete active space (CAS) SCF method and an extended basis set of atomic natural orbitals. The dipole moment in the lowest¹ B ₁ (0.640 a.u.) and³ B ₁ (0.416 a.u.) states is found to be antiparallel to that in the ground electronic state of H₂O. The shape of the quadrupole moment ellipsoid is significantly modified by the electronic excitation to both states investigated in this paper. All components of the excited state dipole polarizability tensor increase by about an order of magnitude compared to their values in the ground electronic state. The present results are used to discuss some aspects of intermolecular interactions involving molecules in their excited electronic states.     

Time-Dependent Density Functional Theory Study on the Electronic Excited State of Hydrogen-Bonded Clusters Formed by 2-Hydroxybenzonitrile (<Emphasis Type="Italic">o</Emphasis>-Cyanophenol) and Carbon Monoxide

Time-dependent density functional theory (TDDFT) method has been carried out to investigate excited-state hydrogen-bonding dynamics between 2-hydroxybenzonitrile (o-cyanophenol) and carbon monoxide. We have demonstrated that intermolecular hydrogen bond between 2-hydroxybenzonitrile (o-cyanophenol) and C=O group are significantly strengthened in the electronically excited state by theoretically monitoring the changes of the bond lengths of hydrogen bonds and hydrogen-bonding groups in different electronic states. In this study, we firstly analyze frontier molecular orbitals (MOs). Our results are consistent with the intermolecular hydrogen bond strengthening in the electronically excited state of Coumarin 102 in alcoholic solvents, which has been demonstrated for the first time by Zhao and Han. Moreover, the calculated electronic excitation energies of the hydrogen bonding C=O and O–H groups are markedly red-shifted upon photoexcitation, which illustrates the hydrogen bonds strengthen in the electronically excited state again. And the geometric structures in both ground state and the S₁ state of this hydrogen-bonded complex are calculated using the density functional theory (DFT) and TDDFT methods, respectively.     

Polarization of thallium atoms produced in molecular photodissociation: experiment and theory

A study has been made of the oriented ground state Tl(6² P _1/2) atoms produced in the photodissociation of TlBr molecules by circularly polarized 266-nm laser light. A significant degree of atomic orientation (15%) has been measured in the experiment which corresponds to the initial degree of orientation of 37%. A high value of depolarization cross section (210 Ų) for the oriented Tl atoms colliding with TlBr molecules has been also observed. The obtained experimental results have been treated theoretically. We present a general quantum mechanical theory of the orientation phenomenon in which all possible nonadiabatic interactions as well as molecular rotation are properly treated. The application of the theory to the case of TlBr photodissociation allowed to understand the obtained experimental results and to evaluate the probability of the earlier unknown radial nonadiabatic transition in the decaying molecule.