Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (DFT) calculations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phase stability and coercivity in La 2 Fe 14 B magnet

Critical rare-earth free La 2 Fe 14 B (2:14:1) has the potential to be a gap permanent magnet. However, La 2 Fe 14 B decomposes into La, α-Fe, and La 4 Fe 4 phases below 1067 K. The phase stability and coercivity have been studied in La 2 Fe 14 B magnet using first principles DFT (density functional theory) calculation and micromagnetic simulation. For a perfect La 2 Fe 14 B cube (edge length of 256 nm) without any structural defects and soft magnetic secondary phases, the coercivity (8.5 kOe) is reduced to ~40% of its magnetocrystalline anisotropy field ( H A = 20 kOe). Further, the coercivity sharply reduces to 3.2 kOe upon forming a thin layer (2 nm) of α-Fe on the surface of the La 2 Fe 14 B cube particle. The DFT calculations indicate that a partial replacement of La by other rare-earth (RE) elements can enhance the structural stability of 2:14:1. The gains in cohesive energy are 0.75, 0.10, and 0.33 eV per formula unit in (La 0.5 RE 0.5 ) 2 Fe 14 B with RE = Ce, Pr, and Nd, respectively. Stabilizing the 2:14:1 structure and mitigating the formation of soft magnetic structural defects or impurity phases such as α-Fe is necessary to develop La 2 Fe 14 B based magnet, which can be moderately achieved via partial substitution of La by other rare earth elements such as Ce, Pr, and Nd.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Theoretical assessments of CO 2 activation and hydrogenation pathways on transition-metal surfaces

Carbon dioxide (CO 2 ) hydrogenation on transition-metal active sites offers a promising carbon utilization route toward mitigating greenhouse gas emissions. C 1 products are often formed in parallel during CO 2 hydrogenation, prompting investigations into the intrinsic properties of transition metals that drive activity and product selectivity. Here, in this work, close-packed surfaces of a selection of transition-metal catalysts (Ni, Co, Rh, Ru, Pd, and Pt) were studied with density functional theory (DFT) calculations to understand their fundamental reactivities for CO 2 transformation reactions. Results indicate that CO 2 conversion proceeds through CO* formation and hydrogenation to form C 1 products (* denotes an adsorbed species). Ni, Co, Rh, and Ru favor CO/CH 4 formation, while Pd and Pt favor CO/CH 3 OH formation. The ability of a metal to dissociate C-O bonds drives selectivity between CH 4 and CH 3 OH, while competition between CO* desorption and surface hydrogenation describes CO selectivities. The C-O bond dissociation steps often impose the highest barrier along CH 4 formation reaction profiles, suggesting their kinetic relevance for CH 4 formation rates. The provided DFT-derived data sets detail a comprehensive reaction network of elementary steps relevant to C 1 chemistries, ultimately offering a benchmark for insights into design strategies for materials that exploit transition-metal active sites in carbon capture or utilization processes.

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

Site-specific surface reactivity on MgO for atomic layer deposition via selective hydration

Atomic layer deposition (ALD) is a powerful technique for thin film synthesis, offering atomic-scale precision and conformality. While ALD of MgO has been widely studied for applications in energy storage and microelectronics, its potential as surface on which deposition may be selective and defects repaired remains underexplored. Here, we present a combined theoretical and experimental investigation of MgO surface hydration and its implications for targeted ALD growth using water and dimethyl aluminum isopropoxide (DMAI) as reactants. We perform density functional theory (DFT) calculations to examine molecular and dissociative H 2 O adsorption on MgO (100) terraces and step-edge sites, including pristine surfaces and those with Mg/O vacancies. Reaction Gibbs free energies are calculated under various conditions to quantify surface reactivity. Our findings reveal facet- and defect-dependent hydration behaviors that align with experimental ALD growth trends on MgO (100). This study provides a molecular-level understanding of MgO surface chemistry critical for optimizing ALD processes for thin film growth and defect repair.

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

Enhancing Photostability of Sn-Pb Perovskite Solar Cells by an Alkylammonium Pseudo-Halogen Additive

High-performance tin-lead perovskite solar cells (PSCs) are needed for all-perovskite-tandem solar cells. However, iodide related fast photodegradation severely limits the operational stability of Sn-Pb perovskites despite the demonstrated high efficiency and thermal stability. Herein, this work employs an alkylammonium pseudo-halogen additive to enhance the power conversion efficiency (PCE) and photostability of methylammonium (MA)-free, Sn-Pb PSCs. Density functional theory (DFT) calculations reveal that the pseudo-halogen tetrafluoroborate (BF 4 - ) has strong binding capacity with metal ions (Sn 2+ /Pb 2+ ) in the Sn-Pb perovskite lattice, which lowers iodine vacancy formation. Upon combining BF 4 - with an octylammonium (OA + ) cation, the PCE of the device with a built-in light-scattering layer is boosted to 23.7%, which represents a new record for Sn-Pb PSCs. The improved efficiency benefits from the suppressed defect density. Under continuous 1 sun illumination, the OABF 4 embodied PSCs show slower generation of interstitial iodides and iodine, which greatly improves the device photostability under open-circuit condition. Moreover, the device based on OABF 4 retains 88% of the initial PCE for 1000 h under the maximum-power-point tracking (MPPT) without cooling.

14 SOLAR ENERGY↗
DOE OSTI2024

Chemical Environment and Structural Variations in High Entropy Oxide Thin Film Probed with Electron Microscopy

For this work, we employ analytical transmission electron microscopy (TEM) to correlate the structural and chemical environment variations within a stacked epitaxial thin film of the high entropy oxide (HEO) Mg 0.2 Co 0.2 Ni 0.2 Cu 0.2 Zn 0.2 O (J14), with two layers grown at different substrate temperatures (500 and 200 °C) using pulsed laser deposition (PLD). Electron diffraction and atomically resolved STEM imaging reveal the difference in out-of-plane lattice parameters in the stacked thin film, which is further quantified on a larger scale using four-dimensional STEM (4D-STEM). In the layer deposited at a lower temperature, electron energy loss spectroscopy (EELS) mapping indicates drastic changes in the oxidation states and bonding environment for Co ions, and energy-dispersive X-ray spectroscopy (EDX) mapping detects more significant cation deficiency. Ab initio density functional theory (DFT) calculations validate that vacancies on the cation sublattice of J14 result in significant electronic and structural changes. The experimental and computational analyses indicate that low temperatures during film growth result in cation deficiency, an altered chemical environment, and reduced lattice parameters while maintaining a single phase. Our results demonstrate that the complex correlation of configurational entropy, kinetics, and thermodynamics can be utilized for accessing a range of metastable configurations in HEO materials without altering cation proportions, enabling further engineering of functional properties of HEO materials.

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

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Interaction mechanism of transition metal phthalocyanines on transition metal nitride supports

Here, we investigated the electronic interactions between transition metal phthalocyanine (TMPc’s) on a refractory transition metal nitride support, specifically copper phthalocyanine (CuPc) on titanium nitride (TiN). X-ray Photoelectron Spectroscopy (XPS) results suggest a presence of a few nanometer native oxide layer on the surface of the TiN nanoparticles, which consists of TiN, TiO 2 , and Titanium oxynitrides (Ti x O y N z ). A TiNCuPc nanocomposite was synthesized via a simple mixing method due to the strong binding between CuPc and TiN confirmed by density functional theory (DFT) calculations. Both XPS data and DFT calculations revealed an electron transfer from TiN substrate to CuPc molecule. The nature of charge transfer is not influenced by the presence of an oxide layer on the surface of TiN. Substantial deviations are however found between photoelectron emission microscopy (PEEM) measured work function for TiN (4.68 eV) and theoretically calculated work function for pristine stoichiometric TiN (2.63 eV). This behavior is attributed to the presence of an oxide layer on the TiN surface. TiNCuPc composite system has a work function value between those of TiN and CuPc. Our studies open up an opportunity to apply a new class of materials based on transition metal phthalocyanine/transition metal nitride composites to catalysis and optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of Fluoroethylene Carbonate Additives on the Initial Formation of the Solid Electrolyte Interphase on an Oxygen-Functionalized Graphitic Anode in Lithium-Ion Batteries

The formation of a solid electrolyte interphase (SEI) at the electrode/electrolyte interface substantially affects the stability and lifetime of lithium-ion batteries (LIBs). One of the methods to improve the lifetime of LIBs is by the inclusion of additive molecules to stabilize the SEI. To understand the effect of additive molecules on the initial stage of SEI formation, in this work we compare the decomposition and oligomerization reactions of a fluoroethylene carbonate (FEC) additive on a range of oxygen-functionalized graphitic anodes to those of an ethylene carbonate (EC) organic electrolyte. A series of density functional theory (DFT) calculations augmented by ab initio molecular dynamics (AIMD) simulations reveal that EC decomposition on an oxygen-functionalized graphitic ($11\bar{2}0$) edge facet through a nucleophilic attack on an ethylene carbon site (C E ) of an EC molecule (S2 mechanism) is spontaneous during the initial charging process of LIBs. However, decomposition of EC through a nucleophilic attack on a carbonyl carbon (C C ) site (S1 mechanism) results in alkoxide species regeneration that is responsible for continual oligomerization along the graphitic surface. In contrast, FEC prefers to decompose through an S1 pathway, which does not promote alkoxide regeneration. Including FEC as an additive is thus able to suppress alkoxide regeneration and results in a smaller and thinner SEI layer that is more flexible toward lithium intercalation during the charging/discharging process. In addition, we find that the presence of different oxygen functional groups at the surface of graphite dictates the oligomerization products and the LiF formation mechanism in the SEI.

25 ENERGY STORAGE↗
DOE OSTI2026

Pb-Apatite Framework as a Generator of Novel Flat-Band CuO-Based Physics

Based upon density functional theory (DFT) calculations, we present the basic electronic structure of CuPb 9 (PO 4 ) 6 O (Cu-doped lead apatite, aka LK-99), in two scenarios: (1) where the structure is constrained to the P3 symmetry and (2) where no symmetry is imposed. At the DFT level, the former is predicted to be metallic while the latter is found to be a charge-transfer insulator. In both cases the filling of these states is nominally d 9 , consistent with the standard Cu 2+ valence state, and Cu with a local magnetic moment of order 0.7 μ B . In the metallic case we find these states to be unusually flat (∼0.2 eV dispersion), giving a very high density of electronic states (DOS) at the Fermi level that we argue can be a host for novel electronic physics. The flatness of the bands is the likely origin of symmetry-lowering gapping possibilities that would remove the spectral weight from E F . Motivated by some initial experimental observations of metallic or semiconducting behavior, we propose that disorder (likely structural) is responsible for closing the gap. Here, we consider a variety of possibilities that could possibly close the charge-transfer gap but limit consideration to kinds of disorder that preserve electron count. Of the possible kinds we considered (spin disorder, O populating vacancy sites, and Cu on less energetically favorable Pb sites), the local Cu moment, and consequently the charge-transfer gap, remains robust. We conclude that disorder responsible for metallic behavior entails some kind of doping where the electron count changes. Further, we claim that the emergence of the flat bands should be due to weak wave function overlap between the orbitals on Cu and O sites, owing to the directional character of the constituent orbitals. Therefore, finding an appropriate host structure for minimizing hybridization between Cu and O while allowing them to still weakly interact should be a promising route for generating flat bands at E F which can lead to interesting electronic phenomena, regardless of whether LK-99 is a superconductor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

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

Elucidation of Active Sites in Aldol Condensation of Acetone over Single-Facet Dominant Anatase TiO 2 (101) and (001) Catalysts

Aldol condensations of carbonyl compounds for C–C bond formation are a very important class of reactions in organic synthesis and upgrading of biomass-derived feedstocks. However, the atomic level understanding of reaction mechanisms and structure-activity correlation on widely used transition metal oxide catalysts are limited due to the high degree of struc-tural heterogeneity of catalysts such as commercial TiO 2 powders. Here, we provide a deep understanding of the reaction mechanisms, kinetics, and structure–function relationships for vapor phase acetone aldol condensation through the con-trolled synthesis of two catalysts with high surface areas and clean, dominant facets, coupled with detailed characterization and kinetic studies that are further assisted by density functional theory (DFT) calculations. Temperature-dependent diffuse reflectance infrared Fourier transform spectroscopy showed the existence of abundant acetone bonded to surface hydroxyl groups (acetone-O s H) and acetone bonded to Lewis acid sites (acetone-Ti 5c ) on surface of both {101} and {001} facet dom-inant TiO 2 . Intermolecular C–C coupling of enolate intermediate from acetone-Ti 5c and a vicinal acetone-O s H is a kinetically relevant step, which is consistent with kinetic and isotopic studies and DFT calculations. The {001} facet showed a lower apparent activation energy (or higher reactivity) than the {101} facet. This is likely caused by the weaker Lewis acid and Brønsted base strengths of the {001} facet which favors the reprotonation–desorption of the coupled intermediate, making the C–C coupling step more exothermic on the {001} facet and resulting in an earlier transition state with a lower activation barrier. It is also possible that the {001} facet has a smoother surface configuration and less steric hindrance during intermo-lecular C–C bond formation than the {101} facet.

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

Hydrodechlorination of 1,2-Dichloroethane on Platinum Catalysts: Insights from Reaction Kinetics Experiments, Density Functional Theory, and Microkinetic Modeling

Catalytic hydrodechlorination is a promising strategy for treating industrial 1,2-dichloroethane wastes, for which Pt and Pt-based alloy catalysts are widely used. Here, we performed a detailed mechanistic study for 1,2-dichloroethane hydrodechlorination on Pt using a synergistic approach combining density functional theory (DFT) calculations, reaction kinetics experiments, and microkinetic modeling. Using planewave DFT calculations, we evaluated the reaction energy and activation energy barrier of each elementary step involved in the reaction network on Pt(111). The calculated energetics were then incorporated into a comprehensive mean-field microkinetic model accounting for a total of 65 elementary steps. The model-predicted reaction rates were compared with the results from our reaction kinetics experiments on SiO2-supported Pt catalysts. Our results indicated that the hydrodechlorination of 1,2-dichloroethane on Pt(111) starts with a H-removal step; then, it proceeds through a sequence of alternating dechlorination and dehydrogenation steps until vinylidene (CH 2 C*) is formed; finally, CH 2 C* is hydrogenated to the final product, ethane, sequentially via vinyl (CH 2 CH*), ethylene, and ethyl (CH 3 CH 2 *) intermediates. After model parameter adjustments, we achieved good agreement between our theoretical model and experimental results; the adjustments to the calculated parameters are consistent with the typically anticipated coverage effects. Furthermore, our study offers valuable mechanistic insights, which are useful for improving catalysts for this chemistry.

1,2-dichloroethane↗
DOE OSTI2021

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Role of Lithium Codoping in Enhancing the Scintillation Yield of Aluminate Garnets

The aim of this work is to clarify the scintillation-yield enhancement in Lu YAG : Pr scintillators obtained by Li codoping via integrated study of the valence state of activators, the preferential site occupancy of Li codopants, and defect structures from experimental and theoretical insights. Here, with Li codoping, the light yield and energy resolution of 10 × 10 × 10 mm 3 Lu YAG : Pr samples are improved from 156000 to 248000 photons/MeV, and 5.3 to 4.3% at 662 keV,respectively.The optical absorption spectra indicate that Li codoping does not induce conversion of stable Pr 3 + to Pr 4 + in Lu YAG : Pr single crystals. Based on the formation energies of substitutional and interstitial Li sites using density-functional-theory (DFT) calculations and the 7 Li nuclear magnetic resonance results, it is shown that the Li ions prefer to dominantly occupy the fourfold coordinated interstitial sites and fourfold coordinated Al sites. The systematic analysis of thermoluminescence glow curves, positron annihilation lifetime spectroscopies, and defect formation energies derived from DFT calculations reveals that the concentration of isolated Lu and Al vacancies as dominant acceptor defects is reduced by Li codoping, whilst the shallow Li i interstitial defects and the deep V O oxygen vacancies are introduced simultaneously. We propose that the lowering of hole trapping at defects resulting from Li codoping contributes to the scintillation-yield enhancement.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Coupling computational thermodynamics with density-function-theory based calculations to design L12 precipitates in Fe Ni based alloys

Achieving a high-volume fraction of thermodynamically stable L12-type precipitates that are resistant to coarsening is of great importance for the development of low-cost FeNi based austenitic steels. With the aid of computational thermodynamics, this work designed two model alloys: Fe-37.4Ni-6.1Al-2.9Ti (FNAT) and Fe-45.2Ni-5.9Al-8.5Si (FNAS). Both alloys were designed to contain a similar amount of L1 2 precipitate in Fe-Ni matrix without forming other precipitates. Density-Function-Theory (DFT) calculation was coupled with computational thermodynamics to predict the critical radius at which the precipitates change shape from spherical to cuboidal. The calculation results suggest that critical radius for the FNAT alloy is about two orders of magnitude larger than that for the FeNiAlSi alloy. Phase stability and morphology of the L1 2 precipitates in these two alloys were experimentally investigated through X-ray diffraction, atom probe tomography, and scanning and transmission electron microscopy. The L1 2 precipitates in the FeNiAlSi system were found to be cuboidal and rod shaped, with much larger size than the spherical ones in the FeNiAlTi system, agreeing with the calculation results. This work suggested that coupling computational thermodynamics with DFT calculations can be reliably used to design L1 2 precipitates in FeNi based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

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

Solution-Phase Growth of Cu Nanowires with Aspect Ratios Greater Than 1000: Multiscale Theory

Penta-twinned metal nanowires are finding widespread application in existing and emerging technologies. However, little is known about their growth mechanisms. Here, we probe the origins of chloride- and alkylamine-mediated, solution-phase growth of penta-twinned Cu nanowires from first principles using multi-scale theory. Using quantum density functional theory (DFT) calculations, we characterize the binding and surface diffusion of Cu atoms on chlorine-covered Cu(100) and Cu(111) surfaces. We find stronger binding and slower diffusion of Cu atoms on chlorinated Cu(111) than on chlorinated Cu(100) – a reversal of the trend for bare Cu surfaces. We also probe inter-facet diffusion and find this proceeds faster from Cu(100) to Cu(111) than the reverse. Using the DFT rates for hopping between individual sites at Ångstrom scales, we calculate coarse-grained, inter-facet rates for nanowires of various lengths – up to hundreds of μm – and diameters in the 10-nm range. We predict nanowires with aspect ratios around 100 based on surface diffusion, alone. We also account for the influence of a self-assembled alkylamine layer that covers most of the {100} facets, but is absent or thin and disordered on the {111} facets and in an “end zone” near the {100}-{111} boundary. With an end zone, we predict a wide range of nanowire aspect ratios in the experimental ranges. Our work reveals the mechanisms by which a halide – chloride – promotes the growth of high-aspect-ratio nanowires.

36 MATERIALS SCIENCE↗