Volume 121, Issue 23, 15 December 2004
 COMMUNICATIONS


Predissociation spectroscopy of the argonsolvated “zundel” cation in the 1000–1900 cm^{−1} region
View Description Hide DescriptionPredissociationspectra of the cluster ions are reported in the 1000–1900 cm^{−1} region. The weakly bound argon atoms enable investigation of the complex in a linear action mode, and the resulting spectra are much simpler than those reported previously in this region [Asmis et al., Science 299, 1375 (2003) and Fridgen et al., J. Phys. Chem. A 108, 9008 (2004)], which were obtained using infrared multiphoton dissociation of the bare complex. The observed spectrum consists of two relatively narrow bands at 1080 and 1770 cm^{−1} that are likely due to excitation of the shared proton and intramolecular bending vibrations of the two water molecules, respectively. The narrow linewidths and relatively small (60 cm^{−1}) perturbation introduced by the addition of a second argon atom indicate that the basic “zundel” character of the ion survives upon complexation.

Nearinfrared electronic spectrum of
View Description Hide DescriptionFour vibronic bands of the electronic transition of have been observed in the near infrared from using a highly sensitive spectrometer based on a Ti:sapphire laser and a positive column plasma. The two states are derived from the lowest state of the linear conformation as a result of a strong Renner–Teller effect. As expected, the rovibronic structure of the spectrum is very similar to that of the isoelectronic radical reported by Herzberg and Johns. Results of a preliminary analysis are given in which the rotational structures of the two states are treated independently, the ground state as a nearprolate asymmetric rotor and the excited state as a linear molecule.

Realtime detection of doorway states in the intramolecular vibrational energy redistribution of the OH/OD stretch vibration of phenol
View Description Hide DescriptionA picosecond timeresolved IRUV pump–probe spectroscopic study was carried out for the intramolecular vibrational energy redistribution of the OH/OD stretching vibration of isolated phenol and its isotopomers in supersonic beams. The time evolution due to IVR showed a significant isotope effect; the OH stretch vibration showed a single exponential decay and its lifetime is greatly lengthened upon the deuterium substitution of the CH group. The OD stretch vibration exhibited prominent quantum beats. Especially, in phenol the electronic transitions from the doorway states were clearly observed. They exhibited an outofphase quantum beat with respect to that of the OD stretch level and disappeared due to further IVR to the dense bath states. The transient spectra as well as the time evolution clearly evidenced the tiermodel of the description of intramolecular vibrational energy redistribution.
 Top

 ARTICLES

 Theoretical Methods and Algorithms

Using preconditioned adaptive step size RungeKutta methods for solving the timedependent Schrödinger equation
View Description Hide DescriptionIf the Hamiltonian is time dependent it is common to solve the timedependent Schrödinger equation by dividing the propagation interval into slices and using an (e.g., split operator, Chebyshev, Lanczos) approximate matrix exponential within each slice. We show that a preconditioned adaptive step size RungeKutta method can be much more efficient. For a chirped laser pulse designed to favor the dissociation of HF the preconditioned adaptive step size RungeKutta method is about an order of magnitude more efficient than the time sliced method.

Curvysteps approach to constraintfree extendedLagrangian ab initio molecular dynamics, using atomcentered basis functions: Convergence toward Born–Oppenheimer trajectories
View Description Hide DescriptionA dynamical extension of the “curvysteps” approach to linearscaling selfconsistent field calculations is presented, which yields an extendedLagrangian formulation of ab initiomolecular dynamics. An exponential parametrization of the oneelectron density matrix, expressed in terms of atomcentered Gaussian basis functions, facilitates propagation along the manifold of density matrices in a geometrically correct fashion that automatically enforces idempotency constraints. The extended Lagrangian itself is constraint free, thus neither density matrix purification nor expensive, iterative solution for Lagrange multipliers is required. Propagation is highly efficient, and time steps compare favorably to those used in Car–Parrinello molecular dynamics simulations. The behavior of the method, especially with regard to the maintenance of adiabatic decoupling of nuclei and electrons, is examined for a sequence of diatomic molecules, and comparison is made to trajectories propagated on the converged Born–Oppenheimer surface. Certain claims to the contrary notwithstanding, our results demonstrate that vibrational frequencies may depend on the value of the fictitious mass parameter, even in an atomcentered basis. Lightatom stretching frequencies can be significantly redshifted, even when the nuclear and electronic energy scales are well separated. With a sufficiently small fictitious mass and a short time step, accurate frequencies can be obtained; we characterize appropriate values of these parameters for a wide range of vibrational frequencies.

A ground state potential energy surface for using Monte Carlo methods
View Description Hide DescriptionUsing variational Monte Carlo and a simple explicitly correlated wave function we have computed the Born–Oppenheimer energy of the ground state at 24 internuclear distances. We have also calculated the diagonal correction to the Born–Oppenheimer approximation and the lowestorder relativistic corrections at each distance using variational Monte Carlo techniques. The nonadiabatic values are evaluated from numerical derivatives of the wave function with respect to the nuclear coordinates. With this potential energy surface we have computed several of the lowest vibrationalrotational energies for this system. Our results are in good agreement with the best values found in the literature.

The curvature of the conical intersection seam: An approximate secondorder analysis
View Description Hide DescriptionWe present a method for analyzing the curvature (second derivatives) of the conical intersection hyperline at an optimized critical point. Our method uses the projected Hessians of the degenerate states after elimination of the two branching space coordinates, and is equivalent to a frequency calculation on a single Born–Oppenheimer potentialenergysurface. Based on the projected Hessians, we develop an equation for the energy as a function of a set of curvilinear coordinates where the degeneracy is preserved to second order (i.e., the conical intersection hyperline). The curvature of the potentialenergysurface in these coordinates is the curvature of the conical intersection hyperline itself, and thus determines whether one has a minimum or saddle point on the hyperline. The equation used to classify optimized conical intersection points depends in a simple way on the first and secondorder degeneracy splittings calculated at these points. As an example, for fulvene, we show that the two optimized conical intersection points of symmetry are saddle points on the intersection hyperline. Accordingly, there are further intersection points of lower energy, and one of symmetry—presented here for the first time—is found to be the global minimum in the intersection space.

Semiclassical Liouville method for the simulation of electronic transitions: Single ensemble formulation
View Description Hide DescriptionIn this paper, we describe a single ensemble implementation of the semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, one ensemble of trajectories supports the evolution of all semiclassical density matrix elements, rather than employing a distinct ensemble for each. The ensemble evolves classically under a single reference Hamiltonian, which is chosen based on physical grounds; for electronic relaxation of an initially excited state, the initially populated upper surface Hamiltonian is the natural choice. Classical trajectories evolving on the reference potential then represent the timedependent upper state population density and also the electronic coherence and the ground state density created by electronic transition. The error made in the classical motion of the trajectories for these latter distributions is compensated for by incorporating the difference between the correct and reference Liouville propagators into the calculation of the coefficients of the individual trajectories. This approach gives very accurate results for a number of model problems and cases describing ultrafast electronic relaxation dynamics.

The approximation in density matrix functional and wave function theory
View Description Hide DescriptionVarious energyfunctionals applying the approximation which leads to twoindex twoelectron integrals instead of fourindex twoelectron integrals in the electron–electron interaction term of the electronic energy are presented. Numerical results of multiconfiguration selfconsistent field calculations for the best possible wave function are compared to those obtained from the pair excitation multiconfiguration selfconsistent (PEMCSCF) method and two versions of density matrix functionaltheory. One of these is derived making explicit use of some necessary conditions for N representability of the secondorder density matrix. It is shown that this method models the energyfunctional based on the best possible wave function with good accuracy. The calculations also indicate that only a minor fraction of the total correlationenergy is incorporated by approaches for larger molecules.

Multiresolution quantum chemistry: Basic theory and initial applications
View Description Hide DescriptionWe describe a multiresolution solver for the allelectron local density approximation KohnSham equations for general polyatomic molecules. The resulting solutions are obtained to a userspecified precision and the computational cost of applying all operators scales linearly with the number of parameters. The construction and use of separated forms for operators (here, the Green’s functions for the Poisson and boundstate Helmholtz equations) enable practical computation in three and higher dimensions. Initial applications include the alkaliearth atoms down to strontium and the water and benzene molecules.

HEAT: High accuracy extrapolated ab initio thermochemistry
View Description Hide DescriptionA theoretical model chemistry designed to achieve high accuracy for enthalpies of formation of atoms and small molecules is described. This approach is entirely independent of experimental data and contains no empirical scaling factors, and includes a treatment of electron correlation up to the full coupledcluster singles, doubles, triples and quadruples approach. Energies are further augmented by anharmonic zeropoint vibrational energies, a scalar relativistic correction, firstorder spin–orbit coupling, and the diagonal Born–Oppenheimer correction. The accuracy of the approach is assessed by several means. Enthalpies of formation (at 0 K) calculated for a test suite of 31 atoms and molecules via direct calculation of the corresponding elemental formationreactions are within 1 kJ mol^{−1} to experiment in all cases. Given the quite different bonding environments in the product and reactant sides of these reactions, the results strongly indicate that even greater accuracy may be expected in reactions that preserve (either exactly or approximately) the number and types of chemical bonds.
 Gas Phase Dynamics and Structure: Spectroscopy, Molecular Interactions, Scattering, and Photochemistry

Photochemistry of ethylene: A multireference configuration interaction investigation of the excitedstate energy surfaces
View Description Hide DescriptionMultireference configuration interaction with singles and doubles (MRCISD) calculations have been performed for the optimization of conical intersections and stationary points on the ethylene excitedstate energy surfaces using recently developed methods for the computation of analytic gradients and nonadiabatic coupling terms. Basis set dependence and the effect of various choices of reference spaces for the MRCISD calculations have been investigated. The crossing seam between the and states has been explored in detail. This seam connects all conical intersections presently known for ethylene. Major emphasis has been laid on the hydrogenmigration path. Starting in the state of twistedorthogonal ethylene, a barrierless path to ethylidene was found. The feasibility of ethylidene formation will be important for the explanation of the relative yield of cis and trans elimination.

Investigating the exciplexes, Density functional approach
View Description Hide DescriptionPotential curves for the ground and the first lowest excited states of the (where exciplexes are calculated using the density functional theory(DFT) formalism. Relativistic calculations are carried out with and without spinorbit (SO) coupling effect, using a zeroth order regular approximation (ZORA) approach. The depth and position parameters of the potential curves for the case without spinorbit effect are presented and compared with other works. Potential curves for LiHe, CsHe, and FrHe without spinorbit effect and CsHe with the spinorbit effect are shown. A bond analysis is presented too, since the central interest of the study of the exciplexes is the nature and the existence of the bonding states that are the origin of emission spectra observed experimentally for these systems.

Study of ab initio molecular data for inelastic and reactive collisions involving the quasimolecule
View Description Hide DescriptionThe lowest two ab initiopotential energy surfaces (PES), and the corresponding nonadiabatic couplings between them, have been obtained for the system; the molecular data are compared to those calculated with the diatomic in molecules (DIM) method. The form of the couplings is discussed in terms of the topology of the molecular structure of the triatomic. The method of Baer is employed to generate “diabatic” states and the residual nonadiabatic couplings are calculated. The ab initio results for these are markedly different from the corresponding DIM data, and show the need to consider the third PES.

Zeeman effect in
View Description Hide DescriptionThe Zeeman effect in the excited state of CaF is measured and analyzed over a wide range of magnetic fields. It is found that the splitting of the Zeeman levels is largely determined by the coupling between different rotational states and there are no lowfield seeking states in the manifold of Zeeman levels at high magnetic fields. A model of the Zeemanspectrum based on the ligandfield theory of CaF is shown to be accurate in the interval of magnetic fields 0–5 Tesla. This demonstrates that the magnetic moment of the molecule is effectively determined by the spin angular momentum of a single electron and the orbital motion of the valence electron around the core. An analysis of the Zeemanspectrum as a function of the molecular rotational constant indicates that molecules should have significant rotational constants (at least as large as twice the rotational constant of CaF) to be amenable to magnetic trapping in high fields.

Photodissociation of laboratory oriented molecules: Revealing molecular frame properties of nonaxial recoil
View Description Hide DescriptionWe report the photodissociation of laboratory oriented OCS molecules. A molecular beam of OCS molecules is hexapole stateselected and spatially oriented in the electric field of a velocity map imaging lens. The oriented OCS molecules are dissociated at 230 nm with the linear polarization set at 45° to the orientation direction of the OCS molecules. The photofragments are quantum stateselectively ionized by the same 230 nm pulse and the angular distribution is measured using the velocity map imaging technique. The observed images are strongly asymmetric and the degree of asymmetry varies with the CO rotational state J. From the observed asymmetry in the laboratory frame we can directly extract the molecular frame angles between the final photofragment recoil velocity and the permanent dipole moment and the transition dipole moment. The data for CO fragments with high rotational excitation reveal that the dissociation dynamics is highly nonaxial, even though conventional wisdom suggests that the nearly limiting β parameter results from fast axial recoil dynamics. From our data we can extract the relative contribution of parallel and perpendicular transitions at 230 nm excitation.

Absolute electronimpact total ionization cross sections of chlorofluoromethanes
View Description Hide DescriptionAn experimental study is reported on the electronimpact total ionization cross sections (TICSs) of and molecules. The kinetic energy of the colliding electrons was in the 10–85 eV range. TICSs were obtained as the sum of the partial ionization cross sections of all fragment ions, measured and identified in a linear double focusing timeofflightmass spectrometer. The resulting TICS profiles—as a function of the electronimpact energy—have been compared both with those computed by ab initio and (semi)empirical methods and with the available experimental data. The computational methods used include the binaryencounterBethe (BEB) modified to include atoms with principal quantum numbers the Deutsch and Märk (DM) formalism, and the modified additivity rule (MAR). It is concluded that both modified BEB and DM methods fit the experimental TICS for and to a high accuracy, in contrast with the poor accord of the MAR method. A discussion on the factors influencing the discrepancies of the fittings is presented.

A density functional study on nitrogendoped carbon clusters
View Description Hide DescriptionUsing molecular graphics software, we designed numerous models of Geometry optimization and calculation on vibration frequency were carried out by the B3LYP density functional method. After comparison of structure stability, we found that the structures of groundstate and are bent chains with a nitrogen atom at either end, whereas when the groundstateclusters show three branches, each with a nitrogen atom located at the end. When the longest branch of is polyacetylenelike. When or 7, the longest branch is connected to the central carbon in a nonlinear manner. The with an even number of carbon atoms are more stable than those with odd numbers, matching the peak pattern observed in laserinduced mass spectra of The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electron affinities, and incremental binding energies.

The triplet state of cytosine and its derivatives: Electron impact and quantum chemical study
View Description Hide DescriptionThe excitation of the lowest electronic states and vibrational excitation of cytosine (C) have been studied using electron energy loss spectroscopy(EELS, 0–100 eV) with angular analysis. The singlet states have been found to be in good agreement with UVVIS absorption results on sublimed films, slightly blueshifted by about 0.1 eV. The EEL spectra recorded at residual energy below 2 eV show clear shoulders at energy losses of 3.50 and 4.25 eV (±0.1 eV). They are assigned to the lowest triplet electronic states of cytosine. Energies and molecular structures of the lowestlying triplet state of C and its methylated and halogenated and substituted derivatives F, Cl, and Br) have been studied using quantum chemical calculations with both molecular orbital and density functional methods, in conjunction with the and augccpVTZ basis sets. The tripletsinglet energy gap obtained using coupledcluster theory [CCSD(T)] and density functional theory(DFT) methods agrees well with those derived from EELS study. The first C’s vertical triplet state is located at 3.6 eV, in good agreement with experiment. The weak band observed at 4.25 eV is tentatively assigned to the second C’s vertical triplet excitation. For the substituted cytosines considered, the vertical triplet state is consistently centered at 3.0–3.2 eV above the corresponding singlet ground state but about 1.0 eV below the first excited singlet state. Geometrical relaxation involving outofplane distortions of hydrogen atoms leads to a stabilization of 0.6–1.0 eV in favor of the equilibrium triplet. The lowestlying adiabatic triplet states are located at 2.3–3.0 eV. Halogen substitution at both C(5) and C(6) positions tends to reduce the tripletsinglet separations whereas methylation tends to enlarge it. The vibrational modes of triplet cytosine and the ionization energies of substituted derivatives were also evaluated.

Lifetime of reactive scattering resonances: Qmatrix analysis and angular momentum dependence for the reaction by the hyperquantization algorithm
View Description Hide DescriptionWe report a study on the behavior with total angular momentumJ of several resonances occurring at collision energies below or slightly above the reaction barrier in the reaction.Resonance positions and widths are extracted from exact timeindependent quantum mechanical calculations using the hyperquantization algorithm and Smith’s Qmatrix formalism which exploits complete Smatrix information. The results confirm previous work but provide much greater insight. Identification of quasibound states responsible for the resonances based on adiabatic models for the longrange atom–molecule interactions both in the entrance and exit channels, is successful except for the feature occurring at the lowest energy, which is found to overlap with an exitchannel resonance for The two features are analyzed as overlapping resonances and their excellent Lorentzian fits, with wellbehaved Jdependences of positions and widths, support the interpretation of the lowenergy feature as a resonance to be associated to the triatomic transition state of the reaction.Resonance role on the reactive observables (integral cross sections and angular distributions) is investigated. The mechanism leading to forward scattering in the reactive differential cross section is commented, while the effects on rate constants, as well as the sensitivity of the resonancepattern to modification of the potential energy surface, are fully discussed elsewhere.