Journal
COMPUTER PHYSICS COMMUNICATIONS
Volume 180, Issue 11, Pages 2025-2033Publisher
ELSEVIER SCIENCE BV
DOI: 10.1016/j.cpc.2009.06.004
Keywords
Schrodinger equation; Molecular dynamics; Projection operators; Finite difference methods; Grid computing
Funding
- European Union ToK [MTKD-CT-2005-029583]
- COST program [D37]
Ask authors/readers for more resources
A method for solving the Schrodinger equation of N-atom molecules in 3N - 3 Cartesian coordinates usually defined by Jacobi vectors is presented. The separation and conservation of the total angular momentum are obtained not by transforming the Hamiltonian in internal curvilinear coordinates but instead, by keeping the Cartesian formulation of the Hamiltonian operator and projecting the initial wavefunction onto the proper irreducible representation angular momentum subspace. The increased number of degrees of freedom from 3N - 6 to 3N - 3, compared to previous methods for solving the Schrodinger equation, is compensated by the simplicity of the kinetic energy operator and its finite difference representations which result in sparse Hamiltonian matrices. A parallel code in Fortran 95 has been developed and tested for model potentials of harmonic oscillators. Moreover, we compare data obtained for the three-dimensional hydrogen molecule and the six-dimensional water molecule with results front the literature. The availability of large clusters of computers with hundreds of CPUs and GBytes of memory. as well as the rapid development of distributed (Grid) computing, make the proposed method. which is unequivocally highly demanding in memory and Computer time, attractive for studying Quantum Molecular Dynamics. (C) 2009 Elsevier B.V. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available