Electronic structure simulations in the cloud computing environment.

The transformative impact of modern computational paradigms and technologies, such as high-performance computing (HPC), quantum computing, and cloud computing, has opened up profound new opportunities for scientific simulations. Scalable computational chemistry is one beneficiary of this technological progress. The main focus of this paper is on the performance of various quantum chemical formulations, ranging from low-order methods to high-accuracy approaches, implemented in different computational chemistry packages and libraries, such as NWChem, NWChemEx, Scalable Predictive Methods for Excitations and Correlated Phenomena, ExaChem, and Fermi-Löwdin orbital self-interaction correction on Azure Quantum Elements, Microsoft's cloud services platform for scientific discovery.

The accuracies of effective interactions in downfolding coupled-cluster approaches for small-dimensionality active spaces.

This paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism.

Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory.

Polariton chemistry has attracted great attention as a potential route to modify chemical structure, properties, and reactivity through strong interactions among molecular electronic, vibrational, or rovibrational degrees of freedom. A rigorous theoretical treatment of molecular polaritons requires the treatment of matter and photon degrees of freedom on equal quantum mechanical footing. In the limit of molecular electronic strong or ultrastrong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of quantum chemistry, yielding an approach we refer to as cavity quantum electrodynamics, where the photon degrees of freedom are treated at the level of cavity quantum electrodynamics.

Quantum Flow Algorithms for Simulating Many-Body Systems on Quantum Computers.

We conducted quantum simulations of strongly correlated systems using the quantum flow (QFlow) approach, which enables sampling large subspaces of the Hilbert space through coupled variational problems in reduced dimensionality active spaces. Our QFlow algorithms significantly reduce circuit complexity and pave the way for scalable and constant-circuit-depth quantum computing. Our simulations show that QFlow can optimize the collective number of wave function parameters without increasing the required qubits using active spaces having an order of magnitude fewer number of parameters.

NWChem: Recent and Ongoing Developments.

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0).

TAMM: Tensor algebra for many-body methods.

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods.

Real-Time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-Body Methods Infrastructure.

We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ 2020, 152, 174113] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions.

Optimized Quantum Phase Estimation for Simulating Electronic States in Various Energy Regimes.

While quantum algorithms for simulations exhibit better asymptotic scaling than their classical counterparts, they currently cannot be accurately implemented on real-world devices. Instead, chemists and computer scientists rely on costly classical simulations of these quantum algorithms. In particular, the quantum phase estimation (QPE) algorithm is among several approaches that has attracted much attention in recent years due to its genuine quantum character.

Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer.

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N - 1 electron systems.

Coupled cluster downfolding methods: The effect of double commutator terms on the accuracy of ground-state energies.

Downfolding coupled cluster techniques have recently been introduced into quantum chemistry as a tool for the dimensionality reduction of the many-body quantum problem. As opposed to earlier formulations in physics and chemistry based on the concept of effective Hamiltonians, the appearance of the downfolded Hamiltonians is a natural consequence of the single-reference exponential parameterization of the wave function. In this paper, we discuss the impact of higher-order terms originating in double commutators.

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations.

For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree-Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals.

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape.

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units.

Index of multi-determinantal and multi-reference character in coupled-cluster theory.

A full configuration interaction calculation (FCI) ultimately defines the innate molecular orbital description of a molecule. Its density matrix and the natural orbitals obtained from it quantify the difference between having N-dominantly occupied orbitals in a reference determinant for a wavefunction to describe N-correlated electrons and how many of those N-electrons are left to the remaining virtual orbitals. The latter provides a measure of the multi-determinantal character (MDC) required to be in a wavefunction.

Toward Quantum Computing for High-Energy Excited States in Molecular Systems: Quantum Phase Estimations of Core-Level States.

This paper explores the utility of the quantum phase estimation (QPE) algorithm in calculating high-energy excited states characterized by the promotion of electrons occupying core-level shells. These states have been intensively studied over the last few decades, especially in supporting the experimental effort at light sources. Results obtained with QPE are compared with various high-accuracy many-body techniques developed to describe core-level states.

Resource-Efficient Chemistry on Quantum Computers with the Variational Quantum Eigensolver and the Double Unitary Coupled-Cluster Approach.

Applications of quantum simulation algorithms to obtain electronic energies of molecules on noisy intermediate-scale quantum (NISQ) devices require careful consideration of resources describing the complex electron correlation effects. In modeling second-quantized problems, the biggest challenge confronted is that the number of qubits scales linearly with the size of the molecular basis. This poses a significant limitation on the size of the basis sets and the number of correlated electrons included in quantum simulations of chemical processes.

Sub-system quantum dynamics using coupled cluster downfolding techniques.

In this paper, we discuss extending the sub-system embedding sub-algebra coupled cluster formalism and the double unitary coupled cluster (DUCC) ansatz to the time domain. An important part of the analysis is associated with proving the exactness of the DUCC ansatz based on the general many-body form of anti-Hermitian cluster operators defining external and internal excitations. Using these formalisms, it is possible to calculate the energy of the entire system as an eigenvalue of downfolded/effective Hamiltonian in the active space, which is identifiable with the sub-system of the composite system.

Quantum simulations of excited states with active-space downfolded Hamiltonians.

Many-body techniques based on the double unitary coupled cluster (DUCC) ansatz can be used to downfold electronic Hamiltonians into low-dimensional active spaces. It can be shown that the resulting dimensionality reduced Hamiltonians are amenable for quantum computing. Recent studies performed for several benchmark systems using phase estimation (PE) algorithms for quantum computers demonstrated that these formulations can recover a significant portion of ground-state dynamical correlation effects that stem from the electron excitations outside of the active space.

The mitochondrial and chloroplast genomes of the green alga Haematococcus are made up of nearly identical repetitive sequences.

The chlamydomonadalean green alga Haematococcus lacustris (strain UTEX 2505) has the largest chloroplast genome on record: 1352 kb with ∼90% non-coding DNA [1,2]. But what of the mitochondrial genome? Here we present sequencing, assembly, and analysis of the mitogenome that shows that it, too, is extremely expanded. What's more, the same repetitive elements have spread throughout the mitochondrial and chloroplast (or plastid) DNA (mtDNA and ptDNA, respectively), resulting in the situation whereby these two distinct organelle genomes are made up of nearly identical sequences.

Downfolding of many-body Hamiltonians using active-space models: Extension of the sub-system embedding sub-algebras approach to unitary coupled cluster formalisms.

In this paper, we discuss the extension of the recently introduced subsystem embedding subalgebra coupled cluster (SES-CC) formalism to unitary CC formalisms. In analogy to the standard single-reference SES-CC formalism, its unitary CC extension allows one to include the dynamical (outside the active space) correlation effects in an SES induced complete active space (CAS) effective Hamiltonian. In contrast to the standard single-reference SES-CC theory, the unitary CC approach results in a Hermitian form of the effective Hamiltonian.

Unrestrained markerless trait stacking in through combined genome editing and marker recycling technologies.

Robust molecular tool kits in model and industrial microalgae are key to efficient targeted manipulation of endogenous and foreign genes in the nuclear genome for basic research and, as importantly, for the development of algal strains to produce renewable products such as biofuels. While Cas9-mediated gene knockout has been demonstrated in a small number of algal species with varying efficiency, the ability to stack traits or generate knockout mutations in two or more loci are often severely limited by selectable agent availability. This poses a critical hurdle in developing production strains, which require stacking of multiple traits, or in probing functionally redundant gene families.

Next-Generation Sequencing of Haematococcus lacustris Reveals an Extremely Large 1.35-Megabase Chloroplast Genome.

is an industrially relevant microalga that is used for the production of the carotenoid astaxanthin. Here, we report the use of PacBio long-read sequencing to assemble the chloroplast genome of strain UTEX:2505. At 1.

Application of the CC(P;Q) Hierarchy of Coupled-Cluster Methods to the Beryllium Dimer.

The performance of coupled-cluster approaches with higher-than-doubly excited clusters, including the CCSD(T), CCSD(2), CR-CC(2,3), CCSD(TQ), and CR-CC(2,4) corrections to CCSD, the active-space CCSDt, CCSDtq, and CCSDTq methods, and the CC(t;3), CC(t,q;3), CC(t,q;3,4), and CC(q;4) corrections to CCSDt, CCSDtq, and CCSDTq resulting from the CC(P;Q) formalism, in reproducing the CCSDT and CCSDTQ potential energy curves and vibrational term values characterizing Be in its electronic ground state is assessed. The correlation-consistent aug-cc-pVnZ and aug-cc-pCVnZ (n = T and Q) basis sets are employed. Among the CCSD-based corrections, the completely renormalized CR-CC(2,3) and CR-CC(2,4) approaches perform the best.

Coupled-cluster interpretation of the photoelectron spectrum of Ag3 (.).

We use the scalar relativistic ionized equation-of-motion coupled-cluster (IP-EOMCC) approaches to investigate the photoelectron spectrum of Ag3 (-), examining the effects of basis set, number of correlated electrons, level of applied theory including up to 3-hole-2-particle terms, and geometry relaxation. By employing an IP-EOMCC-based extrapolation scheme, we are able to provide an accurate interpretation and complete assignment of peaks and other key features in the experimentally observed spectra, including electron binding energies as high as about 6.5 eV.

Communication: Coupled-cluster interpretation of the photoelectron spectrum of Au₃⁻.

We use the scalar relativistic ionized equation-of-motion coupled-cluster approaches, correlating valence and semi-core electrons and including up to 3-hole-2-particle terms in the ionizing operator, to investigate the photoelectron spectrum of Au₃⁻. We provide an accurate assignment of peaks and shoulders in the experimental photoelectron spectrum of Au₃⁻ for the first time.