Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular case study”

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 163 records · Page 9

Optimal adjustment sets for causal query estimation in partially observed biomolecular networks

Abstract Causal query estimation in biomolecular networks commonly selects a ‘valid adjustment set’, i.e. a subset of network variables that eliminates the bias of the estimator. A same query may have multiple valid adjustment sets, each with a different variance. When networks are partially observed, current methods use graph-based criteria to find an adjustment set that minimizes asymptotic variance. Unfortunately, many models that share the same graph topology, and therefore same functional dependencies, may differ in the processes that generate the observational data. In these cases, the topology-based criteria fail to distinguish the variances of the adjustment sets. This deficiency can lead to sub-optimal adjustment sets, and to miss-characterization of the effect of the intervention. We propose an approach for deriving ‘optimal adjustment sets’ that takes into account the nature of the data, bias and finite-sample variance of the estimator, and cost. It empirically learns the data generating processes from historical experimental data, and characterizes the properties of the estimators by simulation. We demonstrate the utility of the proposed approach in four biomolecular Case studies with different topologies and different data generation processes. The implementation and reproducible Case studies are at https://github.com/srtaheri/OptimalAdjustmentSet.

59 BASIC BIOLOGICAL SCIENCES↗

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗

Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. Here, we find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bipolar Redox‐Active Molecules in Non‐Aqueous Organic Redox Flow Batteries: Status and Challenges

Abstract In this Review, recent advances in the development of organic bipolar redox‐active molecules (BRMs) for redox‐flow batteries (RFBs) are discussed, with special emphasis on their utilization in non‐aqueous systems. BRMs are single molecules that can be used as both anolyte and catholyte. Their unique redox chemistry enables identical components to be present in RFBs, which has distinct advantages in storage, mitigating crossover and compound degradation. Recent scientific innovations have uncovered three main strategies in preparing these bipolar species: developing stable radicals that are capable of both positive and negative redox reactions, covalently bridging electron donor‐acceptor motifs on one single molecule, and forming biredox eutectic mixtures through physically mixing active components. These approaches not only enabled symmetric flow battery design, but also resulted in electroactive materials with improved properties such as solubility, stability, and cell voltage. Nevertheless, much is still to be done. Importantly, developing BRMs that are competitive with vanadium RFBs in all performance indicators (e. g. solubility >1.0 M, voltage output >2.0 V, and cycling longevity >10 000 cycles) remains a challenge. Mechanistic studies including molecular dynamics, plausible degradation processes, and electrode/solution interface are necessary to peer into key parameters that dominate their battery performance. Accordingly, we review the current status and suggest future opportunities in these areas.

Li, Min↗

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗

Spectral emission characteristics near 646 nm and plasma properties of laser-induced plasma of gaseous uranium hexafluoride

We present a detailed investigation of the temporal evolution of line and continuum emissions in laser-induced breakdown spectroscopy (LIBS) of gaseous uranium hexafluoride (UF 6 ), focusing on the spectral region near 646 nm. Spectral emission features, signal-to-background ratios (SBRs) of selected uranium lines, and spectral linewidths were examined under varying UF 6 pressures (15–60 Torr) and laser pulse energies (10–60 mJ). Higher pressures and pulse energies enhanced continuum emission and reduced SBRs but did not cause significant spectral congestion. Additional studies with the use of different laser systems, including nanosecond-pulsed Nd:YAG lasers (at fundamental and various harmonics) and a femtosecond-pulsed Ti:sapphire laser, revealed long-lived plasma continua in all cases. This persistent continuum is attributed to a pseudo-continuum from overlapping molecular emissions, as its intensity scales linearly with electron number density, deviating from the expected quadratic dependence of classical ion–electron interactions for free-free and free-bound continuum emission. Based on plasma persistence time, SBR, U II/U I intensity ratios, and electronic excitation temperature, no fundamental advantage was found for femtosecond-pulsed lasers over conventional nanosecond-pulsed ones for UF 6 enrichment assay with direct LIBS measurement.

Laser-induced breakdown spectroscopy↗

A new generation of effective core potentials from correlated and spin–orbit calculations: Selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with 4d, 5d, 6s, and 6p valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin–orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and SO splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin–orbit effects. The subsequent quantum Monte Carlo calculations point out the importance of accurate trial wave functions that, in some cases (mid-series transition elements), require treatment well beyond a single-reference.

Wang, Guangming (ORCID:0000000231327925)↗

Substrate-Dependent Study of Chain Orientation and Order in Alkylphosphonic Acid Self-Assembled Monolayers for ALD Blocking

For years, many efforts in area selective atomic layer deposition (AS-ALD) have focused on trying to achieve high-quality self-assembled monolayers (SAMs), which have been shown by a number of studies to be effective for blocking deposition. In this work, we show that in some cases where a densely packed SAM is not formed, significant ALD inhibition may still be realized. The formation of octadecylphosphonic acid (ODPA) SAMs was evaluated on four metal substrates: Cu, Co, W, and Ru. The molecular orientation, chain packing, and relative surface coverage were evaluated using near-edge X-ray absorption fine structure (NEXAFS), Fourier transform infrared (FTIR) spectroscopy, and electrochemical impedance spectroscopy (EIS). ODPA SAMs formed on Co, Cu, and W showed strong angular dependence of the NEXAFS signal whereas ODPA on Ru did not, suggesting a disordered layer was formed on Ru. Additionally, EIS and FTIR spectroscopy confirmed that Co and Cu form densely packed, “crystal-like” SAMs whereas Ru and W form less dense monolayers, a surprising result since W-ODPA was previously shown to inhibit the ALD of ZnO and Al 2 O 3 best among all the substrates. This work suggests that multiple factors play a role in SAM-based AS-ALD, not just the SAM quality. Therefore, metrological averaging techniques (e.g., WCA and FTIR spectroscopy) commonly used for evaluating SAMs to predict their suitability for ALD inhibition should be supplemented by more atomically sensitive methods. Finally, it highlights important considerations for describing the mechanism of SAM-based selective ALD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent oligomerization pathways facilitated by electrolyte additives during solid-electrolyte interphase formation

The solid–electrolyte interphase (SEI) layer formation is known to play an important role in determining the lifetime of lithium-ion batteries. A thin, stable SEI layer is linked to overall improved battery performance and longevity, however, the factors and mechanisms that lead to optimal SEI morphology and composition are not well understood. Inclusion of electrolyte additives (fluoroethylene carbonate, FEC; and vinylene carbonate, VC) is often necessary for improving SEI characteristics. To understand how these electrolyte additives impact SEI formation, we employed molecular dynamics (MD) and density functional theory (DFT) simulations to study the reaction networks and oligomerization pathways, respectively, for three systems containing ethylene carbonate (EC), a lithium ion, and FEC or VC. MD simulations suggest radical oligomerization pathways analogous to traditional oligomerization with nucleophilic alkoxide species via SN1 reaction mechanisms. Both S N 1 and S N 2 mechanisms were studied for all three systems using DFT. Oligomerization reactions were studied with both a standard alkoxide species and a ring-opened EC radical as the nucleophiles and EC, FEC, and VC as the electrophiles. For all cases, FEC and VC exhibited lower free energy barriers and more stable adducts when compared with EC. We conclude that one of the role of additives is to modify the oligomerization process of EC by introducing branching points (FEC) or termination points (VC).

25 ENERGY STORAGE↗

Structures of elemental potassium at terapascal pressures

Here we investigate the pressure-temperature phase diagram of elemental potassium (K) up to multiterapascal (TPa) pressures using ab initio random structure searching (AIRSS), discovering eight structural phase transitions beyond the already known double hexagonal close-packed (dhcp) structure. Starting at 1.15 TPa, K transitions from the dhcp structure and passes through a variety of close-packed hexagonal and trigonal phases (dhcp → P6 3 /mmc → P3m1 → P$\bar{6}m2$), which differ only in their stacking sequence along the c axis. At 2.55 TPa, K adopts the hexagonal close-packed (hcp) structure, and at 7.21 TPa transitions into a complex orthorhombic phase of Fddd symmetry with 64 atoms in its (conventional) unit cell, before assuming an Ibam structure at 8.60 TPa, and, eventually, the face centered cubic (fcc) structure at 23.6 TPa, which persists well into the petapascal (PPa) regime. Further, we calculate the full pressure-temperature dependence of the melting line of K through extensive molecular dynamics simulations. We study the evolution of the bonding topology of K with pressure, finding that K passes through two quasi-molecular phases featuring diatomic pairs (the Fddd and Ibam structures), before ultimately becoming a high pressure electride (HPE) in the fcc phase. The electron-phonon and superconducting properties of K at these extreme pressures are also investigated, where we find the critical temperature T c rises to a maximum of between 7.65 K and 15.70 K in the hcp phase, before eventually decreasing to essentially zero in the fcc phase. Our results fully elucidate the structural and electronic behavior of K under the most extreme conditions, and provide new case studies for multi-TPa dynamic compression experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reducing the matrix effect in mass spectral imaging of biofilms using flow-cell culture

The interactions between soil microorganisms and soil minerals play a crucial role in the formation and evolution of minerals and the stability of soil aggregates. Due to the heterogeneity and diversity of the soil environment, the under-standing of the functions of bacterial biofilms in soil minerals at the microscale is limited. A soil mineral-bacterial biofilm system was used as a model in this study, and it was analyzed by time-of-flight secondary ion mass spectrometry (ToF-SIMS) to acquire molecular level information. Static culture in multi-wells and dynamic flow-cell culture in microfluidics of biofilms were investigated. Our results show that more characteristic molecules of biofilms can be observed in SIMS spectra of the flow-cell culture. In contrast, biofilm signature peaks are buried under the mineral components in SIMS spectra in the static culture case. Spectral overlay was used in peak selection prior to performing Principal component analysis (PCA). Comparisons of the PCA results between the static and flow-cell culture show more pronounced molecular features and higher loadings of organic peaks of the dynamic cultured specimens. For example, fatty acids secreted from bacterial biofilm extracellular polymeric substance are likely to be responsible for biofilm dispersal due to mineral treatment up to 48 h. Such findings suggest that the use of microfluidic cells to dynamically culture biofilms be a more suitable method for reducing the matrix effect arisen from the growth medium and minerals as a perturbation fac-tor for improved spectral and multivariate analysis of complex mass spectral data in ToF-SIMS. These results show that the interaction mechanism between biofilms and soil minerals at the molecular level can be better studied using the flow-cell culture and advanced mass spectral imaging techniques like ToF-SIMS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Barkas effect in strongly magnetized plasmas

Strongly magnetized plasmas, which are characterized by the particle gyrofrequency exceeding the plasma frequency, exhibit novel transport properties. For example, recent work showed that the friction force on a test charge moving through a strongly magnetized plasma not only consists of the typical stopping power component but also includes components perpendicular to the test charge's velocity. However, these studies only considered test charges that have the same sign as the charge of the plasma particles. Here, we extend these calculations to the case of charges with opposite signs (such as an ion interacting with strongly magnetized electrons). This is done with both a novel generalized Boltzmann kinetic theory and molecular dynamics simulations. It is found that the friction force changes dramatically depending on the sign of the interacting charges. Likewise, the stopping power component for oppositely charged particles decreases in magnitude compared with like-charged particles, and the perpendicular components increase in magnitude. Moreover, the difference between the two cases increases as the gyrofrequency becomes larger compared with the plasma frequency. Finally, the electrical resistivity is calculated from the friction force, where it is found that strong magnetization in conjunction with oppositely charged interactions significantly decreases the parallel resistivity and increases the perpendicular resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A density functional theory analysis of the adsorption and surface chemistry of inorganic iodine species on graphitea

In the event of a nuclear accident, fission products may be released into the environment. The release of 131 I is of particular concern to human health. Iodine can be captured using a number of materials and frequently, this is accomplished with activated carbon impregnated with organic bases. Previous studies have used DFT and the graphite (0001) surface as a surrogate for adsorption, those studies focus on the species I • , I 2 , and CH 3 I. In this work we perform an ab initio study of the adsorption onto the surface of a graphite sheet of I 2 , CH 3 I, and inorganic acidic iodine species (HI, HOI, HIO 2 , and HIO 3 ), which were selected to examine the possible effect of oxidation state on adsorption. The PBE exchange-correlation functional with D3 dispersion was employed. It was found that for molecular iodine, the iodine atoms tended to either situate above the center of a hexagonal site on the graphite or directly atop a carbon atom with the lighter components resting closer to the graphite. For each species the relative binding energies spanned the range of 21–33 kJ mol -1 and graphite-iodine distance was in the range of 3.52–3.93 Å. In all cases we found no significant charge transfer between the iodine species and the graphite, thus we conclude that all the iodine species studied undergo strong physisorption to the graphite.

Ritzmann, Andrew M.↗

Transient Symmetry Controls Photo Dynamics near Conical Intersections

Excited-state chemistry lacks generalised symmetry rules. With many femtochemistry studies focused on individual cases, it is hard to build up the same level of chemical intuition for excited-states as exists for ground states. Here, we unravel the degrees of freedom involved in ultrafast internal conversion (IC) by mapping the vibrational coherence of the initial wavepacket, and the dependence on molecular symmetry in various cyclic tertiary amines. Molecular symmetry plays an important role in the preservation of vibrational coherence in the transit from one electronic state to another. We show here that it is sufficient for the molecule to simply have the possibility of a more symmetric structure to achieve the preservation of vibrational coherence. It can be transient and still lead to preservation. Furthermore, this finding provides an additional angle on how symmetry influences electronic transitions and an additional piece to the puzzle of establishing symmetry-based selection rules for excited-state processes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Unveiling the Interfacial Reconstruction Mechanism Enabling Stable Growth of the Delafossite PdCoO 2 on Al 2 O 3 and LaAlO 3

Delafossites, composed of noble metal (A + ) and strongly correlated sublayers (BO 2 – ), form natural superlattices with highly anisotropic properties. These properties hold significant promise for various applications, but their exploitation hinges on the successful growth of high-quality thin films on suitable substrates. Unfortunately, the unique lattice geometry of delafossites presents a significant challenge to thin-film fabrication. Different delafossites grow differently, even when deposited on the same substrate, ranging from successful epitaxy to complete growth suppression. These variations often lack a clear correlation to obvious causes like lattice mismatch. Unidentified stabilization mechanisms appear to enable growth in certain cases, allowing these materials to form stable thin films or act as buffer layers for subsequent delafossite growth. This study employs advanced scanning transmission electron microscopy techniques to investigate the nucleation mechanism underlying the stable growth of PdCoO 2 films on Al 2 O 3 and LaAlO 3 substrates grown via molecular-beam epitaxy. Our findings reveal the presence of a secondary phase within the substrate surface that stabilizes the films. This mechanism deviates from the conventional understanding of strain relief mechanisms at oxide heterostructure interfaces and differs significantly from those observed for Cu-based delafossites.

36 MATERIALS SCIENCE↗

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A fast computational framework for the design of solvent-based plastic recycling processes

Multicomponent plastics cannot be processed using mechanical recycling technologies, hindering efforts to deal with plastic waste. Multicomponent plastics include multilayer plastic films, which are widely used for food and healthcare packaging. Multilayer films combine several layers (potentially dozens) of different polymers to protect products from external factors (e.g., oxygen, water, temperature, shock, and light). Solvent-based separation processes have emerged as a promising alternative to recycle these complex materials. For instance, the Solvent-Targeted Recovery and Precipitation (STRAP TM ) process uses sequential solvent washes to selectively dissolve and separate constituent polymers from multicomponent plastic waste, including films. STRAP TM process design (separation sequence, type of solvents, and operating conditions) changes significantly depending on the design of the multilayer plastic film (e.g., number, types, and proportions of polymers). The ability to quickly quantify the economic and environmental benefits of diverse STRAP TM process designs is essential to accelerate the development of sustainable recycling processes and more recyclable multilayer film products. In this work, we present a fast computational framework that integrates molecular-scale models, process modeling, and techno-economic and life cycle analysis to quickly evaluate STRAP TM designs. The computational framework is general and can be used to study the processing of complex multilayer plastic waste streams that contain many layers. Furthermore, we highlight the different uses of the framework via targeted case studies.

Computational framework↗

Effect of ionizing radiation on PuO2 and surrogate oxide materials under storage conditions

The stockpile of plutonium separated from used nuclear fuel amounts to many 100s of tonnes. The majority of this plutonium is stored as PuO2 powder packaged in multi-can containers. The conditions during the packaging of reprocessed PuO2 for long term storage are closely controlled to limit water uptake as radiolytic decomposition of any adsorbed water will lead to the formation of both a potentially flammable atmosphere containing molecular hydrogen as well as reactive oxygen species which may be incorporated into the oxide phase. Consequently, the safety case for the long-term storage of PuO2 requires a complete understanding of the fundamental physical, chemical and materials degradation processes occurring at the water-oxide interface inside the storage canisters. This study investigates the effect of ionizing radiation on PuO2 and various surrogate oxides under atmospheres potentially encountered in storage canisters.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗