Search NASA⌕ Search

SEARCH · Search NASA

Results for “Prevalence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

Molecular Determination of Organic Adsorption Sites on Smectite during Fe Redox Processes Using ToF-SIMS Analysis

Turnover of soil organic carbon (SOC) is strongly affected by a balance between mineral protection and microbial degradation. However, the mechanisms controlling the heterogeneous and preferential adsorption of different types of SOC remain elusive, especially for humic substrates (HSs) and microbial carbon (MC), two hotly-debated forms of SOC prevalent in soils. In this work, the heterogeneous adsorption of HSs and MC on a clay mineral (nontronite NAu-2) during microbial-mediated Fe redox cycling was determined with the use of time-of-flight secondary ion mass spectrometry (ToF-SIMS). The results revealed that HSs pre-adsorbed on NAu-2 would partially inhibit the microbial Fe(III) reduction and structural modification of NAu-2, thus retarding the subsequent adsorption of MC. In contrast, NAu-2 without pre-coated HSs was adsorbed by a significant amount of MC from microbial polysaccharides as a result of Fe(III) reduction. This was attributed to the deposition of a thin Al-rich layer on clay surface, which provided active sites for MC adsorption. This study provides direct and detailed molecular evidence for the first time to explain the preferential adsorption of MC over HSs on the surface of clay minerals in the iron redox processes, which could be critical for the preservation of MC in soil. The results also indicate that ToF-SIMS is a unique tool for understanding complex organic-mineral-microbe interactions.

54 ENVIRONMENTAL SCIENCES↗

Elucidating the Role of Single-Atom Pd for Electrocatalytic Hydrodechlorination

In this study, we loaded Pd catalysts onto a reduced graphene oxide (rGO) support in an atomically dispersed fashion [i.e., Pd single-atom catalysts (SACs) on rGO or Pd 1 /rGO] via a facile and scalable synthesis based on anchor-site and photoreduction techniques. The as-synthesized Pd 1 /rGO significantly outperformed the Pd nanoparticle (Pdnano) counterparts in the electrocatalytic hydrodechlorination of chlorinated phenols. Downsizing Pd nano to Pd 1 leads to a substantially higher Pd atomic efficiency (14 times that of Pdnano), remarkably reducing the cost for practical applications. The unique single-atom architecture of Pd 1 additionally affects the desorption energy of the intermediate, suppressing the catalyst poisoning by Cl – , which is a prevalent challenge with Pdnano. Characterization and experimental results demonstrate that the superior performance of Pd 1 /rGO originates from (1) enhanced interfacial electron transfer through Pd–O bonds due to the electronic metal–support interaction and (2) increased atomic H (H*) utilization efficiency by inhibiting H 2 evolution on Pd 1 . Furthermore, this work presents an important example of how the unique geometric and electronic structure of SACs can tune their catalytic performance toward beneficial use in environmental remediation applications.

36 MATERIALS SCIENCE↗

Mercury Reduction, Uptake, and Species Transformation by Freshwater Alga Chlorella vulgaris under Sunlit and Dark Conditions

As a major entry point of mercury (Hg) to aquatic food webs, algae play an important role in taking up and transforming Hg species in aquatic ecosystems. However, little is known how and to what extent Hg reduction, uptake, and species transformations are mediated by algal cells and their exudates, algal organic matter (AOM), under either sunlit or dark conditions. In this study, using Chlorella vulgaris (CV) as one of the most prevalent freshwater model algal species, we show that solar irradiation could enhance the reduction of mercuric Hg(II) to elemental Hg(0) by both CV cells and AOM. AOM reduced more Hg(II) than algal cells themselves due to cell surface adsorption and uptake of Hg(II) inside the cells under solar irradiation. Synchrotron radiation X-ray absorption near-edge spectroscopy (SR-XANES) analyses indicate that sunlight facilitated the transformation of Hg to less bioavailable species, such as β-HgS and Hg-phytochelatins, compared to Hg(Cysteine) 2 -like species formed in algal cells in the dark. These findings highlight important functional roles and potential mechanisms of algae in Hg reduction and immobilization under varying lighting conditions and how these processes may modulate Hg cycling and bioavailability in the aquatic environment.

54 ENVIRONMENTAL SCIENCES↗

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗

Manganese Oxidation States in Volcanic Soils across Annual Rainfall Gradients

Manganese (Mn) exists as Mn(II), Mn(III), or Mn(IV) in soils, and the Mn oxidation state controls the roles of Mn in numerous environmental processes. However, the variations of Mn oxidation states with climate remain unknown. Here we determined the Mn oxidation states in highly weathered bulk volcanic soils (primary minerals free) across two rainfall gradients covering mean annual precipitation (MAP) of 0.25–5 m in the Hawaiian Islands. With increasing MAP, the soil redox conditions generally shifted from oxic to suboxic and to anoxic despite fluctuating at each site; concurrently, the proportions of Mn(IV) and Mn(II) decreased and increased, respectively. Mn(III) was low at both low and high MAP, but accumulated substantially, up to 80% of total Mn, in soils with prevalent suboxic conditions at intermediate MAP. Mn(III) was likely hosted in Mn(III,IV) and iron(III) oxides or complexed with organic matter, and its distribution among these hosts varied with soil redox potentials and soil pH. Soil redox conditions and rainfall-driven leaching jointly controlled exchangeable Mn(II) in soils, with its concentration peaking at intermediate MAP. The Mn redox chemistry was at disequilibrium, with the oxidation states correlating with long-term average soil redox potentials better than with soil pH. The soil redox conditions likely fluctuated between oxic and anoxic conditions more frequently at intermediate than at low and high MAP, creating biogeochemical hot spots where Mn, Fe, and other redox-sensitive elements may be actively cycled.

36 MATERIALS SCIENCE↗

Suspect and Nontarget Screening for Contaminants of Emerging Concern in an Urban Estuary

This study used suspect and nontarget screening with high-resolution mass spectrometry to characterize the occurrence of contaminants of emerging concern (CECs) in the nearshore marine environment of Puget Sound (WA). In total, 87 non-polymeric CECs were identified; those confirmed with reference standards (45) included pharmaceuticals, herbicides, vehicle-related compounds, plasticizers, and flame retardants. Eight polyfluoroalkyl substances were detected; perfluorooctanesulfonic acid (PFOS) concentrations were as high as 72-140 ng/L at one location. Low levels of methamphetamine were detected in 41% of the samples. Transformation products of pesticides were tentatively identified, including two novel transformation products of tebuthiuron. While a hydrodynamic simulation, analytical results, and dilution calculations demonstrated the prevalence of wastewater effluent to nearshore marine environments, the identity and abundance of selected CECs revealed the additional contributions from stormwater and localized urban and industrial sources. For the confirmed CECs, risk quotients were calculated based on concentrations and predicted toxicities, and eight CECs had risk quotients >1. Dilution in the marine estuarine environment lowered the risks of most wastewater-derived CECs, but dilution alone is insufficient to mitigate risks of localized inputs. These findings highlighted the necessity of suspect and nontarget screening and revealed the importance of localized contamination sources in urban marine environments.

Tian, Zhenyu↗

Cis-Divacant Octahedral Fe(II) in a Dimensionally Reduced Family of 2-(Pyridin-2-yl)pyrrolide Complexes

Four-coordinate transition metal complexes can adopt a diverse array of coordination geometries, with square planar and tetrahedral coordination being the most prevalent. Previously, we reported the synthesis of a trinuclear Fe(II) complex, Fe 3 TPM 2 , supported by a three-fold symmetric 2-pyridylpyrrolide ligand (i.e., tris(5-(pyridin-2-yl)-1H-pyrrol-2-yl)methane), that featured a rare cis-divacant octahedral (CDO) geometry at each Fe(II) center. Here, a series of truncated 2-pyridylpyrrolide ligands is described that support mono- and binuclear Fe(II) complexes that also exhibit CDO geometries. Metallation of tetradentate ligand bis(5-(pyridin-2-yl)-1H-pyrrol-2-yl)methane (H 2 BPM) in THF results in a binuclear complex Fe 2 (BPM) 2 (THF) 2 in which both Fe(II) ions are octahedrally coordinated. The coordinated THF solvent ligands are labile: THF dissociation leads to Fe 2 (BPM) 2 , which features five-coordinate Fe(II) ions. The Fe–Fe distance in these binuclear complexes can be elongated by ligand methylation. Metalation of bis(5-(6-methylpyridin-2-yl)-1H-pyrrol-2-yl)methane (H 2 BPM Me ) in THF leads to the formation of four-coordinate, CDO Fe(II) centers in Fe(BPM Me ) 2 . Further ligand truncation affords bidentate ligands 2-(1H-pyrrol-2-yl)pyridine (PyrPyrrH) and 2-methyl-6-(1H-pyrrol-2-yl)pyridine (Pyr Me PyrrH). Metalation of these ligands in THF affords six-coordinate complexes Fe(PyrPyrr) 2 (THF) 2 and Fe(Pyr Me Pyrr) 2 (THF) 2 . Dissociation of labile solvent ligands provides access to fourcoordinate Fe(II) complexes. Ligand disproportionation at Fe(PyrPyrr) 2 results in the formation of Fe(PyrPyrr) 3 and Fe(0). Ligand methylation suppresses this disproportionation and enables isolation of Fe(Pyr Me Pyrr) 2 , which is rigorously CDO. Complete ligand truncation, by separating the 2-pyridylpyrrolide ligands into the constituent monodentate pyridyl and pyrrolide donors, affords Fe(Pyr) 2 (Pyrr) 2 in which the Fe(II) is tetrahedrally coordinated. Computational analysis indicates that the potential energy surface that dictates the coordination geometry in this family of four-coordinate complexes is fairly flat in the vicinity of CDO coordination. Furthermore, these synthetic studies provide the structural basis to explore the implications of CDO geometry on Fe-catalyzed reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Oxygen Defect Configurations in Single-Phase UO 2.15

Hyperstoichiometric UO 2.15 was characterized by neutron total scattering at high temperature in the single-phase UO 2+x region of the U/O phase diagram. The diffraction data confirmed a single-phase fluorite structure at high temperature. Analysis of the short-range data showed that the same structural model does not fit the pair distribution functions well. Instead, structural models containing specific configurations of oxygen defect clusters best represent the local atomic arrangement. In conclusion, prevalent defect clusters previously proposed were fit to the experimental data, and moderately distorted oxygen cuboctahedra hypothesized by recent molecular dynamics simulations fit the data most accurately.

neutron pair-distribution function analysis↗

A Computational Assessment of Actinide Dioxide Cations AnO 2 2+ for An = U to Lr: The Limited Stability Range of the Hexavalent Actinyl Moiety, [O=An=O] 2+

The isolated gas-phase actinide dioxide dications, AnO 2 2+ , were evaluated by DFT and coupled cluster CCSD(T) calculations for 12 actinides, An = U–Lr. CASSCF calculations were used to define the orbitals for the CCSD(T) calculations. The characteristic linear [O$=$An$=$O] 2+ hexavalent actinyl(VI) was found to be the lowest energy structure for An = U, Np, and Pu, which also form stable actinyl(VI) species in solution and possibly for Am when spin–orbit effects are included. For Am, there is a divalent [An II (O 2 )] 2+ structure where the dioxygen is an end-on physisorbed η 1 - 3 O 2 2 kcal/mol above the actinyl when spin–orbit effects are included which lower the energy of the actinyl structure. For An = Cm, Bk, and Lr, the lowest energy structure is trivalent [An III (O 2 – )] 2+ where the dioxygen is a side-on superoxide, η 2 -O 2 – . For Cm, the actinyl is close in energy to the ground state when spin–orbit effects are included. For An = Cf, Es, Fm, Md, and No, the lowest energy structure is divalent [An II (O 2 )] 2+ where the dioxygen is an end-on physisorbed η 1 - 3 O 2 . The relative energies suggest that curyl(VI) and berkelyl(VI), like well-known americyl(VI), might be stabilized by coordinating ligands in condensed phases. The results further indicate that for californyl and beyond, the actinyl(VI) moieties will probably be elusive even using strong donor ligands. The prevalence of low oxidation states (OSs) An(II) and An(III) for transplutonium actinides reflects stabilization of the 5f orbitals and validates established trends, including the remarkably high stability of divalent No. Finally, bond distances and other parameters suggest maximum bond covalency around plutonyl(VI), with a particularly substantial decrease in bond strength between americyl(VI) and curyl(VI).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantitative Account of the Bonding Properties of a Rubredoxin Model Complex [Fe(SCH3)4]q, q = -2, -1, +2, +3

Iron-sulfur clusters play important roles in biology as parts of electron transfer chains and catalytic cofactors. Here, we report a detailed computational analysis of a structural model of the simplest natural iron-sulfur cluster of rubredoxin and its cationic counterparts. Specifically, we report results for the ground and low-lying electronically excited states of the complex [Fe(SCH3)4]2-/1-/2+/3+, using Multi-Reference (CASSCF, MRCISD), and Coupled Cluster [CCSD(T)] methodology in order to provide accurate adiabatic reduction energies, dissociation energies and insights into the bonding analysis. The nature of the Fe-S chemical bond and the magnitude of the ionization potentials in the anionic and cationic [Fe(SCH3)4] complexes offer a physical rationale for the relative stabilization, structure and speciation of these complexes. Anionic and cationic complexes present different types of chemical bonds: prevalently ionic in [Fe(SCH3)4]2-/1- complexes and covalent in [Fe(SCH3)4]2+/3+ complexes. The ionic bonds result in an energy gain for the transition [Fe(SCH3)4]2-®[Fe(SCH3)4]- (i.e., FeII®FeIII) of 1.5 eV, while the covalent bonds result in an energy loss for the transition [Fe(SCH3)4]2+®[Fe(SCH3)4]3+ of 16.6 eV, almost half of the IP of Fe2+. The ionic vs covalent bond character influences the Fe-S bond strength and length, i.e., ionic Fe-S bonds are longer than covalent ones by about 0.2 Å (for FeII) and 0.04 Å (for FeII). Finally, the average Fe-S heterolytic bond strength is 6.7 eV (FeII) and 14.6 (FeIII) eV at the RCCSD(T) level of theory.

Tzeli, Demeter↗

Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations

Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. Here, in this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing the Limitations of Self-Interaction-Corrected Functionals for Describing the Hydrated Electron

Simulating the hydrated electron using density functional theory is challenging due to the prevalence of self-interaction error in standard functionals. Hybrid functionals like PBEh(40) can reasonably describe the chemistry of an excess electron in water and partially mitigate self-interaction error by incorporating exact Hartree–Fock exchange, but they are computationally expensive making them impractical for large-scale and long-time ab initio molecular dynamics simulations. Explicit self-interaction correction schemes that are applied on an orbital-by-orbital basis offer a potential alternative when the correction is limited to the singly occupied molecular orbital obtained with a generalized gradient approximation functional. Here, we examine whether the Perdew–Zunger self-interaction correction scheme applied to the revPBE functional can provide a computationally efficient and physically sensible alternative to PBEh(40) for the hydrated electron. We find that functionals incorporating a self-interaction correction scheme should be viewed with caution when applied to the hydrated electron and its reactivity. Furthermore, we show that it is critical to consider extensive sampling and diverse chemical environments when validating their performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of Pyridopyrimidinones that Selectively Inhibit the H1047R PI3Kα Mutant Protein

The H1047R mutation of PIK3CA is highly prevalent in breast cancers and other solid tumors. Selectively targeting PI3Kα H1047R over PI3Kα WT is crucial due to the role that PI3Kα WT plays in normal cellular processes, including glucose homeostasis. Currently, only one PI3Kα H1047R -selective inhibitor has progressed into clinical trials, while three pan mutant (H1047R, H1047L, H1047Y, E542K, and E545K) selective PI3Kα inhibitors have also reached the clinical stage. Herein, we report the design and discovery of a series of pyridopyrimidinones that inhibit PI3Kα H1047R with high selectivity over PI3Kα WT , resulting in the discovery of compound 17. When dosed in the HCC1954 tumor model in mice, 17 provided tumor regressions and a clear pharmacodynamic response. X-ray cocrystal structures from several PI3Kα inhibitors were obtained, revealing three distinct binding modes within PI3Kα H1047R including a previously reported cryptic pocket in the C-terminus of the kinase domain wherein we observe a ligand-induced interaction with Arg1047.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature Dependence of Short and Intermediate Range Order in Molten MgCl 2 and Its Mixture with KCl

This article addresses the non-Debye-Waller temperature behavior in the intermediate range order for molten MgCl 2 and its mixtures with KCl from a theory, Molecular Dynamics and experimental X-ray scattering perspective and puts these findings in the context of discussions and controversies extending at least four decades. We find that these liquids are defined by two structural motifs. The first motif is associated with chains of positive-negative charge alternation; the second motif, which results in a prepeak in the structure function S(q), is associated with the interaction of Mg 2+ and Cl - ions that do not belong to the same charge alternation chain or aggregate. Our complementary X-ray scattering and computational results provide a quantitative explanation for the increase in intensity of the prepeak with temperature as opposed to the behavior of other peaks following normal Debye-Waller behavior. We note temperature has opposite effects on the prevalence of each of the two structural motifs, and the enhancement of one pattern appears to be at the detriment of the other. Whereas the intensity in S(q) associated with the charge alternation motif is diminished at higher temperature the opposite is true for the prepeak associated with the intermediate range order due to the second structural motif.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗