Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microscopic simulation”

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 127 records · Page 7

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic Monte Carlo Simulations of Electrochemical Oxidation and Reduction of Pt(111)

Kinetic Monte Carlo simulations of electrochemical oxidation and reduction are presented that match many features of the experimentally observed electrochemical and morphological response of Pt(111). Included in the simulation are all relevant microscopic transitions, including the formation of Pt-OH and Pt-O from Pt, surface diffusion of all three species, as well as an effective place exchange diffusion at high potential. A detailed description of this approach to modeling such a complex surface is also presented. Overall, it is found that many features of the Pt(111) CV, including hydroxylation, hysteresis, and surface roughening, can be correlated to events associated with n -coordinated surface species, such as the hydroxylation wave corresponding to a one-electron oxidation of 9-coordinated terrace sites. Oxidation to Pt-O species at potentials above 1.0 V are shown to correlate to the presence of growing surface roughness, and the simulations suggest the onset of Pt-O formation in steady-state cyclic voltammetry is dominated by the oxidation of 8-coordinated step edges rather than terrace sites. Implications for the stability of Pt(111) catalysts after thousands of voltammetric cycles are discussed.

36 MATERIALS SCIENCE↗

Electron optics for 4 MeV ultrafast electron microscope - Final Report

Calabazas Creek Research, Inc. (CCR), in collaboration with Northern Illinois University (NIU), proposed to leverage CCR’s experience in electron beam generation and focusing and NIU’s experience in electron microscopy to develop a 4 MeV UEM with sub-ns temporal resolution and 0.3 nm spatial resolution with a relative energy spread of 10-5. Important design tools for the program were CCR’s 3D, finite element code Beam Optics Analyzer (BOA) and the photocathode design code PHOTOQE. However, early in the program it became clear that information necessary to use PHOTOQE was not available, and the program focused on the modifying to BOA to provide the required functionality.

beam optics simulation↗

Sources of Emittance in RF Photocathode Injectors

Advances in electron beam technology have been central to creating the current generation of x-ray free electron lasers and ultra-fast electron microscopes. These once exotic devices have become essential tools for basic research and applied science. One important beam technology for both is the electron source which, for many of these instruments, is the photocathode RF gun. The invention of the photocathode gun and the concepts of emittance compensation and beam matching in the presence of space charge and RF forces have made these high-quality beams possible. Achieving even brighter beams requires a taking a finer resolution view of the electron dynamics near the cathode during photoemission and the initial acceleration of the beam. In addition, the high brightness beam is more sensitive to degradation by the optical aberrations of the gun’s RF and magnetic lenses. This paper discusses these topics including the beam properties due to fundamental photoemission physics, space charge effects close to the cathode, and optical distortions introduced by the RF and solenoid fields. Analytic relations for these phenomena are derived and compared with numerical simulations.

43 PARTICLE ACCELERATORS↗

Conservative tightly-coupled simulations of stochastic multiscale systems

Multiphysics problems often involve components whose macroscopic dynamics is driven by microscopic random fluctuations. The fidelity of simulations of such systems depends on their ability to propagate these random fluctuations throughout a computational domain, including subdomains represented by deterministic solvers. When the constituent processes take place in nonoverlapping subdomains, system behavior can be modeled via a domain-decomposition approach that couples separate components at the interfaces between these subdomains. Its coupling algorithm has to maintain a stable and efficient numerical time integration even at high noise strength. We propose a conservative domain-decomposition algorithm in which tight coupling is achieved by employing either Picard's or Newton's iterative method. Coupled diffusion equations, one of which has a Gaussian white-noise source term, provide a computational testbed for analysis of these two coupling strategies. Fully-converged (“implicit”) coupling with Newton's method typically outperforms its Picard counterpart, especially at high noise levels. This is because the number of Newton iterations scales linearly with the amplitude of the Gaussian noise, while the number of Picard iterations can scale superlinearly. At large time intervals between two subsequent inter-solver communications, the solution error for single-iteration (“explicit”) Picard's coupling can be several orders of magnitude higher than that for implicit coupling. Increasing the explicit coupling's communication frequency reduces this difference, but the resulting increase in computational cost can make it less efficient than implicit coupling at similar levels of solution error, depending on the communication frequency of the latter and the noise strength. This trend carries over into higher dimensions, although at high noise strength explicit coupling may be the only computationally viable option.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microscopic force driving the photoinduced ultrafast phase transition: Time-dependent density functional theory simulations of IrTe 2

Photoinduced phase transitions can have complex and intriguing behaviors more than material ground-state dynamics. Understanding the underlying mechanism can help us to design new ways to manipulate the materials. A variety of mechanisms has been proposed to explain the photoinduced phase transitions of IrTe 2 , but a consensus has yet to be reached. Here, we study the photo-induced phase transitions of IrTe 2 by performing the real-time time-dependent density functional theory (rt-TDDFT) simulations in combination with the occupation constrained DFT method. We reveal that the microscopic driving force for the photo-induced phase transitions arises from the tendency to lower the energy levels occupied by the excited carriers, through the increase or decrease of the associated atomic pair distances, depending on whether the newly occupied states are antibonding or bonding states, respectively. The geometric constraints between different bonds represented by the Poisson ratio can bring together different tendencies from different atomic pairs, thus forming a complex intriguing dynamic picture depending on the intensity of the excitation. Furthermore, we also find that phonons don't play a primary role, but can assist the phase transition. These findings imply that one can control the structural phase transitions by selectively exciting photocarriers into designated atomic states using appropriate photon sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

When SLIPS meets TIPS: An endogenous lubricant-infused surface by taking the diluent as the lubricant

Slippery lubricant-infused porous surfaces (SLIPSs) hold great promise in fields requiring adhesion resistance such as self-cleaning, anti-fouling, anti-icing, and anti-scaling. However, their road to practical application is still blocked by the complicated fabrication process and the restricted coating area. To address these challenges, we propose a new concept of endogenous SLIPS in which the porous structure grows from the solution coating of polypropylene/caprylyl methicone undergoing a thermally induced phase separation (TIPS) process. Caprylyl methicone serves as the diluent during TIPS as well as the lubricant in the SLIPS so that the coating can be fabricated in one step and unrestrictedly. Furthermore, the dissolution and phase separation mechanisms have been revealed by molecular dynamic simulation and in situ microscopic observation. Such SLIPSs exhibits impressive sliding property even for sticky and viscous liquids and exhibit extremely low ice adhesion in anti-icing tests.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic-Scale Dynamics of Five-Fold Twin Mediated Coalescence: Pathway-Dependent and Defect-Governed Nonclassical Growth Mechanisms

Defective crystals with distinct properties have been discovered in many systems. However, the growth mechanism of defective crystals is still poorly understood. Here, in this work, using a 5-fold twinned gold (Au) nanocrystal (NC) as a model system, three new coalescence pathways involving detwinning or twinning have been uncovered through atomic-scale dynamic observations in an aberration-corrected transmission electron microscope coupled with atomistic simulations. This demonstrates that beyond crystal size, coalescence growth dynamics involving 5-fold twins (5-FTs) are highly dependent on crystal defect density and the approach pathways of the crystals. When a 5-FT encounters a smaller 5-FT or a smaller NC in a face-to-face way, a new, larger 5-FT is produced at a relatively fast coalescence growth rate; while in a corner-to-corner way, the coalescence dynamics are more retarded and sluggish, which is conducive to the formation of complex multitwined structures rather than 5-FTs. This highlights that the planar defect density and crystal approach pathway influence the coalescence dynamics of crystals containing 5-FT. Moreover, a column-by-column grain boundary (GB) migration mechanism, which results in bent GBs, was also unveiled during the crystal coalescence process. These results enrich the general understanding of the crystal growth theory and provide new insights into the controllable fabrication of 5-FTs by crystal coalescence mechanisms.

crystallization↗

Giant electric field-induced second harmonic generation in polar skyrmions

Electric field-induced second harmonic generation allows electrically controlling nonlinear light-matter interactions crucial for emerging integrated photonics applications. Despite its wide presence in materials, the figures-of-merit of electric field-induced second harmonic generation are yet to be elevated to enable novel device functionalities. Here, we show that the polar skyrmions, a topological phase spontaneously formed in PbTiO3/SrTiO3 ferroelectric superlattices, exhibit a high comprehensive electric field-induced second harmonic generation performance. The second-order nonlinear susceptibility and modulation depth, measured under non-resonant 800 nm excitation, reach ~54.2 pm V −1 and ~664% V −1 , respectively, and high response bandwidth (higher than 10 MHz), wide operating temperature range (up to ~400 K) and good fatigue resistance (>10 10 cycles) are also demonstrated. Through combined in-situ experiments and phase-field simulations, we establish the microscopic links between the exotic polarization configuration and field-induced transition paths of the skyrmions and their electric field-induced second harmonic generation response. Our study not only presents a highly competitive thin-film material ready for constructing on-chip devices, but opens up new avenues of utilizing topological polar structures in the fields of photonics and optoelectronics.

Nonlinear optics↗

Observation of the spiral spin liquid in a triangular-lattice material

The spiral spin liquid (SSL) is a highly degenerate state characterized by a continuous contour or surface in reciprocal space spanned by a spiral propagation vector. Although the SSL state has been predicted in a number of various theoretical models, very few materials are so far experimentally identified to host such a state. Via combined single-crystal wide-angle and small-angle neutron scattering, we report observation of the SSL in the quasi-two-dimensional delafossite-like AgCrSe 2 . We show that it is a very close realization of the ideal Heisenberg J 1 –J 2 –J 3 frustrated model on the triangular lattice. By supplementing our experimental results with microscopic spin-dynamics simulations, we demonstrate how such exotic magnetic states are driven by thermal fluctuations and exchange frustration.

36 MATERIALS SCIENCE↗

Plug-and-Play Priors for Bright Field Electron Tomography and Sparse Interpolation

Many material and biological samples in scientific imaging are characterized by nonlocal repeating structures. These are studied using scanning electron microscopy and electron tomography. Sparse sampling of individual pixels in a two-dimensional image acquisition geometry, or sparse sampling of projection images with large tilt increments in a tomography experiment, can enable high speed data acquisition and minimize sample damage caused by the electron beam. In this paper, we present an algorithm for electron tomographic reconstruction and sparse image interpolation that exploits the nonlocal redundancy in images. We adapt a framework, termed plug-and-play priors, to solve these imaging problems in a regularized inversion setting. The power of the plug-and-play approach is that it allows a wide array of modern denoising algorithms to be used as a “prior model” for tomography and image interpolation. We also present sufficient mathematical conditions that ensure convergence of the plug-and-play approach, and we use these insights to design a new nonlocal means denoising algorithm. Finally, we demonstrate that the algorithm produces higher quality reconstructions on both simulated and real electron microscope data, along with improved convergence properties compared to other methods.

47 OTHER INSTRUMENTATION↗

Quantitative Nonlinear Optical Polarimetry with High Spatial Resolution

Nonlinear optical microscopy such as in the optical second-harmonic generation (SHG) modality has become a popular tool today for probing materials in the physical and biological sciences. While imaging and spectroscopy are widely used in the microscopy mode, nonlinear polarimetry, which can shed light on materials’ symmetry and microstructure, is relatively underdeveloped. This is partly because quantitative analytical modeling of the optical SHG response for anisotropic crystals and films largely assumes low-numerical aperture (NA) focusing of light, where the plane-wave approximation is sufficient. Tight focusing provides unique benefits in revealing out-of-plane polarization responses, which cannot be detected by near-plane-wave illumination at normal incidence. Here, we outline a method for quantitatively analyzing SHG polarimetry measurements obtained under high-NA focusing within a microscope geometry. Experiments and simulations of a variety of standard samples, from single crystals to thin films, are in good agreement, including measured and simulated spatial SHG maps of ferroelectric domains. A solution to the inverse problem is demonstrated, where the spatial distribution of an SHG tensor with unknown tensor coefficient magnitudes is determined by experimentally measured polarimetry. The ability to extract the out-of-plane component of the nonlinear polarization in normal incidence is demonstrated, which can be valuable for high-resolution polarimetry of 2D materials, thin films, heterostructures, and uniaxial crystals with a strong out-of-plane response.

36 MATERIALS SCIENCE↗

Final Scientific/Technical Report/DE-EE0006832

Lithium-sulfur batteries utilize a lithium metal anode and a sulfur-carbon composite cathode electrode. They are promising due to their theoretical high energy and power density. However in practice they suffer of relatively short lifetimes due to rapid degradation of the battery components. This project aimed to develop a better understanding of the issues and possible solutions that could result in an improvement of the performances of the lithium-sulfur batteries. In particular, the project focused on developing a better understanding of the mechanisms of battery discharge and charge, and of the chemistry taking place during these events. An integrated theoretical-experimental approach was utilized, where atomistic simulations informed about the microscopic behavior of the system, intermediate time and length scales were investigated with mesoscopic models, and experiments were carried out focusing on the synthesis and characterization of the cathode materials as well as buildup and tests of the complete electrochemical cell. This combined theoretical-experimental strategy saves costs by utilizing predictions from computational analyses to guide experimentation.

25 ENERGY STORAGE↗

Understanding and Strategies for Controlled Interfacial Phenomena in Lithium-Ion Batteries and Beyond

Electrolyte chemistry and properties and electrode structure and chemical properties particularly at the interfaces are crucial for the development of advanced battery components. This project focused on elucidating the roles of the electrolyte and electrode on the formation and evolution of the SEI layer and cell electrochemical performance on silicon and on Li metal anodes. For Si nano/microstructures we evaluated lithiation, volume changes, reactivity, and chemo-mechanical transformations as functions of nanoparticle size, shape, presence of coatings, electrolyte composition, electron leakage to the electrolyte, and cycling. For Li metal anodes we provided detailed analyses of the electrode and electrolyte effects on interfacial reactivity and as a function of microstructural evolution and underlying stochasticity as dendrites nucleate and grow. We thoroughly investigated and simulated how the chemistry of the various components of the electrolyte and the electrode architecture may affect the electrochemical reactions as well as cell degradation. A multiscale modeling approach was utilized, where atomistic simulations informed about the microscopic behavior of the system, and intermediate time and length scales were investigated with mesoscopic models. Results were also tested against selected experiments from our collaborators. The combined theoretical-experimental strategy saves costs by utilizing predictions from computational analyses to guide experimentation.

25 ENERGY STORAGE↗

Deep Learning-Based Queue-Aware Eco-Approach and Departure System for Plug-In Hybrid Electric Buses at Signalized Intersections: A Simulation Study

Eco-Approach and Departure (EAD) has been considered as a promising eco-driving strategy for vehicles traveling in an urban environment, where information such as signal phase and timing (SPaT) and geometric intersection description is well utilized to guide vehicles passing through intersections in the most energy-efficient manner. Previous studies formulated the optimal trajectory planning problem as finding the shortest path on a graphical model. While this method is effective in terms of energy saving, its computation efficiency can be further enhanced by adopting machine learning techniques. In this paper, we propose an innovative deep learning-based queue-aware eco-approach and departure (DLQ-EAD) system for a plug-in hybrid electric bus (PHEB), which is able to provide an online optimal trajectory for the vehicle considering both the downstream traffic condition (i.e. traffic lights, queues) and the vehicle powertrain efficiency. Based on optimal solutions obtained from the graph-based trajectory planning algorithm (GTPA), a deep neural network (DNN) is developed to learn the optimal vehicle speed for the next time step given its current state. It is demonstrated that the trained DNN can provide energy-efficient trajectories with high computational efficiency and high flexibility adopting to dynamic changes in the surrounding environment. To address the impact of downstream traffic, a queue prediction model is further developed using data from radars and connected vehicles (CVs), as well as signal timing data from SPaT messages. A comprehensive simulation study in the microscopic traffic modeling software PTV VISSIM shows that the proposed DLQ-EAD can achieve 18.7%-24.0% energy efficiency improvements for a single PHEB on various traffic congestion levels. The proposed queue prediction model can be of practical significance even at low penetration rates of CVs. Specifically, additional energy savings of 2.0%-8.2% can be further achieved with 20% vehicles in the network.

Ye, Fei↗

Characterization of coarse bainite transformation in low carbon steel during simulated welding thermal cycles

Coarse austenite to bainite transformation in low carbon steel under simulated welding thermal cycles was morphologically and crystallographically characterized by means of optical microscope, transmission electron microscope and electron backscattered diffraction technology. The results showed that the main microstructure changes from a mixture of lath martensite and bainitic ferrite to granular bainite with the increase in cooling time. The width of bainitic laths also increases gradually with the cooling time. For a welding thermal cycle with relatively short cooling time (e.g. t{sub 8/5} is 30 s), the main mode of variant grouping at the scale of individual prior austenite grains changes from Bain grouping to close-packed plane grouping with the progress of phase transformation, which results in inhomogeneous distribution of high angle boundaries. As the cooling time is increased, the Bain grouping of variants becomes predominant mode, which enlarges the effective grain size of product phase. - Highlights: • Main microstructure changes and the width of lath structure increases with cooling time. • Variant grouping changes from Bain zone to close-packed plane grouping with the transformation. • The change of variant grouping results in uneven distribution of high angle grain boundary. • Bain grouping is main mode for large heat input, which lowers the density of high angle boundary.

36 MATERIALS SCIENCE↗

Berni Julian Alder, theoretical physicist and inventor of molecular dynamics, 1925–2020

Berni Julian Alder, one of the leading figures in the invention of molecular dynamics simulations used for a wide array of problems in physics and chemistry, died on September 7th, 2020. His career, spanning more than 65 years, transformed statistical mechanics, many body physics, the study of chemistry and the microscopic dynamics of fluids, by making atomistic computational simulation (in parallel with traditional theory and experiment) a new pathway to unexpected discoveries. Among his many honors, the CECAM prize, recognizing exceptional contributions to the simulation of the microscopic properties of matter is named for him. He was awarded the National Medal of Science by President Obama in 2008.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗