Search NASASearch

SEARCH · Search NASA

Results for “time scale separation”

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

Hierarchical-embedding autoencoder with a predictor as efficient architecture for learning time-evolution in multi-scale turbulent flows

We introduce a scale-aware, data-driven deep learning modeling framework for accurately predicting the time evolution of multi-scale turbulent plasma and liquid flows. The approach is motivated by the idea of scale separation. Structures of vastly different length scales emerge in these systems, and interactions between these structures occur only locally. To exploit this structure, the flow state is transformed by a hierarchical, fully convolutional autoencoder, not into a single embedding layer as in conventional convolutional surrogate models, but into a series of embedding layers. A stepwise training strategy ensures that fine-scale features are encoded on a high-resolution grid, while larger structures are represented on progressively coarser layers. The time evolution predictor advances all embedding layers in sync, capturing local interactions between features at the same scale as well as between all scales. This approach enables efficient modeling of multi-scale systems since negligible interactions between distant, small-scale structures do not need to be directly modeled. Our hierarchical-embedding autoencoder with a predictor framework is evaluated on canonical examples of multi-scale turbulence: two-dimensional Kolmogorov flow and Hasegawa–Wakatani plasma turbulence. In both cases, the proposed framework significantly improves predictive accuracy relative to conventional convolutional network architectures. A significant improvement in prediction accuracy was observed for crucial statistical characteristics of the Hasegawa–Wakatani plasma as well as for individual trajectories of the Kolmogorov flow turbulence. Importantly, the model's rollout for the Hasegawa–Wakatani problem demonstrates a four-order-of-magnitude speedup compared to traditional numerical solvers.

Khrabry, Alexander I. [Princeton Univ., NJ (United

Hot-Carrier Injection and Millisecond Charge Separation from a Robust Heteroleptic Iron(II) Chromophore Immobilized on TiO 2

The synthesis, spectroscopic characterization, computational analysis, and photoelectrochemical behavior of a new iron-based chromophore, [(Cpy) 2 Fe- (deeb)](PF 6 ) 2 (Fe(Cpy) 2 (deeb)), where Cpy is 1-methyl-3-(2-pyridyl)imidazole and deeb is 4,4’-(CO 2 CH 2 CH 3 ) 2 -2,2’-bipyridine, is reported. Electrochemically reversible waves assigned to a metal-centered E o (Fe III/II ) = +0.48 and a ligand-centered E o (Fe 2+/+ ) = −1.47 V vs Fc +/0 reduction were evident in cyclic voltammetry measurements. The combination of a strong σ-donor and a π-acceptor lowered the energy of the metal-toligand charge-transfer (MLCT) excited state relative to the metal-centered state. Two MLCT transitions appear in the visible region at 424 and 580 nm. TDDFT calculations revealed that the lower-energy band was well formulated as Fe(II)→deeb, and the higher-energy transition was charge transfer to both the deeb and Cpy ligands. Resonance Raman spectroscopy supports these findings showing enhanced deeb vibrational modes with 532 nm excitation, both deeb and Cpy modes with 473 nm excitation, and exclusively Cpy with 405 nm excitation. Ultrafast spectroscopy reveals a short-lived (∼2 ps) MLCT excited state and a longer-lived (∼20 ps) metal-centered state. Efficient methods to deprotect the ester groups and anchor the complex to mesoporous TiO 2 (anatase) thin films in high surface coverages, Fe(Cpy) 2 (dcb)|TiO 2 σ = 3 × 10 −8 mol/cm 2 , were established. Pulsed light excitation of Fe(Cpy) 2 (dcb)|TiO 2 resulted in rapid excited state injection (k inj > 10 8 s −1 ) and formation of a charge-separated state, Fe III (Cpy) 2 (dcb)|TiO 2 (e), which persists on the millisecond time scale before returning cleanly to the ground state with secondorder kinetics. Injection yields measured 50 ns after light excitation were found to double from Φ = 0.15 with green (532 nm) light to 0.30 with blue (457 nm) light excitation. Incident photon-to-current efficiency (% IPCE) measurements as a function of excitation wavelength in a 0.5 M LiI/I 2 /CH 3 CN electrolyte provide clear evidence for band-selective “hot carrier” injection from the remote Cpy-localized excited state. Collectively, the spectroscopic and photoelectrochemical data indicate that a semiconductor can intercept hot electrons from iron chromophores even when the excited-state dipole is oriented away from the surface-anchoring ligand.

Charge transfer

Experimental and theoretical evidence of universality in superfluid vortex reconnections

The minimum separation between reconnecting vortices in fluids and superfluids obeys a universal scaling law with respect to time. The prereconnection and the postreconnection prefactors of this scaling law are different, a property related to irreversibility and to energy transfer and dissipation mechanisms. In the present work, we determine the temperature dependence of these prefactors in superfluid helium from experiments and a numeric model which fully accounts for the coupled dynamics of the superfluid vortex lines and the thermal normal fluid component. At all temperatures, we observe a pre- and postreconnection asymmetry similar to that observed in other superfluids and in classical viscous fluids, indicating that vortex reconnections display a universal behavior independent of the small-scale regularizing dynamics. We also numerically show that each vortex reconnection event represents a sudden injection of energy in the normal fluid. Finally we argue that in a turbulent flow, these punctuated energy injections can sustain the normal fluid in a perturbed state, provided that the density of superfluid vortices is large enough.

Reconnections

Cross-Linker Selection Controls Glass Transition Elevation or Reduction in Dynamic Covalently Bonded Polymer Networks

Introducing cross-links is a powerful approach to improve polymeric material performance relevant to controlling viscoelasticity, thermal and creep resistance, degradability, and efficient membrane separations. The chemically specific glass transition temperature T g is of fundamental importance in determining the time scales of key dynamical processes and physical state of the material in such applications. Here, we study experimentally how the introduction of relatively large cross-linking molecules in slowly exchanging dynamic bond-forming polymers (vitrimers) impacts vitrification for diverse polymer chemistries and a wide range of cross-link fractions. We find T g can increase, decrease, or even remain essentially unchanged, in qualitative contrast to the generic elevation of T g in traditional permanent polymer networks. We formulate an effective terpolymer network model to understand this rich behavior, which emerges as a consequence of a competition between pure cross-linking and generalized plasticization effects. The latter is associated with the tunable cross-linker size and intrinsic dynamic mobility that can offset slowing down due to traditional permanent cross-linking constraints. Here, a new strategy for functional polymer network design is suggested based on adjusting the relative importance of the two competing physical effects, which potentially can significantly enhance energy savings in applications while retaining other intrinsic properties germane to advanced materials performance.

Copolymers

Experiments and gyrokinetic simulations of the nonlinear interaction between spinning magnetized plasma pressure filaments

A set of experiments using controlled, skin depth-sized plasma pressure filaments in close proximity have been carried out in a large linear magnetized plasma device. Two- and three-filament configurations have been used to determine the scale of cross field nonlinear interaction. When the filaments are separated by a distance of approximately five times the size of a single filament or less, a significant transfer of charge and energy occurs, leading to the generation of inter-filament electric fields. This has the effect of rotating the filaments and influencing the merging dynamics. Nonlinear gyrokinetic simulations using seeded filaments confirm the presence of unstable drift-Alfvén modes driven by the steep electron temperature gradient. When the filaments are within a few collisionless electron skin depths (separations twice the size of a single filament), the unstable perturbations drive the convective mixing of the density and temperature and rearrange the gradients such that they maximize in the region surrounding the filament bundle.

Sydora, R. D. (ORCID:0000000192543149)

Detecting dark matter substructures on small scales with fast radio bursts

The matter power spectrum is only weakly constrained on subgalactic scales, while physics beyond the Standard Model can leave unique imprints, especially on sub-parsec scales. We propose measuring the arrival-time difference of fast radio bursts (FRBs) along two adjacent sightlines as a new probe to dark matter substructures on scales down to ∼ 1 AU . We discuss two observational scenarios in which it may be possible to place interesting constraints on such models through the monitoring of repeating FRB sources: (i) By sending radio receivers to space to form a baseline of tens of AU or more and measuring the temporal variation of the arrival-time difference between receivers. (ii) By measuring the temporal variation of the arrival-time difference between two lensed images of one strongly lensed repeater. In both scenarios, obtaining interesting constraints requires correlating the voltage time series to measure the radio signal arrival time to sub-nanosecond precision. We find that two radio dishes separated by 20 AU may be sensitive to the enhancement of small-scale structures at ∼ 10 - 8 M ⊙ masses in the QCD axion dark matter scenario, or from an early epoch of matter domination with a reheating temperature up to 60 MeV. Other dark matter models, such as those composed of ∼ 10 - 13 M ⊙ primordial black holes produced during inflation, would also be probed by this method. We further show that a strong lensing situation of multiple images provides an equivalent ∼ 2000 AU ( σ v / 10 3 km s - 1 ) ( δ t / 10 yr ) baseline, for a typical velocity of dark matter substructures σ v and an observational time span δ . This is much more sensitive, but with the uncertainty that intervening decoherence from the interstellar medium may degrade the timing precision, and that spatial variation in the FRB emission spot may result in confounding signals. We show that the lensing magnifications of Type Ia supernovae constrain a similar quantity to such FRB timing, with present limits being equivalent to ruling out the same parameter space that would be probed by a 0.14 AU baseline.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Modeling and Scaling up of Membrane Modules Leveraging Dimensional Analysis

In the current analysis, we study the separation performance using membrane modules through dimensional analysis (DA). We formulate the main process equations to identify relevant dimensionless numbers inherent in the physics. In particular, we identify that the critical step in the separation process is mass transfer through the selective layer. Remarkably, the dimensionless feed flow (DFfeed) emerges as a crucial factor in describing this process. Not only does DFfeed directly appear in the governing equations, but it also holds a physical significance associated with the time scales for the mass transfer across the feed side and through the selective layer. Regarding the output performance variables, we consider the recovery, stage cut, productivity and purity. In this context, we profit from experimental data and CFD simulations to evaluate the separation performance of the modules when varying the input flowrate, the scale of the module, and the CO2 permeance. These datasets enable us to establish correlations between performance metrics and the dimensionless feed flow. Using simple power functions of DFfeed, we obtain R2 coefficients exceeding 0.99, indicating the accuracy of the correlations built in the present work. In the future, we wish to use DA to understand key transport mechanisms, predict and control module performance, and challenge the universality of these findings by testing various gas separations across different membrane modules beyond our case study.

Pedrozo, Hector A.

Characterization of DESI fiber assignment incompleteness effect on 2-point clustering and mitigation methods for DR1 analysis

We present an in-depth analysis of the fiber assignment incompleteness in the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1). This incompleteness is caused by the restricted mobility of the robotic fiber positioner in the DESI focal plane, which limits the number of galaxies that can be observed at the same time, especially at small angular separations. As a result, the observed clustering amplitude is suppressed in a scale-dependent manner, which, if not addressed, can severely impact the inference of cosmological parameters. We discuss the methods adopted for simulating fiber assignment on mocks and data. In particular, we introduce the fast fiber assignment (FFA) emulator, which was employed to obtain the power spectrum covariance adopted for the DR1 full-shape analysis. We present the mitigation techniques, organised in two classes: measurement stage and model stage. We then use high fidelity mocks as a reference to quantify both the accuracy of the FFA emulator and the effectiveness of the different measurement-stage mitigation techniques. This complements the studies conducted in a parallel paper for the model-stage techniques, namely the θ-cut approach. We find that pairwise inverse probability (PIP) weights with angular upweighting recover the “true” clustering in all the cases considered, in both Fourier and configuration space. Notably, we present the first ever power spectrum measurement with PIP weights from real data.

cosmological simulations

Probing Curved Spacetime with a Distributed Atomic Processor Clock

Quantum dynamics on curved spacetime has never been directly probed beyond the Newtonian limit. Although we can describe such dynamics theoretically, experiments would provide empirical evidence that quantum theory holds even in this extreme limit. The practical challenge is the minute spacetime curvature difference over the length scale of the typical extent of quantum effects. Here, we propose a quantum network of alkaline earth (like) atomic processors for constructing a distributed quantum state that is sensitive to the differential proper time between its constituent atomic processor nodes, implementing a quantum observable that is affected by post-Newtonian curved spacetime. Conceptually, we propose to delocalize one clock between three locations by encoding the presence or absence of a clock into the state of the local atoms. By separating three atomic nodes over approximately kilometer-scale elevation differences and distributing one clock between them via a 𝑊 state, we demonstrate that the curvature of spacetime is manifest in the interference of the three different proper times that give rise to three distinct beat notes in our nonlocal observable. We further demonstrate that 𝑁-atom entanglement within each node enhances the interrogation bandwidth by a factor of 𝑁. We discuss how our proposed system can probe new facets of fundamental physics, such as the linearity, unitarity, and probabilistic nature of quantum theory on curved spacetime. Our protocol combines several recent advances with neutral atom and trapped ions to realize a novel quantum probe of gravity uniquely enabled by quantum networks.

Covey, Jacob P. [Univ. of Illinois at Urbana-Champ

Fermionic mean-field dynamics for spin systems beyond free fermions

We introduce the fermionized time-dependent Hartree–Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions, and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

Dutta, Rishab

Community structure and function during periods of high performance and system upset in a full-scale mixed microalgal wastewater resource recovery facility

Microalgae have the potential to exceed current nutrient recovery limits from wastewater, enabling water resource recovery facilities (WRRFs) to achieve increasingly stringent effluent permits. The use of photobioreactors (PBRs) and the separation of hydraulic retention and solids residence time (HRT/SRT) further enables increased biomass in a reduced physical footprint while allowing operational parameters (e.g., SRT) to select for desired functional communities. However, as algal technology transitions to full-scale, there is a need to understand the effect of operational and environmental parameters on complex microbial dynamics among mixotrophic microalgae, bacterial groups, and pests (i.e., grazers and pathogens) and to implement robust process controls for stable long-term performance. Here, we examine a full-scale, intensive WRRF utilizing mixed microalgae for tertiary treatment in the US (EcoRecover, Clearas Water Recovery Inc.) during a nine-month monitoring campaign. We investigated the temporal variations in microbial community structure (18S and 16S rRNA genes), which revealed that stable system performance of the EcoRecover system was marked by a low-diversity microalgal community (D INVSIMPSON = 2.01) dominated by Scenedesmus sp. (MRA = 55 %-80 %) that achieved strict nutrient removal (effluent TP < 0.04 mg·L -1 ) and steady biomass concentration (TSS monthly avg . = 400–700 mg·L −1 ). Operational variables including pH, alkalinity, and influent ammonium (NH 4 + ), correlated positively (p < 0.05, method = Spearman) with algal community during stable performance. Further, the use of these parameters as operational controls along with N/P loading and SRT allowed for system recovery following upset events. Importantly, the presence or absence of bacterial nitrification did not directly impact algal system performance and overall nutrient recovery, but partial nitrification (potentially resulting from NO 2 − accumulation) inhibited algal growth and should be considered during long-term operation. The microalgal communities were also adversely affected by zooplankton grazers (ciliates, rotifers) and fungal parasites (Aphelidium), particularly during periods of upset when algal cultures were experiencing culture turnover or stress conditions (e.g., nitrogen limitation, elevated temperature). Altogether, the active management of system operation in order to maintain healthy algal cultures and high biomass productivity can result in significant periods (>4 months) of stable system performance that achieve robust nutrient recovery, even in winter months in northern latitudes (WI, USA).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

HOLISMOKES XIX

We present imaging and spectroscopic observations of supernova SN 2025wny, associated with the lens candidate PS1 J0716+3821. Photometric monitoring from the Lulin and Maidanak observatories confirms multiple point-like images, consistent with SN 2025wny being strongly lensed by two foreground galaxies. Optical spectroscopy of the brightest image with the Nordic Optical Telescope and the University of Hawaii 88-inch Telescope allowed us to determine the redshift to be z SN = 2.008 ± 0.001, based on narrow absorption lines originating in the interstellar medium of the supernova host galaxy. At this redshift, SN 2025wny shows a very high rest-frame UV flux and broad spectral features even weeks after the explosion, which is consistent with superluminous supernovae of Type I. We find a high ejecta temperature and depressed spectral lines compared to other similar objects. We also measured, for the first time, the redshift of the fainter of the two lens galaxies (the ‘perturber’) to be z p = 0.375 ± 0.001, which is fully consistent with the DESI spectroscopic redshift of the main deflector at z d = 0.3754. Thus, SN 2025wny represents the first confirmed galaxy-scale strongly lensed supernova with time delays likely in the range of days to weeks, as judged from the image separations. This makes SN 2025wny suitable for cosmography, offering a promising new system for independent measurements of the Hubble constant. Following a tradition in the field of strongly lensed supernovae, we give SN 2025wny the nickname SN Winny.

Taubenberger, Stefan

Characterization and controllability of radiated power via extrinsic impurity seeding in strongly negative triangularity plasmas in DIII-D

Experiments with extrinsic impurity seeding in strongly negative triangularity shapes in DIII-D achieved radiated power fractions (relative to input power) of up to ≈85% total radiation and ≈55% core radiation in steady-operating conditions. The relationship between core and total radiation was sensitive to impurity species and input power. Attempts to reach higher radiation levels via higher impurity flows resulted in radiative collapse disruptions. Nitrogen, neon, argon, and krypton were tested. Injection was by gas puffing, usually controlled by feeding back real-time estimates for core or total radiating fractions ($P_\textrm{rad} / P_\textrm{input}$). Argon and krypton were controlled by feeding back total radiating fraction, whereas neon was controlled by feeding back the core radiating fraction, due to the lower efficiency of neon as a divertor radiator. Nitrogen flows were pre-programmed. Reasonable total $P_\textrm{rad}$ control target following was achieved with argon or krypton. Control was more challenging with the neon/core radiation configuration, which was more prone to slow response and overshooting of the control target. Poor particle removal contributed to the control challenge: neon particle inventory within the last closed flux surface was roughly constant for up to 1 s (the longest duration tested) after neon injection was halted. However, separate experiments with laser-blow off of non-recycling impurities measured a short impurity confinement time, on the scale of the energy confinement time of ∼100 ms. Modeling with the Aurora impurity transport simulation matched experimental neon density profiles with full recycling (R = 1.0) and weak pumping but predicted rapid decreases in neon inventory if pumping were increased or recycling decreased. This indicates that changes outside the confined plasma (adding a well-placed pump) would improve controllability for all highly recycling species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

High-fidelity remote entanglement of trapped atoms mediated by time-bin photons

Abstract Photonic interconnects between quantum processing nodes are likely the only way to achieve large-scale quantum computers and networks. The bottleneck in such an architecture is the interface between well-isolated quantum memories and flying photons. We establish high-fidelity entanglement between remotely separated trapped atomic qubit memories, mediated by photonic qubits stored in the timing of their pulses. Such time-bin encoding removes sensitivity to polarization errors, enables long-distance quantum communication, and is extensible to quantum memories with more than two states. Using a measurement-based error detection process and suppressing a fundamental source of error due to atomic recoil, we achieve an entanglement fidelity of 97% and show that fundamental limits due to atomic recoil still allow fidelities in excess of 99.9%.

Science & Technology - Other Topics

Characterizing leaf-scale fluorescence with spectral invariants

Sun-induced chlorophyll fluorescence (SIF) is increasingly recognized as a non-destructive probe for tracking terrestrial photosynthesis. Emerging developments in spectral invariants theory provide an innovative and efficient approach for representing SIF radiative transfer processes at the canopy scale. However, modeling leaf-scale fluorescence based on the spectral invariants properties (SIP) remains underexplored. In this study, the spectral invariants theory is employed for the first time to model the leaf-scale total, backward and forward fluorescence (leaf-SIP SIF). The leaf-SIP SIF model separates the leaf-scale radiative transfer process into two distinct components: the wavelength-dependent one associated with leaf biochemical properties, and the wavelength-independent component linked to leaf structural characteristics. The leaf structure-related effects are characterized by two spectrally invariant parameters: the photon recollision probability (p) and the scattering asymmetry parameter (q), which are parameterized using the directly measurable leaf dry matter. Evaluation against field measurements shows that the proposed leaf-SIP SIF model has a good performance, with coefficient of determination (R 2 ) of 0.89, 0.89, 0.90 and root mean squared errors (RMSE) of 1.28, 0.69, 0.74 Wm -2 µm -1 sr -1 , respectively for the total, backward, and forward fluorescence (660–800 nm). The leaf-SIP SIF model with a more concise formulation demonstrates comparable performance with the widely used Fluspect model. Further, the leaf-SIP SIF model provides a simple and efficient approach for simulating leaf-scale fluorescence, with the potential to be integrated into a unified SIP-based model framework for simulating the radiative transfer processes across the soil-leaf-canopy-atmosphere continuum.

59 BASIC BIOLOGICAL SCIENCES

Development of a bench-scale dissolution concept for the direct extraction of nuclear fuel

Current used nuclear fuel reprocessing efforts utilize a hydrometallurgical approach in which the UNF is dissolved in hot nitric acid followed by solvent extraction into an organic solvent to harvest target nuclides. Previous studies have shown that the dissolution and loading process could be combined into a single organic dissolution/extraction step, producing loaded organic in a single step process. This single step process also includes the advantage of selectively targeting key nuclides in the dissolution while leaving undesirable constituents as part of the undissolved solids. This process is referred to hereafter as direct extraction. Ongoing research from multiple national labs has proven the effectiveness of this technique at research scale. Therefore, potential methods to implement direct extraction at both bench and industrial scale have been developed. The key features of potential dissolver system designs were identified via the team at Pacific Northwest National Laboratory and several designs based on industrial counterparts were assessed for feasibility. This report summarizes the advantages and disadvantages of multiple methods, concluding with a path forward to create multiple unique dissolver designs. The first design will be a single stage recirculating eductor mixer. The second design recommendation is a stator rotor static mixing flow loop design. Each dissolver could be utilized separately, simultaneously, or in series to answer questions surrounding reaction kinetics including residence time, provide a proof of concept for targeted extractions of specific nuclides, and inform needs for industrial scale implementation

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Isolating the Vibrational Spectra of the Red Chlorophylls in Photosystem I with Multispectral Two-Dimensional Spectroscopy

Photosystem I (PSI) uses an antenna of chlorophyll (Chl) molecules to create a charge separated state with high quantum efficiency. Understanding the charge separation mechanism is currently hindered by spectral overlap between the antenna and reaction center (RC) Chls and the fact that energy transfer and electron transfer occur with similar time scales. Here, we characterize the antenna excited states by applying two-dimensional electronic (2DES) and two-dimensional electronic-vibrational (2DEV) spectroscopy to PSI complexes with closed RCs. Comparison of the 2DES and 2DEV spectra, which evolve with the same kinetics, enabled characterization of the vibrational modes of the antenna during energy equilibration between spectrally distinct Chls. Through global analysis, we learn how energy transfer between the Bulk and Red Chls presents in the 2DEV spectra and we definitively identify vibrations of the cationic components of the mixed exciton and intermolecular charge transfer states associated with the Red Chls. This work enables future studies of the initial charge separation mechanism of PSI by 2DEV spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic