Solar flares and particle acceleration

Energy release in solar flares occurs during the impulsive phase, which is a period of a few to about ten minutes, during which energy is injected into the flare region in bursts with durations of various time scales, from a few tens of seconds down to 0.1 s or even shorter. Non-thermal heating is observed during a short period, not longer than a few minutes, in the very first part of the impulsive phase; in average flares, with ambient particle densities not larger than a few times 10¹⁰ cm^–3 it is due to thick-target electron beam injection, causing chromospheric ablation followed by convection. In flares with larger densities the heating is due to thermal fronts (Section 1). The average energy released in chromospheric regions is a few times 10³⁰ erg, and an average number of 10³⁸ electrons with E 15 keV is accelerated. In subsecond pulses these values are about 10³⁵ electrons and about 10²⁷ erg per subsecond pulse. The total energy released in flares is larger than these values (Section 2). Energization occurs gradually, in a series of fast non-explosive flux-thread interactions, on the average at levels about 10⁴ km above the solar photosphere, a region permeated by a large number ( 10) of fluxthreads, each carrying electric currents of 10¹⁰–10¹¹ A. The energy is fed into the flare by differential motions of magnetic fields driven by photospheric-chromospheric movements (Section 3). In contrast to these are the high-energy flares, characterized by the emission of gamma-radiation and/or very high-frequency (millimeter) radiobursts. Observations of such flares, of the flare neutron emission, as well as the observation of ³He-rich interplanetary plasma clouds from flares all point to a common source, identified with shortlived ( 0.1 s) superhot ( 10⁸ K) flare knots, situated in chromospheric levels (Section 4). Pre-flare phenomena and the existence of homologous flares prove that flare energization can occur repeatedly in the same part of an active region: the consequent conclusions are that only seldom the full energy of an active region is exhausted in one flare, or that the flare energy is generated anew between homologous flares; this latter case looks more probable (Section 5). Flare energization requires the formation of direct electric fields, in value comparable with, or somewhat smaller than the Dreicer field (Section 6). Such fields originate by current-thread reconnection in a regime in which the current sheet is thin enough to let resistive instability originate (Section 7). Particle acceleration occurs by fast reconnection in magnetic fields 100 G and electric fields exceeding about 0.3 times the Dreicer field at fairly low particle densities ( 10¹⁰ cm^–3); for larger densities plasma heating is expected to occur (Section 8). Transport of accelerated particles towards interplanetary space demands a field-line configuration open to space. Such a configuration originates mainly after the gradual gamma-ray/proton flares, and particularly after two-ribbon flares; these flares belong to the dynamic flares in Sturrock and vestka's flare classification. Acceleration to GeV energies occurs subsequently in shock waves, probably by first-order Fermi acceleration (Section 9).     

COS-B views on the diffuse galactic gamma-ray emission and some point sources

The correlation between diffuse galactic gamma rays and gas tracers is studied using the final COS-B database and H i and CO surveys covering the entire galactic plane. A good quantitative fit to the gamma rays is obtained, with a small galacto-centric gradient in the gamma-ray emissivity per hydrogen atom. The average ratio of H₂ column density to integrated CO temperature is determined, the best estimate being (2.3 ± 0.3) × 10² molecules cm^–2 (K km s^–1)^–1. Strictly taken, this value is an upper limit. The corresponding mass of molecular hydrogen in the inner galaxy, derived using both 1st and 4th quadrants, is 1.0 × 10⁹ M .The softer gamma-ray spectrum towards the inner galaxy found in previous work can be attributed to a steeper emissivity gradient at low energies and/or to a softer gamma-ray spectrum of the emission distributed like molecular gas. A steeper emissivity gradient at low energies could be related to cosmic-ray spectral variations in the Galaxy, to different distributions of cosmic-ray electrons and nuclei, or to a contribution from discrete sources. A softer spectrum for the emission associated with molecular clouds may be physically related to the clouds themselves (i.e., cosmic-ray spectral variations) or to an associated discrete source distribution.New results on the temporal and spectral characteristics of the high-energy (50 MeV to 5 GeV) gammaray emission from the Vela pulsar are presented. The whole pulsed flux is found to exhibit long-term variability. Five discrete emission regions within the pulsar lightcurve have been identified, with the spectral characteristics and long-term behaviour being different. These characteristics differ significantly from those reported earlier for the Crab pulsar. However, geometrical pulsar models have been proposed (e.g., Morini, 1983; Smith, 1986) which could explain many of these features.     

The local bubble

Recently, observations with the rosat PSPC instrument and the spectrometers onboard the euve satellite have given new detailed information on the structure and physical conditions of the Local Bubble. From the early rocket experiments, and in particular from the WISCONSIN Survey, the existence of a diffuse hot gas in the vicinity of the solar system, extending out to about 100 pc, has been inferred in order to explain the emission below 0.3 keV. The higher angular resolution and sensitivity of rosat made it possible to use diffuse neutral clouds as targets for shadowing the soft X-ray background. Thus, in some directions, more than half of the flux in the 0.25 keV band appears to come from outside the Local Bubble. Further, measurements of the diffuse EUV in the LISM, show surprisingly few emission lines. These findings are in conflict with the standard LHB model, which assumes a local hot (T 10⁶ K) plasma in CIE. Model calculations, based on the non-equilibrium cooling of an expanding plasma, show a promising way of reconciling all available observations. Thus the present temperature within the LB may be as low as 4 × 10⁴ K and its number density as large as 2 × 10^–2 cm ^–3, giving a total pressure that is roughly in agreement with the Local Cloud.Abbreviations CIE collisional ionization equilibrium - ISM Interstellar Medium - LHB Local Hot Bubble - LB Local Bubble - LISM Local ISM - SB superbubble - SXR soft X-ray - SXRB SXR Background - VLISM Very Local ISM Heisenberg Fellow     

Starburst galaxies

Bursts of massive star formation can temporarily dominate the luminosity of galaxies spanning a wide range of morphological types. This review is concerned primarily with such events in the central 1 kpc region of spiral galaxies which result from bar driven inflows of gas triggered by interactions or mergers. Most of the stellar radiant luminosity of such bursts is absorbed by dust and re-emitted in the far-infrared and is accompanied by radio and X-ray emission from supernova remnants which can also act collectively to drive galaxy scale outflows. Both evolutionary stellar models and estimates of the gas depletion times are consistent with typical burst durations of 10^7–8 yr. Spatially-resolved studies of nearby starburst galaxies reveal that the activity is distributed over many individual star forming complexes within rings and other structures organized by interactions between bars and the disc over a range of scales. More distant and extreme examples associated with mergers of massive spirals have luminosities > 10¹³ L and molecular gas masses > 10¹⁰ M implying star formation rates > 1000 M yr^–1 which can only be sustained for 10⁷ yr. In the most luminous merging systems, however, the relative importance of starburst and AGN activity and their possible evolutionary connection is still a hotly debated issue. Also controversial are suggestions that starbursts in addition to a black hole are required to account for the properties of AGNs or that starbursts alone may be sufficient under certain conditions. In a wider context, starbursts must clearly have played an important role in galaxy formation and evolution at earlier times. Recent detections of high redshift galaxies show that star formation was underway at z 4 but do not support a continuing increase of the strong evolution in the co-moving star formation density seen out to z l. Primeval starburst pre-cursors of spheroidal systems also remain elusive. The most distant candidates are radio galaxies and quasars at z = 4–5 and a possible population of objects responsible for an isotropic sub-mm wave background tentatively claimed to have been detected by the COBE satellite.     

The flux of interstellar dust observed by Ulysses and Galileo

Interstellar dust detected by the dust sensor onboard Ulysses was first identified after the Jupiter flyby when the spacecraft's trajectory changed dramatically (Grün et al., 1994). Here we report on two years of Ulysses post-Jupiter data covering the range of ecliptic latitudes from 0° to –54° and distances from 5.4 to 3.2 AU. We find that, over this time period, the flux of interstellar dust particles with a mean mass of 3·10^–13 g stays nearly constant at about 1·10^–4, m^–2 s^–1 ( sr)^–1, with both ecliptic latitude and heliocentric distance.Also presented are 20 months of measurements from the identical dust sensor onboard the Galileo spacecraft which moved along an in-ecliptic orbit from 1.0 to 4.2 AU. From the impact direction and speeds of the measured dust particles we conclude that Galileo almost certainly sensed interstellar dust outside 2.8 AU; interstellar particles may also account for part of the flux seen between 1 and 2.8 AU.     

Characteristics of nearby interstellar matter

There is a warm tenuous partially ionized cloud (T10⁴ K,n(HI)0.1 cm^–3,n(Hii 0.22–0.44 cm^–3) surrounding the solar system which regulates the environment of the solar system, determines the structure of the heliopause region, and feeds neutral interstellar gas into the inner solar system. The velocity (V–20 km s^–1 froml335°,b0° in the local standard of rest) and enhanced Caii and Feii abundances of this cloud suggest an origin as evaporated gas from cloud surfaces in the Scorpius-Centaurus Association. Although the soft X-ray emission attributed to the Local Bubble is enigmatic, optical and ultraviolet data are consistent with bubble formation caused by star formation epochs in the Scorpius-Centaurus Association as regulated by the nearby spiral arm configuration. The cloud surrounding the solar system (the local fluff) appears to be the leading region of an expanding interstellar structure (the squall line) which contains a magnetic field causing polarization of the light of nearby stars, and also absorption features in nearby upwind stars. The velocity vectors of the solar system and local fluff are perpendicular in the local standard of rest. Combining this information with the low column densities seen towards Sirius in the anti-apex direction, and the assumption that the cloud velocity vector is parallel to the surface normal, suggests that the Sun entered the local fluff within the historical past (less than 10 000 years ago) and is skimming the surface of the cloud. Comparison of magnesium absorption lines towards Sirius and anomalous cosmic-ray data suggest the local fluff is in ionization equilibrium.Reason has moons, but moons not hers, Lie mirror'd on her sea, Confounding her astronomers, But, O! delighting me.Ralph Hodgson     

New Experimental Data and What It Tells Us About the Sources and Acceleration of Cosmic Rays

This paper summarizes new data in several fields of astronomy that relate to the origin and acceleration of cosmic rays in our galaxy and similar nearby galaxies. Data from radio astronomy shows that supernova remnants, both in our galaxy and neighboring galaxies, appear to be the sources of most of the accelerated electrons observed in these galaxies. -ray measurements also reveal several strong sources associated with supernova remnants in our galaxy. These sources have -ray spectra that are suggestive of the acceleration of cosmic-ray nuclei. Cosmic-ray observations from the Voyager and Ulysses spacecraft suggest a source composition that is very similar to the solar composition but with distinctive differences in the ⁴He, ¹²C,¹⁴ N and ²²Ne abundances that are the imprint of giant W-R star nucleosynthesis. Injection effects which depend on the first ionization potential (FIP) of the elements involved are also observed, in a manner similar to the fractionization observed between the solar photosphere and corona and also analogous to the preferential acceleration observed for high FIP elements at the heliospheric solar wind termination shock. Most of the ⁵⁹Ni produced in the nucleosynthesis of Fe peak nuclei just prior to a SN explosion appears to have decayed to ⁵⁹Co before the cosmic rays have been accelerated, suggesting that the⁵⁹ Ni is accelerated at least 10⁵ yr after it is produced. The decay of certain K capture isotopes produced during cosmic-ray propagation has also been observed for the first time. These measurements suggest that re-acceleration after an initial principal acceleration cannot be large. The high energy spectral indices of cosmic-ray nuclei show a significant charge dependent trend with the index of hydrogen being -2.76 and that of Fe -2.61. The escape length dependence of cosmic rays from our galaxy can now be measured up to ~300 GeV nucl^-1 using the Fe sec/Fe ratio. This escape length is P ^-0.05 above 10 GeV nucl^-1 leading to a typical source spectral index of (2.70±0.10) -0.50 = -2.20 for nuclei. This is similar to the source index of -2.3 inferred for electrons within the errors of ±0.1 in the index for both components. Spacecraft measurements in the outer heliosphere suggest that the local cosmic-ray energy density is ~2eV cm^-3 – larger than previously assumed. Gamma-ray measurements of electron bremsstrahlung below 50 MeV from the Comptel experiment on CGRO show that fully 20–30% of this energy is in electrons, several times that previously assumed. New estimates of the amount of matter traversed by cosmic rays using measurements of the B/C ratio are also higher than earlier estimates – this value is now ~10 g cm^-2 at 1 GeV nucl^-1. Thus altogether cosmic rays are energetically a more important component of our galaxy than previously assumed. This has implications both for the types of sources that are capable of accelerating cosmic rays and also for the role that cosmic rays may play in ionizing the diffuse interstellar medium.     

Spectral observation of the composite supernova remnant G 29.7-0.3

The X-ray properties of the supernova remnant G 29.7-0.3 are discussed based on spectral data from the EXOSAT satellite. In the 2 to 10 keV range a featureless power-law spectrum is obtained, the best-fit parameters being: energy spectral index =-0.77, hydrogen column density on the line of sight N_H=2.3.10²² cm^–2. The incident X-ray flux from the source is (3.6±0.1) 10¹¹ erg cm^–2 s^–1 in the 2 to 10 keV range corresponding to an intrinsic luminosity of about 2. 10³⁶ erg s^–1 for a distance of 19 kpc. The source was not seen with the imaging instrument thus constraining the hydrogen column density to be N_H=(3.3 ±0.3) 10²² cm^–2 and the energy spectral index =1.0±0.15. This new observation is consistent with emission by a synchroton nebula presumably fed by an active pulsar. An upper limit of 1.5% for the pulsed fraction in the range of periods 32ms to 10⁴ s has been obtained.     

Recent results on the parameters of the interstellar helium from the ULYSSES/GAS experiment

Velocity and direction of the flow of the interstellar helium and its temperature and density have been determined from the measurements of the ULYSSES/GAS experiment for two different epochs: during the in-ecliptic path of ULYSSES, representing solar maximum conditions, and during the south to the north pole transition (11/94-6/95), close to the solar minimum conditions. Within the improved error bars the values are consistent with results published earlier.The determination of the density n of the interstellar helium at the heliospheric boundary from observations in the inner solar system requires knowledge about the loss processes experienced by the particles on their way to the observer. The simultaneous observation of the helium particles arriving on direct and indirect orbits at the observer provides a tool to directly determine the effects of the loss processes assumed to be predominantly photoionization and — for particles travelling close to the Sun — electron impact ionization by high-energy solar wind electrons.Such observations were obtained with the ULYSSES/GAS instrument in February 1995, before the spaceprobe passed its perihelion. From these measurements values for the loss rates and the interstellar density could be derived. Assuming photoionization to be the only loss process reasonable fits to the observations were obtained for an ionization rate = 1.1 · 10^–7 s^–1 and a density n 1.7 · 10^–2 cm^–3. Including, in addition, electron impact ionization, a photoionization = 0.6 · 10^–7 s^–1 was sufficient to fit both observations, resulting in a density n 1.4 · 10^–2 cm^–3.On leave from Space Research Centre, Warsaw, Poland.     

Large scale structure of the universe

The current situation with the cosmological model and fundamental constants is briefly reviewed. Here, we concentrate on evolutionary effects of large-scale structure formation, in particular, the relationship with the quasar distribution and dynamics is discussed. We argue that groups of bright quasars with few or more than dozen of members within regions l _LS(100–150)h ^–1 Mpc found atz<2 may belong to concentrations of young rich clusters of galaxies, and thus be distant Great Attractors like the local GA or the Shapley concentration. These early large-scale galactic structures (i) provide a natural way to bias the distribution of Abell clusters, and (ii) suggest that the spectrum of primordial density perturbations is nearly flat on scales encompassing both the cluster and GAs,l=k ^–1(10,100)h ^–1 Mpc: _k ² k ³ P(k) k , =1 _–0.4 ^+0.6 , whereP(k) is the power spectrum of density perturbations.     

In situ measurements of interstellar dust with the Ulysses and Galileo spaceprobes

Interstellar dust was first identified by the dust sensor onboard Ulysses after the Jupiter flyby in February 1992. These findings were confirmed by the Galileo experiment on its outbound orbit from Earth to Jupiter. Although modeling results show that interstellar dust is also present at the Earth orbit, a direct identification of interstellar grains from geometrical arguments is only possible outside of 2.5 AU. The flux of interstellar dust with masses greater than 6 · 10^–14 g is about 1 · 10^–4 m ^–2 s ^–1 at ecliptic latitudes and at heliocentric distances greater than 1AU. The mean mass of the interstellar particles is 3 · 10^–13 g. The flux arrives from a direction which is compatible with the influx direction of the interstellar neutral Helium of 252° longitude and 5.2° latitude but it may deviate from this direction by 15 – 20°.     

LISM structure — Fragmented superbubble shell?

Small scale structure in local interstellar matter (LISM) is considered. Overall morphology of the local cloud complex is inferred from Ca II absorption lines and observations of H I in white dwarf stars. Clouds with column densities ranging from 2–100 × 10¹⁷ cm^–2 are found within 20 pc of the Sun. Cold (50 K) dense (10⁵ cm^–3) small (5–10 au) clouds could be embedded and currently undetected in the upwind gas. The Sun appears to be embedded in a filament of gas with thickness 0.7 pc, and cross-wise column density 2 × 10¹⁷ cm^–2. The local magnetic field direction is parallel to the filament, suggesting that the physical process causing the filamentation is MHD related. Enhanced abundances of refractory elements and LISM kinematics indicate outflowing gas from the Scorpius-Centaurus Association. The local flow vector and Sco data are consistent with a 4,000,000 year old superbubble shell at –22 km s^–1, which is a shock front passing through preshock gas at –12 km s^–1, and yielding cooled postshock gas at –26 km s^–1in the upwind direction. A preshock magnetic field strength of 1.6 G, and postshock field strength of 5.2 G embedded in the superbubble shell, are consistent with the data.Abbreviations LISM Local ISM - SIC Surrounding Interstellar Cloud - LIC Local Interstellar Cloud     

Early nucleosynthesis: The present situation

The most recent developments in the very light element abundance observations, especially the determination of a high (D/H) ratio 3 10^–4 on a remote QSO line of sight, are analyzed and some of the consequences on the Big Bang nucleosynthesis (BBN) models, on the predicted baryonic density and on the galactic evolution schemes are reviewed. Should this (D/H) determination be confirmed by further observations, the simplest standard BBN model would be the most successful in these abundance predictions: the baryon density should be close to the luminous mass density ^B (5–8)10^–3 if H=100 km s^–1 Mpc^–1; the challenge is now to devise galactic evolution models able to account for a large D destruction and avoiding to overproduce³He and metals.     

Tailored analyses of 24 Galactic WN stars

The fundamental properties of 24 Galactic WN stars are determined from analyses of their optical, UV and IR spectra using sophisticated model atmosphere codes (Hillier, 1987, 1990). Terminal velocities, stellar luminosities, temperatures, mass loss rates and abundances of hydrogen, helium, carbon, nitrogen and oxygen are determined. Stellar parameters are derived using diagnostic lines and interstellar reddenings found from fitting theoretical continua to observed energy distributions.Our results confirm that the parameters of WN stars span a large range in temperature (T_*=30–90,000 K), luminosity (log L_*/L=4.8–5.9), mass loss (M=0.9–12×10^–5 M yr^–1) and terminal velocity (v =630–3300 km s^–1). Hydrogen abundances are determined, and found to be low in WNEw and WNEs stars (<15% by mass) and considerable in most WNL stars (1–50%). Metal abundances are also determined with the nitrogen content found to lie in the range N/He=1–5×10^–3 (by number) for all subtypes, and C/N 0.02 in broad agreement with the predictions of Maeder (1991). Enhanced O/N and O/C is found for HD 104994 (WN3p) suggesting a peculiar evolutionary history. Our results suggest that single WNL+abs stars may represent an evolutionary stage immediately after the Of phase. Since some WNE stars exist with non-negligible hydrogen contents (e.g. WR136) evolution may proceed directly from WNL+abs to WNE in some cases, circumventing the luminous blue variable (LBV) or red supergiant (RSG) stage.     

Fast Dust in the Heliosphere

The dynamics of dust particles in the solar system is dominated by solar gravity, by solar radiation pressure, or by electromagnetic interaction of charged dust grains with the interplanetary magnetic field. For micron-sized or bigger dust particles solar gravity leads to speeds of about 30 to 40 km s^–1 at the Earths distance. Smaller particles that are generated close to the Sun and for which radiation pressure is dominant (the ratio of radiation pressure force over gravity F _rad/F _grav is generally termed ) are driven out of the solar system on hyperbolic orbits. Such a flow of -meteoroids has been observed by the Pioneer 8, 9 and Ulysses spaceprobes. Dust particles in interplanetary space are electrically charged to typically +5 V by the photo effect from solar UV radiation. The dust detector on Cassini for the first time measured the dust charge directly. The dynamics of dust particles smaller than about 0.1 m is dominated by the electromagnetic interaction with the ambient magnetic field. Effects of the solar wind magnetic field on interstellar grains passing through the solar system have been observed. Nanometer sized dust stream particles have been found which were accelerated by Jupiters magnetic field to speeds of about 300 km s^–1.     

Galactic X-rays observed with X-ray astronomy satellite ‘Hakucho’

Highlights of the results obtained with Japanese X-ray astronomy satellite Hakucho are reviewed. After a brief account of instrumentation (Section 2), some new features of non-bursting, non-pulsating objects are presented (Sections 3–5). The main part of the present review is devoted for X-ray bursts which are found more complex than one might have thought (Sections 6–11). The observation of X-ray pulsar, including a change of spin rate of Vela X-1, is described (Section 12). The main results obtained in the first two years are summarized in Section 13.     

Some problems of gamma-astronomy

The modern state of gamma-ray astronomy is reviewed, the paper being mainly devoted to the theoretical models that describe generation of gamma-ray emission under astrophysical conditions. Basic information on the processes of generation and absorption of gamma-rays, as well as the results of observations for various gamma-ray photon energies are reported.In the region of soft gamma-ray emission (i.e., for energies less than tens of MeV), where emission in gamma-ray lines dominates, we also discuss the nature of gamma-ray bursts, the origin of gamma-ray emission from the galactic centre, etc.Discrete sources and, in particular, the mysterious source Cyg X-3 are discussed in the region of very high (E > 10¹² eV) and ultra-high (E > 10¹⁵ eV) energy gamma-ray emission.A larger portion of the review is devoted to the analysis of cosmic-ray origin on the basis of the available gamma-ray data in the region from several tens of MeV to several GeV. The peculiarity of this energy range is, in particular, in the fact that the diffuse galactic emission was observed mainly there. We also discuss the problem of determination of the cosmic-ray density gradient from the gamma-ray data.The origin of high-latitude gamma-ray emission, the problem of galactic gamma-ray halo, etc., are discussed.The theoretical models explaining the nature of unidentified gamma-ray sources, as well as the results of measurements and theoretical estimations of a gamma-ray flux from SN1987A are analysed.List of Notations m electron mass, m = 9.108 × 10^–28 g, - M proton mass, M = 1.672 × 10^–24 g, - e electron charge, e = 4.803 × 10^–10 CGS - c velocity of light, c = 2.9979 × 10¹⁰ cm s^–1, - k Boltzmann constant, k = 1.380 × 10^–16 erg grad^–1, - e electron - p proton - gamma-ray photon - p antiproton - ^{0 0}-meson - -lepton - e ⁺ positron - r, , x radio-frequency, gamma-ray, and X-ray emission bands - E total energy of a particle - E _k kinetic energy - p particle momentum - spectral index for particles - spectral index for emission - n particle density (concentration) - H magnetic field strength - T temperature - _ph energy of low-energy photons - emission frequency - r _H Larmor radius of relativistic particles - k wave number - , z cylindric coordinates, in this case the coordinate (radius) along the galactic disk, z perpendicular to the galactic disk - M solar mass, M = 1.99 × 10³³ g.     

Ultra high-energy gamma-ray observations of Geminga

Results of the observations of Geminga (2CG 195 + 4) in the energy range E 10¹² eV, carried out in 1979, 1981, and 1983 with the Tien Shan high-altitude facility for recording the erenkov flashes of extensive air showers are reported. The mean flux density averaged over the whole protracted data is (5.7 ± 2.5) × 10^–11 quanta cm^–1 s^–1. The flux is variable with a period 59 s. The character of the period variation with time is hard to be reconciled with earlier findings by other authors. The importance of further simultaneous observations at various energies is indicated.     

Cygnus X-1 revisited

We review the X- and gamma-ray observations of Cygnus X-1 and their theoretical interpretations, with emphasis on new developments since the mid-1970's. The overall data base at present is most consistent with the inverse Compton model by hot thermal electrons of T _e 10⁹ K, for the hard X-ray luminosity (10–200 keV). However, the origin of the soft X-rays ( 10 keV) in high states and gamma rays (> 200 keV) remain unsettled.Operated under DOE Contract W-7405-Eng-48.Partially supported by NASA Grant NGR 05-020-668.NRC/NRL Research Associate.