Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 415 records · Page 23

Evaluating the Transport Mechanism for Lignocellulosic Inhibitors in Zymomonas mobilis

The conversion of lignocellulosic biomass to biofuels and bioproducts is limited by small inhibitory molecules generated during the pretreatment process. The diverse compounds within the cell wall hydrolysate including aldehydes, carboxylic acids, phenolics, and alcohols inhibit microbial growth by disrupting microbial membranes. Although Zymomonas mobilis combines high ethanol productivity with hopanoid-rich membranes that confer solvent tolerance, hydrolysate compounds still inhibit microbial growth. One proposed mitigation strategy is to pump these molecules out through active transport processes, reducing direct toxicity. However, if passive permeation rates are high, active transport would create a futile cycle where the exported molecules diffuse back in, creating a net drag on fitness. To investigate this hypothesis, we measure the passive permeation for thirty-three lignocellulose derived compounds using atomistic molecular dynamics simulations with enhanced sampling. Permeability coefficients determined spanned more than eight orders of magnitude driven largely by the balance between hydrophobicity and polarity. Hydrophobic inhibitors such as aromatic acids, aldehydes, and short chain alcohols readily partitioned into the membrane and showed high permeability consistent with rapid passive influx and strong toxicity. Most polar molecules have lower but still appreciable permeation rates, with only the most polar sugars having low permeabilities. By integrating these predictions together with hydrolysate composition information, we find that most of these molecules are so permeable that futile cycles where passive permeation overwhelms efflux is highly likely. These findings provide a predictive framework for strain engineering, highlighting that, for highly permeable inhibitors, reducing extracellular inhibitor concentrations through hydrolysate detoxification, pretreatment optimization, or dilution may be necessary to mitigate futile cycling and improve microbial productivity in industrial biorefineries.

Zymomonas mobilis↗

Autonomous reinforcement learning agent for stretchable kirigami design of 2D materials

Abstract Mechanical behavior of 2D materials such as MoS 2 can be tuned by the ancient art of kirigami. Experiments and atomistic simulations show that 2D materials can be stretched more than 50% by strategic insertion of cuts. However, designing kirigami structures with desired mechanical properties is highly sensitive to the pattern and location of kirigami cuts. We use reinforcement learning (RL) to generate a wide range of highly stretchable MoS 2 kirigami structures. The RL agent is trained by a small fraction (1.45%) of molecular dynamics simulation data, randomly sampled from a search space of over 4 million candidates for MoS 2 kirigami structures with 6 cuts. After training, the RL agent not only proposes 6-cut kirigami structures that have stretchability above 45%, but also gains mechanistic insight to propose highly stretchable (above 40%) kirigami structures consisting of 8 and 10 cuts from a search space of billion candidates as zero-shot predictions.

36 MATERIALS SCIENCE↗

Probing flow-induced nanostructure of complex fluids in arbitrary 2D flows using a fluidic four-roll mill (FFoRM)

Abstract Engineering flow processes to direct the microscopic structure of soft materials represents a growing area of materials research. In situ small-angle neutron scattering under flow (flow-SANS) is an attractive probe of fluid microstructure under simulated processing conditions, but current capabilities require many different sample environments to fully interrogate the deformations a fluid experiences in a realistic processing flow. Inspired by recent advances in microfluidics, we present a fluidic four-roll mill (FFoRM) capable of producing tunable 2D flow fields for in situ SANS measurements, that is intended to allow characterization of complex fluid nanostructure under arbitrary complex flows within a single sample environment. Computational fluid dynamics simulations are used to design a FFoRM that produces spatially homogeneous and sufficiently strong deformation fields. Particle tracking velocimetry experiments are then used to characterize the flows produced in the FFoRM for several classes of non-Newtonian fluids. Finally, a putative FFoRM-SANS workflow is demonstrated and validated through the characterization of flow-induced orientation in a semi-dilute cellulose nanocrystal dispersion under a range of 2D deformations. These novel experiments confirm that, for steady state straining flows at moderate strain rates, the nanocrystals orient along the principal strain-rate axis, in agreement with theories for rigid, rod-like Brownian particles in a homogeneous flow.

36 MATERIALS SCIENCE↗

A bond-order potential for the Al–Cu–H ternary system

Al-Based Al–Cu alloys have a very high strength to density ratio, and are therefore important materials for transportation systems including vehicles and aircrafts. These alloys also appear to have a high resistance to hydrogen embrittlement, and as a result, are being explored for hydrogen related applications. To enable fundamental studies of mechanical behavior of Al–Cu alloys under hydrogen environments, we have developed an Al–Cu–H bond-order potential according to the formalism implemented in the molecular dynamics code LAMMPS. Our potential not only fits well to properties of a variety of elemental and compound configurations (with coordination varying from 1 to 12) including small clusters, bulk lattices, defects, and surfaces, but also passes stringent molecular dynamics simulation tests that sample chaotic configurations. Careful studies verified that this Al–Cu–H potential predicts structural property trends close to experimental results and quantum-mechanical calculations; in addition, it properly captures Al–Cu, Al–H, and Cu–H phase diagrams and enables simulations of H 2 dissociation, chemisorption, and absorption on Al–Cu surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Selectivity of nitrate and chloride ions in microporous carbons: the role of anisotropic hydration and applied potentials

Understanding ion transport in porous carbons is critical for a wide range of technologies, including supercapacitors and capacitive deionization for water desalination, yet many details remain poorly understood. For instance, an atomistic understanding of how ion selectivity is influenced by the molecular shape of ions, morphology of the micropores and applied voltages is largely lacking. In this work, we combined molecular dynamics simulations with enhanced sampling methods to elucidate the mechanism of nitrate and chloride selectivity in subnanometer graphene slit-pores. We demonstrate that nitrate is preferentially adsorbed over chloride in the slit-like micropores. This preferential adsorption was found to stem from the weaker hydration energy and unique anisotropy of the ion solvation of nitrate. Beside the effects of ion dehydration, we found that applied potential plays an important role in determining the ion selectivity, leading to a lower selectivity of nitrate over chloride at a high applied potential. We conclude that the measured ion selectivity results from a complex interplay between voltage, confinement, and specific ion effects-including ion shape and local hydration structure.

36 MATERIALS SCIENCE↗

Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication

Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.

36 MATERIALS SCIENCE↗

Solvent effects on extractant conformational energetics in liquid–liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid–liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant–metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. Here, we calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plasma flow reactor for steady state monitoring of physical and chemical processes at high temperatures

Here, we present the development of a steady state plasma flow reactor to investigate gas phase physical and chemical processes that occur at high temperature (1000 < T < 5000 K) and atmospheric pressure. The reactor consists of a glass tube that is attached to an inductively coupled argon plasma generator via an adaptor (ring flow injector). We have modeled the system using computational fluid dynamics simulations that are bounded by measured temperatures. In situ line-of-sight optical emission and absorption spectroscopy have been used to determine the structures and concentrations of molecules formed during rapid cooling of reactants after they pass through the plasma. Emission spectroscopy also enables us to determine the temperatures at which these dynamic processes occur. A sample collection probe inserted from the open end of the reactor is used to collect condensed materials and analyze them ex situ using electron microscopy. The preliminary results of two separate investigations involving the condensation of metal oxides and chemical kinetics of high-temperature gas reactions are discussed.

36 MATERIALS SCIENCE↗

Time filtering of event based neutron scattering data: A pathway to study the dynamic structural responses of materials

Time-resolved diffraction has become a vital tool for probing dynamic responses to an applied stimulus. Such experiments traditionally use hardware solutions to histogram measured data into their respective bin. We will show that a major advantage of event-based data acquisition, which time-stamps measured diffraction data with 100 ns accuracy, is much preferred over hardware histogramming of the data by enabling postprocessing for advanced custom binning using a software solution. This approach is made even more powerful by coupling measured diffraction data with metadata about the applied stimuli and material response. In this work, we present a time-filter approach that leverages the power of event-based diffraction collection to reduce stroboscopic data measured over many hours into equally weighted segments that represent subsets of the response to a single cycle of the applied stimulus. We demonstrate this approach by observing ferroelectric/ferroelastic domain wall motion during electric field cycling of BaTiO 3 . Furthermore, the developed approach can readily be expanded to investigate other dynamic phenomena using complex sample environments.

47 OTHER INSTRUMENTATION↗

Persistent Octrees for Parallel Mesh Refinement Through Non-Volatile Byte-Addressable Memory

We report that octree-based mesh adaptation has enabled simulations of complex physical phenomena. Existing meshing algorithms were proposed with the assumption that computer memory is volatile. Consequently, for failure recovery, the in-core algorithms need to save memory states as snapshots with slow file I/Os. The out-of-core algorithms store octants on disks for persistence. However, neither of them was designed to leverage unique characteristics of non-volatile byte-addressable memory (NVBM). In this paper, we propose a novel data structure Distributed Persistent Merged octree (DPM-octree) for both meshing and in-memory storage of persistent octrees using NVBM. It is a multi-version data structure and can recover from failures using its earlier persistent version stored in NVBM. In addition, we design a feature-directed sampling approach to help dynamically transform the DPM-octree layout for reducing NVBM-induced memory write latency. DPM-octree uses parity trees which are created using erasure coding and stored in NVBM to support low-latency in-memory octant recovery after data loss. DPM-octree has been successfully integrated with Gerris software for simulation of fluid dynamics. Finally, our experimental results with real-world scientific workloads show that DPM-octree scales up to 1.1 billion mesh elements with 1000 processors on the Titan supercomputer.

97 MATHEMATICS AND COMPUTING↗

Role of Mg{sup 2+}, Sr{sup 2+}, and F{sup –} ions in octacalcium phosphate crystallization

We have synthesized octacalcium phosphate (OCP) in the presence of inorganic additives (magnesium, strontium, and fluoride ions) and studied the composition, morphology, thermal stability, and dynamic dissolution of the samples thus obtained. It has been shown that, in addition to OCP, magnesium and strontium ions favor the formation of brushite and hydroxyapatite (HA), whereas fluoride ions favor the formation of HA and fluorohydroxyapatite (FHA). We have proposed a process for the preparation of powder materials whose resorption kinetics in corrosive liquid media are corrected by adding dopants capable of activating the dissolution process.

36 MATERIALS SCIENCE↗

Coupled Electronic–Nuclear Magnetostatic Oscillations in Magnetic Materials

Magnetization dynamics in spheroidal ferromagnetic samples is studied theoretically. It is shown that, in a magnetostatic approximation, electronic-nuclear magnetostatic modes with a discrete spectrum of eigen oscillations exist in such samples. The structure and the field dependence of the frequencies of these oscillations are substantially dependent on the parameter of the spheroid shape that is the ratio of it’s axes and also on the external magnetic field value. In a certain region of values of these parameters, the dependence of the eigenfrequencies of the system on them becomes nontrivial. The external field strength and the sample shape determine not only the eigenfrequencies in the system, but also the number of the eigenmodes. In addition, for each of the eigenmodes, there exists a “forbidden” region of magnetic fields and shape parameters in which this mode cannot be observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of γ-Radiation and Addition of Silicone Dioxide on the Properties of Polytetrafluoroethylene

We perform a comparative study of melting of polytetrafluoroethylene (PTFE) in its unmodified form, after irradiation, and after adding plant-derived silicon dioxide as an additive (content, 1%). Calculation of peak shapes on thermal capacity profiles within the theory of smeared phase transitions enable us to identify the structural changes and their features resulting from irradiation of PTFE and incorporation of the additive in it. α-relaxation at the glass transition in our samples is studied by dynamic mechanical analysis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kraken2 Metagenomic Virus Database

The Database: Kraken2 [1] database built from a classification tree containing over 700k metagenomic viruses from JGI IMG/VR [2]. (1) Wood, D. E., Lu, J., and Langmead, B. (2019). Improved metagenomic analysis with Kraken 2. Genome Biol., 20(1), 1–13. doi: 10.1186/s13059-019-1891-0 (2) Paez-Espino D, Chen I-MA, Palaniappan K, Ratner A, Chu K, Szeto E, et al. IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses. Nucleic Acids Res. 2017;45:D457-65. For Paper: Title: A k-mer based approach for virus classification in metatranscriptomic and metagenomic samples identifies viral associations in the Populus phytobiome and autism brains Abstract Background Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an important role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which further limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples. Methods To identify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses. We then integrated the viral classification tree with the NCBI taxonomy for use with ParaKraken, a metagenomic/transcriptomic classifier. Results To illustrate the breadth of our utility for classifying viruses with ParaKraken, we analyzed data from a plant metagenome study identifying the differences between two Populus genotypes in three different compartments and on a human metatranscriptome study identifying the differences between Autism Spectrum Disorder patients and controls in post mortem brain biopsies. In the Populus study, we identified genotype and compartment specific viral signatures, while in the Autism study we identified a significant increased abundance of eight viral sequences in Autism brain biopsies. Conclusion Viruses represent an important aspect of the microbiome. The ability to classify viruses represents the first step in being able to better understand their role in the microbiome. The viral classification method presented here allows for more complete identification of viral sequences for use in identifying associations between viruses and the host and viruses and other microbiome members. Acknowledgements and Funding This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725. This research was also supported by the Plant-Microbe Interfaces Scientific Focus Area in the Genomic Science Program, the Office of Biological and Environmental Research (BER) in the U.S. Department of Energy Office of Science, and by the Department of Energy, Laboratory Directed Research and Development funding (ProjectID 8321), at the Oak Ridge National Laboratory. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the US DOE under contract DE-AC05-00OR22725. This research used resources of the Compute and Data Environment for Science (CADES).

59 BASIC BIOLOGICAL SCIENCES↗

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗

A Convex Data-Driven Approach for Nonlinear Control Synthesis

We consider a class of nonlinear control synthesis problems where the underlying mathematical models are not explicitly known. We propose a data-driven approach to stabilize the systems when only sample trajectories of the dynamics are accessible. Our method is built on the density-function-based stability certificate that is the dual to the Lyapunov function for dynamic systems. Unlike Lyapunov-based methods, density functions lead to a convex formulation for a joint search of the control strategy and the stability certificate. This type of convex problem can be solved efficiently using the machinery of the sum of squares (SOS). For the data-driven part, we exploit the fact that the duality results in the stability theory can be understood through the lens of Perron–Frobenius and Koopman operators. This allows us to use data-driven methods to approximate these operators and combine them with the SOS techniques to establish a convex formulation of control synthesis. The efficacy of the proposed approach is demonstrated through several examples.

97 MATHEMATICS AND COMPUTING↗

Development of Masitinib Derivatives with Enhanced Mpro Ligand Efficiency and Reduced Cytotoxicity

Recently, a high-throughput screen of 1900 clinically used drugs identified masitinib, an orally bioavailable tyrosine kinase inhibitor, as a potential treatment for COVID-19. Masitinib acts as a broad-spectrum inhibitor for human coronaviruses, including SARS-CoV-2 and several of its variants. In this work, we rely on atomistic molecular dynamics simulations with advanced sampling methods to develop a deeper understanding of masitinib’s mechanism of M pro inhibition. To improve the inhibitory efficiency and to increase the ligand selectivity for the viral target, we determined the minimal portion of the molecule (fragment) that is responsible for most of the interactions that arise within the masitinib-Mpro complex. We found that masitinib forms highly stable and specific H-bond interactions with Mpro through its pyridine and aminothiazole rings. Importantly, the interaction with His 163 is a key anchoring point of the inhibitor, and its perturbation leads to ligand unbinding within nanoseconds. Based on these observations, a small library of rationally designed masitinib derivatives (M1–M5) was proposed. Our results show increased inhibitory efficiency and highly reduced cytotoxicity for the M3 and M4 derivatives compared to masitinib.

60 APPLIED LIFE SCIENCES↗

Giant Metrewave Radio Telescope Observations of Head–Tail Radio Galaxies

We present results from a study of seven large known head–tail radio galaxies based on observations using the Giant Metrewave Radio Telescope at 240 and 610 MHz. These observations are used to study the radio morphologies and distribution of the spectral indices across the sources. The overall morphology of the radio tails of these sources is suggestive of random motions of the optical host around the cluster potential. The presence of multiple bends and wiggles in several head–tail sources is possibly due to the precessing radio jets. We find steepening of the spectral index along the radio tails. The prevailing equipartition magnetic field also decreases along the radio tails of these sources. These steepening trends are attributed to the synchrotron aging of plasma toward the ends of the tails. The dynamical ages of these sample sources have been estimated to be ∼10{sup 8} yr, which is a factor of six more than the age estimates from the radiative losses due to synchrotron cooling.

79 ASTRONOMY AND ASTROPHYSICS↗