Search

Photoelectron spectroscopy (PES) is a standard experimental method for material characterization, but its interpretation can be hampered by its reliance on standard materials. To facilitate the study of unknown systems, theoretical methods are desirable. Here, we present a real-time equation-of-motion coupled cluster (RT-EOM-CC) approach for valence PES, extending our core-level development. We demonstrate that RT-EOM-CC yields ionization energies and spectral functions in good agreement with experimental and CI-based results, even for some more correlated cases.

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ J. Chem. Phys. 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. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams-Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photoemission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.

Many-body excitations in X-ray photoemission spectra have been difficult to simulate from first principles. We have recently developed a cumulant-based one-electron Green's function method using the real-time coupled-cluster-singles equation-of-motion approach (RT-EOM-CCS) that provides a general framework for treating these problems. Here we extend this approach to include double excitations in the ground-state energy and in the coupled cluster amplitudes, which have been implemented using subroutines generated by the Tensor Contraction Engine (TCE). As in the case of the singles approximation, RT-EOM-CCSD yields a nonperturbative cumulant form of the Green's function in terms of the time-dependent cluster amplitudes, adding nonlinear corrections to the traditional cumulant forms. The extended approach is applied to the core-hole spectral function for small molecular systems. We find that, when core-optimized basis sets are used, the doubles contributions reduce the mean absolute errors in the core binding energies of the 10 e systems from 0.8 to 0.3 eV. They also significantly improve the quasiparticle-satellite gap by reducing its overestimation from about 3-5 to about 0-1 eV in CH 4 , NH 3 , and H 2 O, and also improving the overall shape of the satellite features. Finally, we demonstrate the application of the new implementation to the larger, classical XPS ESCA series of molecules and show that the singles approximation can be paired with a modest basis set to study carbon speciation.

Green's function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Here, we develop the cumulant form of the one-electron Green's function using a real-time coupled-cluster equation-of-motion approach, in an extension of our previous study (Rehr J.; et al. J. Chem. Phys. 2020, 152, 174113). The approach yields a nonperturbative expression for the cumulant in terms of the solution to a set of coupled first-order, nonlinear differential equations. The method thereby adds nonlinear corrections to traditional cumulant methods, which are linear in the self-energy. The approach is applied to the core-hole Green's function and is illustrated for a number of small molecular systems. For these systems, we find that the nonlinear contributions yield significant improvements, both for quasiparticle properties such as core-level binding energies and for inelastic losses that correspond to satellites observed in photoemission spectra.