Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “IONIZATION”

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 541 records · Page 30

Spectra of High-Ionization Seyfert 1 Galaxies: Implications for the Narrow-Line Region

We present line profiles and profile parameters for the Narrow-Line Regions (NLRs) of six Seyfert 1 galaxies with high-ionization lines: MCG 8-11-11, Mrk 79, Mrk 704, Mrk 841, NGC 4151, and NGC 5548. The sample was chosen primarily with the goal of obtaining high-quality [Fe VII] lambda6087 and, when possible, [Fe X] lambda6374 profiles to determine if these lines are more likely formed in a physically distinct 'coronal line region' or are formed throughout the NLR along with lines of lower critical density (n(sub cr)) and/or Ionization Potential (IP). We discuss correlations of velocity shift and width with n(sub cr) and IP. In some objects, lines of high IP and/or n(sub cr) are systematically broader than those of low IP/n(sub cr). Of particular interest, however, are objects that show no correlations of line width with either IP or n(sub cr). In these objects, lines of high and low IP/n(sub cr), are remarkably similar, which is difficult to reconcile with the classical picture of the NLR, in which lines of high and low IP/n(sub cr) are formed in physically distinct regions. We argue for similar spatial extents for the flux in lines with similar profiles. Here, as well as in a modeling-oriented companion paper, we develop further an idea suggested by Moore & Cohen that objects that do and do not show line width correlations with IP/n(sub cr) can both be explained in terms of a single NLR model with only a small difference in the cloud column density distinguishing the two types of object. Overall, our objects do not show correlations between the Full Width at Half-Maximum (FWHM) and IP and/or n(sub cr). The width must be defined by a parameter that is sensitive to extended profile wings in order for the correlations to result. We present models in which FWHM correlations with IP and/or n(sub cr) result only after simulating the lower spectral resolution used in previous observational studies. The models that simulate the higher spectral resolution of our observational study produce line width correlations only if the width is defined by a parameter that is more sensitive to extended profile wings than is the FWHM. Our sample of six objects is in effect augmented by incorporating the larger sample (16 objects) of Veilleux into some of our discussion. This paper focuses on new interpretations of NLR emission-line spectra and line profiles that stem directly from the observations. Paper 2 focuses on modeling and complements this paper by illustrating explicitly the effects that spatial variations in electron density, ionization parameter, and column density have on model profiles. By comparing model profiles with the observed profiles presented here, as well as with those presented by Veilleux, Paper 2 yields insight into how the electron density, ionization parameter, and column density likely vary throughout the NLR.

Moore, David↗

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗

Semiclassical Study of the Wave Vector Dependence of the Interband Impact Ionization Rate in Bulk Silicon

We present calculations of the interband impact ionization rate calculated using a wave vector dependent (k-dependent) semiclassical formulation of the transition rate. The transition rate is determined using Fermi's golden rule from a two-body screened Coulomb interaction assuming energy and momentum conservation. The transition rate is calculated for the first two conduction bands of silicon by numerically integrating over the full Brillouin zone. The overlap integrals in the expression for the transition rate are determined numerically using a 15 band k-p calculation. It is found that the transition rate depends strongly on the initiating electron wave vector (k vector) and that the transition rate is greatest for electrons originating within the second conduction band than the first conduction band. An ensemble Monte Carlo simulation, which includes the numerically determined ionization transition rate as well as the full details of the first two conduction bands, is used to calculate the total impact ionization rate in bulk silicon. Good agreement with the experimentally determined electron ionization rate data is obtained.

Wang, Yang↗

Numerical Study of the Wave-Vector Dependence of the Electron Interband Impact Ionization Rate in Bulk GaAs

Ensemble Monte Carlo calculations of the electron interband impact ionization rate in bulk GaAs are presented using a wave-vector (k)-dependent formulation of the ionization transition rate. The transition rate is evaluated through use of numerically generated wavefunctions determined via a k-p calculation within the first two conduction bonds at numerous points within a finely spaced three-dimensional grid in k space. The transition rate is determined to be greatest for states within the second conduction band. It is found that the interband impact ionization transition rate in bulk GaAs is best characterized as having an exceedingly "soft" threshold energy. As a consequence, the dead space, defined as the distance over which the ionization probability for a given carrier is assumed to be zero, is estimated to be-much larger than that estimated using a "harder" threshold. These results have importance in the design of multiquantum-well avalanche photodiodes.

Wang, Yang↗

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗

Viscosity Coefficient Curve Fits for Ionized Gas Species Grant Palmer

Viscosity coefficient curve fits for neutral gas species are available from many sources. Many do a good job of reproducing experimental and computational chemistry data. The curve fits are usually expressed as a function of temperature only. This is consistent with the governing equations used to derive an expression for the neutral species viscosity coefficient. Ionized species pose a more complicated problem. They are subject to electrostatic as well as intermolecular forces. The electrostatic forces are affected by a shielding phenomenon where electrons shield the electrostatic forces of positively charged ions beyond a certain distance. The viscosity coefficient for an ionized gas species is a function of both temperature and local electron number density. Currently available curve fits for ionized gas species, such as those presented by Gupta/Yos, are a function of temperature only. What they did was to assume an electron number density. The problem is that the electron number density they assumed was unrealistically high. The purpose of this paper is two-fold. First, the proper expression for determining the viscosity coefficient of an ionized species as a function of both temperature and electron number density will be presented. Then curve fit coefficients will be developed using the more realistic assumption of an equilibrium electron number density. The results will be compared against previous curve fits and against highly accurate computational chemistry data.

Palmer, Grant↗

Hydrodynamic Models of Line-Driven Accretion Disk Winds III: Local Ionization Equilibrium

We present time-dependent numerical hydrodynamic models of line-driven accretion disk winds in cataclysmic variable systems and calculate wind mass-loss rates and terminal velocities. The models are 2.5-dimensional, include an energy balance condition with radiative heating and cooling processes, and includes local ionization equilibrium introducing time dependence and spatial dependence on the line radiation force parameters. The radiation field is assumed to originate in an optically thick accretion disk. Wind ion populations are calculated under the assumption that local ionization equilibrium is determined by photoionization and radiative recombination, similar to a photoionized nebula. We find a steady wind flowing from the accretion disk. Radiative heating tends to maintain the temperature in the higher density wind regions near the disk surface, rather than cooling adiabatically. For a disk luminosity L (sub disk) = solar luminosity, white dwarf mass M(sub wd) = 0.6 solar mass, and white dwarf radii R(sub wd) = 0.01 solar radius, we obtain a wind mass-loss rate of M(sub wind) = 4 x 10(exp -12) solar mass yr(exp -1) and a terminal velocity of approximately 3000 km per second. These results confirm the general velocity and density structures found in our earlier constant ionization equilibrium adiabatic CV wind models. Further we establish here 2.5D numerical models that can be extended to QSO/AGN winds where the local ionization equilibrium will play a crucial role in the overall dynamics.

Pereyra, Nicolas Antonio↗

Dissociative Ionization of Benzene by Electron Impact

We report a theoretical study of the dissociative ionization (DI) of benzene from the low-lying ionization channels. Our approach makes use of the fact that electron motion is much faster than nuclear motion and DI is treated as a two-step process. The first step is electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model. For the unimolecular dissociation step, we study the steepest descent reaction path to the minimum of the ion potential energy surface. The path is used to analyze the probability of unimolecular dissociation and to determine the product distributions. Our analysis of the dissociation products and the thresholds of the productions are compared with the result dissociative photoionization measurements of Feng et al. The partial oscillator strengths from Feng et al. are then used in the iBED cross section calculations.

Huo, Winifred↗

The Ionized Gas and Nuclear Environment in NGC 3783: Variability and Modeling of the 900 ks CHANDRA Spectrum - IV

We present a detailed analysis of the 900 ks spectrum of NGC3783 obtained by Chandra in 2000-2001 (Kaspi et al. 2002). We split the data in various ways to look for time dependent and luminosity dependent spectral variations. This analysis, the measured equivalent widths of a large number of X-ray lines, and our photoionization calculations, lead us to the following conclusions: 1) NGC 3783 fluctuated in luminosity, by a factor N 1.5, during individual 170 ks observations. The fluctuations were not associated with significant spectral variations. 2) On a longer time scale, of 20-120 days, we discovered two very different spectral shapes that are noted the high state and the low state spectra. The observed changes between the two can be described as the appearance and disappearance of a soft continuum component. The spectral variations are not related, in a simple way, to the brightening or the fading of the short wavelength continuum, as observed in other objects. NGC3783 seems to be the first AGN to show this unusual behavior. 3) The appearance of the soft continuum component is consistent with beeing the only spectral variation and there is no need to invoke changes in the absorber s opacity. In particular, all absorption lines with reliable measurements show the same equivalent width, within the observational uncertainties, during high and low states. 4) Photoionization model calculations show that a combination of three ionization components, each split into two kinematic components, explain very well the intensity of almost all absorption lines and the bound-free absorption. The components span a large range of ionization and a total column of about 3 x 10(exp 22) per square centimeter Moreover, all components are thermally stable and are situated on the vertical branch of the stability curve.. This means that they are in pressure equilibrium and perhaps occupy the same volume of space. This is the first detection of such a multi-component equilibrium gas in AGN. 5) The only real discrepancy between the model and the observations is the wavelength of the iron M-shell UTA feature. This is most likely due to an underestimation of the dielectronic recombination O VI and discuss its possible origin. 6) The lower limit on the distance of the absorbing gas in NGC3783 is between 0.2 and 3.2 pc, depending of the specific ionization component. The constant pressure assumption imposes an upper limit of about 25 pc on the distance of the least ionized gas from the central sourec.

Netzer, Hagai↗

Dissociative Ionization and Product Distributions of Benzene and Pyridine by Electron Impact

We report a theoretical study of the dissociative ionization (DI) and product distributions of benzene (C6H6) and pyridine (C5H5N) from their low-lying ionization channels. Our approach makes use of the fact that electronic motion is much faster than nuclear motion allowing DI to be treated as a two-step process. The first step is the electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step, the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model [W.M. Huo, Phys. Rev. A64,042719-I (2001)]. For the unimolecular dissociation, we use multiconfigurational self-consistent field (MCSCF) methods to determine the steepest descent pathways to the possible product channels. More accurate methods are then used to obtain better energetics of the paths which are used to determine unimolecular dissociation probabilities and product distributions. Our analysis of the dissociation products and the thresholds of their productions for benzene are compared with the recent dissociative photoionization meausurements of benzene by Feng et al. [R. Feng, G. Cooper, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,211 (2002)] and the dissociative photoionization measurements of pyridine by Tixier et al. [S. Tixier, G. Cooper, R. Feng, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,185 (2002)] using dipole (e,e+ion) coincidence spectroscopy.

Dateo, Christopher E.↗

Dissociative Ionization of Aromatic and Heterocyclic Molecules

Space radiation poses a major health hazard to humans in space flight. The high-energy charged particles in space radiation ranging from protons to high atomic number, high-energy (HZE) particles, and the secondary species they produce, attack DNA, cells, and tissues. Of the potential hazards, long-term health effects such as carcinogenesis are likely linked to the DNA lesions caused by secondary electrons in the 1 - 30 eV range. Dissociative ionization (DI) is one of the electron collision processes that can damage the DNA, either directly by causing a DNA lesion, or indirectly by producing radicals and cations that attack the DNA. To understand this process, we have developed a theoretical model for DI. Our model makes use of the fact that electron motion is much faster than nuclear motion and assumes DI proceeds through a two-step process. The first step is electron-impact ionization resulting in a particular state of the molecular ion in the geometry of the neutral molecule. In the second step the ion undergoes unimolecular dissociation. Thus the DI cross section sigma(sup DI)(sub a) for channel a is given by sigma(sup DI)(sub a) = sigma(sup I)(sub a) P(sub D) with sigma(sup I)(sub a) the ionization cross section of channel a and P(sub D) the dissociation probability. This model has been applied to study the DI of H2O, NH3, and CH4, with results in good agreement with experiment. The ionization cross section sigma(sup I)(sub a) was calculated using the improved binary encounter-dipole model and the unimolecular dissociation probability P(sub D) obtained by following the minimum energy path determined by the gradients and Hessians of the electronic energy with respect to the nuclear coordinates of the ion. This model is used to study the DI from the low-lying channels of benzene and pyridine to understand the different product formation in aromatic and heterocyclic molecules. DI study of the DNA base thymine is underway. Solvent effects will also be discussed.

Huo, Winifred M.↗

Hard X-Ray Emission and the Ionizing Source in LINERs

We report X-ray luminosities of 21 LINERs (low-ionization nuclear emission-line regions) and 17 low-luminosity Seyferts obtained with ASCA and discuss the ionizing source in LINERs. Most LINERs with broad H-alpha emission in their optical spectra (LINER 1s) have a compact hard X-ray source and their 2-10 keV X-ray luminosities (LX) are proportional to their H alpha luminosities (L-H-alpha). This correlation strongly supports the hypothesis that the dominant ionizing source in LINER 1s is photoionization by hard photons from low-luminosity AGNs. Although some LINERs without broad H-alpha emission (LINER 2s) have X-ray properties similar to LINER 1s, the X-ray luminosities of many LINER 2s in our sample are lower than LINER 1s at a given H-alpha luminosity. The observed X-ray luminosities in these objects are insufficient to power their H-alpha luminosities, suggesting that their primary ionizing source is something other than an AGN, or that an AGN, if present, is obscured even at energies above 2 keV. LINER 2s having small LX/LH-alpha occupy a localized region with small [OI]/H-alpha on the excitation diagram. Such LINER spectra can be reproduced by photoionization by very hot stars.

Terashima, Y.↗

The Properties and the Evolution of the Highly Ionized Gas in MR 2251-178

We present the first XMM-Newton observations of the radio-quiet quasar MR 2251-178 obtained in 2000 and 2002. The EPIC-pn spectra show a power-law continuum with a slope of Gamma = 1.6 at high energies absorbed by at least two warm absorbers (WAs) intrinsic to the source. The underlying continuum in the earlier observation shows a soft excess at low X-ray energies which can be modeled as an additional power-law with Gamma = 2.9. The spectra also show a weak narrow iron K alpha emission line. The high-resolution grating spectrum obtained in 2002 shows emission lines from N VI, O VII, O VIII, Ne IX, and Ne X, as well as absorption lines from the low-ionization ions of O III, O IV, and O V, and other confirmed and suspected weaker absorption lines. The O III - O V lines are consistent with the properties of the emission line gas observed as extended optical (O III) emission in this source. The signal-to-noise of the 2000 grating data is too low to detect any lines. We suggest a model for the high-resolution spectrum which consist of two or three warm-absorber (WA) components. The two-components model has a high-ionization WA with a column density of 10(exp 21.5)-10 (exp 21.8) sq cm and a low-ionization absorber with a column density of 10(exp 20.3) sq cm. In the three-components model we add a lower ionization component that produces the observed iron M-shell absorption lines. We investigate the spectral variations in MR 2251-178 over a period of 8.5 years using data from ASCA, BeppoSAX, and XMM-Newton. All X-ray observations can be fitted with the above two power laws and the two absorbers. The observed luminosity variations seems to correlate with variations in the soft X-ray continuum. The 8.5 year history of the source suggests a changing X-ray absorber due to material that enters and disappears from the line-of-sight on timescales of several months. We also present, for the first time, the entire FUSE spectrum of MR 2251-178. We detect emission from N III, C III, and O VI and at least 4 absorption systems in C III, H I, and O VI, one at -580 km/s and at least 3 others which are blended together and form a wide trough covering the velocity range of 0 to -500 km/s. The general characteristics of the UV and X-ray absorbers are consistent with an origin in the same gas.

Kaspi, Shai↗

Density, Velocity and Ionization Structure in Accretion-Disc Winds

This was a project to exploit the unique capabilities of FUSE to monitor variations in the wind- formed spectral lines of the luminous, low-inclination, cataclysmic variables(CV) -- RW Sex. (The original proposal contained two additional objects but these were not approved.) These observations were intended to allow us to determine the relative roles of density and ionization state changes in the outflow and to search for spectroscopic signatures of stochastic small-scale structure and shocked gas. By monitoring the temporal behavior of blue-ward extended absorption lines with a wide range of ionization potentials and excitation energies, we proposed to track the changing physical conditions in the outflow. We planned to use a new Monte Carlo code to calculate the ionization structure of and radiative transfer through the CV wind. The analysis therefore was intended to establish the wind geometry, kinematics and ionization state, both in a time-averaged sense and as a function of time.

Sonneborn, George↗

Spectroscopy of Jet-Cooled Neutral and Ionized PAHs: Implications for Interstellar Dust

We present the gas-phase spectroscopy of neutral and ionized polycyclic aromatic hydrocarbons (PAHs) measured in the W-Visible-NIR range in an astrophysically relevant environment. These measurements provide data on PAHs and nanometer sized particles that can now be directly compared to astronomical observations. The harsh physical conditions of the IS medium - characterized by a low temperature, an absence of collisions and strong VUV radiation fields - are simulated in the laborat'ory by associating a molecular beam with an ionizing discharge to generate a cold plasma expansion. PAH ions are formed from the neutral precursors in an isolated environment at low temperature ($\sim lOO$-K). The spectra of neutral and ionized PAHs are measured using the high sensitivity methods of cavity ring down spectroscopy (CRDS) and multiplex integrated cavity output spectroscopy (MICOS). These experiments provide unique information on the spectra of free, cold large carbon molecules and ions in the gas phase. The electronic bands measured for ionized PAH are found to be intrinsically broad ($\geq$20 cm$^{-l}$) while the bands associated with the neutral precursors are narrower (of the order of 2 - 10 cm$^{-l}$). The laboratory data are discussed and compared with recent astronomical spectra of large and narrow DIBs and with the spectra of circumstellar environments of selected carbon stars (see contribution of Cami et al.) and the implications for the interstellar PAH population are derived. Preliminary results also show that carbon nanoparticles are formed during the short residence time of the precursors in the plasma. This finding holds great potential for understanding the formation process of interstellar grains.

Salama, F.↗

Electron-impact total ionization cross sections of DNA sugar-phosphate backbone and an additivity principle

The improved binary-encounter dipole (iBED) model [W.M. Huo, Phys. Rev. A64, 042719-1 (2001)l is used to study the total ionization cross sections of the DNA sugar-phosphate backbone by electron impact. Calculations using neutral fragments found that the total ionization cross sections of C3' - and C5', -deoxyribose-phospate, two conformers of the sugar-phosphate backbone, are close to each other. Furthermore, the sum of the ionization cross sections of the separate deoxyribose and phosphate fragments is in close agreement with the C3' - and C5" -deoxyribose-phospate cross sections, differing by less than 10%. The result implies that certain properties of the-DNA, like the total singly ionization cross section, are localized properties and a building-up or additivity principle may apply. This allows us to obtain accurate properties of larger molecular systems built up from the results of smaller subsystem fragments. Calculations are underway using a negatively charged sugar-phosphate backbone with a metal counter-ion.

Huo, Winifred M.↗

Collisional Ionization Equilibrium for Optically Thin Plasmas

Reliably interpreting spectra from electron-ionized cosmic plasmas requires accurate ionization balance calculations for the plasma in question. However, much of the atomic data needed for these calculations have not been generated using modern theoretical methods and their reliability are often highly suspect. We have utilized state-of-the-art calculations of dielectronic recombination (DR) rate coefficients for the hydrogenic through Na-like ions of all elements from He to Zn. We have also utilized state-of-the-art radiative recombination (RR) rate coefficient calculations for the bare through Na-like ions of all elements from H to Zn. Using our data and the recommended electron impact ionization data of Mazzotta et al. (1998), we have calculated improved collisional ionization equilibrium calculations. We compare our calculated fractional ionic abundances using these data with those presented by Mazzotta et al. (1998) for all elements from H to Ni, and with the fractional abundances derived from the modern DR and RR calculations of Gu (2003a,b, 2004) for Mg, Si, S, Ar, Ca, Fe, and Ni.

Bryans, P.↗

Space Weather Nowcasting of Atmospheric Ionizing Radiation for Aviation Safety

There is a growing concern for the health and safety of commercial aircrew and passengers due to their exposure to ionizing radiation with high linear energy transfer (LET), particularly at high latitudes. The International Commission of Radiobiological Protection (ICRP), the EPA, and the FAA consider the crews of commercial aircraft as radiation workers. During solar energetic particle (SEP) events, radiation exposure can exceed annual limits, and the number of serious health effects is expected to be quite high if precautions are not taken. There is a need for a capability to monitor the real-time, global background radiations levels, from galactic cosmic rays (GCR), at commercial airline altitudes and to provide analytical input for airline operations decisions for altering flight paths and altitudes for the mitigation and reduction of radiation exposure levels during a SEP event. The Nowcast of Atmospheric Ionizing Radiation for Aviation Safety (NAIRAS) model is new initiative to provide a global, real-time radiation dosimetry package for archiving and assessing the biologically harmful radiation exposure levels at commercial airline altitudes. The NAIRAS model brings to bear the best available suite of Sun-Earth observations and models for simulating the atmospheric ionizing radiation environment. Observations are utilized from ground (neutron monitors), from the atmosphere (the METO analysis), and from space (NASA/ACE and NOAA/GOES). Atmospheric observations provide the overhead shielding information and the ground- and space-based observations provide boundary conditions on the GCR and SEP energy flux distributions for transport and dosimetry simulations. Dose rates are calculated using the parametric AIR (Atmospheric Ionizing Radiation) model and the physics-based HZETRN (High Charge and Energy Transport) code. Empirical models of the near-Earth radiation environment (GCR/SEP energy flux distributions and geomagnetic cut-off rigidity) are benchmarked against the physics-based CMIT (Coupled Magnetosphere- Ionosphere-Thermosphere) and SEP-trajectory models.

Mertens, Christopher J.↗