Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular driving forces”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Replacing Pyridine with Pyrazine in Molecular Cobalt Catalysts: Effects on Electrochemical Properties and Aqueous H2 Generation

Four new molecular Co(II)tetrapyridyl complexes were synthesized and evaluated for their activity as catalysts for proton reduction in aqueous environments. The pyridine groups around the macrocycle were substituted for either one or two pyrazine groups. Single crystal X-ray analysis shows that the pyrazine groups have minimal impact on the Co(II)–N bond lengths and molecular geometry in general. X-band EPR spectroscopy confirms the Co(II) oxidation state and the electronic environment of the Co(II) center are only very slightly perturbed by the substitution of pyrazine groups around the macrocycle. The substitution of pyrazine groups has a substantial impact on the observed metal- and ligand-centered reduction potentials as well as the overall H2 catalytic activity in a multimolecular system using the [Ru(2,2′-bipyridine)3]Cl2 photosensitizer and ascorbic acid as a sacrificial electron donor. The results reveal interesting trends between the H2 catalytic activity for each catalyst and the driving force for electron transfer between either the reduced photosensitizer to catalyst step or the catalyst to proton reduction step. The work presented here showcases how even the difference of a single atom in a molecular catalyst can have an important impact on activity and suggests a pathway to optimize the photocatalytic activity and stability of molecular systems.

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

Superhydrophilicity of $α$-alumina surfaces results from tight binding of interfacial waters to specific aluminols

Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. Here, contact angles at the alumina (0001) and ($11\bar{2}0$) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

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

Trends in Siting of Metals in Heterometallic Nd–Yb Metal–Organic Frameworks and Molecular Crystals

Several studies suggest that metal ordering within metal–organic frameworks (MOFs) is important for understanding how MOFs behave in relevant applications; however, these siting trends can be difficult to determine experimentally. To garner insight into the energetic driving forces that may lead to nonrandom ordering within heterometallic MOFs, we employ density functional theory (DFT) calculations on several bimetallic metal–organic crystals composed of Nd and Yb metal atoms. We also investigate the metal siting trends for a newly synthesized MOF. Our DFT-based energy of mixing results suggest that Nd will likely occupy sites with greater access to electronegative atoms and that local homometallic domains within a mixed-metal Nd–Yb system are favored. We also explore the use of less computationally extensive methods such as classical force fields and cluster expansion models to understand their feasibility for large system sizes. Here, this study highlights the impact of metal ordering on the energetic stability of heterometallic MOFs and crystal structures.

energy of mixing↗
DOE OSTI2024

Exploration of Tertiary Structure in Sequence-Defined Polymers Using Molecular Dynamics Simulations

Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. Here, we present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.

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

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

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

High‐Performance Fluorinated Fused‐Ring Electron Acceptor with 3D Stacking and Exciton/Charge Transport

Abstract A new fluorinated electron acceptor (FINIC) based on 6,6,12,12‐tetrakis(3‐fluoro‐4‐hexylphenyl)‐indacenobis(dithieno[3,2‐ b ;2 ′ ,3 ′ ‐ d ]thiophene) as the electron‐donating central core and 5,6‐difluoro‐3‐(1,1‐dicyanomethylene)‐1‐indanone as the electron‐deficient end groups is rationally designed and synthesized. FINIC shows similar absorption profile in dilute solution to the nonfluorinated analogue INIC. However, compared with INIC, FINIC film shows red‐shifted absorption, down‐shifted frontier molecular orbital energy levels, enhanced crystallinity, and more ordered molecular packing. Single‐crystal structure data show that FINIC molecules pack into closer 3D “network” motif through H‐bonding and π–π interaction, while INIC molecules pack into incompact “honeycomb” motif through only π–π stacking. Theoretical calculations reveal that FINIC has stronger electronic coupling and more molecular interactions than INIC. FINIC has higher electron mobilities in both horizontal and vertical directions than INIC. Moreover, FINIC and INIC support efficient 3D exciton transport. PBD‐SF/FINIC blend has a larger driving force for exciton splitting, more efficient charge transfer and photoinduced charge generation. Finally, the organic solar cells based on PBD‐SF/FINIC blend yield power conversion efficiency of 14.0%, far exceeding that of the PBD‐SF/INIC‐based devices (5.1%).

Dai, Shuixing↗
DOE OSTI2026

Mo 4 FeGa 17.25– x Ge x : Complementary Point Substitutions, Buffering Frameworks, and Merging of the 18-n and Octet Bonding Schemes

We present the discovery of Mo 4 FeGa 17.25–x Ge x (x ∼ 2.1(4), based on determination of Ge content), a complex gallide featuring a variety of point substitution phenomena. Its crystal structure is derived from the Ti 2 Ni type, in which a diamond network of face-sharing octahedra is interpenetrated by a second diamond network of vertex-sharing stella quadrangula. However, while in the ideal Ti 2 Ni type these two frameworks run uninterrupted through the crystal, Mo 4 FeGa 17.25–x Ge x shows three variations. First, the inner tetrahedron of every other stella quadrangula is replaced with a main group atom, creating tetrahedra reminiscent of the Zintl phase NaTl. Next, a selection of Ga 6 octahedra are filled with Fe atoms, in a manner analogous to the stuffed AuCu 3 -type phases. Finally, the refined crystal structure shows that in each unit cell a single Ga2 atom in the octahedral network is substituted with a dumbbell of Ga/Ge atoms. Electronic structure calculations on ordered models of Mo 4 FeGa 17.25–x Ge x reveal a narrow band near the Fermi energy, which is explained in terms of the 18-n and octet bonding schemes using reversed approximation Molecular Orbital analysis. A DFT-chemical pressure analysis connects the tetrahedron/atom and atom/dumbbell substitutions to the relief of atomic packing tensions and highlights soft atomic motions within the octahedral framework and driving forces for atom/dumbbell substitution. The combination of soft vibrational modes, disorder, and a narrow band gap could make this phase of interest for potential thermoelectric properties.

Chemical structure↗
DOE OSTI2023

Energetics of high temperature degradation of fentanyl into primary and secondary products

Fentanyl is a synthetic opioid used for managing chronic pain. Due to its higher potency (50–100×) than morphine, fentanyl is also an abused drug. A sensor that could detect illicit fentanyl by identifying its thermally degraded fragments would be helpful to law enforcement. While experimental studies have probed the thermal degradation of fentanyl, little theoretical work has been done to understand the mechanism. Here, we studied the thermal degradation pathways of fentanyl using extensive ab initio molecular dynamics simulations combined with enhanced sampling via multiple-walker metadynamics. We calculated the free energy profile for each bond suggested earlier as a potential degradation point to map the thermodynamic driving forces. In conclusion, we also estimated the forward attempt rate of each bond degradation reaction to gain information about degradation kinetics.

Poudel, Bharat↗
DOE OSTI2021

Specific Ion Solvation and Pairing Effects in Glycerol Carbonate

Identifying the driving forces behind the solvation of inorganic salts by nonaqueous solvents is an important step in the development of green solvents. Here we focus on one promising solvent: glycerol carbonate (GC). Using ab initio molecular dynamics simulations, we build upon our previous work by detailing glycerol carbonate’s interactions with a series of anions, a lithium ion, and the LiF ion pair. Through these investigations, we highlight the changes in solvation behavior as the anion size increases, the competition of binding shown by lithium for the oxygens of GC, and the behavior of the LiF ion pair in a GC solution. These results in this work indicate the importance of the cation’s identity in ion-pairing structure and dynamics and lend insight into the key factors behind the specific ion effects seen in GC.

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

Crystals of Benzamide, the First Polymorphous Molecular Compound, Are Helicoidal

Abstract The growth of spontaneously twisted crystals is a common but poorly understood phenomenon. An analysis of the formation of twisted crystals of a metastable benzamide polymorph (form II ) crystallizing from highly supersaturated aqueous and ethanol solutions is given here. Benzamide, the first polymorphic molecular crystal reported (1832), would have been the first helicoidal crystal observed had the original authors undertaken an analysis by light microscopy. Polymorphism and twisting frequently concur as they are both associated with high thermodynamic driving forces for crystallization. Optical and electron microscopies as well as electron and powder X‐ray diffraction reveal a complex lamellar structure of benzamide form II needle‐like crystals. The internal stress produced by the overgrowth of lamellae is shown to be able to create a twist moment that is responsible for the observed non‐classical morphologies.

Shtukenberg, Alexander G.↗
DOE OSTI2021

Atomistic simulations of dislocation mobility in refractory high-entropy alloys and the effect of chemical short-range order

Abstract Refractory high-entropy alloys (RHEAs) are designed for high elevated-temperature strength, with both edge and screw dislocations playing an important role for plastic deformation. However, they can also display a significant energetic driving force for chemical short-range ordering (SRO). Here, we investigate mechanisms underlying the mobilities of screw and edge dislocations in the body-centered cubic MoNbTaW RHEA over a wide temperature range using extensive molecular dynamics simulations based on a highly-accurate machine-learning interatomic potential. Further, we specifically evaluate how these mechanisms are affected by the presence of SRO. The mobility of edge dislocations is found to be enhanced by the presence of SRO, whereas the rate of double-kink nucleation in the motion of screw dislocations is reduced, although this influence of SRO appears to be attenuated at increasing temperature. Independent of the presence of SRO, a cross-slip locking mechanism is observed for the motion of screws, which provides for extra strengthening for refractory high-entropy alloy system.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Gibbs-Thomson Law for Singular Step Segments: Thermodynamics Versus Kinetics

Classical Burton-Cabrera-Frank theory presumes that thermal fluctuations are so fast that at any time density of kinks on a step is comparable with the reciprocal intermolecular distance, so that the step rate is about isotropic within the crystal plane. Such azimuthal isotropy is, however, often not the case: Kink density may be much lower. In particular, it was recently found on the (010) face of orthorhombic lysozyme that interkink distance may exceed 500-600 intermolecular distances. Under such conditions, Gibbs-Thomson law (GTL) may not be applicable: On a straight step segment between two corners, communication between the comers occurs exclusively by kink exchange. Annihilation between kinks of opposite sign generated at the comers results in the grain in step energy entering GTL. If the step segment length l much greater than D/v, where D and v are the kink diffusivity and propagation rate, respectively, the opposite kinks have practically no chance to annihilate and GTL is not applicable. The opposite condition of the GTL applicability, l much less than D/v, is equivalent to the requirement that relative supersaturation Delta(sub mu)/kT much less than alpha/l, where alpha is molecular size. Thus, GTL may be applied to a segment of 10(exp 3)alpha approx. 3 x 10(exp -5)cm approx 0.3 micron only if supersaturation is less than 0.1%, while practically used driving forces for crystallization are much larger. Relationships alternative to the GTL for different, but low, kink density have been discussed. They confirm experimental evidences that the Burton-Cabrera-Frank theory of spiral growth is growth rates twice as low as compared to the observed figures. Also, application of GTL results in unrealistic step energy while suggested kinetic law give reasonable figures.

Chernov, A. A.↗
NASA NTRSDate not supplied

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer↗
NASA NTRSDate not supplied

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer↗
DOE OSTI2021

Determination of Thermal Expansion, Defect Formation Energy, and Defect-Induced Strain of α-U Via ab Initio Molecular Dynamics

Uranium (U) is often alloyed with molybdenum (Mo) or zirconium (Zr) in order to stabilize its high-temperature body-centered cubic phase for use in nuclear reactors. However, in all metallic fuel forms, the α phase of U remains in some fraction. This phase decomposition due to temperature or compositional variance can play an outsized role on fuel performance and microstructural evolution. Relatively little is known about fundamental point defect properties in α -U at non-zero temperatures, from either computational or experimental studies. This work performs the first thorough evaluation of the α phase of U via ab initio molecular dynamics (AIMD). A number of thermophysical properties are calculated as a function of temperature, including equilibrium lattice parameters, thermal expansion, and heat capacity. These results indicate a two-region behavior, with the transition at 400 K. The thermal expansion/contraction in the a / b direction occurs rapidly from 100 up to 400 K, after which a more linear and gradual change in the lattice constant takes place. The volumetric expansion matches experiments quantitatively, but the individual lattice constant expansion only matches experiments qualitatively. Point defect formation energies and induced lattice strains are also determined as a function of temperature, providing insight on defect populations and the fundamentals of irradiation growth in α -U. Interstitials induce significantly more strain than vacancies, and the nature of that strain is highly dependent on the individual lattice directions. The direction of point defect-induced lattice strain is contrary to the irradiation growth behavior of α -U. This work shows that isolated point defects cannot be the primary driving force responsible for the significant irradiation-induced growth of α -U observed experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Finding predictive models for singlet fission by machine learning

Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.

36 MATERIALS SCIENCE↗