SAC-CI general- R study of the ionization spectrum of HCl

Reference Energies for Valence Ionizations and Satellite Transitions

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

Fine Theoretical Spectroscopy Using SAC-CI General-R Method: Outer- and Inner-Valence Ionization Spectra of N2O and HN3

Fine theoretical spectroscopy for the outer- and inner-valence ionization spectra has been presented by using the SAC-CI (symmetry adapted cluster-configuration interaction) general-R method applied to N2O and HN3. The SAC-CI general-R method accurately simulated the experimental spectra and the detailed assignments of the satellite peaks were proposed. The continuous peaks I ∼ V of N2O observed by the dipole (e,2e) spectroscopy were finely reproduced. In particular, the low-lying satellites were calculated in good agreement with the experiment; two 2Π shake-up states at the foot of the C state, 2Σ and 2Π states for the peak I. For HN3, new interpretation was proposed for the outer-valence region, namely, peaks 3 ∼ 5 are composed of the mixture of the single-electron main peaks and the two-electron shake-up peaks. Inner-valence ionization spectrum of HN3 was theoretically predicted for which no experimental spectrum has been reported.

Symmetry-adapted-cluster configuration interaction study of the doublet states of HCl+: Potential energy curves, dipole moments, and transition dipole moments

The electronic structure of the HCl(+) molecular ion has been calculated using the general-R symmetry-adapted-cluster configuration interaction (SAC-CI) method. The authors present the potential energy curves, dipole moments, and transition dipole moments for a series of doublet states. The data are compared with the previous CASSCF and MCSCF calculations. The SAC-CI results reproduce quite well the data available in literature and extend the knowledge on the HCl(+) electronic structure for several higher states. The calculated R-dependent behavior of both dipole moments and transition dipole moments for a series of bound and unbound states reveals an intricate dissociation process at intermediate distances (R>R(e)). The pronounced maxima in transition dipole moment (TDM) describing transitions into high electronic states (X (2)Pi-->3 (2)Pi, X (2)Pi-->3 (2)Sigma, 2 (2)Pi-->3 (2)Pi, 3 (2)Pi-->4 (2)Pi) occur at different interatomic separations. Such TDM features are promising for selection of excitation pathways and, consequently, for an optimal control of the dissociation products.

Explicitly intruder-free valence-universal multireference coupled cluster theory as applied to ionization spectroscopy

Although it is quite promising to compute the spectroscopic energies [say, ionization potential (IP)] via the traditional valence-universal multireference coupled cluster (VUMRCC) method based on the description of the complete model space being seriously plagued by the perennial intruder state problem, the eigenvalue independent partitioning (EIP) based VUMRCC (coined as EIP-MRCC) method is quite effective to predict the spectroscopic energies in an intruder-free manner. Hence, the EIP-MRCC method is suitable for generating both the principal IPs and the satellite IPs of the inner-valence region. An EIP strategy converts the nonlinear VUMRCC equations for M(m,n) dimensional model space of m hole and n particle to a non-Hermitian eigenproblem of larger dimension whose M(m,n) roots are only physically meaningful. To increase the quality of the computed energy differences in the sense of chemical accuracy and to locate the correct position of it in the spectrum, the inclusion of higher-body cluster operators on top of all the standard singles-doubles is not the only pivotal issue, the effect of the size of the basis set is also equally important. This paper illustrates these issues by calculating the principal and satellite IPs of HF and HCl molecules using various basis sets (viz., Dunning's cc-pVDZ, cc-pVTZ, and cc-pVQZ) via EIP-MRCC method with full inclusion of triples (abbreviated as EIP-MRCCSDT). The results seem quite encouraging in comparison with the experimental values. The controversial 2Pi satellite at 28.67 eV of HCl of Svensson et al. [J. Chem. Phys. 89, 7193 (1988)] is also reported.

Analytical energy gradient of the symmetry-adapted-cluster configuration-interaction general-R method for singlet to septet ground and excited states

A method of calculating analytical energy gradients of the singlet and triplet excited states, ionized states, electron-attached states, and high-spin states from quartet to septet states by the symmetry-adapted-cluster configuration-interaction general-R method is developed and implemented. This method is a powerful tool in the studies of geometries, dynamics, and properties of the states of molecules in which not only one-electron processes but also two- and multielectron processes are involved. The performance of the present method was confirmed by calculating the geometries and the spectroscopic constants of the diatomic and polyatomic molecules in various electronic states involving the ground state and the one- to three-electron excited states. The accurate descriptions were obtained for the equilibrium geometries, vibrational frequencies, and adiabatic excitation energies, which show the potential usefulness of the present method. The particularly interesting applications were to the C' 1Ag state of acetylene, the A 2Deltau and B 2Sigmau+ states of CNC and the 4B1 and a 4Piu states of N3 radical.

Theoretical fine spectroscopy with symmetry-adapted-cluster configuration-interaction method: Outer- and inner-valence ionization spectra of furan, pyrrole, and thiophene

Theoretical fine spectroscopy has been performed for the valence ionization spectra of furan, pyrrole, and thiophene with the symmetry-adapted-cluster configuration-interaction general-R method. The present method described that the pi(1) state interacts with the pi(3) (-2)pi*, pi(2) (-2)pi*, and pi(2) (-1)pi(3) (-1)pi* shake-up states providing the split peaks and the outer-valence satellites, both of which are in agreement with the experiments. The intensity distributions were analyzed in detail for the inner-valence region. In particular, for furan, theoretical intensities were successfully compared with the intensity measured by the electron momentum spectroscopy. The interactions of the 3b(2) and 5a(1) states with the shake-up states were remarkable for furan and pyrrole, while the 4b(2) state of thiophene had relatively large intensity.