Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular”

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 667 records · Page 37

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↗

Biological evaluation of Keggin‐type polyoxometalates on tyrosinase: Kinetics and molecular modeling

Abstract Abnormal overexpression of tyrosinase activity can lead to the production of hyperpigmentation in human skin and enzymatic browning in fruits and vegetables. Herein, the inhibition and mechanism of the H 3 PMo 12 O 40 and two transition metal‐substituted Keggin‐type polyoxometalates (Na 7 PMo 11 CoO 40 and Na 7 PMo 11 ZnO 40 ) on tyrosinase were studied by kinetics and molecular modeling. Kinetic studies indicated that all compounds had more potent inhibitory activities than standard arbutin, and H 3 PMo 12 O 40 (IC 50 = 0.443 ± 0.006 m m ) is ~15‐fold stronger inhibition than arbutin. Additionally, all compounds inhibited tyrosinase in a reversible competitive manner. Intriguingly, molecular modeling elucidated that three compounds competitively bind to tyrosinase mainly through more interactions with Cu 2+ ions and the amino acid residue capable of forming cation‐π and hydrogen bonding, forming a reversible non‐covalent complex. Molecular simulation study correlated well with the biological activity of three compounds in vitro. This work provided new insights into design and synthesis of polyoxometalates as tyrosinase inhibitors in the field of medicine, cosmetic, and food.

Chi, Guoxiang↗

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗

Accumulation of collagen molecular unfolding is the mechanism of cyclic fatigue damage and failure in collagenous tissues

Overuse injuries to dense collagenous tissues are common, but their etiology is poorly understood. The predominant hypothesis that micro-damage accumulation exceeds the rate of biological repair is missing a mechanistic explanation. Here, we used collagen hybridizing peptides to measure collagen molecular damage during tendon cyclic fatigue loading and computational simulations to identify potential explanations for our findings. Our results revealed that triple-helical collagen denaturation accumulates with increasing cycles of fatigue loading, and damage is correlated with creep strain independent of the cyclic strain rate. Finite-element simulations demonstrated that biphasic fluid flow is a possible fascicle-level mechanism to explain the rate dependence of the number of cycles and time to failure. Molecular dynamics simulations demonstrated that triple-helical unfolding is rate dependent, revealing rate-dependent mechanisms at multiple length scales in the tissue. The accumulation of collagen molecular denaturation during cyclic loading provides a long-sought “micro-damage” mechanism for the development of overuse injuries.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular mechanism of trehalose 6-phosphate inhibition of the plant metabolic sensor kinase SnRK1

SUCROSE-NON-FERMENTING1-RELATED PROTEIN KINASE1 (SnRK1), a central plant metabolic sensor kinase, phosphorylates its target proteins, triggering a global shift from anabolism to catabolism. Molecular modeling revealed that upon binding of KIN10 to GEMINIVIRUS REP-INTERACTING KINASE1 (GRIK1), KIN10’s activation T-loop reorients into GRIK1’s active site, enabling its phosphorylation and activation. Trehalose 6-phosphate (T6P) is a proxy for cellular sugar status and a potent inhibitor of SnRK1. T6P binds to KIN10, a SnRK1 catalytic subunit, weakening its affinity for GRIK1. Here, we investigate the molecular details of T6P inhibition of KIN10. Molecular dynamics simulations and in vitro phosphorylation assays identified and validated the T6P binding site on KIN10. Under high-sugar conditions, T6P binds to KIN10, blocking the reorientation of its activation loop and preventing its phosphorylation and activation by GRIK1. Under these conditions, SnRK1 maintains only basal activity levels, minimizing phosphorylation of its target proteins, thereby facilitating a general shift from catabolism to anabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. Furthermore, we identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ion Pairing, Clustering and Transport in a LiFSI-TMP Electrolyte as Functions of Salt Concentration using Molecular Dynamics Simulations

Battery capacity is highly related to ion-pairing mechanisms in electrolytes, since a cluster formation can lead to dead Li formation, reducing the number of charge carriers and leading to capacity fading. We use molecular dynamics simulations to model an electrolyte comprising trimethyl phosphate (TMP) solvent and a lithium bis(fluorosulfonyl)imide (LiFSI) salt, exploring effects of salt concentration on solvation and ion-transport. We simulate the LiFSI-TMP electrolyte for salt concentrations of 0.7, 1.43 and 3.82 molar. A statistical analysis was performed to study ion-pairing, clustering, diffusivity, conductivity, and coordination of Li-ions, providing insights into relations between molecular structures and transport properties. Molecular structure of ionic components changes as concentration increases, from a predominant solvent separated ion pair (SSIP) and contact ion pair (CIP) to aggregate salt (AGG) and ionic cluster formation. Given the formation of the ionic cluster, the diffusion mechanism followed by Li-ions changes from a hopping/exchange to a vehicular mechanism as concentration increases; this is reflected in a decrease of ionic conductivities. Ionicity was also calculated to reveal how the ionic motion changes from an uncorrelated to a correlated one as the salt concentration increases. Furthermore, we also compared our results with experimental calculations performed for similar electrolyte systems

25 ENERGY STORAGE↗

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗

Dashboard for Visualizing Molecular Property Prediction Machine Learning Results

This is a dashboard for exploring the results of machine learning models for predicting molecular properties from molecular structure. It includes tools for: 1. Modifying molecules to observe the change in predicted properties 2. Exploring the relationship between molecular structure and predicted properties 3. Recommending structurally similar molecules with improved properties 4. Exploring the impact of data subsampling on model performance metrics

Xu, Audrey↗

Molecular basis for the distinct cellular functions of the Lsm1–7 and Lsm2–8 complexes

Eukaryotes possess eight highly conserved Lsm (like Sm) proteins that assemble into circular, heteroheptameric complexes, bind RNA, and direct a diverse range of biological processes. Among the many essential functions of Lsm proteins, the cytoplasmic Lsm1–7 complex initiates mRNA decay, while the nuclear Lsm2–8 complex acts as a chaperone for U6 spliceosomal RNA. It has been unclear how these complexes perform their distinct functions while differing by only one out of seven subunits. Here, we elucidate the molecular basis for Lsm-RNA recognition and present four high-resolution structures of Lsm complexes bound to RNAs. The structures of Lsm2–8 bound to RNA identify the unique 2',3' cyclic phosphate end of U6 as a prime determinant of specificity. In contrast, the Lsm1–7 complex strongly discriminates against cyclic phosphates and tightly binds to oligouridylate tracts with terminal purines. Lsm5 uniquely recognizes purine bases, explaining its divergent sequence relative to other Lsm subunits. Lsm1–7 loads onto RNA from the 3' end and removal of the Lsm1 carboxy-terminal region allows Lsm1–7 to scan along RNA, suggesting a gated mechanism for accessing internal binding sites. These data reveal the molecular basis for RNA binding by Lsm proteins, a fundamental step in the formation of molecular assemblies that are central to eukaryotic mRNA metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Plasmon-phonon coupling between mid-infrared chiral metasurfaces and molecular vibrations

Plasmon-phonon coupling between metamaterials and molecular vibrations provides a new path for studying mid-infrared light-matter interactions and molecular detection. So far, the coupling between the plasmonic resonances of metamaterials and the phonon vibrational modes of molecules has been realized under linearly polarized light. Here, mid-infrared chiral plasmonic metasurfaces with high circular dichroism (CD) in absorption over 0.65 in the frequency range of 50 to 60 THz are demonstrated to strongly interact with the phonon vibrational resonance of polymethyl methacrylate (PMMA) molecules at 52 THz, under both left-handed and right-handed circularly polarized (LCP and RCP) light. The mode splitting features in the absorption spectra of the coupled metasurface-PMMA systems under both circular polarizations are studied in PMMA layers with different thicknesses. The relation between the mode splitting gap and the PMMA thickness is also revealed. The demonstrated results can be applied in areas of chiral molecular sensing, thermal emission, and thermal energy harvesting.

42 ENGINEERING↗

Deciphering the microbial and molecular responses of geographically diverse Setaria accessions grown in a nutrient-poor soil

The microbial and molecular characterization of the ectorhizosphere is an important step towards developing a more complete understanding of how the cultivation of biofuel crops can be undertaken in nutrient poor environments. The ectorhizosphere of Setaria is of particular interest because the plant component of this plant-microbe system is an important agricultural grain crop and a model for biofuel grasses. Importantly, Setaria lends itself to high throughput molecular studies. As such, we have identified important intra- and interspecific microbial and molecular differences in the ectorhizospheres of three geographically distant Setaria italica accessions and their wild ancestor S . viridis . All were grown in a nutrient-poor soil with and without nutrient addition. To assess the contrasting impact of nutrient deficiency observed for two S . italica accessions, we quantitatively evaluated differences in soil organic matter, microbial community, and metabolite profiles. Together, these measurements suggest that rhizosphere priming differs with Setaria accession, which comes from alterations in microbial community abundances, specifically Actinobacteria and Proteobacteria populations. When globally comparing the metabolomic response of Setaria to nutrient addition, plants produced distinctly different metabolic profiles in the leaves and roots. With nutrient addition, increases of nitrogen containing metabolites were significantly higher in plant leaves and roots along with significant increases in tyrosine derived alkaloids, serotonin, and synephrine. Glycerol was also found to be significantly increased in the leaves as well as the ectorhizosphere. These differences provide insight into how C 4 grasses adapt to changing nutrient availability in soils or with contrasting fertilization schemas. Gained knowledge could then be utilized in plant enhancement and bioengineering efforts to produce plants with superior traits when grown in nutrient poor soils.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular basis of P[II] major human rotavirus VP8* domain recognition of histo-blood group antigens

Initial cell attachment of rotavirus (RV) to specific cell surface glycan receptors, which is the essential first step in RV infection, is mediated by the VP8* domain of the spike protein VP4. Recently, human histo-blood group antigens (HBGAs) have been identified as receptors or attachment factors for human RV strains. RV strains in the P[4] and P[8] genotypes of the P[II] genogroup share common recognition of the Lewis b (Le b ) and H type 1 antigens, however, the molecular basis of receptor recognition by the major human P[8] RVs remains unknown due to lack of experimental structural information. Here, we used nuclear magnetic resonance (NMR) spectroscopy-based titration experiments and NMR-derived high ambiguity driven docking (HADDOCK) methods to elucidate the molecular basis for P[8] VP8* recognition of the Leb (LNDFH I) and type 1 HBGAs. We also used X-ray crystallography to determine the molecular details underlying P[6] recognition of H type 1 HBGAs. Unlike P[6]/P[19] VP8*s that recognize H type 1 HBGAs in a binding surface composed of an α-helix and a β-sheet, referred as the “βα binding site”, the P[8] and P[4] VP8*s bind Le b HBGAs in a previously undescribed pocket formed by the edges of two β-sheets, referred to as the “ββ binding site”. Importantly, the P[8] and P[4] VP8*s retain binding capability to non-Le b type 1 HBGAs using the βα binding site. The presence of two distinct binding sites for Le b and non-Le b HBGA glycans in the P[8] and P[4] VP8* domains suggests host-pathogen co-evolution under structural and functional adaptation of RV pathogens to host glycan polymorphisms. Assessment and understanding of the precise impact of this co-evolutionary process in determining RV host ranges and cross-species RV transmission should facilitate improved RV vaccine development and prediction of future RV strain emergence and epidemics.

59 BASIC BIOLOGICAL SCIENCES↗

Artificial Photosynthesis with Next Generation Molecular Catalysts

The objective of this project was to develop a new class of molecular photocatalysts for CO 2 reduction that will enable breakthrough advances in artificial photosynthesis. A major limitation of existing photocatalytic assemblies is that only a small fraction of absorbed photons are converted into fuel due to inefficiencies in electron transport between the different system components. In the case of molecular catalysts for reduction of CO 2 , highly reducing electron equivalents are required for catalysis. As a result, the electron flux between the photosensitizer and CO 2 reduction catalyst is heavily biased towards the photosensitizer. One method to direct the electron flux towards productive photocatalysis is to design molecular catalysts for CO 2 reduction that operate at more favorable potentials than existing catalysts. The conventional approach is to tune the potential of the catalyst through selection of the supporting ligands, however, this only leads to incremental improvements due to the existence of an unfavorable scaling relationship between the rate and potential of catalysis. Our hypothesis is that the catalytic rate and potential can be simultaneously enhanced using bimetallic catalysts in which one metal improves the potential of catalysis and the second metal increases the rate of catalysis. During the course of this project, a few studies reported some success in improving photocatalytic CO 2 reduction with heterobimetallic complexes, demonstrating the potential of this approach to catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the electron-water interaction at the molecular level: Integrating theory and experiment in the cluster regime

In this final technical report, we summarize our major accomplishments resulting from a combined experimental-theoretical approach that exploits the unique properties of molecular clusters to elucidate how local interactions at the molecular level drive the macroscopic behavior of materials and chemical transformations in liquids. The specific targets for our efforts have need to elucidate how charged entities are accommodated by hydrogen-bonding networks and to establish how the behavior of room temperature ionic liquids can be tailored for specific chemical applications. Because of their central importance in aqueous electrolyte chemistry, we have placed particular emphasis understanding the accommodation of excess electrons and protons by water. These surprisingly complex fundamental electrical charges play a fundamental role in a wide range of chemical and biological processes, and yet even a qualitative picture of their speciation in water and at the air-water interface has proven elusive. The main experimental diagnostic in our approach is to measure and analyze vibrational spectra of temperature and composition controlled molecular aggregates that isolate and amplify the spectral markers for local interactions that are often obscured in measurement on bulk or interfacial water. Our combined experimental-theoretical program has led to a greatly enhanced understanding of the nature of excess electrons and protons in water and in the process developed new and powerful methods of chemical analysis that are currently in use across the chemical sciences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Magnetism in North America Conference

The report provides information on the Molecular Magnetism in North America (MAGNA) conference that the DOE funding supported. It was held 2/21/2020 - 2/24/2020 at the King and Prince Hotel on St. Simon’s Island in southern Georgia. It brought together workers from academic institutions and government laboratories working in the field of molecular magnetism, including synthetic chemistry of molecular and 1D-3D molecule-based materials, and the various physical and spectroscopic methods for their study. The oral program spanned 3½ days and comprised 32 talks (30 mins each) by faculty and government scientists, and 12 talks (15 mins each) by junior people (graduate students and postdoctorals). The poster session took place at 6 pm on Saturday, Feb 22nd, and comprised 21 posters. Total registrations for the meeting were 68.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Aspects of Transport in Thin Films of Controlled Architecture (Final Technical Report)

The scientific and technological objectives of this project are motivated by problems in energy technologies that utilize differential transport in nanoscale confined volumes to achieve molecular separations. The core scientific phenomena that control observable behavior - wetting/dewetting, hydrophobicity, stochastic fluctuations in fluid flow, electrokinetics, etc. - are fundamentally different on the nanoscale, because the underlying forces that direct molecular motion are altered upon scaling to nanometer dimensions. We proposed to address this overarching problem by pursuing two objectives. Objective 1 targeted the development of single molecule-based electrokinetic, electrochemical and spectroelectrochemical probes of transport in one-dimensional (1D) nanostructures. Objective 2 focused on the interaction between wetting phenomena and nanoconfined flows under active (electrical and/or electrochemical) control. The studies were intended to develop the design rules, structural motifs, and operating principles needed to achieve singular control over molecular transport in nano-confined volumes, thereby impacting a broad spectrum of energy-relevant separations processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗