320 likes | 418 Views
Lanczos Representation Methods in Application to Rovibrational Spectroscopy Calculations. Hong Zhang and Sean Smith Quantum & Molecular Dynamics Group Center for Computational Molecular Science The University of Queensland, Australia. 1. Lanczos representation methods.
E N D
Lanczos Representation Methods in Application to Rovibrational Spectroscopy Calculations Hong Zhang and Sean Smith Quantum & Molecular Dynamics Group Center for Computational Molecular ScienceThe University of Queensland, Australia
1. Lanczos representation methods 1. Lanczos representation filter diagonalisation in application to the calculations of ro-vibrational spectroscopy and unimolecular dissociation resonances, which includes recent J > 0 calculations and parallel computing implementation. 2. TI Lanczos subspace wavepacket method; real Lanczos artificial boundary inhomogeneity (LABI); real single Lanczos propagation ABI in bimolecular reactions. 3. Comparison with real Chebyshev wavepacket method.
2. Methodologies 1. Transform the primary representation (DVR or FBR) of the fundamental equations (e.g., TI Schrödinger equation or TI wavepacket-Lippmann-Schwinger equation) into the tridiagonal Lanczos representation. 2. Solve the eigen-problem orlinear systemwithin the subspace to obtain either eigen-pairs or subspace wavefunctions. 3. Extract physical information, e.g., bound state energies, resonance energies and widths, S matrix elements in reactive scattering.
Lanczos representation transformation [Kouri, Arnold & Hoffman, CPL, 203: 166, 1993; Zhang H, Smith SC, JCP, 116: 2354, 2002]
Algorithm (a) Choose Mth element of (Ej) to be arbitrary (but non-zero) and calculate (b) For k=M-1, M-2, …, 2, update scalar k-1: (c) Determine constant (true= c) by normalization or by
Real Chebyshev wavepacket propagation (1) [Zhang H, Smith SC, JCP, 117: 5174, 2002]
Generalized eigenvalue equation: WB=SB Sjj’=((Ej)|(Ej’)) Wjj’=((Ej)|TM|(Ej’)) Lanczos Subspace Filter Diagonalization
Time evolution operator Low Storage Filter Diagonalization
J > 0 Hamiltonian: expansion Multiplying the Hamiltonian on the left side by and utilizing some integral and differential formula, one can obtain the coupled equations of motion in internal coordinates
J > 0 Hamiltonian: coupled equations m = 1 for m = 0 for
J > 0 Hamiltonian: angular DVR 1) Diagonalizing the coordinate operator to obtain the DVR points and transformation matrix: 2) Gauss-Jacobi quadrature
Representation Parallel computing Final analysis Conceptually difficult Propagation Computationally difficult: cpu time and memory In propagation, the most time consuming part is the matrix-vector multiplication. We use Message-passing interface (MPI) to perform parallel computation.
a) Master processor (ID = 0) Perform the main propagation; write auto-correlation functions or , elements; all other related works except the matrix-vector multiplications. b) Working processors Perform the matrix-vector multiplications for each component. Processor assignment
Communications According to the Coriolis coupling rules, only two nearest neighbouring components need to communicate. Load balancing jmin is different for each component, but jmax is the same, i.e, the DVR size for angle is changeable. Thus some processors might need to wait for others. For the highest or the lowest components, only one Coriolis coupling required.
Timing • Due to the communications and loading balance issues, the model doesn’t scale ideally with (J+1) for even spectroscopic symmetry or J for odd spectroscopic symmetry. • However, one can achieve wall clock times (e.g., for even symmetry J = 6 HO2 case) that are within about a factor of 2 of J = 0 calculations. For non parallel computing, the wall clock times will approximately be a factor of 7 of J = 0 calculations.
3. Recent Results • HO2 (J = 0-50): using both Lanczos and Chebyshev method from parallel computing. • HOCl/DOCl (J = 0-30): using parallel Lanczos method.
n LFD CFD Ref 1 .000271 .000270 .000270 10 .295249 .295249 .295272 20 .447875 .447876 .447903 30 .574994 .574994 .575033 Table 1 Selected HO2 bound state energies in eV for J = 1 (even symmetry) for comparison. (Ref: X. T. Wu and E. F. Hayes, JCP, 1997,107, 2705, LFD: Zhang H, Smith SC, JCP, 2003, 118, 10042)
Fig. 1 J = 6 (even symmetry) bound state energies for HO2 using both Lanczos FD and Chebyshev FD method. Circles from Lanczos FD method while squares from Chebyshev FD method.
Table 2 Selected HO2 bound state energies for J = 30 (even symmetry) for comparison.All energy units are in eV (Zhang & Smith, JPCA, 110: 3246, 2006 ).
Table 3 Calculated HOCl ro-vibrational state energies in cm-1 with spectroscopic assignments for J = 20 (Zhang, Smith, Nanbu & Nakamura,JPCA, 2006, 110: 5468).
Table 4 Comparison of experiments and calculations for selected HOCl far infrared transitions in cm-1 (Zhang, Smith, Nanbu & Nakamura, JPCA, 2006, , 110: 5468).
Fig. 2 Parallel computation of resonances for J = 20, 30, 40 and 50 cases in HO2 dissociation. The quantum average rates (green solid line) are compared with statistical/classical results of Troe. et. al.
Development of Lanczos representation methods; Design of a parallel computing model; Combination of both has made rigorous quantum calculations possible for challenging J > 0 applications. 4. Conclusions
Acknowledgements • Prof. Nanbu, Prof. Nukamura and CC members; • Prof. Sean Smith and CCMS members; • Prof. R. Kosloff and Prof. V. A. Mandelshtam; • JSPS Invitation Fellowship; Australian Research Council; • Supercomputer time from APAC, UQ, CC, IMS.