Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

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 667 records · Page 37

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗

Polyamorphism in a solute-lean Al–Ce metallic glass

Polyamorphism discovered in lanthanide-rich metallic glasses (MGs) has been attributed to the electronic transition of the lanthanide element as a solvent element. In this work, we report that pressure-induced polyamorphism still exists in a Ce-poor Al 93 Ce 7 binary MG where the 4f electron element serves as a solute and solute–solute avoidance is expected. The polyamorphic transition, observed by in situ high-pressure synchrotron x-ray diffraction, is accompanied by a volume collapse of ~0.78% and occurs over a narrow pressure range from ~0.8 to ~1.8 GPa. Further synchrotron Ce L 3 -edge x-ray absorption spectroscopy measurements reveal that pressure-induced 4f electron delocalization underlies the polyamorphic transition. Molecular dynamics simulations confirm that the Ce atoms in the MG are completely isolated by the solvent Al atoms. This result demonstrates that 4f element-bearing alloys with extremely dilute concentrations can also exhibit polyamorphic states originating from electronic transitions, extending the compositional space of polyamorphism of MGs into very dilute regions. Our work suggests that tunable properties under compressive stress could be achieved in MGs by even minor doping of elements prone to electronic transitions.

36 MATERIALS SCIENCE↗

Computer simulation of carbonization and graphitization of coal

Abstract This study describes computer simulations of carbonization and graphite formation, including the effects of hydrogen, nitrogen, oxygen, and sulfur. We introduce a novel technique to simulate carbonization, ‘Simulation of Thermal Emission of Atoms and Molecules (STEAM)’, designed to elucidate volatile outgassing and density variations in the intermediate material during carbonization. The investigation analyzes the functional groups that endure through high-temperature carbonization and examines the graphitization processes in carbon-rich materials containing non-carbon impurity elements. The physical, vibrational, and electronic attributes of impure amorphous graphite are analyzed, and the impact of nitrogen on electronic conduction is investigated.

Ugwumadu, C. (ORCID:0000000199207594)↗

Exploring dielectric properties in atomistic models of amorphous boron nitride

Abstract We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about 100 atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electronic model, we analyse the dielectric functions for complex three dimensional models (containing about 10.000 atoms) embedding varying concentrations of sp 1 , sp 2 and sp 3 bonds between B and N atoms. Within the limits of these methodologies, the resulting value of the zero-frequency dielectric constant is shown to be influenced by the population density of such mid-gap states and their localization characteristics. We observe nontrivial correlations between the structure-induced electronic fluctuations and the resulting dielectric constant values. Our findings are however just a first step in the quest of accessing fully accurate dielectric properties of as-grown amorphous BN of relevance for interconnect technologies and beyond.

Materials Science↗

Mechanistic studies of small molecule ligands selective to RNA single G bulges

Abstract Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated molecular dynamics simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity relationship studies and transverse relaxation measurements by nuclear magnetic resonance spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Our work deepened the understanding of RNA-small molecule interactions and integrated a new framework for the rational design of selective small-molecule RNA binders.

Biochemistry & Molecular Biology↗

Direct observation of anisotropic growth of water films on minerals driven by defects and surface tension

Knowledge of the occurrences of water films on minerals is critical for global biogeochemical and atmospheric processes, including element cycling and ice nucleation. The underlying mechanisms controlling water film growth are, however, misunderstood. Using infrared nanospectroscopy, amplitude-modulated atomic force microscopy, and molecular simulations, we show how water films grow from water vapor on hydrophilic mineral nanoparticles. We imaged films with up to four water layers that grow anisotropically over a single face. Growth usually begins from the near edges of a face where defects preferentially capture water vapor. Thicker films produced by condensation cooling completely engulf nanoparticles and form thicker menisci over defects. The high surface tension of water smooths film surfaces and produces films of inhomogeneous thickness. Nanoscale topography and film surface energy thereby control anisotropic distributions and thicknesses of water films on hydrophilic mineral nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vibrational control of the reaction pathway in the H + CHD 3 → H 2 + CD 3 reaction

An accurate full-dimensional quantum state-to-state simulation of the six-atom title reaction based on first-principles theory is reported. Counterintuitive effects are found: Increasing the energy in the reactant’s CD 3 umbrella vibration reduces the energy in the corresponding product vibration. An in-depth analysis reveals the crucial role of the effective dynamical transition state: Its geometry is controlled by the vibrational states of the reactants and subsequently controls the quantum state distribution of the products. This finding enables generalizing the concept of transition state control of chemical reactions to the quantum state–specific level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport Phenomena in Low Temperature Lithium-Ion Battery Electrolytes

Lithium-ion batteries face low temperature performance issues, limiting the adoption of technologies ranging from electric vehicles to stationary grid storage. This problem is thought to be exacerbated by slow transport within the electrolyte, which in turn may be influenced by ion association, solvent viscosity, and cation transference number. How these factors collectively impact low temperature transport phenomena, however, remains poorly understood. Here we show using all-atom classical molecular dynamics (MD) simulations that the dominant factor influencing low temperature transport in LP57 (1 M LiPF 6 in 3:7 ethylene carbonate (EC)/ethyl methyl carbonate (EMC)) is solvent viscosity, rather than ion aggregation or cation transference number. We find that ion association decreases with decreasing temperature, while the cation transference number is positive and roughly independent of temperature. In an effort to improve low temperature performance, we introduce γ -butyrolactone (GBL) as a low viscosity co-solvent to explore two alternative formulations: 1 M LiPF 6 in 15:15:70 EC/GBL/EMC and 3:7 GBL/EMC. While GBL reduces solution viscosity, its low dielectric constant results in increased ion pairing, yielding neither improved bulk ionic conductivity nor appreciably altered ion transport mechanisms. We expect that these results will enhance understanding of low temperature transport and inform the development of superior electrolytes.

25 ENERGY STORAGE↗

mpfit: a robust method for fitting atomic resolution images with multiple Gaussian peaks

The standard technique for sub-pixel estimation of atom positions from atomic resolution scanning transmission electron microscopy images relies on fitting intensity maxima or minima with a two-dimensional Gaussian function. While this is a widespread method of measurement, it can be error prone in images with non-zero aberrations, strong intensity differences between adjacent atoms or in situations where the neighboring atom positions approach the resolution limit of the microscope. Here we demonstrate mpfit, an atom finding algorithm that iteratively calculates a series of overlapping two-dimensional Gaussian functions to fit the experimental dataset and then subsequently uses a subset of the calculated Gaussian functions to perform sub-pixel refinement of atom positions. Based on both simulated and experimental datasets presented in this work, this approach gives lower errors when compared to the commonly used single Gaussian peak fitting approach and demonstrates increased robustness over a wider range of experimental conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leiomodin creates a leaky cap at the pointed end of actin-thin filaments

Improper lengths of actin-thin filaments are associated with altered contractile activity and lethal myopathies. Leiomodin, a member of the tropomodulin family of proteins, is critical in thin filament assembly and maintenance; however, its role is under dispute. Using nuclear magnetic resonance data and molecular dynamics simulations, we generated the first atomic structural model of the binding interface between the tropomyosin-binding site of cardiac leiomodin and the N-terminus of striated muscle tropomyosin. Our structural data indicate that the leiomodin/tropomyosin complex only forms at the pointed end of thin filaments, where the tropomyosin N-terminus is not blocked by an adjacent tropomyosin protomer. This discovery provides evidence supporting the debated mechanism where leiomodin and tropomodulin regulate thin filament lengths by competing for thin filament binding. Data from experiments performed in cardiomyocytes provide additional support for the competition model; specifically, expression of a leiomodin mutant that is unable to interact with tropomyosin fails to displace tropomodulin at thin filament pointed ends and fails to elongate thin filaments. Together with previous structural and biochemical data, we now propose a molecular mechanism of actin polymerization at the pointed end in the presence of bound leiomodin. In the proposed model, the N-terminal actin-binding site of leiomodin can act as a “swinging gate” allowing limited actin polymerization, thus making leiomodin a leaky pointed-end cap. Results presented in this work answer long-standing questions about the role of leiomodin in thin filament length regulation and maintenance.

2D-NMR↗

Atomization of liquid sheets in high pressure airflow

An investigation of liquid sheet atomization is made with combustor simulated inlet air pressures varied from 0.10 to 2.1 MPa. Mean drop diameters are measured with an improved scanning radiometer and correlated with the liquid and air stream Reynolds numbers, RE(1) and RE(A) and the airstream pressure sensitive group GC(2). These data are used in the modeling of the combustion process.

Ingebo, R. D.↗

Electron stimulated desorption of atomic oxygen from silver

The electron stimulated desorption (ESD) of neutral oxygen atoms from polycrystalline silver and of oxygen ions from Ag(110) has been studied. Polycrystalline Ag charged with (16)O2 and (18)O2 and bombarded by low-energy electrons (approx 100 eV) under ultrahigh vacuum (UHV) conditions emitted O atom flux levels of 1 x 10 to the 12th power/sq cm/s at a Ag temperature of 300 C. The flux was detected with a quadrupole mass spectrometer operating in the appearance potential mode. The neutral cross section at about 100 C was determined to be 7 x 10 to the -19 sq cm. Ancillary experiments conducted in a UHV chamber equipped with a cylindrical mirror analyzer and rigged for ion energy distribution and ion angular distribution were used to study O ions desorbed from Ag(110). Two primary O(+) energies of 2.4 and 5.4 eV were detected from the Ag(110) after having been dosed with 2500 L of (16)O2. It also appears that in both experiments there was strong evidence for directionality of the emitted flux. The results of this study serve as a proof of concept for the development of a laboratory atomic oxygen beam generator that simulates the gas flux environment experienced by orbiting vehicles.

Outlaw, R. A.↗

Molecular arrangements and conformations of liquid n-tridecane chains confined between two hard walls

Monte Carlo simulations based on united CH2 atoms linked with fixed bond lengths and angles, and continuously varying torsional angles subject to appropriate interaction energies, were used to study liquid n-tridecane chains confined between two parallel hard walls with 5 mm separation. The molecules located in the central region are noted to be in very close agreement with ideal unperturbed chains. Nearly all the chain units in the first segmental layer adjacent to the walls belong to two-dimensional chains, and exhibit considerable orientational correlations between neighboring segments of different chains.

Vacatello, Michele↗

Thermodynamic Behavior of Nano-sized Gold Clusters on the (001) Surface

We have studied thermal expansion of the surface layers of the hexagonally reconstructed Au (001) surface using a classical Molecular Dynamics (MD) simulation technique with an Embedded Atomic Method (EAM) type many-body potential. We find that the top-most hexagonal layer contracts as temperature increases, whereas the second layer expands or contracts depending on the system size. The magnitude of expansion coefficient of the top layer is much larger than that of the other layers. The calculated thermal expansion coefficients of the top-most layer are about -4.93 x 10(exp -5)Angstroms/Kelvin for the (262 x 227)Angstrom cluster and -3.05 x 10(exp -5)Angstroms/Kelvin for (101 x 87)Angstrom cluster. The Fast Fourier Transform (FFT) image of the atomic density shows that there exists a rotated domain of the top-most hexagonal cluster with rotation angle close to 1 degree at temperature T less than 1000Kelvin. As the temperature increases this domain undergoes a surface orientational phase transition. These predictions are in good agreement with previous phenomenological theories and experimental studies.

Paik, Sun M.↗

Enabling Computational Nanotechnology through JavaGenes in a Cycle Scavenging Environment

A genetic algorithm procedure is developed and implemented for fitting parameters for many-body inter-atomic force field functions for simulating nanotechnology atomistic applications using portable Java on cycle-scavenged heterogeneous workstations. Given a physics based analytic functional form for the force field, correlated parameters in a multi-dimensional environment are typically chosen to fit properties given either by experiments and/or by higher accuracy quantum mechanical simulations. The implementation automates this tedious procedure using an evolutionary computing algorithm operating on hundreds of cycle-scavenged computers. As a proof of concept, we demonstrate the procedure for evaluating the Stillinger-Weber (S-W) potential by (a) reproducing the published parameters for Si using S-W energies in the fitness function, and (b) evolving a "new" set of parameters using semi-empirical tightbinding energies in the fitness function. The "new" parameters are significantly better suited for Si cluster energies and forces as compared to even the published S-W potential.

Globus, Al↗

NASA Tech Briefs, May 2008

Topics covered inclde: Deployable Wireless Camera Penetrators; Hand-Held Units for Short-Range Wireless Biotelemetry; Wearable Wireless Telemetry System for Implantable BioMEMS Sensors; Electronic Escape Trails for Firefighters; Architecture for a High-to-Medium-Voltage Power Converter; 24-Way Radial Power Combiner/Divider for 31 to 36 GHz; Three-Stage InP Submillimeter-Wave MMIC Amplifier; Fast Electromechanical Switches Based on Carbon Nanotubes; Solid-State High-Temperature Power Cells; Fast Offset Laser Phase-Locking System; Fabricating High-Resolution X-Ray Collimators; Embossed Teflon AF Laminate Membrane Microfluidic Diaphragm Valves; Flipperons for Improved Aerodynamic Performance; System Estimates Radius of Curvature of a Segmented Mirror; Refractory Ceramic Foams for Novel Applications; Self-Deploying Trusses Containing Shape-Memory Polymers; Fuel-Cell Electrolytes Based on Organosilica Hybrid Proton Conductors; Molecules for Fluorescence Detection of Specific Chemicals; Cell-Detection Technique for Automated Patch Clamping; Redesigned Human Metabolic Simulator; Compact, Highly Stable Ion Atomic Clock; LiGa(OTf)(sub 4) as an Electrolyte Salt for Li-Ion Cells; Compact Dielectric-Rod White-Light Delay Lines; Single-Mode WGM Resonators Fabricated by Diamond Turning; Mitigating Photon Jitter in Optical PPM Communication; MACOS Version 3.31; Fiber-Optic Determination of N2, O2, and Fuel Vapor in the Ullage of Liquid-Fuel Tanks; Spiking Neurons for Analysis of Patterns; Symmetric Phase-Only Filtering in Particle-Image Velocimetry; Efficient Coupler for a Bessel Beam Dispersive Element; and Attitude and Translation Control of a Solar Sail Vehicle.

Source record↗

GPS Based Autonomous Navigation Study for the Lunar Gateway

This paper describes and predicts the performance of a conceptual autonomous GPS-based navigation system for NASA's planned lunar Gateway. This system is based on the flight-proven Magnetospheric Multiscale (MMS) GPS navigation system, augmented with an earth-pointed high-gain antenna and with an option for an atomic clock. High-fidelity simulations, calibrated against MMS flight data and making use of GPS transmitter patterns from the GPS Antenna Characterization Experiment (ACE) project are developed for operation of the system in the Gateway Near-Rectilinear Halo Orbit (NRHO). The results indicate that GPS can provide an autonomous, realtime navigation capability with comparable, or superior, performance to traditional Deep Space Network approach with eight hours of tracking per day.

Winternitz, Luke B.↗