Search NASA⌕ Search

SEARCH · Search NASA

Results for “Li diffusion”

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 163 records · Page 9

Model Stellar Atmospheres and Real Stellar Atmospheres and Status of the ATLAS12 Opacity Sampling Program and of New Programs for Rosseland and for Distribution Function Opacity

I discuss errors in theory and in interpreting observations that are produced by the failure to consider resolution in space, time, and energy. I discuss convection in stellar model atmospheres and in stars. Large errors in abundances are possible such as the factor of ten error in the Li abundance for extreme Population II stars. Finally I discuss the variation of microturbulent velocity with depth, effective temperature, gravity, and abundance. These variations must be dealt with in computing models and grids and in any type of photometric calibration. I have also developed a new opacity-sampling version of my model atmosphere program called ATLAS12. It recognizes more than 1000 atomic and molecular species, each in up to 10 isotopic forms. It can treat all ions of the elements up through Zn and the first 5 ions of heavier elements up through Es. The elemental and isotopic abundances are treated as variables with depth. The fluxes predicted by ATLAS12 are not accurate in intermediate or narrow bandpass intervals because the sample size is too small. A special stripped version of the spectrum synthesis program SYNTHE is used to generate the surface flux for the converged model using the line data on CD-ROMs 1 and 15. ATLAS12 can be used to produce improved models for Am and Ap stars. It should be very useful for investigating diffusion effects in atmospheres. It can be used to model exciting stars for H II regions with abundances consistent with those of the H II region. These programs and line files will be distributed on CD-ROMs.

Kurucz, Robert L.↗

Induced junction solar cell and method of fabrication

An induced junction solar cell is fabricated on a p-type silicon substrate by first diffusing a grid of criss-crossed current collecting n+ stripes and thermally growing a thin SiO2 film, and then, using silicon-rich chemical vapor deposition (CVD), producing a layer of SiO2 having inherent defects, such as silicon interstices, which function as deep traps for spontaneous positive charges. Ion implantation increases the stable positive charge distribution for a greater inversion layer in the p-type silicon near the surface. After etching through the oxide to parallel collecting stripes, a pattern of metal is produced consisting of a set of contact stripes over the exposed collecting stripes and a diamond shaped pattern which functions as a current collection bus. Then the reverse side is metallized.

Maserjian, J.↗

Lower-Hybrid Wave Structures and Interactions With Electrons Observed in Magnetotail Reconnection Diffusion Regions

We investigate waves close to the lower-hybrid frequency in 12 magnetotail reconnection electron diffusion region (EDR) events with guide field levels of near-zero to 30%. In about half of the events, the wave vector has a small component along the current sheet normal, consistent with known lower-hybrid drift wave properties, but the perpendicular magnetic field fluctuations can be comparable or greater than the parallel component, a feature unique to the waves inside and adjacent to EDRs. Another new wave property is that the wave vector has a significant component along the current sheet normal in some events and completely along the normal for one event. In 1/4of the events, the ∇∙𝑷(sub 𝑒) term has a significant contribution to the wave electric field, possibly a feature of lower-hybrid waves more likely to exist in the diffusion region than further away from the X-line. Electron temperature variations are correlated with the wave potential, due to wave electric field acceleration and crossings at the corrugated separatrix region with different amounts of mixing between reconnection inflowing and outflowing populations. The latter also leads to the anti-correlation between parallel and perpendicular temperature components. Using four-spacecraft measurements, the magnetic field line twisting is demonstrated by the correlated fluctuations in (∇×𝑽(sub 𝐸×𝐵))||and (∇×𝐁)||. The lower-hybrid wave in the EDR of weak guide field reconnection may be generated near separatrices and penetrate to the mid-plane or locally generated, and the latter possibility is beyond the prediction of previous reconnection simulations.

Magnetic reconnection↗

Erasing the Milky Way: New Cleaning Technique Applied to GBT Intensity Mapping Data

We present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15- and 1-h field data of the GBT observations previously presented in Masui et al. (2013) and Switzer et al. (2013), covering about 41 square degrees at 0.6 less than z is less than 1.0, for which cross-correlations may be measured with the galaxy distribution of the WiggleZ Dark Energy Survey. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contamination using an independent component analysis technique (fastica), and develop a Fourier-based optimal estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission through the non-Gaussian nature of their probability distributions. The temperature power spectra of the intensity maps is dominated by instrumental noise on small scales which fastica, as a conservative sub-traction technique of non-Gaussian signals, can not mitigate. However, we determine similar GBT-WiggleZ cross-correlation measurements to those obtained by the Singular Value Decomposition (SVD) method, and confirm that foreground subtraction with fastica is robust against 21cm signal loss, as seen by the converged amplitude of these cross-correlation measurements. We conclude that SVD and fastica are complementary methods to investigate the foregrounds and noise systematics present in intensity mapping datasets.

Wolz, L.↗

GJ 3929: High-precision Photometric and Doppler Characterization of an Exo-Venus and Its Hot, Mini-Neptune-mass Companion

We detail the follow-up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b (TOI-2013b), and its nontransiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star’s Venus zone ( P b = 2.616272 ± 0.000005 days; S b = 17.3 +0.8 -0.7 S ⊕ ) orbiting a nearby M dwarf. GJ 3929c is most likely a nontransiting sub-Neptune. Using the new, ultraprecise NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak National Observatory, we are able to modify the mass constraints of planet b reported in previous works and consequently improve the significance of the mass measurement to almost 4σ confidence ( M b = 1.75 ± 0.45 M ⊕ ). We further adjust the orbital period of planet c from its alias at 14.30 ± 0.03 days to the likely true period of 15.04 ± 0.03 days, and we adjust its minimum mass to m sin i = 5.71 ± 0.92 M ⊕ . Using the diffuser-assisted ARCTIC imager on the ARC 3.5 m telescope at Apache Point Observatory, in addition to publicly available TESS and LCOGT photometry, we are able to constrain the radius of planet b to R p = 1.09 ± 0.04 R ⊕ . GJ 3929b is a top candidate for transmission spectroscopy in its size regime (TSM = 14 ± 4), and future atmospheric studies of GJ 3929b stand to shed light on the nature of small planets orbiting M dwarfs.

Exoplanets↗

Dynamics of Disorder-Order Transitions in Hard Sphere Colloidal Dispersions in micro-g

We performed a series of experiments on 0.518 millimeter PMMA spheres suspended in an index matching mixture of decalin and tetralin the microgravity environment provided by the Shuttle Columbia on mission STS-73. The samples ranged in concentration from 0.49 to 0.62. volume fraction (phi) of spheres, which covers the range in which liquid, coexistence, solid and glass phases are expected from Earth bound experiments. Light scattering was used to probe the static structure, and the particle dynamics. Digital and 35 mm photos provided information on the morphology of the crystals. In general, the crystallites grew considerably larger (roughly an order of magnitude larger) than the same samples with identical treatment in 1 g. The dynamic light scattering shows the typical short time diffusion and long time caging effects found in 1 g. The surprises that were encountered in microgravity include the preponderance of random hexagonal close packed (RHCP) structures and the complete absence of the expected face centered cubic (FCC) structure, existence of large dendritic crystals floating in the coexistence samples (where liquid and solid phases coexist) and the rapid crystallization of samples which exist only in glass phase under the influence of one g. These results suggest that colloidal crystal growth is profoundly effected by gravity in yet unrecognized ways. We suspect that the RCHP structure is related to the nonequilibrium growth that is evident from the presence of dendrites. An analysis of the dendritic growth instabilities is presented within the framework of the Ackerson-Schatzel equation.

Zhu, J. X.↗

Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple

The development of more efficient thermoelectric couple technology capable of operating with high-grade heat sources up to 1,275 K is key to improving the performance of radioisotope thermoelectric generators. Lanthanum telluride La3-xTe4 and 14-1-11 Zintls (Yb14MnSb11) have been identified as very promising materials. The fabrication of advanced high-temperature thermoelectric couples requires the joining of several dissimilar materials, typically including a number of diffusion bonding and brazing steps, to achieve a device capable of operating at elevated temperatures across a large temperature differential (up to 900 K). A thermoelectric couple typically comprises a heat collector/ exchanger, metallic interconnects on both hot and cold sides, n-type and ptype conductivity thermoelectric elements, and cold-side hardware to connect to the cold-side heat rejection and provide electrical connections. Differences in the physical, mechanical, and chemical properties of the materials that make up the thermoelectric couple, especially differences in the coefficients of thermal expansion (CTE), result in undesirable interfacial stresses that can lead to mechanical failure of the device. The problem is further complicated by the fact that the thermoelectric materials under consideration have large CTE values, are brittle, and cracks can propagate through them with minimal resistance. The inherent challenge of bonding brittle, high-thermal-expansion thermoelectric materials to a hot shoe material that is thick enough to carry the requisite electrical current was overcome. A critical advantage over prior art is that this device was constructed using all diffusion bonds and a minimum number of assembly steps. The fabrication process and the materials used are described in the following steps: (1) Applying a thin refractory metal foil to both sides of lanthanum telluride. To fabricate the n-type leg of the advanced thermoelectric couple, the pre-synthesized lanthanum telluride coupon was diffusion bonded to the metal foil using a thin adhesion layer. (2) Repeating a similar process for the 14-1-11 Zintl p-type leg of the advanced thermoelectric couple. (3) Bonding thick CTE-matched metal plates on the metallized lanthanum telluride and Yb14MnSb11 to form the hot and cold sides of the thermoelectric couple. The calculated conversion efficiency of such an advanced couple would be about 10.5 percent, about 35 percent better than heritage radioisotope thermoelectric technology that relies on Si-Ge alloys. In addition, unlike Si-Ge alloys, these materials can be combined with many other thermoelectric materials optimized for operation at lower temperatures to achieve conversion efficiency in excess of 15 percent (a factor of 2 increase over heritage technology).

Ravi, Vilupanur A.↗

Ion Behaviors in the Reconnection Diffusion Region of a Corrugated Magnetotail Current Sheet

The current sheet structure and ion behaviors in a magnetotail reconnection diffusion regionare investigated. The multispacecraft analysis suggests a corrugated current sheet structure, interpretedas due to a apping motion that propagates along geocentric solar magnetospheric along the +y direction inthe Geocentric Solar Magnetospheric (GSM) coordinate. The electric eld (E) and ion distributions havesimilarities with those in a planar current sheet. Energetic ions move along the current direction, suggestingthe acceleration by the observed reconnection E during the meandering motion. Counterstreaming ionsalong the current sheet normal suggest the acceleration by the Hall E that is observed to be the dominantcomponent. However, at certain locations, E and counterstreaming ions signicantly deviate from the localnormal direction, and more than one pair of counterstreaming populations exist, possibly because thecorrugated current sheet enables ions entering the current sheet at different locations with differentvelocities to mix together

Wang, Shan↗

Thermal Diffusivity for III-VI Semiconductor Melts at Different Temperatures

The change of the thermal properties of semiconductor melts reflects the structural changes inside the melts, and a fundamental understanding of this structural transformation is essential for high quality semiconductor crystal growth process. This paper focused on the technical development and the measurement of thermal properties of III-VI semiconductor melts at high temperatures. Our previous work has improved the laser flash method for the specialized quartz sample cell. In this paper, we reported the results of our recent progress in further improvements of the measurement system by minimizing the free convection of the melt, adding a front IR detector, and placing the sample cell in a vacuum environment. The results for tellurium and selenium based compounds, some of which have never been reported in the literature, were obtained at different temperatures as a function of time. The data were compared with other measured thermophysical properties to shed light on the structural transformations of the melt.

Ban, H.↗

Normal freezing of ideal ternary systems of the pseudobinary type

Perfect liquid mixing but no solid diffusion is assumed in normal freezing. In addition, the molar compositions of the freezing solid and remaining liquid, respectively, follow the solidus and liquidus curves of the constitutional diagram. For the linear case, in which both the liquidus and solidus are perfectly straight lines, the normal freezing equation giving the fraction solidified at each melt temperature and the solute concentration profile in the frozen solid was determined as early as 1902, and has since been repeatedly published. Corresponding equations for quadratic, cubic or higher-degree liquidus and solidus lines have also been obtained. The equation of normal freezing for ideal ternary liquid solutions solidified into ideal solid solutions of the pseudobinary type is given. Sample computations with the use of this new equation were made and are given for the Ga-Al-As system.

Li, C. H.↗

On the Electron Diffusion Region in Asymmetric Reconnection with a Guide Magnetic Field

Particle-in-cell simulations in a 2.5-D geometry and analytical theory are employed to study the electron diffusion region in asymmetric reconnection with a guide magnetic field. The analysis presented here demonstrates that similar to the case without guide field, in-plane flow stagnation and null of the in-plane magnetic field are well separated. In addition, it is shown that the electric field at the local magnetic X point is again dominated by inertial effects, whereas it remains dominated by nongyrotropic pressure effects at the in-plane flow stagnation point. A comparison between local electron Larmor radii and the magnetic gradient scale lengths predicts that distribution should become nongyrotropic in a region enveloping both field reversal and flow stagnation points. This prediction is verified by an analysis of modeled electron distributions, which show clear evidence of mixing in the critical region.

reconnection↗

A study of plasmaspheric density distributions for diffusive equilibrium conditions

The plasmaspheric density distribution has been modeled for a range of solar cycle, seasonal and diurnal conditions with a magnetic flux tube dependent diffusive equilibrium model by using experimentally determined values of ionospheric parameters at 675 km as boundary conditions. Data is presented in terms of plasmaspheric H(+) and He(+) density contours, total flux tube content, and equatorial plasma density for a range of L-values from 1.15 to 3.0. The variation of equatorial density with L-value shows good agreement with the 1/L exp 4 dependence observed experimentally. The results show that the model predicts larger solar cycle and diurnal variation in equatorial plasma density than observed using whistler techniques. However, the whistler method requires a model to deduce the equatorial density and is therefore open to interpretation. Seasonal variations are rather artificial since in this general model no attempt has been made to match equatorial densities for flux tubes emanating from the winter and summer hemispheres.

Li, W.↗

Cold Ion Demagnetization near the X-line of Magnetic Reconnection

Although the effects of magnetic reconnection in magnetospheres can be observed at planetary scales, reconnection is initiated at electron scales in a plasma. Surrounding the electron diffusion region, there is an Ion-Decoupling Region (IDR) of the size of the ion length scales (inertial length and gyroradius). Reconnection at the Earths magnetopause often includes cold magnetospheric (few tens of eV), hot magnetospheric (10 keV), and magnetosheath (1 keV) ions, with different gyroradius length scales. We report observations of a subregion inside the IDR of the size of the cold ion population gyroradius (approx. 15 km) where the cold ions are demagnetized and accelerated parallel to the Hall electric field. Outside the subregion, cold ions follow the E x B motion together with electrons, while hot ions are demagnetized. We observe a sharp cold ion density gradient separating the two regions, which we identify as the cold and hot IDRs.

Toledo-Redondo, Serio↗

Energetic Particles of Cosmic Accelerators I: Galactic Accelerators

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Our pursuit of more than a century to uncover the origins and fate of these cosmic energetic particles has given rise to some of the most interesting and challenging questions in astrophysics. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. For instance, GCRs act in concert with galactic magnetic fields to support its disk against its own weight. GCR ionization and heating are essential ingredients in promoting and regulating the formation of stars and protostellar disks. GCR ionization also drives astrochemistry, leading to the build up of complex molecules in the interstellar medium. GCR transport throughout the galaxy generates and maintains turbulence in the interstellar medium, alters its multi-phase structure, and amplifies magnetic fields. GCRs could even launch galactic winds that enrich the circumgalactic medium and alter the structure and evolution of galactic disks. As crucial as they are for many of the varied phenomena in our galaxy, there is still much we do not understand about GCRs. While they have been linked to supernova remnants (SNRs), it remains unclear whether these objects can fully account for their entire population, particularly at the lower (approximately less than 1 GeV per nucleon) and higher (~PeV) ends of the spectrum. In fact, it is entirely possible that the SNRs that have been found to accelerate CRs merely re-accelerate them, leaving the origins of the original GCRs a mystery. The conditions for particle acceleration that make SNRs compelling source candidates are also likely to be present in sources such as protostellar jets, superbubbles, and colliding wind binaries (CWBs), but we have yet to ascertain their roles in producing GCRs. For that matter, key details of diffusive shock acceleration (DSA) have yet to be revealed, and it remains to be seen whether DSA can adequately explain particle acceleration in the cosmos. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. For the case of GCRs, MeV astronomy will: 1) Search for fresh acceleration of GCRs in SNRs; 2) Test the DSA process, particularly in SNRs and CWBs; 3) Search for signs of CR acceleration in protostellar jets and superbubbles.

Venters, Tonia M.↗

Mineralogical Composition of the Different Types of Bright Deposits on Vesta

VIR-MS, Dawn's Visible and Infrared Mapping Spectrometer, obtained hyperspectral images of a wide part of Vesta's surface at a variety of spatial resolutions [1]. Vesta spectra are similar to those of the howardite-eucrite-diogenite (HED) meteorites. Moreover, they are characterized by the two iron-bearing pyroxene bands at 0.9 (band I) and 1.9 microns (band II). Vesta surface's is dominated by eucrite/howardite with some diogenitic regions situated in the southern hemisphere near the Rheasilvia basin [2]. The surface is heavily craterized and impacts can expose fresh material, thus generating the Bright Material Deposits (BMD) observed within and surrounding certain craters. BMD can be classified into six different types based on their morphological characteristics: Crater Wall/Scarp Material (CWM), Radial Material (RM), Slope Material (SM), Patchy Material (PM), Spot Material (SpM) and Diffuse Plains Material (DPM) [3]. The most widespread BMD are CWM, SM and RM. CWM, SM, RM originate from impacts. CWM is situated on the edge of the craters. Mass wasting from the crater walls and generates the SM, while RM is associated with the ejecta of the craters [4]. BMD are characterized by albedo greater than that of the vestan average, 0.38 [5]. Therefore the different types of deposits present distinct levels of reflectance respect to the Surrounding Regions (SR), in particular: the CWM and SM is approx.40% brighter, the RM is approx.30- 40% brighter; the SpM is about 20-25% brighter and the PM is about 20% brighter. Near the edge of the Rheasilvia basin it is possible to find some extremely bright areas ~80% brighter than the vestan average [6].

Zambon, F.↗

Photometric and Spectroscopic Properties of Type Ia Supernova 2018oh with Early Excess Emission from the Kepler 2 Observations

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically confirmed Type Ia supernova (SN Ia) observed in the Kepler field. The Kepler data revealed an excess emission in its early light curve, allowing us to place interesting constraints on its progenitor system. Here we present extensive optical, ultraviolet, and nearinfrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3±0.3 days and Δ(m15)(B)=0.96±0.03 mag, but it seems to have bluer B−V colors. We construct the “UVOIR” bolometric light curve having a peak luminosity of 1.49×10(Exp 43) erg/s, from which we derive a nickel mass as 0.55±0.04M(ʘ) by fitting radiation diffusion models powered by centrally located 56Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of 56Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a nondegenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia but is characterized by prominent and persistent carbon absorption features. The CII features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in an SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

Thomas Barclay↗

Benefits of a 4th Ice Class in the Simulated Radar Reflectivities of Convective Systems Using a Bulk Microphysics Scheme

Numerous cloud microphysical schemes designed for cloud and mesoscale models are currently in use, ranging from simple bulk to multi-moment, multi-class to explicit bin schemes. This study details the benefits of adding a 4th ice class (hail) to an already improved 3-class ice bulk microphysics scheme developed for the Goddard Cumulus Ensemble model based on Rutledge and Hobbs (1983,1984). Besides the addition and modification of several hail processes from Lin et al. (1983), further modifications were made to the 3-ice processes, including allowing greater ice super saturation and mitigating spurious evaporationsublimation in the saturation adjustment scheme, allowing graupelhail to become snow via vapor growth and hail to become graupel via riming, and the inclusion of a rain evaporation correction and vapor diffusivity factor. The improved 3-ice snowgraupel size-mapping schemes were adjusted to be more stable at higher mixing rations and to increase the aggregation effect for snow. A snow density mapping was also added. The new scheme was applied to an intense continental squall line and a weaker, loosely-organized continental case using three different hail intercepts. Peak simulated reflectivities agree well with radar for both the intense and weaker case and were better than earlier 3-ice versions when using a moderate and large intercept for hail, respectively. Simulated reflectivity distributions versus height were also improved versus radar in both cases compared to earlier 3-ice versions. The bin-based rain evaporation correction affected the squall line case more but did not change the overall agreement in reflectivity distributions.

bulk microphysics schemes↗

SAO and Kelvin Waves in the EuroGRIPS GCMS and the UK Meteorological Offices Analyses

This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.

Amodei, M.↗