Search NASA⌕ Search

SEARCH · Search NASA

Results for “force method”

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 343 records · Page 19

PySAGES: flexible, advanced sampling methods accelerated with GPUs

Abstract Molecular simulations are an important tool for research in physics, chemistry, and biology. The capabilities of simulations can be greatly expanded by providing access to advanced sampling methods and techniques that permit calculation of the relevant underlying free energy landscapes. In this sense, software that can be seamlessly adapted to a broad range of complex systems is essential. Building on past efforts to provide open-source community-supported software for advanced sampling, we introduce PySAGES, a Python implementation of the Software Suite for Advanced General Ensemble Simulations (SSAGES) that provides full GPU support for massively parallel applications of enhanced sampling methods such as adaptive biasing forces, harmonic bias, or forward flux sampling in the context of molecular dynamics simulations. By providing an intuitive interface that facilitates the management of a system’s configuration, the inclusion of new collective variables, and the implementation of sophisticated free energy-based sampling methods, the PySAGES library serves as a general platform for the development and implementation of emerging simulation techniques. The capabilities, core features, and computational performance of this tool are demonstrated with clear and concise examples pertaining to different classes of molecular systems. We anticipate that PySAGES will provide the scientific community with a robust and easily accessible platform to accelerate simulations, improve sampling, and enable facile estimation of free energies for a wide range of materials and processes.

Chemistry↗

Efficient molecular dynamics using geodesic integration and solvent–solute splitting

We present an approach to Langevin dynamics in the presence of holonomic constraints based on decomposition of the system into components representing geodesic flow, constrained impulse and constrained diffusion. We show that a particular ordering of the components results in an integrator that is an order of magnitude more accurate for configurational averages than existing alternatives. Moreover, by combining the geodesic integration method with a solvent–solute force splitting, we demonstrate that stepsizes of at least 8 fs can be used for solvated biomolecules with high sampling accuracy and without substantially altering diffusion rates, approximately increasing by a factor of two the efficiency of molecular dynamics sampling for such systems. The methods described in this article are easily implemented using the standard apparatus of modern simulation codes.

97 MATHEMATICS AND COMPUTING↗

Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application

Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Extension of the sasCIF format and its applications for data processing and deposition

Recent advances in small-angle scattering (SAS) experimental facilities and data analysis methods have prompted a dramatic increase in the number of users and of projects conducted, causing an upsurge in the number of objects studied, experimental data available and structural models generated. To organize the data and models and make them accessible to the community, the Task Forces on SAS and hybrid methods for the International Union of Crystallography and the Worldwide Protein Data Bank envisage developing a federated approach to SAS data and model archiving. Within the framework of this approach, the existing databases may exchange information and provide independent but synchronized entries to users. At present, ways of exchanging information between the various SAS databases are not established, leading to possible duplication and incompatibility of entries, and limiting the opportunities for data-driven research for SAS users. In this work, a solution is developed to resolve these issues and provide a universal exchange format for the community, based on the use of the widely adopted crystallographic information framework (CIF). The previous version of the sasCIF format, implemented as an extension of the core CIF dictionary, has been available since 2000 to facilitate SAS data exchange between laboratories. The sasCIF format has now been extended to describe comprehensively the necessary experimental information, results and models, including relevant metadata for SAS data analysis and for deposition into a database. Processing tools for these files ( sasCIFtools ) have been developed, and these are available both as standalone open-source programs and integrated into the SAS Biological Data Bank, allowing the export and import of data entries as sasCIF files. Software modules to save the relevant information directly from beamline data-processing pipelines in sasCIF format are also developed. Lastly, this update of sasCIF and the relevant tools are an important step in the standardization of the way SAS data are presented and exchanged, to make the results easily accessible to users and to promote further the application of SAS in the structural biology community.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Leveraging neural control variates for enhanced precision in lattice field theory

Results obtained with stochastic methods have an inherent uncertainty due to the finite number of samples that can be achieved in practice. In lattice QCD this problem is particularly salient in some observables like, for instance, observables involving one or more baryons and it is the main problem preventing the calculation of nuclear forces from first principles. The method of control variables has been used extensively in statistics and it amounts to computing the expectation value of the difference between the observable of interest and another observable whose average is known to be zero but is correlated with the observable of interest. Recently, control variates methods emerged as a promising solution in the context of lattice field theories. In our current study, instead of relying on an educated guess to determine the control variate, we utilize a neural network to parametrize this function. Using 1 + 1 dimensional scalar field theory as a testbed, we demonstrate that this neural network approach yields substantial improvements. Notably, our findings indicate that the neural network ansatz is particularly effective in the strong coupling regime. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Assessing Gaseous Hydrogen Assisted Fatigue Crack Growth Susceptibility of Pipeline Steel Weld Fusion Zones and Heat Affected Zones

The objective of this work was twofold: (1) measure reliable fatigue crack growth relationships for X65 steel and its girth weld in high-pressure hydrogen gas to enable structural integrity assessments of hydrogen pipelines, and (2) evaluate the hydrogen accelerated fatigue crack growth susceptibility of the weld fusion zone and heat-affected zone relative to the base metal. Fatigue crack growth relationships (da/dN versus ΔK) were measured for girth welded X65 pipeline steel in high pressure hydrogen gas (21 MPa) and in air. Hydrogen assisted fatigue crack growth was observed for the base metal (BM), fusion zone (FZ), and heat-affected zone (HAZ), and was manifested through crack growth rates reaching nearly an order of magnitude acceleration over rates in air. Addditionally, at higher ΔK values, crack growth rates of BM, FZ, and HAZ were coincident; however, at lower ΔK, the fatigue crack growth relationships exhibited some divergence with the fusion zone having the highest crack growth rates. These relative fatigue crack growth rates in the BM, FZ, and HAZ were provisional, however, since both crack closure and residual stress contributed to the crack-tip driving force in specimens extracted from the HAZ. Despite the relatively high applied R-ratio (R = 0.5), crack closure was detected in the heat affected zone tests, in contrast to the absence of crack closure in the base metal tests. Crack closure corrections were performed using the adjusted compliance ratio method and the effect of residual stress on K max was determined by the crack-compliance method. Crack-tip driving forces that account for closure and residual stress effects were quantified as a weighted function of ΔK and K max (i.e., K norm ), and the resulting da/dN versus K norm relationships showed that the HAZ exhibited higher hydrogen accelerated fatigue crack growth rates than the BM at lower K norm values.

36 MATERIALS SCIENCE↗

Facile and quantitative estimation of strain in nanobubbles with arbitrary symmetry in 2D semiconductors verified using hyperspectral nano-optical imaging

When layers of van der Waals materials are deposited via exfoliation or viscoelastic stamping, nanobubbles are sometimes created from aggregated trapped fluids. Though they can be considered a nuisance, nanobubbles have attracted scientific interest in their own right owing to their ability to generate large in-plane strain gradients that lead to rich optoelectronic phenomena, especially in the semiconducting transition metal dichalcogenides. Determination of the strain within the nanobubbles, which is crucial to understanding these effects, can be approximated using elasticity theory. However, the Föppl–von Kármán equations that describe strain in a distorted thin plate are highly nonlinear and often necessitate assuming circular symmetry to achieve an analytical solution. Here, we present an easily implemented numerical method to solve for strain tensors of nanobubbles with arbitrary symmetry in 2D crystals. The method only requires topographic information from atomic force microscopy and the Poisson ratio of the 2D material. We verify that this method reproduces the strain for circularly symmetric nanobubbles that have known analytical solutions. Finally, we use the method to reproduce the Grüneisen parameter of the E' mode for 1L-WS2 nanobubbles on template-stripped Au by comparing the derived strain with measured Raman shifts from tip-enhanced Raman spectroscopy, demonstrating the utility of our method for estimating localized strain in 2D crystals.

2D materials↗

Locating the Source of Forced Oscillations in Transmission Power Grids

A forced oscillation event in power grids refers to a state where malfunctioning or abnormally operating equipment causes persisting periodic disturbances in the system. While power grids are designed to damp most perturbations during standard operations, some of them can excite normal modes of the system and cause significant energy transfers across the system, creating large oscillations thousands of miles away from the source. Localization of the source of such disturbances remains an outstanding challenge due to a limited knowledge of the system parameters outside of the zone of responsibility of system operators. Here, we propose a new method for locating the source of forced oscillations that addresses this challenge by performing a simultaneous dynamic model identification using a principled maximum likelihood approach. We illustrate the validity of the algorithm on a variety of examples where forcing leads to resonance conditions in the system dynamics. Our results establish that an accurate knowledge of system parameters is not required for a successful inference of the source and frequency of a forced oscillation. We anticipate that our method will find a broader application in general dynamical systems that can be well described by their linearized dynamics over short periods of time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Peridynamic Model for Single-Layer Graphene Obtained from Coarse Grained Bond Forces

An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substantial reductions in computational cost compared with molecular dynamics. The calibrated peridynamic model is compared to experimental data on the deflection and perforation of a graphene monolayer by an atomic force microscope probe.

36 MATERIALS SCIENCE↗

Binding partner- and force-promoted changes in αE-catenin conformation probed by native cysteine labeling

Adherens Junctions (AJs) are cell-cell adhesion complexes that sense and propagate mechanical forces by coupling cadherins to the actin cytoskeleton via β-catenin and the F-actin binding protein αE-catenin. When subjected to mechanical force, the cadherin•catenin complex can tightly link to F-actin through αE-catenin, and also recruits the F-actin-binding protein vinculin. In this study, labeling of native cysteines combined with mass spectrometry revealed conformational changes in αE-catenin upon binding to the E-cadherin•β-catenin complex, vinculin and F-actin. A method to apply physiologically meaningful forces in solution revealed force-induced conformational changes in αE-catenin when bound to F-actin. Comparisons of wild-type αE-catenin and a mutant with enhanced vinculin affinity using cysteine labeling and isothermal titration calorimetry provide evidence for allosteric coupling of the N-terminal β-catenin-binding and the middle (M) vinculin-binding domain of αE-catenin. Cysteine labeling also revealed possible crosstalk between the actin-binding domain and the rest of the protein. The data provide insight into how binding partners and mechanical stress can regulate the conformation of full-length αE-catenin, and identify the M domain as a key transmitter of conformational changes.

59 BASIC BIOLOGICAL SCIENCES↗

The Effect of Ion Irradiation on Nanocrystallization and Surface Relief of a Ribbon from Fe{sub 72.5}Cu{sub 1}Nb{sub 2}Mo{sub 1.5}Si{sub 14}B{sub 9} Alloy

Using the methods of X-ray diffraction and atomic force microscopy, the process of crystallization of an amorphous Fe{sub 72.5}Cu{sub 1}Nb{sub 2}Mo{sub 1.5}Si{sub 14}B{sub 9} alloy irradiated with accelerated Ar{sup +} ions is investigated. It is found out that an irradiation by the Ar{sup +} ions with the energy 30 keV at the ion current density 300 μA/cm{sup 2} (fluence 3.75·10{sup 15} cm{sup –2}, irradiation time ~2 s, ion-beam short-duration heating up to 350°С, which is 150°С lower than the thermal crystallization threshold) results in a complete crystallization of this amorphous alloy (throughout the bulk of a 25 μm ribbon) followed by precipitation of solid solution crystals of α-Fe(Si), close in its composition to Fe{sub 80}Si{sub 20}, stable phase of Fe{sub 3}Si, and metastable hexagonal phases. By the methods of atomic force and scanning tunneling microscopy it is shown that nanocrystallization caused by ion irradiation is accompanied by surface relief changes both on the irradiated and unirradiated sides of the Fe{sub 72.5}Cu{sub 1}Nb{sub 2}Mo{sub 1.5}Si{sub 14}B{sub 9} alloy ribbon at the depth exceeding by a factor of 10{sup 3} that of the physical ion penetration for this material. The data obtained, taking into account a significant temperature decrease and multiple acceleration of the crystallization process, serve an evidence of the radiation-dynamic influence of accelerated ions on the metastable amorphous medium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Wavelets as basis functions to represent the coarse-graining potential in multiscale coarse graining approach

In this paper, we apply Multiresolution Analysis (MRA) to develop sparse but accurate representations for the Multiscale Coarse-Graining (MSCG) approximation to the many-body potential of mean force. We rigorously framed the MSCG method into MRA so that all the instruments of this theory become available together with a multitude of new basis functions, namely the wavelets. The coarse-grained (CG) force field is hierarchically decomposed at different resolution levels enabling to choose the most appropriate wavelet family for each physical interaction without requiring an a priori knowledge of the details localization. The representation of the CG potential in this new efficient orthonormal basis leads to a compression of the signal information in few large expansion coefficients. The multiresolution property of the wavelet transform allows to isolate and remove the noise from the CG force-field reconstruction by thresholding the basis function coefficients from each frequency band independently. We discuss the implementation of our wavelet-based MSCG approach and demonstrate its accuracy using two different condensed-phase systems, i.e. liquid water and methanol. Simulations of liquid argon have also been performed using a one-to-one mapping between atomistic and CG sites. The latter model allows to verify the accuracy of the method and to test different choices of wavelet families. Furthermore, the results of the computer simulations show that the efficiency and sparsity of the representation of the CG force field can be traced back to the mathematical properties of the chosen family of wavelets. This result is in agreement with what is known from the theory of multiresolution analysis of signals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GeoVision Analysis: Exploration Task Force Report

The GeoVision study effort included ground-breaking, detailed research on current and future market conditions and geothermal technologies in order to forecast and quantify the electric and non-electric deployment potentials under a range of scenarios, in addition to their impacts on the Nation’s jobs, economy and environment. Coordinated by the U.S. Department of Energy’s (DOE’s) Geothermal Technologies Office (GTO), the GeoVision study development relied on the collection, modeling, and analysis of robust datasets through seven national laboratory partners, which were organized into eight technical Task Force groups. The purpose of this report is to provide a central repository for the research conducted by the Exploration Task Force. The Exploration Task Force consists of four individuals representing three national laboratories: Patrick Dobson (task lead) and Christine Doughty of Lawrence Berkeley National Laboratory, Anna Wall of National Renewable Energy Laboratory, Travis McLing of Idaho National Laboratory, and Chester Weiss of Sandia National Laboratories. As part of the GeoVision analysis, our team conducted extensive scientific and financial analyses on a number of topics related to current and future geothermal exploration methods. The GeoVision Exploration Task Force complements the drilling and resource technology investigations conducted as part of the Reservoir Maintenance and Development Task Force. The Exploration Task Force however has focused primarily on early stage R&D technologies in exploration and confirmation drilling, along with an evaluation of geothermal financing challenges and assumptions, and innovative “blue-sky” technologies. This research was used to develop geothermal resource supply curves (through the use of GETEM) for use in the ReEDS capacity expansion modeling that determines geothermal technology deployment potential. It also catalogues and explores the large array of early-stage R&D technologies with the potential to dramatically reduce exploration and geothermal development costs, forming the basis of the GeoVision Technology Improvement (TI) scenario. These modeling topics are covered in detail in Potential to Penetration task force report. Most of the research contained herein has been published in peer-reviewed papers or conference proceedings and are cited and referenced accordingly. The sections that follow provide a central repository for all of the research findings of the Exploration and Confirmation Task Force. In summary, it provides a comprehensive discussion of Engineered Geothermal Systems (EGS) and associated technology challenges, the risks and costs of conducting geothermal exploration, a review of existing government efforts to date in advancing early-stage R&D in both exploration and EGS technologies, as well as a discussion of promising and innovative technologies and implementation of blue-sky concepts that could significantly reduce costs, lower risks, and shorten the time needed to explore and develop geothermal resources of all types.

15 GEOTHERMAL ENERGY↗

Piezoelectric property of PZT nanofibers characterized by resonant piezo-force microscopy

Nano-piezoelectric materials have drawn tremendous research interest. However, characterization of their piezoelectric properties, especially measuring the piezoelectric strain coefficients, remains a challenge. Normally, researchers use an AFM-based method to directly measure nano-materials’ piezoelectric strain coefficients. But, the extremely small piezoelectric deformation, the influence from the parasitic electrostatic force, and the environmental noise make the measurement results questionable. In this paper, a resonant piezo-force microscopy method was used to accurately measure the piezoelectric deformation from 1D piezoelectric nanofibers. During the experiment, the AFM tip was brought into contact with the piezoelectric sample and set to work at close to its first resonant frequency. A lock-in amplifier was used to pick up the sample’s deformation signal at the testing frequency. By using this technique, the piezoelectric strain constant d 33 of the Lead Zirconate Titanate (PZT) nanofiber with a diameter of 76 nm was measured. The result showed that d 33 of this PZT nanofiber was around 387 pm/V. Meanwhile, by tracking the piezoelectric deformation phase image, domain structures inside PZT nanofibers were identified.

36 MATERIALS SCIENCE↗

Upconverting Nanoparticles as Optical Sensors of Nano- to Micro-Newton Forces

Mechanical forces affect a myriad of processes, from bone growth to material fracture to touch-responsive robotics. While nano- to micro-Newton forces are prevalent at the microscopic scale, few methods have the nanoscopic size and signal stability to measure them in vivo or in situ. Here, we develop an optical force-sensing platform based on sub-25 nm NaYF 4 nanoparticles (NPs) doped with Yb 3+ , Er 3+ , and Mn 2+ . The lanthanides Yb 3+ and Er 3+ enable both photoluminescence and upconversion, while the energetically coupled d-metal Mn 2+ adds force tunability through its crystal field sensitivity. IN using a diamond anvil cell to exert up to 3.5 GPa pressure or ~10 μN force per particle, we track stress-induced spectral responses. The red (660 nm) to green (520, 540 nm) emission ratio varies linearly with pressure, yielding an observed color change from orange to red for α-NaYF 4 and from yellow–green to green for d-metal optimized β-NaYF 4 when illuminated in the near infrared. We record consistent readouts over multiple pressure cycles and hours of illumination. With the nanoscopic size, a dynamic range of 100 nN to 10 μN, and photostability, these nanoparticles lay the foundation for visualizing dynamic mechanical processes, such as stress propagation in materials and force signaling in organisms.

36 MATERIALS SCIENCE↗

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of SnO2 Processing on Ionic Distribution in Double-Cation-Double Halide Perovskites

Moving toward a future of efficient, accessible, and less carbon-reliant energy devices has been at the forefront of energy research innovations for the past 30 years. Metal-halide perovskite (MHP) thin films have gained significant attention due to their flexibility of device applications and tunable capabilities for improving power conversion efficiency. Serving as a gateway to optimize device performance, consideration must be given to chemical synthesis processing techniques. Therefore, how does common substrate processing techniques influence the behavior of MHP phenomena such as ion migration and strain? Here, we demonstrate how a hybrid approach of chemical bath deposition (CBD) and nanoparticle SnO2 substrate processing significantly improves the performance of (FAPbI3)0.97(MAPbBr3)0.03 by reducing micro-strain in the SnO2 lattice, allowing distribution of K+ from K-Cl treatment of substrates to passivate defects formed at the interface and produce higher current in light and dark environments. X-ray diffraction reveals differences in lattice strain behavior with respect to SnO2 substrate processing methods. Through use of conductive atomic force microscopy (c-AFM), conductivity is measured spatially with MHP morphology, showing higher generation of current in both light and dark conditions for films with hybrid processing. Additionally, time-of-flight secondary ionization mass spectrometry (ToF-SIMS) observed the distribution of K+ at the perovskite/SnO2 interface, indicating K+ passivation of defects to improve the power conversion efficiency (PCE) and device stability. We show how understanding the role of ion distribution at the SnO2 and perovskite interface can help reduce the creating of defects and promote a more efficient MHP device.

conductive atomic force microscopy↗

Crack propagation in silica from reactive classical molecular dynamics simulations

Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far–field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J –integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (K IC ), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up–scaled atomic–level energies and stresses through the J –integral. In conclusion, the application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.

36 MATERIALS SCIENCE↗