Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogenization theory”

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 217 records · Page 12

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗

Suppressing Phase Segregation and Improving Stability in Mixed-Halide Perovskites through Spinel Oxide-Directed Epitaxy

Mixed-cation mixed-halide perovskite compositions are essential for achieving the required bandgaps for high-efficiency multijunction photovoltaics, yet their stability remains limited by interfacial defects, phase segregation, and degradation. Here, we introduce spinel oxides as a new family of lattice-matched substrates that enable crystalline, phase-pure, compositionally-uniform, bromide-rich perovskite film growth. The effect of spinel oxides is two-fold: reducing defects at the bottom interface by templating film growth and inducing beneficial compressive strain through mismatch-dependent substrate-perovskite lattice coupling. Spinel oxide substrates facilitate growth of highly crystalline films and eliminate detrimental secondary phases across thicknesses. Using grazing incidence X-ray diffraction, X-ray fluorescence, cathodoluminescence–scanning electron microscopy, cryogenic photoluminescence, and density functional theory, we reveal that Mg-halide bonds at the bottom interface induce lattice mismatch-dependent compressive strain that suppresses halide segregation and further reduces defect formation. In addition, films grown on spinel oxides maintain over 87% of the perovskite phase after 12 h under 100% relative humidity, as monitored by in situ grazing incidence wide-angle X-ray scattering (GIWAXS), compared to less than 70% for control samples. This work extends lattice matching from vapor-deposited epitaxial semiconductors to solution-processed halide perovskites to establish a broadly applicable strategy for defect suppression, phase homogenization, and long-term stability. Based on the fundamental science explored here, we set the stage for the development of lattice-matched spinel oxide charge transport layers to be integrated into perovskite solar cells and other optoelectronic devices.

36 MATERIALS SCIENCE↗

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗

A Micro-Mechanical Constitutive Model to Predict Hygrothermal Aging of Cross-Linked Polymers

Abstract A multi-physics material model is presented to describe the effects of temperature, oxygen, and humidity on the constitutive response of cross-linked polymers. The effect of hygrothermal damage on the mechanical integrity of the polymer matrix can be considered as the result of damage accumulation of two independent aging mechanism namely, i) thermo-oxidative, and ii) hydrolytic aging. In order to capture the mutual effects of thermo-oxidative and hydrolytic aging, an assumption has been made that each of the aging phenomenon can be superposed to each other. In fact, each of them works independently and as a result, they can compete with each other. Utilizing the theory of network decomposition, all phenomena and their correlation were modeled and thus, the strain energy function of the polymer matrix is written with respect to four independent mechanisms, i) the shrinking original matrix that has neither been attacked by water nor oxygen, ii) conversion of the first network to two new network due to the reduction and formation of cross-links, and iii) energy loss from network degradation due to attack of the water molecules to polymer active agents. Moreover, the proposed model is micro-mechanically based and is mainly relevant on thin samples due to our underlying assumption of homogeneous diffusion of oxygen and water throughout the matrix. The model has been validated against extensive data-sets obtained from experiments we specifically designed for concept validation.

Bahrololoumi, Amir↗

An LES-based airborne Doppler lidar simulator and its application to wind profiling in inhomogeneous flow conditions

Wind profiling by Doppler lidar is common practice and highly useful in a wide range of applications. Airborne Doppler lidar can provide additional insights relative to ground-based systems by allowing for spatially distributed and targeted measurements. Providing a link between theory and measurement, a first large eddy simulation (LES)-based airborne Doppler lidar simulator (ADLS) has been developed. Simulated measurements are conducted based on LES wind fields, considering the coordinate and geometric transformations applicable to real-world measurements. The ADLS provides added value as the input truth used to create the measurements is known exactly, which is nearly impossible in real-world situations. Thus, valuable insight can be gained into measurement system characteristics as well as retrieval strategies. As an example application, airborne Doppler lidar wind profiling is investigated using the ADLS. For commonly used airborne velocity azimuth display (AVAD) techniques, flow homogeneity is assumed throughout the retrieval volume, a condition which is violated in turbulent boundary layer flow. Assuming an ideal measurement system, the ADLS allows to isolate and evaluate the error in wind profiling which occurs due to the violation of the flow homogeneity assumption. Overall, the ADLS demonstrates that wind profiling is possible in turbulent wind field conditions with reasonable errors (root mean squared error of 0.36 m s –1 for wind speed when using a commonly used system setup and retrieval strategy for the conditions investigated). Nevertheless, flow inhomogeneity, e.g., due to boundary layer turbulence, can cause an important contribution to wind profiling error and is non-negligible. Results suggest that airborne Doppler lidar wind profiling at low wind speeds (<5 m s –1 ) can be biased, if conducted in regions of inhomogeneous flow conditions.

17 WIND ENERGY↗

Quantifying the hierarchy of structural and mechanical length scales in granular systems

Continuum modeling of granular media is made possible by the existence of a length scale at and above which grain-resolved properties can be meaningfully homogenized. Progress has been made in identifying such length scales relevant to local structural properties such as porosity. However, a systematic analysis of scales above which different mechanical properties can be homogenized has yet to emerge. Here, X-ray tomography and 3D X-ray diffraction data are examined to identify such length scales. The data was obtained in-situ in compressed granular materials with rigid and flexible confinement. The experimental data are supplemented with validated discrete element simulations which examine different system sizes and different boundary conditions. Overall, our study reveals a hierarchy in the length scales of granular solids, with lengths governing structural variables being the shortest, lengths of stress variables being intermediate, and lengths of energy dissipation being the longest. All structural and mechanical length scales obey a power law based on the theory of Geostatistics, implying that the length scales can be found by analyzing samples significantly smaller than the length scales themselves. The length scales are also found to be sensitive to boundary conditions, implying that they are extrinsic features of granular media.

36 MATERIALS SCIENCE↗

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

Supported Single‐Atom Manganese Catalysts for the Trimerization of Ethylene

Selective ethylene oligomerization via oxidative cyclization, forming metallacyclic intermediates, is typically catalyzed by molecular titanium and chromium complexes to produce butenes, hexenes, or octenes, depending on the supporting ligand framework. However, this mechanism requires significant electron density at the metal active site and is not known to be generalizable to other first-row transition metals. In this work, we computationally investigate the electronic modulation of five transition metals (Mn, Fe, Co, Ni, and Cu) supported on titania (TiO₂) through reductive lithium intercalation to promote selective oligomerization via oxidative cyclization, using density functional theory (DFT). Our findings predict that Mn/LiTiO₂ exhibits high catalytic activity due to the exergonic nature of oxidative cyclization with two ethylene molecules. Additionally, lithium titanate (LiTiO₂) supports enhance catalytic performance compared to TiO₂. Experimental validation confirms that Mn/LiTiO₂ achieves higher conversion rates and improved selectivity toward hexene (C₄:C₆ = 1:2.6). The enhanced activity is attributed to lithiation, which alters the electronic environment around Mn active sites. Mechanistic studies reveal that the formation of a seven-membered ring, a key intermediate for hexene formation, is more favorable on LiTiO₂ than TiO₂. This work provides the first evidence of Mn catalyzing selective ethylene oligomerization via oxidative cyclization in either homogeneous or heterogeneous catalysis.

Kim, Yu Lim [Argonne National Laboratory (ANL), Ar↗

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗

A Novel NASICON‐Type Na 4 MnCr(PO 4 ) 3 Demonstrating the Energy Density Record of Phosphate Cathodes for Sodium‐Ion Batteries

Abstract Sodium‐ion batteries (SIBs) have attracted incremental attention as a promising candidate for grid‐scale energy‐storage applications. To meet practical requirements, searching for new cathode materials with high energy density is of great importance. Herein, a novel Na superionic conductor (NASICON)‐type Na 4 MnCr(PO 4 ) 3 is developed as a high‐energy cathode for SIBs. The Na 4 MnCr(PO 4 ) 3 nanoparticles homogeneously embedded in a carbon matrix can present an extraordinary reversible capacity of 160.5 mA h g −1 with three‐electron reaction at ≈3.53 V during the Na + extraction/insertion process, realizing an unprecedentedly high energy density of 566.5 Wh kg −1 in the phosphate cathodes for SIBs. It is intriguing to reveal the underlying mechanism of the unique Mn 2+ /Mn 3+ , Mn 3+ /Mn 4+ , and Cr 3+ /Cr 4+ redox couples via X‐ray absorption near‐edge structure spectroscopy. The whole electrochemical reaction undergoes highly reversible single‐phase and biphasic transitions with a moderate volume change of 7.7% through in situ X‐ray diffraction and ex situ high‐energy synchrotron X‐ray diffraction. Combining density functional theory (DFT) calculations with the galvanostatic intermittent titration technique, the superior performance is ascribed to the low ionic‐migration energy barrier and desirable Na‐ion diffusion kinetics. The present work can offer a new insight into the design of multielectron‐reaction cathode materials for SIBs.

Zhang, Jian↗

Scale-free, programmable design of morphable chain loops of kilobots and colloidal motors

Micron-scale robots require systems that can morph into arbitrary target configurations controlled by external agents such as heat, light, electricity, and chemical environment. Achieving this behavior using conventional approaches is challenging because the available materials at these scales are not programmable like their macroscopic counterparts. To overcome this challenge, we propose a design strategy to make a robotic machine that is both programmable and compatible with colloidal-scale physics. Our strategy uses motors in the form of active colloidal particles that constantly propel forward. We sequence these motors end-to-end in a closed chain forming a two-dimensional loop that folds under its mechanical constraints. We encode the target loop shape and its motion by regulating six design parameters, each scale-invariant and achievable at the colloidal scale. We demonstrate the plausibility of our design strategy using centimeter-scale robots called kilobots . We use Brownian dynamics simulation to explore the large design space beyond that possible with kilobots, and present an analytical theory to aid the design process. Multiple loops can also be fused together to achieve several complex shapes and robotic behaviors, demonstrated by folding a letter shape “M,” a dynamic gripper , and a dynamic pacman . The material-agnostic, scale-free, and programmable nature of our design enables building a variety of reconfigurable and autonomous robots at both colloidal scales and macroscales.

36 MATERIALS SCIENCE↗

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗

Photoinduced electron transfer from zinc meso -tetraphenylporphyrin to a one-dimensional perylenediimide aggregate: Probing anion delocalization effects

Organic photovoltaics incorporating non-fullerene acceptors based on perylenediimide (PDI) now rival fullerene acceptor-based devices in performance, although the mechanisms of charge generation in PDI-based devices are not yet fully understood. Fullerene-based systems are proposed to undergo electron transfer directly from the photoexcited donor into a band of delocalized acceptor states, thus increasing charge generation efficiency. Similarly, anion delocalization has been shown to enhance the rate of electron transfer from a photoexcited donor to two electronically coupled PDI acceptors. Here we investigate how additional electron acceptors may further increase the rate of electron transfer from the donor zinc meso-tetraphenylporphyrin (ZnTPP) to an aggregate of PDI acceptors (PDI[Formula: see text]. Femtosecond transient visible and mid-infrared absorption spectroscopies show that the rate of electron transfer from 1 *ZnTPP to the PDI assembly ZnTPP 2 -PDI 3 is statistically identical to that of the previously examined ZnTPP-PDI 2 . A Marcus theory analysis indicates that the parameters governing electron transfer are nearly identical for the two molecules, suggesting that the maximum electron transfer rate enhancement has been achieved in a cofacial PDI dimer because the ZnTPP directly couples to the first two PDI acceptors whereas the coupling to the third PDI is too weak.

bio-inspired↗

Chapter 21. Atomic-Level Dislocation Dynamics in Irradiated Metals

Primary damage and microstructure evolution in structural nuclear materials operating under conditions of a high flux of energetic atomic particles and high temperature and stress lead to the formation of a high concentration, non-homogeneous distribution of defect clusters in the form of dislocation loops, voids, gas-filled bubbles and radiation-induced precipitates of nanometer scale. They cause changes in many material properties. Being obstacles to dislocation glide, they strongly affect the mechanical properties, with an increase in yield and flow stresses and a reduction in ductility. Atomic-scale computer simulations can provide details of how these effects are influenced by the obstacle structure, applied stress, strain rate and temperature. Processes such as obstacle cutting, transformation, absorption and drag are observed. Some recent results for body-centered and face-centered cubic metals are described in this review and, where appropriate, comparisons are drawn with predictions based on the elasticity theory of crystal defects. Perspectives on how to use this information at higher scales and in particular in mesoscale, dislocation dynamics simulations are also discussed.

Osetskiy, Yury N.↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective field theory for quasicrystals and phasons dynamics

We build an effective field theory (EFT) for quasicrystals -- aperiodic incommensurate lattice structures -- at finite temperature, entirely based on symmetry arguments and a well-define action principle. By means of Schwinger-Keldysh techniques, we derive the full dissipative dynamics of the system and we recover the experimentally observed diffusion-to-propagation crossover of the phason mode. From a symmetry point of view, the diffusive nature of the phason at long wavelengths is due to the fact that the internal translations, or phason shifts, are symmetries of the system with no associated Noether currents. The latter feature is compatible with the EFT description only because of the presence of dissipation (finite temperature) and the lack of periodic order. Finally, we comment on the similarities with certain homogeneous holographic models and we formally derive the universal relation between the pinning frequency of the phonons and the damping and diffusion constant of the phason.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗