Insights into photodissociation dynamics of propionyl chloride from ab initio calculations and molecular dynamics simulations

The potential energy surfaces of isomerization, dissociation, and elimination reactions for CH3CH2COCl in the S0 and S1 states have been mapped with the different ab initio calculations. Mechanistic photodissociation of CH3CH2COCl at 266 nm has been characterized through the computed potential energy surfaces, the optimized surface crossing structure, intrinsic reaction coordinate, and ab initio molecular dynamics calculations. Photoexcitation at 266 nm leads to the CH3CH2COCl molecules in the S1 state. From this state, the C-Cl bond cleavage proceeds in a time scale of picosecond in the gas phase. The barrier to the C-Cl bond cleavage on the S1 surface is significantly increased by effects of the matrix and the internal conversion to the ground state prevails in the condensed phase. The HCl eliminations as a result of internal conversion to the ground state become the dominant channel upon photodissociation of CH3CH2COCl in the argon matrix at 10 K.     

Insights into mechanistic photodissociation of acetyl chloride by ab initio calculations and molecular dynamics simulations

The potential energy surfaces of dissociation and elimination reactions for CH(3)COCl in the ground (S0) and first excited singlet (S1) states have been mapped with the different ab inito calculations. Mechanistic photodissociation of CH(3)COCl has been characterized through the intrinsic reaction coordinate and ab initio molecular dynamics calculations. The alpha-C-C bond cleavage along the S1 pathway leads to the fragments of COCl((2)A' ') and CH(3) ((2)A') in an excited electronic state and a high barrier exists on the pathway. This channel is inaccessible in energy upon photoexcitation of the CH(3)COCl molecules at 236 nm. The S1 alpha-C-Cl bond cleavage yields the Cl((2)P) and CH(3)CO(X(2)A') fragments in the ground state and there is very small or no barrier on the pathway. The S1 alpha-C-Cl bond cleavage proceeds in a time scale of picosecond in the gas phase, followed by CH(3)CO decomposition to CH(3) and CO. The barrier to the C-Cl bond cleavage on the S1 surface is significantly increased by effects of the argon matrix. The S1 alpha-C-Cl bond cleavage in the argon matrix becomes inaccessible in energy upon photoexcitation of CH(3)COCl at 266 nm. In this case, the excited CH(3)COCl(S1) molecules cannot undergo the C-Cl bond cleavage in a short period. The internal conversion from S1 to S0 becomes the dominant process for the CH(3)COCl(S1) molecules in the condensed phase. As a result, the direct HCl elimination in the ground state becomes the exclusive channel upon 266 nm photodissociation of CH(3)COCl in the argon matrix at 11 K.     

Direct ab initio molecular dynamics study of the two photodissociation channels of formic acid

A total of ∼1200 trajectories have been integrated for the two photodissociation channels of formic acid, HCOOH → H₂O + CO (1) and HCOOH → CO₂ + H₂ (2), which occur with 248 and 193 nm photons, using the direct ab initio molecular dynamics method at the RMP2(full)/cc-pVDZ level of theory. It was found that the percentage of the energy distributed to a relative translational mode in reaction (2) is much larger than that in reaction (1). This is mainly due to the difference in the geometry of transition state (TS); the H₂O geometry in the TS of reaction (1) was predicted to significantly deviate from the equilibrium one, whereas the CO₂ and H₂ geometries in the TS of reaction (2) were found to be more similar to their equilibrium ones. It was also found that the product diatomic molecules, CO and H₂, are both vibrationally and rotationally excited. The calculated relative population of the vibrationally excited CO for the 248 nm photodissociation was consistent with experiment.     

Multireference theoretical investigation on selectivity of the bond fissions in photodissociation of acetyl cyanide

The selectivity of the C-CH(3) and C-CN bond fissions upon excitation of acetyl cyanide at 193 nm has been investigated at the theoretical level of multistate complete active space self-consistent field second order perturbation. The calculated results indicated that the initially excited S(3) state relaxes to S(2) via ultrafast internal conversion. The S(2) state could dissociate via two pathways. One, adiabatically dissociates to CH(3)CO(X)+CN(A). The other one internally converts to S(1) before S(1) intersystem crossing to T(1). The T(1) state subsequently dissociates to two groups of products: CH(3)(X)+OCCN(X) and CH(3)CO(X)+CN(X). The experimentally observed preference branching of CN elimination over CH(3) one and bond selectivity are the results of the competition between the adiabatic and nonadiabatic dynamics of the S(2) state.     

Nonadiabatic effects in C-Br bond scission in the photodissociation of bromoacetyl chloride

Bromoacetyl chloride photodissociation has been interpreted as a paradigmatic example of a process in which nonadiabatic effects play a major role. In molecular beam experiments by Butler and co-workers [J. Chem. Phys. 95, 3848 (1991); J. Chem. Phys. 97, 355 (1992)], BrCH2C(O)Cl was prepared in its ground electronic state (S0) and excited with a laser at 248 nm to its first excited singlet state (S1). The two main ensuing photoreactions are the ruptures of the C-Cl bond and of the C-Br bond. A nonadiabatic model was proposed in which the C-Br scission is strongly suppressed due to nonadiabatic recrossing at the barrier formed by the avoided crossing between the S1 and S2 states. Recent reduced-dimensional dynamical studies lend support to this model. However, another interpretation that has been given for the experimental results is that the reduced probability of C-Br scission is a consequence of incomplete intramolecular energy redistribution. To provide further insight into this problem, we have studied the energetically lowest six singlet electronic states of bromoacetyl chloride by using an ab initio multiconfigurational perturbative electronic structure method. Stationary points (minima and saddle points) and minimum energy paths have been characterized on the S0 and S1 potential energy surfaces. The fourfold way diabatization method has been applied to transform five adiabatic excited electronic states to a diabatic representation. The diabatic potential energy matrix of the first five excited singlet states has been constructed along several cuts of the potential energy hypersurfaces. The thermochemistry of the photodissociation reactions and a comparison with experimental translational energy distributions strongly suggest that nonadiabatic effects dominate the C-Br scission, but that the reaction proceeds along the energetically allowed diabatic pathway to excited-state products instead of being nonadiabatically suppressed. This conclusion is also supported by the low values of the diabatic couplings on the C-Br scission reaction path. The methodology established in the present study will be used for the construction of global potential energy surfaces suitable for multidimensional dynamics simulations to test these preliminary interpretations.     

Nonadiabatic interactions in wave packet dynamics of the bromoacetyl chloride photodissociation

The competitive photodissociation of bromoacetyl chloride BrCH2COCl in the first 1A" state (S1) by 248 nm photons is investigated by nonadiabatic wave packet simulations. We show that the preferential breaking of the stronger C-Cl bond (alpha to the excited carbonyl) over the weaker C-Br bond (beta) could be explained by a diabatic trapping or nonadiabatic recrossing as previously proposed. Our energy resolved flux analysis agrees fairly well with the experimental branching ratio (C-Cl:C-Br=1.0:0.4). Even if this does not prove the mechanism, this at least prevents to discard it. A reduced dimensionality approach based on constrained Hamiltonian is used. The nonadiabatic dissociation is studied in the two C-O/C-X (X=Br, Cl) subspaces to emphasize the role of the C-O vibration upon [nO-->piCO*] excitation. The internal torsion and wagging dihedral angles are frozen at their Franck-Condon value, according to preliminary dynamical tests. The other inactive coordinates are optimized at the trans and Cs constrained geometry in the first excited state. Corresponding 2D cuts in the potential energy surfaces have been computed at the CASSCF level. The nonadiabatic kinetic couplings are highly peaked along an avoided crossing seam in both cases. A two-state diabatic model with a constant potential coupling is proposed in the two C-O/C-X subspaces. The inclusion of the C-O stretching in the active coordinates improves the value of the branching ratio over our previous 1D computation.     

Controlling the shape and flexibility of arylamides: a combined ab initio, ab initio molecular dynamics, and classical molecular dynamics study

Using quantum chemistry plus ab initio molecular dynamics and classical molecular dynamics methods, we address the relationship between molecular conformation and the biomedical function of arylamide polymers. Specifically, we have developed new torsional parameters for a class of these polymers and applied them in a study of the interaction between a representative arylamide and one of its biomedical targets, the anticoagulant drug heparin. Our main finding is that the torsional barrier of a C(aromatic)-C(carbonyl) bond increases significantly upon addition of an o-OCH2CH2NH3+ substituent on the benzene ring. Our molecular dynamics studies that are based on the original general AMBER force field (GAFF) and GAFF modified to include our newly developed torsional parameters show that the binding mechanism between the arylamide and heparin is very sensitive to the choice of torsional potentials. Ab initio molecular dynamics simulation of the arylamide independently confirms the degree of flexibility we obtain by classical molecular dynamics when newly developed torsional potentials are used.     

Bifurcated hydrogen bond in lithium nitrate trihydrate probed by ab initio molecular dynamics

The hydrogen-bond dynamics of lithium nitrate trihydrate has been studied by a combined approach based on ab initio molecular dynamics simulations and wavelet analysis. The simultaneous bifurcated interaction between one hydrogen atom of water molecules and two oxygen atoms of nitrate ions is the pivotal feature of the crystal structure: this bifurcated interaction has deep effects on the O-H stretching region of the vibrational spectrum. The structural, dynamic, spectroscopic, and electronic properties of the bifurcated hydrogen bond have been investigated computationally, elucidating at the molecular level the differences with weak and strong hydrogen bonds present in the crystal. These studies corroborate the very recent IR experiments performed on the lithium nitrate trihydrate crystal, offering new perspectives to interpreting the vibrational spectra. In fact, this approach allows obtaining two-dimensional plots, which summarize the essential features of both the hydrogen-bond network and IR spectra, resulting in a peculiar "signature" of the bifurcated interaction.     

Ab initio quantum dynamical study of the multi-state nonadiabatic photodissociation of pyrrole

There has been a substantial amount of theoretical investigations on the photodynamics of pyrrole, often relying on surface hopping techniques or, if fully quantal, confining the study to the lowest two or three singlet states. In this study we extend ab initio based quantum dynamical investigations to cover simultaneously the lowest five singlet states, two π-σ? and two π-π? excited states. The underlying potential energy surfaces are obtained from large-scale MRCI ab initio computations. These are used to extract linear and quadratic vibronic coupling constants employing the corresponding coupling models. For the N-H stretching mode Q(24) an anharmonic treatment is necessary and also adopted. The results reveal a sub-picosecond internal conversion from the S(4) (π-π?) state, corresponding to the strongly dipole-allowed transition, to the S(1) and S(2) (π-σ?) states and, hence, to the ground state of pyrrole. The significance of the various vibrational modes and coupling terms is assessed. Results are also presented for the dissociation probabilities on the three lowest electronic states.     

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first-principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.     

Interconversion behavior of the C-H bond in the CH4+ radical cation: ab initio molecular dynamics study

Thermal motion of CH4+ is investigated by performing an ab initio molecular dynamics method with the second-order M?ller-Plesset (MP2)/6-311G** force field. In the trajectories obtained at 400 K, we have observed rapid interconversion behavior of the geometrical parameters of CH4+ with the frequency of 0.6/ps, where the C-H pair forming the small angle around 55 degrees is switched to another pair on subpicosecond time scale. The switching patterns are found to be classified into the following two types. Type 1: one C-H of the small angled C-H pair is switched to one C-H of the other C-H pair. Type 2: the small angled C-H pair is switched to the other C-H pair, which has been newly observed in the present ab initio MD calculation. The four C-H bonds of CH4+ are characterized by the long and short C-H bonds in a time region of the trajectories, and also for the time-evolution of C-H bonds such interconversion behavior is observed. The switching patterns of the geometrical parameters are compared with those in the interconversion scheme between six equivalent C2v symmetry structures of CH4+ [Paddon-Row, M. N. et al., J Am Chem Soc 1985, 107, 7696]. We have also investigated the electronic energy fluctuation due to thermal motion of CH4+. The standard deviation of total electronic energy at 400 K is evaluated to be 1.2 kcal/mol.     

Combined ab initio and classical molecular dynamics simulations of alkyl-lithium aggregates in ethereal solutions

Molecular dynamics simulations of organolithium aggregates in solution are reported for the first time. We use a combined quantum/classical force field (the so-called QM/MM approach) and study ethyl-lithium aggregates in dimethyl ether (DME) solvent. The solutes are described at the Density Functional Theory level while solvent molecules are described using molecular mechanics. NVT Molecular Dynamics simulations at 200 K are carried out in the Born–Oppenheimer approximation. After equilibration, the production phase was run for 80 ps (monomer), 40 ps (dimer) and 26 ps (tetramer). The analysis of the results focuses on Li coordination as a function of aggregate size and we show that the total Li coordination number is always 4. No decoordination has been observed along the simulations. Fluctuations of the structures are predicted to be large in some cases and possible implications on reactivity are discussed.     

Insights into photodissociation dynamics of benzamide and formanilide from ab initio calculations

In the present study, the five lowest electronic states that control the UV photodissociation of formanilide and benzamide have been characterized using the complete active space self-consistent field theory. The mechanisms for the initial relaxation and subsequent dissociation processes have been determined on the basis of the calculated potential energy surfaces and their intersections. It was found that the S(1)/T(1)/T(2) three-surface intersection plays an important role in the photodissociation processes of benzamide. However, the dissociation behavior of formanilide and benzamide was found to be quite different from that for aliphatic amides. The present study provides several insights into the photodissociation dynamics of formanilide and benzamide.     

An ab initio molecular dynamics and density functional theory study of the formation of phosphate chains from metathiophosphates

The reactions that occur between metathiophosphate (MTP) molecules are identified and examined through ab initio molecular dynamics simulations and static quantum chemical calculations at the density functional level of theory. The simulations show that certain types of MTPs can react to yield phosphate chains, while others only dimerize. These differences are rationalized in terms of reaction energies and the electronic structures of these molecules. In the reaction leading to the formation of phosphate chains, the reactive center, a tri-coordinate phosphorus atom, is continually regenerated. A polymerization mechanism linking MTPs to phosphate chains is developed on the basis of these results. This information sheds light on the underlying processes that may be responsible for the formation of phosphates under high-temperature conditions and may prove useful in the development of protocols for the rational synthesis of complex phosphate structures.     

An ab initio molecular dynamics study on hydrogen bonds between water molecules

The quantitative estimation of the total interaction energy of a molecular system containing hydrogen bonds (H bonds) depends largely on how to identify H bonding. The conventional geometric criteria of H bonding are simple and convenient in application, but a certain amount of non-H bonding cases are also identified as H bonding. In order to investigate the wrong identification, we carry out a systematic calculation on the interaction energy of two water molecules at various orientation angles and distances using ab initio molecular dynamics method with the dispersion correction for the Becke-Lee-Yang-Parr (BLYP) functionals. It is shown that, at many orientation angles and distances, the interaction energies of the two water molecules exceed the energy criterion of the H bond, but they are still identified as H-bonded by the conventional "distance-angle" criteria. It is found that in these non-H bonding cases the wrong identification is mainly caused by short-range interaction between the two neighbouring water molecules. We thus propose that, in addition to the conventional distance and angle criteria of H bonding, the distance d(H···H) between the two neighbouring hydrogen atoms of the two water molecules should also be taken as a criterion, and the distance r(O···H) between the hydrogen atom of the H-bond donor molecule and the oxygen atom of the acceptor molecule should be restricted by a lower limit. When d(H···H) and r(O···H) are small (e.g., d(H···H) < 2.0 ? and r(O···H) < 1.62 ?), the repulsion between the two neighbouring atoms increases the total energy of the two water molecules dramatically and apparently weakens the binding of the water dimer. A statistical analysis and comparison of the numbers of the H bonds identified by using different criteria have been conducted on a Car-Parrinello ab initio molecular dynamics simulation with dispersion correction for a system of 64 water molecules at near-ambient temperature. They show that the majority of the H-bonds counted by using the conventional criteria combined with the d(H···H) criterion and the restriction of r(O···H) match what is identified by the binding energy criteria (e.g., E ≤ -10 kJ/mol), while some of them still have a binding energy that exceeds the energy criterion, indicating that the complicated quantum effects in H bonding can only be described by the three geometric parameters to a certain extent.     

Vibrational dynamics of the hydrogen bond in H(2)S-HF: Fourier-transform-infrared spectra and ab initio theory

The rotationally resolved infrared spectrum of the hydrogen bonded complex H(2)S-HF and of its isotopomer D(2)S-DF in the HF/DF stretching range have been observed in a supersonic jet Fourier-transform infrared (FTIR) experiment and indicate a predissociation lifetime of 130 ps for H(2)S-HF. Complementary spectra taken at a temperature of 190 K in a cell without resolved rotational structure indicate the presence of strong anharmonic couplings between low frequency intermolecular modes and the HF donor stretch mode previously observed in other complexes with heavier acceptor molecules without rotational fine structure. The anharmonic analysis of the hot band progressions and of the rotational data confirm the coupling mechanism. The coupling constants and the absolute frequency of the hydrogen bonded stretch mode are in excellent agreement with theoretical predictions based on adiabatic variational calculations on potential surfaces computed at MP2 and CCSD(T) level. Complementary calculations with a perturbational approach further confirm the coupling model.     

Energy-flow dynamics in the molecular channel of propanal photodissociation, C2H5CHO --> C2H6 + CO: direct ab initio molecular dynamics study

Direct ab initio molecular dynamics calculations have been carried out for the molecular channel of the photodissociation of propanal, C2H5CHO --> C2H6 + CO, at the RMP2(full)/cc-pVDZ level of ab initio molecular orbital theory. The initial conditions were generated using the microcanonical sampling to put the excess energy randomly into all vibrational modes of the TS. Starting from the TS, a total of approximately 700 trajectories were numerically integrated for 100 fs. The obtained final energy distributions for the C2H6 and CO fragments and their relative translational motion were found to be quite similar to those obtained for the acetaldehyde reaction, CH3CHO --> CH4 + CO, in our previous study (Chem. Phys. Lett. 2006, 421, 549) despite the fact that the number of degree of freedom for C2H6 is larger than that for CH4. The coupling between the intrinsic reaction coordinate and one of the generalized normal modes orthogonal to it was predicted substantially strong around s = 1.4 amu(1/2) bohr, and it is expected that the energy flow out of C2H6 proceeds through this coupling. However, the obtained energy distributions strongly suggest that the coupling among the modes in C2H6 is quite small and the intramolecular energy redistribution does not occur efficiently in this molecule.