SEARCH · Search NASA

Results for “Gamma Distribution Function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Neutron knockout from 68,70 Ni ground and isomeric states

Neutron-rich isotopes are an important source of new information on nuclear physics. Specifically, the spin-isospin components in the nucleon-nucleon (NN) interaction, e.g., the proton-neutron tensor force, are expected to modify shell structure in exotic nuclei. These potential changes in the intrinsic shell structure are of fundamental interest. The study of the excitation energy of states corresponding to specific configurations in even-even isotopes, together with the single-particle character of the first excited states of odd-A, neutron-rich Ni isotopes, probes the evolution of the neutron orbitals around the Fermi surface as a function of the neutron number a step forward in the understanding of the region and the nature of the NN interaction at large N/Z ratios. In an experiment carried out at the National Superconducting Cyclotron Laboratory, new spectroscopic information was obtained for 68 Ni and the distribution of single-particle strengths in 67,69 Ni was characterized by means of single-neutron knockout from 68,70 Ni secondary beams. The spectroscopic strengths, deduced from the measured partial cross sections to the individual states tagged by their de-exciting gamma rays, is used to identify and quantify configurations that involve neutron excitations across the N = 40 harmonic oscillator shell closure. Finally, the de-excitation rays were measured with the GRETINA tracking array.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Gamma ray Imaging Based Methods to Improve the Accuracy of Uranium Holdup Quantification

Deposits in commercial uranium processing facilities of interest to nuclear safeguards are predominantly low enriched uranium (LEU). The shapes and sizes of holdup deposits can vary significantly, thus making it difficult to model the deposits accurately for calibration purposes. This presents a challenge to currently employed methods such as the Generalized Geometry Holdup (GGH) that rely on quantifying 235U mass by simplifying deposit shapes as a point, a line, or an area. In this work, gamma ray imaging using high energy resolution Germanium Gamma Imagers (GeGI™) is employed to determine the distribution of uranium inside the source containment. Two types of imaging methods are employed; coded-aperture imaging and Compton imaging. Gamma-ray emissions seen in different regions of the image are quantified using an inverse gamma-ray transport solver being developed at ORNL. Using measured data, the inverse method will be used to solve for unknown source parameters such as source matrix thickness, density, and attenuation due to container wall and shielding. The intrinsic efficiency of the detector is determined based on a library of built-in response functions and key dimensions of the detector (e.g., thickness and radius of the detector crystal). Measurements are performed in both imaging modalities using uranium sources of well-known masses and enrichments configured inside mocked-up holdup fixtures. Results from the two imaging modalities are intercompared. The paper presents preliminary results of the imaging measurements and the progress made in the computational methods to quantify uranium mass.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

NASA NTRS2003

Clump Giants in the Hyades

The project is entitled 'Clump Giants in the Hyades.' This observation of one of the late-type Hyades giants (Gamma Tau) has implications for understanding the formation of late-type stellar coronae as a function of the evolutionary state of the star. The Hyades giants are interesting because they are all clump giants in the Helium burning phase, similar to the cool primary of Capella. The Hyades giants show significantly more magnetic activity than expected from their state of evolution (and slowed-down rotation). Thus these systems provide an important clue to dynamo action. The data were obtained by the satellite on 13 March 2001 for a total RGS exposure of 58220 seconds. These data were delivered to the PI on 7 August 2001. The data could not be reprocessed until SAS Version 5.3.3 which became available 7 June 2002. Although the guidelines for assessing background rates suggested that half the data were contaminated, it does not appear that the spectral region of the RGS was adversely affected by unusually high background. The spectra show strong lines of Fe XVII and XVIII, O VII and VIII, Ne IX and X, along with numerous weaker lines. The emission measure distribution is highly reminiscent of Capella; if anything, the emission measure distribution is steeper at 6 million K than for Capella. Gamma Tau is the second brightest of the Hyades clump giants. Pallavicini et al. have shown that the luminosity of the brightest Hyades giant (Theta Tau) is remarkably similar to its luminosity as measured by Einstein. Short-term variability is also modest. We are addressing the variability issue now for Gamma Tau. Initial results were reported at the 2003 Seattle AAS meeting. A paper is in preparation for submission to the Astrophysical Journal.

Mushotzky, Richard F.↗

DOE OSTI2020

A quantitative study of retention and release of deuterium and tritium during irradiation of y-LiAlO2 pellets

This study reports on the quantification of deuterium in ion-irradiated gamma-LiAlO2 pellets as a function of dose and temperature. The gamma-LiAlO2 pellets were sequentially irradiated with He+ and D2+ ions to the same fluences of 5E16, 1E17 and 2E17 He++D+/cm2 at 188 K. Additional irradiation was performed to 1E17, 2E17 and 4E17 He++D+/cm2 at 573 K. A set of the pellets irradiated at 188 K was shipped and stored at low temperatures from 80 to 132 K and characterized using time-of-flight secondary ion mass spectrometry at ~173 K. The deuterium depth profiles show a Gaussian-like distribution in the low-temperature pellets. The total deuterium retention is found to be directly proportional to the ion fluence. About 27 at.% of the implanted deuterium atoms were released from the pellet irradiated to 2E17 He++D+/cm2 at 188 K during storage at room temperature for ~1 month. Retention of the trapped or bound deuterium during ion irradiation at 573 K increases initially with ion fluence and tends to saturate at a high fluence. The amount of the released deuterium is observed to be quantitatively consistent with that of the released tritium from similar standard pellets during neutron irradiation at 573 K.

Deuterium retention and release, ion irradiation, ↗

NASA NTRS1983

Luminosity function and cosmological evolution of X-ray selected quasars

The preliminary analysis of a complete sample of 55 X-ray sources is presented as part of the Medium Sensitivity Survey of the Einstein Observatory. A pure luminosity evolutionary law is derived in order to determine the uniform distribution of the sources and the rates of evolution for Active Galactic Nuclei (AGNs) observed by X-ray and optical techniques are compared. A nonparametric representation of the luminosity function is fitted to the observational data. On the basis of the reduced data, it is determined that: (1) AGNs evolve cosmologically; (2) less evolution is required to explain the X-ray data than the optical data; (3) the high-luminosity portion of the X-ray luminosity can be described by a power-law with a slope of gamma = 3.6; and (4) the X-ray luminosity function flattens at low luminosities. Some of the implications of the results for conventional theoretical models of the evolution of quasars and Seyfert galaxies are discussed.

Maccacaro, T.↗

DOE OSTI2019

Convenient analytical solution for vibrational distribution function of molecules colliding with a wall

We study formation of the Vibrational Distribution Function (VDF) in a molecular gas at low pressure, when vibrational levels are excited by electron impact and deactivated in collisions with walls and show that this problem has a convenient analytical solution that can be used to obtain VDF and its dependence on external parameters. The VDF is determined by excitation of vibrational levels by an external source and deactivation in collisions with the wall. Deactivation in wall collisions is little known process. However, we found that the VDF is weakly dependent on the functional form of the actual form of probability ${\gamma }_{\upsilon ^{\prime} \to \upsilon }$ for a vibrational number $\upsilon ^{\prime} $ to transfer into a lower level $\upsilon $ at the wall. Because for a given excitation source of vibrational states, the problem is linear the solution for VDF involves solving linear matrix equation. The matrix equation can be easily solved if we approximate probability, in the form: ${\gamma }_{\upsilon ^{\prime} \to \upsilon }=(1/\upsilon ^{\prime} )\theta (\upsilon ^{\prime} -\upsilon ).$ In this case, the steady-state solution for VDF($\upsilon $) simply involves a sum of source rates for levels above $\upsilon ,$ with a factor of $1/(\upsilon +1).$ As an example of application, we study the vibrational kinetics in a hydrogen gas and verify the analytical solution by comparing with a full model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2005

Testing the Gamma-Ray Burst Energy Relationships

Building on Nakar & Piran's analysis of the Amati relation relating gamma-ray burst peak energies E(sub p) and isotropic energies E(sub iso ) we test the consistency of a large sample of BATSE bursts with the Amati and Ghirlanda (which relates peak energies and actual gamma-ray energies E(sub gamma)) relations. Each of these relations can be exp ressed as a ratio of the different energies that is a function of red shift (for both the Amati and Ghirlanda relations) and beaming fraction f(sub B) (for the Ghirlanda relation). The most rigorous test, whic h allows bursts to be at any redshift, corroborates Nakar & Piran's r esult - 88% of the BATSE bursts are inconsistent with the Amati relat ion - while only l.6% of the bursts are inconsistent with the Ghirlan da relation if f(sub B) = 1. Modelling the redshift distribution resu lts in an energy ratio distribution for the Amati relation that is sh ifted by an order of magnitude relative to the observed distributions; any sub-population satisfying the Amati relation can comprise at mos t approx. 18% of our burst sample. A similar analysis of the Ghirland a relation depends sensitively on the beaming fraction distribution f or small values of f(sub B); for reasonable estimates of this distrib ution about a third of the burst sample is inconsistent with the Ghir landa relation. Our results indicate that these relations are an artifact of the selection effects of the burst sample in which they were f ound; these selection effects may favor sub-populations for which the se relations are valid.

Band, David L.↗

DOE OSTI2025

The Psyche Gamma-Ray and Neutron Spectrometer

A Gamma-Ray and Neutron Spectrometer (GRNS) instrument has been developed as part of the science payload for NASA’s Discovery Program Psyche mission to the M-class asteroid (16) Psyche. The GRNS instrument is designed to measure the elemental composition of Psyche with the goal to understand the origin of this mysterious, potentially metal-rich planetary body. The GRNS will measure the near-surface abundances for the elements Ni, Fe, Si, K, S, Al, and Ca, as well as the spatial distribution of Psyche’s metal-to-silicate fraction (or metal fraction). These measurements address three of the five Psyche mission science objectives: determine if Psyche is a core; determine whether small metal bodies incorporate light elements into the metal phase; and determine whether Psyche was formed under reducing conditions. The Gamma-Ray Spectrometer (GRS) uses a cryocooled, high-purity Ge (HPGe) sensor to detect cosmic-ray generated gamma rays in the 60 to 9000-keV energy range. The HPGe sensor is surrounded by a borated plastic anticoincidence shield that provides three functions: active background rejection from charged particle interactions in the HPGe sensor; fast neutron measurements; and direct measurements of the incident galactic cosmic ray flux. The Neutron Spectrometer (NS) uses three 3 He gas proportional sensors, each with different material wraps to measure thermal (<0.4 eV), low-energy epithermal (0.4 eV to 1 keV), and high-energy epithermal (up to 100 keV) neutrons. This paper provides an overview of the Psyche GRNS, including: its science and measurement objectives; the design of the instrument hardware, software, and operation; pre-launch performance measurements and its initial performance in space; and an overview of its data products and expected operation for different Psyche mission phases.

Engineering - Instrumentation related to nuclear s↗

DOE OSTI2019

Multiple Components in the Broadband γ -Ray Emission of the Short GRB 160709A

GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The γ-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole γ-ray spectrum. While the thermal component dominates during ~0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of ~340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, ${R}_{0}={3.8}_{-1.8}^{+5.9}$ × 10 8 cm, the size of the photosphere, ${R}_{\mathrm{ph}}={7.4}_{-1.2}^{+0.8}$ × 10 10 cm, and a bulk Lorentz factor, ${\rm{\Gamma }}={728}_{-93}^{+75}$, assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV–40 MeV and 100 MeV–10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, p = 2.5 ± 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2006

The Supercritical Pile Model for GRBs

We present the spectral and temporal radiative signatures expected within the Supercritical Pile model of Gamma Ray Bursts (GRB). This model is motivated by the need for a process that provides the dissipation necessary in GRB and presents a well defined scheme for converting the energy stored in the relativistic protons of the Relativistic Blast Waves (RBW) associated with GRB into radiation; at the same time it leads to spectra which exhibit a peak in the burst nuF(sub nu) distribution at an energy E(sub p) approximately equal to 1 MeV in the observer s frame, in agreement with observation and largely independent of the Lorentz factor GAMMA of the associated relativistic outflow. Furthermore, this scheme does not require (but does not preclude) acceleration of particles at the shock other than that provided by the isotropization of the flow bulk kinetic energy on the RBW frame. In the present paper we model in detail the evolution of protons, electrons and photons from a RBW to produce detailed spectra of the prompt GRB phase as a function of time from across a very broad range spanning roughly 4 log10 GAMMA decades in frequency. The model spectra are in general agreement with observations and provide a means for the delineating of the model parameters through direct comparison with trends observed in GRB properties.

Kazanas, Demos↗

NASA NTRS2011

Fermi/GBM Observations of SGRJ0501 + 4516 Bursts

We present our temporal and spectral analyses of 29 bursts from SGRJ0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T(sub 90) durations of the bursts can be fit with a log-normal distribution with a mean value of approx. 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T(sub 90)s estimated in count space (following a log-normal distribution with a mean value of approx. 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that E(sub peak) decreases with energy flux (and fluence) to a minimum of approx. 30 keV at F = 8.7 x 10(exp -6)erg/sq cm/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550 - 5418 and 1806 - 20. The isotropic luminosity, L(sub iso), corresponding to these flux values is roughly similar for all sources (0.4 - l.5 x 10(exp 40) erg/s.

Lin, Lin↗

DOE OSTI2020

Discovery of a Gamma-Ray Black Widow Pulsar by GPU-accelerated Einstein@Home

We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations—whose periodic brightness modulations suggested an orbit—no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multidimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with B surf ≈ 4 × 10 7 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass gsim 2 M ⊙ . The companion is lightweight with mass ~0.01 M ⊙ , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A hierarchical analysis of the impact of methodological decisions on statistical downscaling of daily precipitation and air temperatures

Despite the widespread application of statistical downscaling tools, uncertainty remains regarding the role of model formulation in determining model skill for daily maximum and minimum temperature (Tmax and Tmin), and precipitation occurrence and intensity. Impacts of several key aspects of statistical transfer function form on model skill are evaluated using a framework resistant to model overspecification. We focus on: (a) model structure: simple (generalized linear models, GLMs) versus complex (artificial neural networks, ANNs) models. (b) Predictor selection: Fixed number of predictors chosen a priori versus stepwise selection of predictors and inclusion of grid point values versus predictors derived from application of principal components analysis (PCA) to spatial fields. We also examine the influence of domain size on model performance. For precipitation downscaling, we consider the role of the threshold used to characterize a wet day and apply three approaches (Poisson and Gamma distributions in GLM and ANN) to downscale wet-day precipitation amounts. While no downscaling formulation is optimal for all predictands and at 10 locations representing diverse U.S. climates, and due to the exclusion of variance inflation all of the downscaling formulations fail to reproduce the range of observed variability, models with larger suites of prospective predictors generally have higher skill. For temperature downscaling, ANNs generally outperform GLM, with greater improvements for T min than T max . Furthermore, use of PCA-derived predictors does not systematically improve model skill, but does improve skill for temperature extremes. Model skill for precipitation occurrence generally increases as the wet-day threshold increases and models using PCA-derived predictors tend to outperform those based on grid cell predictors. Each model for wet-day precipitation intensity overestimates annual total precipitation and underestimates the proportion derived from extreme precipitation events, but ANN-based models and those with larger predictor suites tend to have the smallest bias.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A semi-supervised learning method to produce explainable radioisotope proportion estimates for NaI-based synthetic and measured gamma spectra

Quantifying the radioactive sources present in gamma spectra is an ever-present and growing national security mission and a time-consuming process for human analysts. While machine learning models exist that are trained to estimate radioisotope proportions in gamma spectra, few address the eventual need to provide explanatory outputs beyond the estimation task. In this work, we develop two machine learning models for a NaI detector measurements: one to perform the estimation task, and the other to characterize the first model’s ability to provide reasonable estimates. To ensure the first model exhibits a behavior that can be characterized by the second model, the first model is trained using a custom, semi-supervised loss function which constrains proportion estimates to be explainable in terms of a spectral reconstruction. The second auxiliary model is an out-of-distribution detection function (a type of meta-model) leveraging the proportion estimates of the first model to identify when a spectrum is sufficiently unique from the training domain and thus is out-of-scope for the model. In demonstrating the efficacy of this approach, we encourage the use of meta-models to better explain ML outputs used in radiation detection and increase trust.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2003

Chandra X-Ray Observatory Observations of the Globular Cluster M28 and its Millisecond Pulsar PSR B1821-24

We report here the results of the first Chandra X-Ray Observatory observations of the globular cluster M28 (NGC 6626). We detect 46 X-ray sources of which 12 lie within one core radius of the center. We show that the apparently extended X-ray core emission seen with the ROSAT HRI is due to the superposition of multiple discrete sources for which we determine the X-ray luminosity function down to a limit of about 6 x 10(exp 30) erg/s. We measure the radial distribution of the X-ray sources and fit it to a King profile finding a core radius of r(sub c,x) approx. 11 sec. We obtain the best-fit mass of the X-ray sources to be M(sub x) approx. 1.9 solar masses. We measure for the first time the unconfused phase-averaged X-ray spectrum of the 3.05-ms pulsar B1821-24 and find it best described by a power law with photon-index Gamma approx. equal to 1.2. We find marginal evidence of an emission line centered at 3.3 kev in the pulsar spectrum, which could be interpreted as cyclotron emission from a corona above the pulsar's polar cap if the the magnetic field is strongly different from a centered dipole. The unabsorbed pulsar flux in the 0.5-8.0 keV band is approx. 3.5 x 10(exp -13) ergs/s/sq cm. We present spectral analyses of the 5 brightest unidentified sources. Based on the spectral parameters of the brightest of these sources, we suggest that it is a transiently accreting neutron star in a low-mass X-ray binary, in quiescence. Fitting its spectrum with a hydrogen neutron star atmosphere model yields the effective temperature T(sup infinity)(sub eff) = 90(sup +30)(sub -10) eV and the radius R(sup infinity)(sub NS) = 14.5(sup +6.9)(sub -3.8) km. In addition to the resolved sources, we detect fainter, unresolved X-ray emission from the central core. Using the Chandra-derived positions, we also present a preliminary report on the result of searching archival Hubble Space Telescope data for possible optical counterparts.

Becker, Werner↗

NASA NTRS1993

Structure in gamma ray burst time profiles: Correlations with other observables

One of the current debates raging in the world of gamma-ray burst physics is whether the sources of these enigmatic bursts arise from a single or from multiple distributions. Several authors contend that the histograms of GRB observables imply the latter. The two most-likely candidate components are galactic and cosmological. For example, some researchers claim that a dip in the V/V max distribution is a result of such a two-component source distribution. Others have used a parameter called the 'burst variability' calculated by dividing the maximum count rate on the 64-msec timescale by that from the 1024-msec timescale to show that a correlation of this parameter with bursts brightness implies a two-component model. This method has met vigorous criticism. We have developed two parameters that measure the variability or structure in the time profiles of BATSE gamma-ray bursts. Both parameters ('structure' and 'spike height') are based on the statistics of 'runs up' and 'runs down.' In short, the structure parameter is the observed number of runs (at several lengths) minus the number expected in a chance distribution. The 'spike height' is the sum of all run heights minus the expected sum. These two are straight-forward to calculate, robust, and measure the variability over the complete profile--not just at the peak. We have applied this algorithm to the profiles of 156 GRB's. In this paper we present graphs of the two parameters as functions of the following: (1) burst duration, (2) burst hardness ratio, (3) V/V max, (4) source galactic longitude, and (5) source galactic latitude. We seek correlations as well as groupings in the data that might indicate a multi-component source distribution.

Lestrade, John Patrick↗

NASA NTRS1992

BATSE analysis techniques for probing the GRB spatial and luminosity distributions

The Burst And Transient Source Experiment (BATSE) has measured homogeneity and isotropy parameters from an increasingly large sample of observed gamma-ray bursts (GRBs), while also maintaining a summary of the way in which the sky has been sampled. Measurement of both of these are necessary for any study of the BATSE data statistically, as they take into account the most serious observational selection effects known in the study of GRBs: beam-smearing and inhomogeneous, anisotropic sky sampling. Knowledge of these effects is important to analysis of GRB angular and intensity distributions. In addition to determining that the bursts are local, it is hoped that analysis of such distributions will allow boundaries to be placed on the true GRB spatial distribution and luminosity function. The technique for studying GRB spatial and luminosity distributions is direct. Results of BATSE analyses are compared to Monte Carlo models parameterized by a variety of spatial and luminosity characteristics.

Hakkila, Jon↗

NASA NTRS1994

Removing cosmic-ray hits from multiorbit HST Wide Field Camera images

We present an optimized algorithm that removes cosmic rays ('CRs') from multiorbit Hubble Space Telescope (HST) Wide Field/Planetary Camera ('WF/PC') images. It computes the image noise in every iteration from the WF/PC CCD equation. This includes all known sources of random and systematic calibration errors. We test this algorithm on WF/PC stacks of 2-12 orbits as a function of the number of available orbits and the formal Poissonian sigma-clipping level. We find that the algorithm needs greater than or equal 4 WF/PC exposures to locate the minimal sky signal (which is noticeably affected by CRs), with an optimal clipping level at 2-2.5 x sigma(sub Poisson). We analyze the CR flux detected on multiorbit 'CR stacks,' which are constructed by subtracting the best CR filtered images from the unfiltered 8-12 orbit average. We use an automated object finder to determine the surface density of CRS as a function of the apparent magnitude (or ADU flux) they would have generated in the images had they not been removed. The power law slope of the CR 'counts' (gamma approximately = 0.6 for N(m) m(exp gamma)) is steeper than that of the faint galaxy counts down to V approximately = 28 mag. The CR counts show a drop off between 28 less than or approximately V less than or approximately 30 mag (the latter is our formal 2 sigma point source sensitivity without spherical aberration). This prevents the CR sky integral from diverging, and is likely due to a real cutoff in the CR energy distribution below approximately 11 ADU per orbit. The integral CR surface density is less than or approximately 10(exp 8)/sq. deg, and their sky signal is V approximately = 25.5-27.0 mag/sq. arcsec, or 3%-13% of our NEP sky background (V = 23.3 mag/sq. arcsec), and well above the EBL integral of the deepest galaxy counts (B(sub J) approximately = 28.0 mag/sq. arcsec). We conclude that faint CRs will always contribute to the sky signal in the deepest WF/PC images. Since WFPC2 has approximately 2.7x lower read noise and a thicker CCD, this will result in more CR detections than in WF/PC, potentially affecting approximately 10%-20% of the pixels in multiorbit WFPC2 data cubes.

Windhorst, Rogier A.↗