Molecular dynamics simulations, conventional and metadynamics, were performed to determine the interaction of model protein Gb1 over kaolinite (001), Na(+)-montmorillonite (001), Ca(2+)-montmorillonite (001), goethite (100), and Na(+)-birnessite (001) mineral surfaces. Gb1, a small (56 residue) protein with a well-characterized solution-state nuclear magnetic resonance (NMR) structure and having α-helix, 4-fold β-sheet, and hydrophobic core features, is used as a model protein to study protein soil mineral interactions and gain insights on structural changes and potential degradation of protein. From our simulations, we observe little change to the hydrated Gb1 structure over the kaolinite, montmorillonite, and goethite surfaces relative to its solvated structure without these mineral surfaces present. Over the Na(+)-birnessite basal surface, however, the Gb1 structure is highly disturbed as a result of interaction with this birnessite surface. Unraveling of the Gb1 β-sheet at specific turns and a partial unraveling of the α-helix is observed over birnessite, which suggests specific vulnerable residue sites for oxidation or hydrolysis possibly leading to fragmentation.

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.

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.