Complex generalized minimal residual algorithm for iterative solution of quantum-mechanical reactive scattering equations

A complex GMRes (generalized minimum residual) algorithm is presented and used to solve dense systems of linear equations arising in variational basis‐set approaches to quantum‐mechanical reactive scattering. The examples presented correspond to physical solutions of the Schrodinger equation for the reactions O+HD→OH+D, D+H2→HD+H, and H+H2→H2+H. It is shown that the computational effort for solution with GMRes depends upon both the dimension of the linear system and the total energy of the reaction. In several cases with dimensions in the range 1110–5632, GMRes outperforms the LAPACK direct solver, with speedups for the linear equation solution as large as a factor of 23. In other cases, the iterative algorithm does not converge within a reasonable time. These convergence differences can be correlated with ‘‘indices of diagonal dominance,’’ which we define in detail and which are relatively easy to compute. Furthermore, we find that for a given energy, the computational effort for GMRes can vary with dime...

[1]  D. Kouri,et al.  Direct calculation of the reactive transition matrix by L2 quantum mechanical variational methods with complex boundary conditions , 1989 .

[2]  David E. Manolopoulos,et al.  Quantum scattering via the log derivative version of the Kohn variational principle , 1988 .

[3]  B. Schneider,et al.  A direct iterative-variational method for solving large sets of linear algebraic equations , 1989 .

[4]  D. Kouri,et al.  Quantum dynamics of chemical reactions by converged algebraic variational calculations , 1990 .

[5]  J. H. Zhang.  Progress of basis optimization techniques in variational calculation of quantum reactive scattering , 1991 .

[6]  D. Truhlar,et al.  Improved techniques for outgoing wave variational principle calculations of converged state-to-state transition probabilities for chemical reactions , 1991 .

[7]  D. Kouri,et al.  L2 amplitude density method for multichannel inelastic and rearrangement collisions , 1988 .

[8]  R. Wyatt,et al.  Translational basis set contraction in variational reactive scattering , 1990 .

[9]  W. Lester,et al.  Theoretical study of inelastic scattering of H2 by Li+ on SCF and CI potential energy surfaces , 1975 .

[10]  R. Wyatt,et al.  THE NEWTON VARIATIONAL FUNCTIONAL FOR THE LOG-DERIVATIVE MATRIX : USE OF THE REFERENCE ENERGY GREEN'S FUNCTION IN AN EXCHANGE PROBLEM , 1990 .

[11]  L. Thomas.  Solution of the coupled equations of inelastic atom–molecule scattering for a single initial state. II. Use of nondiagonal matrix Green functions , 1982 .

[12]  D. Kouri,et al.  A comparative analysis of variational methods for inelastic and reactive scattering , 1990 .

[13]  B. C. Garrett,et al.  2?* Calculations of Accurate Quantal-dynamical Reactive Scattering Transition Probabilities and their Use to test Semiclassical Applications , 1987 .

[14]  R. Wyatt,et al.  Multilevel adaptive technique for quantum reactive scattering , 1987 .

[15]  M. Baer.  Arrangement channel approach to atom–diatom reactive systems: Theory and accurate three‐dimensional probabilities for the H+H2 system , 1989 .

[16]  O. Widlund.  A Lanczos Method for a Class of Nonsymmetric Systems of Linear Equations , 1978 .

[17]  R. Wyatt,et al.  H+H2(0,0)→H2(v’, j’)+H integral cross sections on the double many body expansion potential energy surface , 1990 .

[18]  Antonio Laganà,et al.  Supercomputer algorithms for reactivity, dynamics and kinetics of small molecules , 1989 .

[19]  William H. Miller,et al.  Quantum scattering via the S‐matrix version of the Kohn variational principle , 1988 .

[20]  R. Friedman,et al.  High-energy state-to-state quantum dynamics for D+H2 (v=j=1) → HD (v′=1, j′) + H , 1992 .

[21]  G. Golub,et al.  Iterative solution of linear systems , 1991, Acta Numerica.

[22]  D. Kouri,et al.  Variational basis-set calculations of accurate quantum mechanical reaction probabilities , 1987 .

[23]  Weitao Yang,et al.  Block Lanczos approach combined with matrix continued fraction for the S‐matrix Kohn variational principle in quantum scattering , 1989 .

[24]  Omar A. Sharafeddin,et al.  Spectroscopic analysis of transition state energy levels: Bending–rotational spectrum and lifetime analysis of H3 quasibound states , 1989 .

[25]  D. Kouri,et al.  ℒ2 Solution of the quantum mechanical reactive scattering problem. The threshold energy for D + H2(v = 1) → HD + H , 1986 .

[26]  D. Kouri,et al.  Scattered wave variational principle for atom—diatom reactive scattering: hybrid basis set calculations , 1991 .

[27]  B. C. Garrett,et al.  Semiclassical and Quantum Mechanical Calculations of Isotopic Kinetic Branching Ratios for the Reactionof O(3P) with HD , 1989 .

[28]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[29]  Donald G. Truhlar,et al.  Iterative methods for solving the non-sparse equations of quantum mechanical reactive scattering , 1989 .

[30]  A. Arthurs,et al.  The theory of scattering by a rigid rotator , 1960, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[31]  O. Crawford.  Calculation of Chemical Reaction Rates by the R‐Matrix Method , 1971 .

[32]  W. Miller,et al.  3D quantum scattering calculations of the reaction He + H+2 → HeH+ + H for total angular momentum J = 0 , 1990 .

[33]  D. Kouri,et al.  Quantum mechanical algebraic variational methods for inelastic and reactive molecular collisions , 1988 .

[34]  Donald G. Truhlar,et al.  A double many‐body expansion of the two lowest‐energy potential surfaces and nonadiabatic coupling for H3 , 1987 .

[35]  R. Wyatt,et al.  Calculations relating to the experimental observation of resonances in the H+H2 reaction , 1989 .

[36]  R. Wyatt,et al.  Lanczos recursion, continued fractions, Padé approximants, and variational principles in quantum scattering theory , 1988 .

[37]  W. Miller,et al.  New method for quantum reactive scattering, with applications to the 3-D H+H2 reaction , 1987 .

[38]  R. Glowinski,et al.  Computing Methods in Applied Sciences and Engineering , 1974 .

[39]  D. Kouri,et al.  Exact quantum dynamics and tests of the distorted-wave approximation for the O(3P)+ HD reaction , 1990 .

[40]  D. Truhlar,et al.  Variational reactive scattering calculations: computational optimization strategies , 1991 .

[41]  D. Truhlar,et al.  Benchmark calculations of thermal reaction rates. I. Quantal scattering theory , 1991 .

[42]  R. Wyatt,et al.  QUANTUM REACTIVE SCATTERING VIA THE LOG DERIVATIVE VERSION OF THE KOHN VARIATIONAL PRINCIPLE - GENERAL-THEORY FOR BIMOLECULAR CHEMICAL-REACTIONS , 1989 .

[43]  R. Wyatt,et al.  Recursive generation of individual S-matrix elements: Application to the collinear H + H2 reaction☆ , 1988 .