Advanced calculations of X-ray spectroscopies with FEFF10 and Corvus.

The real-space Green's function code FEFF has been extensively developed and used for calculations of X-ray and related spectra, including X-ray absorption (XAS), X-ray emission (XES), inelastic X-ray scattering, and electron energy-loss spectra. The code is particularly useful for the analysis and interpretation of the XAS fine-structure (EXAFS) and the near-edge structure (XANES) in materials throughout the periodic table. Nevertheless, many applications, such as non-equilibrium systems, and the analysis of ultra-fast pump-probe experiments, require extensions of the code including finite-temperature and auxiliary calculations of structure and vibrational properties. To enable these extensions, we have developed in tandem a new version FEFF10 and new FEFF-based workflows for the Corvus workflow manager, which allow users to easily augment the capabilities of FEFF10 via auxiliary codes. This coupling facilitates simplified input and automated calculations of spectra based on advanced theoretical techniques. The approach is illustrated with examples of high-temperature behavior, vibrational properties, many-body excitations in XAS, super-heavy materials, and fits of calculated spectra to experiment.

Equation of motion coupled-cluster cumulant approach for intrinsic losses in x-ray spectra.

We present a combined equation of motion coupled-cluster cumulant Green's function approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS) and x-ray photoemission spectra. The method is based on a factorization of the transition amplitude in the time domain, which leads to a convolution of an effective one-body absorption spectrum and the core-hole spectral function. The spectral function characterizes intrinsic losses in terms of shake-up excitations and satellites using a cumulant representation of the core-hole Green's function that simplifies the interpretation. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.

Corvus: a framework for interfacing scientific software for spectroscopic and materials science applications.

Corvus, a Python-based package designed for managing workflows of physical simulations that utilize multiple scientific software packages, is presented. Corvus can be run as an executable script with an input file and automatically generated or custom workflows, or interactively, in order to build custom workflows with a set of Corvus-specific tools. Several prototypical examples are presented that link density functional, vibrational and X-ray spectroscopy software packages and are of interest to the synchrotron community. These examples highlight the simplification of complex spectroscopy calculations that were previously limited to expert users, and demonstrate the flexibility of the Corvus infrastructure to tackle more general problems in other research areas.

Dynamic structural disorder in supported nanoscale catalysts.

We investigate the origin and physical effects of "dynamic structural disorder" (DSD) in supported nano-scale catalysts. DSD refers to the intrinsic fluctuating, inhomogeneous structure of such nano-scale systems. In contrast to bulk materials, nano-scale systems exhibit substantial fluctuations in structure, charge, temperature, and other quantities, as well as large surface effects. The DSD is driven largely by the stochastic librational motion of the center of mass and fluxional bonding at the nanoparticle surface due to thermal coupling with the substrate. Our approach for calculating and understanding DSD is based on a combination of real-time density functional theory/molecular dynamics simulations, transient coupled-oscillator models, and statistical mechanics. This approach treats thermal and dynamic effects over multiple time-scales, and includes bond-stretching and -bending vibrations, and transient tethering to the substrate at longer ps time-scales. Potential effects on the catalytic properties of these clusters are briefly explored. Model calculations of molecule-cluster interactions and molecular dissociation reaction paths are presented in which the reactant molecules are adsorbed on the surface of dynamically sampled clusters. This model suggests that DSD can affect both the prefactors and distribution of energy barriers in reaction rates, and thus can significantly affect catalytic activity at the nano-scale.

A transferable H2O interaction potential based on a single center multipole expansion: SCME.

A transferable potential energy function for describing the interaction between water molecules is presented. The electrostatic interaction is described rigorously using a multipole expansion. Only one expansion center is used per molecule to avoid the introduction of monopoles. This single center approach turns out to converge and give close agreement with ab initio calculations when carried out up to and including the hexadecapole. Both dipole and quadrupole polarizability are included. All parameters in the electrostatic interaction as well as the dispersion interaction are taken from ab initio calculations or experimental measurements of a single water molecule. The repulsive part of the interaction is parametrized to fit ab initio calculations of small water clusters and experimental measurements of ice Ih. The parametrized potential function was then used to simulate liquid water and the results agree well with experiment, even better than simulations using some of the point charge potentials fitted to liquid water. The evaluation of the new interaction potential for condensed phases is fast because point charges are not present and the interaction can, to a good approximation, be truncated at a finite range.

Nonstoichiometric intensities in core photoelectron spectroscopy.

X-ray photoemission spectroscopy is used in a great variety of research fields; one observable is the sample's stoichiometry. The stoichiometry can be deduced based on the expectation that the ionization cross sections for innershell orbitals are independent of the molecular composition. Here we used chlorine-substituted ethanes in the gas phase to investigate the apparent carbon stoichiometry. We observe a nonstoichiometric ratio for a wide range of photon energies, the ratio exhibits x-ray-absorption fine structure spectroscopy (EXAFS)-like oscillations and hundreds of eV above the C1s ionization approaches a value far from 1. These effects can be accounted for by considering the scattering of the outgoing photoelectron, which we model by multiple-scattering EXAFS calculations, and by considering the effects of losses due to monopole shakeup and shakeoff and to intramolecular inelastic scattering processes.

Real-Time Equation-of-Motion CCSD Cumulant Green's Function.

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 10e 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 CH4, NH3, and H2O, 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.

Equation-of-Motion Coupled-Cluster Cumulant Green's Function for Excited States and X-Ray Spectra.

Green's function methods provide a robust, general framework within many-body theory for treating electron correlation in both excited states and x-ray spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green's function approach has been introduced, where the cumulant is obtained as the solution to a system of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green's function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function for the core Green's function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then calculated using a convolution of the core spectral function and an effective, one-body XAS. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green's function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. The method is illustrated for the XPS and XAS of NH3.

Real-Time Coupled-Cluster Approach for the Cumulant Green's Function.

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.

The local electronic structure of alpha-Li3N.

New theoretical and experimental investigations of the occupied and unoccupied local electronic densities of states (DOS) are reported for alpha-Li(3)N. Band-structure and density-functional theory calculations confirm the absence of covalent bonding character. However, real-space full-multiple-scattering (RSFMS) calculations of the occupied local DOS find less extreme nominal valences than have previously been proposed. Nonresonant inelastic x-ray scattering, RSFMS calculations, and calculations based on the Bethe-Salpeter equation are used to characterize the unoccupied electronic final states local to both the Li and N sites. There is a good agreement between experiment and theory. Throughout the Li 1s near-edge region, both experiment and theory find strong similarities in the s-and p-type components of the unoccupied local final DOS projected onto an orbital angular momentum basis (l-DOS). An unexpected, significant correspondence exists between the near-edge spectra for the Li 1s and N 1s initial states. We argue that both spectra are sampling essentially the same final DOS due to the combination of long core-hole lifetimes, long photoelectron lifetimes, and the fact that orbital angular momentum is the same for all relevant initial states. Such considerations may be generally applicable for low atomic number compounds.

Quantitative first-principles calculations of valence and core excitation spectra of solid C60.

We present calculated valence and C 1s near-edge excitation spectra of solid C60 and experimental results measured with high-resolution electron energy-loss spectroscopy. The near-edge calculations are carried out using three different methods: solution of the Bethe-Salpeter equation (BSE) as implemented in the OCEAN suite (Obtaining Core Excitations with ab initio methods and the NIST BSE solver), the excited-electron core-hole approach (XCH), and the constrained-occupancy method using the Stockholm-Berlin core-excitation code, StoBe. The three methods give similar results and are in good agreement with experiment, though the BSE results are the most accurate. The BSE formalism is also used to carry out valence level calculations using the NIST Bethe-Salpeter Equation solver (NBSE). Theoretical results include self-energy corrections to the band gap and band widths, lifetime-damping effects, and Debye-Waller effects in the core-excitation case. A comparison of spectral features to those observed experimentally illustrates the sensitivity of certain features to computational details, such as self-energy corrections to the band structure and core-hole screening.

Real-time time-dependent density functional theory approach for frequency-dependent nonlinear optical response in photonic molecules.

We present ab initio calculations of frequency-dependent linear and nonlinear optical responses based on real-time time-dependent density functional theory for arbitrary photonic molecules. This approach is based on an extension of an approach previously implemented for a linear response using the electronic structure program SIESTA. Instead of calculating excited quantum states, which can be a bottleneck in frequency-space calculations, the response of large molecular systems to time-varying electric fields is calculated in real time. This method is based on the finite field approach generalized to the dynamic case. To speed the nonlinear calculations, our approach uses Gaussian enveloped quasimonochromatic external fields. We thereby obtain the frequency-dependent second harmonic generation beta(-2omega;omega,omega), the dc nonlinear rectification beta(0;-omega,omega), and the electro-optic effect beta(-omega;omega,0). The method is applied to nanoscale photonic nonlinear optical molecules, including p-nitroaniline and the FTC chromophore, i.e., 2-[3-Cyano-4-(2-{5-[2-(4-diethylamino-phenyl)-vinyl]-thiophen-2-yl}-vinyl)-5,5-dimethyl-5H-furan-2-ylidene]-malononitrile, and yields results in good agreement with experiment.

Comparison of static first hyperpolarizabilities calculated with various quantum mechanical methods.

The prediction of nonlinear electro-optic (EO) behavior of molecules with quantum methods is the first step in the development of organic-based electro-optic devices. Typical EO molecules may require calculations with several hundred electrons, which prevents all but the fastest methods (semiempirical and density functional theory (DFT)) from being used for EO estimation. To test the reliability of these methods, we compare dipole moments, polarizabilities, and first-order hyperpolarizabilities for a wide range of structures of experimental interest with Hartree-Fock (HF), intermediate neglect of differential overlap (INDO), and DFT methods. The relative merits of molecules are consistently predictable with every method.