Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initio molecular dynamics”

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 325 records · Page 18

Nitrogen: A promising doping strategy for high-performance ovonic threshold switching selectors

The Ovonic Threshold Switching (OTS) selector serves as an essential component in the development of three-dimensional high-density memory integration technology. Nevertheless, the state-of-the-art high-performance OTS materials usually contain toxic elements such as arsenic (As), posing significant risks to both environmental and human health. Nitrogen (N), which belongs to the same group as arsenic (As), has emerged as a highly promising alternative for As doping. However, the underlying mechanisms that govern N-based OTS materials have not yet been extensively investigated. In this study, we delve into the effects of N doping on the structural, bonding, and electronic properties of amorphous GeSe (a-GeNSe) by ab initio molecular dynamics simulations to bridge the knowledge gap. Our findings indicate that upon N doping in a-GeSe, the formation of robust Ge-N bonds, along with N-centered tetrahedral and triangular structures, resulting in the sluggish atomic movement that enhances the thermal stability and endurance of a-GeNSe. The OTS characteristics are significantly influenced by the material’s electronic band structure, and thus the relatively slow performance drift can be attributed to the stabilization of mid-gap states, a result of N doping which effectively slows down the aging process of chalcogenide glass. Moreover, the increased mobility gap in a-GeNSe raises the threshold voltage (V th ), making it more compatible with commercially available phase-change memory materials. Furthermore, our findings reveal the extensive impact of the N element on a typical OTS material and offer valuable perspectives for alternative doping strategies that could potentially supplant As practices.

36 MATERIALS SCIENCE↗

Adsorption of hydroxamic acid ligands for improved extraction of rare earth elements from monazite ores

Efficient separation of rare earth element (REE) ores via froth flotation requires the development of novel ligands with enhanced adsorption capacity and selectivity. To realize these advances, understanding the mechanisms underlying interactions between the ligand and mineral surfaces is essential. This study systematically evaluates the adsorption behavior of alkyl- and aromatic alkyl-substituted hydroxamic acid ligands on monazite surfaces using complementary spectroscopic techniques, including UV–visible (UV–vis) spectroscopy, Raman spectroscopy, infrared spectroscopy, and vibrational sum frequency generation (SFG) spectroscopy, together with the ab initio molecular dynamics (AIMD) simulations. Among the studied ligands, octanohydroxamic acid (OHA) and 4-ethoxy-N,2-dihydroxybenzamide (EDHBA) exhibit high adsorption capacity under basic pH (8–10) by forming multilayers on the surface. OHA has a higher equilibrium adsorption capacity compared to EDHBA, but it forms a less stable multilayer susceptible to disruption in the presence of interfering ions. AIMD results show that OHA adopts a single stable chelating geometry, while EDHBA exhibits multiple binding modes involving distinct interactions with La surface atoms and phosphate-bound oxygens, resulting in more complex adsorption kinetics. The variations in surface binding and intermolecular interactions observed between alkyl and aromatic molecules influence the differences in adsorption kinetics, equilibrium adsorption capacities on the mineral surface, and their flotation performance. This work provides valuable insight into the adsorption mechanism of ligands at mineral interfaces, which is crucial for guiding the design of new ligands with enhanced separation performance.

Zhou, Muchu [ORNL] (ORCID:0000000182650215)↗

Structural and electronic characteristics of amorphous Ge 8 Sb 2 Te 11

GeTe-rich phase-change materials can be utilized in rewriteable optical memory due to the large contrast in reflectivity between amorphous and crystalline phases. Here we explored the structure and electronic properties of amorphous Ge 8 Sb 2 Te 11 using ab initio molecular dynamics simulations. The results indicate that amorphous Ge 8 Sb 2 Te 11 is dominantly composed of Ge-, Sb- and Te-centered octahedrons with distortions, while 30.4 % of Ge-centered clusters are in the form of tetrahedrons which are randomly distributed. The 5-fold rings possess a large proportion, and the Gesingle bondTe and Sbsingle bondTe bonds present larger formation energies than other bonds, leading to the ABAB bonding sequence (A: Ge and Sb, B: Te). The lone pair electrons locating at the opposite direction of bonds possess a large fraction of 14.8 %, which may enhance the distortions of local clusters. Importantly, these structural properties lead to the high stability of amorphous Ge 8 Sb 2 Te 11 and thus long data retention in the information storage.

36 MATERIALS SCIENCE↗

First-principles investigation of the thermophysical properties of NaCl, PuCl 3 , and NaCl-PuCl 3 Molten salts

Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl 3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl 3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl 3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. Here, the heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl 3 .

36 MATERIALS SCIENCE↗

Enhanced ionic conductivity and lack of paddle-wheel effect in pseudohalogen-substituted Li argyrodites

Superionic conductors are key to the development of safe and high-energy-density, all-solid-state batteries. Using a combined theoretical and experimental approach, we explore the feasibility of increasing the ionic conductivity through pseudohalogen substitution in the Li argyrodite structure. Under the guidance of calculated thermodynamic stability, BH 4 -substituted Li argyrodite, Li 5.91 PS 4.91 (BH 4 ) 1.09 , was successfully synthesized via a mechanochemical method. As-synthesized BH 4- substituted Li argyrodite displays an ionic conductivity of 4.8 mS/cm at 25°C. Ab initio molecular dynamics simulation trajectory analysis was used to investigate how BH 4 facilitates Li-ion diffusion and indicates only a weak correlation with the B–H bond motion. Here, we find that the enhanced conductivity mainly originates from the weak interaction between Li and BH 4 and find no evidence of a paddle-wheel effect from the polyanion. This work provides insight into how cluster ions enhance Li diffusion and systematically describes how to explore superionic conductors with pseudohalogen substitution.

25 ENERGY STORAGE↗

Fluorinated co-solvent promises Li-S batteries under lean-electrolyte conditions

Practical lithium-sulfur batteries stipulate the use of a lean-electrolyte and a high Coulombic efficiency of the lithium metal anode. In this work, we employ 1,1,2,2-tetrafluoroethyl 2,2,2-trifluoroethyl ether as a co-solvent in the electrolyte of Li-S batteries to meet the demands. The co-solvent fosters the formation of a LiF-rich solid electrolyte interphase on lithium metal anode, as revealed by Ab initio molecular dynamics and Femtosecond Stimulated Raman spectroscopy. The co-solvent results in an average Coulombic efficiency of 99.4% for lithium plating/stripping cycles. Full cells with a capacity ratio of 2 between lithium anode and the S@PAN cathode (3 gs cm –2 ) exhibited a stable cycle life over 100 cycles at an electrolyte/sulfur ratio of 2 µL mg –1 , validating the high Coulombic efficiency of the lithium anode and demonstrating the compatibility of the electrolyte with both electrodes. To enhance the energy density, we prepared a hybrid cathode composed of 45 wt.% VS 2 mixed with ZnS-coated Li 2 S@graphene as the cathode. Based on the mass of both electrodes and the electrolyte, the full cell delivers an energy density of 483 Wh kg –1 , which demonstrates viable Li-S batteries with the lean-electrolyte.

25 ENERGY STORAGE↗

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis and thermophysical property determination of NaCl-PuCl3 salts

Currently, a knowledge gap exists in the available data and understanding of thermophysical properties relating to fresh fuel salts, especially those containing plutonium. These data are necessary for designing, constructing, and licensing future commercial molten-salt reactors. Thermophysical properties, such as melting temperature, salt stability, density, and heat capacity were ascertained using NaCl-PuCl3 (36 mol% PuCl3) and a more sodium rich composition containing 25 mol% PuCl3. The NaCl-PuCl3 salt mixture was synthesized for this study and contained 63.4 mol% NaCl, 36.3 mol% PuCl3 and was 99.7% pure. Upon heating, the NaCl–PuCl3 eutectic was stable at temperatures up to 800 °C. The onset of melting occurred at 451 ± 3 °C, and the enthalpy of fusion was determined to be 23.0 ± 1.4 kJ/mol. Heat capacity measurements in the liquid phase ranged from 107.7 to 91.3 J/mol.K, with an average value of 104.6 ± 11.4 J/mol.K between 500 and 650 °C. Three independent trials of the molten NaCl-PuCl3 salt found the density to be ?(T) = 3.8589 – 9.5342·10-4 T(°C). In addition, ab initio molecular dynamic simulations to calculate density and heat capacity values are included.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Temperature as a control knob on spin-orbit coupling

Spin-orbit coupling (SOC) governs many physical phenomena. Through ab initio molecular dynamics simulations, Lu and Sun demonstrated that structural disorder at elevated temperatures substantially reduces the effective SOC contribution that stabilizes band inversion, driving a topological-to-normal-insulator transition in Bi2Se3. Their work identifies temperature as a meaningful control parameter for SOC-mediated topology and other properties.

Liang, Liangbo [ORNL] (ORCID:0000000311990049)↗

The stability of FeH x and hydrogen transport at Earth’s core mantle boundary

Iron hydride in Earth’s interior can be formed by the reaction between hydrous minerals (water) and iron. Studying iron hydride improves our understanding of hydrogen transportation in Earth’s interior. Our high-pressure experiments found that face-centered cubic (fcc) FeH x (x≤1) is stable up to 165 GPa, and our ab initio molecular dynamics simulations predicted that fcc FeH x transforms to a superionic state under lower mantle conditions. In the superionic state, H-ions in fcc FeH become highly diffusive-like fluids with a high diffusion coefficient of ~3.7 × 10 -4 cm 2 s -1 , which is comparable to that in the liquid Fe-H phase. The densities and melting temperatures of fcc FeH x were systematically calculated. Similar to superionic ice, the extra entropy of diffusive H-ions increases the melting temperature of fcc FeH. Further, the wide stability field of fcc FeH enables hydrogen transport into the outer core to create a potential hydrogen reservoir in Earth’s interior, leaving oxygen-rich patches (ORP) above the core mantle boundary (CMB).

58 GEOSCIENCES↗

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molten-salt-mediated carbon dioxide capture and superequilibrium utilization with ethane oxidative dehydrogenation

Existing CO 2 -mediated oxidative dehydrogenation (CO 2 -ODH) of ethane has yet to demonstrate >60% single-pass CO yield due to the intrinsic equilibrium limitations. We report a unique approach with mixed molten carbonates as a reaction medium for CO 2 -ODH, which strategically partitions the CO 2 -ODH reactions into gas and molten-salt phases and facilitates integrated CO 2 capture from power plant flue gases. An 89% CO yield was achieved at 770°C, doubling the equilibrium limitation of conventional CO 2 -ODH. The high CO yield in turn enhances ethylene formation. Further characterizations confirmed that molten-salt mediated ODH (MM-ODH) proceeds through a gas-phase cracking and molten-salt mediated reverse water-gas-shift reaction pathway. Based on this understanding, thermodynamic analysis and ab initio molecular dynamics simulations were conducted to develop general principles to optimize the molten-salt reaction medium. Process analyses confirm that MM-ODH has the potential to be significantly more efficient for CO 2 capture and utilization than conventional CO 2 -ODH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An ab initio free energy study of the reaction mechanism and rate-limiting steps of CO 2 capture by aqueous glycine

Aqueous amino acids are promising absorbents for direct air capture (DAC) of CO 2 . Herein, we investigate the possibility of kinetic control of CO 2 absorption with aqueous anionic glycine (GLY - ) by employing extensive ab initio molecular dynamics simulations, free energy analysis, and reaction rate theory. We find that first GLY - binds to CO 2 by overcoming a barrier (7.4 kcal/mol) to form a zwitterion intermediate, which then releases a proton by overcoming a similar barrier. Despite the similarity in the barrier, zwitterion formation appears to be the rate-limiting step because it is two orders of magnitude slower (microseconds) than the proton release step. This is predominantly due to stronger nonequilibrium solvent effects for the former that cause many barrier-recrossing events and effectively slow down the reaction rate. Such a detailed fundamental understanding of the amino acid-based CO 2 -absorption mechanism and rates is key to improving the kinetic efficiency of DAC technology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of proton transport in aqueous acid solutions

The fundamental understanding of proton transport, specifically the Grotthuss-like mechanism, is critical for many technological applications. The present study presents a comprehensive analysis of proton transport in aqueous solutions of sulfuric and phosphoric acids using quasielastic neutron and light scattering, broadband dielectric spectroscopy, rheology, pulsed-field gradient-NMR measurements, and ab initio molecular dynamic simulations. Results show that in all systems, proton transport occurs through short “jumps” of ~0.5 Å that are faster than structural relaxation. However, proton hopping appears to be coupled to structural relaxation in aqueous solutions of sulfuric acid, while these processes are decoupled in phosphoric acid. Neutron scattering indicates that all protons have the same fast mobility in phosphoric acid systems, while at least one proton per sulfuric molecule remains slower than other protons. Analysis reveals correlated proton jumps, but these correlations suppress conductivity, suggesting that expected Grotthuss-like enhancement of conductivity is unlikely in bulk liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of Cation Coordination in the Electrical and Optical Properties of Amorphous Transparent Conducting Oxides

Amorphous oxide semiconductor materials have demonstrated numerous advantages without compromise of electrical properties as compared to their crystalline counterparts, yet understanding of the fundamental principles allowing this has remained elusive. To study the origins of enhanced optoelectronic properties, we apply high-throughput, combinatorial sputtering, structural and spectral mapping, and computationally intensive ab initio molecular dynamics simulations with density functional theory to a ternary, post-transition metal oxide system, namely, zinc tin oxide. The deposited thin films exhibit a high figure of merit, achieving carrier densities in the range of 1019 to 1020 cm–3 and carrier mobilities up to 35 cm2/Vs. These results highlight the role of local distortions and cation coordination in determining the microscopic origins of carrier generation and transport. In particular, we identify the strong likelihood of Sn undercoordination in both Zn-poor and Zn-rich phases leading to the high carrier concentrations observed. This not only diverges from the still widespread historical indictment of oxygen vacancies controlling carrier population in crystalline oxides but also provides a comprehensive framework to describe the unique structure–property relationships using specific structural and electronic descriptors in disordered phase materials.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

Engineering Electrical Conductivity in Stable Zirconium-Based PCN-222 MOFs with Permanent Mesoporosity

Electrically conductive metal–organic frameworks (MOFs) featuring zirconium-based nodes are of great interest for electrochemical and optoelectronic applications due to their exceptional thermal and chemical stability, although the number of such MOFs remains limited. Here in this paper, we demonstrate that electron deficient molecules such as transition metal(IV) bis(dicarbollide) (TM IV (C 2 B 9 H 11 ) 2 ; TM IV CB for short, where TM = Ni, Pd, and Pt) or C 60 , spatially infiltrated in the microporous channels of the zirconium-based porphyrinic MOFs PCN-222 and PCN-222-Zn, generate highly stable and electrically conductive frameworks in which the mesopores remain accessible to other guests. Solid-state density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations indicate incorporation of Ni IV CB and C 60 in the MOF micropores is energetically and structurally feasible. Interestingly, in contrast to the free species, strong donor–acceptor interactions between the MOF and Ni IV CB restrain it to a single conformation. Calculated electronic structures and charge-hopping conduction probabilities illustrate that efficient charge-transfer (CT) from photoexcited linkers to the guest molecules facilitates charge hopping in the framework, making the MOFs electrically conductive. The donor–acceptor conjugates also enhance exciton dissociation at their heterojunctions, fostering the formation of long-lived electron-trapped states with potential utility for photo- and electrochemical devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anionic Oxygen Redox in the High-Lithium Material Li 8 SnO 6

Lithiated transition metal oxide cathode materials that combine oxygen-anion redox [O 2 – /(O 2 ) n– ] with cation redox offer substantial capacities and higher voltage than cathodes without oxygen redox, but oxygen-anion redox is usually accompanied by the formation of peroxo or superoxo species that may cause oxygen release. To clarify the relationship between the oxygen release and peroxide and superoxide dimers without the complication of transition-metal redox, we performed density functional theory calculations to study the layered Li 8 SnO 6 cathode because it has only oxygen-anion redox during delithiation. Features with O–O distances of ~1.5 Å (corresponding to the formation of a peroxide) were observed at 75% lithium concentration, and further delithiation to 62.5% lithium concentration shortened the O–O bond distance to ~1.3 Å (corresponding to the formation of superoxide). The solely anionic oxygen redox stabilizes the delithiated structure of Li 8–x SnO 6 , and it does not cause structural disorder or release oxygen gas. Ab initio molecular dynamics calculations in the supercell at 300 K indicate that oxygen dimerization may occur at 75% lithium concentration. On the (0001) surface with a high lithium concentration (98%), oxygen dimerization is unfavorable. The calculations by HSE06 show that oxygen redox in Li 8 SnO 6 offers an initial voltage (vs Li/Li + ) greater than 4.2 V. The GGA+U calculations indicate that the Li + -polaron in Li 95 Sn 12 O 72 migrates with a migration barrier of only 0.43 eV. Here, this work also quantifies a solely anionic oxygen redox mechanism in Li 8 SnO 6 , and it provides a deeper understanding of anionic oxygen redox in Li-excess oxide cathode materials.

25 ENERGY STORAGE↗