Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular 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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Origin of short- and medium-range order in supercooled liquid Ge 3 Sb 2 Te 6 from ab initio molecular dynamics simulations

Phase-change materials such as Ge–Sb–Te compounds have attracted much attention due to their potential value in electrical data storage. In contrast to the amorphous and crystalline phases, supercooled liquids are far from being deeply understood despite their inevitable role in both amorphization and crystallization processes. To this end, we have studied the dynamics properties and structural characteristics of liquid and supercooled liquid Ge 3 Sb 2 Te 6 during the fast cooling process. As the temperature decreases, chemical bonds become more homogeneous, but coordination numbers of Ge, Sb and Te atoms change very little. Meanwhile, the structural order of short-range configuration is obviously enhanced. Further studies suggest that Ge-centered, Sb-centered and Te-centered configurations change to the more ordered defective octahedrons mainly by adjusting the bond-angle relationship and bond length, rather than just by changing the coordination environment. It is the more ordered octahedrons that promote the formation of medium-range order. Our findings provide a deep insight into the origin of local structural order in supercooled liquid Ge 3 Sb 2 Te 6 , which is of great importance for the comprehensive understanding of amorphization and crystallization processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Analysis of an all-solid state nanobattery using molecular dynamics simulations under an external electric field

Present Li-ion battery (LIB) technology requires strong improvements in performance, energy capacity, charging-time, and cost to expand their application to e-mobility and grid storage. Li-metal is one of the most promising materials to replace commercial anodes such as graphite because of its 10 times higher specific capacity. However, Li-metal has high reactivity with commercial liquid electrolytes; thus, new solid materials are proposed to replace liquid electrolytes when Li-metal anodes are used. We present a theoretical analysis of the charging process in a full nanobattery, containing a LiCoO 2 cathode, a Li 7 P2S 8 I solid-state electrolyte (SSE), a Li-metal anode as well as Al and Cu collectors for the cathode and anode, respectively. In addition, we added a Li 3 P/Li 2 S film as a solid electrolyte interphase (SEI) layer between the Li-anode and SSE. Thus, we focus this study on the SEI and SSE. We simulated the charging of the nanobattery with an external voltage by applying an electric field. We estimated temperature profiles within the nanobattery and analyzed Li-ion transport through the SSE and SEI. Here, we observed a slight temperature rise at the SEI due to reactions forming $PS_{3}^{–}$ and $P_{2}S_{7}^{4}$$^{–}$ fragments at the interfaces; however, this temperature profile changes due to the charging current under the presence of the external electric field ε = 0.75 V Å –1 . Without the external field, the calculated open-circuit voltage (OCV) was 3.86 V for the battery, which is within the range of values of commercial cobalt-based LIBs. This voltage implies a spontaneous fall of available Li-ions from the anode to the cathode (during discharge). The charge of this nanobattery requires overcoming the OCV plus an additional voltage that determines the charging current. Thus, we applied an external potential able to neutralize the OCV, plus an additional 1.6 V to induce the transport of Li + from the cathode up to the anode. Several interesting details about Li + transport paths through the SSE and SEI are discussed.

25 ENERGY STORAGE↗
DOE OSTI2022

High-throughput screening of tribological properties of monolayer films using molecular dynamics and machine learning

Monolayer films have shown promise as a lubricating layer to reduce friction and wear of mechanical devices with separations on the nanoscale. These films have a vast design space with many tunable properties that can affect their tribological effectiveness. For example, terminal group chemistry, film composition, and backbone chemistry can all lead to films with significantly different tribological properties. This design space, however, is very difficult to explore without a combinatorial approach and an automatable, reproducible, and extensible workflow to screen for promising candidate films. Here, using the Molecular Simulation Design Framework (MoSDeF), a combinatorial screening study was performed to explore 9747 unique monolayer films (116 964 total simulations) and a machine learning (ML) model using a random forest regressor, an ensemble learning technique, to explore the role of terminal group chemistry and its effect on tribological effectiveness. The most promising films were found to contain small terminal groups such as cyano and ethylene. The ML model was subsequently applied to screen terminal group candidates identified from the ChEMBL small molecule library. Approximately 193 131 unique film candidates were screened with approximately a five order of magnitude speed-up in analysis compared to simulation alone. The ML model was thus able to be used as a predictive tool to greatly speed up the initial screening of promising candidate films for future simulation studies, suggesting that computational screening in combination with ML can greatly increase the throughput in combinatorial approaches to generate in silico data and then train ML models in a controlled, self-consistent fashion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

Here, we develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗
DOE OSTI2025

Multiscale Molecular Dynamics Simulations: Accelerating Conformational Sampling of Biomolecular Systems by Iterating All-Atom and Coarse-Grained Simulations

We developed the atomistic-coarse-grained multiscale MD simulation method in the OpenMM simulation package by iterating between the all-atom (AA) and coarse-grained (CG) MD simulations to enhance the sampling of biomolecular conformations. As the free energy surfaces are flattened during CG MD simulations, we can accelerate the transitions between different low-energy conformations. The AA-CG-AA cycles are repeated, facilitating the accelerated sampling of biomolecular conformations at a CG level, while the finer atomistic interactions are refined with AA simulators.

Do, Hung Nguyen↗