Search

The potential energy surface of H(5)(+) is characterized using density functional theory. The hypersurface is evaluated at selected configurations employing different functionals, and compared with results obtained from ab initio CCSD(T) calculations. The lowest ten stationary points (minima and saddle-points) on the surface are located, and the features of the short-, intermediate-, and long-range intermolecular interactions are also investigated. A detailed analysis of the surface's topology, and comparisons with extensive CCSD(T) results, as well as a recent ab initio analytical surface, shows that density functional theory calculations using the B3(H) functional represent very well all aspects studied on the H(5)(+) potential. These include the tiny energy difference between the minimum at 1-C(2v) configuration and the 2-D(2d) one corresponding to the transition state for the proton transfer between the two equivalent C(2v) minima, and also the correct asymptotic behavior of the long-range interactions. The calculated binding energy and dissociation enthalpies compare very well with previous benchmark coupled-cluster ab initio data, and with experimental data available. Based on these results the use of such approach to perform first-principles molecular dynamics simulations could provide reliable information regarding the dynamics of protonated hydrogen clusters.

The Ar(3) system has been studied between T=0 K and T=40 K by means of a path-integral Monte Carlo (PIMC) method. The behavior of the average energy in terms of the temperature has been explained by comparison with results obtained with the thermal averaged rovibrational spectra estimated via: (i) a quantum mechanical method based on distributed Gaussian functions for the interparticle distances and (ii) an analytical model which precisely accounts for the participation of the dissociative continua Ar(2)+Ar and Ar+Ar+Ar. Beyond T approximately 20 K, the system explores floppier configurations than the rigid equilateral geometry, as linear and Ar-Ar(2)-like arrangements, and fragmentates around T approximately 40 K. A careful investigation of the specific heat in terms of a confining radius in the PIMC calculation seems to discard a proper phase transition as in larger clusters, in apparent contradiction with previous reports of precise values for a liquid-gas transition. The onset of this noticeable change in the dynamics of the trimer occurs, however, at a remarkably low value of the temperature in comparison with Ar(n) systems formed with more Ar atoms. Quantum mechanical effects are found of relevance at T

Recently, the predissociation lifetimes of the NeBr(2)(B) complex for different initial vibrational excitation (10

Pump-probe results are reported for NeCl(2) excited to the Cl(2) B state, undergoing vibrational predissociation, and then probed via E <-- B transitions. Intensities, lifetimes and product vibrational branching ratios are reported for 16 < or = v' < or = 19 Cl(2) stretching quanta. The intensity of the signal rapidly decreases above v' = 17. Detailed wave packet calculations of the vibrational predissociation dynamics are performed to determine if the experimental results can be explained by the onset of IVR dynamics. The calculations and the experiment are in close accord for low vibrational levels. For higher levels, some, but not all, of the loss of experimental signal can be attributed to IVR. To test whether electronic relaxation dynamics are important for NeCl(2) and ArCl(2), excited state potential surfaces that incorporate spin orbit coupling effects are calculated. These surfaces are then used in a wave packet calculation that includes both vibrational predissociation and electronic predissociation dynamics. The results show that electronic predissociation is important for ArCl(2) levels above v' = 12. For NeCl(2) the calculation suggests that the onset of electronic predissociation should occur for levels as low as v' = 13 but may not contribute markedly to the observed loss of signal above v' = 17. Suggestions are made for further studies of this puzzling problem.

A global potential energy surface is obtained for the ground state of the endoergic Au((2)S)+H(2)(X (1)Sigma(g) (+))-->AuH((1)Sigma(+))+H((2)S) reaction. The global potential is obtained by fitting highly correlated ab initio calculations on the system, using relativistic pseudopotential for the gold atom. Several electronic states are calculated correlating with Au((2)S)+H(2), Au((2)D)+H(2), and H(2), Au((2)P)+H(2) asymptotes. These states show several conical intersections and curve crossings along the minimum energy reaction path which are analyzed in detail. One of them gives rise to an insertion well in which there are important contributions from the Au((2)D) and Au((2)P) states of gold, which is interesting because it is analog to the deep chemisorption well appearing in larger gold clusters. Quantum wave packet and quasiclassical trajectory dynamical calculations performed for the reaction at zero total angular momentum are in good agreement, provided that a Gaussian binning method is used to account for the zero-point energy of products. Finally, integral and differential cross sections are calculated for the reaction with quasiclassical trajectories. Two different reaction mechanisms are found, one direct and the second indirect, in which the Au atom inserts in between the two hydrogen atoms because of the existence of the insertion well discussed above.

A new method is proposed to partition the density of a system in two portions. The density on each subsystem is the solution of a Fock equation modified by the addition of an embedding potential. This embedding potential is obtained iteratively by minimizing the difference between the electronic densities of the total system and the sum of the subsystems. Thus, the electronic density partition and the embedding potential are obtained at the same time within the procedure, guaranteeing the v-representability of the densities partitioned. This fact is a considerable improvement of a recently proposed embedding potential inversion technique, [O. Roncero, M. P. de Lara-Castells, P. Villarreal, F. Flores, J. Ortega, M. Paniagua, and A. Aguado, J. Chem. Phys. 129, 184104 (2008)], in which the embedding potential is obtained once the electronic density is previously partitioned. The method is first applied to a linear H(10) chain to illustrate how it works. The orbitals obtained are localized on each subsystem, and can be used to include local electronic correlation with currently available ab initio programs. Finally, the method is applied to include the electronic correlation needed to describe the van der Waals interaction between H(10) chains and H(2) molecules, of approximately 12 meV, giving very accurate results.

A detailed study of the rovibrational spectrum of the Ar trimer is performed by means of an exact hyperspherical coordinate (HC) method and a variational approach based on distributed Gaussian functions (DGFs) to describe the interparticle distances. The good agreement observed between the energy levels obtained with both procedures for high values of the total angular momentum (J=15 and 20) reveals the quality of the DGF method to describe the rotation of the title system. Rotational constants for the lowest bound states, obtained as averages for each vibrational state, have been obtained and compared to previous results. A detailed analysis of density probability functions obtained by means of the HC approach for rovibrational states at J=0 and 20 shows close similitudes thus supporting the vibration-rotation separation adopted within the DGF scheme for the Ar(3) system.

A new embedding method to include local correlation in large systems is proposed. In this method the density of the whole system, calculated via density functional theory approaches, is partitioned in two pieces, one corresponding to the subsystem of interest and the rest to the environment. In the second step, an embedding potential is obtained iteratively using as a driving force the self-repulsion due to the density difference, in a similar form as proposed by Zhao et al. [Phys. Rev. A 50, 2138 (1994)], to obtain the "exact" exchange-correlation functional. Such potential is added to the Fock equation to build the localized molecular orbitals which are further used to include the local electronic correlation in the subsystem of interest. This method is an alternative to the previous DFT-based embedding methods first proposed by Wesolowski and Washell [J. Phys. Chem. 97, 8050 (1993)] and after enhanced by Govind et al. [J. Chem. Phys. 110, 7677 (1999)] and adapted to metal extended systems, which use density functionals to describe the kinetic energy contribution to the embedding potential, whose precise form has been largely treated in the literature and its crucial role is discussed here. The method is applied to hydrogen chains and its van der Waals interaction with H(2). The results obtained are in very good agreement with exact calculations performed on the whole system, which demonstrates that the method proposed is a very promising route to introduce correlation in large systems.

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.

We report the first spectroscopic observation of and vibrational predissociation dynamics for ArBr(2). Measurements are reported for the linear and T-shaped isomers with time and frequency resolution near the Fourier limit of a 15 ps pulse. For the T-shaped isomer, the Ar-Br(2) bond energy D(0) for the B state, nu(')=19, is found to be 200 cm(-1), yielding a D(0)(X) value of 213 cm(-1). Product appearance rates are determined for nu(')=16-25 of the B state of ArBr(2). While the rate generally increases with increasing vibrational quantum number, the trend is not monotonic. Also, obtaining reproducible rates for any given vibrational level requires very careful control of the experimental conditions. The data suggest that ArBr(2) undergoes vibrational predissociation (VP) in the sparse intramolecular vibrational relaxation regime. These observations are consistent with theoretical calculations that show that T-shaped ArBr(2) undergoes VP in the sparse regime, such that lifetimes are strongly dependent on both the vibrational and rotational quantum numbers. As for ArI(2), a linear isomer of ArBr(2) is found to contribute a quasicontinuous background to the excitation spectra. Direct excitation of this continuum results in a very broad product vibrational distribution.

An energy-based method is proposed for the diabatization of the OH(2Pi)+F(2P)-->O(3P)+HF(1Sigma+) reaction. It is demonstrated that the diabatic representation obtained is regularized, i.e., the residual derivative couplings do not present singularities at the conical intersections appearing along the reaction path. This method only requires the knowledge of the 1,2 3A" and 1 3A' eigenvalues and does not require any adjustable parameter. Thus, many convergence problems arising in other derivative-based diabatization methods are avoided, and the description of the configuration space along the reaction path is enormously simplified. Three-dimensional coupled diabatic energy surfaces are obtained by an interpolation procedure using approximately 4000 accurate ab initio points. The angular resolved photodetachment cross sections are obtained in the diabatic and adiabatic representations using a wave packet method. An excellent agreement is obtained with recent experimental data [D. M. Neumark, Phys. Chem. Chem. Phys. 7, 433 (2005)] for high electron kinetic energies where only the triplet electronic states contribute.

The H+ + H2 exchange reaction has been studied theoretically by means of a different variety of methods as an exact time independent quantum mechanical, approximate quantum wave packet, statistical quantum, and quasiclassical trajectory approaches. Total and state-to-state reaction probabilities in terms of the collision energy for different values of the total angular momentum obtained with these methods are compared. The dynamics of the reaction is extensively studied at the collision energy of E(coll)=0.44 eV. Integral and differential cross sections and opacity functions at this collision energy have been calculated. In particular, the fairly good description of the exact quantum results provided by the statistical quantum method suggests that the dynamics of the process is governed by an insertion mechanism with the formation of a long-lived collision complex.

A procedure for the transformation from reactant to product Jacobi coordinates is proposed, which is designed for the extraction of state-to-state reaction probabilities using a time-dependent method in a body-fixed frame. The method consists of several steps which involve a negligible extra computational time as compared with the propagation. Several intermediate coordinates are used, in which the efficiency depends on the masses of the atoms involved in the reaction. A detailed study of the relative efficiency of using reactant and product Jacobi coordinates is presented for several systems, and simple arguments are found depending on the masses of the atoms involved in the reaction. It is found that the proposed method is, in general, more efficient than the use of product Jacobi coordinates, specially for nonzero total angular momentum. State-to-state reaction probabilities are obtained for Li+FH-->LiF+H and F+HO-->FH+O collisions for several total angular momenta.

The H(+)+D(2) and D(+)+H(2) reactive collisions are studied using a recently proposed adiabatic potential energy surface of spectroscopic accuracy. The dynamics is studied using an exact wave packet method on the adiabatic surface at energies below the curve crossing occurring at approximately 1.5 eV above the threshold. It is found that the reaction is very well described by a statistical quantum method for a zero total angular momentum (J) as compared with the exact ones, while for higher J some discrepancies are found. For J >0 different centrifugal sudden approximations are proposed and compared with the exact and statistical quantum treatments. The usual centrifugal sudden approach fails by considering too high reaction barriers and too low reaction probabilities. A new statistically modified centrifugal sudden approach is considered which corrects these two failures to a rather good extent. It is also found that an adiabatic approximation for the helicities provides results in very good agreement with the statistical method, placing the reaction barrier properly. However, both statistical and adiabatic centrifugal treatments overestimate the reaction probabilities. The reaction cross sections thus obtained with the new approaches are in rather good agreement with the exact results. In spite of these deficiencies, the quantum statistical method is well adapted for describing the insertion dynamics, and it is then used to evaluate the differential cross sections.

Global three-dimensional adiabatic potential-energy surfaces for the excited 2(3)A" and 1(3)A' triplet states of OHF are obtained to study the F(2P)+OH(2pi)-->O(3P)+HF(1sigma+) reaction. Highly accurate ab initio calculations are obtained for the two excited electronic states and fitted to analytical functions with small deviations. The reaction dynamics is studied using a wave-packet treatment within a centrifugal sudden approach, which is justified by the linear transition state of the two electronic states studied. The reaction efficiency presents a marked preference for perpendicular orientation of the initial relative velocity vector and the angular momentum of the OH reagent, consistent in the body-fixed frame used with an initial collinear geometry which facilitates the access to the transition state. It is also found that the reaction cross section presents a rather high threshold so that, in an adiabatic picture, the two excited triplet states do not contribute to the rate constant at room temperature. Thus, only the lowest triplet state leads to reaction under these conditions and the simulated rate constants are too low as compared with the experimental ones. Such disagreement is likely to be due to nonadiabatic transitions occurring at the conical intersections near the transition state for this reaction.

Ca+HCl(upsilon,j) reactive collisions were studied for different rovibrational states of the HCl reactant using wave-packet calculations in reactant Jacobi coordinates. A recently proposed potential-energy surface was used with a barrier of approximately 0.4 eV followed by a deep well. The possibility of an insertion mechanism due to this last well has been analyzed and it was found that once the wave packet passes over the barrier most of it goes directly to CaCl+H products, which shows that the reaction dynamics is essentially direct. It was also found that there is no significant change in the reaction efficiency as a function of the initial HCl rovibrational state, because CaHCl at the barrier has an only little elongated HCl bond. Near the threshold for reaction with HCl(upsilon=0), however, the reaction shows significant steric effects for j > 0. In a complementary study, the infrared excitation from the Ca-HCl van der Waals well was simulated. The spectrum thus obtained shows several series of resonances which correspond to quasibound states correlating to excited HCl(upsilon) vibrations. The Ca-HCl binding energies of these quasibound states increase dramatically with upsilon, from 75 to 650 cm(-1), because the wave function spreads increasingly over larger HCl bond lengths. Thus it explores the region of the barrier saddle point and the deep insertion well. Although also the charge-transfer contribution increases with upsilon, the reaction probability for resonances of the upsilon=2 manifold, which are well above the reaction threshold, is still negligible. This explains the relatively long lifetimes of these upsilon=2 resonances. The reaction probability becomes significant at upsilon=3. Our simulations have shown that an experimental study of this type will allow a gradual spectroscopic probing of the barrier for the reaction.

The potential-energy surface of the ground electronic state of CaHCl has been obtained from 6400 ab initio points calculated at the multireference configuration-interaction level and represented by a global analytical fit. The Ca+HCl-->CaCl+H reaction is endothermic by 5100 cm(-1) with a barrier of 4470 cm(-1) at bent geometry, taking the zero energy in the Ca+HCl asymptote. On both sides of this barrier are potential wells at linear geometries, a shallow one due to van der Waals interactions in the entrance channel, and a deep one attributed to the H(-)Ca(++)Cl(-) ionic configuration. The accuracy of the van der Waals well depth, approximately 200 cm(-1), was checked by means of additional calculations at the coupled-cluster singles and doubles with perturbative triples level and it was concluded that previous empirical estimates are unrealistic. Also, the electric dipole function was calculated, analytically fitted in the regions of the two wells, and used to analyze the charge shifts along the reaction path. In the insertion well, 16,800 cm(-1) deep, the electric dipole function confirmed the ionic structure of the HCaCl complex and served to estimate effective atomic charges. Finally, bound rovibrational levels were computed both in the van der Waals well and in the insertion well, and the infrared-absorption spectrum of the insertion complex was simulated in order to facilitate its detection.

The competition between vibrational and electronic predissociations of the ArI2(B) van der Waals complex has been studied using several dynamical computational methods: exact quantum wave-packet propagation, time-dependent golden rule, and quasiclassical trajectory with quantum jumps model. Five electronic states are considered using recent three-dimensional coupled surfaces obtained with a perturbative diatoms-in-molecules method. Final vibrational and electronic populations, predissociation rates, and absorption spectra have been computed for I2(B,v=18-24)<--I2(X,v=0) excitations within the complex. The contribution of vibrational predissociation into the total decay oscillates as a function of vibrational excitation due to intramolecular vibrational relaxation in a sparse-intermediate regime, which induces irregular variations of the total decay rate. Franck-Condon oscillations control the branching ratios of the individual electronic predissociation channels. However, since these oscillations are out of phase as a function of vibrational excitation, they have limited effect on the oscillatory behavior of the total predissociation rate. Comparison between exact quantum and perturbative golden rule calculations shows that vibrational predissociation has some impact on the electronic predissociation process and affects the final electronic distributions. On the contrary, vibrational product distributions are not significantly affected by the electronic predissociation. A classical description of the ArI2 dynamics provides an averaged picture of the competing predissociation processes, being better adapted for treating intermolecular vibrational relaxation in the statistical limit., ((c) 2005 American Institute of Physics.)

Wave packet simulations of the photodetachment spectrum of OHF(-) are performed on several electronic adiabatic states, two triplets and four singlets of neutral OHF. The transition moments to these six states have been approximated using the ab initio electronic wave functions of OHF(-) and OHF calculated at the equilibrium configuration of the parent anion. In a first step, two-dimensional simulations of the spectrum are performed on new two-dimensional potential energy surfaces (PESs) of the neutral in a OHF collinear geometry. The resulting simulated spectrum is in rather good agreement with the experimental one, reproducing all the structures from 0 to 2.5 eV electron kinetic energies. At energies below 0.5 eV, all calculated states, singlets and triplets, contribute to the total spectrum. At higher energies, however, only the triplet states participate. In a second step, to improve the description of the spectrum, three-dimensional wave packet simulations of the spectrum are performed, getting an excellent agreement with the experiment. The collinear (3)Sigma(-) and (3)Pi states split in two (3)A(") and one (3)A('). New adiabatic PESs are used in this work for the 2 (3)A(") and 1 (3)A(') states, while the one recently proposed was used for the ground 1 (3)A("). It is found that the minimum energy paths of the (3)Sigma(-) and (3)Pi states cross twice at collinear geometry, so that at the transition state the ground state corresponds to (3)Pi, while (3)Sigma(-) is the lowest state otherwise. Such conical intersections are expected to give rise to important Sigma-Pi vibronic effects, requiring a complete three-dimensional model of coupled diabatic states to improve our understanding of the reaction dynamics in this kind of systems., (Copyright 2004 American Institute of Physics.)

A theoretical study of the F(2P) + OH(2Pi) --> HF(1Sigma+) + O(3P) reactive collisions is carried out on a new global potential energy surface (PES) of the ground 3A" adiabatic electronic state. The ab initio calculations are based on multireference configuration interaction calculations, using the aug-cc-pVTZ extended basis sets of Dunning et al. A functional representation of the PES shows no nominal barrier to reaction, contrary to previous results by others. Wave packet and quasiclassical trajectory calculations have been performed for this PES to study the F + OH(v = 0,j) reactive collision. The comparison was performed at fixed and constant values of the total angular momentum from 0 to 110 and relative translational energy up to 0.8 eV. The reaction presents a dynamical barrier, essentially due to the zero-point energy for the bending vibration near the saddle point. This determines two different reaction mechanisms. At energies higher than approximately 0.125 eV the reaction is direct, while below that value it is indirect and mediated by heavy-light-heavy resonances. Such resonances, also found in the simulations of the photodetachment spectrum of the triatomic anion, manifest themselves in the quasiclassical simulations, too, where they are associated to periodic orbits., (Copyright 2004 American Institute of Physics)

In this work a simulation of the OHF(-) photodetachment spectrum is performed in a three-dimensional potential energy surface recently developed for OHF((3)A(")). The ground (2)A(') state potential of the anion is calculated in three dimensions based on accurate ab initio calculations and the reaction dynamics is studied using a wave packet method. The calculated spectrum shows a sequence of bands associated to vibrational HF(v) up to v=3. Each band is formed by a continuous spectrum and resonant structures. These resonances are associated to the OH-F channel well of the (3)A(") PES, in which fragmentation occurs through vibrational predissociation. Above the OH(v=0) threshold a new resonant pattern appears corresponding to heavy-light-heavy resonances. Special attention is paid to the assignment of these resonances because they mediate the reaction dynamics in the OH+F collision at low kinetic energies. The sequence of bands is in rather good agreement with that appearing in the experimental spectrum, especially at higher electron kinetic energies. At low kinetic energies, however, some other electronic states may contribute. The resonance structures might be washed out by the rotational average and the relatively low energy resolution of the experiment., ((c) 2004 American Institute of Physics.)

Aim: A Spanish Consensus Conference had been arranged by the Spanish Club for the Study of Helicobacter pylori to encourage the use of eradication regimens in appropriate indications. The aim of our study was to assess whether the publication of these recommendations in November 1999 induced a change on the indications of eradication in the referring primary care district., Method: Patients who had undergone Helicobacter pylori eradication therapy indicated by a gastroenterologist and referred to perform the 13C-labelled marked urea breath test were evaluated. Indications of eradication therapy were analysed and divided in: a) use of Conference recommended procedures (duodenal or gastric ulcer, erosive duodenitis, gastric MALT lymphoma and after resection of gastric adenocarcinoma); and b) other procedures not included in the recommendations (the rest). Indications established during former and latter years of publication of the Club's recommendations were compared., Results: A total of 659 established eradications in 1999 against 537 in 2000, were studied. Regarding established eradications in 1999, 418 out of 659 (63%) followed the recommendations made by the Spanish Club, while 241 (37%) did not. In 2000, the proportion of eradication therapy administered using the recommendations made by the Spanish Club was similar to previous year: 338 (63%) followed the recommendations, against 199 (37%) that did not (p = not significant). In both years, considering those indications not recommended, patients either with superficial chronic gastritis with or without dyspepsia or with gastrooesophageal reflux disease were the main disorders to indicate eradication therapy., Conclusions: The recommendations of the Spanish Helicobacter pylori Study Club have not triggered a significant change after their publication on the indications for Helicobacter pylori eradication therapy in the setting studied.