Search NASA⌕ Search

SEARCH · Search NASA

Results for “Soft computing”

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 1,045 records · Page 58

SU-F-BRF-07: Impact of Different Patient Setup Strategies in Adaptive Radiation Therapy with Simultaneous Integrated Volume-Adapted Boost of NSCLC

Purpose: To evaluate the potential impact of several setup error correction strategies on a proposed image-guided adaptive radiotherapy strategy for locally advanced lung cancer. Methods: Daily 4D cone-beam CT and weekly 4D fan-beam CT images were acquired from 9 lung cancer patients undergoing concurrent chemoradiation therapy. Initial planning CT was deformably registered to daily CBCT images to generate synthetic treatment courses. An adaptive radiation therapy course was simulated using the weekly CT images with replanning twice and a hypofractionated, simultaneous integrated boost to a total dose of 66 Gy to the original PTV and either a 66 Gy (no boost) or 82 Gy (boost) dose to the boost PTV (ITV + 3mm) in 33 fractions with IMRT or VMAT. Lymph nodes (LN) were not boosted (prescribed to 66 Gy in both plans). Synthetic images were rigidly, bony (BN) or tumor and carina (TC), registered to the corresponding plan CT, dose was computed on these from adaptive replans (PLAN) and deformably accumulated back to the original planning CT. Cumulative D98% of CTV of PT (ITV for 82Gy) and LN, and normal tissue dose changes were analyzed. Results: Two patients were removed from the study due to large registration errors. For the remaining 7 patients, D98% for CTV-PT (ITV-PT for 82 Gy) and CTV-LN was within 1 Gy of PLAN for both 66 Gy and 82 Gy plans with both setup techniques. Overall, TC based setup provided better results, especially for LN coverage (p = 0.1 for 66Gy plan and p = 0.2 for 82 Gy plan, comparison of BN and TC), though not significant. Normal tissue dose constraints violated for some patients if constraint was barely achieved in PLAN. Conclusion: The hypofractionated adaptive strategy appears to be deliverable with soft tissue alignment for the evaluated margins and planning parameters. Research was supported by NIH P01CA116602.

60 APPLIED LIFE SCIENCES↗

Efficient Solution of Three-Dimensional Problems of Acoustic and Electromagnetic Scattering by Open Surfaces

We present a computational methodology (a novel Nystrom approach based on use of a non-overlapping patch technique and Chebyshev discretizations) for efficient solution of problems of acoustic and electromagnetic scattering by open surfaces. Our integral equation formulations (1) Incorporate, as ansatz, the singular nature of open-surface integral-equation solutions, and (2) For the Electric Field Integral Equation (EFIE), use analytical regularizes that effectively reduce the number of iterations required by iterative linear-algebra solution based on Krylov-subspace iterative solvers.

sound-soft acoustic scattering↗

System and component design and test of a 10 hp, 18,000 rpm AC dynamometer utilizing a high frequency AC voltage link, part 1

Hard and soft switching test results conducted with one of the samples of first generation MOS-controlled thyristor (MCTs) and similar test results with several different samples of second generation MCT's are reported. A simple chopper circuit is used to investigate the basic switching characteristics of MCT under hard switching and various types of resonant circuits are used to determine soft switching characteristics of MCT under both zero voltage and zero current switching. Next, operation principles of a pulse density modulated converter (PDMC) for three phase (3F) to 3F two-step power conversion via parallel resonant high frequency (HF) AC link are reviewed. The details for the selection of power switches and other power components required for the construction of the power circuit for the second generation 3F to 3F converter system are discussed. The problems encountered in the first generation system are considered. Design and performance of the first generation 3F to 3F power converter system and field oriented induction moter drive based upon a 3 kVA, 20 kHz parallel resonant HF AC link are described. Low harmonic current at the input and output, unity power factor operation of input, and bidirectional flow capability of the system are shown via both computer and experimental results. The work completed on the construction and testing of the second generation converter and field oriented induction motor drive based upon specifications for a 10 hp squirrel cage dynamometer and a 20 kHz parallel resonant HF AC link is discussed. The induction machine is designed to deliver 10 hp or 7.46 kW when operated as an AC-dynamo with power fed back to the source through the converter. Results presented reveal that the proposed power level requires additional energy storage elements to overcome difficulties with a peak link voltage variation problem that limits reaching to the desired power level. The power level test of the second generation converter after the addition of extra energy storage elements to the HF link are described. The importance of the source voltage level to achieve a better current regulation for the source side PDMC is also briefly discussed. The power levels achieved in the motoring mode of operation show that the proposed power levels achieved in the generating mode of operation can also be easily achieved provided that no mechanical speed limitation were present to drive the induction machine at the proposed power level.

Lipo, Thomas A.↗

Chemical and Molecular Descriptors for the Reactivity of Amines with CO{sub 2}

Amine-based solvents are likely to play an important role in CO{sub 2} capture applications in the future, and the identification of amines with superior performance will facilitate their use in CO{sub 2} capture. While some improvements in performance will be achieved through process modifications, modifying the CO{sub 2} capture performance of an amine also implies in part an ability to modify the reactions between the amine and CO{sub 2} through development of new functionalized amines. We present a computational study of trends in the reactions between CO{sub 2} and functionalized amines with a focus on identifying molecular descriptors that determine trends in reactivity. We examine the formation of bicarbonate and carbamate species on three classes of functionalized amines: alkylamines, alkanolamines, and fluorinated alkylamines including primary, secondary and tertiary amines in each class. These functional groups span electron-withdrawing to donating behavior, hydrogen-bonding, extent of functionalization, and proximity effects of the functional groups. Electron withdrawing groups tend to destabilize CO{sub 2} reaction products, whereas electron-donating groups tend to stabilize CO{sub 2} reaction products. Hydrogen bonding stabilizes CO{sub 2} reaction products. Electronic structure descriptors based on electronegativity were found to describe trends in the bicarbonate formation energy. A chemical correlation was observed between the carbamate formation energy and the carbamic acid formation energy. The local softness on the reacting N in the amine was found to partially explain trends carbamic acid formation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasi-deuteron model at low renormalization group resolution

The quasi-deuteron model introduced by Levinger is used to explain cross sections for knocking out high-momentum protons in photoabsorption on nuclei. This is within a framework we characterize as exhibiting high renormalization group (RG) resolution. Assuming a one-body reaction operator, the nuclear wave function must include two-body short-range correlations (SRCs) with deuteronlike quantum numbers. In Phys. Rev. C 104, 034311 (2021), we showed that SRC physics can be naturally accounted for at low RG resolution. We describe the quasi-deuteron model at low RG resolution and determine the Levinger constant, which is proportional to the ratio of nuclear photoabsorption to that for photodisintegration of a deuteron. We extract the Levinger constant based on the ratio of momentum distributions at high relative momentum. We compute momentum distributions evolved under similarity RG (SRG) transformations where the SRC physics is shifted into the operator as a universal two-body term. The short-range nature of this operator motivates using local-density approximations with uncorrelated wave functions in evaluating nuclear matrix elements, which greatly simplifies the analysis. The operator must be consistently matched to the RG scale and scheme of the interaction for a reliable extraction. We apply SRG transformations to different nucleon-nucleon (NN) interactions and use the deuteron wave functions and Weinberg eigenvalues to determine approximate matching scales. We predict the Levinger constant for several NN interactions and a wide range of nuclei comparing to experimental extractions. The predictions at low RG resolution are in good agreement with experiment when starting with a hard NN interaction and the initial operator. Similar agreement is found using soft NN interactions when the additional two-body operator induced by evolution from hard to soft is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neural Network and Regression Soft Model Extended for PAX-300 Aircraft Engine

In fiscal year 2001, the neural network and regression capabilities of NASA Glenn Research Center's COMETBOARDS design optimization testbed were extended to generate approximate models for the PAX-300 aircraft engine. The analytical model of the engine is defined through nine variables: the fan efficiency factor, the low pressure of the compressor, the high pressure of the compressor, the high pressure of the turbine, the low pressure of the turbine, the operating pressure, and three critical temperatures (T(sub 4), T(sub vane), and T(sub metal)). Numerical Propulsion System Simulation (NPSS) calculations of the specific fuel consumption (TSFC), as a function of the variables can become time consuming, and numerical instabilities can occur during these design calculations. "Soft" models can alleviate both deficiencies. These approximate models are generated from a set of high-fidelity input-output pairs obtained from the NPSS code and a design of the experiment strategy. A neural network and a regression model with 45 weight factors were trained for the input/output pairs. Then, the trained models were validated through a comparison with the original NPSS code. Comparisons of TSFC versus the operating pressure and of TSFC versus the three temperatures (T(sub 4), T(sub vane), and T(sub metal)) are depicted in the figures. The overall performance was satisfactory for both the regression and the neural network model. The regression model required fewer calculations than the neural network model, and it produced marginally superior results. Training the approximate methods is time consuming. Once trained, the approximate methods generated the solution with only a trivial computational effort, reducing the solution time from hours to less than a minute.

Patnaik, Surya N.↗

A Tunable, Software-based DRAM Error Detection and Correction Library for HPC

Proposed exascale systems will present a number of considerable resiliency challenges. In particular, DRAM soft-errors, or bit-flips, are expected to greatly increase due to the increased memory density of these systems. Current hardware-based fault-tolerance methods will be unsuitable for addressing the expected soft error frequency rate. As a result, additional software will be needed to address this challenge. In this paper we introduce LIBSDC, a tunable, transparent silent data corruption detection and correction library for HPC applications. LIBSDC provides comprehensive SDC protection for program memory by implementing on-demand page integrity verification. Experimental benchmarks with Mantevo HPCCG show that once tuned, LIBSDC is able to achieve SDC protection with 50\% overhead of resources, less than the 100\% needed for double modular redundancy.

99 GENERAL AND MISCELLANEOUS↗

An Ionomeric Renewable Thermoplastic from Lignin-Reinforced Rubber

An ionomeric, leathery thermoplastic with high mechanical strength is prepared by a new thermal processing method from a soft, melt-processable rubber. Compositions made by incorporation of equal-mass lignin, a renewable oligomeric feedstock, in an acrylonitrile-butadiene rubber often yield weak rubbers with large lignin domains (1–2 µm). The addition of zinc chloride (ZnCl2) in such a composition based on sinapyl alcohol-rich lignin during a solvent-free synthesis induces a strong interfacial crosslinking between lignin and rubber phases. This compositional modification results in finely interspersed lignin domains (<100 nm) that essentially reinforce the rubbery matrix with a 10–22 °C rise in the glassy-to-rubbery transition temperature. The ion-modified polymer blends also show improved materials properties, like a 100% increase in ultimate tensile strength and an order of magnitude rise in Young's modulus. Coarse-grained molecular dynamics (MD) simulations verify the morphology and dynamics of the ionomeric material. The computed result also confirms that the ionomers have glassy characteristics.

36 MATERIALS SCIENCE↗

3D exemplar-based random walks for tooth segmentation from cone-beam computed tomography images

Purpose: Tooth segmentation is an essential step in acquiring patient-specific dental geometries from cone-beam computed tomography (CBCT) images. Tooth segmentation from CBCT images is still a challenging task considering the comparatively low image quality caused by the limited radiation dose, as well as structural ambiguities from intercuspation and nearby alveolar bones. The goal of this paper is to present and discuss the latest accomplishments in semisupervised tooth segmentation with adaptive 3D shape constraints. Methods: The authors propose a 3D exemplar-based random walk method of tooth segmentation from CBCT images. The proposed method integrates semisupervised label propagation and regularization by 3D exemplar registration. To begin with, the pure random walk method is to get an initial segmentation of the teeth, which tends to be erroneous because of the structural ambiguity of CBCT images. And then, as an iterative refinement, the authors conduct a regularization by using 3D exemplar registration, as well as label propagation by random walks with soft constraints, to improve the tooth segmentation. In the first stage of the iteration, 3D exemplars with well-defined topologies are adapted to fit the tooth contours, which are obtained from the random walks based segmentation. The soft constraints on voxel labeling are defined by shape-based foreground dentine probability acquired by the exemplar registration, as well as the appearance-based probability from a support vector machine (SVM) classifier. In the second stage, the labels of the volume-of-interest (VOI) are updated by the random walks with soft constraints. The two stages are optimized iteratively. Instead of the one-shot label propagation in the VOI, an iterative refinement process can achieve a reliable tooth segmentation by virtue of exemplar-based random walks with adaptive soft constraints. Results: The proposed method was applied for tooth segmentation of twenty clinically captured CBCT images. Three metrics, including the Dice similarity coefficient (DSC), the Jaccard similarity coefficient (JSC), and the mean surface deviation (MSD), were used to quantitatively analyze the segmentation of anterior teeth including incisors and canines, premolars, and molars. The segmentation of the anterior teeth achieved a DSC up to 98%, a JSC of 97%, and an MSD of 0.11 mm compared with manual segmentation. For the premolars, the average values of DSC, JSC, and MSD were 98%, 96%, and 0.12 mm, respectively. The proposed method yielded a DSC of 95%, a JSC of 89%, and an MSD of 0.26 mm for molars. Aside from the interactive definition of label priors by the user, automatic tooth segmentation can be achieved in an average of 1.18 min. Conclusions: The proposed technique enables an efficient and reliable tooth segmentation from CBCT images. This study makes it clinically practical to segment teeth from CBCT images, thus facilitating pre- and interoperative uses of dental morphologies in maxillofacial and orthodontic treatments.

60 APPLIED LIFE SCIENCES↗

The X-ray spectra of the Vela and Puppis supernova remnants and the shock-wave model of supernova remnants

The structure and spectrum of the Vela and Puppis supernova remnants (SNRs) were observed in soft (0.1 to 1.5 keV) X rays from a rocket. The spectral data of both objects below 1.5 keV are consistent with thermal radiation from a hot plasma and are not consistent with a simple power-law function. Fitting the data to a Tucker-Koren model results in values of T = (4.3 plus or minus 0.3) x 1,000,000 K for Vela and T = (7 plus or minus 1) x 1,000,000 K for Puppis. The present data, earlier data from the Cygnus Loop, and results on Cas A and Tycho reported by others are considered within the context of a simple model in which X rays are thermally produced as a result of a shock wave expanding into the interstellar medium. X-ray data on the temperature, intrinsic diameter, and intrinsic luminosity of the five SNRs are used to compute the age, energy release, and initial interstellar density for each one. From measurements of X-ray absorption along the line of sight, and also from shock-wave model interpretations of temperature-diameter-luminosity relations, it is concluded that the average density of interstellar material is about 0.1 per cu cm between the sun and the nearest X-ray emitting SNR, Vela.

Gorenstein, P.↗

X-ray refelection from photoionized media in active galactic nuclei

We calculate the spectrum of X-ray radiation and reprocessed by a partly ionized optically thick medium in an active galactic nucleus. We self-consistently calculate the ionization balance and thermal balance in the medium along with the distribution of X-ray intensity with optical depth. In addition to absorption or scattering of the incident X-rays, we also compute the spectrum of X-rays emitted by the material, including lines, edges, and bremsstrahlung. The albedo of the medium depends primarily on the X-ray ionization parameter (ratio of incident flux to gas density, zeta(sub Chi), and secondarily on the UV flux generated by dissipation inside the disk; we locate the critical range of zeta(sub Chi) over which the albedo increases from small to nearly unity. While the continuum reflection is very weak below 10 keV when zeta(sub Chi) is small, significnat fluxes are emitted in atomic lines and edges in this energy range. In the limit of large zeta(sub Chi), the albedo below 10 keV increases, but reflection in this band is never gray: some photoelectric absorption remains up to rather large values of zeta(sub Chi), while at still higher values, inverse Compton scattering amplifies the soft X-ray flux. These features are sufficiently sharp that current and near-future X-ray experiments should permit diagnostic measures of zeta(sub Chi).

Zycki, Piotr T.↗

Energy levels, lifetimes and radiative data of Ba XXVI

Highlights: • Lowest 110 levels for Ba XXVI are calculated. • Radiative data for E1, E2, M1 and M2 transitions within lowest 110 levels are presented. • Lifetimes for lowest 110 fine structure levels are provided. We report an extensive and an elaborate theoretical study of atomic data for Ba XXVI by considering Singlet, Doublet and Triplet (SDT) electron excitations within N-shell and single excitations from N-shell to O-shell. We have calculated energy levels and lifetimes for lowest 110 fine structure levels by using Multi-configuration Dirac–Fock method (MCDF). We have also considered Quantum Electrodynamics (QED) and Breit corrections in our calculations. We have presented the radiative data for electric and magnetic dipole (E1, M1) and quadrupole (E2, M2) transitions among lowest 110 levels. We have made comparisons of our calculated excitation energies and EUV (Extreme Ultraviolet) transition wavelengths with experimentally observed energy levels and wavelengths and achieved good agreement. We have also computed energy levels by performing similar relativistic distorted wave calculations using Flexible Atomic Code (FAC). Additionally, we have provided new atomic data for Ba XXVI which are not published elsewhere in the literature. We believe that our results may be beneficial in fusion plasma research and astrophysical investigations and applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Freely cooling granular gases with short-ranged attractive potentials

We treat the case of an undriven gas of inelastic hard-spheres with short-ranged attractive potentials via an extension of the pseudo-Liouville operator formalism. New evolution equations for the granular temperature and coordination number are obtained. The granular temperature exhibits deviation from both Haff’s law and the case of long-ranged potentials. We verify this departure using soft-sphere discrete element method simulations. Excellent agreement is found for the duration of the simulation even beyond where exclusively binary collisions are expected. Simulations show the emergence of strong spatial-velocity correlations on the length scale of the last peak in the pair-correlation function but do not show strong correlations beyond this length scale. We argue that molecular chaos may remain an adequate approximation if the system is modelled as a Smoluchowski type equation with aggregation and break-up processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Class II algorithm for charged-particle transport simulation - Paper 142

Simulation of the transport of high-energy charged particles is difficult because of the large number of interactions experienced by these particles before being brought to rest. Possible strategies to deal with this difficulty were analyzed by Berger. Most modern Monte Carlo codes use a mixture of complete-grouping (class I) and mixed (class II) schemes. In class I simulations, the global effect of all the interactions that occur along a trajectory segment of a given length s is described by means of precalculated distributions derived from multiple-scattering theories. In class II simulations, hard (or catastrophic) interactions (with energy loss W or polar scattering angle {theta} larger than given cutoff values) are sampled individually, while the effect of soft interactions (with energy losses and scattering angles smaller than the respective cutoffs) is described using a multiple-scattering approach. Although Berger effectively employed class II simulation only for energy-loss events, he did recognize the advantages of class II over class I also for multiple elastic scattering. On the one hand, because hard interactions are described exactly according to the adopted differential cross sections (DCSs), class II schemes are more accurate. On the other hand, multiple-scattering approximations, which are still required to describe the cumulative effect of soft interactions, are simpler and more accurate when individual interactions are limited to small energy transfers and angular deflections. For instance, for path lengths that are not too small, but still much less than the CSDA range, the energy-loss distribution due to multiple soft events is approximately Gaussian, by virtue of the central limit theorem. Similarly, the angular distribution caused by multiple small-angle deflections is also approximately Gaussian. However, class II schemes are more difficult to implement because of their dependence on specificities of the DCSs. In this article we briefly present the essentials of the class II algorithm for charged particles implemented in the last version of the penelope code, and provide evidence of its accuracy and stability with the simulation parameters. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

WE-G-BRD-03: Development of a Real-Time Optical Tracking Goggle System (OTGS) for Intracranial Stereotactic Radiotherapy

Purpose: Optical tracking systems (OTS) are an acceptable alternative to frame-based stereotactic radiotherapy (SRT). However, current surface-based OTS lack the ability to target exclusively rigid/bony anatomical features. We propose a novel marker-based optical tracking goggle system (OTGS) that provides real-time guidance based on the nose/facial bony anatomy. This ongoing study involves the development and characterization of the OTGS for clinical implementation in intracranial stereotactic radiotherapy. Methods: The OTGS consists of eye goggles, a custom thermoplastic nosepiece, and 6 infrared markers pre-attached to the goggles. A phantom and four healthy volunteers were used to evaluate the calibration/registration accuracy, intrafraction accuracy, interfraction reproducibility, and end-to-end accuracy of the OTGS. The performance of the OTGS was compared with that of the frameless SonArray system and cone-beam computed tomography (CBCT) for volunteer and phantom cases, respectively. The performance of the OTGS with commercial immobilization devices and under treatment conditions (i.e., couch rotation and translation range) was also evaluated. Results: The difference in the calibration/registration accuracy of 24 translations or rotation combinations between CBCT and in-house OTS software was within 0.5 mm/0.4°. The mean intrafraction and interfraction accuracy among the volunteers was 0.004+/−0.4mm with −0.09+/−0.5° (n=6,170) and −0.26+/−0.8mm with 0.15+/0.8° (n=11), respectively. The difference in end-to-end accuracy between the OTGS and CBCT was within 1.3 mm/1.1°. The predetermined marker pattern (1) minimized marker occlusions, (2) allowed for continuous tracking for couch angles +/− 90°, (3) and eliminated individual marker misplacement. The device was feasible with open and half masks for immobilization. Conclusion: Bony anatomical localization eliminated potential errors due to facial hair changes and/or soft tissue deformation. The OTGS offers a workflow-friendly, patient-friendly solution for intracranial SRT, while being comparable to other real-time options. The minimum rotation uncertainty of the OTGS can be combined with CBCT to ensure maximum accuracy for high-precision SRT.

60 APPLIED LIFE SCIENCES↗

Experimental and computational assessment of U Si N ternary phases

Uranium nitride-silicide composites are being considered as a high-density and high thermal conductivity fuel option for light water reactors. During development, chemical interactions were observed near the silicide melting point which resulted in formation of an unknown U-Si-N ternary phase. Here in the present work, U-Si-N composite samples were produced by arc-melting U 3 Si 2 under an argon-nitrogen atmosphere to form the ternary phase. The resulting samples were characterized by SEM/EDS-EPMA and XRD, and demonstrated an equilibrium between U 3 Si 2 , UN, USi and a U-Si-N phase with a distinct crystallographic structure. Rietveld refinement of the ternary structure was performed, considering the ternary structures existent in the analogue U-Si-C system, and a good fit was obtained for the hexagonal U 20 Si 16 N 3 phase. DFT + U calculations were performed in parallel to evaluate the thermodynamic and dynamic stability of the ternaries U 20 Si 16 N 3 and U 3 Si 2 N 2 . The calculated enthalpy of formation and phonon dispersion support the existence of stable U 20 Si 16 N 3 and U 3 Si 2 N 2 , although some soft modes in the U 20 Si 16 N 3 phase phonons are observed. The results presented here thus demonstrate the occurrence of at least one ternary phase in the U-Si-N system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mechanical specific energy versus depth of cut in rock cutting and drilling

The relationship between Mechanical Specific Energy (MSE) and the Rate of Penetration (ROP), or equivalently the depth of cut per revolution, provides an important measure for strategizing a drilling operation. This study explores how MSE evolves with depth of cut, and presents a concerted effort that encompasses analytical, computational and experimental approaches. A simple model for the relationship between MSE and cutting depth is first derived with consideration of the wear progression of a circular cutter. This is an extension of Detournay and Defourny's phenomenological cutting model. Wear is modeled as a flat contact area at the bottom of a cutter referred to as a wear flat, and that wear flat in the past is often considered to be fixed during cutting. But during a drilling operation by a full bit that consists of multiple circular cutters, the wear flat length may increase because of various wear mechanisms involved. The wear progression of cutters generally results in reduced efficiency with either increased MSE or decreased ROP. Also, an accurate estimate of removed rock volume is found important for the evaluation of MSE. The derived model is compared with experiment results from a single circular cutter, for cutting soft rock under ambient pressure with actual depth measured through a micrometer, and for cutting high strength rock under high pressure with actual cutting area measured by a confocal microscope. Lastly, the model is employed to interpret the evolution of MSE with depth of cut for a full drilling bit under confining pressure. The general form of equation of the developed model is found to describe well the experiment data and can be applied to interpret the drilling data for a full bit.

42 ENGINEERING↗

Edwards statistical mechanics for jammed granular matter

In 1989, Sir Sam Edwards made the visionary proposition to treat jammed granular materials using a volume ensemble of equiprobable jammed states in analogy to thermal equilibrium statistical mechanics, despite their inherent athermal features. Since then, the statistical mechanics approach for jammed matter—one of the very few generalizations of Gibbs-Boltzmann statistical mechanics to out-of-equilibrium matter—has garnered an extraordinary amount of attention by both theorists and experimentalists. Its importance stems from the fact that jammed states of matter are ubiquitous in nature appearing in a broad range of granular and soft materials such as colloids, emulsions, glasses, and biomatter. Indeed, despite being one of the simplest states of matter—primarily governed by the steric interactions between the constitutive particles—a theoretical understanding based on first principles has proved exceedingly challenging. Here a systematic approach to jammed matter based on the Edwards statistical mechanical ensemble is reviewed. The construction of microcanonical and canonical ensembles based on the volume function, which replaces the Hamiltonian in jammed systems, is discussed. The importance of approximation schemes at various levels is emphasized leading to quantitative predictions for ensemble averaged quantities such as packing fractions and contact force distributions. An overview of the phenomenology of jammed states and experiments, simulations, and theoretical models scrutinizing the strong assumptions underlying Edwards approach is given including recent results suggesting the validity of Edwards ergodic hypothesis for jammed states. A theoretical framework for packings whose constitutive particles range from spherical to nonspherical shapes such as dimers, polymers, ellipsoids, spherocylinders or tetrahedra, hard and soft, frictional, frictionless and adhesive, monodisperse, and polydisperse particles in any dimensions is discussed providing insight into a unifying phase diagram for all jammed matter. Furthermore, the connection between the Edwards ensemble of metastable jammed states and metastability in spin glasses is established. This highlights the fact that the packing problem can be understood as a constraint satisfaction problem for excluded volume and force and torque balance leading to a unifying framework between the Edwards ensemble of equiprobable jammed states and out-of-equilibrium spin glasses.

97 MATHEMATICS AND COMPUTING↗