Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio thermodynamics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗
DOE OSTI2026

Surface Reordering During Layer‐by‐Layer Growth on SrTiO 3

With the development of layer‐by‐layer growth techniques such as molecular beam epitaxy (MBE), it is now possible to construct materials with atomic‐level precision, with the layer sequence primarily determined by the growth recipe. However, materials can still restructure in the high‐temperature growth environment with surface thermodynamics strongly influencing growth behavior. Here, we demonstrate that growth on ${\rm SrTiO}_{3}$ (001), the principal platform for oxide electronics, does not occur in a layer‐by‐layer fashion but follows a more complex process in which a ${\rm TiO}_{2}$ plane continually diffuses toward the growth surface. Employing in situ synchrotron X‐ray scattering combined with ab initio thermodynamic calculations, we discover the existence of a stable ${\rm TiO}_{2}$ double‐layer structure on the pristine substrate and the occurrence of dynamic layer rearrangement during homoepitaxial growth by oxide MBE. Our findings suggest that the current methods used to precisely control ${\rm SrTiO}_{3}$ surfaces are limited, as well as our ability to dictate the composition of ultrathin films, resulting in important ramifications regarding the surface reactivity of perovskite materials grown on ${\rm SrTiO}_{3}$.

Coherent Bragg Rod Analysis↗
DOE OSTI2020

Spherically expanding flame in silane–hydrogen–nitrous oxide–argon mixtures

The effect of silane addition on the laminar flame speed ($S^0_u$) of flames propagating in hydrogen–nitrous oxide–argon mixtures has been investigated experimentally for the first time using the spherically expanding flame technique in a constant volume combustion chamber. Replacing hydrogen by silane and maintaining the equivalence ratio constant, much higher flame speeds, explosion peak pressures, and pressure rise coefficients were measured. A previously developed detailed reaction model has been updated based on ab initio thermodynamic properties calculations and collision limit violation analysis. The improved reaction model demonstrates encouraging performance in predicting the flame speed, with a mean absolute error below 11%. To explain the effect of silane addition on the flame dynamics, a number of parameters have been calculated including OH and H rate of production, heat release rate per reaction, and sensitivity coefficient on . The dynamics of freely propagating flames in SiH 4 –H 2 –N 2 O–Ar mixtures is essentially controlled by reactions of the H–O–N chemical system: N 2 O+H=N 2 +OH, OH+H 2 =H 2 O+H, and N 2 O(+M)=N 2 +O(+M). Whereas silane addition does not influence much the rate of production of OH, it significantly modifies that of H with a number of pyrolytic chemical pathways of silicon hydrides, such as SiH+H 2 =SiH 2 +H and Si+H 2 =SiH+H, which act as sink of H atom as they proceed in the backward direction. The reactions forming SiO(s) and SiO 2 (s), such as SiO+OH=SiO 2 (s)+H and 2SiO=2SiO(s), are exothermic and significantly contribute to the temperature increase. The adiabatic, constant pressure flame temperature for mixture containing silane is significantly higher, up to several 100’s K. The increase of induced by silane addition seems to be mostly related to the large increase of the flame temperature which leads to higher energy release rate.

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

Computational descriptor for electrochemical currents of carbon dioxide reduction on Cu facets

Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO 2 R) reduction that produce valuable hydrocarbons and oxygenates. In this work, leveraging density functional theory calculations for the CO adsorption energy ΔE CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE CO and the experimentally measured CO 2 R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining ab initio thermodynamics of early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically controlling step. Beyond the general CO adsorption energy that only shows a linear trend with CO 2 R activity, we show that the reaction free energy of CO* → CHO* is the descriptor for the overall CO 2 R activity for Cu facets, as it displays a volcano relationship with the experimental current. Importantly, we show that high step and kink density of the Cu terminations not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. In addition, we explain that the high activity of (741) is due to its relatively low hydrogen evolution reaction activity compared with the other Cu surfaces.

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

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2025

First-Principles Study of Molecular Hydrogen Activation by Defects in Boron Nitride

Here, we used density functional theory simulations in combination with ab initio thermodynamics to determine the H 2 partial pressure (p H2 )-dependent energetics associated with H 2 activation and recovery at various defect sites in hexagonal boron nitride (h-BN). We found that some defects are very reactive with hydrogen, thereby definitely trapping hydrogen in defective h-BN. However, depending on hydrogen partial pressure, less reactive defect sites can be populated. Because of the lower binding capability of these sites, they would allow hydrogen to be recycled and recovered. For small defect sizes, we found that hydrogen preferentially binds to nitrogen sites by forming N–H bonds, and if no N sites are available then boron sites would be the next to bind hydrogen. Hydrogen dissociation via frustrated Lewis pair is found to be more favorable than forming only N–H bonds but only if the defect size is large enough to accommodate steric effects. For specific conditions such as T = 400 K, p H2 = 1 bar, and only considering one molecular H 2 per defect, three defects, namely, the N monovacancy, 3V(1B2N), and hexagonal 6V(3B3N) could play a role in both the activation and recycling of H 2 as they would be reacting enough to allow a favorable splitting of H 2 while not binding too strongly to allow its recovery. More broadly, a range of p H2 and hydrogen loading conditions were investigated for different types of defects and the finding suggests that p H2 could be used to fine-tune the Gibbs free energy of hydrogenation, thereby allowing several types of defects at different hydrogen loading contents to play a role in the activation/recovery process of H 2 in defective h-BN.

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

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗
DOE OSTI2022

Dynamical Study of Adsorbate-Induced Restructuring Kinetics in Bimetallic Catalysts Using the PdAu(111) Model System

The dynamic restructuring of bimetallic catalysts plays a crucial role in their catalytic activity and selectivity. In particular, catalyst pretreatment with species such as carbon monoxide and oxygen has been shown to be an effective strategy for tuning the surface composition and morphology. Mechanistic and kinetic understanding of such restructuring are fundamental to the chemistry and engineering of surface active sites but have remained challenging due to the large structural, chemical, and temporal degrees of freedom. Here, we combine time-resolved temperature-programmed infrared reflection absorption spectroscopy, ab initio thermodynamics, and machine-learning molecular dynamics to uncover previously unidentified timescale and kinetic parameters of in situ restructuring in Pd/Au(111), a highly relevant model system for dilute Pd-in-Au nanoparticle catalysts. The key innovation lies in utilizing CO not only as a chemically sensitive probe of surface Pd, but also as an agent that induces restructuring of the surface. Upon annealing in vacuum, as-deposited Pd islands became encapsulated by Au and partially dissolved into the subsurface, leaving behind isolated Pd monomers on the surface. Subsequent exposure to 0.1 mbar CO enabled Pd monomers to repopulate the surface up to 373 K, above which complete Pd dissolution occurred by 473 K, with apparent activation energies of 0.14 and 0.48 eV, respectively. Furthermore, these restructuring processes occurred over the span of ~1000 s at a given temperature. Such a minute-timescale dynamics not only elucidates the fluxional nature of alloy catalysts but also presents an opportunity to fine-tune the surface at moderate temperature and pressure conditions.

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

Computational and Experimental Characterization of the Ligand Environment of a Ni-Oxo Catalyst Supported in the Metal–Organic Framework NU-1000

Heterogeneous catalysts exhibit significant changes in composition due to the influence of operating conditions, and these compositional changes can have dramatic effects on catalytic performance. For traditional bulk metal heterogeneous catalysts, relationships between composition and catalytic operating conditions are well documented. However, the influence of operating conditions on the compositions of single-site heterogeneous catalysts remains largely unresolved. To address this, we report a combined computational and experimental characterization of a Ni oxo catalyst under catalytic hydrogenation conditions. Specifically, pair distribution function (PDF) analysis is combined with ab initio thermodynamic modeling to investigate ligand environments present on a Ni oxo cluster supported in the metal–organic framework NU-1000. Comparisons of the experimentally observed and simulated Ni–O coordination numbers and Ni–O, Ni···Ni, and Ni···Zr distances provide insight into the Ni ligand environment under H 2 (g). These comparisons suggest significant OH and H 2 O content and, further, that different Ni ions within the cluster and/or NU-1000 structure may comprise subtly different numbers of these ligands. Further, the observation of significant H 2 O content under H 2 (g) suggests that the NU-1000 support supplies H 2 O to the cluster. Examples of ligand environments that could lead to the observed PDFs are provided. Furthermore, the combination of simulations and experiments provides new insights into the ligand environment for Ni-NU-1000 catalysts that will be useful for understanding the ligand environments of other single-site Ni catalysts as well.

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

Hydrogenation of Covalent Organic Framework Induces Conjugated π Bonds and Electronic Topological Transition to Enhance Hydrogen Evolution Catalysis

Recently, many topological materials have been discovered as promising electrocatalysts in chemical conversion processes and energy storage. However, it remains unclear how the topological electronic states specifically modulate the catalytic reaction. Here, the two-dimensional metal phthalocyanine-based covalent organic framework (MP c -COF) is studied by ab initio thermodynamic calculations to clearly reveal the promotional effect on the electrochemical hydrogen evolution reaction (HER) induced by topological gapless bands (TGBs). We find that the prehydrogenated (and fluorinated) H 4 CdPc-COF(F) shows the best HER performance, with 0.016 V (near zero) overpotential. By tracking changes to the electronic structure and free energy as the prehydrogenation and HER processes occur, we are able to separately attribute the high HER efficiency in part due to the increase of the electron bath by donating electrons to the conjugated π bonds and also to the existence of TGBs. Specifically, the significant catalytic promotion by TGBs is proven to decrease the free energy by 0.218 eV to near zero. When the TGBs are destroyed, e.g., by replacing N with P and opening a band gap, the HER efficiency is reduced. As a result, this study opens avenues for deterministically harnessing topological band features to improve electrocatalysis.

25 ENERGY STORAGE↗
DOE OSTI2020

Oxidative dehydrogenation of propane on the oxygen adsorbed edges of boron nitride nanoribbons

Metal-free boron nitride has recently been reported to exhibit high catalytic activity and selectivity towards oxidative dehydrogenation of propane (ODHP). Despite several experimental studies, the exact nature and function of the active sites in this emerging catalyst are unknown. In the present work, density functional theory calculations are combined with microkinetic modeling to systematically explore the ODHP on oxygen passivated boron nitride nanoribbons (BNNRs). Here, the relative stabilities of different edge structures of BNNR were examined using ab initio thermodynamic analysis, and the most stable oxygen passivated edge structure was used for a mechanistic study. Microkinetic analysis revealed that a N 2 O or NO x -type active site is active and selective for the conversion of propane to propene. In addition to a heterogeneous catalytic cycle, the proposed N 2 O/NO x -type active sites are expected to be able to generate gas phase C 3 H 7 ˙, which in turn can trigger gas phase reactions of ODHP as experimentally speculated.

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

Adsorption of ethylenediamine on Cu surfaces: attributes of a successful capping molecule using first-principles calculations

The shape-controlled synthesis of Cu nanocrystals can benefit a wide range of applications, though challenges exist in achieving high and selective yields to a particular shape. Capping agents play a pivotal role in controlling shape, but their exact role remains ambiguous. In this study, the adsorption of ethylenediamine (EDA) on Cu(100) and Cu(111) was investigated with quantum density functional theory (DFT) to reveal the complex roles of EDA in promoting penta-twinned Cu nanowire growth. We find EDA has stronger binding on Cu(100) than on Cu(111), which agrees the general expectation that penta-twinned Cu nanowires express facets with stronger capping-molecule binding. Despite this stronger binding, ab initio thermodynamics reveals the surface energy of EDA-covered Cu(111) is lower than that EDA-covered Cu(100) at all solution-phase EDA chemical potentials, so there is no thermodynamic driving force for penta-twinned nanowires. We also investigated the capability of EDA to protect Cu surfaces from oxidation in water by quantifying energy barriers for a water molecule to diffuse through EDA layers on Cu(100) and Cu(111). The energy barrier on Cu(100) is significantly lower, which supports observations of faster oxidation of Cu(100) in electrochemical experiments. Thus, we elucidate another possible function of a capping agent – to enable selective oxidation of crystal facets. This finding adds to the general understanding of successful attributes of capping agents for shape-selective nanocrystal growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lithium cuprate, a multifunctional material for NO selective catalytic reduction by CO with subsequent carbon oxide capture at moderate temperatures

Li 2 CuO 2 was evaluated as a possible catalyst for the NO selective catalytic reduction (NO SCR) by CO. The initial idea was supported by ab initio thermodynamic analysis. This analysis showed that NO SCR was favored by using Li 2 CuO 2 as a catalyst, based on the calculated heat of reaction and Gibbs free energy. The NO SCR reaction was evaluated under dynamic, isothermal and cyclic conditions using different NO : CO ratios (2 : 1, 1 : 1 and 1 : 2). The dynamic and isothermal experiments showed that Li 2 CuO 2 does reduce NO in the presence of CO to reach higher conversions (>90%) between 400 and 450 °C. However, it must be mentioned that the catalytic process was reduced once some of the produced CO 2 was captured on the ceramic. At that point, Li 2 CuO 2 acts as a bifunctional material, catalyzing the NO SCR reaction and capturing some of the produced CO 2 . Moreover, cyclic experiments evidenced high NO reduction. Finally, a different set of experiments showed that this catalytic reduction performed differently if NO and CO were previously fed. Finally, when lithium cuprate is previously carbonated (initial CO feed), the catalytic reaction does not complete. In contrast, when there was a previous nitrate formation (initial NO feed), the reaction reached 100% efficiency.

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

Elucidating the role of reactive nitrogen intermediates in hetero-cyclization during hydrothermal liquefaction of food waste

Hydrothermal liquefaction (HTL) has emerged as a promising strategy for converting abundant, water-rich organic streams into an energy-dense, sustainable, biocrude. Despite major strides in improving oil yields and process intensification, a key underlying challenge persists in relating the composition of the feedstock to the ultimate fate of nitrogen, which is especially important for biocrude quality. To elucidate how nitrogen partitions into the respective gas, aqueous, oil and solid hydrothermal products, food waste was treated under HTL conditions (15 wt% solids, 575 K reaction temperature, and 1 hour reaction time) and elemental, chromatographic and mass spectroscopy analyses were performed on the products to establish dominant product classes and to close mass, carbon, and nitrogen balances. Liquid products (aqueous and biocrude) were found to contain nearly 60% of the nitrogen, with a majority of the biocrude-phase nitrogen in the form of amides. Ab initio thermodynamic and kinetic simulations on model reactants were used to evaluate potential reaction pathways involving reactive nitrogen intermediates. A subsequent kinetic model was evaluated at reaction temperatures ranging 300–1000 K and for a range of feedstock compositions to identify primary reaction pathways. The Maillard reaction is revealed to be the dominant pathway for converting reactive nitrogen intermediates into observed products, especially those in which nitrogen is present as a heterocyclic aromatic. Subsequent product analysis by Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) confirms strong agreement among the experimental data, kinetic model and proposed Maillard pathways.

09 BIOMASS FUELS↗
DOE OSTI2023

Theoretical assessments of Pd–PdO phase transformation and its impacts on H 2 O 2 synthesis and decomposition pathways

The direct synthesis of H 2 O 2 from O 2 and H 2 provides a green pathway to produce H 2 O 2 , a popular industrial oxidant. Here, in this study, we theoretically investigate the effects of Pd oxidation states, coordination environments, and particle sizes on primary H 2 O 2 selectivities, assessed by calculating the ratio of rate constants for the formation of H 2 O 2 (via OOH* reduction; k O–H ) and the decomposition of OOH* (via O–O cleavage; k O–O ). For Pd metals, the k O–H /k O–O ratio decreased from 10 -4 for Pd(111) to 10 -10 for the Pd 13 cluster at 300 K, indicating poorer H 2 O 2 selectivity as Pd particle size decreases and low primary selectivities for H 2 O 2 overall. As the oxygen chemical potential increases and metals form surface and bulk oxides, the perturbation of Pd–Pd ensemble sites by lattice O atoms results in selectivities that become dramatically higher than unity. For instance, at 300 K, the k O–H /k O–O ratio increases significantly from 10 -4 to 10 9 to 10 16 as Pd(111) oxidizes to Pd 5 O 4 /Pd(111) and to PdO(100), respectively. In contrast, such selectivity enhancements are not observed for surface and bulk oxides that persistently contain rows of more metallic, undercoordinated Pd–Pd ensemble sites, such as PdO(101)/Pd(100) and PdO(101). These Pd–Pd ensembles are also absent when smaller Pd nanoparticles fully oxidize, indicating that smaller PdO clusters can be more selective for H 2 O 2 synthesis. These trends for primary H 2 O 2 selectivities were found to inversely correlate with trends for H 2 O 2 decomposition rates via O–O bond cleavage, demonstrating that catalysts with high primary H 2 O 2 selectivity can also hinder H 2 O 2 decomposition. Ab initio thermodynamic calculations are used to estimate the thermodynamically favored phase among Pd, PdO/Pd and PdO in O 2 , H 2 O 2 /H 2 O, and O 2 /H 2 environments. These results are combined to show that smaller Pd nanoparticles are more prone to be oxidized at lower oxygen chemical potentials, upon which they become more selective than larger Pd particles for H 2 O 2 synthesis.

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

A first principles analysis of potential-dependent structural evolution of active sites in Fe-N-C catalysts

Fe-N-C (iron–nitrogen–carbon) electrocatalysts have emerged as potential alternatives to precious metal-based materials for the oxygen reduction reaction (ORR). However, the structure of these materials under electrochemical conditions is not well understood, and their poor stability in acidic environments poses a formidable challenge for successful adoption in commercial fuel cells. To provide molecular-level insights into these complex phenomena, we combine periodic density functional theory (DFT) calculations, exhaustive treatment of coadsorption effects for ORR reaction intermediates, including O and OH, and comprehensive analysis of solvation stabilization effects to construct voltage-dependent ab initio thermodynamic phase diagrams that describe the in situ structure of the active sites. These structures are further linked to activity and stability descriptors that can be compared with experimental parameters such as the half-wave potential for ORR and the onset potential for carbon corrosion and CO 2 evolution. The results indicate that pyridinic Fe sites at zigzag carbon edges, as well as other edge sites, exhibit high activity for ORR compared to sites in the bulk. However, edges neighboring the active sites are prone to instability via overoxidation and consequent site loss. The results suggest that it could be beneficial to synthesize Fe-N-C catalysts with small sizes and large perimeter edge lengths to enhance ORR activity, while voltage fluctuations should be limited during fuel cell operation to prevent carbon corrosion of overoxidized edges.

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

Surface science of shape-selective metal nanocrystal synthesis from first-principles: Growth of Cu nanowires and nanocubes

The authors present results from first-principles density functional theory aimed at understanding the aqueous solution-phase growth of fivefold twinned copper nanowires and single-crystal nanocubes capped by hexadecylamine (HDA). The role of solution-phase chloride, present in the Cu salt or as an additive, is emphasized. Using ab initio thermodynamics, the authors delineate the range of solution-phase conditions, characterized by the chemical potentials of chloride and HDA, under which Cu nanowires and nanocubes can be grown. The authors discuss the likelihood of thermodynamic and/or kinetic nanostructures for various solution-phase concentrations. Their results are in good agreement with experiments and indicate that methods and insights developed for surface science in gas-phase or vacuum conditions can yield much insight into liquid-phase systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LANL Contributions to pyiron: An integrated development environment for computational materials science

To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron (http://pyiron.org). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

Janssen, Jan↗