Search

Spectroscopy simulations are of paramount importance for the interpretation of experimental electronic spectra, the disentangling of overlapping spectral features, and the tracing of the microscopic origin of the observed signals. Linear and nonlinear simulations are based on the results drawn from electronic structure calculations that provide the necessary parameterization of the molecular systems probed by light. Here, we investigate the applicability of excited-state properties obtained from linear-response time-dependent density functional theory (TDDFT) in the description of nonlinear spectra by employing the pseudowavefunction approach and compare them with benchmarks from highly accurate RASSCF/RASPT2 calculations and with high temporal resolution experimental results. As a test case, we consider the prediction of femtosecond transient absorption and two-dimensional electronic spectroscopy of a perylene bisimide dye in solution. We find that experimental signals are well reproduced by both theoretical approaches, showing that the computationally cheaper TDDFT can be a suitable option for the simulation of nonlinear spectroscopy of molecular systems that are too large to be treated with higher-level RASSCF/RASPT2 methods.

Resonant stimulated X-ray Raman spectroscopy of the bimetallic [Mn III Mn IV (μ-O) 2 (μ-OAC)(tacn) 2 ] 2+ manganese complex is investigated in a simulation study. Essential biological processes, including water oxidation in photosynthesis, involve charge transfer between manganese sites of different oxidation states. We study a prototypical binuclear mixed-valence transition-metal complex with two Mn atoms in different oxidation states surrounded by ligand structures and employ a pump-probe sequence of resonant X-ray Raman excitations to follow the charge transfer occurring in the molecule. This allows us to generate and monitor valence-electron wave packets at selected regions in the molecule by exploiting element-specific core-excited states. A two-color protocol is presented, with pump and probe pulses tuned to the Mn and N K-edges. A natural orbital decomposition allows the visualization of the electron dynamics underlying the signal.

A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [ J. Chem. Phys. 1986 , 85 , 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.

We present a time-dependent density functional theory (TDDFT) approach to compute the light-matter couplings between two different manifolds of excited states relative to a common ground state in the context of 4d transition metal systems. These quantities are the necessary ingredients to solve the Kramers-Heisenberg (KH) equation for resonant inelastic X-ray scattering (RIXS) and several other types of two-photon spectroscopies. The procedure is based on the pseudo-wavefunction approach, where the solutions of a TDDFT calculation can be used to construct excited-state wavefunctions, and on the restricted energy window approach, where a manifold of excited states can be rigorously defined based on the energies of the occupied molecular orbitals involved in the excitation process. Thus, the present approach bypasses the need to solve the costly TDDFT quadratic-response equations. We illustrate the applicability of the method to 4d transition metal molecular complexes by calculating the 2p4d RIXS maps of three representative ruthenium complexes and comparing them to experimental results. The method can capture all the experimental features in all three complexes to allow the assignment of the experimental peaks, with relative energies correct to within ∼0.6 eV at the cost of two independent TDDFT calculations.

Psi4NumPy demonstrates the use of efficient computational kernels from the open-source Psi4 program through the popular NumPy library for linear algebra in Python to facilitate the rapid development of clear, understandable Python computer code for new quantum chemical methods, while maintaining a relatively low execution time. Using these tools, reference implementations have been created for a number of methods, including self-consistent field (SCF), SCF response, many-body perturbation theory, coupled-cluster theory, configuration interaction, and symmetry-adapted perturbation theory. Furthermore, several reference codes have been integrated into Jupyter notebooks, allowing background, underlying theory, and formula information to be associated with the implementation. Psi4NumPy tools and associated reference implementations can lower the barrier for future development of quantum chemistry methods. These implementations also demonstrate the power of the hybrid C++/Python programming approach employed by the Psi4 program.

The Hartree-Fock problem was recently recast as a semidefinite optimization over the space of rank-constrained two-body reduced-density matrices (RDMs) [ Phys. Rev. A 2014 , 89 , 010502(R) ]. This formulation of the problem transfers the nonconvexity of the Hartree-Fock energy functional to the rank constraint on the two-body RDM. We consider an equivalent optimization over the space of positive semidefinite one-electron RDMs (1-RDMs) that retains the nonconvexity of the Hartree-Fock energy expression. The optimized 1-RDM satisfies ensemble N-representability conditions, and ensemble spin-state conditions may be imposed as well. The spin-state conditions place additional linear and nonlinear constraints on the 1-RDM. We apply this RDM-based approach to several molecular systems and explore its spatial (point group) and spin ( Ŝ 2 and Ŝ 3 ) symmetry breaking properties. When imposing Ŝ 2 and Ŝ 3 symmetry but relaxing point group symmetry, the procedure often locates spatial-symmetry-broken solutions that are difficult to identify using standard algorithms. For example, the RDM-based approach yields a smooth, spatial-symmetry-broken potential energy curve for the well-known Be-H 2 insertion pathway. We also demonstrate numerically that, upon relaxation of Ŝ 2 and Ŝ 3 symmetry constraints, the RDM-based approach is equivalent to real-valued generalized Hartree-Fock theory.

An explicitly time-dependent (TD) approach to equation-of-motion (EOM) coupled-cluster theory with single and double excitations (CCSD) is implemented for simulating near-edge X-ray absorption fine structure in molecular systems. The TD-EOM-CCSD absorption line shape function is given by the Fourier transform of the CCSD dipole autocorrelation function. We represent this transform by its Padé approximant, which provides converged spectra in much shorter simulation times than are required by the Fourier form. The result is a powerful framework for the blackbox simulation of broadband absorption spectra. K-edge X-ray absorption spectra for carbon, nitrogen, and oxygen in several small molecules are obtained from the real part of the absorption line shape function and are compared with experiment. The computed and experimentally obtained spectra are in good agreement; the mean unsigned error in the predicted peak positions is only 1.2 eV. We also explore the spectral signatures of protonation in these molecules.

We report an explicitly time-dependent approach to the generation of linear absorption spectra for molecular systems within the framework of equation-of-motion (EOM) coupled-cluster (CC) theory. While most time-dependent CC approaches consider the perturbation and time-evolution of a CC wave function, the present work considers the time-evolution of a CC dipole function. The dipole function formalism introduces no approximations and requires the evolution of only one time-dependent quantity, either the left or right dipole function. This time-dependent framework can be used to compute linear absorption spectra for molecules with a high density of states over a broad spectral range, a case for which conventional frequency-domain computations may become impractical. We validate the approach by comparing absorption spectra for small molecules computed at EOM second-order approximate CC (CC2) and time-dependent EOM-CC2 (TD-EOM-CC2) levels of theory. TD-EOM-CC2 computations are also used to predict extreme ultraviolet absorption spectra for third-row ions that are in reasonable agreement with experiment.

We present a parametrized coupled-pair functional, optimized to describe molecular interactions (PCPF-MI). The method has the same computational cost as singles and doubles configuration interaction, is size extensive, and yields energies that are stationary with respect to variations in all of its excitation coefficients. This last property facilitates the construction of density matrices and the evaluation of one- and two-electron properties. For the S22, HSG, S66, and A24 databases of van der Waals dimers, PCPF-MI computations yield interaction energies with mean unsigned errors of 0.326, 0.149, 0.214, and 0.044 kcal mol(-1), respectively, relative to benchmark computations. PCPF-MI interaction energies consistently improve upon those obtained from several other coupled-pair methods, including the averaged coupled-pair functional (ACPF), the coupled electron-pair approximation [CEPA(n), n = 0,1,3], and averaged quadratic coupled cluster (AQCC). Optimal parameters for spin-component-scaled (SCS) variants of each of these methods are also presented. SCS-CEPA(0), SCS-CEPA(1), SCS-ACPF, and SCS-PCPF-MI all perform similarly, with average errors for the four databases of roughly 0.1 kcal mol(-1). The PCPF-MI method, without additional SCS parametrization, has an average error of 0.192 kcal mol(-1) over the four databases and, when compared to the SCS-parametrized coupled-pair methods, has the added benefit that the energy is stationary with respect to variations in all its excitation coefficients.