X-ray absorption spectroscopy and theoretical investigations of the effect of extended ligands in potassium organic matter interaction.

Potassium (K) is an essential nutrient for plant growth, and despite its abundance in soil, most of the K is structurally bound in minerals, limiting its bioavailability and making this soil K reservoir largely inaccessible to plants. Microbial biochemical weathering has been shown to be a promising pathway to sustainably increase plant available K. However, the mechanisms underpinning microbial K uptake, transformation, storage, and sharing are poorly resolved. To better understand the controls on microbial K transformations, we performed K K-edge x-ray absorption near-edge structure (XANES) spectroscopy on K-organic salts, including acetate, citrate, nitrate, oxalate, and tartrate, which are frequently observed as low molecular weight organic acids secreted by soil microbes, as well as humic acid, which acts as a proxy for higher molecular weight organic acids. The organic salts display feature-rich K XANES spectra, each demonstrating numerous unique features spanning ∼13 eV range across the absorption edge. In contrast, the spectra for humic acid have one broad, wide feature across the same energy range. We used a combination of time-dependent density functional theory and the Bethe-Salpeter equation based approach within the OCEAN code to simulate the experimental spectra for K-nitrate (KNO3) and K-citrate [K3(C6H5O7)·H2O] to identify the electronic transitions that give rise to some of the outlying and unique spectral features in the organic salts. KNO3 has both the lowest and highest lying energy features, and K3(C6H5O7)·H2O is produced by several soil microbes and is effective at mineral weathering. Our results analyze the K-organic salt bonding in detail to elucidate why the spectral shapes differ and indicate that the K K-edge XANES spectra are associated with the entire ligand despite similar first-shell bonding environments around the K center. The improved understanding of K bonding environments with organic ligands and their use for interpretation of the K-XANES spectra provides an important toolkit to understand how K is transformed by microbial processes and made bioavailable for plant uptake.

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties.

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute-solvent interactions in modulating the excited-state dynamics. To tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV-vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC)5FeIII(µ-CN)RuII(pyridine)4(µ-NC)FeIII(CN)5]4-). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute-solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm-Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

Atomistic simulations for investigation of substrate effects on lipid in-source fragmentation in secondary ion mass spectrometry.

In beam-based ionization methods, the substrate plays an important role on the desorption mechanism of molecules from surfaces. Both the specific orientation that a molecule adopts at a surface and the strength of the molecule-surface interaction can greatly influence desorption processes, which in turn will affect the ion yield and the degree of in-source fragmentation of a molecule. In the beam-based method of secondary ion mass spectrometry (SIMS), in-source fragmentation can be significant and molecule specific due to the hard ionization method of using a primary ion beam for molecule desorption. To investigate the role of the substrate on orientation and in-source fragmentation, we have used atomistic simulations-molecular dynamics in combination with density functional theory calculations-to explore the desorption of a sphingolipid (palmitoylsphingomyelin) from a model surface (gold). We then compare SIMS data from this model system to our modeling findings. Using this approach, we found that the combined adsorption and binding energy of certain bonds associated with the headgroup fragments (C3H8N+, C5H12N+, C5H14NO+, and C5H15PNO4 +) was a good predictor for fragment intensities (as indicated by relative ion yields). This is the first example where atomistic simulations have been applied in beam-based ionization of lipids, and it presents a new approach to study biointerfacial lipid ordering effects on SIMS imaging.

SOMAS: a platform for data-driven material discovery in redox flow battery development.

Aqueous organic redox flow batteries offer an environmentally benign, tunable, and safe route to large-scale energy storage. The energy density is one of the key performance parameters of organic redox flow batteries, which critically depends on the solubility of the redox-active molecule in water. Prediction of aqueous solubility remains a challenge in chemistry. Recently, machine learning models have been developed for molecular properties prediction in chemistry and material science. The fidelity of a machine learning model critically depends on the diversity, accuracy, and abundancy of the training datasets. We build a comprehensive open access organic molecular database "Solubility of Organic Molecules in Aqueous Solution" (SOMAS) containing about 12,000 molecules that covers wider chemical and solubility regimes suitable for aqueous organic redox flow battery development efforts. In addition to experimental solubility, we also provide eight distinctive quantum descriptors including optimized geometry derived from high-throughput density functional theory calculations along with six molecular descriptors for each molecule. SOMAS builds a critical foundation for future efforts in artificial intelligence-based solubility prediction models.

Femtosecond X-ray Spectroscopy Directly Quantifies Transient Excited-State Mixed Valency.

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)5FeIICNRuIII(NH3)5]-, 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.

Graphical Gaussian process regression model for aqueous solvation free energy prediction of organic molecules in redox flow batteries.

The solvation free energy of organic molecules is a critical parameter in determining emergent properties such as solubility, liquid-phase equilibrium constants, pKa and redox potentials in an organic redox flow battery. In this work, we present a machine learning (ML) model that can learn and predict the aqueous solvation free energy of an organic molecule using the Gaussian process regression method based on a new molecular graph kernel. To investigate the performance of the ML model for electrostatic interaction, the nonpolar interaction contribution of the solvent and the conformational entropy of the solute in the solvation free energy, three data sets with implicit or explicit water solvent models, and contribution of the conformational entropy of the solute are tested. We demonstrate that our ML model can predict the solvation free energy of molecules at chemical accuracy with a mean absolute error of less than 1 kcal mol-1 for subsets of the QM9 dataset and the Freesolv database. To solve the general data scarcity problem for a graph-based ML model, we propose a dimension reduction algorithm based on the distance between molecular graphs, which can be used to examine the diversity of the molecular data set. It provides a promising way to build a minimum training set to improve prediction for certain test sets where the space of molecular structures is predetermined.

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer.

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute-solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (FeiiRuiii) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120-200 cm-1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

Quantification of High-Temperature Transition Al2 O3 and Their Phase Transformations*.

High-temperature treatment of γ-Al2 O3 can lead to a series of polymorphic transformations, including the formation of δ-Al2 O3 and θ-Al2 O3 . Quantification of the microstructure in the range where δ- and θ-Al2 O3 are formed represents a formidable challenge, as both phases accommodate a high degree of structural disorder. In this work, we explore the use of an XRD recursive-stacking formalism for the quantification of high-temperature transition aluminas. We formulate the recursive-stacking methodology for modelling of disorder in δ-Al2 O3 and twinning in θ-Al2 O3 and show that explicitly accounting for the disorder is necessary to reliably model the XRD patterns of high-temperature transition alumina. We also use the recursive stacking approach to study phase transformation during high-temperature (1050 °C) treatment. We show that the two different intergrowth modes of δ-Al2 O3 have different transformation characteristics and that a significant portion of δ-Al2 O3 is stabilized with θ-Al2 O3 even after prolonged high-temperature exposures.

Elucidation of the photoaquation reaction mechanism in ferrous hexacyanide using synchrotron x-rays with sub-pulse-duration sensitivity.

Ligand substitution reactions are common in solvated transition metal complexes, and harnessing them through initiation with light promises interesting practical applications, driving interest in new means of probing their mechanisms. Using a combination of time-resolved x-ray absorption spectroscopy and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations and x-ray absorption near-edge spectroscopy calculations, we elucidate the mechanism of photoaquation in the model system iron(ii) hexacyanide, where UV excitation results in the exchange of a CN- ligand with a water molecule from the solvent. We take advantage of the high flux and stability of synchrotron x-rays to capture high precision x-ray absorption spectra that allow us to overcome the usual limitation of the relatively long x-ray pulses and extract the spectrum of the short-lived intermediate pentacoordinated species. Additionally, we determine its lifetime to be 19 (±5) ps. The QM/MM simulations support our experimental findings and explain the ∼20 ps time scale for aquation as involving interconversion between the square pyramidal (SP) and trigonal bipyramidal pentacoordinated geometries, with aquation being only active in the SP configuration.

IRMOF-74(n)-Mg: a novel catalyst series for hydrogen activation and hydrogenolysis of C-O bonds.

Metal-Organic Frameworks (MOFs) that catalyze hydrogenolysis reactions are rare and there is little understanding of how the MOF, hydrogen, and substrate molecules interact. In this regard, the isoreticular IRMOF-74 series, two of which are known catalysts for hydrogenolysis of aromatic C-O bonds, provides an unusual opportunity for systematic probing of these reactions. The diameter of the 1D open channels can be varied within a common topology owing to the common secondary building unit (SBU) and controllable length of the hydroxy-carboxylate struts. We show that the first four members of the IRMOF-74(Mg) series are inherently catalytic for aromatic C-O bond hydrogenolysis and that the conversion varies non-monotonically with pore size. These catalysts are recyclable and reusable, retaining their crystallinity and framework structure after the hydrogenolysis reaction. The hydrogenolysis conversion of phenylethylphenyl ether (PPE), benzylphenyl ether (BPE), and diphenyl ether (DPE) varies as PPE > BPE > DPE, consistent with the strength of the C-O bond. Counterintuitively, however, the conversion also follows the trend IRMOF-74(III) > IRMOF-74(IV) > IRMOF-74(II) > IRMOF-74(I), with little variation in the corresponding selectivity. DFT calculations suggest the unexpected behavior is due to much stronger ether and phenol binding to the Mg(ii) open metal sites (OMS) of IRMOF-74(III), resulting from a structural distortion that moves the Mg2+ ions toward the interior of the pore. Solid-state 25Mg NMR data indicate that both H2 and ether molecules interact with the Mg(ii) OMS and hydrogen-deuterium exchange reactions show that these MOFs activate dihydrogen bonds. The results suggest that both confinement and the presence of reactive metals are essential for achieving the high catalytic activity, but that subtle variations in pore structure can significantly affect the catalysis. Moreover, they challenge the notion that simply increasing MOF pore size within a constant topology will lead to higher conversions.

Comprehensive Experimental and Computational Spectroscopic Study of Hexacyanoferrate Complexes in Water: From Infrared to X-ray Wavelengths.

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.

Modeling Optical Spectra of Large Organic Systems Using Real-Time Propagation of Semiempirical Effective Hamiltonians.

We present an implementation of a time-dependent semiempirical method (INDO/S) in NWChem using real-time (RT) propagation to address, in principle, the entire spectrum of valence electronic excitations. Adopting this model, we study the UV/vis spectra of medium-sized systems such as P3B2 and f-coronene, and in addition much larger systems such as ubiquitin in the gas phase and the betanin chromophore in the presence of two explicit solvents (water and methanol). RT-INDO/S provides qualitatively and often quantitatively accurate results when compared with RT- TDDFT or experimental spectra. Even though we only consider the INDO/S Hamiltonian in this work, our implementation provides a framework for performing electron dynamics in large systems using semiempirical Hartree-Fock Hamiltonians in general.

Surface Interactions and Confinement of Methane: A High Pressure Magic Angle Spinning NMR and Computational Chemistry Study.

Characterization and modeling of the molecular-level behavior of simple hydrocarbon gases, such as methane, in the presence of both nonporous and nanoporous mineral matrices allows for predictive understanding of important processes in engineered and natural systems. In this study, changes in local electromagnetic environments of the carbon atoms in methane under conditions of high pressure (up to 130 bar) and moderate temperature (up to 346 K) were observed with 13C magic-angle spinning (MAS) NMR spectroscopy while the methane gas was mixed with two model solid substrates: a fumed nonporous, 12 nm particle size silica and a mesoporous silica with 200 nm particle size and 4 nm average pore diameter. Examination of the interactions between methane and the silica systems over temperatures and pressures that include the supercritical regime was allowed by a novel high pressure MAS sample containment system, which provided high resolution spectra collected under in situ conditions. For pure methane, no significant thermal effects were found for the observed 13C chemical shifts at all pressures studied here (28.2, 32.6, 56.4, 65.1, 112.7, and 130.3 bar). However, the 13C chemical shifts of resonances arising from confined methane changed slightly with changes in temperature in mixtures with mesoporous silica. The chemical shift values of 13C nuclides in methane change measurably as a function of pressure both in the pure state and in mixtures with both silica matrices, with a more pronounced shift when meso-porous silica is present. Molecular-level simulations utilizing GCMC, MD, and DFT confirm qualitatively that the experimentally measured changes are attributed to interactions of methane with the hydroxylated silica surfaces as well as densification of methane within nanopores and on pore surfaces.

Theoretical Modeling of (99)Tc NMR Chemical Shifts.

Technetium-99 (Tc) displays a rich chemistry due to its wide range of accessible oxidation states (from -I to +VII) and ability to form coordination compounds. Determination of Tc speciation in complex mixtures is a major challenge, and (99)Tc nuclear magnetic resonance (NMR) spectroscopy is widely used to probe chemical environments of Tc in odd oxidation states. However, interpretation of (99)Tc NMR data is hindered by the lack of reference compounds. Density functional theory (DFT) calculations can help to fill this gap, but to date few computational studies have focused on (99)Tc NMR of compounds and complexes. This work evaluates the effectiveness of both pure generalized gradient approximation and their corresponding hybrid functionals, both with and without the inclusion of scalar relativistic effects, to model the (99)Tc NMR spectra of Tc(I) carbonyl compounds. With the exception of BLYP, which performed exceptionally well overall, hybrid functionals with inclusion of scalar relativistic effects are found to be necessary to accurately calculate (99)Tc NMR spectra. The computational method developed was used to tentatively assign an experimentally observed (99)Tc NMR peak at -1204 ppm to fac-Tc(CO)3(OH)3(2-). This study examines the effectiveness of DFT computations for interpretation of the (99)Tc NMR spectra of Tc(I) coordination compounds in high salt alkaline solutions.

Protein-Mineral Interactions: Molecular Dynamics Simulations Capture Importance of Variations in Mineral Surface Composition and Structure.

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.

Charge retention of soft-landed phosphotungstate Keggin anions on self-assembled monolayers.

Soft landing of mass-selected ions onto surfaces often results in partial loss of charge that may affect the structure and reactivity of deposited species. In this study, Keggin phosphotungstate anions in two selected charge states, PW12O40(3-) (WPOM(3-)) and PW12O40(2-) (WPOM(2-)), were soft-landed onto different self-assembled monolayer (SAM) surfaces and examined using in situ infrared reflection absorption spectroscopy (IRRAS) and density functional theory (DFT) calculations. Partial retention of the 3- charge was observed when WPOM(3-) was soft-landed onto the fluorinated SAM (FSAM), while the charge state distribution was dominated by the 2- charge after both WPOM(3-) and WPOM(2-) were deposited onto a hydrophilic alkylthiol SAM terminated with cationic NH3(+) functional groups (NH3(+)SAM). We found that during the course of the soft landing of WPOM(3-), the relative abundance of WPOM(3-) on FSAM decreased while that of WPOM(2-) increased. We propose that the higher stability of immobilized WPOM(2-) in comparison with WPOM(3-) makes it the preferred charge state of WPOM on both the FSAM and NH3(+)SAM. We also observe weaker binding of WPOM anions to SAMs in comparison with phosphomolybdate ions (MoPOM) reported previously (J. Phys. Chem. C, 2014, 118, 27611-27622). The weaker binding of WPOM to SAMs is attributed to the lower reactivity of WPOM reported in the literature. This study demonstrates that both the charge retention and the reactivity of deposited anionic POM clusters on surfaces are determined by the type of addenda metal atoms in the cluster.

Self-learning kinetic Monte Carlo simulations of Al diffusion in Mg.

Vacancy-mediated diffusion of an Al atom in the pure Mg matrix is studied using the atomistic, on-lattice self-learning kinetic Monte Carlo (SLKMC) method. Activation barriers for vacancy-Mg and vacancy-Al atom exchange processes are calculated on the fly using the climbing image nudged-elastic-band method and binary Mg-Al modified embedded-atom method interatomic potential. Diffusivities of an Al atom obtained from SLKMC simulations show the same behavior as observed in experimental and theoretical studies available in the literature; that is, an Al atom diffuses faster within the basal plane than along the c-axis. Although the effective activation barriers for an Al atom diffusion from SLKMC simulations are close to experimental and theoretical values, the effective prefactors are lower than those obtained from experiments. We present all the possible vacancy-Mg and vacancy-Al atom exchange processes and their activation barriers identified in SLKMC simulations. A simple mapping scheme to map an HCP lattice onto a simple cubic lattice is described, which enables simulation of the HCP lattice using the on-lattice framework. We also present the pattern recognition scheme which is used in SLKMC simulations to identify the local Al atom configuration around a vacancy.

Highly Selective Colorimetric and Luminescence Response of a Square-Planar Platinum(II) Terpyridyl Complex to Aqueous TcO4(-).

Molecular recognition of an aqueous pertechnetate (TcO4(-)) anion is fundamentally challenging partly due to the charge-diffuse nature of this anion, which hampers design of new technologies for its separation and detection. To address this gap, simple salts of transition metal complexes that undergo a distinct spectroscopic change upon exposure to aqueous anions were explored. The Pt(II) complex [Pt(tpy)Br]SbF6 (tpy = 2,2';6',2″-terpyridine) undergoes a dramatic color change and intense luminescence response upon TcO4(-) uptake due to concomitant enhancement of Pt···Pt interactions. The spectroscopic response was highly selective and quantitative for aqueous TcO4(-) among other competing anions. Complementary Raman spectroscopy and microscopy techniques, structural determination, and theoretical methods were employed to elucidate the mechanism of this response at the molecular level.