Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ionizing radiation”

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 505 records · Page 28

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The effect of metastable atomic levels in neutral tungsten for erosion and impurity transport applications

Long-lived metastable states in neutral tungsten are shown to potentially impact measurements of tungsten erosion from plasma facing components. Time-dependent collisional radiative modeling of neutral tungsten is used to analyze the role of these states in tungsten emission and ionization. The large number of non-quasistatic atomic states in neutral tungsten can take on the order of milliseconds to reach equilibrium, depending on plasma conditions, causing erosion measurements to be affected by the metastable populations. Previous measurements using the 400.88 nm tungsten emission line could be affected by these non-quasistatic metastable effects. Therefore, a scheme for measuring the relative metastable fractions is proposed through simultaneous observation of multiple ultraviolet spectral lines of neutral W. The accuracy of gross erosion measurements could potentially be increased by inclusion of these previously unconsidered metastable effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Experimental Observation of the Tetrabromine Cluster Br 4 by Synchrotron Photoionization Mass Spectrometry

We present a first spectroscopic characterization of the homoatomic polyhalogen tetrabromine, Br 4 , in the gas phase. Photolysis of CHBr 3 at 248 nm is used to generate atomic bromine radicals in a flow tube reactor. Resulting combination products are detected by photoionization mass spectrometry at the Advanced Light Source of the Lawrence Berkeley National Laboratory. Interpretation of the experimental mass spectra is informed by calculated adiabatic ionization energies carried out at the CCSD(T)/aug-cc-pVTZ//M06-2X/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ//cam-B3LYP/6-311++g** levels of theory. Tunable VUV synchrotron radiation enables the collection of the mass-selected photoionization spectra by which Br 4 is assigned using Franck-Condon simulations of a Br 2 dimer with a stretched tetrahedral geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Amplified spontaneous emission in the extreme ultraviolet by expanding xenon clusters

Focused short-wavelength free-electron laser (FEL) pulses interacting with gas phase samples can induce by inner-shell ionization a short-lived population inversion, followed by coherent collective emission of directed, short, and strong radiation bursts. We extend our studies into the warm-dense matter (WDM) regime by investigating the nanoplasmas produced in an ensemble of nanometer-sized clusters by FEL irradiation. Here, additional pathways can also lead to strong, laserlike emission: Electron-ion collisions can yield a long-lived population inversion, and subsequent amplified spontaneous emission. We observe amplified spontaneous emission (ASE) in the extreme ultraviolet in xenon clusters excited by soft x-ray FEL pulses, we diagnose the generated nanoplasmas by fluorescence spectroscopy, and we study under various cluster and FEL parameters the directed ASE from the Xe 2+ 65 nm line. We show its exponential increase as a function of FEL irradiation power, and an accompanying collisional broadening of the emission spectra. These findings are corroborated by extensive numerical simulations based on theory, combining detailed hydrodynamic and kinetic simulations with time-dependent calculations of radiation transport, amplification, and collective emission in the WDM nanoplasma. Furthermore, our theoretical findings underline that population inversion is due to electron-ion collisions and that the observed decoherence processes can be empirically characterized by a phenomenological decoherence time in the range of 100–200 fs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pressurized Gas Beam Monitor for Extremely High Intensities – Phase II

A novel pressurized gas-filled multi-RF-cavity beam profile monitor has been studied that is simple and robust in high-radiation environments. Charged particles passing through each RF-cavity in the monitor produce intensity-dependent ionized plasma, which changes the gas permittivity. Standard RF techniques to measure the change in quality factor (Q) and frequency (f) as a function of time are then used to determine the change in permittivity and corresponding beam intensity in each cavity in the profile monitor. The sensitivity to beam intensity is adjustable using gas pressure and RF gradient. The performance of the gas-filled beam profile monitor has been numerically simulated to evaluate the sensitivity of permittivity measurements and the effectiveness of calibration strategies. An entire multi-RF-cavity system has been designed and a demonstration test has been prepared.

43 PARTICLE ACCELERATORS↗

Review of the 1st EUV Light Sources Code Comparison Workshop

We review the results of the 1 st Extreme Ultraviolet (EUV) Light Sources Code Comparison Workshop. The goal of this workshop was to provide a platform for specialists in EUV light source plasma modeling to benchmark and validate their numerical codes using well-defined case studies. Detailed consideration of a plethora of atomic collisional and radiative processes is required for modeling EUV light source plasmas. Eight institutions spanning four countries contributed data to the workshop. Two topics were addressed, namely (i) the atomic kinetics and radiative properties of tin plasmas under EUV-generating conditions and (ii) laser absorption in a fully ionized, one-dimensional hydrogen plasma. In this paper, we summarize the key findings of the workshop and outline plans for future iterations of the code comparison activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Examining Interactions of Uranyl(VI) Ions with Amino Acids in the Gas Phase

Gas-phase experiments, using electrospray ionization quadrupole ion trap mass spectrometry (ESI-QIT/MS), were conducted to probe basic interactions of the uranyl(VI) ion, UO 2 2+ , with selected natural amino acids, namely, L-cysteine (Cys), L-histidine (His), and L-aspartic acid (Asp), which strongly bind to metal ions. The simplest amino acid, glycine (Gly), was also studied for comparison. Cys, His, and Asp have additional potentially coordinating groups beyond the amino and carboxylic acid functional groups, specifically thiol in Cys, imidazole in His, and a second carboxylate in Asp. Gas-phase experiments comprised collision-induced dissociation (CID) of uranyl–amino acid complexes and competitive CID to assess the relative binding strength of different amino acids in the same uranyl complex. Reactivity of selected uranyl–amino acid complexes with water provided further insights into relative stabilities. In positive ion mode, CID and ensuing reactions with water suggested that uranyl–neutral AA binding strength decreased in the order His > Asp > Cys > Gly, which is similar to amino acid proton affinities. In negative ion mode, CID revealed a decreasing dissociation tendency in the order Gly >> His ≈ Cys > Asp, presumably reflecting a reverse enhanced binding to uranyl of the doubly deprotonated amino acids formed in CID.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Quantum-Path-Resolved Attosecond High-Harmonic Spectroscopy

Strong-field ionization of molecules releases electrons which can be accelerated and driven back to recombine with their parent ion, emitting high-order harmonics. This ionization also initiates attosecond electronic and vibrational dynamics in the ion, evolving during the electron travel in the continuum. Revealing this subcycle dynamics from the emitted radiation usually requires advanced theoretical modeling. Here we show that this can be avoided by resolving the emission from two families of electronic quantum paths in the generation process. The corresponding electrons have the same kinetic energy, and thus the same structural sensitivity, but differ by the travel time between ionization and recombination-the pump-probe delay in this attosecond self-probing scheme. We measure the harmonic amplitude and phase in aligned CO 2 and N 2 molecules and observe a strong influence of laser-induced dynamics on two characteristic spectroscopic features: a shape resonance and multichannel interference. This quantum-path-resolved spectroscopy thus opens wide prospects for the investigation of ultrafast ionic dynamics, such as charge migration.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resonance ionization of zirconium

Abstract Measuring the isotopic composition of trace Zr in presolar stardust grains allows us to study the environment of slow neutron-capture nucleosynthesis in asymptotic giant branch stars. Here, we present a newly characterized Zr resonance ionization scheme that can be saturated with state-of-the-art titanium-sapphire lasers and yields a useful yield of (5.4 ± 0.4)%. Resonance ionization is achieved in two steps: after first being excited with a photon at 319.215 nm, neutral Zr atoms are then ionized with a photon at 388.988 nm, where both wavelengths are reported as in vacuum.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins (0.5-1, 1-4, and 13-16Ry) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors ≳10 times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗

Temporal evolution of the light emitted by a thin, laser-ionized plasma source

We present an experimental and simulation-based investigation of the temporal evolution of light emission from a thin, laser-ionized helium plasma source. We demonstrate an analytic model to calculate the approximate scaling of the time-integrated, on-axis light emission with the initial plasma density and temperature, supported by the experiment, which enhances the understanding of plasma light measurement for plasma wakefield accelerator (PWFA) plasma sources. Our model simulates the plasma density and temperature using a split-step Fourier code and a particle-in-cell code. A fluid simulation is then used to model the plasma and neutral density, and the electron temperature as a function of time and position. We then show the numerical results of the space-and-time-resolved light emission and that collisional excitation is the dominant source of light emission. We validate our model by measuring the light emitted by a laser-ionized plasma using a novel statistical method capable of resolving the nanosecond-scale temporal dynamics of the plasma light using a cost-effective camera with microsecond-scale timing jitter. This method is ideal for deployment in the high radiation environment of a particle accelerator that precludes the use of expensive nanosecond-gated cameras. Our results show that our models can effectively simulate the dynamics of a thin, laser-ionized plasma source. In addition, this work provides a detailed understanding of the plasma light measurement, which is one of the few diagnostic signals available for the direct measurement of PWFA plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions computed by Cloudy in the presence of a generalized incident local radiation field. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors $\gtrsim 10$ times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics of Product Branching from Radical Intermediates in Elementary Combustion Reactions (Final Technical Report)

While the total rate constant for many elementary reactions is well-characterized, understanding the product branching in complex reactions presents a formidable challenge. To gain an incisive probe of such reactions, our experiments investigate the dynamics of the product channels that arise from transient radical intermediates along the bimolecular reaction coordinates. Our recent published work uses the methodology developed in my group in the last fifteen years, using both imaging and scattering apparatuses. The experiments generate a particular isomeric form of an unstable radical intermediate along a bimolecular reaction coordinate and study the branching between the ensuing product channels of the energized radical as a function of its internal rotational and vibrational energy under collision-less conditions. The experiments use a combination of: 1) measurement of product velocity and angular distributions in a crossed laser-molecular beam apparatus, with electron bombardment detection in my lab in Chicago or 2) with tunable vacuum ultraviolet photoionization detection at Taiwan's National Synchrotron Radiation Research Center (NSRRC), and 3) velocity map imaging using state-selective REMPI and single photon VUV ionization of radical intermediates and reaction products. We also employ tunable VUV photoionization detection in our imaging apparatus, using difference frequency four-wave mixing to produce photoionization light tunable from 8 to 10.8 eV. This year of support funded our work to photolytically generate the radical intermediate important in the OH + propene reaction when OH adds to an end carbon. Our other major result this past year was on the photoionization detection of vinoxy radicals. Our experiments at the NSRRC definitively show that vinoxy dissociatively ionizes to both CH 3 + and HCO + , and establishes the absolute partial photoionization cross section of vinoxy to CH 3 + at 10.5 and 11.44 eV. This should allow experimentalists detecting vinoxy radicals produced in bimolecular reactions such as OH + propene to get good branching fractions to that product channel. Finally we initiated studies of the photodissociation of 1-bromo-2-chloroethane and 1,1-bromochloroethane at 193 nm at the NSRRC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic performance of Low-Gain Avalanche Diodes for different types of radiation

LGADs (Low-Gain Avalanche Diodes or Detectors) are a type of silicon Avalanche Photo-Diodes originally developed for the fast detection of minimum ionizing particles in high-energy physics experiments. Thanks to their fast timing performance, the LGAD paradigm enables detectors to accurately measure minimum ionizing particles with a timing resolution of a few tens of picoseconds. Such a performance is due to a thin substrate and the presence of a moderate signal gain. This internal gain of a few tens is enough to compensate for the reduced charge deposition in the thinner substrate and the noise of fast read-out systems. While LGADs are optimized for the detection of minimum ionizing particles for high-energy particle detectors, it is critical to study their performance for the detection of different types of particle, such as X-rays, gamma-rays, or alphas. In this paper, we evaluate the gain of three types of LGADs: two devices with different geometries and doping profiles fabricated by Brookhaven National Laboratory, and one fabricated by Hamamatsu Photonics with a different process. Since the gain in LGADs depends on the bias voltage applied to the sensor, pulse-height spectra have been acquired for bias voltages spanning from the depletion voltage up to the breakdown voltage. Finally, the signal-to-noise ratio of the generated signals and the shape of their spectra allow us to probe the underlying physics of the multiplication process.

47 OTHER INSTRUMENTATION↗

A 100-Mrad (Si) JFET-based Sensing and Communications System for Extreme Nuclear Instrumentation Environments

Dry cask storage is one of two storage methods approved by the U.S. Nuclear Regulatory Commission for spent fuel after removal from reactor cores. Dry casks consist of a stainless steel canister enclosed in a concrete overpack to contain the hazardous radioactive spent fuel rods and provide radiation shielding. Monitoring spent fuel storage casks is desired to ensure the safe containment of the enclosed spent fuel, but is very difficult due to the related harsh temperature and radiation environment. The sensors and associated electronics to monitor temperature, pressure, and/or radiation need to survive high temperatures and radiation doses for extended time periods. For this reason, there is a severe need for radiation-hardened electrical systems that survive well beyond the existing capabilities of commercially available radiation-rated electronic components, which have primarily been developed for space applications. Junction-gate field-effect transistor (JFET) devices are inherently radiation hardened [exceeding 100 Mrad (Si)]. Additionally, when JFETs are used as building blocks for sensing and communication electronics (i.e., oscillators, amplifiers, filters, and mixers), inherently radiation-hardened circuits can be achieved. To this end, JFET-based radiation-hardened electronics interfacing with cask-embedded sensors capable of driving modulated sensor signals through a stainless steel barrier were designed and tested at a dose rate of approximately 500 krad/h (Si) to beyond a 200-Mrad (Si) total ionizing dose. After 200 Mrad (Si), the sensor and communication circuit signals were correctly decoded at the receiver despite oscillator drift. The results from this experiment demonstrate the potential for creating more complex radiation-hardened JFET-based electrical systems for nuclear environments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of an Enhanced Radiation Physics Toolset for Modeling Spectral and Imaging Signatures in the Warm Dense Matter Experiments

Radiative and atomic processes in plasmas play a critical role in a wide variety of high energy density laboratory plasma (HEDLP) experiments. The emission, absorption, and transport of radiation can strongly affect the overall energetics and evolution of such plasmas. In addition, radiation-based diagnostics – including imaging, spectroscopy, and absolute flux measurements – are widely used to determine key features of HEDLPs. To advance our understanding of HEDLP science, it is vital to have high-fidelity computational physics tools that have well-tested radiation physics modeling, and that are readily accessible to researchers in the HEDLP community. Simulations play an extremely important role for planning and designing the experiments, as well as for post-experiment data analysis. Prism Computational Sciences develops software that is used by National Laboratories and universities (including five members of LaserNetUS network). Prominent examples of such research efforts include z-pinch and short-pulse laser experiments designed to study the basic physics of photoionized plasmas and photoionization fronts, as well as their application to astrophysical plasmas. The main effort was dedicated to the development of non-equilibrium equation-of-state (EOS) models within the HELIOS-CR code, a hydrodynamics code with inline collisional-radiative atomic kinetics. Gas cell experiments on Z and Omega demonstrated the importance of non-equilibrium effects on atomic kinetics in photoionized plasmas. Recent proof-of-principle experiments on Omega EP confirmed the advantages of using a short-pulse laser to create an intense radiation drive, leading to additional experiments being proposed. Photoionization front experiments at LLE also emphasizes the importance of radiation and atomic physics. In both studies, HELIOS-CR simulations played a crucial role in computing non-equilibrium opacities and ionization distributions. A newly developed non-LTE EOS model will help addressing possible non-equilibrium effects, on for example specific heat, and their influence in plasma evolution. Prism also implemented support for open-source atomic data generated by the Flexible Atomic Code. This allows researchers to generate custom atomic tables and use them within the complex framework of simulation tools developed by Prism. The ability to use open-source atomic data would be extremely valuable for hydrodynamics and spectroscopic simulations that include high-Z materials, e.g., picosecond x-ray pulse generation experiments. Support for new atomic structures was fully implemented, and the data can be used by all simulation tools developed at Prism: radiation-hydrodynamics, imaging and spectroscopy, EOS and opacity. The development resulted in a significant fidelity enhancement to the simulations tools developed by Prism that are currently used in other cutting-edge experiments including: opacity measurement experiments performed to both understand the basic radiative and atomic properties of plasmas as well as provide data for more accurately modeling the internal structure of the Sun and other stars, high-intensity short-pulse laser experiments performed to develop short-wavelength light sources for use as backlighters and to investigate fast ignition concepts for inertial fusion energy; capsule implosion experiments designed to develop inertial fusion as an energy source, etc.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Chemical abrasion: the mechanics of zircon dissolution

Chemical abrasion is a technique that combines thermal annealing and partial dissolution in hydrofluoric acid (HF) to selectively remove radiation-damaged portions of zircon crystals prior to U–Pb isotopic analysis, and it is applied ubiquitously to zircon prior to U–Pb isotope dilution thermal ionization mass spectrometry (ID-TIMS). The mechanics of zircon dissolution in HF and the impact of different leaching conditions on the zircon structure, however, are poorly resolved. We present a microstructural investigation that integrates microscale X-ray computed tomography (µCT), scanning electron microscopy, and Raman spectroscopy to evaluate zircon dissolution in HF. We show that µCT is an effective tool for imaging metamictization and complex dissolution networks in three dimensions. Acid frequently reaches crystal interiors via fractures spatially associated with radiation damage zoning and inclusions to dissolve soluble high-U zones, some inclusions, and material around fractures, leaving behind a more crystalline zircon residue. Other acid paths to crystal cores include the dissolution of surface-reaching inclusions and the percolation of acid across zones with high defect densities. In highly crystalline samples dissolution is crystallographically controlled with dissolution proceeding almost exclusively along the c axis. Increasing the leaching temperature from 180 to 210 °C results in deeper etching textures, wider acid paths, more complex internal dissolution networks, and greater volume losses. How a grain dissolves strongly depends on its initial radiation damage content and defect distribution as well as the size and position of inclusions. As such, the effectiveness of any chemical abrasion protocol for ID-TIMS U–Pb geochronology is likely sample-dependent. We also briefly discuss the implications of our findings for deep-time (U-Th)/He thermochronology.

58 GEOSCIENCES↗