G PU A ccelerated IN tensities MPI (GAIN-MPI): A new method of computing Einstein-A coefficients

Abstract Calculating dipole transition intensities or the related Einstein A coefficients can dominate the computer usage for large line lists of transitions such as those being computed to model radiative transport through hot atmospheres. An algorithm for the efficient computation of line strengths is presented based on the use of the half-linestrength. This is implemented on GPUs that are shown to give up to a thousandfold speed-up compared to calculations on conventional computers. This algorithm is implemented in the program GAIN which was developed as part of the TROVE nuclear motion program, but can be adapted for use by other similar programs in a straightforward fashion. Program summary Program title: GAIN-MPI Program Files doi: http://dx.doi.org/10.17632/4x75jsphc6.1 Licensing provisions: MIT licence. Programming language: C++ 99, CUDA C and Fortran 95. Nature of problem: Computation of linestrengths using GPU hardware Solution method : Split the linestrength into smaller blocks and compute them in the GPU in parallel Restrictions: The current version is restricted to separable rovibrational basis sets Unusual features: Can be extended by user supplied concrete C++ classes for the MPI version

[1]  M. Mladenović.  Discrete variable approaches to tetratomic molecules: part I: DVR(6) and DVR(3) + DGB methods. , 2002, Spectrochimica acta. Part A, Molecular and biomolecular spectroscopy.

[2]  T. Carrington,et al.  Discrete‐Variable Representations and their Utilization , 2007 .

[3]  J. Tennyson,et al.  A global, high accuracy ab initio dipole moment surface for the electronic ground state of the water molecule. , 2011, The Journal of chemical physics.

[4]  S. Tashkun,et al.  A room temperature CO2 line list with ab initio computed intensities , 2016, 1601.05334.

[5]  Ralf Schneider,et al.  Ab initio modeling of molecular IR spectra of astrophysical interest: Application to CH4 , 2009 .

[6]  J. Muckerman,et al.  A General Variational Algorithm to Calculate Vibrational Energy Levels of Tetraatomic Molecules , 2002 .

[7]  Andr'e Nauts,et al.  Exact numerical computation of a kinetic energy operator in curvilinear coordinates , 2002 .

[8]  Anthony Papagiannis.  Intern , 2010, BMJ : British Medical Journal.

[9]  V. Szalay,et al.  Variational vibrational calculations using high-order anharmonic force fields , 2004 .

[10]  T. Carrington,et al.  Computing rovibrational levels of methane with curvilinear internal vibrational coordinates and an Eckart frame. , 2013, The Journal of chemical physics.

[11]  Joel M. Bowman,et al.  The adiabatic rotation approximation for rovibrational energies of many-mode systems: Description and tests of the method , 1998 .

[12]  J. Tennyson,et al.  ExoMol line lists – VII. The rotation–vibration spectrum of phosphine up to 1500 K , 2014, 1410.2917.

[13]  T. Carrington,et al.  A multi-dimensional Smolyak collocation method in curvilinear coordinates for computing vibrational spectra. , 2015, The Journal of chemical physics.

[14]  J. Hodges,et al.  High-Accuracy CO(2) Line Intensities Determined from Theory and Experiment. , 2015, Physical review letters.

[15]  J. Tennyson,et al.  A variationally computed line list for hot NH3 , 2010, 1011.1569.

[16]  T. Carrington,et al.  A discrete variable representation method for studying the rovibrational quantum dynamics of molecules with more than three atoms. , 2009, The Journal of chemical physics.

[17]  J. Tennyson,et al.  Hybrid variational–perturbation method for calculating ro-vibrational energy levels of polyatomic molecules , 2014, 1411.6098.

[18]  Hua Guo.  Recursive Solutions to Large Eigenproblems in Molecular Spectroscopy and Reaction Dynamics , 2007 .

[19]  F. Gatti,et al.  Exact and constrained kinetic energy operators for polyatomic molecules: The polyspherical approach , 2009 .

[20]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[21]  J. Tennyson,et al.  The ab initio calculation of the vibrational‐rotational spectrum of triatomic systems in the close‐coupling approach, with KCN and H2Ne as examples , 1982 .

[22]  S. Yurchenko,et al.  Automatic differentiation method for numerical construction of the rotational-vibrational Hamiltonian as a power series in the curvilinear internal coordinates using the Eckart frame. , 2015, The Journal of chemical physics.

[23]  B. M. Fulk.  MATH , 1992 .

[24]  Joel M Bowman,et al.  Calculations of rovibrational energies and dipole transition intensities for polyatomic molecules using MULTIMODE. , 2009, The Journal of chemical physics.

[25]  M. Mladenović.  Rovibrational Hamiltonians for general polyatomic molecules in spherical polar parametrization. I. Orthogonal representations , 2000 .

[26]  P. Jensen,et al.  Theoretical ROVibrational Energies (TROVE) : A robust numerical approach to the calculation of rovibrational energies for polyatomic molecules , 2007 .

[27]  P. Jensen,et al.  Rotation-vibration motion of pyramidal XY₃ molecules described in the Eckart frame: The calculation of intensities with application to NH₃ , 2005 .

[28]  J. Tennyson,et al.  Stark coefficients for highly excited rovibrational states of H2O. , 2012, The Journal of chemical physics.

[29]  M. Mladenović.  Discrete variable approaches to tetratomic molecules: part II: application to H2O2 and H2CO. , 2002, Spectrochimica acta. Part A, Molecular and biomolecular spectroscopy.

[30]  P. Jensen.  A new morse oscillator-rigid bender internal dynamics (MORBID) Hamiltonian for triatomic molecules , 1988 .

[31]  A. Császár,et al.  Toward black-box-type full- and reduced-dimensional variational (ro)vibrational computations. , 2009, The Journal of chemical physics.

[32]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[33]  A. Nikitin,et al.  An efficient method for energy levels calculation using full symmetry and exact kinetic energy operator: tetrahedral molecules. , 2015, The Journal of chemical physics.

[34]  Jonathan Tennyson,et al.  New vibration-rotation code for tetraatomic molecules exhibiting wide-amplitude motion: WAVR4 , 2004, Comput. Phys. Commun..

[35]  J. Tennyson,et al.  A hybrid line list for CH4 and hot methane continuum. , 2017, Astronomy and astrophysics.

[36]  T. Carrington,et al.  Using a nondirect product discrete variable representation for angular coordinates to compute vibrational levels of polyatomic molecules. , 2008, The Journal of chemical physics.

[37]  T. Carrington,et al.  CALCULATING ROVIBRATIONAL ENERGY LEVELS OF A TRIATOMIC MOLECULE WITH A SIMPLE LANCZOS METHOD , 1999 .

[38]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[39]  J. Tennyson,et al.  Communication: Visible line intensities of the triatomic hydrogen ion from experiment and theory. , 2014, The Journal of chemical physics.

[40]  P. Bunker,et al.  Molecular symmetry and spectroscopy , 1979 .

[41]  J. Tennyson,et al.  Spectrum of hot methane in astronomical objects using a comprehensive computed line list , 2014, Proceedings of the National Academy of Sciences.

[42]  Jonathan Tennyson,et al.  DVR3D: a program suite for the calculation of rotation-vibration spectra of triatomic molecules , 2004, Comput. Phys. Commun..

[43]  J. Tennyson,et al.  A variationally calculated room temperature line-list for H2O2 , 2015 .

[44]  H. Meyer,et al.  Rovibrational spectroscopy using a kinetic energy operator in Eckart frame and the multi-configuration time-dependent Hartree (MCTDH) approach. , 2014, The Journal of chemical physics.

[45]  F. Gatti,et al.  Exact and constrained kinetic energy operators for polyatomic molecules: The polyspherical approach , 2008 .

[46]  D. Schwenke.  Variational Calculations of Rovibrational Energy Levels and Transition Intensities for Tetratomic Molecules , 1996 .

[47]  T. Carrington,et al.  Discrete variable representations of complicated kinetic energy operators , 1994 .

[48]  N. Handy,et al.  Extensions and tests of “multimode”: a code to obtain accurate vibration/rotation energies of many-mode molecules , 1998 .

[49]  S. Carter,et al.  MULTIMODE: A code to calculate rovibrational energies of polyatomic molecules , 2003 .

[50]  J. Tennyson,et al.  ExoMol line lists – XV. A new hot line list for hydrogen peroxide , 2016, 1607.00498.

[51]  T. Carrington.  Two new methods for computing vibrational energy levels , 2015 .

[52]  Jonathan Tennyson,et al.  The ExoMol project: Software for computing large molecular line lists , 2016, 1607.01220.

[53]  S. Carter,et al.  Large vibrational variational calculations using ‘multimode’ and an iterative diagonalization technique , 2004 .

[54]  Sergei N. Yurchenko,et al.  ExoMol line lists IV: The rotation-vibration spectrum of methane up to 1500 K , 2014, 1401.4852.

[56]  J. Bloino.  A VPT2 route to near-infrared spectroscopy: the role of mechanical and electrical anharmonicity. , 2015, The journal of physical chemistry. A.

[57]  J. Tennyson,et al.  ExoMol line lists – VIII. A variationally computed line list for hot formaldehyde , 2015, 1506.00172.

[58]  J. Tennyson,et al.  A hybrid variational-perturbation calculation of the ro-vibrational spectrum of nitric acid. , 2014, The Journal of chemical physics.

[59]  A. Nikitin,et al.  THEORETICAL HOT METHANE LINE LISTS UP TO T = 2000 K FOR ASTROPHYSICAL APPLICATIONS , 2014 .