Search NASA⌕ Search

SEARCH · Search NASA

Results for “Random Interactions”

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 451 records · Page 25

Development of a coarse-grained molecular dynamics model for poly(dimethyl- co -diphenyl)siloxane

Polydimethylsiloxane is an important polymeric material with a wide range of applications. However, environmental effects like low temperature can induce crystallization in this material with resulting changes in its structural and dynamic properties. The incorporation of phenyl-siloxane components, e.g., as in a poly(dimethyl-co-diphenyl)siloxane random copolymer, is known to suppress such crystallization. Molecular dynamics (MD) simulations can be a powerful tool to understand such effects in atomistic detail. Unfortunately, all-atomistic molecular dynamics (AAMD) is limited in both spatial dimensions and simulation times it can probe. Here, to overcome such constraints and to extend to more useful length- and time-scales, we systematically develop a coarse-grained molecular dynamics (CGMD) model for the poly(dimethyl-co-diphenyl)siloxane system with bonded and non-bonded interactions determined from all-atomistic simulations by the iterative Boltzmann inversion (IBI) method. Additionally, we propose a lever rule that can be useful to generate non-bonded potentials for such systems without reference to the all-atomistic ground truth. Our model captures the structural and dynamic properties of the copolymer material with quantitative accuracy and is useful to study long-time dynamics of highly-entangled systems, sequence-dependent properties, phase behaviour, etc.

36 MATERIALS SCIENCE↗

Block copolymer self-assembly to pattern gold nanodots for site-specific placement of DNA origami and attachment of nanomaterials

Directed placement of DNA origami could play a key role in future integrated nanoelectronic devices. Here we demonstrated the site-selective attachment of DNA origami on gold dots formed using a pattern transfer method through block copolymer self-assembly. First, a random copolymer brush layer is grafted on the Si surface and then poly (styrene-b-methylmethacrylate) block copolymer is spin-coated to give a hexagonal nanoarray after annealing. UV irradiation followed by acetic acid etching is used to remove the PMMA, creating cylindrical holes and then oxygen plasma etching removes the random copolymer layer inside those holes. Next, metal evaporation, followed by lift-off creates a gold dot array. We evaluated different ligand functionalization of Au dots, as well as DNA hybridization to attach DNA origami to the nanodots. DNA-coated Au nanorods are assembled on the DNA origami as a step towards creating nanowires and to facilitate electron microscopy characterization of the attachment of DNA origami on these Au nanodots. The DNA hybridization approach showed better DNA attachment to Au nanodots than localization by electrostatic interaction. Here, this work contributes to the understanding of DNA-templated assembly, nanomaterials, and block copolymer nanolithography. Furthermore, the work shows potential for creating DNA-templated nanodevices and their placement in ordered arrays in future nanoelectronics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Machine learning using host/guest energy histograms to predict adsorption in metal–organic frameworks: Application to short alkanes and Xe/Kr mixtures

A machine learning (ML) methodology that uses a histogram of interaction energies has been applied to predict gas adsorption in metal–organic frameworks (MOFs) using results from atomistic grand canonical Monte Carlo (GCMC) simulations as training and test data. In this work, the method is first extended to binary mixtures of spherical species, in particular, Xe and Kr. In addition, it is shown that single-component adsorption of ethane and propane can be predicted in good agreement with GCMC simulation using a histogram of the adsorption energies felt by a methyl probe in conjunction with the random forest ML method. Here, the results for propane can be improved by including a small number of MOF textural properties as descriptors. We also discuss the most significant features, which provides physical insight into the most beneficial adsorption energy sites for a given application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modal dynamics of high-frequency transverse combustion instabilities

This paper investigates the modal dynamics of self-excited high frequency transverse instabilities in a multi-nozzle, swirl-stabilized can combustor. It utilizes results from multiple pressure sensors to decompose the disturbance into clockwise (CW) and counter-clockwise (CCW) waves, enabling a reconstruction of the pressure amplitude, phase, anti-nodal line location, and spin ratio. While these types of decompositions have been demonstrated on small scale combustors, particularly annular combustors, this is the first study to analyze these results for a high power, multi-nozzle can combustor. Its particular focus is in characterizing the system evolution in phase space and the nature of the system attractors under linearly unstable conditions; results are presented in spin ratio-phase difference space, so that limit cycle oscillations appear as fixed points. Results are shown for conditions where three different attractors dominate (standing, CW, and CCW) but are taken for very long time periods (~50,000 cycles) so that the system's interaction with these attractors can be observed. For the standing mode-dominated case, the phase portrait shows spiral trajectories converging to a single fixed point. For the spinning mode dominated cases, however, multiple fixed points and saddle points were observed, and the trajectories showed a more complicated structure. It was also shown that even though the CW mode is dominant, this mode intermittently switches to standing or CCW modes for short period. These results suggest that the strength of each attractor depends on test conditions and that random noise stochastically perturbs the system so that large regions of the phase space are explored. Lastly, one interesting observation is that the angular velocity of the anti-nodal line linearly varies with the spin ratio under standing-mode dominant conditions.

42 ENGINEERING↗

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

Sequence Programmable Order–Disorder Transitions in Supramolecular Assembly of Peptide Nanofibers

Protein–protein interactions determine the assembly of complexes that are responsible for numerous key biological processes. The assembly of many natural protein complexes is mediated by post-translational structural changes and environmental stimuli. In this study, we show that incorporation of adjacent lysine residues results in the pH-tunable stability of peptide secondary structure and assembly, allowing for the incorporation of complementary order-inducing motifs. The strategic placement of cysteine pairs in the same peptide sequence results in redox-dependent disulfide staple formation, inducing a transition from random coil to β-sheet conformation and subsequent supramolecular nanofiber assembly from otherwise disordered peptide monomers. Spectroscopic, imaging, molecular dynamics, and kinetic studies highlight the critical role of sequence motif location, oligomerization, and the competitive interplay between intra- and interpeptide disulfide bonding in determining assembly outcomes. We extend this approach to demonstrate phosphorylation-dependent assembly from the design of the same parent peptide sequence, suggesting a general approach to the design of diverse stimulus-responsive peptide sequences for supramolecular assembly. Furthermore, these findings also provide a framework for investigating sequence-dependent pathways in amyloid fiber formation with potential implications for neurodegenerative disease research.

Disulfides↗

Electron capture of superheavy nuclei with realistic lepton wave functions

The superheavy nuclei push the periodic table of the elements and the chart of the nuclides to their limits, providing a unique laboratory for studies of the electron-nucleus interactions. The most important weak decay mode in known superheavy nuclei is electron capture (EC). In the standard calculations of EC, the lepton wave functions are usually considered in the lowest-order approximation. In this work, we investigate the sensitivity of EC rates on the choice of the electron wave functions by (i) assuming the single-particle approximation for the electron wave functions, and (ii) carrying out Dirac-Hartree-Fock (DHF) calculations. The nuclear response is generated based on the state-of-the-art quasiparticle random phase approximation employing relativistic nuclear energy density functional theory. Here, we show that using the improved lepton wave functions reduces the EC rates up to 40% in the superheavy nucleus oganesson (𝑍=118). Interestingly, because of screening effects, the difference between the EC rates obtained with the DHF and single-particle calculations is fairly small.

Atomic orbital↗

Machine Learning for Neutron Resonance Evaluations [Slides]

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit a resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this presentation, we describe the application of machine learning to automate the quantum number assignments. Scikit-learn classifiers were trained on simulated resonance data whose statistical properties were chosen to mimic real data. We explored the use of several physics (and random matrix theory)-motivated features for training the classifiers, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results demonstrated that we can determine resonance spin groups somewhat reliably. We are now investigating the application of our approach to 52 Cr resonance data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Securing Grid-interactive Efficient Buildings (GEB) through Cyber Defense and Resilient System (CYDRES)

The DOE CYDRES project is driven by the urgent need to address critical research gaps in the domain of cyber-physical security of smart buildings, including Grid-interactive Efficient Buildings (GEBs). CYDRES, a real-time advanced building resilient platform, aims to enhance the cyber-attack-immune capabilities of buildings through multi-layered prevention, detection, and adaptation mechanisms. CYDRES consists of five key modules: a multi-layer network analyzer, an Automatic Fault Detection, Diagnosis, and Prognosis (AFDDP) framework, an intelligent mode selector, a cyber-resilient control framework, and a situation awareness platform. The Network Analyzer employs a data-driven framework that includes a protocol state learning tool and a CRF (Conditional Random Field) command validator. In Hardware-In-the-Loop (HIL) testbeds, it achieved 100% detection accuracy with a false alarm rate of 3%, validating its efficacy in identifying selected cyber-attacks. The AFDDP framework leverages pattern matching, PCA (Principal Component Analysis)-based strategies, and a DBN (Dynamic Bayesian Network)-based fault diagnosis approach to pinpoint the causes of physical system abnormalities using Building Automation System (BAS) data. In HIL experiments, the AFDDP module attained a detection accuracy of over 95% with a false alarm rate below 7%. Additionally, the fault detector utilized machine learning (Random Forest) and deep learning (Multi-Layer Perceptron) methods with acoustic sensor data to achieve a 100% fault detection accuracy in Heating, Ventilation, and Air-Conditioning (HVAC) equipment. The Mode Selector offered real-time impact analysis, allowing immediate actions to protect BASs in the face of emerging threats. The cyber-resilient control framework included an adaptive Model Predictive Control (MPC) and a measurement compensator, reducing temperature violations by up to 94% and improving the total demand flexibility by up to 70% in HIL experiments. Such HIL experiments covered a cyber-attack case and a physical fault case, showcasing CYDRES’ efficiency in maintaining operational continuity during threats. The situation awareness platform in Grafana enhanced real-time threat detection and response visualization, augmenting the operational awareness for building operators. CYDRES demonstrated high technical effectiveness in various test scenarios, particularly in HIL environments. The project's phased development approach ensured efficient use of resources, highlighting its practical feasibility and readiness for commercialization. By enhancing the security and resilience of building operations, CYDRES represents a significant advance in mitigating risks associated with cyber-physical systems, thereby enhancing public confidence in the safety of modern building infrastructure. Future directions for the project include expanding testing protocols, refining AFDDP methodologies, exploring more comprehensive resilient control strategies, and testing in real commercial buildings.

42 ENGINEERING↗

Polymer Conformations and Diffusion Through a Monolayer of Confining Nanoparticles

We present coarse-grained molecular dynamics simulations to probe chain-scale polymer conformations and diffusion between confining nanoparticles (NPs). By constructing a monolayer of hexagonally-packed NPs in a polymer melt with athermal interactions, we observe the magnitude and length-scale over which homogeneously confining NPs impact the polymer behavior, which provides fundamental insights into more complex polymer nanocomposites. We show that polymer conformations are more perturbed under strong confinement (i.e. when the interparticle distance, ID, is less than twice the polymer radius of gyration, 2R g ) as compared to around an isolated NP, and the effect depends on the ratio of R NP /R g rather than either independently. In fact, these conformations can be quantitatively replicated by executing a simple random walk in a similarly confining environment. We also show that polymer diffusion is slowed by the presence of NPs and the slowing persists far beyond the length-scale over which polymer conformations are perturbed. Although the slowing is strongest ~R g from the NPs, the diffusion coefficient is slower even beyond ~5R g from the NPs. Furthermore, by analyzing the directional van Hove distributions, we show polymer diffusion away from the NP monolayer is bulk-like while diffusion through the monolayer is slower with increasing NP confinement. Furthermore, these molecular dynamics simulations provide fundamental insights into the temporal and spatial effect of confining, athermal NPs on chain-scale polymer conformations and diffusion.

36 MATERIALS SCIENCE↗

Interfacial triferroicity in monolayer chromium dihalide

Couplings between different ferroics in two-dimensional (2D) materials, with atomic thickness and van der Waals surface, have been long sought but not yet realized. We show the first-principles evidence of unique triferroic couplings in recently synthesized CrI 2 monolayer. Its ferroelasticity with a low switching barrier stems from the Jahn-Teller effect. When coupled with the interaction of a substrate, a type of interfacial ferroelectricity with direction determined by the ferroelastic states can emerge which can be switched following the same low-barrier pathway of ferroelastic switching. The direction of its striped antiferromagnetism is also governed by its ferroelastic states, leading to a coupling of ferroelasticity, ferroelectricity, and magnetism. Such mechanism of multiferroic couplings, denoted as interfacial Jahn-Teller triferroicity, can be applied to other 2D materials such as CuC l2 and their Janus monolayers, enabling electrical manipulation of spintronics for efficient nonvolatile random-access memories.

36 MATERIALS SCIENCE↗

Origins of magnetic memory and strong exchange bias bordering magnetic compensation in mixed-lanthanide systems

Here, the unexpected physical phenomena resulting from the seemingly inconsequential substitutions of chemically similar lanthanide elements in the Pr 1-x Gd x ScGe system are exploited to further the understanding of rare-earth magnetism and inform materials design. By directly probing magnetic moments of crystallographically indistinguishable Pr and Gd we solve the puzzles of how an unusual magnetic memory and strong exchange bias emerge at specific, easily predictable chemistries. Both effects are rooted in a robust antiparallel arrangement of large 4f magnetic moments of light and heavy lanthanides. This enables precise control of nearly zero net magnetization either opposed to, or aligned with, the external magnetic field that persists over a wide range of temperatures and fields. Further, spontaneous perturbations in the random distribution of lanthanide ions makes strong exchange bias possible in bulk single-phase compounds bordering magnetic compensation, consequently expanding the materials base beyond artificial magnetic multilayers and broadening the range of potential applications of the phenomenon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transcripts and genomic intervals associated with variation in metabolite abundance in maize leaves under field conditions

Abstract Plants exhibit extensive environment-dependent intraspecific metabolic variation, which likely plays a role in determining variation in whole plant phenotypes. However, much of the work seeking to use natural variation to link genes and transcript’s impacts on plant metabolism has employed data from controlled environments. Here, we generated and analyzed data on the variation in the abundance of 26 metabolites across 660 maize inbred lines under field conditions. We employ these data and previously published transcript and whole plant phenotype data reported for the same field experiment to identify both genomic intervals (through genome-wide association studies (GWAS)) and transcripts (using both transcriptome-wide association studies (TWAS) and an explainable artificial intelligence (AI) approach based on random forest (RF)) associated with variation in metabolite abundance. Both genome-wide association and random forest-based methods identified substantial numbers of significant associations including genes with plausible links to the metabolites they are associated with. In contrast, the transcriptome-wide association identified only six significant associations. In three cases, genetic markers associated with metabolic variation in our study colocalized with markers linked to variation in non-metabolic traits scored in the same experiment. We speculate that the poor performance of transcriptome-wide association studies in identifying transcript-metabolite associations may reflect a high prevalence of non-linear interactions between transcripts and metabolites and/or a bias towards rare transcripts playing a large role in determining intraspecific metabolic variation.

Mathivanan, Ramesh Kanna↗

Unraveling the dislocation–precipitate interactions in high-entropy alloys

The precipitates play a significant role in not only enhancing the strength, but also maintaining the high toughness in alloys. However, the interactions of the nanoscale precipitates with dislocations in the high entropy alloys (HEAs) are difficult to observe directly by in-situ TEM experiments due to the limits of the resolution and time. Here, using atomic simulations we report the synergistic strengthening of the coherent precipitate and atomic-scale lattice distortion in the HEAs at cryogenic/elevated temperatures. The effects of temperature, chemical disorder, precipitate spacing, precipitate size, elemental segregation, and dislocation-cutting number on the critical stress for the dislocation to overcome a row of precipitates are studied. A random stacking fault energy landscape along the slip plane, the lattice distortion at different temperatures, and the interface/surface energy at various crystallographic orientations are obtained. Compared with the traditional metals and alloys, HEAs have the severe atomic-scale lattice distortions to generate the local high tensile/compressive stress fields. This complex stress causes the dislocation line to bend, and thus improves the dislocation slip resistance, resulting in the strong solid-solution strengthening. The stacking fault strengthening induced by the obvious difference of the stacking fault energies between the HEA matrix and precipitate (within the inner of the HEA matrix), and the formation of the antiphase domain boundary contribute to the high strength. The precipitate embedded by the solute atoms produces the strong lattice distortion to enhance the dislocation slip resistance at high temperatures. Hence, the current results provide the mechanistic insight into the phenomenon that the coherent precipitate combined with the severe atomic-scale lattice distortion can enhance the strength at cryogenic/elevated temperatures to further broaden the scope of applications of advanced HEAs.

36 MATERIALS SCIENCE↗

Storylines for the 1997 New Year’s Flood: The role of watershed antecedent conditions and future warming in shaping discharge in the Truckee River watershed

The 1997 New Year’s flood was among the most devastating floods in the Truckee River watershed located in western Nevada. This event resulted from complex interactions of flood drivers, such as extreme precipitation, wet antecedent watershed conditions, warm temperatures and rapid snowmelt. We leveraged simulated forcings from the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) and a process-based hydrological model to recreate the 1997 New Year’s flood for the Truckee River watershed across four climate warming levels ranging from the current temperatures to + 4° C. For each scenario, we conducted ensemble simulations with the same forcing but with 100 different seasonal watershed antecedent conditions, which were randomly sampled from long-term hydrological simulations. The results show that the 1997 New Year’s flood can be reproduced or exceeded consistently only when the antecedent watershed conditions are wet, specifically when streamflows are above the 75th percentile of the climatological value. There is negligible change in ensemble mean peakflows for Truckee River near Reno; however, there are increases of 18% and 14% under the warming levels of + 3° C and + 4° C, respectively. The increases in peakflows under future climate warming are attributed to wetter antecedent watershed conditions and enhanced snowmelt. Furthermore, the largest increases in peakflows occur at small, high-elevation headwater basins along the Sierra Nevada crest. This study highlights that changes in extreme flood events will result from the complex interplay of multiple flood drivers. It also demonstrates the potential of storyline approaches to analyze future realizations of these extreme events under different climate scenarios.

Climate change impact study↗

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗

Heisenberg-limited Hamiltonian learning for interacting bosons

We develop a protocol for learning a class of interacting bosonic Hamiltonians from dynamics with Heisenberg-limited scaling. For Hamiltonians with an underlying bounded-degree graph structure, we can learn all parameters with root mean square error ϵ using ${\mathcal{O}}(1/\epsilon )$ total evolution time, which is independent of the system size, in a way that is robust against state-preparation and measurement error. In the protocol, we only use bosonic coherent states, beam splitters, phase shifters, and homodyne measurements, which are easy to implement on many experimental platforms. A key technique we develop is to apply random unitaries to enforce symmetry in the effective Hamiltonian, which may be of independent interest.

computer science↗

Nonlocal correlations in iron pnictides and chalcogenides

Deviations of low-energy electronic structurse of iron-based superconductors from density-functional-theory predictions have been parametrized in terms of band- and orbital-dependent mass renormalizations and energy shifts. The former have typically been described in terms of a local self-energy within the framework of dynamical mean field theory, while the latter appears to require nonlocal effects due to interband scattering. By calculating the renormalized band structure in both random phase approximation (RPA) and the two-particle self-consistent approximation (TPSC), we show that correlations in pnictide systems like LaFeAsO and LiFeAs can be described rather well by a nonlocal self-energy. In particular, Fermi pocket shrinkage as seen in experiments occurs due to repulsive interband finite-energy scattering. For the canonical iron chalcogenide system FeSe in its bulk tetragonal phase, the situation is, however, more complex since even including momentum-dependent band renormalizations cannot explain experimental findings. We propose that the nearest-neighbor Coulomb interaction may play an important role in band-structure renormalization in FeSe. Finally, we further compare our evaluations of nonlocal quasiparticle scattering lifetime within RPA and TPSC with experimental data for LiFeAs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗