PubMed Health⌕ Search

Biomedical subjects

Pierre-Nicholas Roy

Publications and source records attributed to Pierre-Nicholas Roy.

13 recordsLinked to original sources
PubMed2007

Geometric constraints in semiclassical initial value representation calculations in Cartesian coordinates: excited states.

The authors show that a recently proposed approach [J. Chem. Phys. 123, 084103 (2005)] for the inclusion of geometric constraints in semiclassical initial value representation calculations can be used to obtain excited states of weakly bound complexes. Sample calculations are performed for free and constrained rare gas clusters. The results show that the proposed approach allows the evaluation of excited states with reasonable accuracy when compared to exact basis set calculations.

Journal Article↗
PubMed2006

Inclusion of inversion symmetry in centroid molecular dynamics: a possible avenue to recover quantum coherence.

Inversion symmetry is included in the operator formulation of the centroid molecular dynamics (CMD). This work involves the development of a symmetry-adapted CMD (SA-CMD), here particularly for symmetrization and antisymmetrization projections. A symmetry-adapted quasidensity operator, as defined by Blinov and Roy [J. Chem. Phys. 115, 7822 (2001)], is employed to obtain the centroid representation of quantum mechanical operators. Numerical examples are given for a single particle confined to one-dimensional symmetric quartic and symmetric double-well potentials. Two SA-CMD simulations are performed separately for both projections, and centroid position autocorrelation functions are obtained. For each projection, the quality of the approximation as well as the accuracy are similar to those of regular CMD. It is shown that individual trajectories from two separate SA-CMD simulations can be properly combined to recover trajectories for Boltzmann statistics. Position autocorrelation functions are compared to the exact quantum mechanical ones. This explicit account of inversion symmetry provides a qualitative improvement on the conventional CMD approach and allows the recovery of some quantum coherence.

Journal Article↗
PubMed2006

Path integral ground state study of finite-size systems: application to small (parahydrogen)N (N=2-20) clusters.

We use the path integral ground state method to study the energetic and structural properties of small para-H2 clusters of sizes ranging from 2 to 20 molecules. A fourth order formula is used to approximate the short imaginary-time propagator and two interaction potentials are considered. Our results are compared to those of exact basis set calculations and other quantum Monte Carlo methods when available. We find that for all cluster sizes considered, our results show a lower ground state energy than literature values obtained by diffusion Monte Carlo and variational Monte Carlo. For the dimer and trimer, ground state energies are in good agreement with exact results obtained using the discrete variable representation. Structural properties are found to be insensitive to the choice of interaction potential. We explore the use of Pekeris coordinates to analyze the importance of linear arrangement in trimers and for trimers within clusters of larger size.

Journal Article↗
PubMed2006

Excited States of weakly bound bosonic clusters: discrete variable representation and quantum Monte Carlo.

We compare two approaches for the accurate calculation of the energy levels of weakly bound boson trimers. The first approach is based on correlation function Monte Carlo employing optimized trial functions, while the second approach is based on the discrete variable representation. A trimer with atoms of half the mass of neon is used as a test problem for benchmark calculations. The two approaches yield identical results, within error bars, for all the J = 0 energy levels below the dissociation threshold. The relative merits of the two techniques are discussed, and a perspective is given for extension to larger clusters.

Journal Article↗
PubMed2005

Geometric constraints in semiclassical initial value representation calculations in Cartesian coordinates: accurate reduction in zero-point energy.

An approach for the inclusion of geometric constraints in semiclassical initial value representation calculations is introduced. An important aspect of the approach is that Cartesian coordinates are used throughout. We devised an algorithm for the constrained sampling of initial conditions through the use of multivariate Gaussian distribution based on a projected Hessian. We also propose an approach for the constrained evaluation of the so-called Herman-Kluk prefactor in its exact log-derivative form. Sample calculations are performed for free and constrained rare-gas trimers. The results show that the proposed approach provides an accurate evaluation of the reduction in zero-point energy. Exact basis set calculations are used to assess the accuracy of the semiclassical results. Since Cartesian coordinates are used, the approach is general and applicable to a variety of molecular and atomic systems.

Journal Article↗
PubMed2005

Path integral ground state with a fourth-order propagator: application to condensed helium.

Ground state properties of condensed helium are calculated using the path integral ground state (PIGS) method. A fourth-order approximation is used as short (imaginary) time propagator. We compare our results with those obtained with other quantum Monte Carlo (QMC) techniques and different propagators. For this particular application, we find that the fourth-order propagator performs comparably to the pair product approximation, and is far superior to the primitive approximation. Results obtained for the equation of state of condensed helium show that PIGS compares favorably to other QMC methods traditionally utilized for this type of calculation.

Journal Article↗
PubMed2004

Rotational spectrum, potential energy surface, and bound states of the weakly bound complex He-N2O.

Pure rotational transitions of the weakly bound complex He-N(2)O and three minor isotopomers (He-(14)N(15)NO, He-(15)N(14)NO, and He-(15)N(15)NO) were measured in the frequency region from 6 to 20 GHz. Predictions for the microwave transition frequencies were based on the infrared work by Tang and McKellar [J. Chem. Phys. 117, 2586 (2002)]. In the case of (14)N containing isotopomers, nuclear quadrupole hyperfine structure of the rotational transitions was observed and analyzed. The resulting spectroscopic parameters were used to determine geometrical and dynamical information about the complex. An ab initio potential energy surface was calculated at the coupled cluster level of theory with single and double excitations and perturbative inclusion of triple excitations. This surface was constructed using the augmented correlation consistent polarized valence triple zeta basis set for all atoms with the inclusion of bond functions for the van der Waals bond. Bound state calculations were done to determine the energies of low-lying rovibrational levels that are supported by the potential energy surface. The resulting transition energies agree with the experimental values to 1% or better.

Journal Article↗
PubMed2004

Energy levels and wave functions of weakly-bound 4Hex 20NeyH (x+y=2) systems using Pekeris coordinates and a symmetry-adapted Lanczos approach.

Energy levels and wave functions of floppy triatomic rare gas hydrides are calculated using a Pekeris coordinate system and the importance of various triangular configurations is assessed through the calculation of reduced distribution functions and relative weights. The calculations are performed using a symmetry-adapted Lanczos recursion within the discrete variable representation. For the 4He2H- anion, the present results are compared with those obtained from calculations based on other methods, and the accuracy of the present method is discussed. Calculations are also performed for the case of 4He2H and 20Ne2H, as well as for the mixed 4He 20NeH neutrals. Our results show that no bound states are found for 4He2H while only one bound state is found for both the 20Ne2H and 4He 20NeH complexes. Interestingly, a very important and common property of these systems is that there is a significant contribution from linear configurations to their bound states.

Journal Article↗
PubMed2004

Quantum Monte Carlo study of helium clusters doped with nitrous oxide: quantum solvation and rotational dynamics.

Dynamical and structural properties of small (4)He(N)-N(2)O complexes have been analyzed using ground-state and finite-temperature Monte Carlo simulations. The effective rotational constants resulting from the ground-state calculations are in excellent agreement with the results of a recent spectroscopic study [Y. Xu et al., Phys. Rev. Lett. 91, 163401 (2003)]. After an initial decrease for cluster sizes up to N=8, the rotational constant increases, signaling a transition from a molecular complex to a quantum-solvated system. Such a turnaround is not present in the rotational constants extracted from the finite-temperature Monte Carlo calculations, performed for Boltzmann statistics, thus highlighting the importance of exchange effects to explain the decoupling between a solvated dopant and the helium motion.

Journal Article↗
PubMed2004

Path integral Monte Carlo approach for weakly bound van der Waals complexes with rotations: algorithm and benchmark calculations.

A path integral Monte Carlo technique suitable for the treatment of doped helium clusters with inclusion of the rotational degrees of freedom of the dopant is introduced. The extrapolation of the results to the limit of infinite Trotter number is discussed in detail. Benchmark calculations for small weakly bound (4)He(N)--OCS clusters are presented. The Monte Carlo results are compared with those of basis set calculations for the He--OCS dimer. A technique to analyze the orientational imaginary time correlation function is suggested. It allows one to obtain information regarding the effective rotational constant for a doped helium cluster based on a model for the rotational Hamiltonian. The renormalization of the effective rotational constant for (4)He(N)--OCS clusters derived from the orientational imaginary time correlation function is in good agreement with experimental results.

Journal Article↗
PubMed2004

On the calculation of single-particle time correlation functions from Bose-Einstein centroid dynamics.

The calculation of single-particle time correlation functions using the Bose-Einstein centroid dynamics formalism is discussed. A new definition of the quasidensity operator is used to calculate the centroid force on a given particle for an anharmonic system. The force includes correlation effects due to quantum statistics and is used for the calculation of the classical-like dynamics of phase-space centroid variables within the centroid molecular dynamics approximation. Time correlation functions are then obtained for single-particle quantities. These correspond to the double-Kubo transform of exact quantum-mechanical correlation functions. The centroid dynamics results are compared to those of exact basis-set calculations and a good agreement is found. The level of accuracy is in fact the same as what was observed earlier for the calculation of center-of-mass correlation functions for Fermi-Dirac and Bose-Einstein statistics, and for any correlation function for Boltzmann statistics. These results show that it is now possible to use Bose-Einstein centroid molecular dynamics to calculate single-particle correlation functions for systems where quantum exchange effects are present.

Journal Article↗
PubMed2004

Connection between the observable and centroid structural properties of a quantum fluid: application to liquid para-hydrogen.

It is shown that the discrepancy between path integral Monte Carlo [M. Zoppi et al., Phys. Rev. B 65, 092204 (2002)] and path integral centroid molecular dynamics [F. J. Bermejo et al., Phys. Rev. Lett. 84, 5359 (2000)] calculations of the static structure factor of liquid para-hydrogen can be explained based on a deconvolution equation connecting centroid and physical radial distribution functions. An explicit expression for the kernel of the deconvolution equation has been obtained using functional derivative techniques. In the superposition approximation, this kernel is given by the functional derivative of the effective potential with respect to the pairwise classical potential. Results of path integral Monte Carlo calculations for the radial distribution function and the static structure factor of liquid para-hydrogen are presented.

Journal Article↗