Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized continuum”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2023

Infrared spectroscopy of nearby radio active early-type galaxies - II: spectral atlas

ABSTRACT We present a near-infrared spectroscopic atlas of nearby, bright early-type galaxies with radio emission, containing 163 galaxies observed by the Palomar 200″ TripleSpec instrument, measuring the emission line fluxes for H, He, [Fe ii], and H2 and determined the nuclear excitation mechanisms. By stacking spectra, we deduced the H2 excitation temperature (1957 ± 182 K) and dominant excitation mechanism (thermal and shock heating combined) from the K-band emission line sequence. Stacking also produces an ‘average’ spectrum of absorption features and spectral indices from the literature; the CO12 absorption line index versus J − K colour shows a trend of stronger nuclear activity producing a weaker CO12 index and a redder (flatter) continuum. The correlations between the radio and the emission-line luminosities finds a trend with radio power; however, the large scatter in the upper limits shows that the two are not directly coupled and the duty cycles of SF and AGN activity are not synchronized.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Radiative Electronic Bound States in the Continuum from Defects in Semiconductors

Continuum-buried defect states in semiconductors are generally expected to be optically inactive because of their strong coupling to continuum bands. Here, we show that such defects can instead host radiative electronic bound states in the continuum (BICs) using the silicon G center as a prototypical example. Hybrid functional first-principles calculations with a Hubbard U correction reveal that a localized defect state, initially buried below the valence band maximum (VBM) in the ground state, undergoes exchange-driven energy-level reordering under optical excitation and shifts above the VBM. This exchange-induced transition suppresses nonradiative decay and enables a robust radiative emission. By computing temperature-dependent nonradiative lifetimes and comparing them to experimental photoluminescence (PL) lifetimes, we quantitatively reproduce the observed temperature dependence of the emission. These results uncover a stabilization mechanism for continuum-embedded defect states and establish electronic BICs as a general paradigm for designing defect-based optical systems, including quantum emitters and qubits.

Hong, SeongYun↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast hole relaxation dynamics in quantum dots revealed by two-dimensional electronic spectroscopy

Elucidating the population dynamics of correlated electron-hole pairs (bound excitons) in semiconducting quantum dots (QDs) is key for developing our fundamental understanding of nanoscale photophysics as well as for the optimal design of devices, such as lasers. For decades, it was assumed that holes did not contribute to band edge bleach signals in QDs. Here, we employ two-dimensional electronic spectroscopy to monitor electron and hole dynamics in both CdSe and CdSe/CdS/ZnS QDs to probe electron and hole dynamics. Based on a combination of time and frequency resolution, we observe a previously unresolved bleaching signal in CdSe QDs on timescales faster than 30 fs due to hole cooling. Atomistic semiempirical pseudopotential calculations are used to rationalize the order of magnitude difference in the observed hole dynamics in CdSe and CdSe/CdS/ZnS QDs. This picture advances our understanding of QD excitonics past the prevailing continuum effective mass theories generally used to describe QD electronic structure and dynamics.

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

Field-tuned quantum renormalization of spin dynamics in the honeycomb lattice Heisenberg antiferromagnet YbCl 3

The basis for our understanding of quantum magnetism has been the study of elegantly simple model systems. However, even for the antiferromagnetic honeycomb lattice with isotropic spin interactions–one of the simplest model systems–a detailed understanding of quantum effects is still lacking. Here, using inelastic neutron scattering measurements of the honeycomb lattice material YbCl 3 , we elucidate how quantum effects renormalize the single-magnon and multimagnon excitations and how this renormalization can be tuned and ultimately driven to the classical limit by applying a magnetic field. Additionally, our work reveals that the quantum effects tuned by the magnetic field not only renormalize the magnetic excitations but also induce a distinctive sharp feature inside the multimagnon continuum. From a more general perspective, this result demonstrates that structures within magnetic continua can occur over a wide experimental parameter space and can be used as a reliable means of identifying quantum phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Abell 1430: A merging cluster with exceptional diffuse radio emission

Diffuse radio emission has been found in many galaxy clusters, predominantly in massive systems which are in the state of merging. The radio emission can usually be classified as relic or halo emission, which are believed to be related to merger shocks or volume-filling turbulence, respectively. Recent observations have revealed radio bridges for some pairs of very close galaxy clusters. The mechanisms that may allow one to explain the high specific density of relativistic electrons, which are necessary to explain the radio luminosity of these bridge regions, have been poorly explored until now. When inspecting the first data release of the LOFAR Two-Metre Sky Survey (LoTSS), we discovered diffuse radio emission in the galaxy cluster Abell 1430. Here, for this work, we aim to determine the dynamical state of the cluster and characterise the diffuse radio emission. We analysed the LoTSS data in detail and complemented them with recent Karl G. Jansky Very Large Array observations in the L-band. To study the dynamical state of the cluster, we analysed XMM-Newton data, Chandra data, and Sloan Digital Sky Survey data. Moreover, we compared our results to clusters extracted from THE THREE HUNDRED PROJECT cosmological simulation. We find that Abell 1430 consists of two components, namely A1430-A and A1430-B, with a mass ratio of about 2:1. The massive component shows diffuse radio emission which can be classified as radio halo which shows a low radio power at 1.4 GHz with respect to the mass of the cluster. Most interestingly, there is extended diffuse radio emission in the following dubbed as the ‘Pillow’ according to its morphology, which is apparently related to A1430-B and which is neither typical halo nor typical relic emission. The origin of this emission is puzzling. We speculate that the two components of Abell 1430 undergo an off-axis merger. In this scenario, A1430-B is moving towards the main cluster component and may have compressed and stirred the medium in the filament between the two cluster components. We have discovered evidence for diffuse radio emission related to the low-density intracluster or intergalactic medium in Abell 1430. To date, only a few examples of emission originating from such regions are known. These discoveries are crucial to constrain possible acceleration mechanisms which may allow us to explain the presence of relativistic electrons in these regions. In particular, our results indicate a spectral index of α 144 MHz 1.5 GHz = -1.4±0.5 for the Pillow. If upcoming observations confirm a slope as flat as -1.4 or even flatter, this would pose a challenge for the electron acceleration scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

How plane are plane shock waves in solids

The notion of plane shock waves is a macroscopic, very fruitful idealization of near discontinuous disturbance propagating at supersonic speed. Such a picture is comparable to the picture of shorelines seen from a very high altitude. When viewed at the grain scale where the structure of solids is inherently heterogeneous and stochastic, features of shock waves are non-laminar and field variables, such as particle velocity and pressure, fluctuate. This paper reviews select aspects of such fluctuating nonequilibrium features of plane shock waves in solids with focus on grain scale phenomena and raises the need for a paradigm change to achieve a deeper understanding of plane shock waves in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Temporal dynamics of the alignment of the turbulent stress and strain rate

We report that the flux of energy between scales in turbulence can be cast as the result of the interaction between a turbulent stress and a rate of strain. Efficient energy transfer requires that these two tensors be oriented properly relative to each other, but previous work has shown that the instantaneous alignment between them is poor. Here, we consider the temporal dynamics of this alignment in a direct numerical simulation of isotropic turbulence. We show that the orientation of the stress lags behind that of the strain rate, both at a single location and along trajectories. However, the timescale of the reorientation of the stress in these cases does not follow the expected dynamical scaling with length scale. To capture the proper dynamical scaling, we reformulate the energy flux between scales using the right Cauchy-Green strain tensor and the second Piola-Kirchhoff stress tensor. Our results highlight the key role played by the deformation of fluid elements in the physics of the energy cascade, and suggest that their irreversible deformation is a Lagrangian manifestation of the cascade's broken time-reversal symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Density matrix renormalization group description of the island of inversion isotopes F 28 – 33

Recent experiments have confirmed that the neutron-rich isotopes 28,29 F belong to the so-called island of inversion (IOI), a region of the nuclear chart around Z = 10 and N = 20 where nuclear structure deviates from the standard shell model predictions due to deformation and continuum effects. However, while the general principles leading to the IOI are relatively well understood, the details of the low-lying structure of the exotic fluorine isotopes 28–33 F are basically unknown. In this study, we perform large-scale shell model calculations including continuum states to investigate the properties of the neutron-rich isotopes 25–33 F, from a core of 24 O and using an effective two-body interaction with a small number of adjustable parameters in the central and tensor channels. We develop two models adjusted on experimentally confirmed states in 25,26 O and 25–27 F based on different assumptions concerning the positions of the neutron 0d 3/2 and 1p 3/2 shells, and solve the many-body problem using the density matrix renormalization group (DMRG) method for open quantum systems in an sd–fp model space. We obtain the first detailed spectroscopy of 25–33F in the continuum and show how the interplay between continuum effects and deformation explains the recent data on 28,29 F. Several deformed one- and two-neutron halo states are predicted in 29,31 F, and we provide some information about the possible structure of the heaviest fluorine isotopes. We also suggest several experimental studies of interest to constraint models and test the present predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

MIGHTEE: total intensity radio continuum imaging and the COSMOS/XMM-LSS Early Science fields

ABSTRACT MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectropolarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image ∼20 deg2 over the COSMOS, E-CDFS, ELAIS-S1, and XMM-Newton Large Scale Structure field (XMM-LSS) extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength data sets from numerous existing and forthcoming observational campaigns. Here, we describe and validate the data processing strategy for the total intensity continuum aspect of MIGHTEE, using a single deep pointing in COSMOS (1.6 deg2) and a three-pointing mosaic in XMM-LSS (3.5 deg2). The processing includes the correction of direction-dependent effects, and results in thermal noise levels below 2 $\mathrm{\mu }$Jy beam−1 in both fields, limited in the central regions by classical confusion at ∼8 arcsec angular resolution, and meeting the survey specifications. We also produce images at ∼5 arcsec resolution that are ∼3 times shallower. The resulting image products form the basis of the Early Science continuum data release for MIGHTEE. From these images we extract catalogues containing 9896 and 20 274 radio components in COSMOS and XMM-LSS, respectively. We also process a close-packed mosaic of 14 additional pointings in COSMOS and use these in conjunction with the Early Science pointing to investigate methods for primary beam correction of broad-band radio images, an analysis that is of relevance to all full-band MeerKAT continuum observations, and wide-field interferometric imaging in general. A public release of the MIGHTEE Early Science continuum data products accompanies this article.

Heywood, I. (ORCID:0000000168645057)↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

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

Holographic codes from hyperinvariant tensor networks

Holographic quantum-error correcting codes are models of bulk/boundary dualities such as the anti-de Sitter/conformal field theory (AdS/CFT) correspondence, where a higher-dimensional bulk geometry is associated with the code’s logical degrees of freedom. Previous discrete holographic codes based on tensor networks have reproduced the general code properties expected from continuum AdS/CFT, such as complementary recovery. However, the boundary states of such tensor networks typically do not exhibit the expected correlation functions of CFT boundary states. In this work, we show that a new class of exact holographic codes, extending the previously proposed hyperinvariant tensor networks into quantum codes, produce the correct boundary correlation functions. This approach yields a dictionary between logical states in the bulk and the critical renormalization group flow of boundary states. Furthermore, these codes exhibit a state-dependent breakdown of complementary recovery as expected from AdS/CFT under small quantum gravity corrections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Field-driven ion pairing dynamics in concentrated electrolytes

Here, we investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5M LiPF 6 in water and acetonitrile under uniform electric fields, we compute transition path theory observables, including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of 40% at 50 mV/Å in acetonitrile, compared to that of less than 10% in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

Moon, Seokjin [University of California, Berkeley,↗
DOE OSTI2023

Dissociation dynamics of ion-pair states accessed by low-energy electron collisions

Abstract Ion-pair (IP) states are the superexcited states of the neutral molecule that dissociate as a pair of positive and negative ions. These states are generally present near the ionization continuum of the molecule. IP states can be accessed by using photoexcitation or by using the electron collision technique. Different photoexcitation techniques are widely used over the years to study the threshold (threshold IP production spectroscopy) and the dynamics (IP imaging spectroscopy) of the IP states. However, the electron collision technique is ignored over the years and only a few studies are available. In this review, we will discuss different experimental techniques to probe IP states by using electron collision and also the dynamics of the IP states that are accessed by the electron collision.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Connecting radio emission to AGN wind properties with broad absorption line quasars

ABSTRACT Broad absorption line quasars (BALQSOs) show strong signatures of powerful outflows, with the potential to alter the cosmic history of their host galaxies. These signatures are only seen in ∼10 per cent of optically selected quasars, although the fraction significantly increases in IR and radio selected samples. A proven physical explanation for this observed fraction has yet to be found, along with a determination of why this fraction increases at radio wavelengths. We present the largest sample of radio matched BALQSOs using the LOFAR Two-metre Sky Survey Data Release 2 and employ it to investigate radio properties of BALQSOs. Within the DR2 footprint, there are 3537 BALQSOs from Sloan Digital Sky Survey DR12 with continuum signal-to-noise ≥5. We find radio-detections for 1108 BALQSOs, with an important subpopulation of 120 LoBALs, an unprecedented sample size for radio matched BALQSOs given the sky coverage to date. BALQSOs are a radio-quiet population that show an increase of ×1.50 radio-detection fraction compared to non-BALQSOs. LoBALs show an increase of ×2.22 that of non-BALQSO quasars. We show that this detection fraction correlates with wind-strength, reddening, and C iv emission properties of BALQSOs and that these features may be connected, although no single property can fully explain the enhanced radio detection fraction. We create composite spectra for subclasses of BALQSOs based on wind strength and colour, finding differences in the absorption profiles of radio-detected and radio-undetected sources, particularly for LoBALs. Overall, we favour a wind-ISM interaction explanation for the increased radio-detection fraction of BALQSOs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗