Search NASA⌕ Search

SEARCH · Search NASA

Results for “NEUTRON DISTRIBUTION”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Theoretical description of fission yields: Toward a fast and efficient global model

A recently proposed approach to global modeling of fission fragment distributions is extended to account for odd-even staggering in charge yields and for neutron evaporation. We show that the proposed hybrid model can reproduce experimental mass and charge fragment yields, including the odd-even staggering for a wide range of fissioning nuclei. Experimental isotopic yields can be described within a simple neutron evaporation scheme. Here, we also explore fission fragment distributions of exotic neutron-rich and superheavy systems and compare our predictions with other state-of-the art global calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Sample spinning to mitigate polarization artifact and interstitial-vacancy imbalance in ion-beam irradiation

Accelerator-based ion-beam irradiation has been widely used to mimic the effects of neutron radiation damage in nuclear reactors. However, ion radiation is most often monodisperse in the incoming ions’ momentum direction, leading to excessive polarization in defect distribution, while the scattering under neutron irradiation is often more isotropic and has less radiation-induced polarization. Mitigation of the excess-polarization as well as the damage non-uniformity artifact might be crucial for making the simulation of neutron radiation by ion-beam radiation more realistic. In this work, a general radiation polarization theory in treating radiation as external polar stimuli is established to understand the natural material responses in different contexts, and the possibility to correct the defect polarization artifact in ion-beam irradiation. Inspired by Magic Angle Spinning in Nuclear Magnetic Resonance, we present a precise sample spinning strategy to reduce the point-defect imbalance effect in ion-beam irradiation. It can be seen that with optimized surface inclination angle and the axis of sample rotation, the vacancy-interstitial population imbalance, as well as the damage profile non-uniformity in a designated region in the target are both reduced. It is estimated that sample spinning frequency on the order of kHz should be sufficient to scramble the ion momentum monodispersity for commonly taken ion fluxes and dose rates, which is experimentally feasible.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Inference of Multichannel r -process Element Enrichment in the Milky Way Using Binary Neutron Star Merger Observations

Observations of GW170817 strongly suggest that binary neutron star (BNS) mergers produce rapid neutron-capture nucleosynthesis ( r -process) elements. However, it remains an open question whether these mergers can account for all the r -process element enrichment in the Milky Way’s history. Here, we constrain the contributions of the BNS channel using astrophysical neutron star observations. The rate and mass distributions are constrained by LIGO/Virgo/Kagra through the latest catalog GWTC-3, the neutron star equation of state by gravitational-wave, radio, and X-ray observations, and the delay time distribution by short gamma-ray burst (GRB) host galaxy associations. We present a Bayesian framework to consistently combine these observations with abundance information to quantify the contribution and uncertainties of single and multiple astrophysical enrichment sources, and obtain a distribution of per-event BNS r -process element yields consistent with geophysical and astrophysical abundance constraints. We then adopt a Galactic chemical evolution model assuming an instantaneous and fixed amount of Fe enrichment from core-collapse supernovae, and show that BNS-only enrichment scenarios remain inconsistent with the observed r-process abundance trend of disk stars in the Galaxy even with the uncertainties in BNS merger observations. Using stellar abundance observations instead of the short GRB constraints, we can infer a shorter BNS delay time distribution with power-law index α ≤ −2.0 and minimum delay time ${t}_{{\rm{\min }}}\leqslant 40$ Myr at 90% confidence, consistent with detailed Galactic chemical evolution models. Such delay times are in tension with those predicted by standard BNS formation models. Alternatively, we confirm that a two-channel scenario, in which the second channel tracks the star formation history without significant delay, can account for both Galactic stellar and short GRB observations. We estimate that 45%–90% of the r -process abundance in the Milky Way today would have been produced by this star formation-tracking channel, rather than BNS mergers with significant delay times.

gravitational wave astronomy↗
DOE OSTI2021

Nucleon off-shell structure and the free neutron valence structure from A=3 inclusive electron scattering measurements

Understanding the differences between the distribution of quarks bound in protons and neutrons is key for constraining the mechanisms of SU(6) spin-flavor symmetry breaking in Quantum Chromodynamics (QCD). While vast amounts of proton structure measurements were done, data on the structure of the neutron is much more spars as experiments typically extract the structure of neutrons from measurements of light atomic nuclei using model-dependent corrections for nuclear effects. Recently the MARATHON collaboration performed such an extraction by measuring inclusive deep-inelastic electron-scattering on helium-3 and tritium mirror nuclei where nuclear effects are expected to be similar and thus be suppressed in the helium-3 to tritium ratio. Here we evaluate the model dependence of this extraction by examining a wide range of models including the effect of using instant-form and light-cone nuclear wave functions and several different parameterizations of nucleon modification effects, including those with and without isospin dependence. We find that, while the data cannot differentiate among the different models of nuclear structure and nucleon modification, they consistently prefer a neutron-to-proton structure function ratio of at $x_B \rightarrow 1$ of $\sim 0.4$ with a typical uncertainty ($1\sigma$) of $\sim0.05$ and $\sim0.10$ for isospin-independent and isospin-dependent modification models, respectively. While strongly favoring SU(6) symmetry breaking models based on perturbative QCD and the Schwinger-Dyson equation calculation, the MARATHON data do not completely rule out the scalar di-quark models if an isospin-dependent modification exist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Measured Thermal And Fast Neutron Fluence Rates ATR Cycle 166B

This report contains the thermal (2200 m/s) and fast (E>1MeV) neutron fluence rate data for ATR Cycle 166B which were measured by the Radiation Measurements Laboratory (RML) as requested by the Power Reactor Programs (ATR Experiments) Radiation Measurements Work Order. This report contains fluence rate values corresponding to the particular elevations (relative to the 80-ft. core elevation) where the measurements were taken. The data in this report consists of (1) a table of the ATR power history and distribution, (2) a hard copy listing of all thermal and fast neutron fluence rates, (3) plots of both the thermal and fast neutron fluence rates, and (4) an optical record (compact disk) containing a listing of the thermal and the fast neutron fluence rates, their assigned elevations and proper header identification of all monitor positions contained herein. The fluence rates reported are for the average power levels given in the table of power history and distribution. No neutron monitors were located in the N-2 position during this cycle.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Measured Thermal And Fast Neutron Fluence Rates ATR Cycle 169A

This report contains the thermal (2200 m/s) and fast (E>1MeV) neutron fluence rate data for ATR Cycle 169A which were measured by the Radiation Measurements Laboratory (RML) as requested by the Power Reactor Programs (ATR Experiments) Radiation Measurements Work Order. This report contains fluence rate values corresponding to the particular elevations (relative to the 80-ft. core elevation) where the measurements were taken. The data in this report consists of ( 1) a table of the ATR power history and distribution, (2) a hard copy listing of all thermal and fast neutron fluence rates, (3) plots of both the thermal and fast neutron fluence rates, and (4) an optical record (compact disk) containing a listing of the thermal and the fast neutron fluence rates, their assigned elevations and proper header identification of all monitor positions contained herein. The fluence rates reported are for the average power levels given in the table of power history and distribution. No neutron monitors were located in the N-2 position during this cycle. All "SR" holder monitor wires for this cycle are 55 inches long. This length allows measurements from 24.73 inches above core midplane to 25.25 inches below core midplane. The distance from the end of the wires to the first count position was 4.25 inches for all wires counted from this cycle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Two-neutron knockout as a probe of the composition of states in Mg 22 , Al 23 , and Si 24

Simpson and Tostevin proposed that the width and shape of exclusive parallel momentum distributions of the A - 2 residue in direct two-nucleon knockout reactions carry a measurable sensitivity to the nucleon single-particle configurations and their couplings within the wave functions of exotic nuclei. We further report here on the first benchmarks and use of this new spectroscopic tool. Exclusive parallel momentum distributions for states in the neutron-deficient nuclei Mg 22 , Al 23 , and Si 24 populated in such direct two-neutron removal reactions were extracted and compared to predictions combining eikonal reaction theory and shell-model calculations. For the well-known Mg 22 and Al 23 nuclei, measurements and calculations were found to agree, supporting the dependence of the parallel momentum distribution width on the angular momentum composition of the shell-model two-neutron amplitudes. In Si 24 , a level at 3439(9) keV, of relevance for the important Al 23 ( p , γ ) Si 24 astrophysical reaction rate, was confirmed to be the 2 2 + state, whereas the 4 1 + state, expected to be strongly populated in two-neutron knockout, was not observed. This puzzle is resolved by theoretical considerations of the Thomas-Ehrman shift, which also indicates that a previously reported 3471-keV state in Si 24 is, in fact, the ( 0 2 + ) level with one of the largest experimental mirror-energy shifts ever observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Transient convective delayed neutron precursor of {sup 235}U for the molten salt reactor experiment

This work investigates the delayed neutron precursors from the fission of {sup 235}U in molten salt reactors. The six delayed neutron groups predict the spatial distribution of the advective transport for radioactive mass transfer in the molten-salt reactor experiment (MSRE). The Mole code, using the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE), was used to compute the delayed neutron concentration of the whole system. This approach can be used for analysis of reactivity in circulating conditions. The effect of delayed neutron concentrations on the static or dynamic behavior of the system in MSRE was analyzed. Using this approach, the pattern and validation of delayed neutron for a range of parameters of interest are reported and discussed herein. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Measurement of the neutron charge radius and the role of its constituents

The neutron is a cornerstone in our depiction of the visible universe. Despite the neutron zero-net electric charge, the asymmetric distribution of the positively- (up) and negatively-charged (down) quarks, a result of the complex quark-gluon dynamics, lead to a negative value for its squared charge radius, $\langle r^2_n\rangle$. The precise measurement of the neutron’s charge radius thus emerges as an essential part of unraveling its structure. Here we report on a $\langle r^2_n\rangle$ measurement, based on the extraction of the neutron electric form factor, $G^n_E$, at low four-momentum transfer squared (Q 2 ) by exploiting the long known connection between the N → Δ quadrupole transitions and the neutron electric form factor. Our result, $\langle r^2_n\rangle$=-0.110±0.008(fm 2 ), addresses long standing unresolved discrepancies in the $\langle r^2_n\rangle$ determination. The dynamics of the strong nuclear force can be viewed through the precise picture of the neutron’s constituent distributions that result into the non-zero $\langle r^2_n\rangle$ value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Wireless Instrumented RB Experiment Preliminary Design and Analysis

The ability to deploy new nuclear fuels for current or future reactor concepts requires carefully designed experiments to generate data to support fuel qualification. Ideally these experiments would include state of-the-art sensing to maximize the amount of in situ data that can be collected during operation. Furthermore, advanced reactor systems can take advantage of integrated in-core sensing technologies to maximize fuel utilization, reduce unnecessary conservativism in design margins, and improve operator’s understanding of limiting peaking factors. Before any novel sensing technologies can be readily adopted for nuclear applications, they must first demonstrate acceptable performance in test reactors. This report summarizes the preliminary design and analysis of the most highly instrumented irradiation experiment ever performed in the removable beryllium (RB) positions of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). The Wireless Instrumented RB Experiment 2021 (WIRE-21) will test a wide range of sensors including wireless sensors being developed by Westinghouse Electric Company (WEC) that could provide in situ measurements of peak fuel temperatures and fuel rod pressurization due to fission gas release. The ability to wirelessly transmit a signal through the fuel rod’s cladding is critical to improving fuel monitoring capabilities without requiring signal penetrations through the cladding pressure boundary, which would significantly impact fuel fabrication, handling, and operation. Other sensors that will be tested in WIRE-21 include an array of thermocouples, self-powered neutron detectors (SPNDs), and spatially distributed fiber-optic temperature sensors. More generally, WIRE-21 will establish a flexible irradiation vehicle design to allow accelerated, economical testing of advanced sensor technologies while leveraging the extremely high neutron flux that is available in HFIR. This report summarizes the mechanical design for WIRE-21, the experimental test matrix, initial neutronic and thermal design analyses, and the active monitoring and control system enhancements necessary to support testing of advanced sensor technologies. The containment for WIRE-21 is similar to previous RB irradiation vehicles but includes a few modifications, most notably the use of integrated compression seals to pass a larger number of sensor leads through the experiment’s pressure boundary. In addition to the sensor leads, inert gas lines are passed into the experiment to enable active temperature control and the ability to pneumatically actuate a bellows-driven pressure sensor. WIRE-21 is targeting temperatures (300–350°C) and neutron fluence levels (~10 22 n/cm 2 ) relevant to light water reactors (LWRs), but the flexible design of the experiment vehicle allows much higher operating temperatures (>1,100°C). Neutronic calculations determine the neutron flux conditions as well as the nuclear heating within the experiments. These results are used as inputs to detailed thermal finite element calculations, which are required to evaluate the complex, three-dimensional heat transfer that occurs within WEC’s wireless sensor enclosures. Initial results show that the temperatures of the sensors’ enclosures and the metal bellows can be operated near the temperature range of LWR coolants and cladding while simultaneously increasing the temperature of a surrogate fuel material to values in the range of 800–1200°C to simulate centerline fuel temperatures during LWR operation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Mass-asymmetric fission of 205,207,209 Bi at energies close to the fission barrier using proton bombardment of 204,206,208 Pb

Background: Recent observation of mass-asymmetric fission in neutron-deficient Hg and Pt nuclei has reignited interest in fission fragment mass distributions close to Pb. Investigations at energies close to the fission barrier, where mass-asymmetric fission is expected to be most obvious and the sensitivity to shell effects is maximized, are limited in this mass region. Purpose: To measure fission mass distributions for 205,207,209 Bi nuclei at the lowest possible excitation energies to determine how the mass distributions change with excitation energy and the neutron number of the compound nucleus. Method: Proton beams bombarding targets of 204,206,208 Pb were used to study the fission of 205,207,209 Bi at energies from just above to 10 MeV above their fission barriers. Fission fragments were measured using the CUBE fission spectrometer. Fission fragment mass distributions were determined using a newly developed time difference analysis method. Mass distributions were characterized by triple-Gaussian fits to determine the systematic trends across each isotope with excitation energy. Results: Measured mass distributions of all three Bi isotopes exhibit a component of mass-asymmetric fission at all energies studied. The probability of mass-asymmetric fission decreases significantly with increasing excitation energy, from ≈70 to ≈40% over a 10-MeV range. Comparisons between the three Bi isotopes hint at an increase in the mass-symmetric fission yield with increasing neutron number, which could be due to a decrease in the difference between the symmetric and asymmetric fission barriers. The centroids of the mass-asymmetric peaks suggest that several deformed shell gaps in the fission fragments could be contributing to the presence of the mass-asymmetric fission mode with Z light ≃ 38, Z heavy ≃ 45, and N light ≃ 56 all present in the fission fragments. Conclusions: Measurements of fission mass distributions at the lowest possible excitation energies above the fission barrier provide an excellent platform to investigate the origins of the mass-asymmetric fission mode. Finally, further systematic measurements at these energies offer an opportunity to rigorously test new models of fission in this mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

New Measurements of the Deuteron-to-Proton 𝐹 2 Structure-Function Ratio

Nucleon structure functions, as measured in lepton-nucleon scattering, have historically provided a critical observable in the study of partonic dynamics within the nucleon. However, at very large parton momenta, it is both experimentally and theoretically challenging to extract parton distributions due to the probable onset of nonperturbative contributions and the unavailability of high-precision data at critical kinematics. Extraction of the neutron structure and the d quark distribution have been further challenging because of the necessity of applying nuclear corrections when utilizing scattering data from a deuteron target to extract the free neutron structure. However, a program of experiments has been carried out recently at the energy-upgraded Jefferson Lab electron accelerator aimed at significantly reducing the nuclear correction uncertainties on the d quark distribution function at large partonic momentum. This allows leveraging the vast body of deuterium data covering a large kinematic range to be utilized for d quark parton distribution function extraction. In this Letter, we present new data from experiment E12-10-002, carried out in Jefferson Lab Experimental Hall C, on the deuteron to proton cross section ratio at large Bjorken 𝑥. These results significantly improve the precision of existing data and provide a first look at the expected impact on quark distributions extracted from parton distribution function fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Application of recoil-imaging time projection chambers to directional neutron background measurements in the SuperKEKB accelerator tunnel

Gaseous time projection chambers (TPCs) with high readout segmentation are capable of reconstructing detailed 3D ionization distributions of nuclear recoils resulting from elastic neutron scattering. Here, using a system of six compact TPCs with pixel ASIC readout, filled with a 70:30 mixture of He:CO 2 gas, we analyze the first directional measurements of beam-induced neutron backgrounds in the tunnel regions surrounding the Belle II detector at the SuperKEKB e + e – collider. With the use of 3D recoil tracking, we show that these TPCs are capable of maintaining nearly 100% nuclear recoil purity to reconstructed ionization energies (E reco ) as low as 5 keV ee . Using a large sample of Monte-Carlo (MC)-simulated 4 He, 12 C, and 16 O recoil tracks, we find consistency between predicted and measured recoil energy spectra in five of the six TPCs, providing useful validation of the neutron production mechanisms modeled in simulation. Restricting this sample to 4 He recoil tracks with E reco > 40 keV ee , we further demonstrate axial angular resolutions within 8° and we introduce a procedure that under suitable conditions, correctly assigns the vector direction to 91% of these simulated 4 He recoils. Applying this procedure to assign vector directions to measured 4 He recoil tracks, we observe consistency between the angular distributions of observed and simulated recoils, providing first experimental evidence of localized neutron “hotspots” in the accelerator tunnel. Observed rates of nuclear recoils in these TPCs suggest that simulation overestimates the neutron flux from these hotspots. Despite this, we estimate these hotspots to produce the majority of neutron backgrounds in the accelerator tunnel at SuperKEKB’s target luminosity of 6.3 x 10 35 cm –2 s –1 , making them important regions to continue to monitor.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A liquid hydrogen tube moderator arrangement for SNS second target station

Liquid hydrogen filled tubes arranged in a triangular shape surrounded by light-water premoderators were investigated as cold moderators coupled to a neutron production zone of a tungsten target fed by a 1.3 GeV proton beam in a short-pulse mode. A moderator concept optimized for the tube length, premoderator thickness, and target position was found to deliver the highest pulse-integrated neutron brightness suitable to send cold neutron beams for neutron scattering instrumentation geared toward medium resolution and small samples. Reducing the tube diameter from 50 to 10 mm boosts the cold neutron brightness (E < 10 meV) by a factor of 2.25, also with enhancements seen by downsizing the height of cylindrical moderators in earlier studies. Attempts were made to characterize this novel moderator system with regard to pulse shapes, local brightness distribution, and angular distribution of emission at the neutron emission port.

Gallmeier, Franz X. (ORCID:0000000171500764)↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Multi-scale investigation of heterogeneous swift heavy ion tracks in stannate pyrochlore

Er 2 Sn 2 O 7 pyrochlore was irradiated with swift heavy Au ions (2.2 GeV), and the induced structural modifications were systematically examined using complementary characterization techniques including transmission electron microscopy (TEM), X-ray diffraction (XRD), and neutron total scattering with pair distribution function (PDF) analysis. Each technique probes different aspects and length scales of the transformed material regions. TEM revealed a core–shell ion track structure—an amorphous core surrounded by a disordered, anion-deficient fluorite shell—which was confirmed by XRD. Neutron total scattering, with sensitivity to the oxygen sublattice, provided relative fractions of amorphous and disordered fluorite phases and confirmed the presence of a defective pyrochlore phase, which largely maintains its structural ordering but is clearly distinct from the pristine pyrochlore matrix. This defect-rich pyrochlore phase forms a halo extending radially beyond the well-characterized core–shell track morphology observed in electron micrographs. Despite their differing long-range periodicity, the short-range structures of the amorphous, disordered, and defective pyrochlore phases are all modeled well with a weberite-type configuration. Evolution of the phase fractions with increasing ion fluence was examined to ascertain the phase-to-phase pathways that occur during primary and secondary ion impact. Furthermore, this approach extends knowledge about the multi-scale response of stannate pyrochlores to swift heavy ion irradiation in the electronic energy loss regime and improves existing track-overlap models.

36 MATERIALS SCIENCE↗