Accurate ab initio potential energy curve of F2. II. Core-valence correlations, relativistic contributions, and long-range interactions

Coupled Electron Pair-Type Approximations for Tensor Product State Wave Functions

Size extensivity, defined as the correct scaling of energy with system size, is a desirable property for any many-body method. Traditional configuration interaction (CI) methods are not size extensive, hence the error increases as the system gets larger. Coupled electron pair approximation (CEPA) methods can be constructed as simple extensions of a truncated CI that ensures size extensivity. One of the major issues with the CEPA and its variants is that singularities arise in the amplitude equations when the system starts to be strongly correlated. In this work, we extend the traditional Slater determinant based coupled electron pair approaches like CEPA-0, averaged coupled-pair functional, and average quadratic coupled-cluster to a new formulation based on tensor product states (TPS). We show that a TPS basis can often be chosen such that it removes the singularities that commonly destroy the accuracy of CEPA based methods. A suitable TPS representation can be formed by partitioning the system into separate disjoint clusters and forming the final wave function as the tensor product of the many body states of these clusters. We demonstrate the application of these methods on simple bond breaking systems such as CH₄ and F₂ where determinant based CEPA methods fail. We further apply the TPS-CEPA approach to stillbene isomerization and few planar π-conjugated systems. Overall, the results show that the TPS-CEPA method can remove the singularities and provide improved numerical results compared to common electronic structure methods.

Relativistic four-component MRCISD+Q calculations of the six lowest valence states of molecular F[Formula: see text] anion including Breit interactions

In this study, the potential energy curves of the ground and the excited states of molecular fluorine anion (F[Formula: see text]) were investigated at multireference configuration interaction (MRCISD) with Davidson size-extensivity correction (denoted as +Q) within fully relativistic four-component relativistic framework including Breit interaction. Spectroscopic constants (R_e, ω_e, ω_ex_e, ω_ey_e, D_e,D₀,B_e, α_e, β_e, γ_e ), accurate extended Rydberg analytical form and rovibrational levels for ground state X:[Formula: see text] are presented, as well as spectroscopic constants for non dissociative excited states. For most states these spectroscopic constants are presented for the first time in literature and they are of interest for experimental studies, specially regarding electron attachment of F₂. Results suggest that inclusion of relativistic effects at 4-component level and correlation effects treated at MRCISD+Q level are needed to obtain reliable results, which we report for X:[Formula: see text] ground state's R_e, ω_e and D_e the values of 1.999 Å, 391 cm^- 1 and 1.22 eV, respectively.

Multiconfiguration Density-Coherence Functional Theory

This paper presents a new theory called multiconfiguration density coherence functional theory (MC-DCFT). This theory provides a new route to define density functionals for multiconfiguration wave functions, in particular by using the one-particle density matrix in the coordinate representation. The theory is illustrated by calculating the dissociation curve of four heteronuclear and homonuclear diatomic molecules, namely, H₂, F₂, N₂, and HF, using density coherence functionals converted from PBE, BLYP, and PW91. By introducing two parameters in the converted density functionals, we are able to calculate bond dissociation energies of comparable accuracy as those calculated by multiconfiguration pair-density functional theory (MC-PDFT) and complete active space second-order perturbation theory (CASPT2). This demonstrates that it would be possible to build a successful multiconfiguration density functional theory based on density coherence.

Accurate ground state potential of Cu2 up to the dissociation limit by perturbation assisted double-resonant four-wave mixing

Perturbation facilitated double-resonant four-wave mixing is applied to access high-lying vibrational levels of the X ¹Σ_g ⁺ (0_g ⁺) ground state of Cu₂. Rotationally resolved transitions up to v^″ = 102 are measured. The highest observed level is at 98% of the dissociation energy. The range and accuracy of previous measurements are significantly extended. By applying the near dissociation equation developed by Le Roy [R. J. Le Roy, J. Quant. Spectrosc. Radiat. Transfer 186, 197 (2017)], a dissociation energy of D_e = 16 270(7) hc cm^-1 is determined, and an accurate potential energy function for the X ¹Σ_g ⁺ (0_g ⁺) ground state is obtained. Molecular constants are determined from the measured transitions and by solving the radial Schrödinger equation using this function and are compared with results from earlier measurements. In addition, benchmark multi-reference configuration interaction computations are performed using the Douglas-Kroll-Hess Hamiltonian and the appropriate basis of augmented valence quadruple ζ type. Coupled-cluster single, double, and perturbative triple calculations were performed for comparison.

Reproducing benchmark potential energy curves of molecular bond dissociation with small complete active space aided with density and density-matrix functional corrections

Various effects of electron correlation accompany molecular bond dissociation, which makes the efficient calculation of potential energy curves a notoriously difficult problem. In an attempt to reliably reproduce both absolute energies and shapes of the benchmark dissociation curves, calculations with the combined CASΠDFT method are carried out for the prototype molecules H₂, BH, F₂, and N₂. The complete active space (CAS) part of CASΠDFT accounts for long-range nondynamic correlation, while short-range dynamic correlation is accounted for with the corrected Lee-Yang-Parr correlation functional of density functional theory (DFT). The correction represents the suppression of dynamic correlation with nondynamic correlation, and it is a function of the ratio x(r) between the conditional and conventional densities obtained with the CAS on-top pair density Π(r). For the single-bonded molecules H₂, BH, and F₂, CASΠDFT succeeds in reproducing the shapes and absolute energies (for H₂ and BH) of the benchmark curves, while for the triple-bonded N₂ molecule, the addition to CASΠDFT of a multibond correction is required. It accounts for the middle-range dynamic correlation of the same-spin electrons in the (symmetrized) high-spin atomic electron configurations of the dissociating N₂.

Benchmark of correlation matrix renormalization method in molecule calculations

We report benchmark calculations of the correlation matrix renormalization (CMR) approach for 23 molecules in the well-established G2 molecule set. This subset represents molecules with spin-singlet ground state in a variety of chemical bonding and coordination environments. The QUAsi-atomic minimal basis-set orbitals (QUAMBOs) are used as local orbitals in both CMR and full configuration interaction (FCI) calculations for comparison. The results obtained from the calculations are also compared with available experimental data. It is shown that the CMR method produces binding and dissociation energy curves in good agreement with the QUAMBO-FCI calculations as well as experimental results. The CMR benchmark calculations yield a standard deviation of 0.09 Å for the equilibrium bond length and 0.018 Hartree/atom for the formation energy, with a gain of great computational efficiency which scales like Hartree-Fock method.

Term rules for simple metal clusters

Hund’s term rules are only valid for isolated atoms, but have no generalization for molecules or clusters of several atoms. We present a benchmark calculation of Al₂ and Al₃, for which we find the high and low-spin ground states ³Π_u and , respectively. We show that the relative stabilities of all the molecular terms of Al₂ and Al₃ can be described by simple rules pertaining to bonding structures and symmetries, which serve as guiding principles to determine ground state terms of arbitrary multi-atom clusters.

Relativistic state-specific multireference perturbation theory incorporating improved virtual orbitals: Application to the ground state single-bond dissociation

Using four-component (4c) relativistic spinors, we present a computationally economical relativistic ab initio method for molecular systems employing our recently proposed second-order state-specific multireference perturbation theory (SSMRPT) incorporating the improved virtual orbital-complete active space configuration interaction (IVO-CASCI) reference wavefunction. The resulting method, 4c-IVO-SSMRPT [calculate one state at a time] is tested in pilot calculations on the homonuclear dimers including Li(2), Na(2), K(2), Rb(2), F(2), Cl(2), and Br(2) through the computations of the ground state potential energy curves (PECs). As SSMRPT curbs intruder effects, 4c-IVO-SSMRPT is numerically stable. To our knowledge, the SSMRPT in the 4c relativistic framework has not been explored in the past. Selective spectroscopic constants that are closely related to the correct shape and accuracy of the energy surfaces have been extracted from the computed PECs. For the halogen molecules, a relativistic destabilization of the bond has been found. Relativistic and electron correlation effects need to be incorporated to get reliable estimates. Our results are in good accordance with reference theoretical and experimental data which manifests the computational accuracy and efficiency of the new 4c-IVO-SSMRPT method. The method opens for an improved description of MR systems containing heavy elements. The inexpensiveness of IVO-CASCI makes 4c-IVO-SSMRPT method promising for studies on large systems of heavy elements.

Intramolecular halogen-halogen bonds?

By analysing the properties of the electron density in the structurally simple perhalogenated ethanes, X3C-CY3 (X, Y = F, Cl), a previously overlooked non-covalent attraction between halogens attached to opposite carbon atoms is found. Quantum chemical calculations extrapolated towards the full solution of the Schrödinger equation reveal the complex nature of the interaction. When at least one of the halogens is a chlorine, the strength of the interaction is comparable to that of hydrogen bonds. Further analysis shows that the bond character is quite different from standard non-covalent halogen bonds and hydrogen bonds; no bond critical points are found between the halogens, and the σ-holes of the halogens are not utilised for bonding. Thus, the nature of the intramolecular halogen···halogen bonding studied here appears to be of an unusually strong van der Waals type.

A companion perturbation theory for state-specific multireference coupled cluster methods

A partitioning scheme is applied to the state-specific Mukherjee multireference coupled cluster method to derive a companion perturbation theory (Mk-MRPT2). A production-level implementation of Mk-MRPT2 is reported. The effectiveness of the Mk-MRPT2 method is demonstrated by application to the classic F(2) dissociation problem and the lowest-lying electronic states of meta-benzyne, including computations with up to 766 atomic orbitals. We show that Mk-MRPT2 theory is particularly useful in multireference focal point extrapolations to determine ab initio limits.

Multipole polarizabilities and long-range interactions of the fluorine atom

The polarizabilities of the ground state of fluorine are determined by the box-based B-spline configuration-interaction method. The ground-state static and tensor polarizabilities are calculated as alpha(d) = 3.49a(0)(3) and alpha(2) = 0.303a(0)(3), respectively. The resulting excitation matrix elements are employed to compute the interatomic dispersion interactions for the HF and F(2) molecules.

Accurate ab initio potential energy curve of F2. III. The vibration rotation spectrum

An analytical expression is found for the accurate ab initio potential energy curve of the fluorine molecule that has been determined in the preceding two papers. With it, the vibrational and rotational energy levels of F(2) are calculated using the discrete variable representation. The comparison of this theoretical spectrum with the experimental spectrum, which had been measured earlier using high-resolution electronic spectroscopy, yields a mean absolute deviation of about 5 cm(-1) over the 22 levels. The dissociation energy with respect to the lowest vibrational energy is calculated within 30 cm(-1) of the experimental value of 12 953+/-8 cm(-1). The reported agreement of the theoretical spectrum and dissociation energy with experiment is contingent upon the inclusion of the effects of core-generated electron correlation, spin-orbit coupling, and scalar relativity. The Dunham analysis [Phys. Rev. 41, 721 (1932)] of the spectrum is found to be very accurate. New values are given for the spectroscopic constants.