Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) 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

Polymers in Deep Eutectic Solvents

The project investigated the behavior of polymers in ionic liquids such as deep eutectic solvents using an array of techniques. These included the development of atomistic force fields, coarse graining these force fields to the united atom level, large scale molecular dynamics (MD) simulations using new thermostat algorithms, and machine learning (ML) methods for the phase behavior. The project demonstrated the feasibility and accuracy of first principles force fields for ionic liquids, deep-eutectic solvents, and urea-water mixtures. A novel hierarchical coarse graining method was then used to develop accurate and efficient united-atom models, using which microsecond simulations were performed for polymers in ionic liquids. These simulations were in quantitative agreement with experiment, thus resolving previous controversies. Methods were also developed to obtain the potential of mean force between complex ions in solution. Finally, supervised ML methods were developed for the phase behavior of polymers in ionic liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rational Design of High-Performing Electrodes in Energy Storage Devices

Our approach in this project is to control conductivities of different ions by modulating the pore sizes in MOFs as sulfur hosts in cathodes. To validate our hypothesis that controlling the pore size of MOFs can effectively prevent LiPSs shuttling while allowing sufficient Li ion transport, which will ensure electrochemical cycling performance, we performed molecular dynamics (MD) simulations using a series of UiO MOFs with different linker sizes including LiPSs and electrolyte components (1M LiTFSI in DME/DOL 1:1 volume) in pores. Specifically, UiO-66, UiO-67, and UiO-68 are Zr-based MOFs with a linker length of 5.67 Å, 10.07 Å, and 14.35 Å, respectively, leading to different pore diameters of 8.4 Å, 12.8 Å, and 17.8 Å. We calculated the mean-squared displacements of Li and LiPS from MD simulations at 300 K for 50 ns or over to assess their mobility in MOFs.

25 ENERGY STORAGE↗

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

Establishing defect-property relationships for 2D-nanomaterials (Final technical report)

Studies of new families of two-dimensional nanomaterials (2DNMs) have established that they possess unique properties that diverge from those of their bulk counterparts. This project aimed to determine how tolerant 2DNMs are to extreme photon and particle fluxes and to identify the mechanisms governing their radiation response. Early work indicated that graphene is not as representative of other 2DNMs as previously assumed; however, the origin of this difference remained unclear. To address these knowledge gaps, this project investigated the physical processes occurring at multiple length scales in transition metal dichalcogenides (TMDs) and quantified their structural stability and property evolution under far-from-equilibrium conditions. In-situ and ex-situ ion and electron irradiations were performed to directly control and monitor defect formation in TMDs. Complementary density functional theory (DFT) and molecular dynamics (MD) simulations were employed to elucidate the mechanisms of defect generation and evolution. The outcomes of this work established mechanistic understanding of irradiation-induced defect formation in 2DNMs, quantified the radiation tolerance of TMDs, and built a fundamental knowledge base for correlating defect structures with material properties. Collectively, these results provide new insights into the stability of low-dimensional materials in extreme environments and enable the predictive design of radiation-tolerant 2DNMs.

2D materials↗

The Role of Surface Hydrophobicity on the Structure and Dynamics of CO2 and CH4 Confined in Silica Nanopores

Advancing a portfolio of technologies that range from the storage of excess renewable natural gas for distributed use to the capture and storage of CO 2 in geological formation are essential for meeting our energy needs while responding to challenges associated with climate change. Delineating the surface interactions and the organization of these gases in nanoporous environments is one of the less explored approaches to ground advances in novel materials for gas storage or predict the fate of stored gases in subsurface environments. To this end, the molecular scale interactions underlying the organization and transport behavior of CO 2 and CH 4 molecules in silica nanopores need to be investigated. To probe the influence of hydrophobic surfaces, a series of classical molecular dynamics (MD) simulations are performed to investigate the structure and dynamics of CO 2 and CH 4 confined in OH-terminated and CH 3 -terminated silica pores with diameters of 2, 4, 6, 8, and 10 nm at 298 K and 10 MPa. Higher adsorption extents of CO 2 compared to CH 4 are noted on OH-terminated and CH 3 -terminated pores. The adsorbed extents increase with the pore diameter. Further, the interfacial CO 2 and CH 4 molecules reside closer to the surface of OH-terminated pores compared to CH 3 -terminated pores. The lower adsorption extents of CH 4 on OH-terminated and CH 3 -terminated pores result in higher diffusion coefficients compared to CO 2 molecules. The diffusivities of both gases in OH-terminated and CH 3 -terminated pores increase systematically with the pore diameter. The higher adsorption extents of CO 2 on OH-terminated and CH 3 -terminated pores are driven by higher van der Waals and electrostatic interactions with the pore surfaces, while CH 4 adsorption is mainly due to van der Waals interactions with the pore walls. These findings provide the interfacial chemical basis underlying the organization and transport behavior of pressurized CO 2 and CH 4 gases in confinement.

Mohammed, Sohaib↗

The Role of Hydrophobic Nodes in the Dynamics of Class A β-Lactamases

Class A β-lactamases are known for being able to rapidly gain broad spectrum catalytic efficiency against most β-lactamase inhibitor combinations as a result of elusively minor point mutations. The evolution in class A β-lactamases occurs through optimisation of their dynamic phenotypes at different timescales. At long-timescales, certain conformations are more catalytically permissive than others while at the short timescales, fine-grained optimisation of free energy barriers can improve efficiency in ligand processing by the active site. Free energy barriers, which define all coordinated movements, depend on the flexibility of the secondary structural elements. The most highly conserved residues in class A β-lactamases are hydrophobic nodes that stabilize the core. To assess how the stable hydrophobic core is linked to the structural dynamics of the active site, we carried out adaptively sampled molecular dynamics (MD) simulations in four representative class A β-lactamases (KPC-2, SME-1, TEM-1, and SHV-1). Using Markov State Models (MSM) and unsupervised deep learning, we show that the dynamics of the hydrophobic nodes is used as a metastable relay of kinetic information within the core and is coupled with the catalytically permissive conformation of the active site environment. Finally, our results collectively demonstrate that the class A enzymes described here, share several important dynamic similarities and the hydrophobic nodes comprise of an informative set of dynamic variables in representative class A β-lactamases.

59 BASIC BIOLOGICAL SCIENCES↗

Elucidating the Mechanical Energy for Cyclization of a DNA Origami Tile

DNA origami has emerged as a versatile method to synthesize nanostructures with high precision. This bottom-up self-assembly approach can produce not only complex static architectures, but also dynamic reconfigurable structures with tunable properties. While DNA origami has been explored increasingly for diverse applications, such as biomedical and biophysical tools, related mechanics are also under active investigation. Here we studied the structural properties of DNA origami and investigated the energy needed to deform the DNA structures. We used a single-layer rectangular DNA origami tile as a model system and studied its cyclization process. This origami tile was designed with an inherent twist by placing crossovers every 16 base-pairs (bp), corresponding to a helical pitch of 10.67 bp/turn, which is slightly different from that of native B-form DNA (~10.5 bp/turn). We used molecular dynamics (MD) simulations based on a coarse-grained model on an open-source computational platform, oxDNA. We calculated the energies needed to overcome the initial curvature and induce mechanical deformation by applying linear spring forces. We found that the initial curvature may be overcome gradually during cyclization and a total of ~33.1 kcal/mol is required to complete the deformation. These results provide insights into the DNA origami mechanics and should be useful for diverse applications such as adaptive reconfiguration and energy absorption.

59 BASIC BIOLOGICAL SCIENCES↗

Microphysics of liquid water in sub-10 nm ultrafine aerosol particles

Abstract. Ultrafine aerosol particles with sizes smaller than 50 nm have been shown in recent studies to serve as a large source of cloud condensation nuclei that can promote additional cloud droplet formation under supersaturation conditions. Knowledge of the microphysics of liquid water in these droplets remains limited, particularly in the sub-10 nm particle size range, due to experimental and theoretical challenges associated with the complexity of aerosol components and the small length scales of interest (e.g., difficulty of precisely sampling the liquid–air interface, questionable validity of mean-field theoretical representations). Here, we carried out molecular dynamics (MD) simulations of aerosol particles with diameters between 1 and 10 nm and characterized atomistic-level structure and water dynamics in well-mixed and phase-separated systems with different particle sizes, NaCl salinities, and pimelic acid (PML) organic surface loadings as a function of distance from the time-averaged Gibbs dividing interface or instantaneous water–air interface. We define a sphericity factor (ϕ) that can shed light on the phase-mixing state of nanodroplets, and we reveal an unexpected dependence of mixing state on droplet size. Our results also evidence an ion concentration enhancement in ultrafine aerosols, which should modulate salt nucleation kinetics in sub-10 nm droplets, and provide detailed characterization of the influence of droplet size on surface tension and on water self-diffusivity near the interface. Analysis of water evaporation free energy and water activity demonstrates the validity of the Kelvin equation and Köhler theory at droplet sizes larger than 4 nm under moderate salinities and organic loadings and the need for further extension to account for ion concentration enhancement in sub-10 nm aerosols, droplet-size-dependent phase separation effects, and a sharp decrease in the cohesiveness of liquid water in sub-4 nm droplets. Finally, we show that an idealized fractional surface coating factor (fs) can be used to categorize and reconcile water accommodation coefficients (α*) observed in MD simulations and experimental results in the presence of organic coatings, and we resolve the droplet size dependence of α*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Dynamics study of Pb overlayer on Cu(100)

Isothermal-isobaric Molecular Dynamics (MD) simulation of a submonolayer Pb film in c(2x2) ordered structure adsorbed on a Cu(100) substrate showed retention of order to high T. The Embedded Atom Method (EAM) calculated the energy of atoms of overlayer and substrate. The time-averaged squared modulus of the two dimensional structure factor for the Pb overlayer measured the order of the overlayer. The results are for increasing T only, and require verification by simulated cooling.

Karimi, M.↗

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.↗

Thermal Expansion and Diffusion Coefficients of Carbon Nanotube-Polymer Composites

Classical molecular dynamics (MD) simulations employing Brenner potential for intra-nanotube interactions and van der Waals forces for polymer-nanotube interface have been used to investigate thermal expansion and diffusion characteristics of carbon nanotube-polyethylene composites. Addition of carbon nanotubes to polymer matrix is found to significantly increase the glass transition temperature Tg, and thermal expansion and diffusion coefficients in the composite above Tg. The increase has been attributed to the temperature dependent increase of the excluded volume for the polymer chains, and the findings could have implications in the composite processing, coating and painting applications.

Wei, Chengyu↗

Observation and Modeling of Single Wall Carbon Nanotube Bend Junctions

Single wall carbon nanotube (SWNT) bends, with diameters from approx. 1.0 to 2.5 nm and bend angles from 18 deg. to 34 deg., are observed in catalytic decomposition of hydrocarbons at 600 - 1200 C. An algorithm using molecular dynamics simulation (MD) techniques is developed to model these structures that are considered to be SWNT junctions formed by topological defects (i.e. pentagon-heptagon pairs). The algorithm is used to predict the tube helicities and defect configurations for bend junctions using the observed tube diameters and bend angles. The number and arrangement of the defects at the junction interfaces are found to depend on the tube helicities and bend angle. The structural and energetic calculations using the Brenner potential show a number of stable junction configurations for each bend angle with the 34 deg. bends being more stable than the others. Tight binding calculations for local density of state (LDOS) and transmission coefficients are carried out to investigate electrical properties of the bend junctions.

Han, Jie↗

Analysis of Carbon Nanotube Pull-out from a Polymer Matrix

Molecular dynamics (MD) simulations of carbon nanotube (NT) pull-out from a polymer matrix are carried out. As the NT pull-out develops in the simulation, variations in the displacement and velocities of the NT are monitored. The existence of a carbon-ring-based period in NT sliding during pull-out is identified. Linear trends in the NT velocity-force relation are observed and used to estimate an effective viscosity coefficient for interfacial sliding at the NT/polymer interface. As a result, the entire process of NT pull-out is characterized by an interfacial friction model that is based on a critical pull-out force, and an analog of Newton's friction law used to describe the NT/polymer interfacial sliding.

Frankland, S. J. V.↗

Nanomechanics of Carbon and CxByNz Nanotubes: Via a Quantum Molecular Dynamics Method

Nanomechanics of single-wall C, BN and BC$_3$ and B doped C nanotubes under axial compression and tension are investigated through a generalized tight-binding molecular dynamics (GTBMD) and {\it ab-initio} electronic structure methods. The dynamic strength of BN, BC$_3$ and B doped C nanotubes for small axial strain are comparable to each other. The main difference is in the critical strain at which structural collapse occurs. For example, even a shallow doping with B lowers the value of critical strain for C nanotubes. The critical strain for BN nanotube is found to be more than that for the similar C nanotube. Once the structural collapse starts to occur we find that carbon nanotubes irreversibly go into plastic deformation regime via the formation of tetrahedral (four-fold coordinated) bonds at the location of sharp pinches or kinks. This finding is considerably different from the classical MD (molecular dynamics) simulation results known so far. The energetics and electronic densities of states of the collapsed structures, investigated with {\it ab-initio) methods, will also be discussed.

Srivastava, Deepak↗

Theory of Transport of Long Polymer Molecules through Carbon Nanotube Channels

A theory of transport of long chain polymer molecules through carbon nanotube (CNT) channels is developed using Fokker-Planck equation and direct molecular dynamics (MD) simulations. The mean transport or translocation time tau is found to depend on the chemical potential energy, entropy and diffusion coefficient. A power law dependence tau approx. N(sup 2)is found where N is number of monomers in a molecule. For 10(exp 5)-unit long polyethylene molecules, tau is estimated to be approx. 1micro-s. The diffusion coefficient of long polymer molecules inside CNTs, like that of short ones, are found to be few orders of magnitude larger than in ordinary silicate based zeolite systems.

Wei, Chenyu↗