Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “filled”

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,297 records · Page 72

Magnetic properties and pairing tendencies of the iron-based superconducting ladder BaFe 2 S 3 : Combined ab initio and density matrix renormalization group study

The recent discovery of superconductivity under high pressure in the two-leg ladder compound BaFe 2 S 3 [H. Takahashi et al., Nat. Mater. 14, 1008 (2015)] opens a broad avenue of research, because it represents the first report of pairing tendencies in a quasi-one-dimensional iron-based high-critical-temperature superconductor. Similarly, as in the case of the cuprates, ladders and chains can be far more accurately studied using many-body techniques and model Hamiltonians than their layered counterparts, particularly if several orbitals are active. In this publication, we derive a two-orbital Hubbard model from first principles that describes individual ladders of BaFe 2 S 3 . The model is studied with the density matrix renormalization group. These first reported results are exciting for two reasons: (i) at half-filling, ferromagnetic order emerges as the dominant magnetic pattern along the rungs of the ladder, and antiferromagnetic order along the legs, in excellent agreement with neutron experiments; and (ii) with hole doping, pairs form in the strong coupling regime, as found by studying the binding energy of two holes doped on the half-filled system. In addition, orbital selective Mott phase characteristics develop with doping, with only oneWannier orbital receiving the hole carriers while the other remains half-filled. Lastly, these results suggest that the analysis of models for iron-based two-leg ladders could clarify the origin of pairing tendencies and other exotic properties of iron-based high-critical-temperature superconductors in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/t from 2 to 8, for a range of densities including half-filling and n = 0.3,0.5,0.6, 0.75, and 0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in this paper in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to 16 × 16 are also provided for benchmark purposes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ferromagnetism beyond Lieb's theorem

The noninteracting electronic structures of tight-binding models on bipartite lattices with unequal numbers of sites in the two sublattices have a number of unique features, including the presence of spatially localized eigenstates and flat bands. When a uniform on-site Hubbard interaction U is turned on, Lieb proved rigorously that at half-filling (ρ = 1) the ground state has a nonzero spin. In this paper we consider a “CuO 2 lattice” (also known as “Lieb lattice,” or as a decorated square lattice), in which “d orbitals” occupy the vertices of the squares, while “ p orbitals” lie halfway between two d orbitals; both d and p orbitals can accommodate only up to two electrons. We use exact determinant quantum Monte Carlo (DQMC) simulations to quantify the nature of magnetic order through the behavior of correlation functions and sublattice magnetizations in the different orbitals as a function of U and temperature; we have also calculated the projected density of states, and the compressibility. Here, we study both the homogeneous (H) case, U d = U p , originally considered by Lieb, and the inhomogeneous (IH) case, U d ≠ U p . For the H case at half-filling, we found that the global magnetization rises sharply at weak coupling, and then stabilizes towards the strong-coupling (Heisenberg) value, as a result of the interplay between the ferromagnetism of like sites and the antiferromagnetism between unlike sites; we verified that the system is an insulator for all U. For the IH system at half-filling, we argue that the case U p ≠ U d falls under Lieb's theorem, provided they are positive definite, so we used DQMC to probe the cases U p = 0, U d = U and Up = U, U d = 0. We found that the different environments of d and p sites lead to a ferromagnetic insulator when U d = 0; by contrast, U p = 0 leads to to a metal without any magnetic ordering. In addition, we have also established that at density ρ = 1/3, strong antiferromagnetic correlations set in, caused by the presence of one fermion on each d site.

36 MATERIALS SCIENCE↗

Shape of magnetic domain walls formed by coupling to mobile charges

Magnetic domain walls, which are crucially important in both fundamental physics and technical applications, often have a preference in their form due to many different origins, such as the crystalline shape, lattice symmetry, and magnetic anisotropy. Here we theoretically investigate yet another origin stemming from the coupling to mobile charges in itinerant magnets. Performing a large-scale numerical simulation in a minimal model for itinerant magnets, i.e., the Kondo lattice model with classical localized spins, we show that the shape of magnetic domain walls depends on the electronic band structure and electron filling. While Néel and 120° antiferromagnetic states do not show a strong preference in the shape of domain walls, noncoplanar spin states with scalar chiral ordering have distinct directional preferences of the domain walls depending on the electron filling. We find that the directional preference is rationalized by the wave-number dependence of the effective magnetic interactions induced by the mobile charges, which are set by the band structure and electron filling. Finally, we also observe that, in the noncoplanar chiral states, an electric current is induced along the domain walls owing to the spin Berry phase mechanism, with different spatial distributions depending on whether the bulk state is metallic or insulating.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cascades between Light and Heavy Fermions in the Normal State of Magic-Angle Twisted Bilayer Graphene

We present a framework for understanding the cascade transitions and the Landau level degeneracies of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a band minimum at the moiré Brillouin zone center with a small mass that compares well with the experiment; if towards the charge neutrality point, then it has a much larger mass and a higher degeneracy. At a nonzero density away from an integer filling the excitations interact. The system on the small mass side has a large bandwidth and forms a Fermi liquid. On the large mass side the bandwidth is narrow, the compressibility is negative and the Fermi liquid is likely unstable. This explains the observed sawtooth features in compressibility, the Landau fans pointing away from charge neutrality and their degeneracies. The framework sets the stage for superconductivity at lower temperatures.

36 MATERIALS SCIENCE↗

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund’s Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

For this work, we apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)], to obtain its Kondo lattice limit. In this limit localized $\mathcal{f}$ electrons on a triangular lattice interact with topological conduction $\mathcal{c}$ electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0, ±1, ±2 are controlled by the heavy $\mathcal{f}$ electrons, while ν=±3 is at the border of a phase transition between two $\mathcal{f}$-electron fillings. For ν=0, ±1, ±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the $\mathcal{f}$ moments in the full model and analytically prove the SU(4) Hund’s rule for the ground state which maintains that two $\mathcal{f}$ electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite “bond-operators” method. We then compute the excitation spectrum of the $\mathcal{f}$ moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $\mathcal{f}$ moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Correlation effects in magic-angle twisted bilayer graphene: An auxiliary-field quantum Monte Carlo study

Magic-angle twisted bilayer graphene (MATBG) presents a fascinating platform for investigating the effects of electron interactions in topological flat bands. The Bistritzer-MacDonald (BM) model provides a simplified quantitative description of the flat bands. Introducing long-range Coulomb interactions leads to an interacting BM (IBM) Hamiltonian, a momentum-space continuum description which offers a very natural starting point for many-body studies of MATBG. Accurate and reliable many-body computations in the IBM model are challenging, however, and have been limited mostly to special fillings or smaller lattice sizes. We employ a state-of-the-art auxiliary-field quantum Monte Carlo (AFQMC) method to study the IBM model, which constrains the sign problem to enable accurate treatment of large system sizes. We determine ground-state properties and quantify errors compared to mean-field theory calculations. Our calculations identify correlated metal states and their competition with the insulating Kramers intervalley-coherent state at both half-filling and charge neutrality. Additionally, we investigate one- and three-quarter fillings, and examine the effect of many-body corrections beyond single Slater determinant solutions. We discuss the effect that details of the IBM Hamiltonian have on the results, including different forms of double-counting corrections, and the need to establish and precisely specify many-body Hamiltonians to allow more direct and quantitative comparisons with experiments in MATBG. Published by the American Physical Society 2025

Xiao, Zhi-Yu (ORCID:0000000219531579)↗

Microscopic origin of temperature-dependent magnetism in spin-orbit-coupled transition metal compounds

A few 4 d and 5 d transition metal compounds with various electron fillings were recently found to exhibit magnetic susceptibilities χ and magnetic moments that deviate from the well-established Kotani model. This model has been considered for decades to be the canonical expression for descriing the temperature dependence of magnetism in systems with nonnegligible spin-orbit coupling effects. In this paper, we uncover the origin of such discrepancies and determine the applicability and limitations of the Kotani model by calculating the temperature dependence of the magnetic moments of a series of 4 d (Ru-based) and 5 d (W-based) systems at different electron fillings. For this purpose, we perform exact diagonalization of -derived relativistic multiorbital Hubbard models on finite clusters and compute their magnetic susceptibilities. Comparison with experimentally measured magnetic properties indicates that contributions such as a temperature-independent χ 0 background, crystal field effects, Coulomb and Hund's couplings, and intersite interactions—not included in the Kotani model—are especially crucial for correctly describing the temperature dependence of χ and magnetic moments at various electron fillings in these systems. Based on our results, we propose a generalized approach beyond the Kotani model to accurately describe their magnetism. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spaced Orthocyclic Winding Pattern for Improved Die Compressed Coils

Winding fill, coil fill, or packing factor, is critical to the performance of many types of electromagnetic devices from electric machines to actuators, transformers, inductors, etc. One technique to increase the packing of conductors in the winding area is to compress round magnet wire conductors in a die, so that they deform and air voids between coil turns are removed. The winding layout of the coil turns on a bobbin before compression has a direct impact on the achievable coil fill and the damage to the coil insulation or conductors after compression. It is thought that winding the turns in an orthocyclic pattern on the bobbin before compression would minimize damage to the turns. This however causes turns to displace during the winding compression process. A new pre-compression winding pattern is proposed, a spaced orthocyclic winding, to minimize turn displacement and damage. A design of experiments study has been carried out to optimize this winding method and to compare it to the traditional orthocyclic winding.

33 ADVANCED PROPULSION SYSTEMS↗

NDX: Neutron Dose Rate Meters with Extended Capabilities

Particle accelerators, irradiation facilities, and nuclear power plants would benefit from using neutron dose rate meters featuring improved energy response covering full range of expected neutron energies, improved dynamic range of measured dose rates while providing the real-time data, and improved ability to operate in harsh radiation environments including heavy photon backgrounds and pulsed beams. The novel design of the neutron moderator for the NDX dose rate meter implements a layer of a Beryllium-loaded medium Z material which, similar to the earlier known WENDI detectors, allows the device to detect higher energy neutron flux and convert it to dose. Using Beryllium helps to fill the known dip in the response functions around 5-50 MeV. Optimized design parameter set, obtained using FLUKA Monte Carlo simulations, features the NDX moderator made of polyethylene embedding a layer of BeCu alloy inside, with the dose rate response constant within ±10% in the neutron energy range from about 100 keV to 500 MeV, and ±25% essentially in the full energy range above thermal. A new sensor technique implemented in the NDX is based on using the Ionization Chamber (IC) filled with 3He gas to capture the moderated neutron fluence, and includes a symmetric IC filled with regular 4He to evaluate the contributions due to the non-neutron ionizing radiation. Advantages of the method include large accessible dynamic range, relative insensitivity of the IC probes to photon radiation fields, and the ability to subtract their contribution. Two prototype NDX detectors were built, calibrated, and operated in various conditions at the CEBAF 12 GeV electron accelerator at JLab. The first results showed reliable neutron dose rate measurements in the presence of overwhelming photon radiation fields, and at practical high dose rate levels of up to approximately 10 Sv/h. Detector design, operation, first results, and further development plans are presented.

Degtiarenko, Pavel V.↗

SU-F-T-141: Proton Dose Validation in a Phantom Beyond TRUFILL N-BCA Embolization Glue

Purpose: To validate the treatment planning system predicted proton dose beyond a heterogeneity (n-BCA glue) by making a measurement in a custom acrylic phantom. Methods: A custom cubic acrylic phantom was designed for this experiment. A container was designed to fit in the phantom. This container was filled with TRUFILL™ n-Butyl Cyanoacrylate(n-BCA) glue. When the container was placed in the phantom, its center was at a distance of 7.4cm from the entrance. This depth allows us to make measurements around the center of modulation of a 126 MeV proton beam with a 3cm spread-out-Bragg peak. To make measurements at other beam energies, additional acrylic can be added in front of the phantom, to adjust the depth of the heterogeneity. A diamond detector was cross calibrated against a standard cylindrical ion chamber in a 126MeV beam. The diamond detector was then used to make dose measurements beyond the inhomogeneity. The measurement was repeated with the container filled with water. Several measurements were made at each setup, to check reproducibility of measurements. Results: For the same number of Tic3R1 counts, the dose measured with the diamond detector beyond n-BCA glue was 1.053 times the dose measured beyond the water filled container. This result is in agreement with the measured stopping power of glue (1.06). These measurements were in agreement with the dose predicted by the treatment planning system when the electron density of the heterogeneity was replaced with 1.06 before the dose calculation. Conclusion: Our initial measurements validate the dose predicted by our treatment plan in the presence of heterogeneity in a phantom. The material tested (n-BCA glue) is commonly used in the treatment of AVM’s prior to an SRS treatment. An error in dose predicted by the treatment plan in the presence of the glue can be detrimental in a single fraction high dose SRS treatment I received the n-BCA liquid embolic system samples from Codman and Shurtleff, Inc.

60 APPLIED LIFE SCIENCES↗

SU-F-T-171: Manufacturing Cost Effective Heterogeneous Phantoms for Use in Proton Therapy

Purpose: To study the feasibility of 3D printing cost effective heterogeneous phantoms for use in proton therapy treatment planning quality assurance. Methods: A desktop 3D printer was utilized to create a series of 2 cm × 2 cm × 4 cm PLA plastic blocks of varying fill materials and hexagonal fill pattern. The blocks were than tested when filled with air, polyurethane foam, paraffin, silicone, and caulk of calcium carbonate – acrylic polymer blend. The blocks were evaluated with a “GE Lightspeed” 16 slice CT scanner and average CT# of the materials’ centers evaluated. Blocks were then placed into a custom aperture fitted to a Mevion Proton system to determine the relative stopping power of each. Scans were performed in water tank with Marcus type parallel plate chamber under a beam with a range of 15 cm and modulation of 2 cm. Shifts in range occurring relative to the 80% distal edge of the open SOBP were evaluated. Results: The CT#s of the blocks were plotted against their measured relative stopping power. This curve was compared to that which is in clinical use. While the trend agrees generally, specific differences between the relative stopping powers were as great as 10%. Conclusion: We have demonstrated that it is possible to utilize different cost effective materials in the manufacturing of phantoms for use in proton therapy. While different materials may provide better agreement to established calibration curves, a custom curve specific to the materials used may be utilized to accurately predict proton treatment dose distributions.

60 APPLIED LIFE SCIENCES↗

Technical Note: Evaluation of a 160-mm/256-row CT scanner for whole-heart quantitative myocardial perfusion imaging

Purpose: The authors investigated the performance of a recently introduced 160-mm/256-row CT system for low dose quantitative myocardial perfusion (MP) imaging of the whole heart. This platform is equipped with a gantry capable of rotating at 280 ms per full cycle, a second generation of adaptive statistical iterative reconstruction (ASiR-V) to correct for image noise arising from low tube voltage potential/tube current dynamic scanning, and image reconstruction algorithms to tackle beam-hardening, cone-beam, and partial-scan effects. Methods: Phantom studies were performed to investigate the effectiveness of image noise and artifact reduction with a GE Healthcare Revolution CT system for three acquisition protocols used in quantitative CT MP imaging: 100, 120, and 140 kVp/25 mAs. The heart chambers of an anthropomorphic chest phantom were filled with iodinated contrast solution at different concentrations (contrast levels) to simulate the circulation of contrast through the heart in quantitative CT MP imaging. To evaluate beam-hardening correction, the phantom was scanned at each contrast level to measure the changes in CT number (in Hounsfield unit or HU) in the water-filled region surrounding the heart chambers with respect to baseline. To evaluate cone-beam artifact correction, differences in mean water HU between the central and peripheral slices were compared. Partial-scan artifact correction was evaluated from the fluctuation of mean water HU in successive partial scans. To evaluate image noise reduction, a small hollow region adjacent to the heart chambers was filled with diluted contrast, and contrast-to-noise ratio in the region before and after noise correction with ASiR-V was compared. The quality of MP maps acquired with the CT system was also evaluated in porcine CT MP studies. Myocardial infarct was induced in a farm pig from a transient occlusion of the distal left anterior descending (LAD) artery with a catheter-based interventional procedure. MP maps were generated from the dynamic contrast-enhanced (DCE) heart images taken at baseline and three weeks after the ischemic insult. Results: Their results showed that the phantom and animal images acquired with the CT platform were minimally affected by image noise and artifacts. For the beam-hardening phantom study, changes in water HU in the wall surrounding the heart chambers greatly reduced from >±30 to ≤ ± 5 HU at all kVp settings except one region at 100 kVp (7 HU). For the cone-beam phantom study, differences in mean water HU from the central slice were less than 5 HU at two peripheral slices with each 4 cm away from the central slice. These findings were reproducible in the pig DCE images at two peripheral slices that were 6 cm away from the central slice. For the partial-scan phantom study, standard deviations of the mean water HU in 10 successive partial scans were less than 5 HU at the central slice. Similar observations were made in the pig DCE images at two peripheral slices with each 6 cm away from the central slice. For the image noise phantom study, CNRs in the ASiR-V images were statistically higher (p < 0.05) than the non-ASiR-V images at all kVp settings. MP maps generated from the porcine DCE images were in excellent quality, with the ischemia in the LAD territory clearly seen in the three orthogonal views. Conclusions: The study demonstrates that this CT system can provide accurate and reproducible CT numbers during cardiac gated acquisitions across a wide axial field of view. This CT number fidelity will enable this imaging tool to assess contrast enhancement, potentially providing valuable added information beyond anatomic evaluation of coronary stenoses. Furthermore, their results collectively suggested that the 100 kVp/25 mAs protocol run on this CT system provides sufficient image accuracy at a low radiation dose (<3 mSv) for whole-heart quantitative CT MP imaging.

60 APPLIED LIFE SCIENCES↗

Poster - 38: On the physical and dosimetric properties of 3D printed electron bolus fabricated using polylactic acid

Purpose: 3D printing technology could simplify and improve electron bolus fabrication. The purpose of this study was to characterize the density, dimensional accuracy, uniformity, and attenuation of PLA boluses fabricated with a low-cost 3D printer. Methods: Several solid square slabs were printed with specific requested dimensions and 100% infill using different fill patterns. These pieces were imaged using an x-ray flat panel imager in order to check for uniformity of the prints. Percentage depth doses (PDDs) were measured downstream of the slabs in solid water using a parallel plate chamber and compared to measurements in water in order to characterize attenuation. The dimensions of the PLA slabs were measured using digital calipers. The slabs were also weighed to find their density. Results: The fill pattern used to create boluses can affect the attenuation of the bolus. Fill patterns should be chosen carefully and quality assurance should be done for each printed piece. PLA causes the electron PDD to shift towards shallower depths, compared to water, by 1.7 mm for each centimeter of PLA. Agreement with design dimensions was within 1 mm in the plane of the printer bed, and within 1/3 of a millimeter (roughly the thickness of a single layer), perpendicular to the printer bed. Average density was in the range 1.20 – 1.22. Conclusions: 3D printing shows great promise for use in fabricating electron bolus. This work indicates that printed PLA can be a suitable material provided the increased attenuation is properly accounted for.

60 APPLIED LIFE SCIENCES↗

Interaction of a normally-incident plane wave with a nonlinear poroelastic fracture

While it has been recognized that a large amplitude incident wave upon a dry fracture can exhibit nonlinear seismic wave scattering due to its stress-dependent mechanical compliance, the impact of pore fluid in the fracture and a fluid-filled poroelastic background medium—features common for fractures in the Earth—are not well understood. As a first step toward an understanding of the nonlinear poroelastic response of elastic waves in fractured media, analytical approximate formulas are used for the amplitude and phase of a normally incident plane wave using a perturbation method, assuming a fluid-filled, highly compliant nonlinear interface embedded in a linear poroelastic solid. The stress-closure behavior of the fracture is modeled by nonlinear, poroelastic displacement-discontinuity boundary conditions (a linear–slip interface). The theory predicts that the static (“Direct current,” or DC) and higher-order-harmonic waves produced by the nonlinear scattering can be greatly reduced by the presence of fluid in the fracture. This, however, depends upon a number of parameters, including fracture compliance, fluid properties (compressibility and viscosity), and the permeability of the background medium, as well as environmental parameters such as the initial fluid pressure and stress acting on the fracture. The static effect produces low-frequency fluid pressure pulses when a finite-duration wave is incident upon the fracture—behavior unique to fluid-filled fractures within a poroelastic medium.

58 GEOSCIENCES↗

Space-tiled colloidal crystals from DNA-forced shape-complementary polyhedra pairing

Generating space-filling arrangements of most discrete polyhedra nanostructures of the same shape is not possible. However, if the appropriate individual building blocks are selected (e.g., cubes), or multiple shapes of the appropriate dimensions are matched (e.g., octahedra and tetrahedra) and their pairing interactions are subsequently forced, space-filled architectures may be possible. With flexible molecular ligands (polyethylene glycol–modified DNA), the shape of a polyhedral nanoparticle can be deliberately altered and used to realize geometries that favor space tessellation. In this work, 10 new colloidal crystals were synthesized from DNA-modified nanocrystal building blocks that differed in shapes and sizes, designed to form space-filling architectures with micron-scale dimensions. The insights and capabilities provided by this new strategy substantially expand the scope of colloidal crystals possible and provide an expanded tool kit for researchers interested in designing metamaterials.

Science & Technology - Other Topics↗

Rapid clay precipitation in explosion-induced fractures

Fractures within the earth control rock strength and fluid flow, but their dynamic nature is not well understood. As part of a series of underground chemical explosions in granite in Nevada, we collected and analyzed microfracture density data sets prior to, and following, individual explosions. Our work shows an ~4-fold increase in both open and filled microfractures following the explosions. Based on the timing of core retrieval, filling of some new fractures occurs in as little as 6 wk after fracture opening under shallow (<100 m) crustal conditions. These results suggest that near-surface fractures may fill quite rapidly, potentially changing permeability on time scales relevant to oil, gas, and geothermal energy production; carbon sequestration; seismic cycles; and radionuclide migration from nuclear waste storage and underground nuclear explosions.

58 GEOSCIENCES↗