NeI2 Photofragmentation Dynamics Through Quasiclassical and Semiclassical Studies.

The triatomic system NeI2 is studied under the consideration that the diatom is found in an excited electronic state (B). The vibrational levels (v=13, …, 23) are considered within two well-known theoretical procedures: quasi-classical trajectories (QCT), where the classical equations of motion for nuclei are solved on a single potential energy surface (PES), and the trajectory surface hopping (TSH) method, where the same are solved in a bunch of crossed vibrational PES (diabatic representation). The trajectory surface hopping fewest switches (TSHFS) is implemented to minimize the number of hoppings, thus allowing the calculations of hopping probability between the different PES's, and the kinetic mechanism to track the dissociation path. From these calculations, several observables such as, the lifetimes, vibrational and rotational energies (I2 ), dissociation channels, are obtained. Our results are compared with previous experimental and theoretical work.

Two-Photon Excited Fluorescence Dynamics in NADH in Water-Methanol Solutions: The Role of Conformation States.

The dynamics of polarized fluorescence in reduced nicotinamide adenine dinucleotide (NADH) at 460 nm under two-photon excitation at 720 nm by femtosecond laser pulses in water-methanol solutions has been studied experimentally and theoretically as a function of methanol concentration. A number of fluorescence parameters have been determined from experiment by means of the global fit procedure and then compared with the results reported by other authors. A comprehensive analysis of experimental errors was made. Ab initio calculations of the structure of NADH in water and methanol and of ß-nicotinamide mononucleotide (NMNH) in vacuum have been carried out for clarifying the role of decay time heterogeneity. The main results obtained are as follows. An explanation of the heterogeneity in the measured fluorescence decay times in NADH has been suggested based on the influence of the internal molecular electric field in the nicotinamide ring on nonradiative decay rates. We suggest that different charge distributions in the cis and trans configurations result in different internal electrostatic field distributions that lead to the decay time heterogeneity. A slight but noticeable rise of the fluorescence decay times τ1 and τ2 with methanol concentration was observed and treated as a minor effect of a nonradiative relaxation slowing due to the decrease in solution polarity. Relative concentrations of the folded and unfolded NADH conformations in solutions have been determined using a new method of analysis of the rotational diffusion time τr as a function of methanol concentration on the basis of the Stokes-Einstein-Debye equation. The analysis of the fluorescence anisotropy parameters obtained under linearly and circularly polarized excitation and the parameter Ω has been carried out and resulted in the determination of the two-photon excitation tensor components and suggested the existence of two excitation channels with comparable intensities. These were the longitudinal excitation channel dominated by the diagonal tensor component Szz and the mixed excitation channel dominated by the off-diagonal tensor components |Sxz2 + Syz2|1/2.

Mechanistic Insights into the Formation of Lithium Fluoride Nanotubes.

Born-Oppenheimer molecular dynamics (BOMD) and periodic density functional theory (DFT) calculations have been applied for describing the mechanism of formation of lithium fluoride (LiF) nanotubes with cubic, hexagonal, octagonal, decagonal, dodecagonal, and tetradecagonal cross-sections. It has been shown that high energy structures, such as nanowires, nanorings, nanosheets, and nanopolyhedra are transient species for the formation of stable nanotubes. Unprecedented (LiF)n clusters (n≤12) were also identified, some of them lying less than 10 kcal mol-[1] above their respective global minima. Such findings indicate that stochastic synthetic techniques, such as laser ablation and chemical vapor deposition, should be combined with a template-driven procedure in order to generate the nanotubes with adequate efficiency. Apart from the stepwise growth of LiF units, the formation of nanotubes was also studied by rolling up a planar square sheet monolayer, which could be hypothetically produced from the exfoliation of the FCC crystal structure. It was shown that both pathways could lead to the formation of alkali halide nanotubes, a still unprecedented set of one-dimensional materials.

Time-Dependent Quantum Wave Packet Study of the Si + OH → SiO + H Reaction: Cross Sections and Rate Constants.

The dynamics of the Si(3P) + OH(X2Π) → SiO(X1Σ+) + H(2S) reaction is investigated by means of the time-dependent wave packet (TDWP) approach using an ab initio potential energy surface recently developed by Dayou et al. ( J. Chem. Phys. 2013 , 139 , 204305 ) for the ground X2A' electronic state. Total reaction probabilities have been calculated for the first 15 rotational states j = 0-14 of OH(v=0,j) at a total angular momentum J = 0 up to a collision energy of 1 eV. Integral cross sections and state-selected rate constants for the temperature range 10-500 K were obtained within the J-shifting approximation. The reaction probabilities display highly oscillatory structures indicating the contribution of long-lived quasibound states supported by the deep SiOH/HSiO wells. The cross sections behave with collision energies as expected for a barrierless reaction and are slightly sensitive to the initial rotational excitation of OH. The thermal rate constants show a marked temperature dependence below 200 K with a maximum value around 15 K. The TDWP results globally agree with the results of earlier quasi-classical trajectory (QCT) calculations carried out by Rivero-Santamaria et al. ( Chem. Phys. Lett. 2014 , 610-611 , 335 - 340 ) with the same potential energy surface. In particular, the thermal rate constants display a similar temperature dependence, with TDWP values smaller than the QCT ones over the whole temperature range.

A theoretical simulation of the resonant Raman spectroscopy of the H2O⋯Cl2 and H2O⋯Br2 halogen-bonded complexes.

The resonant Raman spectra of the H2O⋯Cl2 and H2O⋯Br2 halogen-bonded complexes have been studied in the framework of a 2-dimensional model previously used in the simulation of their UV-visible absorption spectra using time-dependent techniques. In addition to the vibrational progression along the dihalogen mode, a progression is observed along the intermolecular mode and its combination with the intramolecular one. The relative intensity of the inter to intramolecular vibrational progressions is about 15% for H2O⋯Cl2 and 33% for H2O⋯Br2. These results make resonant Raman spectra a potential tool for detecting the presence of halogen bonded complexes in condensed phase media such as clathrates and ice.

Quasi-classical trajectories study of Ne2Br2(B) vibrational predissociation: kinetics and product distributions.

The vibrational predissociation of the Ne(2)Br(2)(B) van der Waals complex has been investigated using the quasi-classical trajectory method (QCT), in the range of vibrational levels v(') = 16-23. Extensive comparison is made with the most recent experimental observations [Pio et al., J. Chem. Phys. 133, 014305 (2010)], molecular dynamics with quantum transitions simulations [Miguel et al., Faraday Discuss. 118, 257 (2001)], and preliminary results from 24-dimensional Cartesian coupled coherent state (CCCS) calculations. A sequential mechanism is found to accurately describe the theoretical dynamical evolution of intermediate and final product populations, and both QCT and CCCS provide very good estimates for the dissociation lifetimes. The capabilities of QCT in the description of the fragmentation kinetics are analyzed in detail by using reduced-dimensionality models of the complexes and concepts from phase-space transport theory. The problem of fast decoupling of the different coherent states in CCCS simulations, resulting from the high dimensionality of phase space, is tackled using a re-expansion scheme. QCT ro-vibrational product state distributions are reported. Due to the weakness of the van der Waals couplings and the low density of vibrational states, QCT predicts a larger than observed propensity for Δv' = -1 and -2 channels for the respective dissociation of the first and second Ne atoms.

A new ab initio potential energy surface for studying vibrational relaxation in NO(v) + NO collisions.

A new ab initio potential energy surface for the ground state of the NO-NO system has been calculated within a reduced dimensionality model. We find an unusually large vibrational dependence of the interaction potential which explains previous spectroscopic observations. The potential can be used to model vibrational energy transfer, and here we perform quantum scattering calculations of the vibrational relaxation of NO(v). We show that the vibrational relaxation for v = 1 is 4 orders of magnitude larger than that for the related O(2)(v) + O(2) system without having to invoke nonadiabatic mechanisms as had been suggested in the past. For highly vibrationally excited states, we predict a strong dependence of the rates on the vibrational quantum number as has been observed experimentally, although there remain important quantitative differences. The importance of a chemically bound isomer on the relaxation mechanism is analyzed, and we conclude it does not play a role for the values of v considered in the experiment. Finally, the intriguing negative temperature dependence of the vibrational relaxation rate constants observed in experiments was studied using an statistical model to include the presence of many asymptotically degenerate spin-orbit states.

Large shift and small broadening of Br2 valence band upon dimer formation with H2O: an ab initio study.

Valence electronic excitation spectra are calculated for the H(2)O···Br(2) complex using highly correlated ab initio potentials for both the ground and the valence electronic excited states and a 2-D approximation for vibrational motion. Due to the strong interaction between the O-Br and the Br-Br stretching motions, inclusion of these vibrations is the minimum necessary for the spectrum calculation. A basis set calculation is performed to determine the vibrational wave functions for the ground electronic state and a wave packet simulation is conducted for the nuclear dynamics on the excited state surfaces. The effects of both the spin-orbit interaction and temperature on the spectra are explored. The interaction of Br(2) with a single water molecule induces nearly as large a shift in the spectrum as is observed for an aqueous solution. In contrast, complex formation has a remarkably small effect on the T = 0 K width of the valence bands due to the fast dissociation of the dihalogen bond upon excitation. We therefore conclude that the widths of the spectra in aqueous solution are mostly due to inhomogeneous broadening.

Cartesian coupled coherent states simulations: Ne(n)Br2 dissociation as a test case.

In this article, we describe coupled coherent states (CCS) simulations of vibrational predissociation of weakly bounded complexes. The CCS method is implemented in the Cartesian frame in a manner that is similar to classical molecular dynamics. The calculated lifetimes of the vibrationally excited Ne-Br(2)(ν) complexes agree with experiment and previous calculations. Although the CCS method is, in principle, a fully quantum approach, in practice it typically becomes a semiclassical technique at long times. This is especially true following dissociation events. Consequently, it is very difficult to converge the quantum calculations of the final Br(2) vibrational distributions after predissociation and of the autocorrelation functions. However, the main advantage of the method is that it can be applied with relative ease to determine the lifetimes of larger complexes and, in order to demonstrate this, preliminary results for tetra- and penta-atomic clusters are reported.

Argon solid response upon Rydberg photoexcitation of the NO chromosphore: case of using ab initio potential energy surfaces and comparison to similar studied systems.

Molecular dynamics simulations of NO-doped Ar solid upon Rydberg photoexcitation of the impurity have been carried out taking into account angular dependent potential energy surfaces (PESs) in the ground and excited states. To go beyond isotropic potentials simulations, the effects of anisotropy of potentials on the structure, dynamics, and energetics are investigated by taking into account two cases, namely, the whole PESs and the isotropic parts. Results have been compared to those obtained in previous works for similar NO-doped rare gas systems. Radial distribution functions (RDF) for the ground and excited state indicate that for both cases the shell structure of the lattice is kept ordered and is characterized by well-defined bands. No influence of the anisotropy of potentials has been detected in the RDFs since the anisotropy is rather manifested at short distances. The well part, which has been proven to be unimportant for the dynamics in previous works, arises here to be important for the right simulation of the spectrum. In general, our results show a reasonable agreement with respect to the experimental values for the dynamics and energetics when ab initio potentials are used, although better results can be obtained if higher level ab initio PESs are used.

On the local relaxation of solid neon upon Rydberg excitation of a NO impurity: the role of the NO(A)-Ne interaction potential and zero-point quantum delocalization.

The local relaxation of solid neon subsequent to the impulsive excitation of the NO chromophore to its A(3s sigma) Ryberg state is investigated using molecular dynamics simulations. This study makes use of empirical NO(X,A)-Ne isotropic pair potentials as well as a recently developed ab initio triatomic potential energy surface for the excited state. The role of these interaction potentials is analyzed, including many-body effects. In particular, empirical potentials, designed to reproduce correctly both the NO X-A steady-state absorption and emission bands, are shown to lead to a good description of the subpicosecond relaxation dynamics. The 600 fs expansion of the electronic bubble fairly agrees with experimental data. This relatively long time scale with respect to solid Argon, which was previously attributed to the range of the NO(A)-Ne interaction, is presumably related to the quantum nature of the medium. The time-resolved local relaxation of the Ne solid is understandably intermediate between that of classical solids (e.g., Ar) and that of quantum solids (e.g., H(2)).

An ab initio study of the Ar-NO(A 2Sigma+) intermolecular potential.

More complete molecular dynamics simulations of NO doped Ar solid upon photoexcitation of the impurity should include effects of angular variations of Ar-NO intermolecular potential. This is the main reason for presenting in this work an ab initio study of the Ar-NO(A (2)Sigma(+)) intermolecular potential. Ab initio calculations were carried out at the level of CASSCF-MRCI, with the aug-cc-pVTZ basis sets. In order to evaluate the influence of the quadruple excitations on the topology of the potential energy surface (PES), two cases were considered, that is, with and without taking into account Davidson's correction for quadruple excitations during the calculations. An analytical representation of the PES has been obtained as a function of the Jacobi coordinates of the system. In general, the PES is repulsive, except for linear directions, where two shallow wells appear. When quadruple excitations are considered, wells are located at 4.2 A (alpha=0 degrees) and 6.08 A (alpha=180 degrees) with energies of -20 and -15 cm(-1), respectively; and when are not considered, wells are located at 6.1 A (alpha=0 degrees) and 6.8 A (alpha=180 degrees) with energies of -15 and -10 cm(-1), respectively. For distances beyond 7 A, it is observed a very low energy decay and a rapid tendency to isotropic interactions.

An ab initio calculation of the valence excitation spectrum of H2O...Cl2: comparison to condensed phase spectra.

Valence electronic excitation spectra are calculated for the H(2)O...Cl(2) dimer using state-of-the art ab initio potentials for both the ground and the valence excited states, a basis set calculation of the ground state nuclear wave function, and a wave packet analysis to simulate the dynamics on the excited state surface. The peak of the H(2)O...Cl(2) dimer spectrum is blue-shifted by 1250 cm(-1) from that of the free Cl(2) molecule. This is less than the value previously estimated from vertical excitation energies but still significantly more than the blue shift in aqueous solution and clathrate-hydrate solid. Seventy percent of the blue shift is attributed to ground state stabilization, the rest to excited state repulsion. Spin-orbit effects are found to be small for this dimer. Homogeneous broadening is found to be slightly smaller for the dimer than for the free Cl(2). The reflection principle and spectator model approximations were tested and found to be quite satisfactory. This is promising for an eventual simulation of the condensed phase spectra.

Ab initio ground and excited state potential energy surfaces for NO-Kr complex and dynamics of Kr solids with NO impurity.

The intermolecular potentials for the NO(X 2Pi)-Kr and NO(A 2Sigma+)-Kr systems have been calculated using highly accurate ab initio calculations. The spin-restricted coupled cluster method for the ground 1 2A' state [NO(X 2Pi)-Kr] and the multireference singles and doubles configuration interaction method for the excited 2 2A' state [NO(A 2Sigma+)-Kr], respectively, were used. The potential energy surfaces (PESs) show two linear wells and one that is almost in the perpendicular position. An analytical representation of the PESs has been constructed for the triatomic systems and used to carry out molecular dynamics (MD) simulations of the NO-doped krypton matrix response after excitation of NO. MD results are shown comparatively for three sets of potentials: (1) anisotropic ab initio potentials [NO molecule direction fixed during the dynamics and considered as a point (its center of mass)], (2) isotropic ab initio potentials (isotropic part in a Legendre polynomial expansion of the PESs), and (3) fitted Kr-NO potentials to the spectroscopic data. An important finding of this work is that the anisotropic and isotropic ab initio potentials calculated for the Kr-NO triatomic system are not suitable for describing the dynamics of structural relaxation upon Rydberg excitation of a NO impurity in the crystal. However, the isotropic ab initio potential in the ground state almost overlaps the published experimental potential, being almost independent of the angle asymmetry. This fact is also manifested in the radial distribution function around NO. However, in the case of the excited state the isotropic ab initio potential differs from the fitted potentials, which indicates that the Kr-NO interaction in the matrix is quite different because of the presence of the surrounding Kr atoms acting on the NO molecule. MD simulations for isotropic potentials reasonably reproduce the experimental observables for the femtosecond response and the bubble size but do not match spectroscopic results. A general overall view of the results suggests that, when the Kr-NO interaction takes place inside the matrix, potentials are rather symmetric and less repulsive than those for the triatomic system.