Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular modeling”

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 109 records · Page 6

A proposal for the proper use of pseudopotentials in molecular orbital cluster model studies of chemisorption

The interaction of CO with Cu5, Ni5, and Al4 are treated as model systems for molecular adsorption on metal surfaces. The effect of the use of pseudopotentials for the metal atoms is studied by considering three types of clusters. In the first case, all of the metal electrons are explicitly included in the wave function; an all electron (AE) treatment. In the second case, the metal atom which directly interacts with the CO is described by AE but the remaining metal atoms include a pseudopotential for their core electrons. Finally, in the third case, all of the metal atoms in the cluster have a pseudopotential treatment for the core electrons. The AE cluster results are taken as reference values for the two pseudopotential treatments. The mixed cluster results are in excellent agreement with those of the all AE clusters; however, the results for the all pseudopotential cluster of Ni5CO or of Cu5CO are qualitatively different. The pseudopotential treatment for all of the metal atoms often leads to results that contain serious errors and it is not a reliable approach.

Bagus, P. S.↗

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will be maintained only if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74%–94% of the variance in aneuploid strains’ growth rates is explained by the cumulative cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of small nucleolar RNAs (snoRNAs) and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

genic load↗

Screening green solvents for multilayer plastic film recycling processes

Multilayer (ML) plastic films are essential packaging materials that help protect products from diverse external factors; however, only 5% of all ML films are recycled in the United States. Solvent-based technologies are a promising alternative for recycling ML films because they enable recovery of constituent polymer resins. For example, the Solvent Targeted Recovery and Precipitation (STRAPTM) process sequentially dissolves and separates polymer components using a series of targeted solvent washes. A crucial design aspect of this process is the impact of selected solvents on human health and on the environment. Here, this work introduces a computational framework that integrates molecular modeling, process modeling, techno-economic analysis (TEA), and life-cycle analysis (LCA) to quickly screen green solvents for solvent-based ML recycling processes. Initial screening for solvents based on selectivity is performed by estimating temperature-dependent solubilities using molecular-scale models. Subsequent screening uses basic estimates of energy use and octanol-water partition coefficients (logP) as key measures of health, safety, and environmental hazards. Detailed process modeling, TEA, and LCA are used on a reduced set of promising solvents identified in early screening steps to more accurately determine how solvent selection and associated operating conditions impact overall economics and environmental impacts. The framework is used for the identification of green solvents (from a database of 1,000 solvents) that separate an industrial ML film composed of polyethylene (PE), ethylene vinyl alcohol (EVOH), and polyethylene terephthalate (PET). Our analysis shows the effectiveness of the framework and reveals fundamental trade-offs between solvent greenness, solubility, and economics. Our work emphasizes the importance of taking a holistic systems view during solvent design and aims to inform the development of new processes for ML film recycling and the identification of new ML films that are easier to recycle.

economics↗

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗

Preparation, Fabrication, and Evaluation of Advanced Polymeric and Composite Materials

The thesis titles are given below: physical and mechanical behavior of amorphous poly(arylene ether-co-imidasole)s and poly(arylene ether-co-imidasole) modification epoxies; the requirements of patentability as applied to the chemical arts; fabrication of thermoplastic polymer composite ribbon; blend of reactive diluents with phenylethynyl-terminated arylene ether oligomers; the synthesis, characterization, and application of ether-containing polyimides; the synthesis of reflective and electrically conductive polyimide films via an in-situ self-metalization procedure using silver (I) complexes; the thermal cure of phenylethynyl terminated polyimides and selected model compounds; and the synthesis, characterization, and molecular modeling of cyclic arylene ether oligomers.

Orwoll, Robert A.↗

A full-coordinate model of the polymerase domain of HIV-1 reverse transcriptase and its interaction with a nucleic acid substrate

We present a full-coordinate model of residues 1-319 of the polymerase domain of HIV-I reverse transcriptase. This model was constructed from the x-ray crystallographic structure of Jacobo-Molina et al. (Jacobo-Molina et al., P.N.A.S. USA 90, 6320-6324 (1993)) which is currently available to the degree of C-coordinates. The backbone and side-chain atoms were constructed using the MAXSPROUT suite of programs (L. Holm and C. Sander, J. Mol. Biol. 218, 183-194 (1991)) and refined through molecular modeling. A seven base pair A-form dsDNA was positioned in the nucleic acid binding cleft to represent the template-primer complex. The orientation of the template-primer complex in the nucleic acid binding cleft was guided by the positions of phosphorus atoms in the crystal structure.

Non-NASA Center↗

All-Atom Modeling and Simulation of Biopolymer Interface: Dual Role of Antifouling Polymer Brushes

Antifouling polymer brushes are well-known for their exceptional resistance to unwanted protein adsorption. While experimental studies have extensively characterized protein–polymer brush interactions, computational investigations remain limited, largely due to the challenges in accurately modeling and integrating these complex, multicomponent systems without resorting to oversimplification. To address this challenge, this study presents one of the most comprehensive and realistic model systems to date, comprising the substrate, grafted polymers, and proteins. Our work interprets the interactions between polymer brush and protein based on realistic modeling without simplification. In particular, this study utilizes molecular modeling and simulation of polycationic and polyzwitterionic brushes─poly(dimethylaminoethyl methacrylate) (PDMAEMA), poly(2-(N-oxide-N,N-dimethylamino)ethyl methacrylate) (PNOMA), and poly(2-(N-3-sulfopropyl-N,N-dimethylammonium)ethyl methacrylate) (PSBMA)─grafted onto α-quartz substrates via polymerization initiator linkers. The brush models were developed to closely replicate experimentally synthesized samples and to provide detailed insights into structural and dynamical changes at the molecular level during protein adsorption. Using steered molecular dynamics simulations, we show that the PSBMA brush, due to its high local density, exhibits the greatest resistance to protein insertion. Cα root-mean-square deviation and interaction pattern analyses further reveal that the PSBMA brush also induces the most significant destabilization of lysozyme, while the PDMAEMA brush enhances protein stability through ion-mediated interactions. The PNOMA brush, while requiring the lowest force for protein adsorption, induces greater protein destabilization than the PDMAEMA brush, primarily due to electrostatic repulsion caused by a short carbon spacer length. Hydration analysis reveals that both the PSBMA brush and the lysozyme interacting with it exhibit the most rapid dehydration, attributed to the brush’s high local chain density, which results in the greatest lysozyme destabilization and the highest adsorption force. These findings highlight the dual role of antifouling polymer brushes: resisting protein adsorption and modulating protein structural dynamics. In conclusion, this study provides valuable insights for the rational design of next-generation antifouling materials and offers a framework for realistic model development in complex multicomponent systems.

Adsorption↗

On the coalescence-dispersion modeling of turbulent molecular mixing

The general coalescence-dispersion (C/D) closure provides phenomenological modeling of turbulent molecular mixing. The models of Curl and Dopazo and O'Brien appear as two limiting C/D models that bracket the range of results one can obtain by various models. This finding is used to investigate the sensitivtiy of the results to the choice of the model. Inert scalar mixing is found to be less model-sensitive than mixing accompanied by chemical reaction. Infinitely fast chemistry approximation is used to relate the C/D approach to Toor's earlier results. Pure mixing and infinite rate chemistry calculations are compared to study further a recent result of Hsieh and O'Brien who found that higher concentration moments are not sensitive to chemistry.

Givi, Peyman↗

Computer display and manipulation of biological molecules

This paper describes a computer model that was designed to investigate the conformation of molecules, macromolecules and subsequent complexes. Utilizing an advanced 3-D dynamic computer display system, the model is sufficiently versatile to accommodate a large variety of molecular input and to generate data for multiple purposes such as visual representation of conformational changes, and calculation of conformation and interaction energy. Molecules can be built on the basis of several levels of information. These include the specification of atomic coordinates and connectivities and the grouping of building blocks and duplicated substructures using symmetry rules found in crystals and polymers such as proteins and nucleic acids. Called AIMS (Ames Interactive Molecular modeling System), the model is now being used to study pre-biotic molecular evolution toward life.

Coeckelenbergh, Y.↗

Thermodynamics of a solar mixture of molecular hydrogen and helium at high pressure

The thermodynamic properties of a model molecular hydrogen and helium mixture are calculated in the strongly interacting region of 0.005 to 0.3 per cu cm for a range of temperatures that are of interest for the envelopes of the Jovian planets. Computed adiabats fit the gravity data and boundary conditions from model atmospheres of Jupiter.

Slattery, W. L.↗

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

Chemistry↗

Key Factors in Semi-Generic Coarse-Grained Modeling of Solid Polymer Electrolytes

Given the long length and time scales of interest and strong interactions present in solid polymer electrolytes, coarse-grained molecular models can be useful in understanding the molecular basis of their structure and transport properties. Rather than representing atoms individually, coarse-grained models use beads representing groups of atoms, making the system simpler and more efficient. Generic coarse-grained models, that are not built based on matching a particular atomistic system, can provide insight into the important physical considerations that apply across different chemical systems, though it can be unclear how to best match to experimental systems and capture relevant experimental behaviors with as simple of a model as possible. Here, we build on generic bead-spring models, but use stiff angle potentials, set different bead properties for different polymer types, and include additional ion parameters, creating a semi-generic coarse-grained model that can map more specifically to polymers and copolymers with different chain architectures, component glass transition temperatures, and ion solvation behaviors. We set most parameters with basic homopolymer data such as glass transition temperatures, Kuhn length, density, and dielectric constant, with further adjustment based on data from the polymer electrolyte system. We specifically model polystyrene-block-poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) with lithium triflate salt, considering the POEM to be made of a poly(methylmethacrylate) backbone with poly(ethylene oxide) side chains. Solvation of ions is accounted for by additional polymer-ion interactions of the form -S/r4, plus additional lithium-polymer Lennard-Jones interactions. We find these potentials have different effects and discuss strategies for setting these parameters.

Zhang, Yuanhao↗

Investigation on the Origins of Comets as Revealed through IR High Resolution Spectroscopy: I. Molecular Abundances

We report and analyze molecular abundances updatedin twenty comets by employingmodern data reduction procedures and molecular models. Using boxplots and scatterplots, we examine how the different molecular species are distributed among the comet population, while by means of pie charts we investigate the relative proportions of these molecular species in each comet. We compare these results with the orbital parametersof the selected targets in order to identify trends related to the dynamical history of each comet. In this way, we identify at least three chemical classes based mainly on relative abundances of CO, CH3OH, CH4, C2H6, HCN and NH3. The combination of relative abundances and orbital parameters is then compared with chemical models of planetary system formation and may offer an alternative approach to investigate the origin and evolution of the material in cometary nuclei.Our study also confirms the need to add more objects to our sample in order to improve the statistics, especially for hyper-volatiles (i.e. CH4and CO) in Jupiter Family comets.

M Lippi↗

High throughput, accurate gene annotation through AI and HPC-enabled structural analysis

With the advances in next generation sequencing technologies, the number of sequenced genomes is growing exponentially, resulting in a technology bottleneck for the translation of sequence information into usable hypotheses about the function of each gene. We have proposed leveraging our leadership high-performance computing (HPC) resources to help break this annotation bottleneck. Here we design an HPC-based framework to infer gene function from gene sequence by incorporating information about protein structure and interactions predicted by deep learning approaches. Accurate functional prediction and gene annotation using computational methods will facilitate breakthroughs in the genomic sciences essential to understanding and harnessing life processes in bacteria, fungi and plants. The development and applications of the state-of-the-art deep neural networks to protein structural modeling, interaction prediction, sequence comparison, and quality assessment of protein structural models will be made possible by leadership computational resources. These HPC-enabled bioinformatics and molecular modeling tools will provide powerful insights into molecular functions of genes.

59 BASIC BIOLOGICAL SCIENCES↗

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

The Earliest Ion Channels in Protocellular Membranes

Cellular membranes with their hydrophobic interior are virtually impermeable to ions. Bulk of ion transport through them is enabled through ion channels. Ion channels of contemporary cells are complex protein molecules which span the membrane creating a cylindrical pore filled with water. Protocells, which are widely regarded as precursors to modern cells, had similarly impermeable membranes, but the set of proteins in their disposal was much simpler and more limited. We have been, therefore, exploring an idea that the first ion channels in protocellular membranes were formed by much smaller peptide molecules that could spontaneously selfassemble into short-lived cylindrical bundles in a membrane. Earlier studies have shown that a group of peptides known as peptaibols is capable of forming ion channels in lipid bilayers when they are exposed to an electric field. Peptaibols are small, non-genetically encoded peptides produced by some fungi as a part of their system of defense against bacteria. They are usually only 14-20 residues long, which is just enough to span the membrane. Their sequence is characterized by the presence of non-standard amino acids which, interestingly, are also expected to have existed on the early earth. In particular, the presence of 2-aminoisobutyric acid (AIB) gives peptaibols strong helix forming propensities. Association of the helices inside membranes leads to the formation of cylindrical bundles, typically containing 4 to 10 monomers. Although peptaibols are excellent candidates for models of the earliest ion channels their structures, which are stabilized only by van der Waals forces and occasional hydrogen bonds between neighboring helices, are not very stable. Although it might properly reflect protobiological reality, it is also a major obstacle in studying channel behavior. For this reason we focused on two members of the peptaibol family, trichotoxin and antiamoebin, which are characterized by a single conductance level. This indicates that their structures are unique and stable. In addition, it is also believed that the trichotoxin channel displays some selectivity between potassium and chloride ions. This makes trichotoxin and antiamoebin ideal models of the earliest ion channels that could provide insight into the origins of ion conductance and selectivity. In the absence of crystal structure of the trichotoxin and antiamoebin channels, we propose their molecular models based on experimentally determined number of monomers forming the bundles. We use molecular dynamics simulations to validate the models in terms of their conductance and selectivity. On the basis of our simulations we show that the emergence of channels built of small, alpha-helical peptides was protobiologically plausible and did not require highly specific amino acid sequences, which is a convenient evolutionary trait. Despite their simple structure, such channels could possess properties that, at the first sight, appear to require markedly larger complexity. To this end, we will discuss how the amino acid sequence and structure of primitive channels give rise to the phenomena of ionic conductance and selectivity across the earliest cell walls, which were essential functions for the emergence and early evolution of protocells. Furthermore, we will argue that even though architectures of membrane proteins are not nearly as diverse as those of water-soluble proteins, they are sufficiently flexible to adapt readily to the functional demands arising during evolution.

Mijajlovic, Milan↗