Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “enhanced sampling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗
DOE OSTI2024

CO 2 Uptake and Stability Enhancement in Vinyltrimethoxysilane‐Treated SBA‐15 Solid Amine‐Based Sorbents

Abstract Silica‐supported amine absorbents, including materials produced by tethering aminosilanes or infusion of poly(ethyleneimine), represent a promising class of materials for CO 2 capture applications, including direct air and point source capture. Various silica surface treatments and functionalization strategies are explored to enhance stability and CO 2 uptake in amine‐based solid sorbent systems. Here, the synthesis and characterization of novel vinyltrimethoxysilane‐treated Santa Barbara Amorphous‐15 (SBA‐15) supports and the corresponding enhancement in CO 2 uptake compared to various SBA‐15‐based control supports are presented. The relationship between CO 2 diffusion and amine efficiency in these systems is explored using a previously reported kinetic model. The synthesized materials are characterized with CO 2 and H 2 O isotherms, diffuse reflectance infrared Fourier transform spectroscopy, 1 H T 1 – T 2 relaxation correlation NMR, and rapid thermal cycling experiments. The novel support materials are shown to enable high amine efficiencies, approaching a fourfold improvement over standard SBA‐15‐supported amines, while simultaneously exhibiting excellent stability when cycled rapidly under humid conditions. As the poly(ethyleneimine) loadings are held constant across the various samples, enhancements in CO 2 uptake are attributed to differences in the way the poly(ethyleneimine) interacts with the support surface.

Vallace, Anthony↗
DOE OSTI2019

Porting Adaptive Ensemble Molecular Dynamics Workflows to the Summit Supercomputer

Molecular dynamics (MD) simulations must take very small (femtosecond) integration steps in simulation-time to avoid numerical errors. Efficient use of parallel programming models and accelerators in state-of-the art MD programs now is pushing Moore’s limit for time-per-MD step. As a result, directly simulating timescales beyond milliseconds will not be attainable directly, even at exascale. However, concepts from statistical physics can be used to combine many parallel simulations to provide information about longer timescales and to adequately sample the simulation space, while preserving details about the dynamics of the system. Implementing such an approach requires a workflow program that allows adaptable steering of task assignments based on extensive statistical analysis of intermediate results. Here we report the implementation of such an adaptable workflow program to drive simulations on the Summit IBM Power System AC922, a pre-exascale supercomputer at the Oak Ridge Leadership Computing Facility (OLCF). We compare to experiences on Titan, Summit’s predecessor, report the performance of the workflow and its components, and describe the porting process. We find that using a workflow program managed by a Mongo database can provide the fault tolerance, scalable performance, task dispatch rate, and reconfigurability required for robust and portable implementation of ensemble simulations such as are used in enhanced-sampling molecular dynamics. This type of workflow generator can also be used to provide adaptive steering of ensemble simulations for other applications in addition to MD.

Ossyra, John R.↗
DOE OSTI2024

High Entropy Protected Sharp Magnetic Transitions in Highly Disordered Spinel Ferrites

How disorder affects magnetic ordering is always an intriguing question, and it becomes even more interesting in the recently rising high entropy oxides due to the extremely high disorder density. However, due to the lack of high-quality single crystal samples, the strong compositional disorder effect on magnetic transition has not been deeply investigated. In this work, we have successfully synthesized high-quality single crystalline high entropy spinel ferrites (Mg 0.2 Mn 0.2 Fe 0.2 Co 0.2 Ni 0.2 ) x Fe 3-x O 4 . Our findings from high-temperature magnetization and neutron diffraction experiments showed ferrimagnetic transitions at 748, 694, and 674 K for x values of 1, 1.5, and 1.8, respectively. Notably, the magnetic transition almost showed no broadening for x values of 1 and 1.5, compared to Fe 3 O 4 . Extended X-ray absorption fine structure measurements provided insights into the elemental distribution among the octahedral and tetrahedral sites. The random distribution of elements across these sites reduced the formation of local clusters and short-range orders, enhancing sample homogeneity and preserving the sharpness of the magnetic transition, despite bond length variation. Here, our study not only marks the first successful synthesis of an HEO bulk single crystal exhibiting long-range magnetic order but also sheds light on the interaction between high configurational entropy and magnetic orderings. This opens new avenues for future research and applications of magnetic high entropy oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The Role of Nonequilibrium Solvent Effects in Enhancing Direct CO 2 Capture at the Air–Aqueous Amino Acid Interface

Direct air capture (DAC) technologies are limited by the poor understanding of the dynamic role of interfaces in modulating the chemisorption of CO 2 from air into solutions. While the reactivity of aqueous amine-based solvents in the bulk environment is strongly inhibited by nonequilibrium solvent effects, promoting DAC at interfaces posits a possibility to reduce the coupling with the solvent and significantly accelerate DAC. Building on an experimentally proven concept to bring an anionic glycine absorbent to the interface through ion-pairing interactions with a positively charged surfactant, we establish the fundamental time scales for key elementary steps involved in DAC with rate theory and enhanced-sampling ab initio molecular dynamics simulations. Here, we elucidate the mechanism by which water influences the free energy barriers and dynamical crossing-recrossing of those barriers, affecting the reaction rates. We find that water reorganizes to partially dehydrate [-NH 2 ], facilitating S N 2-based CO 2 conversion to a zwitterion, which then releases a proton via overhydration of [-NH 2 ]. The low-density interfacial water favors dehydration over overhydration, leading to a comparatively higher barrier (slower kinetics) for proton release. The barrier-recrossing events neutralize this effect, letting both steps occur at the same time scale (sub-microseconds) and making the overall DAC process faster at the interface than in the bulk water. Such an understanding of environment-sensitive solvent effects on the reaction kinetics will help design tailored interfaces for enhanced CO 2 capture kinetics via control of solvation and ion paring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Improved sample utilization in thermal ionization mass spectrometry isotope ratio measurements: refined development of porous ion emitters for nuclear forensic applications

The precise and accurate determination of isotopic composition in nuclear forensic samples is vital for assessing origin, intended use and process history. Thermal ionization mass spectrometry (TIMS) is widely accepted as the gold standard for high performance isotopic measurements and has long served as the workhorse in the isotopic ratio determination of nuclear materials. Nuclear forensic and safeguard specialists have relied heavily on such methods for both routine and atypical e orts. Despite widespread use, TIMS methods for the assay of actinide systems continue to be hindered by poor ionization e ciency, often less than tenths of a percent; the majority of a sample is not measured. This represents a growing challenge in addressing nextgeneration nuclear detection needs by limiting the ability to analyze ultratrace quantities of high priority elements that could potentially provide critical nuclear forensic signatures. Porous ion emitter (PIE) thermal ion sources were developed in response to the growing need for new TIMS ion source strategies for improved ionization e ciency, PIEs have proven to be simple to implement, straightforward approach to boosting ion yield. This work serves to expand the use of PIE techniques for the analysis of trace quantities of plutonium and americium. PIEs exhibited superior plutonium and americium ion yields when compared to direct lament loading and the resin bead technique, one of the most e cient methods for actinide analysis, at similar mass loading levels. Initial attempts at altering PIE composition for the analysis of plutonium proved to enhance sample utilization even further. Preliminary investigations of the instrumental fractionation behavior of plutonium and uranium analyzed via PIE methods were conducted. Data collected during these initial trial indicate that PIEs fractionate in a consistent, reproducible manner; a necessity for high precision isotope ratio measurements. Ultimately, PIEs methods were applied for the age determination of various uranium isotopic standards. PIEs did not exhibit signi cant advantages for the determination of model ages when compared to traditional laments; however, this trial was able to provide valuable insight for guiding future investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Application of the ESMACS Binding Free Energy Protocol to a Multi-Binding Site Lactate Dehydogenase A Ligand Dataset

Over the past two decades, the use of fragment-based lead generation has become a common, mature approach to identify tractable starting points in chemical space for the drug discovery process. This approach naturally involves the study of the binding properties of highly heterogeneous ligands. Such datasets challenge computational techniques to provide comparable binding free energy estimates from different binding modes. The performance of a range of statistically robust ensemble-based binding free energy calculation protocols, called ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent), is evaluated. Ligands designed to target two binding pockets in the lactate dehydogenase, a target protein, which vary in size, charge, and binding mode, are studied. When compared to experimental results, excellent statistical rankings are obtained across this highly diverse set of ligands. In addition, three approaches to account for entropic contributions are investigated: 1) normal mode analysis, 2) weighted solvent accessible surface area (WSAS), and 3) variational entropy. Normal mode analysis and WSAS correlate strongly with each other—although the latter is computationally far cheaper—but do not improve rankings. Variational entropy corrects exaggerated discrimination of ligands bound in different pockets but creates three outliers which reduce the quality of the overall ranking.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Molecular Insights Into the Ionic Assembly of Poly-Galacturonic Acid Oligomers - Impact of Charge, Ionic Radius, and Polymer Functionalization

Pectin, a major class of matrix polysaccharides present in plant cell walls (PCW), contains widespread anionic saccharides that cross-link in the presence of cations. It modulates important functions such as cell-cell adhesion and determines the PCW's biomechanical properties. It is known that mono-, di-, and tri-valent cations facilitate cross-linking; however, significant knowledge gaps remain in understanding the structure and mechanism of pectin cross-linking. In this study, replica-exchange molecular dynamics (REMD) simulations were employed to elucidate the role of ionic charge, ionic radii, and functional groups on the cross-linking of homogalacturonan (HG), the most abundant pectin molecule. Our enhanced sampling approach in fully solvated environments suggests more effective cross-linking with higher-valent and smaller ions, and that the "zipper" conformation is more favorable than the prevalent "egg-box" conformation. These findings advance our fundamental understanding of pectin matrix structure in PCWs and provide a solid foundation to probe structure-property relationships in pectic polysaccharides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

The Effect of SiC Addition on Microstructure and Mechanical Properties of Gas Tungsten Arc-Welded Ti-6Al-4V Alloy

In the present study, SiC powder was added in situ for bead-on-plate weldments and fusion butt joints to modify the microstructure and the corresponding mechanical properties. Cross sections of bead-on-plate welds demonstrated that the SiC addition resulted in significant changes of the titanium microstructure. Specifically, a fine lamellar structure of α-phase with a small amount of intergranular β-titanium within the fusion zone could be established. Moreover, the size of prior β-grains was significantly reduced. After welding, no intermetallic silicides were observed. In the fusion zone of the SiC-enhanced samples, a small increase in hardness resulted, and the samples showed an increase in tensile strength, yield strength and elongation.

36 MATERIALS SCIENCE↗
DOE OSTI2013

miR-206 is down-regulated in breast cancer and inhibits cell proliferation through the up-regulation of cyclinD2

Highlights: ► miR-206 was downexpressed in tumor samples compared with matched normal samples. ► Enhanced expression of miR-206 could inhibit breast cancer growth in vitro. ► Luciferase confirmed miR-206 functions as an anti-oncogene by targeting cyclinD2. ► A reverse correlation between miR-206 and cyclinD2 in breast cancer was found. -- Abstract: MicroRNAs act as important gene regulators in human genomes, and their aberrant expression is linked to many malignancies. Aberrant expression of miR-206 has been frequently reported in cancer studies; however, the role and mechanism of its function in breast cancer remains unclear. Quantitative real-time PCR was performed to detect the relative expression levels of miR-206 in breast cancer and normal breast tissues. Lower expression of miR-206 in breast cancer tissues was associated with larger tumour size and a more advanced clinical stage. Further in vitro observations showed that the enforced expression of miR-206 in MCF-7 breast cancer cells inhibited cell growth by blocking the G1/S transition and suppressed cell proliferation and colony formation, implying that miR-206 functions as a tumour suppressor in the progression of breast cancer. Interestingly, Luciferase assays first revealed that miR-206 inhibited cyclinD2 expression by targeting two binding sites in the 3′-untranslated region of cyclinD2 mRNA. qRT-PCR and Western blot assays verified that miR-206 reduced cyclinD2 expression at both the mRNA and protein levels. A reverse correlation between miR-206 and cyclinD2 expression was noted in breast cancer tissues. Altogether, our results identify a crucial tumour suppressive role of miR-206 in the progression of breast cancer, at least partly via up-regulation of the expression of cyclinD2, and suggest that miR-206 might be a candidate prognostic predictor or an anticancer therapeutic target for breast cancer patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Effects of pH on an IDP conformational ensemble explored by molecular dynamics simulation

The conformational ensemble of intrinsically disordered proteins, such as α-synuclein, are responsible for their function and malfunction. Misfolding of α-synuclein can lead to neurodegenerative diseases, and the ability to study their conformations and those of other intrinsically disordered proteins under varying physiological conditions can be crucial to understanding and preventing pathologies. In contrast to well-folded peptides, a consensus feature of IDPs is their low hydropathy and high charge, which makes their conformations sensitive to pH perturbation. We examine a prominent member of this subset of IDPs, α-synuclein, using a divide-and-conquer scheme that provides enhanced sampling of IDP structural ensembles. We constructed conformational ensembles of α-synuclein under neutral (pH ~ 7) and low (pH ~ 3) pH conditions and compared our results with available information obtained from smFRET, SAXS, and NMR studies. Specifically, α-synuclein has been found to in a more compact state at low pH conditions and the structural changes observed are consistent with those from experiments. We also characterize the conformational and dynamic differences between these ensembles and discussed the implication on promoting pathogenic fibril formation. We find that under low pH conditions, neutralization of negatively charged residues leads to compaction of the C-terminal portion of α-synuclein while internal reorganization allows α-synuclein to maintain its overall end-to-end distance. Here, we also observe different levels of intra-protein interaction between three regions of α-synuclein at varying pH and a shift towards more hydrophilic interactions with decreasing pH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Influence of artificial aging on corrosion of abraded Al-Zn-Mg-Cu alloys

Surface abrasion has been shown to create altered surface layers (ASLs) that are hundreds of nm thick on all types of Al alloys. Such ASLs have lower corrosion resistance than the underlying substrate of Al-Zn-Mg-Cu alloys. Here, we demonstrate how heat treatments affect the ASL microstructure on abraded Al-Zn-Mg-Cu alloys using transmission electron microscopy and atom probe tomography, and also how the corrosion properties change using electrochemical polarization. Additionally, high temperature treatments on the abraded bulk samples enhance η phase precipitation in the ASL, leading to a decreased Zn content in the ASL solid solution, which ennobles the ASL breakdown potential.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗
DOE OSTI2023

Disordered interfaces of alkaline aluminate salt hydrates provide glimpses of Al 3+ coordination changes

Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within a highly disordered quasi-two-dimensional phase at the solid-solution interface. Experiments and simulations: In this work, in situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase. Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH)$_4^–$) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ion linkages. The solution phase Al(OH)$_4^–$ ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (~30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Structural and magnetic characteristics of PVA/CoFe{sub 2}O{sub 4} nano-composites prepared via mechanical alloying method

Highlights: • Single phase CoFe{sub 2}O{sub 4} nano-particles synthesized in one step by mechanical alloying. • PVA/CoFe{sub 2}O{sub 4} magnetic nano-composites were fabricated via mechanical milling. • FTIR confirmed the interaction between PVA and magnetic CoFe{sub 2}O{sub 4} particles. • Increasing in milling time and PVA amount led to well dispersion of CoFe{sub 2}O{sub 4}. - Abstract: In this research, polyvinyl alcohol/cobalt ferrite nano-composites were successfully synthesized employing a two-step procedure: the spherical single-phase cobalt ferrite of 20 ± 4 nm mean particle size was synthesized via mechanical alloying method and then embedded into polymer matrix by intensive milling. The results revealed that increase in polyvinyl alcohol content and milling time causes cobalt ferrite particles disperse more homogeneously in polymer matrix, while the mean particle size and shape of cobalt ferrite have not been significantly affected. Transmission electron microscope images indicated that polyvinyl alcohol chains have surrounded the cobalt ferrite nano-particles; also, the interaction between polymer and cobalt ferrite particles in nano-composite samples was confirmed. Magnetic properties evaluation showed that saturation magnetization, coercivity and anisotropy constant values decreased in nano-composite samples compared to pure cobalt ferrite. However, the coercivity values of related nano-composite samples enhanced by increasing PVA amount due to domain wall mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Expanded Coverage of Phytocompounds by Mass Spectrometry Imaging Using On-Tissue Chemical Derivatization by 4-APEBA

Probing the entirety of any species metabolome is an analytical grand challenge, especially at a cellular scale. Where spatial metabolomics, completed primarily by matrix-assisted laser desorption/ionization (MALDI), has limited molecular coverage for several reasons. To expand the scope of spatial metabolomics, we developed an on-tissue chemical derivatization (OTCD) workflow using 4-APEBA for confident identification of several dozen elusive phytocompounds, including several phytohormones, which have various roles within stress responses and cellular communication. Superiority of 4-APEBA is established in comparison to other derivatization agents with (1) broad specificity towards carbonyls, (2) low background, and (3) introduction of bromine isotopes, where the latter two facilitate confident bioinformatics. In conclusion, the outlined workflow trailblazes a path towards spatial hormonomics within plant samples, enhancing detection of carboxylates, aldehydes, ketones, and plausibly phenols.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Peptoid Backbone Flexibilility Dictates Its Interaction with Water and Surfaces: A Molecular Dynamics Investigation

Peptoids are peptide-mimetic biopolymers that are easy-to-synthesize and adaptable for use in drugs, chemical scaffolds, and coatings. However, there is insufficient information about their structural preferences and interactions with the environment in various applications. We conducted a study to understand the fundamental differences between peptides and peptoids using molecular dynamics simulations with semi-empirical (PM6) and empirical (AMBER) potentials, in conjunction with metadynamics enhanced sampling. From studies of single molecules in water and on surfaces, we found that sarcosine (model peptoid) is much more flexible than alanine (model peptide) in different environments. However, the sarcosine and alanine interact similarly with a hydrophobic or a hydrophilic. Finally, this study highlights the conformational landscape of peptoids and the dominant interactions that drive peptoids towards these conformations. ACKNOWLEDGMENT: MD simulations and manuscript preparation were supported by the MS3 (Materials Synthesis and Simulation Across Scales) Initiative at Pacific Northwest National Laboratory (PNNL), a multi-program national laboratory operated by Battelle for the U.S. Department of Energy. CJM was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences Division of Chemical Sciences, Geosciences, and Biosciences. MDB was supported by the US Department of Energy, Office of Basic Energy Sciences, Biomolecular Materials Program at PNNL. Computing resources were generously allocated by University of Washington's IT department and PNNL's Institutional Computing program. The authors greatly acknowledge conversations with Dr. Kayla Sprenger, Josh Smith, and Dr. Yeneneh Yimer.

alanine↗
DOE OSTI2023

Data-Driven Discovery of Linear Molecular Probes with Optimal Selective Affinity for PFAS in Water

Approaches to tackle the wide and growing variety of highly persistent per- and polyfluoroalkyl substances (PFAS) are of pressing global need because of their detrimental human health effects, such as cancer, birth defects, and hormone imbalance. Sensitive, selective, and easy-to-use real-time sensors to monitor and detect PFAS and sorbents to extract them are critical to meeting government-mandated environmental concentrations. In this work, we combine all-atom molecular dynamics simulations, enhanced sampling, deep representational learning, and Bayesian optimization to perform high-throughput virtual screening for highly sensitive and selective molecular probes. Our molecular design space consists of 3850 linear hydrocarbon chains with varying degrees of halogenation with and without amine- and phosphine-based headgroups. By employing a data-driven search process, we efficiently explore the molecular design space to optimize the sensitivity to perfluorooctanesulfonic acid (PFOS) as a prototypical PFAS analyte and selectivity relative to a sodium dodecyl sulfate (SDS) interferent. We calculate 504 Gibbs free energies of probe-analyte and probe-interferent interactions and identify probes with PFOS association free energies of up to (-ΔG PFOS ) = 9.8 ± 0.2 kJ/mol and selectivities relative to SDS of (-ΔΔG PFOS–SDS ) = 3.1 ± 1.5 kJ/mol. A C 11 Br 23 P(CH 3 ) 2 probe containing 11 backbone brominated carbons and a tertiary phosphine headgroup possesses the most sensitive binding constant to PFOS within the defined search space of K b PFOS = 177.4 ± 12.7, and a semibrominated probe C 5 H 11 C 7 Br 14 N(CH 3 ) 2 containing 12 backbone carbons and a tertiary amine headgroup possesses the highest selectivity relative to SDS of K b PFOS /K b SDS = 4.6 ± 1.7. A retrospective analysis of our data to extract interpretable design rules reveals that the sensitivity of linear hydrogenated probes increases by approximately 1 kJ/mol per C–C bond. The addition or removal of halogen atoms and amine or phosphine headgroups produces nonmonotonic changes in both sensitivity and selectivity with changes to the sensitivity of up to 2.5 kJ/mol. Finally, this work places empirical limitations on the performance of a wide range of linear probes for PFOS detection and offers a generic strategy for high-throughput computational screening to promote selective and sensitive binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗