Numerical Analysis of Fluidized Bed Chemical Vapor Deposition of Pyrolytic Carbon for Tristructural-Isotropic (TRISO) Fuel Particles
Explore the source record for details and available documents.
Dal Forno Chuahy, Flavio↗
SEARCH · Search NASA
Results for “ISOTROPISM”
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.
3D orthorhombic earth model effects on seismic source characterization
Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).
58 GEOSCIENCES↗
R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method
The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Controlling Factors of Chorus Spectral Gaps
Abstract The present study compares a single‐band chorus wave against a banded chorus wave observed by Van Allen Probes at adjacent times, and demonstrates that the single‐band chorus wave is associated with an anisotropic electron population over a broad energy range, while the banded chorus wave is accompanied by an electron phase space density plateau and an electron anisotropy reduction around Landau resonant energies. We further compare banded chorus waves with different spectral gap widths, and show that a wider spectral gap is associated with electron isotropization extending to higher energies with respect to the equatorial Landau resonant energy. We suggest that early generated chorus waves isotropize electrons via Landau resonant acceleration, and the waves that propagate to higher latitudes isotropize electrons at higher energies. The isotropization extending to higher energies leads to a larger spectral gap of new chorus waves after electrons bounce back to the equator.
Li, Jinxing↗
Inferring thermal ion temperature and residual kinetic energy from nuclear measurements in inertial confinement fusion implosions
In inertial con finement fusion implosion experiments, the presence of residual anisotropic fluid motion within the stagnating hot spot leads to signifi cant variations in ion-temperature measurements using neutron time- of-flight detectors along different lines of sight. The minimum ion-temperature measurement is typically used as representative of the thermal temperature. In the presence of isotropic flows, however, even the minimum DT neutron-inferred ion temperature can be well above the plasma thermal temperature. Using both DD and DT neutron-inferred ion-temperature measurements, we show that it is possible to determine the contribution of isotropic flows and infer the DT burn-averaged thermal ion temperature. The contribution of large isotropic flows on driving the ratio of DD to DT neutron-inferred ion temperatures well below unity and approaching the lower bound of 0.8 is demonstrated in multimode simulations. The minimum DD neutron- inferred ion temperature is determined from the velocity variance analysis, accounting for the presence of isotropic flows. Being close to the DT burn-averaged thermal ion temperature, the inferred DD minimum ion temperatures demonstrate a strong correlation with the experimental yields in the OMEGA implosion database. An analytical expression is also derived to explain the effect of mode ℓ = 1 ion-temperature measurement asymmetry on yield degradations caused by the anisotropic flows.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Fusion yield of plasma with velocity-space anisotropy at constant energy
Velocity-space anisotropy can significantly modify fusion reactivity. The nature and magnitude of this modification depends on the plasma temperature, as well as the details of how the anisotropy is introduced. For plasmas that are sufficiently cold compared to the peak of the fusion cross section, anisotropic distributions tend to have higher yields than isotropic distributions with the same thermal energy. At higher temperatures, it is instead isotropic distributions that have the highest yields. However, the details of this behavior depend on exactly how the distribution differs from an isotropic Maxwellian. Therefore, this paper describes the effects of anisotropy on fusion yield for the class of anisotropic distribution functions with the same energy distribution as a 3D isotropic Maxwellian and compares those results with the yields from bi-Maxwellian distributions. In many cases, especially for plasmas somewhat below reactor-regime temperatures, the effects of anisotropy can be substantial.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model
We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Flows into de Sitter space from anisotropic initial conditions: An effective field theory approach
For decades, physicists have analyzed various versions of a “cosmic no-hair” conjecture, to understand under what conditions a spacetime that is initially spatially anisotropic and/or inhomogeneous will flow into an isotropic and homogeneous state. Wald’s theorem, in particular, established that homogeneous but anisotropic spacetimes, if filled with a positive cosmological constant plus additional matter sources that satisfy specific energy conditions, will necessarily flow toward an (isotropic) de Sitter state at late times. Here, in this paper, we study the flow of homogeneous but anisotropic spacetimes toward isotropic states under conditions more general than those to which Wald’s theorem applies. We construct an effective field theory (EFT) treatment for generic “single-clock” systems in anisotropic spacetimes—which are not limited to realizations compatible with scalar-field constructions—and identify fixed points in the resulting phase space. We identify regions of this phase space that flow to isotropic fixed points—including a de Sitter fixed point—even in the absence of a bare cosmological constant, and for matter sources that do not obey the energy conditions required for Wald’s theorem. Such flows into de Sitter reveal the emergence of an effective cosmological constant.
79 ASTRONOMY AND ASTROPHYSICS↗
Coal to Carbon Fiber (C2CF) Continuous Processing for High Value Composites (Final Report)
Coal tar is a condensed and recovered by-product of the coking of metallurgical coal for steel production. The heaviest fraction of distilled coal tar is an isotropic pitch largely used as a binder in the manufacturing of carbonaceous electrodes for primary aluminum smelting and in electric arc furnaces. Coal tar pitch offers high carbon yield upon carbonization. In this project, a process to convert the domestically sourced isotropic coal tar pitch, containing very low particulates (QI = 0.32 wt.%), to form flow-domain mesophase pitch amenable for melt spinning into precursor (green) fibers for carbon fiber, was developed and optimized. The final reproducible processing developed is reviewed in this report, along with several characterizations of the mesophase pitch. A final definition of the characteristics of a ‘spinnable’ mesophase pitch is presented. With this mesophase pitch, reproducible and stable multifilament melt spinning was developed, producing green fiber tows with filament diameters of approximately 20 m. The multifilament melt spinning was the most challenging aspect of the project and required the most effort. The best practices learned from this project for melt spinning are reviewed in this report. Once spun, the green fibers were oxidatively stabilized, carbonized and graphitized under inert gas atmosphere to form the final carbon fibers. Given the relatively high softening point of the mesophase pitch, no issues of interfilament fusion were observed during batch oxidation, and subsequent batch carbonization and graphitization went smoothly in all cases. An ~ 80 wt.% conversion of the green fiber mass to final carbon fiber was achieved. After graphitization, the carbon fibers showed high tensile moduli (most were ~ 600 GPa, or 87 Msi) consistent with commercially-available high-performance pitch-based carbon fiber. However, tensile strength and strain to failure were comparatively low. Further work to reduce defects in and on the fiber surfaces would increase these properties. SEM imaging of the graphitic fiber textures is presented herein. Rudimentary composites were fabricated from the carbon fibers and characterized showing similar modulus to baseline composites fabricated with commercial carbon fiber. Finally, a basic economic analysis was done showing the potential to increase the value of the isotropic coal tar pitch by up to 13.6 to 136 times based on a carbon fiber value of $\$$5/lb to $\$$50/lb, respectively. Moreover, the site case study suggests that the coal tar from the single integrated steel mill could supply production of up to 16 kt/yr of carbon fiber. Finally, a technological gap analysis was done which shed light on remaining technical challenges. These challenges included: recovery and utilization of condensates from the mesophase pitch processing, further advancing and increasing stability of the multifilament melt spinning processing, optimization of the oxidation processing, defect reduction for increased carbon fiber strength, and the development of a weaving process towards carbon fiber fabrics. To maximize the coal value chain, the primary objectives of this project were to (a) develop and scale efficient processing technology for producing melt-spinnable mesophase pitch from isotropic coal tar pitch, (b) clarify and simplify tedious continuous fiber processing technologies (particularly multifilament melt spinning of mesophase pitch) towards the efficient production of high performance carbon fiber, and (c) demonstrate and characterize representative composite parts derived from the final carbon fiber. Immense progress was made on all 3 objectives and is detailed in this report.
01 COAL, LIGNITE, AND PEAT↗
The Linear Point Standard Ruler with DESI DR1 and DR2 Data
The linear point, a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler. Compared to the peak in the correlation function, it is more robust to late-time nonlinear effects at the percent level. In light of improved simulations and high quality data, we revisit the robustness of the linear point and use it as an alternative to template-based fitting approaches typically used in BAO analyses. We present the linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the DESI survey. We convert the linear point into a dimensionless parameter $α_{iso,LP}$, defined as the ratio of the linear point in the fiducial cosmology and the observed value, analogous to the isotropic BAO scaling parameter $α_{iso}$ used in previous BAO measurements. Using the 2nd generation of AbacusSummit mock catalogs, we find that linear point measurements are more precise when calculated in the post-reconstruction regime with 15-60% smaller uncertainties than those pre-reconstruction. We find a systematic shift in the linear point measurements compared against the isotropic BAO measurements in mocks; we attribute this to the isotropic damping parameter responsible for smearing the linear point in the nonlinear regime. We propose a sample-dependent correction that mitigates the impact of late-time nonlinear effects. While this introduces a cosmology dependence in an otherwise model-independent measurement, this is necessary given the sub-percent precision dictated by current cosmological surveys. Comparing $α_{iso,LP}$ with isotropic BAO measurements made on the DESI DR1 and DR2 galaxy samples, we find excellent agreement after applying this correction, particularly post-reconstruction. We discuss future scope regarding cosmological inference with linear point measurements.
Uberoi, N. [Yale U.] (ORCID:0000000275179629)↗
Measurements of multijet event isotropies using optimal transport with the ATLAS detector
A measurement of novel event shapes quantifying the isotropy of collider events is performed in 140 fb -1 of proton-proton collisions with $\sqrt{s}$ = 13 TeV centre-of-mass energy recorded with the ATLAS detector at CERN’s Large Hadron Collider. These event shapes are defined as the Wasserstein distance between collider events and isotropic reference geometries. This distance is evaluated by solving optimal transport problems, using the ‘Energy-Mover’s Distance’. Isotropic references with cylindrical and circular symmetries are studied, to probe the symmetries of interest at hadron colliders. The novel event-shape observables defined in this way are infrared- and collinear-safe, have improved dynamic range and have greater sensitivity to isotropic radiation patterns than other event shapes. The measured event-shape variables are corrected for detector effects, and presented in inclusive bins of jet multiplicity and the scalar sum of the two leading jets’ transverse momenta. The measured distributions are provided as inputs to future Monte Carlo tuning campaigns and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Determining elastic anisotropy of textured polycrystals using resonant ultrasound spectroscopy
Abstract Polycrystalline materials can have complex anisotropic properties depending on their crystallographic texture and crystal structure. In this study, we use resonant ultrasound spectroscopy (RUS) to nondestructively quantify the elastic anisotropy in extruded aluminum alloy 1100-O, an inherently low-anisotropy material. Further, we show that RUS can be used to indirectly provide a description of the material’s texture, which in the present case is found to be transversely isotropic. By determining the entire elastic tensor, we can identify the level and orientation of the anisotropy originated during extrusion. The relative anisotropy of the compressive (c 11 /c 33 ) and shear (c 44 /c 66 ) elastic constants is 1.5% ± 0.5% and 5.7% ± 0.5%, respectively, where the elastic constants (five independent elastic constants for transversely isotropic) are those associated with the extrusion axis that defines the symmetry of the texture. These results indicate that the texture is expected to have transversely isotropic symmetry. This finding is confirmed by two additional approaches. First, we confirm elastic constants and the degree of elastic anisotropy by direct sound velocity measurements using ultrasonic pulse echo. Second, neutron diffraction (ND) data confirm the symmetry of the bulk texture consistent with extrusion-induced anisotropy, and polycrystal elasticity simulations using the elastic self-consistent model with input from ND textures and aluminum single-crystal elastic constants render similar levels of polycrystal elastic anisotropy to those measured by RUS. We demonstrate the ability of RUS to detect texture-induced anisotropy in inherently low-anisotropy materials. Therefore, as many other common materials have intrinsically higher elastic anisotropy, this technique should be applicable for similar levels of texture, providing an efficient general diagnostic and characterization tool.
36 MATERIALS SCIENCE↗
Beyond curvature-driven grain growth: Insights from fully anisotropic Monte Carlo Potts simulations
Grain boundary (GB) motion away from the center of curvature, termed anti-curvature behavior, has recently been observed in 3D experiments but is not predicted by classical grain growth theory. In this study, we investigate this behavior using a novel, fully anisotropic Monte Carlo Potts (MCP) model that incorporates both misorientation and inclination dependencies of GB energy. We perform 3D grain growth simulations with isotropic and anisotropic GB energies to explore the relationship between GB velocity and curvature. Contrary to the classical relation that velocity is a product of reduced mobility and mean curvature, we observe no consistent correlation between velocity and mean curvature for individual GBs, even under isotropic conditions, though there is correlation between the average velocity and curvature for some cases. The 3D simulations exhibit frequent anti-curvature motion regardless of GB energy anisotropy including with isotropic GBs, though larger curvatures occur with anisotropic GB energy functions that promote low-energy GBs. In 2D simulations, anti-curvature behavior only occurs with the anisotropic functions that favor low-energy GBs. This difference between 2D and 3D results suggests that anti-curvature behavior results in part from the increased freedom of motion intrinsic in 3D GB networks. Furthermore, our results support recent experimental observations that demonstrate that simple curvature-driven models are insufficient for describing GB migration in polycrystals.
Anti-curvature↗
The oscillatory anisotropy in the spatially flat cosmological models
Highlights: • Relations among asymptotic anisotropy and matter content. • A new scheme for obtaining physically motivated exact solutions. • The process of isotropization/anisotropization in cosmology. • The conversion among the anisotropy and the energy density into one another. In this manuscript, we investigate the oscillatory behaviour of the anisotropy in the diagonal Bianchi-I spacetimes. Our starting point is a simplification of Einstein’s equations using only observable or physical variables. As a consequence, we are able to: (a) Prove general results concerning the existence of oscillations of the anisotropy in the primordial and the late-time universe. For instance, in the expanding scenario, we show that a past weakly mixmaster behaviour (oscillations as we approach the Kasner solutions) might appear even with no violation of the usual energy conditions, while in the future, the pulsation (oscillations around isotropic solutions) seems to be most favoured; (b) Determine a large scheme for deriving classes of physically motivated exact solutions, and we give some (including the general barotropic perfect fluid and the magnetic one); (c) Understand the physical conditions for the occurrence of the isotropization or anisotropization during the cosmological evolution; (d) Understand how anisotropy and energy density are converted one into another. In particular, we call attention to the presence of a residue in the energy density in a late-time isotropic universe coming from its past anisotropic behaviour.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing
Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.
36 MATERIALS SCIENCE↗
A note on the instability and pattern formation of shrinkage cracks in viscoplastic soils
In this work we present a theoretical study on the conditions for the onset of cracks, as well as the corresponding pattern formation, in saturated viscoplastic soils under isotropic loading (extension). The type of stress applied is left unspecified, to cover a variety of loadings including shrinkage due to desiccation, isotropic thermal expansion, mechanical loading and so forth. By treating the saturated soil as rigid viscoplastic, we obtain a 2D extension of the Cnoidal Waves equations (Veveakis and Regenauer-Lieb, 2015). By numerically solving the corresponding boundary value problem, we retrieve conditions for the onset of cracking instability in 2D loading, and identify the characteristic spacing between cracks to be a length scale combining all the hydro-mechanical parameters of the problem. Finally, we show that in a rectangular slab of clay under isotropic extension, patterns of triangular, rectangular and hexagonal cracks can tessellate the domain, with the hexagonal pattern being the energetically favored, as it minimizes the free energy of the system.
58 GEOSCIENCES↗
Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository
To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Machine Learning Full NMR Chemical Shift Tensors of Silicon Oxides with Equivariant Graph Neural Networks
The nuclear magnetic resonance (NMR) chemical shift tensor is a highly sensitive probe of the electronic structure of an atom and furthermore its local structure. Recently, machine learning has been applied to NMR in the prediction of isotropic chemical shifts from a structure. Current machine learning models, however, often ignore the full chemical shift tensor for the easier-to-predict isotropic chemical shift, effectively ignoring a multitude of structural information available in the NMR chemical shift tensor. Here we use an equivariant graph neural network (GNN) to predict full 29 Si chemical shift tensors in silicate materials. The equivariant GNN model predicts full tensors to a mean absolute error of 1.05 ppm and is able to accurately determine the magnitude, anisotropy, and tensor orientation in a diverse set of silicon oxide local structures. When compared with other models, the equivariant GNN model outperforms the state-of-the-art machine learning models by 53%. The equivariant GNN model also outperforms historic analytical models by 57% for isotropic chemical shift and 91% for anisotropy. The software is available as a simple-to-use open-source repository, allowing similar models to be created and trained with ease.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗