Higher-order Electron Correlation Research Articles

Simulating transient X-ray photoelectron spectra of Fe(CO)5 and its photodissociation products with multireference algebraic diagrammatic construction theory.

Accurate simulations of transient X-ray photoelectron spectra (XPS) provide unique opportunities to bridge the gap between theory and experiment in understanding the photoactivated dynamics in molecules and materials. However, simulating X-ray photoelectron spectra along a photochemical reaction pathway is challenging as it requires accurate description of electronic structure incorporating core-hole screening, orbital relaxation, electron correlation, and spin-orbit coupling in excited states or at nonequilibrium ground-state geometries. In this work, we employ the recently developed multireference algebraic diagrammatic construction theory (MR-ADC) to investigate the core-ionized states and X-ray photoelectron spectra of Fe(CO)5 and its photodissociation products (Fe(CO)4, Fe(CO)3) following excitation with 266 nm light. The simulated transient Fe 3p and CO 3σ XPS spectra incorporating spin-orbit coupling and high-order electron correlation effects are shown to be in a good agreement with the experimental measurements by Leitner et al. [J. Chem. Phys., 2018, 149, 044307]. Our calculations suggest that core-hole screening, spin-orbit coupling, and ligand-field splitting effects are similarly important in reproducing the experimentally observed chemical shifts in transient Fe 3p XPS spectra of iron carbonyl complexes. Our results also demonstrate that the MR-ADC methods can be very useful in interpreting the transient XPS spectra of transition metal compounds.

Ab initio potential energy surface and vibration-rotation energy levels of magnesium monohydroxide revisited

The accurate potential energy surface of magnesium monohydroxide, MgOH, in its ground electronic state X˜2Σ+ has been determined from ab initio calculations using the coupled-cluster approach in conjunction with the correlation-consistent basis sets up to septuple-zeta quality. The core-electron correlation, higher-order electron correlation, scalar relativistic, and adiabatic effects were taken into account. The equilibrium configuration of the MgOH molecule was confirmed to be linear, although with the bending potential energy function being very flat near the minimum. The vibration–rotation-spin energy levels of the MgOH, MgOD, 25MgOH, and 26MgOH isotopologues were predicted using a variational approach. The spectroscopic constants of these isotopologues were determined to high accuracy.

Ab initio calculations of the 2p3/2→2s transition in He-, Li-, and Be-like uranium

The bound-state QED approach is applied to calculations of the $2p_{3/2} \rightarrow 2s$ transition energies in He-, Li-, and Be-like uranium. For U$^{90+}$ and U$^{89+}$, standard perturbation theory for a single level is employed, while the calculations of U$^{88+}$ have required its counterpart for quasidegenerate levels. The utilized approach merges the rigorous QED treatment up to the second order of perturbation theory with the higher-order electron-correlation contributions evaluated within the Breit approximation. The higher-order screened QED effects are estimated by means of the model-QED operator. The nuclear recoil, nuclear polarization, and nuclear deformation effects are taken into account as well. Along with the transition energies, their pairwise differences are calculated. The comprehensive analysis of the uncertainties due to uncalculated effects is carried out, and the most accurate theoretical predictions, which are in perfect agreement with available experimental data, are obtained.

Complex-scaled ab initio QED approach to autoionizing states

Ab initio method based on a complex-scaling approach and aimed at a rigorous QED description of autoionizing states is worked out. The autoionizing-state binding energies are treated nonperturbatively in $\alpha Z$ and include all the many-electron QED contributions up to the second order. The higher-order electron correlation, nuclear recoil, and nuclear polarization effects are taken into account as well. The developed formalism is demonstrated on the $LL$ resonances in heliumlike argon and uranium. The most accurate theoretical predictions for the binding energies are obtained.

Reconciliation of Theoretical Lifetimes of the 5s5p Metastable State for 88Sr with Measurement: The Role of the Blackbody-Radiation-Induced Decay

We conducted measurement and calculation to resolve the long-standing large discrepancy in the metastable state lifetime for the 88Sr atom between theoretical and experimental results. The present lifetime s, measured using the magneto-optical trap as a photon amplifier to detect the weak decay events, is approximately 60% larger than the previous experimental value s. By considering the electron correlation effects in the framework of the multiconfiguration Dirac–Hartree–Fock theory, we obtained a theoretical lifetime of 1079(54) s, which lies in the range of measurements with error bars. Furthermore, we considered the higher-order electron correlation and Breit interaction to control the uncertainty of the theoretical calculation. The significant improvement in the agreement between calculations and measurements is attributed to the updated blackbody radiation-induced decay rate.

Signatures of moiré trions in WSe2/MoSe2 heterobilayers.

Moiré superlattices formed by van der Waals materials can support a wide range of electronic phases, including Mott insulators1-4, superconductors5-10 and generalized Wigner crystals2. When excitons are confined by a moiré superlattice, a new class of exciton emerges, which holds promise for realizing artificial excitonic crystals and quantum optical effects11-16. When such moiré excitons are coupled to charge carriers, correlated states may arise. However, no experimental evidence exists for charge-coupled moiré exciton states, nor have their properties been predicted by theory. Here we report the optical signatures of trions coupled to the moiré potential in tungsten diselenide/molybdenum diselenide heterobilayers. The moiré trions show multiple sharp emission lines with a complex charge-density dependence, in stark contrast to the behaviour of conventional trions. We infer distinct contributions to the trion emission from radiative decay in which the remaining carrier resides in different moiré minibands. Variation of the trion features is observed in different devices and sample areas, indicating high sensitivity to sample inhomogeneity and variability. The observation of these trion features motivates further theoretical and experimental studies of higher-order electron correlation effects in moiré superlattices.

Approaching meV level for transition energies in the radium monofluoride molecule RaF and radium cation Ra+ by including quantum-electrodynamics effects.

Highly accurate theoretical predictions of transition energies in the radium monofluoride molecule, 226RaF, and radium cation, 226Ra+, are reported. The considered transition X2Σ1/2 → A2Π1/2 in RaF is one of the main features of this molecule and can be used to laser-cool RaF for a subsequent measurement of the electron electric dipole moment. For molecular and atomic predictions, we go beyond the Dirac-Coulomb Hamiltonian and treat high-order electron correlation effects within the coupled cluster theory with the inclusion of quadruple and ever higher amplitudes. The effects of quantum electrodynamics (QED) are included non-perturbatively using the model QED operator that is now implemented for molecules. It is shown that the inclusion of the QED effects in molecular and atomic calculations is a key ingredient in resolving the discrepancy between the theoretical values obtained within the Dirac-Coulomb-Breit Hamiltonian and the experiment. The remaining deviation from the experimental values is within a few meV. This is more than an order of magnitude better than the "chemical accuracy," 1 kcal/mol = 43 meV, that is usually considered as a guiding thread in theoretical molecular physics.

Abinitio Calculations of Energy Levels in Be-Like Xenon: Strong Interference between Electron-Correlation and QED Effects.

The strong mixing of close levels with two valence electrons in Be-like xenon greatly complicates abinitio QED calculations beyond the first-order approximation. Because of a strong interplay between the electron-electron correlation and QED effects, the standard single-level perturbative QED approach may fail, even if it takes into account the second-order screened QED diagrams. In the present Letter, the corresponding obstacles are overcome by working out the QED perturbation theory for quasidegenerate states. The contributions of all the Feynman diagrams up to the second order are taken into account. The many-electron QED effects are rigorously evaluated in the framework of the extended Furry picture to all orders in the nuclear-strength parameter αZ. The higher-order electron-correlation effects are considered within the Breit approximation. The nuclear recoil effect is accounted for as well. The developed approach is applied to high-precision QED calculations of the ground and singly excited energy levels in Be-like xenon. The most accurate theoretical predictions for the binding and excitation energies are obtained. These results deviate from the most precise experimental value by 3σ but perfectly agree with a more recent measurement.

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions.

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree-Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

QED theory of the normal mass shift in few-electron atoms

The electron-electron interaction correction of first order in $1/Z$ to the one-electron part of the nuclear recoil effect on binding energies in atoms and ions is considered within the framework of the rigorous QED approach. The calculations to all orders in $\alpha Z$ are performed for the $1s^2$ state in heliumlike ions and the $1s^2 2s$ and $1s^2 2p_{1/2}$ states in lithiumlike ions in the range $Z=5$--$100$. The results obtained are compared with the Breit-approximation values. The performed calculations complete a systematic treatment of the QED nuclear recoil effect up to the first order in $1/Z$. The correction obtained is combined with the previously studied two-electron part as well as the higher-order electron-correlation corrections evaluated within the Breit approximation to get the total theoretical predictions for the mass shifts.

Stochastic many-body perturbation theory for electron correlation energies.

Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone's time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone's time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.

Beyond chemical accuracy in the heavy p-block: The first ionization potentials and electron affinities of Ga-Kr, In-Xe, and Tl-Rn.

A relativistic coupled-cluster version of the Feller-Peterson-Dixon composite method has been used to accurately calculate the first ionization potentials (IPs) and electron affinities (EAs) of the post-d, p-block elements Ga-Rn. Complete basis set extrapolations including outer-core correlation at the CCSD(T) level of theory were combined with contributions from higher order electron correlation up to CCSDTQ, quantum electrodynamic effects (Lamb shift), and spin-orbit (SO) coupling including the Gaunt contribution. Several methods for including SO were investigated, in which all involved the four-component (4c) Dirac-Coulomb (DC) Hamiltonian. The treatment of SO coupling was the contribution that limited the final accuracy of the present results. In the cases where 4c-DC-CCSD(T) could be reliably used for the SO contributions, the final composite IPs and EAs agreed with the available experimental values to within an unsigned average error of just 0.16 and 0.20 kcal/mol, respectively. In all cases, the final IPs and EAs were within 1 kcal/mol of the available experimental values, except for the EAs of the group 13 elements (Ga, In, and Tl), where the currently accepted experimental values appear to be too large by as much as 4 kcal/mol. The values predicted in this work, which have estimated uncertainties of ±0.5 kcal/mol, are 5.25 (Ga), 7.69 (In), and 7.39 (Tl) kcal/mol. For the EAs of Po and At, which do not have experimental values, the current calculations predict values of 34.2 and 55.8 kcal/mol with estimated uncertainties of ±0.6 and ±0.3 kcal/mol, respectively.

On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions.

In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C 6 dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H2O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α iso) and nearly 50% to the anisotropic (α aniso) polarizability values. In doing so, we find that α iso = 9.8307 bohr3, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.

RI-MP3 calculations of biomolecules based on the fragment molecular orbital method.

In this study, the third-order Møller-Plesset perturbation (MP3) theory using the resolution of the identity (RI) approximation was combined with the fragment molecular orbital (FMO) method to efficiently calculate a high-order electron correlation energy of biomolecular systems. We developed a new algorithm for the RI-MP3 calculation, which can be used with the FMO scheme. After test calculations using a small molecule, the FMO-RI-MP3 calculations were performed for two biomolecular systems comprising a protein and a ligand. The computational cost of these calculations was only around 5 and 4 times higher than those of the FMO-RHF calculations. The error associated with the RI approximation was around 2.0% of the third-order correlation contribution to the total energy. However, the RI approximation error in the interaction energy between the protein and ligand molecule was insignificantly small, which reflected the negligible error in the inter fragment interaction energy. © 2018 Wiley Periodicals, Inc.

Polymer Segments at the Folding Limit: Raman Scattering for the Diglyme Benchmark.

Methyl-capped polyethers model flexible polar polymer chains. Their intrinsic all-trans conformational preference gives way to folded conformations if the chain is sufficiently long. We find by Raman cryospectroscopy in the gas phase and quantum chemical calculations that diglyme still prefers the stretched state, although folded variants come very close in energy. Three-body dispersion correction or higher-order electron correlation is important to match experiments, despite the polar character of C-O bonds. This turns oligoglymes into challenging benchmark systems for a balanced description of intra-chain torsion and inter-segment van der Waals attraction.

Ground-state ionization energies of boronlike ions

High-precision QED calculations of the ground-state ionization energies are performed for all boronlike ions with the nuclear charge numbers in the range $16 \leqslant Z\leqslant 96$. The rigorous QED calculations are performed within the extended Furry picture and include all many-electron QED effects up to the second order of the perturbation theory. The contributions of the third- and higher-order electron-correlation effects are accounted for within the Breit approximation. The nuclear recoil and nuclear polarization effects are taken into account as well. In comparison with the previous evaluations of the ground-state ionization energies of boronlike ions the accuracy of the theoretical predictions has been improved significantly.

Plain s-Wave Superconductivity near Magnetic Criticality: Enhancement of Attractive Electron–Boson Coupling Vertex Corrections

Recent experiments revealed that the plain $s$-wave state without any sign-reversal emerges in various metals near the magnetic criticality. To understand this counter-intuitive phenomenon, we study the gap equation for the multiorbital Hubbard-Holstein model, by analyzing the vertex correction (VC) due to the higher-order electron-correlation effects. We find that the phonon-mediated orbital fluctuations are magnified by the VC for the susceptibility ($\chi$-VC). In addition, the charge-channel attractive interaction is enlarged by the VC for the coupling-constant ($U$-VC), which is significant when the interaction has prominent $q$-dependences so the Migdal theorem fails. Due to both $\chi$-VC and $U$-VC, the plain $s$-wave state is caused by the small electron-phonon interaction near the magnetic criticality against the repulsive Coulomb interaction. We find that the direct Coulomb repulsion for the plain $s$-wave Cooper pair is strongly reduced by the "multiorbital screening effect".

Fragment molecular orbital (FMO) calculations on DNA by a scaled third-order Møller-Plesset perturbation (MP2.5) scheme

A series of fragment molecular orbital (FMO) calculations on DNA models was performed to assess the reliability of the scaled third-order Møller-Plesset perturbation method (MP2.5) (Pitoňák et al., Chem. Phys. Chem., 10 (2009) 282) for the incorporation of higher-order electron correlation. The Watson-Crick H-bonding (horizontal interaction) and the 1,2-stacking (vertical interaction) of four combinations of two base pairs (AT:AT, AT:TA, AT:GC, and CG:GC) were calculated at the levels of MP2, MP3 (nonscaled), MP2.5, and coupled cluster singles and doubles with perturbative triples (CCSD(T)). The comparison showed sufficient agreement between MP2.5 and the reference CCSD(T), and indicated that MP2.5 is an inexpensive but precise alternative to costly CCSD(T). The influence of basis sets including the model core potential (MCP) was also examined for the two base pairs. The Watson-Crick H-bonding interaction energies were fairly uninfluenced by these methods. In contrast, the values of the stacking interactions among bases quantitatively varied depending on the calculation options. As a cost-effective procedure, the FMO-MP2.5/MCP calculations were applied to a couple of realistic DNA models explicitly solvated by water: a 12-base-pair model and a 10-base-pair model cross-linked with cis-platin (Pt(NH3)2). The H-bonding and stacking interactions in them were discussed from a biochemical viewpoint.

Ab initio potential energy surface and vibration-rotation energy levels of beryllium monohydroxide.

The accurate potential energy surface of beryllium monohydroxide, BeOH, in its ground electronic state X 2A' has been determined from ab initio calculations using the coupled-cluster approach in conjunction with the correlation-consistent core-valence basis sets up to septuple-zeta quality. The higher-order electron correlation, scalar relativistic, and adiabatic effects were taken into account. The BeOH molecule was confirmed to be bent at equilibrium, with the BeOH angle of 141.2° and the barrier to linearity of 129 cm-1 . The vibration-rotation energy levels of the BeOH and BeOD isotopologues were predicted using a variational approach and compared with recent experimental data. The results can be useful in a further analysis of high-resolution vibration-rotation spectra of these interesting species. © 2016 Wiley Periodicals, Inc.

X-ray fluorescence spectrum of highly charged Fe ions driven by strong free-electron-laser fields

The influence of nonlinear dynamical effects is analyzed on the observed spectra of controversial 3C and 3D astrophysically relevant x-ray lines in neonlike Fe and the A, B, and C lines in natriumlike Fe ions. First, a large-scale configuration-interaction calculation of oscillator strengths is performed with the inclusion of higher-order electron-correlation effects. Also, quantum-electrodynamic corrections to the transition energies are calculated. Further considered dynamical effects provide a possible resolution of the discrepancy between theory and experiment found by recent x-ray free-electron-laser measurements of these controversial lines. We find that, for strong x-ray sources, the modeling of the spectral lines by a peak with an area proportional to the oscillator strength is not sufficient and nonlinear dynamical effects have to be taken into account. Thus, we advocate the use of light–matter-interaction models also valid for strong light fields in the analysis and interpretation of the associated astrophysical and laboratory spectra. We investigate line-strength ratios distinguishing between the coherent and incoherent parts of the emission spectrum. In addition, the spectrum of Fe, an autoionizing ion which was also present in the recent laboratory experiment, is analyzed as well.