摘要
We present the Fortran package H2SOLV for an efficient computation of the nonrelativistic energy levels and the wave functions of diatomic two-electron molecules within the Born–Oppenheimer approximation. The wave function is obtained as a linear combination of the explicitly correlated exponential (Kołos–Wolniewicz) functions. The computations of H2SOLV are performed within the arbitrary-precision arithmetics, where the number of working digits can be adjusted by the user. The key part of H2SOLV is the implementation of the algorithm of an efficient computation of the two-center two-electron integrals for arbitrary values of internuclear distances developed by one of us (Pachucki, 2013). This have been one of the long-standing problems of quantum chemistry. The code is parallelized, suitable for large-scale computations limited only by the computer resources available and can produce highly accurate results. As an example, we report several benchmark results obtained with H2SOLV, including the energy value accurate to 18 decimal digits.Program summaryProgram title: H2SOLVCatalogue identifier: AFBA_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AFBA_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 8740No. of bytes in distributed program, including test data, etc.: 43089Distribution format: tar.gzProgramming language: Fortran 90.Computer: PCs or higher performance computers.Operating system: Linux / Unix.Has the code been vectorized or parallelized?: YesRAM: From several Mbytes to several Gbytes, depending on the size of the basis.Classification: 2.1.External routines: MPFR library [2] and GMP library [3] need to be pre-installed on the computer.Nature of problem:Numerical solution of the two-center two-electron Schrödinger equation within the Born–Oppenheimer approximation using the explicitly correlated basis set of Kołos–Wolniewicz functions.Solution method:The method of solution is based on the algorithm developed in Ref. [1].Running time:From seconds to days, depending on the size of the basis. The tests provided each take between 2 and 15 minutes to run.References:[1]K. Pachucki, Efficient approach to two-center exponential integrals with application to excited states of molecular hydrogen, Phys. Rev. A 88 (2013) 022507.[2]The GNU Multiple Precision Floating-Point Library, http://www.mpfr.org/.[3]The GNU Multiple Precision Arithmetic Library, https://gmplib.org/.