Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular evolution”

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

Model for photothermal ionization and molecular recombination during pulsed ablation of polyethylene

The surface erosion of polyethylene is studied under conditions characteristic of extremely high-rate heating expected in laser-driven and high energy density experiments. A coupled model for photoionization and volumetric photothermal degradation of polyethylene is developed that includes finite rate chemistry and molecular recombination of radicals. First, the model is used to explore the evolution of radicals and hydrocarbons during isochoric heating, for which it is found that polyethylene is not in chemical equilibrium for rates >108K/s. Then, the model is used to explore the cooperativity between photoionization and photothermal ionization during one-dimensional ablation from a pulsed heating source, accounting for coupled energy deposition, thermokinetics, thermochemistry, hydrodynamics, mechanics, and thermal conduction in the ALEGRA multi-physics code. The ablation depth per energy pulse is found to agree with an analytical model for instantaneous energy deposition and absent molecular recombination or photoionization. Otherwise, molecular recombination of radicals reduces the depth of ablation per pulse, owing to thermal stabilization on forming branched molecules. Radical formation through photoionization can be compensated for by recombination, but it accelerates ablation at high enough photoionization rates. Finally, the influence of molecular recombination on the ablation depth is found to decrease for long pulses, owing to decreased recombination reaction rates at lower vapor densities near the ablated surface. These and other findings highlight competition between the processes of scission, recombination, and photoionization during pulsed ablation of polyethylene.

Moore, N. W. (ORCID:0000000272259674)↗
DOE OSTI2020

Acetic Acid Adsorption and Reactions on Ni(110)

Acetic acid adsorption and reactions at multiple surface coverage values on Ni(110) were studied with temperature-programmed desorption (TPD) and infrared reflection absorption spectroscopy (IRAS) at 90–500 K. The experimental measurements were interpreted with density functional theory (DFT) calculations that provided information on adsorbate geometries, energies, and vibrational modes. Below the monolayer saturation coverage of 0.36 ML at 90 K, acetic acid adsorbs mostly molecularly. Above this coverage, a physisorbed layer is formed with dimers and catemers, without detectable monomers. Dimers and catemers desorb as molecular acetic acid at 157 and 172 K, respectively. Between 90 and 200 K, the O–H bond in acetic acid breaks to form bridge-bonded bidentate acetate that becomes the dominant surface species. Desorption-limited hydrogen evolution is observed at 265 K. However, even after the acetate formation, acetic acid desorbs molecularly at 200–300 K due to recombination. Minor surface species observed at 200 K, acetyls or acetates with a carbonyl group, decompose below 350 K and generate adsorbed carbon monoxide. At 350 K, the surface likely undergoes restructuring, the extent of which increases with acetic acid coverage. The initial dominant bridge-bonded bidentate acetate species formed below 200 K remain on the surface, but they now mostly adsorb on the restructured sites. Here, the acetates and all other remaining hydrocarbon species decompose simultaneously at 425 K in a narrow temperature range with concurrent evolution of hydrogen, carbon monoxide, and carbon dioxide. Above 425 K, only carbon remains on the surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Formation of I 1 stacking fault by deformation defect evolution from grain boundaries in Mg

I 1 stacking faults (SFs) in Mg alloys are regarded as the nucleation sites of $\langle \text{c+a} \rangle$ dislocations that are critical for these alloys to achieve high ductility. Previously it was proposed that the formation of I 1 SFs requires the accumulations of a large number of vacancies, which are difficult to achieve at low temperatures. In this study, molecular dynamics (MD) and molecular statics (MS) simulations based on empirical interatomic potentials were applied to investigate the deformation defect evolutions from the symmetric tilt grain boundaries (GBs) in Mg and Mg-Y alloys under external loading along $\langle c \rangle$-axis. The results show the planar faults (PFs) on Pyramidal I planes first appear due to the nucleation and glide of $\langle {\frac {1}{2}}c + p\rangle$ partial dislocations from GBs, where $\langle p\rangle$ = ${\frac {1}{3}} \langle10\bar{1}0\rangle$. These partial dislocations with pyramidal PFs interact with other defects, including pyramidal PFs themselves, GBs, and $\langle p\rangle$ partial dislocations, generating a large amount of I 1 SFs. Detailed analyses show the nucleation and growth of I 1 SFs are achieved by atomic shuffle events and deformation defect reactions without the requirements of vacancy diffusion. Our simulations also suggest the Y clusters at GBs can reduce the critical stress for the formation of pyramidal PFs and I 1 SFs, which provide a possible reason for the experimental observations that Y promotes the $\langle{\text{c+a}}\rangle$ dislocation activities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

BioPhotovoltaics: New paradigm towards high-efficiency and high-stability cells

In this project, we demonstrated significant progress in the development of Bio-Photovoltaic (BioPV) technology, with a particular focus on the transition from the initial success with Artemisinin (ART) to the development of the E1 compound. This journey began with the exploration of less conformationally restricted analogs of ART, leading to the discovery of E1. The initial success in the first quarter with ART set a precedent for the project, guiding our approach in molecular selection and design. Our computational studies provided a solid rationale for selecting specific biomolecules, with density functional theory calculations revealing the potential of certain molecules to form beneficial interactions with perovskite. This was a crucial step in narrowing down the candidate molecules from a broader selection. Subsequently, our approach involved simplifying these molecules to refine their properties and enhance their performance in bioPV applications. The ART-MAPbI3 films, for example, showcased not only high carrier mobility and hydrophobicity but also a significant increase in PCE. The evolution from ART to E1 was marked by a thorough understanding of molecular interactions and their impact on the material’s performance. This progression, from the complexity of lead candidates to the modeling and testing of simplified compounds, has culminated in the development of next-generation biomolecules with vastly improved properties. The link between E1 and ART, through this enhanced understanding, has been compelling and instrumental in achieving the milestones set forth in our project. The success in material and device performance underscores the importance of fundamental molecular design parameters, pointing towards future potential in the field of bioPV technology.

14 SOLAR ENERGY↗
DOE OSTI2024

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

X-ray and molecular dynamics study of the temperature-dependent structure of molten NaF-Zr⁢F 4

The local atomic structure of NaF-Zr⁢F 4 (53–47 mol%) molten system and its evolution with temperature are examined with x-ray scattering measurements which are then used to validate the quality of ab initio and neural network-based molecular dynamics (NNMD) calculations in the temperature range 515–700°⁢C. The machine-learning enhanced NNMD calculations offer improved efficiency while maintaining accuracy at higher distances compared to ab initio calculations. Looking at the evolution of the pair distribution function with increasing temperature, a fundamental change in the liquid structure within the selected temperature range, accompanied by a slight decrease in overall correlation is revealed. NNMD calculations indicate the coexistence of three different fluorozirconate complexes: [Zr⁢F 6 ] 2– , [Zr⁢F 7 ] 3– , and [Zr⁢F 8 ] 4– , with a shift in the dominant coordination state from the 7-coordinated Zr cation toward a 6-coordinated cation with increasing temperature. The study also highlights the metastability of different local coordination structures, with frequent interconversions between the 6- and 7-coordinate states. Analysis of the Zr-F-Zr angular distribution function reveals the presence of both “edge-sharing” and “corner-sharing” fluorozirconate complexes with specific bond angles and distances in accord with previous studies, while the next-nearest-neighbor cation-cation correlations demonstrate a clear preference for unlike cations as nearest-neighbor pairs, emphasizing nonrandom arrangement. Finally, these findings contribute to a comprehensive understanding of the complex local structure of the molten salt, providing insights into temperature-dependent preferences and correlations within the molten system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of Surface Orientation and Termination Plane on Glass‐to‐Crystal Transformation of Lithium Disilicate by Molecular Dynamics Simulations

Glass‐to‐crystal transformation of lithium disilicate is studied using molecular dynamics simulations using an effective partial charge potential. The structural evolution of the interface between glassy and crystalline lithium disilicate is analyzed to simulate crystallization of glass on pre‐existing crystal seeds. Besides previously used atomic number density, the distribution of Q n species (Si tetrahedra with n bridging oxygen) is shown to be an effective parameter for following this transformation quantitatively. The early stages of crystal growth are significantly affected by the orientation and termination of the surface of adjacent crystal, as indicated by calculated atomic density, partial ordering, atomic segregation, and an increase in Q 3 concentration. In particular, under‐coordinated Si within the outer crystal layer is found to be most effective in transforming the amorphous structure toward crystallinity. The increase in Q 3 in the glass close to interface region most clearly shows the initial stage of lithium disilicate crystal growth.

Sun, Wei↗
DOE OSTI2023

MD Simulation of Water Using a Rigid Body Description Requires a Small Time Step to Ensure Equipartition

In simulations of aqueous systems, it is common to freeze the bond vibration and angle bending modes in water to allow for a longer time step δt for integrating the equations of motion. Thus, δt = 2 fs is often used in simulating rigid models of water. We simulate the SPC/E model of water using δt from 0.5 to 3.0 fs and up to 4 fs using hydrogen mass repartitioning. In these simulations, we find that for all but δt = 0.5 fs, equipartition is not obtained between translational and rotational modes, with the rotational modes exhibiting a lower temperature than the translation modes. To probe the reasons for the lack of equipartition, we study the autocorrelation of the translational velocity of the center of mass and the angular velocity of the rigid water molecule, respectively. Here, we find that the rotational relaxation occurs on a timescale comparable to vibrational periods, calling into question the original motivations for freezing the vibrations. Furthermore, a time step with δt ≥ 1 fs is not able to capture accurately the fast rotational relaxation, which reveals its impact as an effective slowing-down of rotational relaxation. The fluctuation–dissipation relation then leads to the conclusion that the rotational temperature should be cooler for δt greater than the reference value of 0.5 fs. Consideration of fluctuation–dissipation in equilibrium molecular dynamics simulations also emphasizes the need to capture the temporal evolution of fluctuations with fidelity and the role of δt in this regard. The time step also influences the solution thermodynamic properties: both the mean system potential energies and the excess entropy of hydration of a soft repulsive cavity are sensitive to δt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Spatiotemporal Formation and Growth Kinetics of Polyelectrolyte Complex Micelles with Millisecond Resolution

We have directly observed the in situ self-assembly kinetics of polyelectrolyte complex (PEC) micelles by synchrotron time-resolved small-angle X-ray scattering, equipped with a stopped-flow device that provides millisecond temporal resolution. A synthesized neutral-charged diblock polycation and homopolyanion that we have previously investigated as a model charge-matched, core-shell micelle system were selected for this work. The initial micellization of the oppositely charged polyelectrolytes was completed within the dead time of mixing of 100 ms, followed by micelle growth and equilibration up to several seconds. By combining the structural evolution of the radius of gyration (R-g) with complementary molecular dynamics simulations, we show how the self-assemblies evolve incrementally in size over time through a two-step kinetic process: first, oppositely charged polyelectrolyte chains pair to form nascent aggregates that immediately assemble into spherical micelles, and second, these PEC micelles grow into larger micellar entities. Furthermore, this work has determined one possible kinetic pathway for the initial formation of PEC micelles, which provides useful physical insights for increasing fundamental understanding self-assembly dynamics, driven by polyelectrolyte complexation that occurs on ultrafast time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Step-by-step state-selective tracking of fragmentation dynamics of water dications by momentum imaging

Abstract The double photoionization of a molecule by one photon ejects two electrons and typically creates an unstable dication. Observing the subsequent fragmentation products in coincidence can reveal a surprisingly detailed picture of the dynamics. Determining the time evolution and quantum mechanical states involved leads to deeper understanding of molecular dynamics. Here in a combined experimental and theoretical study, we unambiguously separate the sequential breakup via D + + OD + intermediates, from other processes leading to the same D + + D + + O final products of double ionization of water by a single photon. Moreover, we experimentally identify, separate, and follow step by step, two pathways involving the b 1 Σ + and a 1 Δ electronic states of the intermediate OD + ion. Our classical trajectory calculations on the relevant potential energy surfaces reproduce well the measured data and, combined with the experiment, enable the determination of the internal energy and angular momentum distribution of the OD + intermediate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Electrical conductivities of (Mg,Fe)O at extreme pressures and implications for planetary magma oceans

During planet formation, planets undergo many impacts that can generate magma oceans. When these crystallize, part of the magma densifies via iron enrichment and migrates to the core–mantle boundary, forming an iron-rich basal magma ocean (BMO). The BMO could generate a dynamo in early Earth and super-Earths if the electrical conductivity of the BMO, which is thought to be sensitive to its Fe content, is sufficiently high. To test this hypothesis, here we conduct laser-driven shock experiments on ferropericlase (Mg x ,Fe 1−x )O (0.95 ≤ x ≤ 1) as an Fe-rich BMO analogue, perform density functional theory molecular dynamics simulations on MgO and calculate the long-term evolution of super-Earths. We find that the d.c. conductivities of MgO and (Mg,Fe)O are indistinguishable between 467 GPa and 1,400 GPa, despite previous predictions. Here, we predict that super-Earths larger than 3–6 Earth masses can produce BMO-driven dynamos that are almost one order of magnitude stronger than core-driven dynamos for several billion years.

Exoplanets↗
DOE OSTI2022

Strong Coulomb coupling influences ion and neutral temperatures in atmospheric pressure plasmas

Molecular dynamics simulations are used to model ion and neutral temperature evolution in partially-ionized atmospheric pressure plasma at different ionization fractions. Results show that ion–ion interactions are strongly coupled at ionization fractions as low as 10 –5 and that the temperature evolution is influenced by effects associated with the strong coupling. Specifically, disorder-induced heating is found to rapidly heat ions on a timescale of the ion plasma period (~10 s ps) after an ionization pulse. This is followed by the collisional relaxation of ions and neutrals, which cools ions and heats neutrals on a longer (~ns) timescale. Slight heating then occurs over a much longer (~100 s ns) timescale due to ion-neutral three-body recombination. An analytic model of the temperature evolution is developed that agrees with the simulation results. Finally, a conclusion is that strong coupling effects are important in atmospheric pressure plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Study of Microbial Physiology Under Microoxic Conditions Is Critical but Neglected

ABSTRACT During the early evolution of life on Earth, the environment was largely free of molecular oxygen, and only anaerobic life existed. With the subsequent oxidation of oceans and the atmosphere, a wide range of environmental niches, ranging from anoxic to microoxic/hypoxic and oxic, developed. Despite this broad range of natural environments, microbiology as a field has focused on the physiology, metabolism, and genetics of aerobic microorganisms, with less attention paid to anaerobes and much less attention paid to microaerophiles. The disparity in studies between aerobic and anaerobic conditions is rampant in host‐associated systems, particularly in human health, and studies of microorganisms in intermediate oxygen conditions between fully aerobic and fully anoxic conditions are exceedingly rare. Studies on the physiological behaviour, metabolism, growth response, and drug susceptibility patterns of commensal and pathogenic organisms are almost totally neglected in microoxic conditions. Furthermore, microorganisms from microaerobic and microoxic ecosystems have been less robustly explored in terms of physiology, growth, and metabolism. In this work, we highlight the importance of understanding the physiological and metabolic behaviours of microorganisms under hypoxic or microoxic conditions.

Environmental Sciences & Ecology↗
DOE OSTI2025

Propulsion of synthetic protocells and coacervates driven by biochemical catalysis

An important area in the design of biomolecular materials is the mimicry of cellular function through the assembly of molecular components into synthetic cells, or protocells. Through compartmentalization and subsequent evolution, biological cells have achieved the goals of specificity, amplification, and response. By assembling cell-like structures that mimic the organization and size of a natural cell but that incorporate novel components and respond to synthetic triggers, we can extend their capabilities beyond what is currently achievable in biology. For applications in energy, these include the ability to convert chemical to mechanical energy, and to enable the response of sensors and responders once a threshold has been reached. This technical report summarizes progress made on this project as of Dec 31, 2023.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deep mutational scanning reveals the molecular determinants of RNA polymerase-mediated adaptation and tradeoffs

RNA polymerase (RNAP) is emblematic of complex biological systems that control multiple traits involving trade-offs such as growth versus maintenance. Laboratory evolution has revealed that mutations in RNAP subunits, including RpoB, are frequently selected. However, we lack a systems view of how mutations alter the RNAP molecular functions to promote adaptation. We, therefore, measured the fitness of thousands of mutations within a region of rpoB under multiple conditions and genetic backgrounds, to find that adaptive mutations cluster in two modules. Mutations in one module favor growth over maintenance through a partial loss of an interaction associated with faster elongation. Mutations in the other favor maintenance over growth through a destabilized RNAP-DNA complex. The two molecular handles capture the versatile RNAP-mediated adaptations. Combining both interaction losses simultaneously improved maintenance and growth, challenging the idea that growth-maintenance tradeoff resorts only from limited resources, and revealing how compensatory evolution operates within RNAP.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Electronic coherences in nonadiabatic molecular photophysics revealed by time-resolved photoelectron spectroscopy

Significance Time-resolved photoelectron spectroscopy (TRPES) is a promising technique for the study of ultrafast molecular processes, such as the nonadiabatic dynamics taking place at conical intersections. Directly accessing the evolution of the coherences generated at the conical intersection should provide most valuable dynamical information. However, the signals are dominated by background contributions due to state populations, and most theoretical treatments completely neglect the role of the coherences. Here we show that distinguishable signatures of molecular coherences appear in TRPES. These can be recorded using currently available ultrashort pulses and unambiguously extracted at the postprocessing stage. The technique thus provides direct access to nonadiabatic coherence dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Light-driven hydrogen evolution via a novel pincer/no pincer mechanism including a possible concerted proton electron transfer

The road to efficient molecular catalysts for artificial photosynthesis requires a great deal of basic research to address the challenges of anthropogenic global warming. An unexpected result of such research is a novel pincer/no pincer mechanism discovered for the DFT-calculated catalytic cycle of PNP-C2 complexes. Structural parameters from theoretical modeling are supported by single crystal X-ray data for the corresponding complexes. Using the examples of [M(II)Cl(PNP-C2-R)] + (M = Pd 2+ , Ni 2+ ; R = CH 3 , H; PNP-C2 = N,N-bis{(di(2-methoxy-phenyl)phosphanyl)ethyl}-N-alkylamine) this pincer/no pincer mechanism is supported by differential pulse voltammetric, cyclic voltammetric, UV-Vis spectroscopic (combined with voltammetry) and 31 P nuclear magnetic resonance spectroscopic studies. The efficiency of hydrogen evolution during artificial photosynthesis depends significantly on the ratio of photosensitiser and water reduction catalysts. At first sight, the results from irradiation experiments (gas chromatography and mass spectrometry) seem obvious, but show a considerable impact of secondary and tertiary amines within the pincer structure. The best result obtained for the palladium pincer complex was with a secondary amine as a proton relay, delivering a turnover number of 2237 after only 2.5 hours. In contrast the analogous water reduction catalyst with a methylated (tertiary) amine instead of a secondary one gives only 972 turnovers under the same experimental conditions, over the same time. For all photocatalytic measurements the amount of hydrogen is approximately 2.5-times higher when a secondary amine is used as the proton relay. The influences of these different amines on the activity of water reduction catalysts during irradiation experiments in H 2 O or D 2 O lead to evidence of a kinetic isotope effect, supporting a possible concerted proton electron transfer. DFT calculations of reaction enthalpies, activation energies and intrinsic reaction coordinates of the catalytic cycle support the high efficiency of the presented catalysts. In particular, results from secondary amines suggest new candidates as examples of the rarely known concerted proton electron transfer.

14 SOLAR ENERGY↗