We present two efficient iterative algorithms for solving the linear response eigenvalue problem arising from the time dependent density functional theory. Although the matrix to be diagonalized is nonsymmetric, it has a special structure that can be exploited to save both memory and floating point operations. In particular, the nonsymmetric eigenvalue problem can be transformed into an eigenvalue problem that involves the product of two matrices M and K. We show that, because MK is self-adjoint with respect to the inner product induced by the matrix K, this product eigenvalue problem can be solved efficiently by a modified Davidson algorithm and a modified locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm that make use of the K-inner product. The solution of the product eigenvalue problem yields one component of the eigenvector associated with the original eigenvalue problem. We show that the other component of the eigenvector can be easily recovered in an inexpensive postprocessing procedure. As a result, the algorithms we present here become more efficient than existing methods that try to approximate both components of the eigenvectors simultaneously. In particular, our numerical experiments demonstrate that the new algorithms presented here consistently outperform the existing state-of-the-art Davidson type solvers by a factor of two in both solution time and storage.

Long-range electron transfer (ET) is a crucial step in many energy conversion processes and biological redox reactions in living organisms. We show that newly developed X-ray pulses can directly probe the evolving oxidation states and the electronic structure around selected atoms with detail not available through conventional time-resolved infrared or optical techniques. This is demonstrated in a simulation study of the stimulated X-ray Raman (SXRS) signals in Re-modified azurin, which serves as a benchmark system for photoinduced ET in proteins. Nonlinear SXRS signals offer a direct novel window into the long-range ET mechanism.

We demonstrate how optical cavities can be exploited to control both valence- and core-excitations in a prototypical model transition metal complex, ferricyanide ([Fe(iii)(CN)₆]^3-), in an aqueous environment. The spectroscopic signatures of hybrid light-matter polariton states are revealed in UV/Vis and X-ray absorption, and stimulated X-ray Raman signals. In an UV/Vis cavity, the absorption spectrum exhibits the single-polariton states arising from the cavity photon mode coupling to both resonant and off-resonant valence-excited states. We further show that nonlinear stimulated X-ray Raman signals can selectively probe the bipolariton states via cavity-modified Fe core-excited states. This unveils the correlation between valence polaritons and dressed core-excitations. In an X-ray cavity, core-polaritons are generated and their correlations with the bare valence-excitations appear in the linear and nonlinear X-ray spectra.

Core-resonant circular dichroism (CD) signals are induced by molecular chirality and vanish for achiral molecules and racemic mixtures. The highly localized nature of core excitations makes them ideal probes of local chirality within molecules. Simulations of the circular dichroism spectra of several molecular families illustrate how these signals vary with the electronic coupling to substitution groups, the distance between the X-ray chromophore and the chiral center, geometry, and chemical structure. Clear insight into the molecular structure is obtained through analysis of the X-ray CD spectra.

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

Molecular clustering and associated dynamic processes of lithium polysulfide species were unraveled using classical molecular dynamics and ab initio metadynamics calculations. The spectroscopic signatures of polysulfide clusters were analyzed using a multimodal analysis including experimental and computational nuclear magnetic resonance (NMR) and X-ray absorption spectroscopies. Lithium polysulfide solutes (Li2S4, Li2S6, and Li2S8) and their mixtures in a 1,3-dioxolane and 1,2-dimethoxyethane (DOL/DME) solvent undergo aggregation driven by intramolecular lithium-sulfur (Li-S) interactions, leading to distributions of cluster sizes, which could critically influence the functioning of lithium-sulfur batteries. Representative polysulfide clusters with systematic increases in molecular size were extracted from the classical molecular dynamics trajectories for subsequent structural and spectroscopic property calculations using density functional theory analysis. Structural analysis of these clusters reveals progressively decreasing solvent involvement in Li+ coordination varying from Li2S4 to Li2S8, with more pronounced variation and changes in DME compared with those in DOL. These observations are reflected in the analysis of the experimental and theoretical 7Li and 17O NMR chemical shifts and pulsed field gradient-NMR diffusion measurements. A comparison of experimental and theoretical S K-edge X-ray absorption near edge spectra shows that relatively large lithium sulfide chain clusters are likely to occur in the DOL/DME-solvated lithium sulfide systems. Ab initio metadynamics simulations and NMR analysis indicate that Li+ solvated by only the solvent can occur through Li+ dissociation from sulfide chains. However, the occurrence of "sulfide-free" Li+ is a minor mechanism compared with the dynamic aggregation and shuttling processes of polysulfide solvates in DOL/DME-based electrolytes of Li-S batteries. Overall, atomistic insights gained about clustering and lithium exchange dynamics will be critical for the predictive understanding of the polysulfide shuttling and nucleation process that dictates the Li-S battery performance.

Mechanisms of CaCO₃ nucleation from solutions that depend on multistage pathways and the existence of species far more complex than simple ions or ion pairs have recently been proposed. Herein, we provide a tightly coupled theoretical and experimental study on the pathways that precede the initial stages of CaCO₃ nucleation. Starting from molecular simulations, we succeed in correctly predicting bulk thermodynamic quantities and experimental data, including equilibrium constants, titration curves, and detailed x-ray absorption spectra taken from the supersaturated CaCO₃ solutions. The picture that emerges is in complete agreement with classical views of cluster populations in which ions and ion pairs dominate, with the concomitant free energy landscapes following classical nucleation theory.

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.

We present a joint experimental and computational study of the hexacyanoferrate aqueous complexes at equilibrium in the 250 meV to 7.15 keV regime. The experiments and the computations include the vibrational spectroscopy of the cyanide ligands, the valence electronic absorption spectra, and Fe 1s core hole spectra using element-specific-resonant X-ray absorption and emission techniques. Density functional theory-based quantum mechanics/molecular mechanics molecular dynamics simulations are performed to generate explicit solute-solvent configurations, which serve as inputs for the spectroscopy calculations of the experiments spanning the IR to X-ray wavelengths. The spectroscopy simulations are performed at the same level of theory across this large energy window, which allows for a systematic comparison of the effects of explicit solute-solvent interactions in the vibrational, valence electronic, and core-level spectra of hexacyanoferrate complexes in water. Although the spectroscopy of hexacyanoferrate complexes in solution has been the subject of several studies, most of the previous works have focused on a narrow energy window and have not accounted for explicit solute-solvent interactions in their spectroscopy simulations. In this work, we focus our analysis on identifying how the local solvation environment around the hexacyanoferrate complexes influences the intensity and line shape of specific spectroscopic features in the UV/vis, X-ray absorption, and valence-to-core X-ray emission spectra. The identification of these features and their relationship to solute-solvent interactions is important because hexacyanoferrate complexes serve as model systems for understanding the photochemistry and photophysics of a large class of Fe(II) and Fe(III) complexes in solution.

Quantifying charge delocalization associated with short-lived photoexcited states of molecular complexes in solution remains experimentally challenging, requiring local element specific femtosecond experimental probes of time-evolving electron transfer. In this study, we quantify the evolving valence hole charge distribution in the photoexcited charge transfer state of a prototypical mixed valence bimetallic iron-ruthenium complex, [(CN)₅Fe^IICNRu^III(NH₃)₅]^-, in water by combining femtosecond X-ray spectroscopy measurements with time-dependent density functional theory calculations of the excited-state dynamics. We estimate the valence hole charge that accumulated at the Fe atom to be 0.6 ± 0.2, resulting from excited-state metal-to-metal charge transfer, on an ∼60 fs time scale. Our combined experimental and computational approach provides a spectroscopic ruler for quantifying excited-state valency in solvated complexes.