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 1,441 records · Page 80

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗

Electronic Properties of Carbon Nanotubes and Junctions

Metallic and semiconducting Single Wall Carbon Nanotubes (CNT) have recently been characterized using scanning tunneling microscopy (STM) and the manipulation of individual CNT has been demonstrated. These developments make the prospect of using CNT as molecular wires and possibly as electronic devices an even more interesting one. We have been modeling various electronic properties such as the density of states and the transmission coefficient of CNT wires and junctions. These studies involve first calculating the stability of junctions using molecular dynamics simulations and then calculating the electronic properties using a pi-electron tight binding Hamiltonian. We have developed the expertise to calculate the electronic properties of both finite-sized CNT and CNT systems with semi-infinite boundary conditions. In this poster, we will present an overview of some of our results. The electronic application of CNT that is most promising at this time is their use as molecular wires. The conductance can however be greatly reduced because of reflection due to defects and contacts. We have modeled the transmission through CNT in the presence of two types of defects: weak uniform disorder and strong isolated scatterers. We find that the conductance is affected in significantly different manners due to these defects Junctions of CNT have also been imaged using STM. This makes it essential to derive rules for the formation of junctions between tubes of different chirality, study their relative energies and electronic properties. We have generalized the rules for connecting two different CNT and have calculated the transmission and density of states through CNT junctions. Metallic and semiconducting CNT can be joined to form a stable junction and their current versus voltage characteristics are asymmetric. CNT are deformed by the application of external forces including interactions with a substrate or other CNT. In many experiments, these deformation are expected to occur naturally. We will present some preliminary results of our calculations of the modification of CNT electronic properties as a result of deformations.

Anantram, M. P.↗

Examining Bioethanol-Producing Ultrastructures with Electron Microscopy and Molecular Dynamics

A longstanding goal in the biofuel sector is to increase the efficiency of the circularization of the use of materials (1). One such instance is the degradation of plant material, such as cellulose, and converting it to biofuel via bacterial digestion (2). The bifunctional aldehyde-alcohol dehydrogenase (AdhE) from the anaerobic thermophile C. thermocellum seems to be vital for the production and cellular tolerance of bioethanol; however, it lacks the efficiency to produce ethanol at industry standards (3). Therefore, to understand how C. thermocellum AdhE functions, we used cryo-electron microscopy (cryo-EM) to obtain a 3.2 A structure of the AdhE complex. When compared to previously published structures from E. coli (4-6), we identified potential regions that control the native conformation of the ultrastructure, as well as located channels that isolate the intermediate aldehyde from the cellular milieu. This high-resolution structure, in combination with molecular dynamics simulations, provides insight into one example of substrate channeling and establishes a basis for future mutagenesis studies.

AdhE↗

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗

Addition of Chloroform in a Solvent-in-Salt Electrolyte: Outcomes in the Microscopic Dynamics in Bulk and Confinement

Solvent-in-salt electrolytes (SISEs) are a promising alternative to the electrolytes currently used in commercial devices. Despite the SISEs’ advantages, their utilization is not yet realized due to the poor mobility of their chemical species. Here, we explore this problem by adding chloroform to a SISE formed by acetonitrile and a Li-salt. First, we performed illustrative cycling experiments to highlight the potential of this approach. Then, we focused on the description of the microscopic dynamics of the electrolytes and exposed the relevant aspects to be considered for their optimal performance. While the conductivity at low temperatures may be enhanced by the addition of chloroform, only subtle changes occur at room temperature. As revealed by molecular dynamics simulations and quasielastic neutron scattering (QENS) experiments, this effect is related to the preservation of the structure expected for a highly concentrated solution and promotion of the formation of ionic aggregates. These outcomes occur despite the increase in the overall mobility of the chemical species. The dynamics of the electrolytes in porous carbon was also investigated using QENS. In these circumstances, low concentrations of chloroform lead to diffusivities of the molecular species higher than those observed for the bulk electrolytes. As chloroform’s concentration increases, no further changes in the diffusivities are observed. Nonetheless, chloroform is mostly immobilized on the carbon surfaces and this behavior may be intensified at compositions closer to the eutectic mixture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗

Functional enzyme–polymer complexes

Significance The use of biological enzyme catalysts could have huge ramifications for chemical industries. However, these enzymes are often inactive in nonbiological conditions, such as high temperatures, present in industrial settings. Here, we show that the enzyme PETase (polyethylene terephthalate [PET]), with potential application in plastic recycling, is stabilized at elevated temperature through complexation with random copolymers. We demonstrate this through simulations and experiments on different types of substrates. Our simulations also provide strategies for designing more enzymatically active complexes by altering polymer composition and enzyme charge distribution.

59 BASIC BIOLOGICAL SCIENCES↗

Solvation Energy of Ions in a Stockmayer Fluid

We calculate the solvation energy of monovalent and divalent ions in various liquids with coarse-grained molecular dynamics simulations. Our theory treats the solvent as a Stockmayer fluid, which accounts for the intrinsic dipole moment of molecules and the rotational dynamics of the dipoles. Despite the simplicity of the model, we obtain qualitative agreement between the simulations and experimental data for the free energy and enthalpy of ion solvation, which indicates that the primary contribution to the solvation energy arises mainly from the first and possibly second solvation shells near the ions. Our results suggest that a Stockmayer fluid can serve as a reference model that enables direct comparison between theory and experiment and may be invoked to scale up electrostatic interactions from the atomic to the molecular length scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enzyme kinetics by GH7 cellobiohydrolases on chromogenic substrates is dictated by non‐productive binding: insights from crystal structures and MD simulation

Cellobiohydrolases (CBHs) in the glycoside hydrolase family 7 (GH7) ( EC3.2.1.176 ) are the major cellulose degrading enzymes both in industrial settings and in the context of carbon cycling in nature. Small carbohydrate conjugates such as p‐ nitrophenyl‐β‐ d ‐cellobioside (pNPC), p‐ nitrophenyl‐β‐ d ‐lactoside (pNPL) and methylumbelliferyl‐β‐ d ‐cellobioside have commonly been used in colorimetric and fluorometric assays for analysing activity of these enzymes. Despite the similar nature of these compounds the kinetics of their enzymatic hydrolysis vary greatly between the different compounds as well as among different enzymes within the GH7 family. Through enzyme kinetics, crystallographic structure determination, molecular dynamics simulations, and fluorometric binding studies using the closely related compound o‐nitrophenyl‐β‐ d ‐cellobioside (oNPC), in this work we examine the different hydrolysis characteristics of these compounds on two model enzymes of this class, TrCel7A from Trichoderma reesei and PcCel7D from Phanerochaete chrysosporium . Protein crystal structures of the E212Q mutant of TrCel7A with pNPC and pNPL, and the wildtype TrCel7A with oNPC, reveal that non‐productive binding at the product site is the dominating binding mode for these compounds. Enzyme kinetics results suggest the strength of non‐productive binding is a key determinant for the activity characteristics on these substrates, with PcCel7D consistently showing higher turnover rates ( k cat ) than TrCel7A, but higher Michaelis–Menten ( K M ) constants as well. Furthermore, oNPC turned out to be useful as an active‐site probe for fluorometric determination of the dissociation constant for cellobiose on TrCel7A but could not be utilized for the same purpose on PcCel7D, likely due to strong binding to an unknown site outside the active site.

59 BASIC BIOLOGICAL SCIENCES↗

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atomic Simulation of Complex DNA DSBs and the Interactions with the Ku70/80 Heterodimer

DNA double strand breaks (DSBs) induced by ionizing radiation (IR) usually contain modified bases such as 8-oxo-7,8-dihydroguanine (8-oxoG) and thymine glycol, apurinic/apyrimidinic (AP) sites, 2-deoxyribonolactone, or single-strand breaks (SSBs). The presence of such lesions in close proximity to the DSB terminus makes the DNA nicks more difficult to repair and rejoin than endogenously induced simple DSBs, and as such a major determinant of the biological effects of high linear energy transfer (LET) radiation as encountered in space travel. In this study we conducted molecular dynamics simulations on a series of DNA duplexes with various complex lesions of 8-oxoG and AP sites, in an effort to investigate the effects of such lesions to the structural integrity and stability of DNA after insulted by IR. We also simulated the interaction of such complex DSBs with the Ku70/80 heterodimer, the first protein in mammalian cells to embark the non-homologous end joining (NHEJ) DNA repair pathway. The results indicate, compared to DNA with simple DSBs, the complex lesions can enhance the hydrogen bonds opening rate at the DNA terminus, and increase the mobility of the whole duplex, thus they present more deleterious effects to the genome integrity if not captured and repaired promptly in cells. Simulations also demonstrate the binding of Ku drastically reduces structural disruption and flexibility caused by the complex lesions, and the interactions of Ku with complex DSBs have a different potential energy landscape from the bound structure with simple DSB. In all complex DSBs systems, the binding of DSB terminus with Ku70 is softened while the binding of the middle duplex with Ku80 is tightened. This energy shift may help the Ku protein to secure at the DSB terminus for a longer time, so that other end processing factors or repair pathways can proceed at the lesions before NHEJ repair process starts. These atomic simulations may provide valuable new insight into the selective action of repair proteins on damaged DNA.

Hu, Shaowen↗

Scaling Theory of Neutral Sequence-Specific Polyampholytes

The primary sequence of charged monomers in a polyampholyte has a profound effect on its conformational properties. In this work, theory and simulations are used to predict how the sequence influences the single-chain behavior of globally neutral polyampholytes (PAs) under salt-free conditions. We consider PAs with Markov statistics of charges, where charge “blockiness” is defined by the correlation λ along the chain. These PAs cover a wide spectrum of primary sequences ranging from alternating (λ = –1) to ideally random (λ = 0) to diblock (λ → 1). In a theta solvent, sufficiently long PAs of any primary sequence form globules, but their internal structure, as well as the position and sharpness of the coil–globule transition, crucially depends on λ. Further, using scaling arguments, we demonstrate that the excess charge of the globule correlation volume (blob) increases with increasing λ, thereby inducing stronger Coulomb attractions between neighboring blobs of opposite charge and leading to higher density PA globules. Scaling predictions are supported by the results of the random phase approximation (RPA) and coarse-grained molecular dynamics simulations. A generalization to good and poor solvents is presented that enables the construction of the resulting scaling diagram of PA conformations in terms of solvent quality and charge blockiness, thereby providing a comprehensive molecular understanding of the role of the PA primary sequence in the chain structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale . Here, our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Biomolecular condensates form spatially inhomogeneous network fluids

The functions of biomolecular condensates are thought to be influenced by their material properties, and these will be determined by the internal organization of molecules within condensates. However, structural characterizations of condensates are challenging, and rarely reported. Here, we deploy a combination of small angle neutron scattering, fluorescence recovery after photobleaching, and coarse-grained molecular dynamics simulations to provide structural descriptions of model condensates that are formed by macromolecules from nucleolar granular components (GCs). We show that these minimal facsimiles of GCs form condensates that are network fluids featuring spatial inhomogeneities across different length scales that reflect the contributions of distinct protein and peptide domains. The network-like inhomogeneous organization is characterized by a coexistence of liquid- and gas-like macromolecular densities that engenders bimodality of internal molecular dynamics. These insights suggest that condensates formed by multivalent proteins share features with network fluids formed by systems such as patchy or hairy colloids.

59 BASIC BIOLOGICAL SCIENCES↗