Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Molecular models for creep in oriented polyethylene fibers

Highly oriented and crystalline polyetheylene (PE) fibers have a large failure stress under rapid tensile loading but exhibit significant creep at much smaller stresses that limits applications. A possible mechanism is slip of chains due to stress-enhanced, thermally activated nucleation of dislocations at chain ends in crystalline regions. Molecular dynamics simulations are used to parameterize a Frenkel–Kontorova model that provides analytic expressions for the limiting stress and activation energy for dislocation nucleation as a function of stress. Results from four commonly used hydrocarbon potentials are compared to show that the qualitative behavior is robust and estimate quantitative uncertainties. In all cases, the results can be described by an Eyring model with values of the zero-stress activation energy Ea0≈1.5 eV and activation volume V* ≈ 45 Å3 that are consistent with the experimental results for increasingly crystalline materials. The limiting yield stress is ∼8 GPa. These results suggest that activated dislocation nucleation at chain ends is an important mechanism for creep in highly oriented PE fibers.

O’Connor, Thomas C. (ORCID:0000000293930295)↗
DOE OSTI2025

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗
NASA NTRS1981

A molecular model for ice nucleation and growth, attachment 2

The central field empirical pair potential model is applied to studying the effects of kinks, ledges, and vacancies on the absorption of water molecules from the vapor. Molecular dynamics simulations indicate that cluster and/or surface modes play a primary role in the absorption process, the flexibility of the hydrogen bond serves to decrease the energy required for structural interconversion, and the rapid distribution of added energy in a hydrogen bonded system lead to aggregate stability which greatly exceeds that predicted by static energy calculations.

Hale, B. N.↗
NASA NTRS1981

A molecular model for ice nucleation and growth, attachment 1

The quantum mechanical technique is used to study ionic, configurational, and impurity defects in the ice surface. In addition to static calculations of the energetics of the water monomer-ice surface interactions, molecular dynamics studies were initiated. The calculations of the monomer-ice surface interaction, molecular dynamics studies were initiated. The calculations of monomer-ice surface interactions indicate that many adsorption sites exist on the ice surfaces and that the barriers between bonding sites are relatively low. Bonding on the prism face of ice is preferentially above lattice sites.

Plummer, P. L. M.↗
DOE OSTI2026

A High-Fidelity Molecular Model of the Cu(111) Repeating Unit

Dynamic processes at surfaces are central to heterogeneous catalysis, but their atomistic mechanism(s) can prove difficult to elucidate due to variations in material structure and the corresponding impact on reactivity. Moreover, disparities between reaction conditions and those employed for spectroscopic characterization at surfaces can inhibit detailed understanding of catalysis-relevant chemistries. Herein, we substantiate the so-called “cluster-surface” analogy by leveraging a low-valent tricopper architecture ( 1 ) as a model system for small molecule activation at Cu(111). Two reaction classes are explored: the adsorption of carbon monoxide (CO) and the dissociative adsorption of dihydrogen (H 2 ). These processes serve as an ideal testbed to compare the reactivity of a molecular cluster ( 1 ) to that of a heterogeneous surface, as both reactions have empirical data from measurements performed on crystalline Cu(111). Cluster 1 reversibly binds CO. Variable temperature NMR analysis with 13 CO reveals a favorable enthalpy but large negative entropy (−5.1 kcal × mol –1 and −22.9 cal × mol –1 × K –1 , respectively) for CO binding, affording a process that is marginally endergonic at room temperature (ΔG ads (298.15 K) = 1.7 ± 0.5 kcal × mol –1 ). Similarly, analogous to a Cu(111) surface, 1 is shown to oxidatively add (chemisorb) H 2 . Kinetic parameters were determined for this process and the activation enthalpy (8.4 ± 0.5 kcal × mol –1 ) closely mirrors that established for H 2 binding at the Cu(111) facet (6.0 to 12.4 kcal × mol –1 ). Together, these results showcase that a trinuclear cluster can reproduce the small molecule binding and activation energetics of a bulk crystalline surface, setting the stage for studying less-defined surface processes in an atomically precise molecular setting.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Molecular Modeling of the Poling of Piezoelectric Polyimides

The computational method described in this paper allows the calculation of the dielectric relaxation strength of an amorphous polymer based solely upon its chemical structure. The 4,4' oxydiphthalic anhydride (ODPA) dianhydride and bis-aminophenoxybenzene (APB) diamine based polyimides, (beta-CN) APB-ODPA and APB-ODPA were studied. Amorphous cells were constructed and then poled using molecular dynamics. Dielectric relaxation strengths of Delta(epsilon) = 17.8 for (beta-CN) APB-ODPA and Delta(epsilon) = 7.7 for APB-ODPA were predicted. These values are in excellent agreement with the experimental values. It was found that both the pendant nitrile dipole and the backbone anhydride residue dipole make significant contributions to the polyimides dielectric response. Specifically, it was shown that the difference in the magnitude of the dielectric relaxations is directly attributable to the nitrile dipole. The size of the relaxations indicate an absence of cooperative dipolar motions, The model was used to explain these results in terms of the average orientation of the nitrile and anhydride dipoles to within 51 deg. and 63 deg., respectively, of the applied electric field.

Young, J. A.↗
DOE OSTI2025

Molecular model of TFIIH recruitment to the transcription-coupled repair machinery

Transcription-coupled repair (TCR) is a vital nucleotide excision repair sub-pathway that removes DNA lesions from actively transcribed DNA strands. Binding of CSB to lesion-stalled RNA Polymerase II (Pol II) initiates TCR by triggering the recruitment of downstream repair factors. Yet it remains unknown how transcription factor IIH (TFIIH) is recruited to the intact TCR complex. Combining existing structural data with AlphaFold predictions, we build an integrative model of the initial TFIIH-bound TCR complex. We show how TFIIH can be first recruited in an open repair-inhibited conformation, which requires subsequent CAK module removal and conformational closure to process damaged DNA. In our model, CSB, CSA, UVSSA, elongation factor 1 (ELOF1), and specific Pol II and UVSSA-bound ubiquitin moieties come together to provide interaction interfaces needed for TFIIH recruitment. STK19 acts as a linchpin of the assembly, orienting the incoming TFIIH and bridging Pol II to core TCR factors and DNA. Molecular simulations of the TCR-associated CRL4CSA ubiquitin ligase complex unveil the interplay of segmental DDB1 flexibility, continuous Cullin4A flexibility, and the key role of ELOF1 for Pol II ubiquitination that enables TCR. Collectively, these findings elucidate the coordinated assembly of repair proteins in early TCR.

Paul, Tanmoy↗
DOE OSTI2020

Sorption-induced polymer rearrangement: approaches from molecular modeling

With the growing need for chemical separation and chemical storage solutions, polymeric adsorbents have emerged as a promising class of candidate materials because of their potentially tunable sorption properties, membrane structure and relatively cost consciousness. Moreover, the developing field of polymeric membrane materials has shown particular success at integrating both experimental and computational studies. However, these material systems are known to suffer from varying degrees of induced membrane structural rearrangement upon adsorbate uptake, and thus many polymeric membrane performance metrics are often considered to degrade with an increasing number of ‘guest’ species. In this mini–review, we highlight methodology tradeoffs and provide insights into atomistic molecular simulations used to study adsorption with flexible frameworks, which have the potential to predict separation, storage or catalytic capabilities a priori to experimental efforts. Specifically, molecular simulation methods that have been applied to provide predictions of polymeric membrane properties that have included consideration for sorbate–induced polymer chain rearrangement, swelling and/or plasticization are reviewed. Here, the examples and methodologies described provide demonstrations of the applicability of simulations as an approach to understand adsorption–based phenomena at an atomistic/molecular level, and as a tool to carry out screening studies aimed at efficiently providing analysis for a diversity of polymeric adsorbent–adsorbate systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Modeling to Increase Kraft Pulp Yield

Kraft pulping is an important component of the pulp and paper industry and is the predominant technology for removing lignin from wood carbohydrates. However, kraft pulping is energy-intensive, expensive, and is limited by the degradation of cellulose and hemicellulose. Pretreatment increases yield by stabilizing cellulose against degradation. However, protection of galactoglucomannan (GGM), the primary hemicellulose component of softwood, is minimal when conventional pretreatments are used. Here we investigate the effectiveness of new pretreatment methods on southern pine wood chips under a range of experimental conditions. If successful, improved pretreatment methods will increase carbohydrate yield, reduce waste, reduce energy use, lower the cost of bleaching, and decrease the cost of air emission controls. The purpose of this CRADA was to combine industrial expertise in wood pulping with national laboratory expertise in high-performance computing, leading to improved understanding of molecular-scale processes that limit carbohydrate yield during pretreatment and pulping. A combined computational and experimental approach was used to investigate pretreatment effectiveness under relevant pulping conditions and then use molecular simulation techniques to provide complementary insight into structural and chemical factors that govern the observed behavior. In this report we summarize the accomplishments of the project.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Predicting Melt Curves of Energetic Materials Using Molecular Models

Here, in this work, the solid–liquid coexistence curves of classical fully flexible atomistic models of α-RDX and β-HMX were calculated using thermodynamically rigorous methodologies that identify where the free energy difference between the phases is zero. The free energy difference between each phase at a given state point was computed using the pseudosupercritical path (PSCP) method, and Gibbs–Helmholtz integration was used to evaluate the solid–liquid free energy difference as a function of temperature. This procedure was repeated for several pressures to determine points along the coexistence curve, which were then fit to the Simon–Glatzel functional form. While effective, this method is computationally expensive. An alternative approach is to compute the melting point at a single pressure via the PSCP method, and then use the Gibbs–Duhem integration technique to trace out the coexistence curve in a more computationally economical manner. Both approaches were used to determine the coexistence curve of α-RDX. The Gibbs–Duhem integration method was shown to generate a melt curve that is in good agreement with the PSCP-derived melt curve, while only costing ~10 % of the computational resources used for the PSCP method. For α-RDX, the predicted melting temperature increases significantly more for a given increase in pressure when compared to available experimental data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Implementation and assessment of an extended hydrogenic molecular model in UEDGE

Abstract Two sets of pure deuterium plasmas are simulated using UEDGE, one including atoms and molecules and the other including atoms only as separate fluid species. Simulation results in one and two dimensions are reported to assess the role of molecules in tokamak plasmas. It is shown that thermal coupling of the molecules to the plasma can be relevant for the onset of detachment under highly collisional conditions in simple geometries, but play a modest role at low collisionality. Ion‐molecule equipartition presents an additional ion energy loss channel, dissipating ion energy that would otherwise heat the electrons by equipartition. The resulting reduced electron temperature increases the ion‐electron recombination by an order of magnitude as the plasma density is increased, yielding earlier and deeper detachment.

Holm, A.↗