Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal plasma”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 505 records · Page 28

Convective instabilities of electromagnetic ion cyclotron waves in the outer magnetosphere

The path-integrated linear growth of electromagnetic ion cyclotron waves in the outer (L is greater than or equal to 7) magnetosphere is investigated using a realistic thermal plasma distribution with an additional anisotropic energetic ring current H(+) to provide free energy for instability. The results provide a realistic simulation of the recent Active Magneto- spheric Particle Tracer Explorers (AMPTE) observations. For conditions typical of the dayside magnetosphere, high plasma beta effects reduce the group velocity and significantly increase the spatial growth rates for left-handed polarized instabilities just below the helium gyrofrequency Omega(sub He(+)), and on the guided mode above Omega(sub He(+)) but below the cross over frequency omega(sub cr). Relatively high densities, typical of the afternoon local time sector, favor these low group velocity effects for predominantly field-aligned waves. Lower densities, typical of those found in the early morning local time sector, increase the group velocity but allow strong convective instabilities at high normalized frequencies well above Omega(sub He(+)). These waves are reflected in the magnetosphere and can exist for several equatorial transits without significant damping. They are left-handed polarized only on the first equatorial crossing and become linearly polarized for the remainder of the ray path. Consequently, these waves should be observed with basically linear polarization at all frequencies and all latitudes in the early morning local time sector. Wave growth below Omega(sub He(+)) is severely limited owing to the narrow bandwidth for instability and the small resonant path lengths. In the afternoon sector, where plasma densities can exceed 10(exp 7)/cu m, intense convective amplification is possible both above and below Omega(sub He(+)). Waves below Omega(sub He(+)) are not subject to reflection when the O(+) concentration is small and therefore should be observed with left-handed polarization near the equator and essentially linear polarization at higher latitudes. Since the He(+) concentration is usually large in the afternoon sector, guided mode waves above Omega(sub He(+)) reflect to form a background distribution with basically linear polarization. We suggest that the strong left-handed polarized emissions observed by AMPTE in the afternoon sector near the equator are probably due to strongly growing low group velocity waves at frequencies just below Omega(sub He(+)), and on the guided mode above Omega(sub He(+)).

Horne, Richard B.↗

Grant Proposal for the Continuation of the Voyager Interstellar Mission: LECP Investigation

This proposal documents the plans of the Low Energy Charged Particle (LECP) investigation team for participation in NASA's Voyager Interstellar Mission (VIM) as the Voyager 1 and 2 spacecraft explore the outer reaches of the heliosphere and search for the termination shock and the heliopause. The proposal covers the four year period from 1 January 1997 to 31 December 2000. The LECP instruments on Voyager 1 and 2 measure in situ intensities of charged particles with energies from about 30 keV to 100 MeV for ions, and about 20 keV to greater than 10 MeV for electrons. The instruments provide detailed spectral, angular, and compositional information about the particles. Composition is available for greater than 200 keV/nuc using multi-parameter measurements. Angular information is obtained by a mechanically scanned platform that rotates at various commanded rates. Measurements of low energy ion and electron intensities versus time and spatial location within the heliosphere contain an abundance of information regarding various transport and acceleration processes on both local (approx. 1 hr, approx. 0.01 AU) and global (approx. 11 yrs, approx. 100 AU) scales. The LECP instruments provide unique observations of such dynamical processes, and we anticipate that it will return critical information regarding the boundaries of the heliosphere. Several recent and exciting discoveries based on LECP measurements emphasize the important role that low energy charged particle distributions play in physical processes in the interplanetary medium. Yet, at the same time, these discoveries also underscore the fact that our understanding of processes in the outer heliosphere is, in most cases, incomplete, and in others, only rudimentary at best. Among the discoveries referred to above are the following: (1) Shocks: Examination of greater than 30 keV ion intensities have revealed: (a) a total absence of acceleration beyond only -100-200 keV at a strong transient shock in May 1991 at 35 AU, despite an enhanced level of seed particles; (b) a large transient shock in September 1991 of global scale, with intensities of shock-accelerated ions greater than or equal to 30 keV to approx. 30 MeV showing complex, highly energy-dependent spatial evolution, and small-scale (approx. few gyroradii), often anisotropic, micro-structures; (c) recurrent intensity increases in greater than or equal to 30 keV to -few MeV ions, with structures that, in some cases, show no correlation with the associated corotating shock. (2) Superthermal ion pressure: A global merged interaction region with a leading shock, downstream of which the superthermal ion (greater than or equal to 30 keV to approx. 4 MeV) pressure is comparable to that of the thermal plasma, and the total particle pressure yields a plasma beta of order unity. (3) Pickup ions: Measurements of the C/O ratio within transient structures at 35-45 AU showing the first clear evidence that transient shocks can pre-accelerate interstellar pickup ions from approx. 1 keV/nuc to at least 1 MeV/nuc. (4) Seed particles: Injection of ions for acceleration to high energies at the termination shock is unlikely to be a problem, since interplanetary transient and recurrent shocks are continually accelerating ions, of solar wind or interstellar origin, to highly superthermal energies. (5) Precursor electrons: Ambient solar electrons (greater than or equal to few tens of keV) that exist in the outer heliosphere ca form a broad precursor, several days wide, that is upstream of the termination shock and potentially observable a few months prior to the shock crossing. (6) Solar wind velocity at Voyager 1: We can use LECP ion data to obtain the solar wind velocity at Voyager 1, enabling us to provide critical measurement of the plasma flow as we approach and encounter the termination shock and other regions (necessary due to the partial failure of the Voyager 1 PLS experiment). The work of the LECP investigator team during the VIM will include: (1) Continuing operations with regard to the receipt, processing, verification, cataloging, display, and distribution of the data from the LECP instruments on Voyager 1 and 2, (2) Monitoring the health and performance of the LECP instruments, and evaluating and characterizing the response of the LECP instruments to various energetic particle and plasma environments, (3) Participating in, and supporting Voyager Project planning exercises and other coordinated activities relevant to exploration of the outer heliosphere, (4) Developing analysis techniques and operational procedures suitable for searching for and characterizing the boundaries and unique regions of the outher heliosphere, (5) Continuing the preparation of data sets appropriate for submission to the National Space Sciences Data Center (NSSDC) and, where appropriate, the Planetary Data System (PDS), (6) Maintaining direct Web access to online LECP data through the JHU/APL Voyager LECP home page, (7) Performing scientific evaluations of the Voyager 1 and 2 LECP data sets in conjunction with other data sets and other investigators, with particular focus on the outer regions of the heliosphere, and (8) Publishing the results of these evaluations in the scientific literature and presenting the results in scientific conferences.

Krimigis, Stamatios M.↗

Empirical Modeling of the Plasmasphere

Over 40 years of ground and spacecraft plasmaspheric measurements have resulted in many statistical descriptions of plasmaspheric properties. In some cases, these properties have been represented as analytical descriptions that are valid for specific regions or conditions. For the most part, what has not been done is to extend regional empirical descriptions or models to the plasmasphere as a whole. In contrast, many related investigations depend on the use of representative plasmaspheric conditions throughout the inner magnetosphere. Wave propagation, involving the transport of energy through the magnetosphere, is strongly affected by thermal plasma density and its composition. Ring current collisional and wave particle losses also strongly depend on these quantities, Plasmaspheric also plays a secondary role in influencing radio signals from the Global Positioning System satellites. The Global Core Plasma Model (GCPM) is an attempt to assimilate previous empirical evidence and regional models for plasmaspheric density into a continuous, smooth model of then-nal plasma density in the inner magnetosphere. In that spirit, the International Reference Ionosphere is currently used to complete the low altitude description of density and composition in the model. The models and measurements on which the GCPM is currently based and its relationship to IRI will be discussed.

Gallagher, Dennis L.↗

Iron K Line Variability in the Low-Luminosity AGN NGC 4579

We present results of new ASCA observations of the low-luminosity AGN (LLAGN) NGC 4579 obtained in 1998 December 18 and 28, and we report on detection of variability of an iron K emission line. The X-ray luminosities in the 2-10 keV band for the two observations are nearly identical, L(sub X) approximately = 2 x 10(exp 4l) ergs/s, but they are approximately 35% larger than that measured in 1995 July by Terashima et al. (1998). An Fe K emission line is detected at 6.39 +/- 0.09 keV (source rest frame) which is lower than the line energy 6.73(sup +0.13, sub -0.12) keV in the 1995 observation. If we fit the Fe lines with a blend of two Gaussians centered at 6.4 keV and 6.73 KeV, the intensity of the 6.7 keV line decreased, while the intensity of the 6.4 keV line increased, within an interval of 3.5 years. This variability rules out thermal plasmas in the host galaxy as the origin of the ionized Fe line in this LLAGN. The detection and variability of the 6.4 keV line suggest that an optically thick standard accretion disk is present and subtends a large solid angle viewed from the nucleus at the Eddington ratio of L(sub Bol)/L(sub Eddington) approximately 2 x 10(exp -3) (Ho 1999). A broad disk-line profile is not clearly seen and the structure of the innermost part of accretion disk remains unclear.

Terashima, Yuichi↗

The Plasmasphere as Seen by the IMAGE EUV Instrument

The IMAGE spacecraft includes the extreme ultraviolet (EUV) imager, which is designed to view He+ plasmaspheric ions in the inner magnetosphere. These images provide the first opportunity to view the global distribution of plasmaspheric ions and their response to solar wind induced convection and the ring current. In addition to the generalized pattern of convection long believed to develop in thermal plasma, the EUV instrument is observing other meso-scale patterns of plasma distribution never before observed or postulated using in situ measurements. This poster presentation will present the "first light" observations of the EUV instrument and their possible interpretations.

Gallagher, D. L.↗

Transport of Photoelectrons in the Nightside Magnetosphere

Kinetic modeling results are analyzed to examine the transport of photoelectrons through the nightside inner magnetosphere. A natural filter exists on the nightside for these particles, with Coulomb collisions eroding the distribution at low energies and low L shells, and magnetospheric convection compressing the flux tubes as they corotate towards dawn. This leads to an electron population with very distinctive features. It is shown that, for low activity levels, a band of photoelectrons forms between L=4 and 6 that extends throughout the nightside local times and into the morning sector. Modeled energy and pitch angle distributions are presented to show the development and features of this electron band, as well as their influences on the thermal plasma. Such a population should be observable in satellite spectrometer measurements.

Khazanov, G. V.↗

The Storm-Time Plasmasphere by IMAGE/EUV

With the availability of routine global images of the plasmasphere from the IMAGE Mission Extreme Ultraviolet Imager (EUV), we now have a growing body of observations that both characterize and quantify the dynamics of the plasmasphere. Direct interactions of the plasmasphere with the ring current and indirect interaction through localized electric fields appear to be the cause of a complex array of structures in the spatial distribution of thermal plasma. This presentation will show those features in detail and discuss possible mechanisms for their cause.

Gallagher, D. L.↗

Extreme Convection Conditions for the Plasmasphere

The IMAGE Extreme Ultraviolet (EUV) imager has now observed the plasmasphere under conditions of extreme erosion. Surprisingly, the plasmasphere is sometimes found to almost disappear. Global EUV images are used together with dynamic plasmasphere modeling to describe the convection electric field necessary to produce observed thermal plasma distributions under extreme conditions. These results will be compared with established measures of subauroral electric fields.

Gallagher, D. L.↗

Understanding the Long-Term Spectral Variability of Cygnus X-1 from BATSE and ASM Observations

We present a spectral analysis of observations of Cygnus X-1 by the RXTE/ASM (1.5-12 keV) and CGRO/BATSE (20-300 keV), including about 1200 days of simultaneous data. We find a number of correlations between intensities and hardnesses in different energy bands from 1.5 keV to 300 keV. In the hard (low) spectral state, there is a negative correlation between the ASM 1.5-12 keV flux and the hardness at any energy. In the soft (high) spectral state, the ASM flux is positively correlated with the ASM hardness (as previously reported) but uncorrelated with the BATSE hardness. In both spectral states, the BATSE hardness correlates with the flux above 100 keV, while it shows no correlation with the flux in the 20-100 keV range. At the same time, there is clear correlation between the BATSE fluxes below and above 100 keV. In the hard state, most of the variability can be explained by softening the overall spectrum with a pivot at approximately 50 keV. The observations show that there has to be another, independent variability pattern of lower amplitude where the spectral shape does not change when the luminosity changes. In the soft state, the variability is mostly caused by a variable hard (Comptonized) spectral component of a constant shape superimposed on a constant soft blackbody component. These variability patterns are in agreement with the dependence of the rms variability on the photon energy in the two states. We interpret the observed correlations in terms of theoretical Comptonization models. In the hard state, the variability appears to be driven mostly by changing flux in seed photons Comptonized in a hot thermal plasma cloud with an approximately constant power supply. In the soft state, the variability is consistent with flares of hybrid, thermal/nonthermal, plasma with variable power above a stable cold disk. Also, based on broadband pointed observations simultaneous with those of the ASM and BATSE, we find the intrinsic bolometric luminosity increases by a factor of approximately 3-4 from the hard state to the soft one, which supports models of the state transition based on a change of the accretion rate.

Zdziarski, Andrzej A.↗

An XMM-Newton Observation of the Lagoon Nebula and the Very Young Open Cluster NGC 6530

We report the results of an XMM-Newton observation of the Lagoon Nebula (M 8). Our EPIC images of this region reveal a cluster of point sources, most of which have optical counterparts inside the very young open cluster NGC 6530. The bulk of these X-ray sources are probably associated with low and intermediate mass pre-main sequence stars. One of the sources experiences a flare-like increase of its X-ray flux making it the second brightest source in M 8 after the O4 star 9 Sgr. The X-ray spectra of most of the brightest sources can be fitted with thermal plasma models with temperatures of kT approximately a few keV. Only a few of the X-ray selected PMS candidates are known to display H(alpha) emission and were previously classified as classical T Tauri stars. This suggests that most of the X-ray emitting PMS stars in NGC 6530 are weak line T Tauri stars. In addition to 9 Sgr, our EPIC field of view contains also a few early-type stars. The X-ray emission from HD 164816 is found to be typical for an O9.5III-IV star. At least one of the known Herbig Be stars in NGC 6530 (LkH(alpha) 115) exhibits a relatively strong X-ray emission, while most of the main sequence stars of spectral type B1 and later are not detected. We also detect (probably) diffuse X-ray emission from the Hourglass Region that might reveal a hot bubble blown by the stellar wind of Herschel 36, the ionizing star of the Hourglass Region.

Rauw, G.↗

Understanding the Long-Term Spectral Variability of Cygnus X-1 with Burst and Transient Source Experiment and All-Sky Monitor Observations

We present a comprehensive analysis of all observations of Cyg X-1 by the Compton Gamma Ray Observatory Burst and Transient Source Experiment (BATSE; 20-300 keV) and by the Rossi X-Ray Timing Explorer all-sky monitor (ASM; 1.5-12 keV) until 2002 June, including approximately 1200 days of simultaneous data. We find a number of correlations between fluxes and hardnesses in different energy bands. In the hard (low) spectral state, there is a negative correlation between the ASM 1.5-12 keV flux and the hardness at any energy. In the soft (high) spectral state, the ASM flux is positively correlated with the ASM hardness but uncorrelated with the BATSE hardness. In both spectral states, the BATSE hardness correlates with the flux above 100 keV, while it shows no correlation with the 20-100 keV flux. At the same time, there is clear correlation between the BATSE fluxes below and above 100 keV. In the hard state, most of the variability can be explained by softening the overall spectrum with a pivot at approximately 50 keV. There is also another, independent variability pattern of lower amplitude where the spectral shape does not change when the luminosity changes. In the soft state, the variability is mostly caused by a variable hard (Comptonized) spectral component of a constant shape superposed on a constant soft blackbody component. These variability patterns are in agreement with the dependencies of the rms variability on the photon energy in the two states. We also study in detail recent soft states from late 2000 until 2002. The last of them has lasted thus far for more than 200 days. Their spectra are generally harder in the 1.5-5 keV band and similar or softer in the 3-12 keV band than the spectra of the 1996 soft state, whereas the rms variability is stronger in all the ASM bands. On the other hand, the 1994 soft state transition observed by BATSE appears very similar to the 1996 one. We interpret the variability patterns in terms of theoretical Comptonization models. In the hard state, the variability appears to be driven mostly by changing flux in seed photons Comptonized in a hot thermal plasma cloud with an approximately constant power supply. In the soft state, the variability is consistent with flares of hybrid, thermal/nonthermal, plasma with variable power above a stable cold disk. The spectral and timing differences between the 1996 and 2000-2002 soft states are explained by a decrease of the color disk temperature. Also, on the basis of broadband pointed observations simultaneous with those of the ASM and BATSE, we find the intrinsic bolometric luminosity increases by a factor of approximately 3-4 from the hard state to the soft one, which supports models of the state transition based on a change of the accretion rate.

Zdziarski, Andrzej A.↗

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The May 2-7, 1998, Storm

Complete description of a self-consistent model for magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves, and back on waves, are considered self-consistently by solving both equations on a global magnetospheric scale under non steady-state conditions. In the paper by Khazanov et al. [2002] this self-consistent model has only been shortly outlined, and discussions of many the model related details have been omitted. For example, in present study for the first time a new algorithm for numerical finding of the resonant numbers for quasilinear wave-particle interaction is described, or it is demonstrated that in order to describe quasilinear interaction in a multi-ion thermal plasma correctly, both e and He(+) modes of electromagnetic ion cyclotron waves should be employed. The developed model is used to simulate the entire May 2-7, 1998 storm period. Trapped number fluxes of the ring current protons are calculated and presented along with their comparison with the data measured by the 3D hot plasma instrument Polar/HYDRA. Examining of the wave (MLT, L shell) distributions produced during the storm progress reveals an essential intensification of the wave emissions in about two days after main phase of storm. This result is well consistent with the earlier ground-based observations. Also the theoretical shapes and the occurrence rates for power spectral densities of electromagnetic ion cyclotron waves are studied. It is found that in about 2 days after the storm main phase on May 4, mainly non Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗

RHESSI Results on Energetic Particles in the Solar Atmosphere

Since launch on 5 February, 2002, the Reuven Ramaty High Energy Solar Spectroscopic Imager has recorded well over 10,000 X-ray flares and at least 6 gamma-ray line flares. I will review some of the major results of the analysis of a small fraction of these observations of particular interest to solar radio astronomers. In particular, I will present X-ray evidence for the generation of current sheets during impulsive X-ray flares and for magnetic reconnection as the main energy release mechanism. I will also demonstrate the close association between several large flares and the accompanying coronal mass ejections. The energetic importance of accelerated electrons and ions as compared to the total energy appearing in the thermal plasma and CME mass motions will be discussed. Possible explanations will be presented for the spatial separation between the hard X-ray footpoints and the centroid of the gamma-ray line source in the July 23,2002, gamma-ray line flare.

Dennis, Brian R.↗

Multi-Epoch XMM Observations of NGC4258

The goal of this project was detection of variability in the X-ray absorption column of the AGN in NGC4258 through monitoring with the XMM satellite. We have accomplished this goal and submitted the results to ApJ for publication in a paper entitled, "X-ray Luminosity and Absorption Column Fluctuations in the H2O Maser Galaxy NGC4258 from Weeks to Years," by Fruscione, A., Greenhill, L.J., Filippedco, A.V., Moran, J.M., Hermstein, J.R., and Galle, E. We have received a favorable referee report and expect the article will appear in 2005. To complete the project, we reduced our four epochs of XMM data for NGC4258, one archival XMM observation, and all existing Chandra datasets for NGC4258 (with the latest calibration tables and a grid of corrections for pileup). Self-consistent reduction of all these data permitted detailed comparison that could not have been accomplished simply by taking published model fits that appear in the literature. To accumulate a broader monitoring record, we combined the Chandra and XMM results with those published for SAX and ASCA. We modeled the Chandra and XMM data self-consistently with partially absorbed, hard power-law, soft thermal plasmas, and soft power-law components. Over nine years, the photo-electric absorbing column exhibited a 40% drop between two ASCA epochs separated by 3 years and a 60% rise between two XMM epochs separated by just 5 months. In contract, uncorrelated factor of of 2-3 changes were seen in absorbed flux on te timescale of years, which suggests intrinsic variability of the central engine. The warped disk that is a known source of H2O maser emission in 4258 is believed to cross the line of sight to the central engine. We have proposed that the variations in absorbing column arise from inhomogeneities in the rotating disk, as they sweep across the line of sight. We estimate from the XMM data that the inhomogeneities are about 1E+15 cm in size at radii greater than 0.27 pc. This is consistent with the estimated radius of the disk crossing estimated entirely independently, and it solidifies evidence that the warped accretion disk is the absorber in this (and possible other) AGN.

Greenhill, Lincoln J.↗

The Physics of Cooling Flow Clusters with Central Radio Sources

Central galaxies in rich clusters are the sites of cluster cooling flows, with large masses of gas cooling through part of the X-ray band. Many of these galaxies host powerful radio sources. These sources can displace and compress the X-ray gas leading to enhanced cooling and star formation. We observed the bright cooling flow Abell 2626 with a strangely distorted central radio source. We wished to understand the interaction of radio and X-ray thermal plasma, and to determine the dynamical nature of this cluster. One aim was to constrain the source of additional pressure in radio "holes" in the X-ray emission needed to support overlying shells of X-ray gas. We also aimed to study the problem of the lack of kT < 1-2 keV gas in cooling flows by searching for abundance inhomogeneities, heating from the radio source, and excess absorption. We also have a Chandra observation of this cluster. There were problems with the pipeline processing of this data due to a telemetry dropout. We are publishing the Chandra and XMM data together. Delays with the Chandra data have slowed up the publication. At the center of the cluster, there is a complex interaction of the odd, Z-shaped radio source, and the X-ray plasma. However, there are no clear radio bubbles. Also, the cluster SO galaxy IC 5337, which is projected 1.5 arcmin west of the cluster center, has unusual tail-like structures in both the radio and X-ray. It appears to be falling into the cluster center. There is a hot, probably shocked region of gas to the southwest, which is apparently due to the merger of a subcluster in this part of the system. There is also a merging subcluster to the northeast. The axes of these two mergers agrees with a supercluster filament structure.

Sarazin, Craig L.↗

Jet-Intracluster Medium Interactions of the Head Tail Radio Galaxy 3C 129

The 50 ksec XMM observations of the galaxy cluster 3C 129 were taken as scheduled, and the data are of good quality. We analyzed the data in the following way. After standard cleaning, we flat-fielded the XMM surface brightness maps. Combining the data from the EPIC MOS and PN Camera CCDs, we performed a cross-correlation analysis of the X-ray surface brightness distribution with the 1.4 GHz VLA radio map. We found evidence for cavities in the X-ray emitting Intra-Cluster Medium (ICM) associated with the radio tail of the head-tail radio galaxy 3C 129. This discovery is very interesting as it excludes the presence of a large fraction of thermal plasma in the radio tail. Together with the observation of an apparent pressure mismatch between the radio plasma and the ICM, and an upper limit on the magnetic field inside the radio tail (from the radio spectral indices map) the observation implies that the tail pressure is dominated either by low-energy electrons/positrons, or, by relativistic protons. Furthermore, we studied the energy spectrum of an X-ray "hot-spot" associated with the head of the radio galaxy 3C 129. It seems likely that the X-ray hot-spot originates from shocked gas in front of the radio galaxy. , The analysis turned out to be much more difficult than anticipated. The main reason is the lack of a comprehensive, publicly available background model that is key for the analysis of extended sources. Small groups like our do not have the man-power to come up with a background model themselves. We used the model from Read & Ponman (A&A 409, 395, 2003). However, the background subtracted X-ray surface brightness maps show a bright ring in the outer 20% of the camera. We tried to get rid of this ring and contacted the XMM helpdesk and Read & Ponman, the authors of the background paper. However, up to this day, we did not entirely succeed to remove the brightness enhancement at the outer parts of the camera. Unfortunately, our results are somewhat sensitive to the uncertainty, as the radio galaxy 3C 129 is very large and occupies a rather large fraction of the XMM field of view. We are now working on a paper describing the results. The paper will include a detailed discussion of the uncertainties associated with the non- perfect background subtraction.

Krawczynski, Henric S. W.↗

Influence of Coronal Abundance Variations

The PI of this project was Jeff Scargle of NASA/Ames. Co-I's were Alma Connors of Eureka Scientific/Wellesley, and myself. Part of the work was subcontracted to Eureka Scientific via SAO, with Vinay Kashyap as PI. This project was originally assigned grant number NCC2-1206, and was later changed to NCC2-1350 for administrative reasons. The goal of the project was to obtain, derive, and develop statistical and data analysis tools that would be of use in the analyses of high-resolution, high-sensitivity data that are becoming available with new instruments. This is envisioned as a cross-disciplinary effort with a number of "collaborators" including some at SA0 (Aneta Siemiginowska, Peter Freeman) and at the Harvard Statistics department (David van Dyk, Rostislav Protassov, Xiao-li Meng, Epaminondas Sourlas, et al). We have developed a new tool to reliably measure the metallicities of thermal plasma. It is unfeasible to obtain high-resolution grating spectra for most stars, and one must make the best possible determination based on lower-resolution, CCD-type spectra. It has been noticed that most analyses of such spectra have resulted in measured metallicities that were significantly lower than when compared with analyses of high- resolution grating data where available (see, e.g., Brickhouse et al., 2000, ApJ 530,387). Such results have led to the proposal of the existence of so-called Metal Abundance Deficient, or "MAD" stars (e.g., Drake, J.J., 1996, Cool Stars 9, ASP Conf.Ser. 109, 203). We however find that much of these analyses may be systematically underestimating the metallicities, and using a newly developed method to correctly treat the low-counts regime at the high-energy tail of the stellar spectra (van Dyk et al. 2001, ApJ 548,224), have found that the metallicities of these stars are generally comparable to their photospheric values. The results were reported at the AAS (Sourlas, Yu, van Dyk, Kashyap, and Drake, 2000, BAAS 196, v32, #54.02), and at the conference on Statistical Challenges in Modem Astronomy (Sourlas, van Dyk, Kashyap, Drake, and Pease, 2003, SCMA 111, Eds. E.D.Feigelson, G.J.Babu, New York:Springer, p489-490). We also described the limitations of one of the most egregiously misused and misapplied statistical tests in astrophysical literature, the F-test for verifying model components (Protassov, van Dyk, Connors, Kashyap, and Siemiginowska, 2002, ApJ, 571,545). Indeed, a search through the ApJ archives turned up 170 papers in the 5 previous years that used the F-test explicitly in some form or the other, and with the vast majority of them not using it correctly! Indeed, looking at just 4 issues of the ApJ in 2001, we found 13 instances of its use, of which nine were demonstrably incorrect. Clearly, it is difficult to understate the importance of this issue. We also worked on speeding up Bayes Blocks and Sparse Bayes Blocks algorithms to make them more tractable for large searches. We also supported staistics students and postdocs in both explicit physics- model-based (spectra with tens of thousands of atomic lines) and "model-free" -- i.e. non-parametric or semi-parametric -- algorithms. Work on using more of the latter is just beginning; while using multi-scale methods for Poisson imaging has come to hition. In fact, "An Image Restoration Technique with Error Estimates", by D. Esch, A. Connors, M. Karovska, and D. van Dyk, was published by ApJ (Esch et a1.2004, ApJ, 610, 1213). The code has been delivered to M. Karovska for CXC; and is available for beta-testing upon request. The other large project we worked on was on the self-consistent modeling of logN-logs curves in the Poisson limit. logN-logs curves are a fundamental tool in the study of source populations, luminosity functions, and cosmological parameters. However, their determination is hampered by statistical effects such as the Eddington bias, incompleteness due to detection efficiency, faint source flux fluctuations, etc. We have develed a new and powerful method using the full Poisson machinery that allows us to model the logN-logs distribution of X-ray sources in a self-consistent manner. Because we properly account for all the above statistical effects, our modeling is valid over the full range of the data, and not just for strong sources, as is normally done. Using a Bayesian approach and modeling the fluxes with known functional forms such as simple or broken power-laws, and conditioning the expected photon counts on the fluxes, the background contamination, effective area, detector vignetting, and detection probability, we can delve deeply into the low counts regime and extend the usefulness of medium sensitivity surveys such as ChAMP by orders of magnitude. The built-in flexibility of the algorithm also allows a simultaneous analysis of multiple datasets. We have applied this analysis to a set a Chandra observations (Sourlas, Kashyap, Zezas, van Dyk, 2004, HEAD #8, #16.32)

Scargle, Jeffrey D.↗

X-Ray and Infrared Observations of Embedded Young Stars in NGC 2264

Images of the NGC 2264 star-forming region, which we have acquired with the XMM-Newton spacecraft, reveal strong X-ray emission from three deeply embedded (Av > 10 mag) young stellar objects in the vicinity of Allen's infrared source (AFGL 989 = IRS 1) and Castelaz & Grasdalen s infrared source (RNO-EW = IRS 2). Thermal plasma models for the brightest source in X-rays, located 11 southwest of Allen's star, yield a quasi-steady luminosity of Lx = 10 ergs s-1 and an extraordinarily high X-ray temperature of 100 MK. The high temperature is consistent with the presence of emission lines of Fe xxv and Fe xxvi at photon energies of 6.7 and 6.9 keV, respectively. An even higher temperature of nearly 140 MK was observed during the rise phase of a powerful impulsive X-ray flare of another young star in the IRS 2 region. Moderate-resolution near-infrared (1-4 um) spectra of the embedded objects, obtained at the NASA Infrared Telescope Facility, exhibit deep water ice absorption bands, as well as a variety of emission and absorption features of H I, CO, and both neutral and ionized metals.

Simon, Theordore↗