Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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

Permeation of niobium through grain boundaries in copper

Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.

36 MATERIALS SCIENCE↗

Permeation of niobium through grain boundaries in copper

Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.

36 MATERIALS SCIENCE↗

Physically Based Modeling of PEMFC Cathode Catalyst Layers: Effective Microstructure and Ionomer Structure–Property Relationship Impacts

Abstract This work presents a pseudo-two-dimensional proton exchange membrane fuel cell (PEMFC) model incorporating Nafion ionomer structure–property relationships in the cathode catalyst layer (CL) to capture and explain losses at low Pt loading. Structural data from neutron reflectometry and thin film Nafion conductivity measurements predict variations in the oxygen diffusion coefficient and ionic conductivity with changing CL ionomer thickness and Pt loading. By including these structure–property relationships, predicted polarization curves agree closely with previously published experimental data from cells with Pt loadings between 0.025 and 0.2 mg/cm2. Results demonstrate that structure–property relationships based on physically measurable ionomer and CL properties provide a feasible interpretation of PEMFC CL phenomena for a range of Pt loadings and help explain previously unaccounted-for losses at low Pt. Results also show that simulations must account for surface species coverage variations in order to properly capture the kinetic losses. Finally, results suggest that an increase in ionomer thickness surrounding the C/Pt surfaces may lead to improved cell performance due to improved ionic conductivity.

Electrochemistry↗

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects to be discovered by future direct imaging surveys. We present a low-resolution, R ∼ 300, M-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K {sub zz}) to estimate the strength of vertical mixing in each object. The K {sub zz} values of cooler gaseous objects are close to their theoretical maximum, and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K {sub zz} values imply that disequilibrium phosphine (PH{sub 3}) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH{sub 3} absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's M-band spectrum by including both CO and water clouds in the atmospheric model.

79 ASTRONOMY AND ASTROPHYSICS↗

Crack formation in strained β-(AlxGa 1−x ) 2 O 3 films grown on (010) β-Ga 2 O 3 substrates

The cracking and local strain relaxation in (010) (Al x Ga 1−x ) 2 O 3 films grown on Ga 2 O 3 substrates are assessed in terms of film composition and thickness. We utilize x-ray diffraction and electron microscopy techniques combined with simulation and modeling to investigate that film cracking on curvature (flatness) has a directly proportional relationship with film thickness and/or aluminum content. Cross section transmission electron microscopy reveals that cracks along both the (001) and (100) cleavage planes penetrate into the substrate. The diffuse scattered intensity observed in reciprocal space maps (RSMs) is directly correlated with the tilt that is introduced due to the change in the deformation conditions near the cracks. While asymmetric RSMs show that the layers are fully strained, the diffuse scattering distribution in reciprocal space can be interpreted to show that the cracking relaxes and locally tilts the lattice ∼350 nm from the crack edges, which is consistent with the larger radius of curvature associated with the films with higher crack densities. For example, a 200 nm (Al 0.13 Ga 0.87 ) 2 O 3 thick film has an average inter-crack spacing of 3.3 µm, so most of the epitaxial layer is fully strained except near the cracks where it deforms elastically and is consistent with the gallium oxide Poisson ratio. A reciprocal space model was developed, which imports the strain and tilt distributions (based on finite element modeling) to match the features observed in the experimental maps. We also note that previous studies involving (Al x Ga 1−x ) 2 O 3 films may show evidence of cracking as observed in their symmetric and asymmetric RSMs.

36 MATERIALS SCIENCE↗

A phase-field study of stainless-steel oxidation from high-temperature carbon dioxide exposure

An electrochemical phase-field model has been developed to investigate the oxidation mechanisms of the 21-2N valve stainless steel alloy exposed to carbon dioxide at 973 K. Three oxide phases observed in oxidation experiments are included in the model: Mn 3 O 4 , Cr 2 O 3 , and MnCr 2 O 4 . Local charge neutrality and the conserved current condition are assumed to include the impact of the electric potential on oxidation, though it is found to be negligible when the system is electrically isolated. The sensitivity of the oxidation processes to the values of the diffusion mobilities is examined, and the oxidation rate is calibrated against experimental data by modifying the sensitive mobilities. We find that both inward oxygen and outward metal diffusion are important for oxidation. Furthermore, by investigating the impact of the order of the initial oxide layers, we find that Cr 2 O 3 serves as a barrier to outward Mn diffusion. Moreover, controlling the initial oxide layer order can be used to limit oxidation.

36 MATERIALS SCIENCE↗

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

58 GEOSCIENCES↗

Dissecting Disorder: Defect-Driven Structural Complexity in Layered Li3InCl6 Solid Electrolyte

Halide solid electrolytes have emerged as promising candidates for solid-state batteries owing to their high oxidative stability and ionic conductivity. Among them, Li3InCl6 (LIC) has attracted significant attention. However, diffraction patterns of LIC synthesized via different methods exhibit distinct differences particularly at low-angle reflectionsindicative of underlying structural disorder. These variations are attributed to deviations from ideal crystallographic order, especially stacking faults, whose impact on structure and ion transport remains poorly understood. Here, we identify and quantify stacking faults in LIC samples prepared under different synthetic conditions. Using X-ray diffraction and time-of-flight neutron diffraction, we construct and refine stacking fault models that accurately reproduce the experimental diffraction features. LIC samples with higher degrees of stacking faults exhibit only negligible differences in ionic conductivities and activation energies. This indicates that stacking faults have a limited impact on altering the Li+ diffusion pathway along the c-axis, likely due to the high concentration of vacancies in the In layers, while Li+ diffusion remains nearly unchanged in the ab-plane. Our results account for the observed differences in diffraction patterns across samples and provide a quantitative assessment of faulting probabilities and stacking sequences. The insights gained from this study are expected to be broadly applicable to other layered halide solid electrolytes and contribute to a deeper understanding of the role of structural disorder in ion transport.

Liu, Jue [ORNL] (ORCID:000000024453910X)↗

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗

Effective Lanthanide Diffusivity through U-Zr Metallic Fuel

Fuel-cladding chemical interaction (FCCI) is one of the main concerns for the performance of U-Zr metallic fuels. Although FCCI results in complex phase transformations and morphological changes in both the fuel and cladding sides of the interface, the formation of a brittle layer in the cladding side of the interface (sometimes referred to as wastage) has one of the most significant impacts on the cladding mechanical integrity. This brittle layer is associated with the formation of intermetallic compounds between the cladding constituents Fe, Cr and the lanthanide fission product species. The formation of this cladding wastage layer has been the focus of recent efforts to develop a mechanistic modeling framework, including production of lanthanides in the fuel, transport through the fuel to the fuel-cladding interface, and phase transformation to form intermetallic phases. To calculate the rate of transport of lanthanides through the fuel, a multi-scale computational approach has been used. Atomistic calculations have been performed to determine the diffusivities of Nd, the most prevalent lanthanide species, through the solid fuel matrix and along pore surfaces. These diffusivities have been used as input parameters for a mesoscale model to determine an effective diffusion coefficient that accounts for the formation of porosity, its interconnection, and infiltration with bond sodium from the fuel-cladding gap. The effective diffusivity has been used in engineering-scale simulations using the BISON fuel performance code, and BISON simulations have been compared to representative validation cases from the EBR-II reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Conserving asphalt resources: Rethinking rejuvenator performance evaluation through peptizing efficiency

Timely rejuvenation and restoration of asphalt are essential conservation practices that help preserve and extend the service life of roads, bridges, driveways, and parking lots. The performance of asphalt rejuvenators is often assessed based on their diffusion rates and softening power, yet these metrics alone fail to capture true rejuvenation potential. Here, this study integrates density functional theory (DFT) modeling with experimental analysis to show that rejuvenation effectiveness is primarily governed by molecular interactions with oxidized asphaltene nanoaggregates. DFT results revealed that amide- and unsaturated-chain compounds, such as hexadecanamide, oleic acid, and 9,17-octadecadienal, act cooperatively to disrupt π–π stacking and exfoliate asphaltene layers, reducing binding strength and enhancing dispersion. Unlike bulky, rigid molecules that remain trapped, these components remain mobile and repeatedly interact with multiple aggregation sites. Experimental validation using rheometry and FTIR confirmed that such molecularly compatible rejuvenators restore the binder flexibility and polydispersity, even when diffusion is relatively slow. Building on this mechanistic foundation, six rejuvenators (A2, A5, A7, A8, A9, and A10) were evaluated using a three-metric performance framework encompassing cracking resistance (Glover–Rowe parameter), UV stability, and surface hydrophobicity. Among the six rejuvenators evaluated, A9 exhibited the highest overall performance, reducing the Glover–Rowe cracking parameter by 85% (from 351 kPa to 53 kPa), demonstrating the greatest resistance to UV-induced aging with a stability index of 4.17 h·kPa⁻¹, and increasing surface hydrophobicity to a contact angle of 103.6°. Its superior performance is primarily attributed to its amide- and unsaturated-chain components, which act cooperatively to disrupt π–π stacking interactions and exfoliate asphaltene layers, thereby promoting molecular deagglomeration and enhancing long-term durability. These results shift the criteria for selecting rejuvenators: effective candidates must pair the electronic capability to unlock aged asphaltenes with sufficient structural stability to resist secondary aging and restore hydrophobicity. Restoring hydrophobicity is critical, as aging reduces the asphalt’s water repellency and increases water diffusion, which in turn accelerates moisture-related damage; consequently, a high-performing rejuvenator must effectively restore the binder’s hydrophobic characteristic. Collectively, these findings provide a framework for the rational design of next-generation of bio-based rejuvenators that enhance pavement longevity and promote long-term sustainability.

Aging↗

Diffusion characteristics of water molecules in a lamellar structure formed by triblock copolymers

The distribution and diffusion of water molecules are playing important roles in determining self-assembly and transport properties of polymeric systems. Small-angle neutron scattering (SANS) experiments and molecular dynamics (MD) simulation have been applied to understand the distribution of water molecules and their dynamics in the lamellar membrane formed by Pluronic L62 block copolymers. Penetration of water molecules into the polyethylene oxide (PEO) layers of the membranes has been estimated using scattering length density (SLD) profiles obtained from SANS measurements, which agree well with the molecular distribution observed from MD simulations. The water diffusion coefficient at different regions of the lamellar membrane was further investigated using MD simulation. The diffusion characteristic shows a transition from normal to anomalous diffusion as the position of the water molecule changes from the bulk to PEO and to the polypropylene oxide (PPO) layer. Here we find that water molecules within the PEO or PPO layers follow subdiffusive dynamics, which can be interpreted by the model of fractional Brownian motion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced battery modeling for interfacial phenomena and optimal charging

Lithium ion batteries are one of the most promising energy storage systems for portable devices, transportation, and renewable grids. To meet the increasing requirements of these applications, higher energy density and areal capacity, long cycle life, fast charging rate and enhanced safety for lithium ion battery (LIBs) are urgently needed. To solve these challenges, the relevant physics at different length scale need to be understood. However, experimental study is time consuming and limited in small scale’s study. Modeling techniques provide us powerful tools to get a deep understanding of the relevant physics and find optimal solutions. This work focuses on studying the mechanism in advanced battery engineering techniques and developing a new charging algorithm by model-based optimization. The research topics are divided into six topics and each topic is reported as a form of journal publication. Paper Ⅰ provides a new aspect of how ALD coating improves the lithium ion diffusion at electrode particles. Paper Ⅱ explains the mechanisms by which 3D electrodes enhance battery performance and reveals guidelines for optimized 3D electrode designs by a 3D electrochemical-mechanical battery model. Paper Ⅲ investigates the electrolyte concentration impact on SEI layer growth and Li plating, especially under high charge rates. Paper Ⅳ proposes an optimized charging protocol for fast charging for reducing the charging time with minimal degradation. Paper Ⅴ reports a comprehensive degradation model for degradation estimation and life predication of energy storage system (ESS). Paper Ⅵ is a study of temperature-dependent state of charge (SOC) estimation for battery pack.

25 ENERGY STORAGE↗

Effects of various parameters of different porous transport layers in proton exchange membrane water electrolysis

Porous transport layers (PTLs) play an important role in proton exchange membrane water electrolysis (PEMWE) cells. The PTL facilitates water and gas transport, as well as thermal and electrical conduction, and is required to sustain good contact with adjacent components. It is expected that using PTLs with variations in material properties such as structure, composition, thickness and wettability results in performance changes of the PEMWE. Here, a general mathematical PEMWE model is developed that separates and analyzes the contributions of ohmic, activation, diffusion and Nernst potentials. For model validation, three inherently different anode PTL structures (carbon paper, sintered titanium particles, and titanium felt) are operated over a range of conditions. Additionally, the effects of PTL wettability were used to verify the model using Polytetrafluoroethylene (PTFE) treated Toray papers with PTFE loading ranging from 0% to 20%. The modeling results of both PTFE treated and untreated materials show good agreement with the experimental data. Mass transport or diffusion loss is the primary reason for performance differences between PTFE treated and untreated PTLs. Sintered titanium PTLs with thicknesses above 1 mm suffer from up to 33% increased ohmic losses without indicating any obvious changes in activation and diffusion losses when compared to untreated PTLs. The losses of the cell increase when using PTFE treated Toray paper. Individual contributions are quantified and assigned to increased ohmic, activation, and diffusion losses. In conclusion, the proposed model offers insights into the overpotential contributions of a PEMWE. It is a useful tool for predicting performance of various PTL materials and can be applied for PTL development and optimization efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal Shape Stability of fcc Metal Nanocrystals Synthesized with Faceted Nonequilibrium Shapes

Highly refined capabilities of the shape-controlled solution-phase synthesis of metal nanocrystals (NCs) allow the generation of NCs with faceted nonequilibrium shapes, which optimize properties for target applications such as catalysis and plasmonics. Often, for such applications and also for TEM analysis, the NCs are removed from the solution-phase environment. We explore the postsynthesis evolution of these metastable NCs in a high-vacuum TEM environment. Specifically, here we analyze their reshaping toward the equilibrium Wulff shapes mediated by surface diffusion, where such reshaping degrades the above-mentioned optimized properties. Typical sizes for these NCs range from 5 to 30 nm or 10 3 –10 6 atoms, and reshaping often occurs on the time scale of minutes for temperatures around, say, 400 °C. We discuss the development of predictive stochastic atomistic-level models for NC evolution with a realistic description of surface diffusion. These models, in contrast to Molecular Dynamics, can naturally address the relevant time and length scales for these systems. KMC simulation results for the stochastic models are described, focusing on the reshaping of slightly elongated nanorods and of mildly truncated octahedra and nanocubes. In addition, we review appropriate theoretical formulations for reshaping, which involves the nucleation and growth on 2D islands or layers on outer facets of the NC. We note the limitations of classical nucleation theory in some scenarios and demonstrate the successes of a more fundamental and general master equation-based analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗