Accelerating Pythonic Coupled-Cluster Implementations: A Comparison Between CPUs and GPUs

In this work, we benchmark several Python routines for time and memory requirements to identify the optimal choice of the tensor contraction operations available. We scrutinize how to accelerate the bottleneck tensor operations of Pythonic coupled-cluster implementations in the Cholesky linear algebra domain, utilizing a NVIDIA Tesla V100S PCIe 32GB (rev 1a) graphics processing unit (GPU). The NVIDIA compute unified device architecture API interacts with CuPy, an open-source library for Python, designed as a NumPy drop-in replacement for GPUs. Due to the limitations of video memory, the GPU calculations must be performed batch-wise. Timing results of some contractions containing large tensors are presented. The CuPy implementation leads to a factor of 10-16 speed-up of the bottleneck tensor contractions compared to computations on 36 central processing unit (CPU) cores. Finally, we compare example CCSD and pCCD-LCCSD calculations performed solely on CPUs to their CPU-GPU hybrid implementation, which leads to a speed-up of a factor of 3-4 compared to the CPU-only variant.

Geminal-Based Strategies for Modeling Large Building Blocks of Organic Electronic Materials

We elaborate on unconventional electronic structure methods based on geminals and their potential to advance the rapidly developing field of organic photovoltaics (OPVs). Specifically, we focus on the computational advantages of geminal-based methods over standard approaches and identify the critical aspects of OPV development. Examples are reliable and efficient computations of orbital energies, electronic spectra, and van der Waals interactions. Geminal-based models can also be combined with quantum embedding techniques and a quantum information analysis of orbital interactions to gain a fundamental understanding of the electronic structures and properties of realistic OPV building blocks. Furthermore, other organic components present in, for instance, dye-sensitized solar cells (DSSCs) represent another promising scope of application. Finally, we provide numerical examples predicting the properties of a small building block of OPV components and two carbazole-based dyes proposed as possible DSSC sensitizers.

Static embedding with pair coupled cluster doubles based methods

Quantum embedding methods have recently been significantly developed to model large molecular structures. This work proposes a novel wave function theory in a density functional theory (WTF-in-DFT) embedding scheme based on pair-coupled cluster doubles (pCCD)-type methods. While pCCD can reliably describe strongly-correlated systems with mean-field-like computational cost, the large extent of the dynamic correlation can be accounted for by (linearized) coupled-cluster corrections on top of the pCCD wave function. Here we focus on the linearized coupled-cluster singles and doubles (LCCSD) ansatz for electronic ground states and its extension to excited states within the equation of motion (EOM) formalism. We test our EOM-pCCD-LCCSD-in-DFT approach for the vertical excitation energies of the hydrogen-bonded water-ammonia complex, micro-solvated thymine, and uranyl tetrahalides (UO2X42-, X = F, Cl, Br). Furthermore, we assess the quality of the embedding potential using an orbital entanglement and correlation analysis. The approximate embedding models successfully capture changes in the excitation energies going from bare fragments to supramolecular structures and represent a promising computational method for excited states in large molecular systems.

The relationship between structure and excited-state properties in polyanilines from geminal-based methods

We employ state-of-the-art quantum chemistry methods to study the structure-to-property relationship in polyanilines (PANIs) of different lengths and oxidation states. Specifically, we focus on leucoemeraldine, emeraldine, and pernigraniline in their tetramer and octamer forms. We scrutinize their structural properties, HOMO and LUMO energies, HOMO-LUMO gaps, and vibrational and electronic spectroscopy using various Density Functional Approximations (DFAs). Furthermore, the accuracy of DFAs is assessed by comparing them to experimental and wavefunction-based reference data. We perform large-scale orbital-optimized pair-Coupled Cluster Doubles (oo-pCCD) calculations for ground and electronically excited states and conventional Configuration Interaction Singles (CIS) calculations for electronically excited states in all investigated systems. The EOM-pCCD+S approach with pCCD-optimized orbitals allows us to unambiguously identify charge transfer and local transitions across the investigated PANI systems-an analysis not possible within a delocalized canonical molecular orbital basis obtained, for instance, by DFAs. We show that the low-lying part of the emeraldine and pernigraniline spectrum is dominated by charge transfer excitations and that polymer elongation changes the character of the leading transitions. Furthermore, we augment our study with a quantum informational analysis of orbital correlations in various forms of PANIs.

Diffuse Basis Functions for Relativistic s and d Block Gaussian Basis Sets

Diffuse s, p, and d functions have been optimized for use with previously reported relativistic basis sets for the s and d blocks of the periodic table. The functions were optimized on the 4:1 weighted average of the s² and p² configurations of the anion, with the d shell in the dⁿ⁺¹ configuration for the d blocks. Exponents were extrapolated for groups 2 and 12, which have unstable or weakly bound anions. The diffuse basis sets have been tested by application to calculations of electron affinities of the group 11 elements (Cu, Ag, and Au), double electron affinities of the group 11 monocations, and potential energy curves of Mg₂ and Ca₂ van der Waals dimers, as well as some response properties of the group 1 anions (Rb^-, Cs^-, and Fr^-), the group 2 elements (Sr, Ba, and Ra), and RbLi, CsLi, and FrLi molecules.

Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry

In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.

Resolving the π-assisted U-N σf-bond formation using quantum information theory

We model the potential energy profiles of the UO₂ (NCO)Cl₂^- → NUOCl₂^- + CO₂ reaction pathway [Y. Gong, V. Vallet, M. del Carmen Michelini, D. Rios and J. K. Gibson, J. Phys. Chem. A, 2014, 118, 325-330] using different pair coupled-cluster doubles (pCCD) methods. Specifically, we focus on pCCD and pCCD-tailored coupled cluster models in predicting relative energies for the various intermediates and transition states along the reaction coordinate. Furthermore, we augment our study on energetics with an orbital-pair correlation analysis of the complete reaction pathway that features two distinct paths. Our analysis of orbital correlations sheds new light on the formation and breaking of respective bonds between the uranium, oxygen, and nitrogen atoms along the reaction coordinates where the "yl" bond is broken and a nitrido compound formed. Specifically, the strengthening of the U-N σ_f-bond is assisted by a π-type interaction that is delocalized over the C-N-U backbone of the UO₂ (NCO)Cl₂^- complex.

Multi-reference ab initio calculations of Hg spectral data and analysis of magic and zero-magic wavelengths

We have identified magic wavelengths for ¹S₀ ↔ ³P_1,2 (m_J = 0) transitions and zero-magic wavelengths for the ³P_1,2 (m_J = 0) states of ²⁰⁰Hg atoms, analysed the robustness of the magic conditions with respect to wavelength and polarization imperfections. We show that the most experimentally feasible magic wavelength for the ¹S₀ ↔ ³P₂ transition is 351.8 nm of π polarized light. Relevant transition wavelengths and transition strengths are calculated using the state-of-the-art Complete Active Space Self-Consistent-Field (CASSCF) method with a perturbative inclusion of spin-orbit coupling. The transition wavelengths are a posteriori corrected for the dynamical energy using the second-order perturbation theory.

Mixed uranyl and neptunyl cation–cation interaction-driven clusters: structures, energetic stability, and nuclear quadrupole interactions

We present a quantum-chemical study of mixed CCI clusters, their structures, energetic stability, and nuclear quadrupole interactions.

Assessing the accuracy of simplified coupled cluster methods for electronic excited states in f0 actinide compounds

We scrutinize the performance of different variants of equation of motion coupled cluster (EOM-CC) methods to predict electronic excitation energies and excited state potential energy surfaces in closed-shell actinide species.

Elucidating cation–cation interactions in neptunyl dications using multi-reference ab initio theory

Understanding the binding mechanism in neptunyl clusters formed due to cation–cation interactions is of crucial importance in nuclear waste reprocessing and related areas of research.

Electron correlation effects of the ThO and ThS molecules in the spinor basis. A relativistic coupled cluster study of ground and excited states properties

We present a comprehensive relativistic coupled cluster study of the electronic structures of the ThO and ThS molecules in the spinor basis. Specifically, we use the single-reference coupled cluster and the multi-reference Fock Space Coupled Cluster (FSCC) methods to model their ground and electronically-excited states. Two variants of the FSCC method have been investigated: (a) one where the electronic spectrum is obtained from sector (1,1) of the Fock space, and (b) another where the excited states come from the doubly attached electronic states to the doubly charged systems (ThO2+ and ThS2+), that is, from sector (0,2) of the Fock space. Our study provides a reliable set of spectroscopic parameters such as bond lengths, excitation energies, and vibrational frequencies, as well as a detailed analysis of the electron correlation effects in the ThO and ThS molecules. Finally, we examine the first ionization potential and electron affinity of the above mentioned molecules.

On the multi-reference nature of plutonium oxides: PuO22+, PuO2, PuO3 and PuO2(OH)2

Actinide-containing complexes present formidable challenges for electronic structure methods due to the large number of degenerate or quasi-degenerate electronic states arising from partially occupied 5f and 6d shells. Conventional multi-reference methods can treat active spaces that are often at the upper limit of what is required for a proper treatment of species with complex electronic structures, leaving no room for verifying their suitability. In this work we address the issue of properly defining the active spaces in such calculations, and introduce a protocol to determine optimal active spaces based on the use of the Density Matrix Renormalization Group algorithm and concepts of quantum information theory. We apply the protocol to elucidate the electronic structure and bonding mechanism of volatile plutonium oxides (PuO₃ and PuO₂(OH)₂), species associated with nuclear safety issues for which little is known about the electronic structure and energetics. We show how, within a scalar relativistic framework, orbital-pair correlations can be used to guide the definition of optimal active spaces which provide an accurate description of static/non-dynamic electron correlation, as well as to analyse the chemical bonding beyond a simple orbital model. From this bonding analysis we are able to show that the addition of oxo- or hydroxo-groups to the plutonium dioxide species considerably changes the π-bonding mechanism with respect to the bare triatomics, resulting in bent structures with a considerable multi-reference character.

Dissecting the cation–cation interaction between two uranyl units

A theoretical study of the CCIs between two bare uranyl units and their spectroscopic characterization.

Singlet ground state actinide chemistry with geminals

We present the first application of the variationally orbital optimized antisymmetric product of 1-reference orbital geminals (vOO-AP1roG) method to singlet-state actinide chemistry.

The effect of nitrido, azide, and nitrosyl ligands on magnetization densities and magnetic properties of iridium PNP pincer-type complexes

The magnetisation density of different iridium PNP pincer-type complexes depends on the chosen quantum chemistry method.

Quantum entanglement in carbon–carbon, carbon–phosphorus and silicon–silicon bonds

We present a quantum entanglement analysis to dissect the bond orders in polyatomic molecules.