The XCHEM code was introduced in 2017 [1] to provide an accurate description of electron correlation and exchange in the electronic continuum of molecules at the same level as complete or restricted active-space self-consistent field (CASSCF or RASSCF) methods. This has allowed for an accurate description of molecular photoionization in the region of Feshbach resonances, shake up processes in which ionization is accompanied by excitation of one or several of the remaining electrons, and interchannel couplings. The success of XCHEM for small molecules has led us to improve its performance in several aspects, which now allows for the description of resonant molecular photoionization in larger systems. In addition, we have incorporated the possibility to calculate photoelectron angular distributions in the laboratory and molecular frames, which are essential to interpret angularly resolved photoionization experiments. Here we show its performance in the N2 and pyrazine molecules. The new version of the code, XCHEM-2.0, is freely available at https://gitlab.com/xchem/xchem_public. Program summaryProgram Title: XCHEM-2.0CPC Library link to program files:https://doi.org/10.17632/t8tbk9gdt2.1Developer's repository link:https://gitlab.com/xchem/xchem_public/-/archive/2.0.0/xchem_public-2.0.0.tar.gzLicensing provisions: LGPLProgramming language: Fortran, pythonNature of problem: Understanding photoionization of molecules is a recurrent goal in Atomic, Molecular and Optical physics, even more nowadays with the advent of XUV and X-ray high-harmonic generation sources and free electron lasers. Photoionization can be employed to induce electron dynamics in molecules, e.g., by using ultrashort XUV or X-ray pulses, or as a probe of these dynamics by recording time-resolved photoelectron spectra. At low photoelectron energies, photoionization may result from various competing processes, such as direct ionization, autoionization of Feshbach and shape resonances, shake up, etc., each of these leaving its trace in the measured spectra. So, interpretation of these spectra can be a difficult task and theoretical simulations are often required to separate the different contributions. However, a method able to do so must be able to describe electron correlation in the electronic continuum beyond the Hartree-Fock (HF) and configuration interaction singles (CIS) approximations. XCHEM-2.0 has been designed to account for all these processes by providing an accurate description of electron correlation in the ionization continuum of molecules.Solution method: XCHEM-2.0 employs a close coupling formalism to describe the electronic continuum of molecules at the level of multiference configuration interaction (MRCI) methods in combination with a hybrid Gaussian-Bspline basis set (the so-called GABs basis). In particular, it makes use of a restricted active space self-consistent field (RASSCF) approximation to evaluate the bound electronic states of the remaining molecular cation using localized Gaussian functions. These localized Gaussians are then supplemented with a single centered Gaussian expansion to allow the electron to leave the nuclear environment. Finally, the basis set includes a set of B-spline functions, centered at the same position as the single-centered Gaussians, which reach the asymptotic region. Thanks to the last B-spline function, which does not vanish at the box boundary, XCHEM-2.0 can determine continuum states fulfilling any required asymptotic behavior. The present version of the code improves on the original version by incorporating a more efficient augmentation procedure to build Ne-electron configurations from (Ne−1) ones, reducing the space to store data, providing a more efficient removal of linear dependencies resulting from the over-completeness of the polycentric+GABS basis and a more efficient solution of the scattering equations, and including new routines to calculate photoelectron angular distributions in the laboratory and molecular frames.Additional comments including restrictions and unusual features: The proposed method focuses on RASSCF methodology, so part of electron correlation is still not included. Also, one- and two electron integrals involving exclusively Gaussian functions are done using a modified version of the OpenMOLCAS software, so the installation of this code is also required and limitation associated with the latter is applicable. The code is limited to include a few ionization channels (around 20) with medium angular momentum (limited to ℓ = 15 as in OpenMOLCAS).
Read full abstract