Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular data”

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,531 records · Page 85

Investigating PVC polymer–plasticizer interactions with atomistic MD simulations and potential of mean force calculations

For this work, atomistic molecular dynamics (MD) simulations coupled with potential of mean force (PMF) calculations were employed to investigate the interactions between PVC polymer chains containing 6–20 repeating units and various plasticizers, with the goal of identifying potential replacements for the toxic plasticizer di(2-ethylhexyl)phthalate (DEHP) used in blood bags. The selected plasticizers belong to various chemical families, such as orthophthalates, citrates, adipates, and the terephthalate DEHT. Both the polymer and the plasticizers lack ionizable groups, and their interactions are primarily governed by van der Waals and electrostatic forces. A model correlating PMF profiles with interaction forces was developed and validated across polyvinyl chloride (PVC) polymers of different lengths and all investigated plasticizers. This model provides insight into how structural variations in plasticizers influence their respective PMF values. The study ranks the investigated plasticizers based on binding affinity, identifying TOTM (tris(2-ethylhexyl) trimellitate, TEHTM), BTHC (butyryl trihexyl citrate, Citroflex B-6), ATHC (acetyl trihexyl citrate, Citroflex A-6, CA-6), and DEHT (bis(2-ethylhexyl)terephthalate, DOTP/DEHTP) as promising alternatives for further investigation. This comprehensive study encompasses PVC polymers of four different lengths and 14 plasticizers, with all data averaged over 30 independent simulations. The approach provides a deeper understanding of molecular interactions, enabling the tailoring of polymer–plasticizer systems for diverse applications, including extractables and leachables, with relevance spanning materials science to biomedical engineering, and also serves as a basis for developing coarse-grained simulation protocols. Overall, this study provides valuable insights for designing safer and more efficient plasticizer substitutes and is well supported by other studies.

Shet, Sai Athmeeya G. [Sri Sathya Sai Institute of↗

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗

Determining the Molecular Growth Mechanisms of Tetragonal Lysozyme Crystals

Studies of the growth of tetragonal lysozyme crystals employing atomic force microscopy (AFM) have shown the advantages of this technique in investigating the growth mechanisms of protein crystals [1]. The resolution of these studies was in the micron range, which revealed surface features such as the occurrence of dislocations and 2D nucleation islands, similar to those found in inorganic systems. They clearly showed that the crystals grew by these surface growth mechanisms. However, the studies also revealed some surprising features, such as bimolecular growth step heights and pronounced growth anisotropies on the (110) face, which could not be explained. In previous studies we employed Periodic Bond Chain (PBC) theory to tetragonal lysozyme crystal growth and found that the crystals were constructed by strongly bonded molecular chains forming helices about the 43 axes [2,3]. The helices were connected to each other with weaker bonds. The growth process was shown to proceed by the formation of these 43 helices, resulting in bimolecular growth steps on the (110) face. It was also shown to explain many other observations on tetragonal lysozyme crystal growth. Although PBC analysis is not a new technique [4], it has not been widely used as the mechanisms predicted from it could not be experimentally verified. In this study the growth process of these crystals was investigated, particularly for the (110) face, employing some newly developed high resolution AFM techniques. These techniques allowed individual lysozyme molecules on the crystal faces to be resolved and predictions from PBC analyses to be tested. The analyses had shown that of the two possible packing arrangements on (110) faces, only one would actually occur. Employing the first of the newly developed techniques, these faces were scanned by high resolution AFM. The resulting images were then compared with the theoretically constructed images for the two possible packing arrangements on the (110) face. The theoretical images were constructed by convolution of the crystal surface shape obtained from crystallographic data with the AFM tip shape. The comparison confirmed the prediction that the molecular packing arrangement of these faces corresponded to that for complete 43 helices. The second AFM technique that was developed was used to follow the growth process by measuring the dimensions of individual growth units on the (110) face. Linescans across a growth step, performed near the saturation limit of the crystals, allowed the growth unit dimensions to be measured. These measurements showed that growth on the (110) face proceeded by the formation of new 43 helices from the addition of at least tetramer units in the [110] direction. In the [001] direction growth proceeded by the addition of various aggregate units corresponding to the 4(sub 3) helices.

Li, Huayu↗

A dynamic subgrid-scale model for compressible turbulence and scalar transport

The dynamic subgrid-scale (SGS) model of Germano et al. (1991) is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

Moin, P.↗

QuantumScents: Quantum-Mechanical Properties for 3.5k Olfactory Molecules

Quantitative structure–odor relationships are critically important for studies related to the function of olfaction. Current literature data sets contain expert-labeled molecules but lack feature data. This paper introduces QuantumScents, a quantum mechanics augmented derivative of the Leffingwell data set. QuantumScents contains 3.5k structurally and chemically diverse molecules ranging from 2 to 30 heavy atoms (CNOS) and their corresponding 3D coordinates, total PBE0 energy, molecular dipole moment, and per-atom Hirshfeld charges, dipoles, and ratios. The authors demonstrate that Hirshfeld charges and ratios contain sufficient information to perform molecular classification by training a Message Passing Neural Network with chemprop (Heid, E.; et al. ChemRxiv, 2023, DOI: 10.26434/chemrxiv-2023-3zcfl) to predict scent labels. Finally, the QuantumScents data set is freely available on Zenodo along with the authors’ code, example models, and data set generation workflow (https://zenodo.org/doi/10.5281/zenodo.8239853).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗

Improved materials characterization for spacecraft applications

A material outgassing test consists of an effusion cell, mass spectrometer, and thermally controlled quartz crystal microbalance (TQCM). The material sample is placed in the effusion cell at a predetermined temperature, and the sample is outgassed in a vacuum of 0.000001 torr range. The total outgassing mass as a function of time is determined with the TQCM which is cooled to LN2 temperatures. Based on the view factor of the TQCM, the total outgassed mass can be determined. The technique can be verified by comparing the results to the actual mass loss of the sample which is determined by hanging the diffusion cell from the Cohn microbalance. The mass spectrometer can be used to determine if there are any low molecular weight components outgassing, such a nitrogen which does not condense on the TQCM. After the material outgasses to a point of saturation, the effusion cell is closed off, and the TQCM is slowly heated to allow the condensed film to be broken down into its components. As the temperature is increased, the components evaporate from the surface at a different rate and can be detected with the mass spectrometer. The relative amount of each component is found by the change in frequency of the TQCM as the component evaporates. Hence, from this test, not only is isothermal kinetic data obtained, but the primary components can also be identified in terms of their molecular weight.

Cull, R.↗

FTICR, NPOC, TN, and Moisture of Variably Inundated Sediment across 48 North American Rivers

This dataset supports a broader study examining drivers of organic matter chemistry in variably inundated hyporheic zone sediments and further linking that chemistry to biogeochemical rates. The dataset provides physical and chemical sediment data collected from 48 rivers across North America. There are related data that will be published, including surficial sediment temperature time series that can be used to infer the dynamics of inundation prior to the collection of sediments. This dataset is comprised of two folders: (1) ECA2_Sediment and (2) ECA2_SitePhotos. The ECA2_Sediment folder contains: (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) non-purgeable organic carbon (NPOC), total nitrogen (TN), and sediment moisture; (5) sediment sampling protocol, (6) readme, and (7) folder of high resolution characterization of organic matter from sediment extracts via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). The FTICR-MS files are .xml, and the data package includes instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All other files are .csv, .txt, or .pdf. The ECA2_SitePhotos folder contains a readme and site photographs taken in the field. All files are .txt or .jpg.

54 ENVIRONMENTAL SCIENCES↗

KRAS interaction with RAF1 RAS-binding domain and cysteine-rich domain provides insights into RAS-mediated RAF activation

The first step of RAF activation involves binding to active RAS, resulting in the recruitment of RAF to the plasma membrane. To understand the molecular details of RAS-RAF interaction, we present crystal structures of wild-type and oncogenic mutants of KRAS complexed with the RAS-binding domain (RBD) and the membrane-interacting cysteine-rich domain (CRD) from the N-terminal regulatory region of RAF1. Our structures reveal that RBD and CRD interact with each other to form one structural entity in which both RBD and CRD interact extensively with KRAS. Mutations at the KRAS-CRD interface result in a significant reduction in RAF1 activation despite only a modest decrease in binding affinity. Combining our structures and published data, we provide a model of RAS-RAF complexation at the membrane, and molecular insights into RAS-RAF interaction during the process of RAS-mediated RAF activation.

60 APPLIED LIFE SCIENCES↗

MASCOT: molecular gas depletion times and metallicity gradients – evidence for feedback in quenching active galaxies

ABSTRACT We present results from the first public data release of the MaNGA-ARO Survey of CO Targets (MASCOT), focusing our study on galaxies whose star formation rates and stellar masses place them below the ridge of the star-forming main sequence. In optically selected type 2 AGN/low-ionization nuclear emission regions (LINERs)/Composites, we find an empirical relation between gas-phase metallicity gradients ∇Z and global molecular gas depletion times $t_{\rm dep} = M_{\rm H_2} /{\rm SFR}$ with ‘more quenched’ systems showing flatter/positive gradients. Our results are based on the O3N2 metallicity diagnostic (applied to star-forming regions within a given galaxy), which was recently suggested to also be robust against emission by diffuse ionized gas (DIG) and LINERs. We conduct a systematic investigation into possible drivers of the observed ∇Z − tdep relation (ouflows, gas accretion, in situ star formation, mergers, and morphology). We find a strong relation between ∇Z or tdep and centralized outflow strength traced by the [O iii] velocity broadening. We also find signatures of suppressed star formation in the outskirts in AGN-like galaxies with long depletion times and an enhancement of metals in the outer regions. We find no evidence of inflows impacting the metallicity gradients, and none of our results are found to be significantly affected by merger activity or morphology. We thus conclude that the observed ∇Z–tdep relation may stem from a combination of metal redistribution via weak feedback, and a connection to in situ star formation via a resolved mass-metallicity–SFR relation.

Astronomy & Astrophysics↗

Limited impact of jet-induced feedback in the multi-phase nuclear interstellar medium of 4C12.50

4C12.50 (IRAS 13451+1232) at z = 0.122 is an ultraluminous infared radio galaxy that has often been proposed as a prime candidate for the link between ultraluminous infared galaxies and young radio galaxies. It is also an interesting target to investigate whether and how radio-induced feedback affects the evolution of galaxies in the early phases of radio activity. We study, in detail for the first time, the hot (≥1500 K) molecular gas in 4C12.50. The potential impact of the radio jet on this gas phase, as well as on the star formation activity, are investigated. We also study the ionised (including coronal) gas as traced by the near-infrared lines. Using near-infrared long slit spectroscopy obtained with EMIR on GTC and X-shooter on VLT, we analyse the emission line spectrum of the ionised, coronal, and, specially, the hot molecular gas in the western nucleus hosting the compact radio jet. Based on high spatial resolution ALMA CO(2–1) data, we also revise the location of 4C12.50 in the Kennicutt-Schmidt diagram in order to investigate whether star formation is suppressed. 4C12.50 hosts (2.1 ± 0.4) × 10 4 M ⊙ of hot molecular gas. An unusually high rotational temperature T rot = 3020 ± 160 K is inferred. The molecular gas mass obeys a power-law temperature distribution, dM H2 /dT ∝ T –5 , from T ~300 K and up to ~3000 K. Both results support the idea that shocks (probably induced by the radio jet) contribute to the heating and excitation of the hot molecular gas. A molecular outflow is not detected. The coupling of the outflowing ionised and neutral outflows with the hot molecular gas is poor. Contrary to other studies, we claim that there is no evidence for star formation suppression in this object. If radio-induced feedback can regulate the star formation activity in galaxies, 4C12.50 is a promising candidate to reveal this phenomenon in action. However, we find no solid evidence for a current or past impact of this mechanism on the evolution of this system, neither by clearing out the dusty central cocoon efficiently, nor by suppressing the star formation activity.

79 ASTRONOMY AND ASTROPHYSICS↗

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗

The hTERT core promoter forms three parallel G-quadruplexes

The structure of the 68 nt sequence with G-quadruplex forming potential within the hTERT promoter is disputed. One model features a structure with three stacked parallel G-quadruplex units, while another features an unusual duplex hairpin structure adjoined to two stacked parallel and antiparallel quadruplexes. We report here the results of an integrated structural biology study designed to distinguish between these possibilities. As part of our study, we designed a sequence with an optimized hairpin structure and show that its biophysical and biochemical properties are inconsistent with the structure formed by the hTERT wild-type sequence. By using circular dichroism, thermal denaturation, nuclear magnetic resonance spectroscopy, analytical ultracentrifugation, small-angle X-ray scattering, molecular dynamics simulations and a DNase I cleavage assay we found that the wild type hTERT core promoter folds into a stacked, three-parallel G-quadruplex structure. The hairpin structure is inconsistent with all of our experimental data obtained with the wild-type sequence. All-atom models for both structures were constructed using molecular dynamics simulations. These models accurately predicted the experimental hydrodynamic properties measured for each structure. We found with certainty that the wild-type hTERT promoter sequence does not form a hairpin structure in solution, but rather folds into a compact stacked three-G-quadruplex conformation.

59 BASIC BIOLOGICAL SCIENCES↗

Genomic clocks and evolutionary timescales

For decades, molecular clocks have helped to illuminate the evolutionary timescale of life, but now genomic data pose a challenge for time estimation methods. It is unclear how to integrate data from many genes, each potentially evolving under a different model of substitution and at a different rate. Current methods can be grouped by the way the data are handled (genes considered separately or combined into a 'supergene') and the way gene-specific rate models are applied (global versus local clock). There are advantages and disadvantages to each of these approaches, and the optimal method has not yet emerged. Fortunately, time estimates inferred using many genes or proteins have greater precision and appear to be robust to different approaches.

Evolution, Molecular↗

Cell and molecular biology of simian virus 40: implications for human infections and disease

Simian virus 40 (SV40), a polyomavirus of rhesus macaque origin, was discovered in 1960 as a contaminant of polio vaccines that were distributed to millions of people from 1955 through early 1963. SV40 is a potent DNA tumor virus that induces tumors in rodents and transforms many types of cells in culture, including those of human origin. This virus has been a favored laboratory model for mechanistic studies of molecular processes in eukaryotic cells and of cellular transformation. The viral replication protein, named large T antigen (T-ag), is also the viral oncoprotein. There is a single serotype of SV40, but multiple strains of virus exist that are distinguishable by nucleotide differences in the regulatory region of the viral genome and in the part of the T-ag gene that encodes the protein's carboxyl terminus. Natural infections in monkeys by SV40 are usually benign but may become pathogenic in immunocompromised animals, and multiple tissues can be infected. SV40 can replicate in certain types of simian and human cells. SV40-neutralizing antibodies have been detected in individuals not exposed to contaminated polio vaccines. SV40 DNA has been identified in some normal human tissues, and there are accumulating reports of detection of SV40 DNA and/or T-ag in a variety of human tumors. This review presents aspects of replication and cell transformation by SV40 and considers their implications for human infections and disease pathogenesis by the virus. Critical assessment of virologic and epidemiologic data suggests a probable causative role for SV40 in certain human cancers, but additional studies are necessary to prove etiology.

Review↗

The 825-1110 angstroms EUV spectrum of Venus

On 15 August 1994 we launched the EUVS sounding rocket payload to observe the 825-1110 angstrom region of Venus's far ultraviolet airglow spectrum. The EUVS telescope/spectrograph obtained good data at five times higher spectral resolution than was previously available in the far ultraviolet. We present these data and compare our results to those obtained by the Galileo UVS and Venera 11/12 UV spectrophotometers. We identify several new spectral emission features, including both singly ionized nitrogen and molecular nitrogen in Venus's spectrum. We also see evidence for electron-impact-induced emission from CO. Finally, the EUVS data indicate that the "Ar" emissions detected in Venus's far ultraviolet spectrum by Venera 11/12 spectrophotometers are in fact not due to argon, thus eliminating the discrepancy between in situ and remote sensing measurements.

Flight Experiment↗

Transcriptomics-Based Points of Departure for Daphnia magna Exposed to 18 Per- and Polyfluoroalkyl Substances

Here, per- and polyfluoroalkyl substances (PFAS) represent a large group of contaminants of concern based on their widespread use, environmental persistence, and potential toxicity. Many traditional models for estimating toxicity, bioaccumulation, and other toxicological properties are not well suited for PFAS. Consequently, there is a need to generate hazard information for PFAS in an efficient and cost-effective manner. In the present study, Daphnia magna were exposed to multiple concentrations of 22 different PFAS for 24 h, in 96-well plate format. Following exposure, whole body RNA was extracted and extracts, each representing five exposed individuals, were subjected to RNA sequencing. Following analytical measurements to verify PFAS exposure concentrations, and quality control on processed cDNA libraries for sequencing, concentration-response modeling was applied to the data sets for 18 of the tested compounds, and the concentration at which a concerted molecular response occurred (transcriptomic point of departure; tPOD) was calculated. The tPODs, based on measured concentrations of PFAS, generally ranged from 0.03-0.58 µM (9.9-350 µg/L; interquartile range). In most cases, these concentrations were two orders of magnitude lower than similarly calculated tPODs for human cell lines exposed to PFAS. They were also lower than apical effect concentrations reported for seven PFAS for which some crustacean or invertebrate toxicity data were available, although there were a few exceptions. Despite being lower than most other available hazard benchmarks, Daphnia magna tPODs were, on average, four orders of magnitude greater than the maximum aqueous concentrations of PFAS measured in Great Lakes tributaries. Overall, this high throughput transcriptomics assay with Daphnia magna holds promise as a component of a tiered hazard evaluation strategy employing new approach methodologies.

59 BASIC BIOLOGICAL SCIENCES↗