ALBERT

Description: Two decades ago, Witten suggested that the ground state of matter might be material of nuclear density made from up, down and strange quarks. Since then, much effort has gone into exploring astrophysical and other implications of this possibility. For example, neutron stars would almost certainly be strange quark stars; dark matter might be strange quark matter. Searches for stable strange quark matter have been made in various mass ranges, with negative, but not conclusive results. Recently, we [D. Anderson, E. Herrin, V. Teplitz, and I. Tibuleac, Bull. Seis. Soc. of Am. 93, 2363 (2003)] reported a positive result for passage through the Earth of a multi-ton "nugget" of nuclear density in a search of about a million seismic reports, to the U.S. Geological Survey for the years 1990-93, not associated with known Earthquakes. I will present the evidence (timing of first signals to the 9 stations involved, first signal directions, and unique waveform characteristics) for our conclusion and discuss potential improvements that could be obtained from exploiting the seismologically quieter environments of the moon and Mars.

Description: Many papers have been published that further elucidate the structure of coronas in cool stars as determined from EUVE, HST, FUSE, Chandra, and XMM-Newton observations. In addition we are exploring the effects of coronas on the He I 1083081 transition that is observed in the infrared. Highlights of these are summarized below including publications during this reporting period and presentations. Ground-based magnetic Doppler imaging of cool stars suggests that active stars have active regions located at high latitudes on their surface. We have performed similar imaging in X-ray to locate the sites of enhanced activity using Chandra spectra. Chandra HETG observations of the bright eclipsing contact binary 44i Boo and Chandra LETG observations for the eclipsing binary VW Cep show X-ray line profiles that are Doppler-shifted by orbital motion. After careful analysis of the spectrum of each binary, a composite line-profile is constructed by adding the individual spectral lines. This high signal-to-noise ratio composite line-profile yields orbital velocities for these binaries that are accurate to 30 km/sec and allows their orbital motion to be studied at higher time resolutions. In conjunction with X-ray lightcurves, the phase-binned composite line-profiles constrain coronal structures to be small and located at high latitudes. These observations and techniques show the power of the Doppler Imaging Technique applied to X-ray line emission.

Description: Dynamic signal fluctuations due to atmospheric scintillations may impair the Ka-band (around 32-GHz) link sensitivities for a low-margin Deep Space Network (DSN) receiving system. The ranges of frequency and power of the fast fluctuating signals (time scale less than 1 min) are theoretically investigated using the spatial covariance and turbulence theory. Scintillation power spectrum solutions are derived for both a point receiver and a finite-aperture receiver. The aperture-smoothing frequency ((omega(sub s)), corner frequency ((omega(sub c)), and damping rate are introduced to define the shape of the spectrum for a finite-aperture antenna. The emphasis is put on quantitatively describing the aperture-smoothing effects and graphically estimating the corner frequency for a large aperture receiver. Power spectral shapes are analyzed parametrically in detail through both low- and high-frequency approximations. It is found that aperture-averaging effects become significant when the transverse correlation length of the scintillation is smaller than the antenna radius. The upper frequency or corner frequency for a finite-aperture receiver is controlled by both the Fresnel frequency and aperture-smoothing frequency. Above the aperture-smoothing frequency, the spectrum rolls off at a much faster rate of exp (-omega(sup 2)/omega(sup 2, sub s), rather than omega(sup -8/3), which is customary for a point receiver. However, a relatively higher receiver noise level can mask the fast falling-off shape and make it hard to be identified. We also predict that when the effective antenna radius a(sub r) less than or = 6 m, the corner frequency of its power spectrum becomes the same as that for a point receiver. The aperture-smoothing effects are not obvious. We have applied these solutions to the scenario of a DSN Goldstone 34-m-diameter antenna and predicted the power spectrum shape for the receiving station. The maximum corner frequency for the receiver (with omega(sub s) = 0.79 omega(sub 0) is found to be 0.44 Hz (or 1.0 omega(sub 0), while the fading rate (or fading slope) is about 0.06 dB/s.

Description: We present the far-ultraviolet spectrum of the Seyfert 1 galaxy 2MASX J21362313-6224008 obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE). The spectrum features absorption from Galactic O VI at two velocities and redshifted H I Ly beta and gamma, C II, CIII, and O VI. The redshifted absorption features represent a single kinematic component blueshifted by approx. 310 km/s relative to the active galactic nucleus. We use photoionization models to derive constraints on the physical parameters of the absorbing gas. An alternative interpretation for the absorption lines is also proposed, wherein the absorbing gas is associated with an intervening galaxy cluster.

Description: Using SIM, we will perform narrow-angle observations of several X-ray binaries to determine their orbits, and we will observe about 50 X-ray binary systems in wide-angle mode to measure their distances and proper motions. Sources with mass estimates for the compact component of greater than 3 solar masses are generally called black hole candidates since this mass is above the theoretical neutron star limit. Narrow-angle observations of these sources provide a direct test of the dynamical mass estimates on which the black hole evidence is based. Better measurements of the black hole masses will provide constraints on possible evolutionary paths that lead to black hole formation. When combined with X-ray data, mass measurements may provide additional constraints on the black hole spin. Precise mass determinations of neutron star systems can address the question of whether neutron stars can be significantly more massive than 1.4 solar masses, which would eliminate soft models of the neutron star equations of state. The wide-angle observations will probe the Galactic distribution of X-ray binaries through parallaxes and proper motions. They will also eliminate the uncertainties in the luminosities of individual sources, which is currently up to a full order of magnitude. This will enable more detailed comparisons of X-ray observations to physical models such as advection-dominated accretion flows (ADAFs). We intend to carry out the following measurements: 1) Determine the orbits of two black hole candidates to measure the black hole masses; 2) Obtain precise mass measurements for two neutron star systems to constrain neutron star equations of state; 3) Determine the distances and thus luminosities of selected representatives of various classes of X-ray binaries (black hole candidates, neutron stars, jet sources); 4) In the process of distance determination, proper motions will also be measured, from which the age of the population can be estimated.

Description: We report a theoretical investigation of the robustness of two-dimensional inviscid MHD flows at low magnetic Reynolds numbers with respect to three-dimensional perturbations. We analyze three model problems, namely flow in the interior of a triaxial ellipsoid, an unbounded vortex with elliptical streamlines, and a vortex sheet parallel to the magnetic field. We demonstrate that motion perpendicular to the magnetic field with elliptical streamlines becomes unstable with respect to the elliptical instability once the velocity has reached a critical magnitude whose value tends to zero as the eccentricity of the streamlines becomes large. Furthermore, vortex sheets parallel to the magnetic field, which are unstable for any velocity and any magnetic field, are found to emit eddies with vorticity perpendicular to the magnetic field and with an aspect ratio proportional to N(sup 1/2). The results suggest that purely two-dimensional motion without Joule energy dissipation is a singular type of flow which does not represent the asymptotic behaviour of three-dimensional MHD turbulence in the limit of infinitely strong magnetic fields.

Description: Observations of the multi-TeV spectra of the Mkn 501 and other nearby BL Lac objects exhibit the high energy cutoffs predicted to be the result of intergalactic annihilation interactions, primarily with IR photons having a flux level as determined by various astronomical observations. After correcting for such intergalactic absorption, these spectra can be explained within the framework of synchrotron self-Compton emission models. Stecker and Glashow have shown that the existence of this annihilation via electron-positron pair production puts strong constraints on Lorentz invariance violation. Such constraints have important implications for some quantum gravity and large extra dimension models. A much smaller amount of Lorentz invariance violation has potential implications for understanding the spectra of ultrahigh energy cosmic rays.

Description: We report the discovery of a new X-ray pulsar, XTE J1810-197, that was serendipitously discovered on 2003 July 15 by the Rossi X-Ray Timing Explorer (RXTE) while observing the soft gamma repeater SGR 1806-20. The pulsar has a 5.54 s spin period, a soft X-ray spectrum (with a photon index of approx. = 4). and is detectable in earlier RXTE observations back to 2003 January but not before. These show that a transient outburst began between 2002 November 17 and 2003 January 23 and that the source's persistent X-ray flux has been declining since then. The pulsar exhibits a high spin-down rate P approx.= l0(exp -11) s/s with no evidence of Doppler shifts due to a binary companion. The rapid spin-down rate and slow spin period imply a supercritical characteristic magnetic field B approx. = 3 x l0(exp 14) G and a young age tau less than or = 7600 yr. Follow-up Chandra observations provided an accurate position of the source. Within its error radius, the 1.5 m Russian-Turkish Optical Telescope found a limiting magnitude R(sub c) = 21.5. All such properties are strikingly similar to those of anomalous X-ray pulsars ad soft gamma repeaters, providing strong evidence that the source is a new magnetar. However, archival ASCA and ROSAT observations found the source nearly 2 orders of magnitude fainter. This transient behavior and the observed long-term flux variability of the source in absence of an observed SGR-like burst activity make it the first confirmed transient magnetar and suggest that other neutron stars that share the properties of XTE 51810- 197 during its inactive phase may be unidentified transient magnetars awaiting detection via a similar activity. This implies a larger population of magnetars than previously surmised and a possible evolutionary connection between magnetars and other neutron star families. Subject headings: pulsars: general -pulsars: individual (XTE 51810- 197) - stars: magnetic fields -

Description: V407 Vu1 (RX J1914.4+2456) is a candidate double-degenerate binary with a putative 1.756 mHz (9.5 min) orbital frequency. In a previous timing study using archival ROSAT and ASCA data we reported evidence for an increase of this frequency at a rate consistent with expectations for gravitational radiation from a detached ultracompact binary system. Here we report the results of new Chandra timing observations which confirm the previous indications of spin-up of the X-ray frequency, and provide much tighter constraints on the frequency derivative, nu (raised dot). We obtained with Chandra a total of 90 ksec of exposure in two epochs separated in time by 11.5 months. The total time span of the archival ROSAT, ASCA and new Chandra data is now approximately equal to 10.5 years. This more than doubles the interval spanned by the ROSAT and ASCA data alone, providing much greater sensitivity to a frequency derivative. With the addition of the Chandra data an increasing frequency is unavoidable, and the mean nu (raised dot) is 7.0 plus or minus 0.8 x l0(exp -18) Hz per second. Although a long-term spin-up trend is confirmed, there is excess variance in the phase timing residuals, perhaps indicative of shorter timescale torque fluctuations or phase instability associated with the source of the X-ray flux. Power spectral searches for periods longward of the 9.5 minute period do not find any significant modulations, however, the sensitivity of searches in this frequency range are somewhat compromised by the dithering of the Chandra attitude. The observed spin-up is of a magnitude consistent with that expected from gravitational radiation decay, however, the factor of approximately equal to 3 variations in flux combined with the timing noise could conceivably result from accretion-induced spin-up of a white dwarf. Continued monitoring to explore correlations of torque with X-ray flux could provide a further test of this hypothesis.

Description: We report on the Chandra observations of the archetypical mixed morphology (or thermal composite) supernova remnant, W44. As with other mixed morphology remnants, W44's projected center is bright in thermal X-rays. It has an obvious radio shell, but no discernable X-ray shell. In addition, X-ray bright knots dot W44's image. The spectral analysis of the Chandra data show that the remnant s hot, bright projected center is metal-rich and that the bright knots are regions of comparatively elevated elemental abundances. Neon is among the affected elements, suggesting that ejecta contributes to the abundance trends. Furthermore, some of the emitting iron atoms appear to be underionized with respect to the other ions, providing the first potential X-ray evidence for dust destruction in a supernova remnant. We use the Chandra data to test the following explanations for W44's X-ray bright center: 1.) entropy mixing due to bulk mixing or thermal conduction, 2.) evaporation of swept up clouds, and 3.) a metallicity gradient, possibly due to dust destruction and ejecta enrichment. In these tests, we assume that the remnant has evolved beyond the adiabatic evolutionary stage, which explains the X-ray dimness of the shell. The entropy mixed model spectrum was tested against the Chandra spectrum for the remnant's projected center and found to be a good match. The evaporating clouds model was constrained by the finding that the ionization parameters of the bright knots are similar to those of the surrounding regions. While both the entropy mixed and the evaporating clouds models are known to predict centrally bright X-ray morphologies, their predictions fall short of the observed brightness gradient. The resulting brightness gap can be largely filled in by emission from the extra metals in and near the remnant's projected center. The preponderance of evidence (including that drawn from other studies) suggests that W44's remarkable morphology can be attributed to dust destruction and ejecta enrichment within an entropy mixed, adiabatic phase supernova remnant. The Chandra data prompts a new question - by what astrophysical mechanisms are the metals distributed so inhomogeneously in the supernova remnant.

Description: The Theoretical Astrophysics Group works on a broad range of topics ranging from string theory to data analysis in the Sloan Digital Sky Survey. The group is motivated by the belief that a deep understanding of fundamental physics is necessary to explain a wide variety of phenomena in the universe. During the three years 2001-2003 of our previous NASA grant, over 120 papers were written; ten of our postdocs went on to faculty positions; and we hosted or organized many workshops and conferences. Kolb and collaborators focused on the early universe, in particular and models and ramifications of the theory of inflation. They also studied models with extra dimensions, new types of dark matter, and the second order effects of super-horizon perturbations. S tebbins, Frieman, Hui, and Dodelson worked on phenomenological cosmology, extracting cosmological constraints from surveys such as the Sloan Digital Sky Survey. They also worked on theoretical topics such as weak lensing, reionization, and dark energy. This work has proved important to a number of experimental groups [including those at Fermilab] planning future observations. In general, the work of the Theoretical Astrophysics Group has served as a catalyst for experimental projects at Fennilab. An example of this is the Joint Dark Energy Mission. Fennilab is now a member of SNAP, and much of the work done here is by people formerly working on the accelerator. We have created an environment where many of these people made transition from physics to astronomy. We also worked on many other topics related to NASA s focus: cosmic rays, dark matter, the Sunyaev-Zel dovich effect, the galaxy distribution in the universe, and the Lyman alpha forest. The group organized and hosted a number of conferences and workshop over the years covered by the grant. Among them were:

Description: We study the rotational distortions of the vacuum dipole magnetic field in the context of geometrical models of the radio emission from pulsars. We find that at low altitudes the rotation deflects the local direction of the magnetic field by at most an angle of the order of r(sup 2 sub n), where r(sub n) = r/R(sub lc), r is the radial distance and R(sub lc) is the light cylinder radius. To the lowest (i.e. second) order in r(sub n) this distortion is symmetrical with respect to the plane containing the dipole axis and the rotation axis ((Omega, mu) plane). The lowest order distortion which is asymmetrical with respect to the (Omega, mu) plane is third order in r(sub n). These results confirm the common assumption that the rotational sweepback has negligible effect on the position angle (PA) curve. We show, however, that the influence of the sweep back on the outer boundary of the open field line region (open volume) is a much larger effect, of the order of r(sup 1/2 sub n). The open volume is shifted backwards with respect to the rotation direction by an angle delta(sub o nu) approx. 0.2 sin alpha r(sup 1/2 sub n) where alpha is the dipole inclination with respect to the rotation axis. The associated phase shift of the pulse profile Delta phi(sub o nu) approx. 0.2 r(sup 1/2 sub n) can easily exceed the shift due to combined effects of aberration and propagation time delays (approx. 2r(sub n)). This strongly affects the misalignment of the center of the PA curve and the center of the pulse profile, thereby modifying the delay radius relation. Contrary to intuition, the effect of sweepback dominates over other effects when emission occurs at low altitudes. For r(sub n) 〈 or approx. 3 x 10(exp -3) the shift becomes negative, i.e. the center of the position angle curve precedes the profile center. With the sweepback effect included, the modified delay-radius relation predicts larger emission radii and is in much better agreement with the other methods of determining r(sub n).

Description: Our group has designed a heterodyne submillimeter receiver that offers a very wide IF bandwidth of 12 GHz, while still maintaining a low noise temperature. The 180-300 GHz double-sideband design uses a single SI5 device excited by a full bandwidth, fixed-tuned waveguide probe on a silicon substrate. The IF output frequency (limited by the MMIC low noise IF preamplifier) is 6-18 GHz. providing an instantaneous RF bandwidth of 24 GHz (double-sideband). Intensive simulations predict that the junction will achieve a conversion loss better than 1-2 dB and a mixer noise temperature of less than 20 K across the band (twice the quantum limit). The single sideband receiver noise temperature goal is 70 K. The wide instantaneous bandwidth and low noise will result in an instrument capable of a variety of important astrophysical and environmental observations beyond the capabilities of current instruments. Lab testing of the receiver will begin this summer, and first light on the CSO should be in the Spring of 2003. At the CSO, we plan to use receiver with WASP2, a wideband spectrometer, to search for spectral lines from SCUBA sources. This approach should allow us to rapidly develop a catalog of redshifts for these objects.

Description: This is my final report on the NASA ATP grant on nucleosynthesis, rotation and magnetism in accreting neutron stars (NAG5-8658). In my last two reports, I summarized the science that I have accomplished, which covered a large range of topics. For this report, I want to point out the graduate students that were partially supported on this grant and where they are now. Andrew Cumming is an Assistant Professor of Physics at McGill University, Greg Ushomirsky is a researcher at MIT s Lincoln Laboratories, Dean Townsley is a postdoctoral researcher at Univ. of Chicago, Chris Deloye is a postdoctoral researcher at Northwestern University. The other two students, Phil Chang and Tony Piro, are still at UCSB and will be completing their PhD s in Summer 05 and Summer 06.

Description: We present results for r-modes of relativistic nonbarotropic stars. We show that the main differential equation, which is formally singular at lowest order in the slow-rotation expansion, can be regularized if one considers the initial value problem rather than the normal mode problem. However, a more physically motivated way to regularize the problem is to include higher order terms. This allows us to develop a practical approach for solving the problem and we provide results that support earlier conclusions obtained for uniform density stars. In particular, we show that there will exist a single r-mode for each permissible combination of 1 and m. We discuss these results and provide some caveats regarding their usefulness for estimates of gravitational-radiation reaction timescales. The close connection between the seemingly singular relativistic r-mode problem and issues arising because of the presence of co-rotation points in differentially rotating stars is also clarified.

Description: We present a detailed study of the effects of mesh refinement boundaries on the convergence and stability of simulations of black hole spacetimes. We find no technical problems. In our applications of this technique to the evolution of puncture initial data, we demonstrate that it is possible to simulaneously maintain second order convergence near the puncture and extend the outer boundary beyond 100M, thereby approaching the asymptotically flat region in which boundary condition problems are less difficult.

Description: We report results of a serendipitous hard X-ray (3-20 keV), nearly all-sky (absolute value of b greater than l0 deg.) survey based on RXTE/PCA observations performed during satellite reorientations in 1996-2002. The survey is 80% (90%) complete to a 4(sigma) limiting flux of approx. = 1.8 (2.5) x 10(exp -l1) erg/s sq cm in the 3-20 keV band. The achieved sensitivity in the 3-8 keV and 8-20 keV subbands is similar to and an order of magnitude higher than that of the previously record HEAO-1 A1 and HEAO-1 A4 all-sky surveys, respectively. A combined 7 x 10(exp 3) sq. deg area of the sky is sampled to flux levels below l0(exp -11) erg/ s sq cm (3-20 keV). In total 294 sources are detected and localized to better than 1 deg. 236 (80%) of these can be confidently associated with a known astrophysical object; another 22 likely result from the superposition of 2 or 3 closely located known sources. 35 detected sources remain unidentified, although for 12 of these we report a likely soft X-ray counterpart from the ROSAT all-sky survey bright source catalog. Of the reliably identified sources, 63 have local origin (Milky Way, LMC or SMC), 64 are clusters of galaxies and 100 are active galactic nuclei (AGN). The fact that the unidentified X-ray sources have hard spectra suggests that the majority of them are AGN, including highly obscured ones (N(sub H) greater than l0(exp 23)/sq cm). For the first time we present a log N-log S diagram for extragalactic sources above 4 x l0(exp -12) erg/ s sq cm at 8-20 keV. Key words. cosmo1ogy:observations - diffuse radiation - X-rays general

Description: The aim of this project is to examine the O VI emission at three positions around the X- ray bright ring of the remnant in order to investigate the relation between the O VI emission, the X-ray O VII and O VIII emission, and the optical [OIII} emission, and how these vary around the rim of the remnant. All three pointings and the background pointing have now been observed; the archive notification for the most recent dataset was Oct 30, 2003. After reprocessing and screening, the net exposure time for the SE exposure is only 54 percent of the approved time (15 kilosec). for the SE exposure, the available statistics are not good enough for analysis. A request for reobservation to make up for the lost time in the SE pointing has been approved. Broad O VI 1032 and O VI 1038 emission is detected with velocity width of at least 800 km/s, and possibly exceeding 1000 km/s. The Flanagan et al. analysis of the Chandra grating data show bulk velocities in the X-ray gas of order +/- 1000 km/s. In the region of the FUSE E0lO2-SE pointing, the Chandra data indicate both blue-shifted and red-shifted emission. Analysis of the velocity structure of the O VI emission will provide additional constraints on the kinematics of the gas: is it emission from a tilted expanding barrel, or a more symmetric expansion? The O VI fluxes are also needed to assess whether the O VI is radiation from recombining O VII or instead fiom cooler gas ionizing toward O VII. The emission is faint, however, which complicates the analysis. Because the lines are so broad, absorption by intervening H2, CII, and foreground OVI must be considered. A number of stars in the SMC have been observed which provide information on foreground OVI absorption. The initial analyses have concentrated on the Li1F channel since the guidance is based on that channel. The Li2F channel is being examined exposure by exposure to see if any of the data can be used to improve the signal to noise in the Li1F data.

Description: The interactions of galactic cosmic rays (GCR) with extraterrestrial bodies produce small amounts of radionuclides and stable isotopes. The production rates of many relatively short-lived radionuclides, including 2.6-year Na-22 and 312-day Mn-54, have been measured in several meteorites collected very soon after they fell. Theoretical models used to calculate production rates for comparison with the measured values rely on input data containing good cross section measurements for all relevant reactions. Most GCR particles are protons, but secondary neutrons make most cosmogenic nuclides. Calculated production rates using only cross sections for proton-induced reactions do not agree well with measurements. One possible explanation is that the contribution to the production rate from reactions initiated by secondary neutrons produced in primary GCR interactions should be included explicitly. This, however, is difficult to do because so few of the relevant cross sections for neutron-induced reactions have been measured.

Description: We present a series of simulations to demonstrate that high-fidelity velocity-delay maps of the emission-line regions in active galactic nuclei can be obtained from time-resolved spectrophotometric data sets like those that will arise from the proposed Kronos satellite. While previous reverberation-mapping experiments have established the size scale R of the broad emission-line regions from the mean time delay tau = R/c between the line and continuum variations and have provided strong evidence for supermassive black holes, the detailed structure and kinematics of the broad-line region remain ambiguous and poorly constrained. Here we outline the technical improvements that will be required to successfully map broad-line regions by reverberation techniques. For typical AGN continuum light curves, characterized by power-law power spectra P (f) is proportional to f(exp -alpha) with a = -1.5 +/- 0.5, our simulations show that a small UV/optical spectrometer like Kronos will clearly distinguish between currently viable alternative kinematic models. From spectra sampled at time intervals Delta t and sustained for a total duration T(sub dur), we can reconstruct high-fidelity velocity-delay maps with velocity resolution comparable to that of the spectra, and delay resolution Delta tau approx. 2 Delta t, provided T(sub dur) exceeds the broad-line region light crossing time by at least a factor of three. Even very complicated kinematical models, such as a Keplerian flow with superimposed spiral wave pattern, are resolved in maps from our simulated Kronos datasets. Reverberation mapping with Kronos data is therefore likely deliver the first clear maps of the geometry and kinematics in the broad emission-line regions 1-100 microarcseconds from supermassive black holes.

Description: We discuss the measurement and associated uncertainties of AGN reverberation-based black-hole masses, since these provide the zero-point calibration for scaling relationships that allow black-hole mass estimates for quasars. We find that reverberation-based mass estimates appear to be accurate to within a factor of about 3.

Description: Despite being the most energetic phenomenon in the known universe, the astrophysics of gamma-ray bursts (GRBs) has still proven difficult to understand. It has only been within the past five years that the GRB distance scale has been firmly established, on the basis of a few dozen bursts with x-ray, optical, and radio afterglows. The afterglows indicate source redshifts of z=1 to z=5, total energy outputs of roughly 10(exp 52) ergs, and energy confined to the far x-ray to near gamma-ray regime of the electromagnetic spectrum. The multi-wavelength afterglow observations have thus far provided more insight on the nature of the GRB mechanism than the GRB observations; far more papers have been written about the few observed gamma-ray burst afterglows in the past few years than about the thousands of detected gamma-ray bursts. One reason the GRB central engine is still so poorly understood is that GRBs have complex, overlapping characteristics that do not appear to be produced by one homogeneous process. At least two subclasses have been found on the basis of duration, spectral hardness, and fluence (time integrated flux); Class 1 bursts are softer, longer, and brighter than Class 2 bursts (with two second durations indicating a rough division). A third GRB subclass, overlapping the other two, has been identified using statistical clustering techniques; Class 3 bursts are intermediate between Class 1 and Class 2 bursts in brightness and duration, but are softer than Class 1 bursts. We are developing a tool to aid scientists in the study of GRB properties. In the process of developing this tool, we are building a large gamma-ray burst classification database. We are also scientifically analyzing some GRB data as we develop the tool. Tool development thus proceeds in tandem with the dataset for which it is being designed. The tool invokes a modified KDD (Knowledge Discovery in Databases) process, which is described as follows.

Description: Being the most massive gravitationally bound objects in the Universe, clusters of galaxies are prime targets for studies of structure formation and evolution. Specifically the comoving space density of virialized clusters of a given mass (or X-ray temperature), but also the frequency and degree of substructure, as well as the shape of the cluster mass profile are quantities whose current values and evolution as a function of lookback time can provide important constraints on the cosmological and physical parameters of structure formation theories. The project funded by NASA grant NAG 5-10041 intended to take such studies to a new level by combining observations of a well-selected cluster sample by three state-of-the-art telescopes: HST, to accurately measure the mass distribution in the cluster core (approx. 0.5 h(sup -1)(sub 50) Mpc) via strong gravitational lensing; CFHT, to measure the large scale mass distribution out to approx. 3 Mpc via weak lensing; and XMM, to measure the gas density and temperature distribution accurately on intermediate scales 〈 1.5 Mpc. XMM plays a pivotal role in this context as the calibration of X-ray mass measurements through accurate, spatially resolved X-ray temperature measurements (particularly in the cosmologically most sensitive range of kT〉 5 keV) is central to the questions outlined above. This set of observations promised to yield the best cluster mass measurements obtained so far for a representative sample, thus allowing us to: 1) Measure the high-mass end of the local cluster mass function; 2) Test predictions of a universal cluster mass profile; 3) calibrate the mass-temperature and temperature-luminosity relations for clusters and the scatter around these relations, which is vital for studies of cluster evolution using the X-ray temperature and X-ray luminosity functions.

Description: One of the best objects for study of the structure, kinematics, and evolutionary status of a middle-aged supernova remnant (SNR) is the VELA SNR, due to its proximity, extensive optical filamentary structure, and an abundance of hot background stars for absorption line research. The VELA remnant is 7.3 degrees in diameter, based on x-ray imagery with ROSAT, with the pulsar nearly centered in the remnant. The western region of the remnant has much lower x-ray surface brightness than the remainder of the remnant and in fact escaped earlier detection with previous instrumentation. The remnant is believed to be about 11,000 years old.

Description: In the last year, we have made an extremely important breakthrough in establishing the relationship between thermonuclear burst oscillations in accreting neutron stars and the stellar spin. More broadly, we have continued t o make significant scientific progress in all four of the key focus areas identified in our original proposal: (1) the disk-magnetosphere interaction in neutron stars, (2) rapid variability in accreting neutron stars, (3) physics of accretion flows, and (4) fundamental properties of neutron stars. A list of all publications that have arising from this work since the start of our program is given.

Description: The scope of the investigation is to extract information on the properties of the bulk solar wind from the minor ion observations that are provided by instruments on board NASA space craft and theoretical model studies. Ion charge states measured in situ in interplanetary space are formed in the inner coronal regions below 5 solar radii, hence they carry information on the properties of the solar wind plasma in that region. The plasma parameters that are important in the ion forming processes are the electron density, the electron temperature and the flow speeds of the individual ion species. In addition, if the electron distribution function deviates from a Maxwellian already in the inner corona, then the enhanced tail of that distribution function, also called halo, greatly effects the ion composition. This study is carried out using solar wind models, coronal observations, and ion calculations in conjunction with the in situ observations.

Description: We use data from Helios, IMP-8, and other spacecraft (e.g. ISEE) to systematically investigate solar energetic particle (SEP) events from different longitudes and distances in the heliosphere. The purpose of the project is to assess empirically the connection between the morphology of the travelling shock and strength with observed enhancements in the flow of energized particles in shock accelerated particle (SEP) events (also often identified as "gradual" solar energetic particle events). Activities during this first year centered on the organization of the SEPs events and their correlation with solar wind observations at multiple spacecraft locations. From an identified list of more than 30 SEPs events at multiple spacecraft locations, currently four single cases for detailed study were selected and are in an advance phase of preparation for publication. Preliminary results of these four cases were presented at AGU Spring and Fall 2003 meetings, and other meetings on SEPs.

Description: Contents include the following: On the Dynamical Evolution of a Nebula and Its Effect on Dust Coagulation and the Formation of Centimeter-sized Particles. The Mineralogy and Grain Properties of the Disk Surfaces in Three Herbig Ae/Be Stars. Astrophysical Observations of Disk Evolution Around Solar Mass Stars. The Systematic Petrology of Chondrites: A Consistent Approach to Assist Classification and Interpretation. Understanding Our Origins: Formation of Sun-like Stars in H II Region Environments. Chondrule Crystallization Experiments. Formation of SiO2-rich Chondrules by Fractional Condensation. Refractory Forsterites from Murchison (CM2) and Yamato 81020 (CO3.0) Chondrites: Cathodoluminescence, Chemical Compositions and Oxygen Isotopes. Apparent I-Xe Cooling Rates of Chondrules Compared with Silicates from the Colomera Iron Meteorite. Chondrule Formation in Planetesimal Bow Shocks: Physical Processes in the Near Vicinity of the Planetesimal. Genetic Relationships Between Chondrules, Rims and Matrix. Chondrite Fractionation was Cosmochemical; Chondrule Fractionation was Geochemical. Chondrule Formation and Accretion of Chondrite Parent Bodies: Environmental Constraints. Amoeboid Olivine Aggregates from the Semarkona LL3.0 Chondrite. The Evolution of Solids in Proto-Planetary Disks. New Nickel Vapor Pressure Measurements: Possible Implications for Nebular Condensates. Chemical, Mineralogical and Isotopic Properties of Chondrules: Clues to Their Origin. Maximal Size of Chondrules in Shock-Wave Heating Model: Stripping of Liquid Surface in Hypersonic Rarefied Gas Flow. The Nature and Origin of Interplanetary Dust: High Temperature Components. Refractory Relic Components in Chondrules from Ordinary Chondrites. Constraints on the Origin of Chondrules and CAIs from Short-lived and Long-lived Radionuclides. The Genetic Relationship Between Refractory Inclusions and Chondrules. Contemporaneous Chondrule Formation Between Ordinary and Carbonaceous Chondrites. Chondrules and Isolated Grains in the Fountain Hills Bencubbinite. Implications of Chondrule Formation in a Gas of Solar Composition. Implications of Meteoritic Cl-36 Abundance for the Origin of Short-lived Radionuclides in the Early Solar System. Size Sorting and the Chondrule Size Spectrum. Comparative Study of Refractory Inclusions from Different Groups of Chondrites. In Situ Investigation of Mg Isotope Distributions in an Allende CAI by Combined LA-ICPMS and SIMS Analyses Photochemical Speciation of Oxygen Isotopes in the Solar Nebula.

Description: Most of the thermal energy in the Galaxy and perhaps most of the baryons in the Universe are found in hot (log T approximately 5.5 - 7) gas. Hot gas is detected in the local interstellar medium, in supernova remnants (SNR), the Galactic halo, galaxy clusters and the intergalactic medium (IGM). In our own Galaxy, hot gas exists in large superbubbles up to several hundred pc in diameter that locally dominate the interstellar medium (ISM) and determine its thermal and dynamic evolution. While X-ray observations using ROSAT, Chandra and XMM have allowed us to make dramatic progress in mapping out the morphology of the hot gas and in understanding some of its spectral characteristics, there remain fundamental questions that are unanswered. Chief among these questions is the way that hot gas interacts with cooler phase gas and the effects these interactions have on hot gas energetics. The theoretical investigations we proposed in this grant aim to explore these interactions and to develop observational diagnostics that will allow us to gain much improved information on the evolution of hot gas in the disk and halo of galaxies. The first of the series of investigations that we proposed was a thorough exploration of turbulent mixing layers and cloud evaporation. We proposed to employ a multi-dimensional hydrodynamical code that includes non-equilibrium ionization (NEI), radiative cooling and thermal conduction. These models are to be applied to high velocity clouds in our galactic halo that are seen to have O VI by FUSE (Sembach et ai. 2000) and other clouds for which sufficient constraining observations exist.

Description: A recent remarkable discovery by Katharina Lodders showed that the abundances of extinct radioactivities (relative to reference nuclei) in chondrites, achondrites, and irons are proportional to the squares of their mean lives, whereas no such abundance-mean life relationship is apparent for the data of calcium-aluminum rich inclusions (CAIs) and a variety of other inclusion types. In this talk I shall interpret these results in terms of galactic and solar nebula processes. The first step in organizing the data for interpretation is to require that both the abundance of the extinct radioactivity (as measured by its decay product abundance) and that of the reference nuclide should have the same nucleosynthesis history, or else that a correction to the ratio be made to compensate for the ratio of the two different production processes. The striking feature of this diagram is that the lower edge of the data (the Lodders Line) is remarkably straight; most of the data are derived from chondrites; and its slope on the log-log diagram is two. So the extinct radioactivities are present in proportion to the square of their mean lives. All the other data lie above the Lodders line in this diagram; the sources of that data are in general somewhat larger particles than those that have contributed to the Lodders Line. A straight relationship between abundances of extinct radioactivities

Description: We present a modification of the relativistic phase shift method of determining the radio emission geometry from pulsar magnetospheres proposed by Gangadhara & Gupta (2001). Our modification provides a method of determining radio emission altitudes which does not depend on the viewing geometry and does not require polarization measurements. We suggest application of the method to the outer edges of averaged radio pulse profiles to identify magnetic field lines associated with'the edges of the pulse and, thereby, to test the geometric method based on the measurement of the pulse width at the lowest intensity level. We show that another relativistic method proposed by Blaskiewicz et al. (1991) provides upper limits for emission altitudes associated with the outer edges of pulse profiles. A comparison of these limits with the altitudes determined with the geometric method may be used to probe the importance of rotational distortions of magnetic field and refraction effects in the pulsar magnetosphere. We provide a comprehensive discussion of the assumptions used in the relativistic methods.

Description: The unexpected discovery of x-ray emission from Comet Hyakutake in March 1996 (Lisse et ai. 1996) has produced a number of questions about the physical mechanism producing the radiation. The original detection and subsequent observations (Dennerl et ai. 1997, Mumma et al. 1997, Krasnopolsky et al. 1998, Owens et al. 1998. Lisse et al. 1999) have shown that the very soft (best fit thermal bremsstrahlung model kT approx. 0.2 keV) emission is due to an interaction between the solar wind and the comet's atmosphere. Using the results fiom the 15 comets detected to date in x-rays, we report on the latest results on cometary x-ray emission, including new results from Chandra and XMM. As-observed morphologies, spectra, and light curves will be discussed. Our emphasis will be on understanding the physical mechanism producing the emission, and using this to determine the nature of the cometary coma, the structure of the solar wind in the heliosphere, and the source of the local soft x-ray background. This work has been graciously supported by grants from the NASA Planetary Astronomy and Astrophysical Data Programs.

Description: We present BVR(sub c)JHK(sub s) photometry of the optical afterglow of the gamma-ray burst GRB 021211 taken at the Magellan, MMT, and WIYN observatories between 0.7 and 50 days after the burst. We find an intrinsic spectral slope at optical and near-infrared wavelengths of 0.69 +/- 0.14 at 0.87 days. The optical decay during the first day is almost identical to that of GRB 990123 except that GRB 021211's optical afterglow was intrinsically approximately 38 times fainter and the transition from the reverse shock to the forward shock may have occurred earlier than it did for GRB 990123. We find no evidence for a jet break or the cooling break passing through optical frequencies during the first day after the burst. There is weak evidence for a break in the J-band decay between 0.89 and 1.87 days which may be due to a jet. The optical and infrared data are consistent with a relativistic fireball where the shocked electrons are in the slow cooling regime and the electron index is 2.3 +/- 0.1. The burst appears to have occurred in a homogeneous ambient medium. Our analysis suggests that the jet of GRB 021211 may have a small opening angle (1.4 deg-4.4 deg) and that the total gamma-ray energy is much less than the canonical value of 1.33 x 10(exp 51) erg. If, this is the case then most of the energy of the burst may be in another form such as a frozen magnetic field, in supernova ejecta, or in a second jet component. The host galaxy of GRB 021211 is subluminous and has a star formation rate of at least 1 solar mass/yr.

Description: A NASA spacecraft set for launch Nov. 17 is designed to help scientists determine the origin of gamma-ray bursts by pinpointing their sources in the sky and quickly focusing optical/ultraviolet and X-ray telescopes on them.

Description: The measurements by neutron detectors on Odyssey have revealed two large poleward regions with large depression of flux of epithermal and high energy neutrons. The flux of neutrons from Mars is known to be produced by the bombardment of the surface layer by galactic cosmic rays. The leakage flux of epithermal and fast neutrons has regional variation by a factor of 10 over the surface of Mars. These variations are mainly produced by variations of hydrogen content in the shallow subsurface. On Mars hydrogen is associated with water. Therefore, the Northern and Southern depressions of neutron emission could be identified as permafrost regions with very high content of water ice. These regions are much larger than the residual polar caps, and could contain the major fraction of subsurface water ice. Here we present the results of HEND neutron data deconvolution for these regions and describe the similarities and differences between them.

Description: We present first INTEGRAL and XMM-Newton ob- servations of a Seyfert galaxy, the type 2 AGN NGC 4388. Several INTEGRAL observations performed in 2003 allow us to study the spectrum in the 20 - 300 keV range. In addition two XMM-Newton observations give detailed insight into the 0.2 - 10 keV emission. Comparison with previous observations by BeppoSAX, SIGMA and CGROIOSSE show that the overall spectrum for soft X-rays up to the gamma-rays can be described by a highly absorbed and variable non-thermal component in addition to constant non-absorbed thermal emission of low abundance (2 - 5%Za), plus a constant Fe K alpha line. The hard X-ray component is well described by a simple power law with a mean photon index of 1 = 1.7. During the INTEGRAL observations the flux at 100 keV increased by a factor of 1.5. The analysis of XMM-Newton data implies that the emission below 3 keV is decoupled from the AGN and probably due to extended emission as seen in Chandra observations. The constant iron line emission is apparently also decoupled from the direct emission of the central engine and likely to be generated in the obscuring material, e.g. in the molecular torus or even further away.

Description: We completed and published two papers in the Astrophysical Journal based on research from grant. In the first paper we analyzed nine X-ray-bright Virgo early-type galaxies observed by both ASCA and ROSAT. Through spatially resolved spectroscopy, we determined the radial temperature profiles and abundances of Mg, Si, and Fe for six galaxies. The temperature profiles are consistent with isothermal temperatures outside of cooler regions at the galaxies' centers. We present new evidence for iron abundance gradients in NGC 4472 and NGC 4649 and confirm the previous results on NGC 4636. Mg and Si abundance gradients on average are flatter than those of iron and correspond to an underabundance of α-process elements at high Fe values, while at low iron the element ratios favor enrichment by Type II supernovae (SNe). We explain the observed trend using the metallicity dependence of SN Ia metal production and present constraints on the available theoretical modeling for low-metallicity inhibition of SNe Ia. In the second paper We analyzed nine X-ray-bright Virgo early-type galaxies observed by both ASCA and ROSAT. Through spatially resolved spectroscopy, we determined the radial temperature profiles and abundances of Mg, Si, and Fe for six galaxies. The temperature profiles are consistent with isothermal temperatures outside of cooler regions at the galaxies' centers. We present new evidence for iron abundance gradients in NGC 4472 and NGC 4649 and confirm the previous results on NGC 4636. Mg and Si abundance gradients on average are flatter than those of iron and correspond to an underabundance of α-process elements at high Fe values, while at low iron the element ratios favor enrichment by Type I1 supernovae (SNe). We explain the observed trend using the metallicity dependence of SN Ia metal production and present constraints on the available theoretical modeling for low-metallicity inhibition of SNe Ia.

Description: We present preliminary results of a XMM-Newton 50 ks observation of the Perseus Cluster that provides an unprecedented view of the central 0.5 Mpc region. The projected gas temperature declines smoothly by a factor of 2 from a maximum value of approx. 7 keV in the outer regions to just above 3 keV at the cluster center. Over this same range, the heavy-element abundance rises slowly from 0.4 to 0.5 solar as the radius decreases from 14 ft. to 5 ft., and then it rises to a peak of almost 0.7 solar at 1&farcm;25 before declining to 0.4 at the center. Th global east-west asymmetry of the gas temperature and surface brightness distributions, approximately aligned with the chain of bright galaxies, suggests an ongoing merger, although the modest degree of the observed asymmetry certainly excludes a major merger interpretation. The chain of galaxies probably traces the filament along which accretion started some time ago and is continuing at the present time. A cold and dense (low-entropy) cluster core like Perseus is probably well "protected" against the penetration of the gas of infalling groups and poor clusters, whereas in non-cooling core clusters such as Coma and A1367, infalling subclusters can penetrate deeply into the core region. In Perseus, gas associated with infalling groups may be stripped completely at the outskirts of the main cluster and only compression waves (shocks) may reach the central regions. We argue, and show supporting simulations, that the passage of such a wave(s) can qualitatively explain the overall horseshoe shaped appearance of the gas temperature map (the hot horseshoe surrounds the colder, low-entropy core) as well as other features of the Perseus Cluster core. These simulations also show that as compression waves traverse the cluster core, they can induce oscillatory motion of the cluster gas that can generate multiple sharp "edges" on opposite sides of the central galaxy. Gas motions induced by mergers may be a natural way to explain the high frequency of "edges" seen in clusters with cooling cores.

Description: Despite the fact that it has been well known since the earliest observations that solar energetic particle events are well associated with solar flares it is often considered that the association is not physically significant. Instead, in large events, the particles are considered to be only accelerated at a shock driven by the coronal mass ejection (CME) that is also always present. If particles are accelerated in the associated flare, it is claimed that such particles do not find access to open field lines and therefore do not escape from the low corona. However recent work has established that long lasting type III radio bursts extending to low frequencies are associated with all prompt solar particle events. Such bursts establish the presence of open field lines. Furthermore, tracing the radio bursts to the lowest frequencies, generated near the observer, shows that the radio producing electrons gain access to a region of large angular extent. It is likely that the electrons undergo cross field transport and it seems reasonable that ions do also. Such observations indicate that particle propagation in the inner heliosphere is not yet fully understood. They also imply that the contribution of flare particles in major particle events needs to be properly addressed.

Description: Our activities here at MIT have largely concentrated on four different binary X-ray pulsars: LMC X-4; 4UO352+3O/XPer; 4U0115+63; and X1908+075. We have also recently initiated a search for millisecond X-ray pulsations in RXTE archival data for several bright LMXBs using a new technique. Since this study is just getting under way, we will not report any results here. Using RXTE timing observations of LMC X-4 we have definitively measured, for the first time, the orbital decay of this high-mass X-ray binary. The e-folding decay time scale is very close to lo6 years, comparable to, but somewhat longer than, the corresponding orbital decay times for SMC X-1 and Cen X-3. We find that the orbital decay in LMC X-4 is likely driven by tidal interactions, where the asynchronism between the orbital motion and the rotation of the companion star is maintained by the evolutionary expansion of the companion. Under NASA grant NAGS7479 we carried out RXTE observations of X Per/4U0352+30 in order to track the pulse phase over a one year interval. This effort was successful in tentatively identifying a N 250-day orbital period. However, due to the fact that the observing interval was only somewhat longer than the orbital period, we asked for the observations of X Per to continue as public, or non-proprietary observations. Dr. Jean Swank kindly agreed to the continuation of the observations and they were carried out on a less frequent basis over the next year and a half. After 72 separate observations of X Per, we have the orbital period and semimajor axis firmly determined. In addition, we were able to measure the orbital eccentricity-which turns out to be remarkably small (e = 0.10) for such a wide binary orbit. This has led us establish the birth of a neutron star with a very small (or zero) natal kick.

Description: We examine a number of evolutionary scenarios for the recently discovered class of accretion-powered millisecond X-ray pulsars in ultracompact binaries, including XTE JO929-314 and XTE J1751-305. These systems have very short orbital periods of Porb = 43.6 and 42.4 minutes, respectively, and extremely small mass functions. We focus on a particular scenario that can naturally explain the present-day properties of these systems. This model invokes a donor star that was either very close to the main-sequence turnoff at the onset of mass transfer or had sufficient time to evolve during the mass-transfer phase. We have run a systematic set of binary evolution calculations with a wide range of initial conditions.

Description: The recent emergence of a new class of accretion-powered, transient, millisecond X-ray pulsars presents some difficulties for the conventional picture of accretion onto rapidly rotating magnetized neutron stars and their spin behavior during outbursts. In particular, it is not clear that the standard paradigm can accommodate the wide range in M(i.e., approx. greater than a factor of 50) over which these systems manage to accrete and the high rate of spindown that the neutron stars exhibit in at least a number of cases. When the accretion rate drops sufficiently, the X-ray pulsar is said to become a "fast rotator," and in the conventional view, this is accompanied by a transition from accretion to "propellering," in which accretion ceases and the matter is ejected from the system. On the theoretical side, we note that this scenario for the onset of propellering cannot be entirely correct because it is not energetically self-consistent. We show that, instead, the transition is likely to take place through disks that combine accretion with spindown and terminate at the corotation radius. We demonstrate the existence of such disk solutions by modifying the Shakura-Sunyaev equations with a simple magnetic torque prescription. The solutions are completely analytic and have the same dependence on M and a (the viscosity parameter) as the original Shakura-Sunyaev solutions, but the radial profiles can be considerably modified, depending on the degree of fastness. We apply these results to compute the torques expected during the outbursts of the transient millisecond pulsars and find that we can explain the large spin-down rates that are observed for quite plausible surface magnetic fields of approx. 10(exp 90 G.

Description: The flow fields of unheated, supersonic free jets from convergent and convergent-divergent nozzles operating at M = 0.99, 1.4, and 1.6 were measured using spectrally resolved Rayleigh scattering technique. The axial component of velocity and temperature data as well as density data obtained from a previous experiment are presented in a systematic way with the goal of producing a database useful for validating computational fluid dynamics codes. The Rayleigh scattering process from air molecules provides a fundamental means of measuring flow properties in a non-intrusive, particle free manner. In the spectrally resolved application, laser light scattered by the air molecules is collected and analyzed using a Fabry-Perot interferometer (FPI). The difference between the incident laser frequency and the peak of the Rayleigh spectrum provides a measure of gas velocity. The temperature is measured from the spectral broadening caused by the random thermal motion and density is measured from the total light intensity. The present point measurement technique uses a CW laser, a scanning FPI and photon counting electronics. The 1 mm long probe volume is moved from point to point to survey the flow fields. Additional arrangements were made to remove particles from the main as well as the entrained flow and to isolate FPI from the high sound and vibration levels produced by the supersonic jets. In general, velocity is measured within +/- 10 m/s accuracy and temperature within +/- 10 K accuracy.

Description: The aim of the proposal is to investigate the absorption properties of a sample of inter-mediate redshift quasars. The main goals of the project are: Measure the redshift and the column density of the X-ray absorbers; test the correlation between absorption and redshift suggested by ROSAT and ASCA data; constrain the absorber ionization status and metallicity; constrain the absorber dust content and composition through the comparison between the amount of X-ray absorption and optical dust extinction. Unanticipated low energy cut-offs where discovered in ROSAT spectra of quasars and confirmed by ASCA, BeppoSAX and Chandra. In most cases it was not possible to constrain adequately the redshift of the absorber from the X-ray data alone. Two possibilities remain open: a) absorption at the quasar redshift; and b) intervening absorption. The evidences in favour of intrinsic absorption are all indirect. Sensitive XMM observations can discriminate between these different scenarios. If the absorption is at the quasar redshift we can study whether the quasar environment evolves with the Cosmic time.

Description: We summarize the major accomplishments of our program to use high angular resolution observations at millimeter wavelengths to probe the structure of protoplanetary disks in nearby regions of star formation. The primary facilities used in this work were the Very Large Array (VLA) of the National Radio Astronomy Observatories (NRAO) located in New Mexico, and the recently upgraded Australia Telescope Compact Array (ATCA), located in Australia (to access sources in the far southern sky). We used these facilities to image thermal emission from dust particles in disks at long millimeter wavelengths, where the emission is optically thin and probes the full disk volume, including the inner regions of planet formation that remain opaque at shorter wavelengths. The best resolution obtained with the VLA is comparable to the size scales of the orbits of giant planets in our Solar System (〈 10 AU).

Description: During the past year we have developed the multispecies code to include large scale flows and shown how those solutions compare with the observations of long lived loops. The time-dependent aspects of this work are under development. We will extend the steady flow investigations to study the effect these flows have on coronal structure as observed with TRACE. Coronal observations derive from heavy-ion emission; thus, we focus on the extent to which flow may modify coronal abundances by examining the heavy-ion abundance stratification within long-lived loops. We discuss the magnitudes of the physical effects modeled and compare simulated results with TRACE observations. These results can have a profound effect on the interpretation of TRACE observations. We will also apply the potential field matching programs developed by Aad Van Ballegooijen and Harvard undergraduate, Ellen Lee, to loops modeled for this project. Having the best possible understanding of the geometric properties of the loops is important for accurate modeling.

Description: This five-year grant involves the study of young neutron stars, particularly those in supernova remnants.In the fourth year of this program, the following studies have been undertaken in support of this effort: 1.CTA 1: Following up on our ROSAT and ASCA studies of this SNR, we obtained observations with the XMM-Newton observatory to investigate the central compact source and surrounding nebula. 2. 3C 58: Based upon our earlier Chandra observations, we submitted a successful Chandra Large Project proposal for a 350 ks observation of this young neutron star and its wind nebula. 3. G347.3 - - 0.5: Our Chandra observations of portions of this SNR were aimed at studying the nonthermal X-ray emission from the remnant shell. 4. Chandra Survey for Compact Objects in Supernova Remnants: We have formed a collaboration to carry out an extensive search for young neutron stars in nearby supernova remnants. Using X-ray observations from an approved Chandra Large Project, as well as from additional approved XMM observations, we are investigating a volume-limited sample of SNRs for which there is currently no evidence of associated neutron stars.

Description: The FUSE observation has been reduced and a paper is in progress. The analysis has been slow because of the very noisy quality of the data. We anticipate submitting a paper in 2 or 3 months. The primary result is that there is no H_2 fluorescent emission. However, there seems to be no accretion disk. The upper limit to the continuum level is well below that observed by IUE, and while O VI lines are present, they are extremely faint and much narrower that the published and archival H alpha profiles. It appears that the accretion has essentially shut off entirely, leaving only very faint O VI and C III emission from the white dwarf. We are considering a magnetic model analogous to that for AE Aqr.

Description: In the completed work listed above, we present an X--ray spectral and timing analysis of 4131543-47 during its 2002 outburst based on 49 pointed observations obtained using the Rossi X-ray Timing Explorer (RXTE). The outburst reached a peak intensity of 4.2 Crab in the 2-12 keV band and declined by a factor of 32 throughout the month-long observation. A 21.9+/-0.6 mJy radio flare was detected at 1026.75 MHz two days before the X-ray maximum; the radio source was also detected late in the outburst, after the X-ray source entered the hard state. The X-ray light curve exhibits the classic shape of a rapid rise and an exponential decay. The spectrum is soft and dominated by emission from the accretion disk. The continuum is fit with a multicolor disk blackbody (kT-max = 1.04 keV) and a power-law (Gamma -2.7). Midway through the decay phase, a strong low-frequency quasi-periodic oscillation (QPO; nu = 7.3-8.1 Hz) was present for several days. The spectra feature a broad Fe K line that is asymmetric, suggesting that the line is due to relativistic broadening rather than Gomptonization. Relativistic Laor models provide much better fits to the line than non-relativistic Gaussian models, particularly near the beginning and end of our observations. The line fits yield estimates for the inner disk radius that are within 6 GM/c 2; this result and additional evidence indicates that this black hole may have a non-zero angular momentum. In current work, we are in the process of analyzing a huge set of RXTE data (-100 pointed observations) for H1743-322, which was in eruption throughout most of 2003. The source s behavior is remarkable. It displays all of the canonical X-ray states, sometimes in unusual ways. In this short report, we feature a single important result, namely, the detection of a pair of high-frequency QPOs. The bright flaring behavior is of special interest. During the flares the shape of the power spectrum and the presence of low-frequency QPOs indicate that the source is in the steep power-law (SPL or very high) state. It is in this state that the high-frequency QPOs, which are believed to emanate from the very inner accretion disk, were oberved. Specifically, we observe a 240 Hz QPO at a significance of 4.7sigma and a 162 Hz QPO at 3.7sigma. The higher frequency QPO was discovered by Homan et al. in an independent data set; they also found marginal evidence (2.5sigma)

Description: We have worked on the analysis of the Chandra and XMM observations of the nearby and distant clusters of galaxies, and on the extended sample of the distant x- ray clusters selected in the archival ROSAT PSPC data. Some of the scientific results are discussed below.

Description: We present the initial results from our study of the nonthermal continuum emission from the supernova remnants Cassiopeia A, MSB 11-61-4, a d CT-4 1. We used the INTEGRAL Core Program data to conduct this study. During the INTEGRAL mission a significant fraction of the total observing time (e.g. 35% in year one) is allocated to the Core Program and is analyzed under the auspices of the INTEGRAL Science Working Team. We report no statistically significant detections thus far but we will continue to analyze the data as more is taken. The results so far are consistent with previous measurements from e.g. RXTE and ASCA.

Description: Responding to the RXTE Cycle 7 NASA Research Announcement, we proposed to use the RXTE X-ray telescopes to intensively observe the TeV Gamma-ray Blazars Markarian 421, Markarian 501, 1ES 1959+650 and 1ES 1426+428, when their X-ray or TeV Gamma-ray fluxes would surpass preset trigger thresholds. In May and June, 2002, the Blazar 1ES 1959+650 (z=0.048) showed a series of spectacular X-ray and gamma-ray flares. Following the detection of a strong Gamma-ray flare on May 16 and 17 with the VERITAS 10 m Cherenkov Telescope, we invoked intensive RXTE observations, as well as complementary radio, optical and GeV/TeV Gamma-ray observations. From May 18 to August 14, more than 150 ksec RXTE observations were taken, yielding a unique data set with simultaneous RXTE and GeV/TeV Gamma-ray coverage.We used the financial support from the ADP program of NASA s Office for Space Science to perform a comprehensive analysis of the RXTE data. We studied in detail the temporal and spectral characteristics of the source. We collected multiwavelength data from a large number of collaborators, and performed a detailed cross-correlation analysis. Eventually, we interpreted the results in the framework of a Synchrotron-Self Compton model. The most important discovery of our research has been the detection of an orphan gamma-ray flare , not associated with an X-ray flare. The discovery showed conclusively that most models invoked to describe the non-thermal emission from blazars are overly simplistic.

Description: During the period September 1, 2002 through August 31, 2003 we have been supported in part by a small NASA 'bridge" grant NAG5-12522 to continue our theoretical investigations of the "Formation, Evolution, and Population Synthesis of Binary Systems Containing Collapsed Stars". This research includes theoretical studies of the formation and evolution of several different types of interacting binary systems containing collapsed stars. Four papers were completed under the auspices of this grant: 1. Theoretical Consideration on the Properties of Accreting Millisecond Pulsars. 2. Accretion Onto Fast X-Ray Pulsars. 3. The Effects of Binary Evolution on the Dynamics of Core Collapse and Neutron-Star.

Description: Overview. The past year for the new SEU Forum has been a highly productive one and has moved us forward on three major objectives: Involving the public in the major scientific discoveries in our SEU theme this past year, especially the exciting discoveries about dark energy, cosmology, and black holes. Implementing the recommendations of the Knappenberger Report in order to strengthen the educational coherence of our collective activities. Developing strategic partnerships with underserved communities and other key customers . Among our activities for the past year are the following: Serving the Informal Science Education Communities: Summary of Goals, Plans and Activities for the Year 9/15/04 - 9/14/05. Interactions with the Space Science Community.

Description: We modeled the long-term and small-scale evolution of molecular clouds using direct 2D and 3D magnetohydrodynamic (MHD) simulations. This work followed up on previous research by our group under auspices of the ATP in which we studied the energetics of turbulent, magnetized clouds and their internal structure on intermediate scales. Our new work focused on both global and smallscale aspects of the evolution of turbulent, magnetized clouds, and in particular studied the response of turbulent proto-cloud material to passage through the Galactic spiral potential, and the dynamical collapse of turbulent, magnetized (supercritical) clouds into fragments to initiate the formation of a stellar cluster. Technical advances under this program include developing an adaptive-mesh MHD code as a successor to ZEUS (ATHENA) in order to follow cloud fragmentation, developing a shearing-sheet MHD code which includes self-gravity and externally-imposed gravity to follow the evolution of clouds in the Galactic potential, and developing radiative transfer models to evaluate the internal ionization of clumpy clouds exposed to external photoionizing UV and CR radiation. Gammie's work at UIUC focused on the radiative transfer aspects of this program.

Description: In this research program, we obtained and analyzed an X-ray observation of the young nearby intermediate mass pre-main sequence star HD 104237 using the XMM-Newton space-based observatory. The observation was obtained on 17 Feb. 2002. This observation yielded high-quality X-ray images, spectra, and timing data which provided valuable information on the physical processes responsible for the X-ray emission. This star is a member of the group of so-called Herbig Ae/Be stars, which are young intermediate mass (approx. 2 - 4 solar masses) pre-main sequence (PMS) stars a few million years old that have not yet begun core hydrogen burning. The objective of the XMM-Newton observation was to obtain higher quality data than previously available in order to constrain possible X-ray emission mechanisms. The origin of the X-ray emission from Herbig Ae/Be stars is not yet known. These intermediate mass PMS stars lie on radiative tracks and are not expected to emit X-rays via solar-like magnetic processes, nor are their winds powerful enough to produce X-rays by radiative wind shocks as in more massive O-type stars. The emission could originate in unseen low-mass companions, or it may be intrinsic to the Herbig stars themselves if they still have primordial magnetic fields or can sustain magnetic activity via a nonsolar dynamo.

Description: We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the gamma-ray burst (GRB). Eight X-ray point sources were localized, three and five each in the original error boxes, S1 and S2, assigned for variable X-ray counteparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by Chandra on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S 1, S 1 a, is coincident with the well-determined radio location of SN 1998bw and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning approx. 1400 days, is obtained by assuming that the burst and supernova were coincident at 35.6 Mpc. When this X-ray light curve is compared with those of the X-ray af "erglows" of ordinary GRBs, X-ray flashes, and ordinary supernovae, evidence emerges for at least two classes of light curves, perhaps bounding a continuum. By 3-10 yr, all these phenomena seem to converge on a common X-ray luminosity, possibly indicative of the supernova underlying them all. This convergence strengthens the conclusion that SN 1998 bw aid GRB 980425 took place in the same object.One possible explanation for the two classes is that a (nearly) standard GRB was observed at different angles, in which case X-ray afterglows with intermediate luminosities should eventually be discovered. Finally, we comment on the contribution of GRB afterglows to the ultraluminous X-ray source population.

Description: Our goal was to use short (10 ksec) observations of selected fields in the halo of M31, to determine the size and characteristics of its X-ray population and to study the connection between globular clusters and X-ray sources. The program of observations has yet to be successfully completed. We received acceptable data from just 2 of the 5 approved fields. Nevertheless, the results were intriguing and we have submitted a paper based on this data to Nature. We find that the X-ray source density is significantly enhanced in the vicinity of one GC, providing the first observational evidence supporting the ejection hypothesis. We also find additional X-ray sources, including some which are very soft, in large enough numbers to suggest that not all could have been formed in GCs. That is, some must be descended from the same primordial halo population that produced any compact stars comprising part of the halo's dark matter. Extrapolating fiom the X-ray source population, we estimate that stellar remnants and dim old stars in the halo could comprise as much as 25% of the estimated mass (approx. 10(exp 12) Solar Mass) of the halo. These results suggest that the other approved fields should be observed soon and also provide strong motivation for the future XMM-Newton programs.

Description: The session"Origin of Planetary Systems" presented the following reports:Planetary Formation in the the Gamma Cephei System by Core:Accretion; Accretion and Heating of Particles by Supersonic Planetesimals; Planetesimal Accretion in Close Binary Systems; Collisions, Gas Flow, and the Formation of Planetesimals; Hydrodynamic Escape of a Proto-Atmosphere Just After a Giant Impact; X-Ray Flare Induced Shock Waves and Chondrule Formation in Upper Solar Nebula; and Making Water Worlds: The Role of 26Al.

Description: The Special Session: Oxygen in the Solar System, II, included the following reports:Evolution of Oxygen Isotopes in the Solar Nebula; Disequilibrium Melting of Refractory Inclusions: A Mechanism for High-Temperature Oxygen; Isotope Exchange in the Solar Nebula; Oxygen Isotopic Compositions of the Al-rich Chondrules in the CR Carbonaceous Chondrites: Evidence for a Genetic Link to Ca-Al-rich Inclusions and for Oxygen Isotope Exchange During Chondrule Melting; Nebular Formation of Fayalitic Olivine: Ineffectiveness of Dust Enrichment; Water in Terrestrial Planets: Always an Oxidant?; Oxygen Barometry of Basaltic Glasses Based on Vanadium Valence Determination Using Synchrotron MicroXANES; A New Oxygen Barometer for Solar System Basaltic Glasses Based on Vanadium Valence; The Relationship Between Clinopyroxene Fe3+ Content and Oxygen Fugacity ; and Olivine-Silicate Melt Partitioning of Iridium.

Description: The Stardust Mission session included the following reports:The Stardust:A Successful Encounter with the Remarkable Comet Wild 2; Stardust Imaging of Comet Wild 2: First Look; Preliminary Results from the Dust Flux Monitoring Instrument During the Encounter of Stardust Spacecraft with Wild-2 Comet; Streaming Clumps Ejection Model and the Heterogeneous Inner Coma of Comet Wild 2; Stardust: First Results from the Cometary and Interstellar Dust Analyzer; STARDUST Sample Collection at Wild 2 and Its Preliminary Examination; Stardust Dynamic Science at Wild 2: First Look; and Preliminary Sample Analysis Plan for the Cometary and Interstellar Samples Being Returned by the Stardust Spacecraft.

Description: The Special Session: Oxygen in the Solar System, I, included the following reports:Oxygen in the Solar System: Origins of Isotopic and Redox Complexity; The Origin of Oxygen Isotope Variations in the Early Solar System; Solar and Solar-Wind Oxygen Isotopes and the Genesis Mission; Solar 18O/17O and the Setting for Solar Birth; Oxygen Isotopes in Early Solar System Materials: A Perspective Based on Microbeam Analyses of Chondrules from CV Carbonaceous Chondrites; Insight into Primordial Solar System Oxygen Reservoirs from Returned Cometary Samples; Tracing Meteorites to Their Sources Through Asteroid Spectroscopy; Redox Conditions Among the Terrestrial Planets; Redox Complexity in Martian Meteorites: Implications for Oxygen in the Terrestrial Planets; Implications of Sulfur Isotopes for the Evolution of Atmospheric Oxygen; Oxygen in the Outer Solar System; and On the Oxidation States of the Galilean Satellites: Implications for Internal Structures.

Description: In this period, the Chandra X-ray Observatory continued to perform exceptionally well, with many scientific observations and spectacular results. The HRMA performance continues to be essentially identical to that predicted from ground calibration data. The Telescope Scientist Team has improved the mirror model to provide a more accurate description to the Chandra observers, enabling them to reduce the systematic errors and uncertainties in their data reduction. There also has been good progress in the scientific program. Using the Telescope Scientist GTO time, we carried out an extensive Chandra program to observe distant clusters of galaxies. The goals of this program were to use clusters to derive cosmological constraints and to investigate the physics and evolution of clusters. A total of 71 clusters were observed with ACIS-I; the last observations were completed in December 2003.

Description: Improved analysis of ultraviolet and optical monitoring data on the Seyfert 1 galaxy NGC 3783 provides evidence for the existence of a supermassive, (8.7 +/- 1.1) x 10(exp 6) solar mass, black hole in this galaxy. By using recalibrated spectra from the International Ultraviolet Explorer satellite and ground-based optical data, as well as refined techniques of reverberation mapping analysis, we have reduced the statistical uncertainties in the response of the emission lines to variations in the ionizing continuum. The different time lags in the emission-line responses indicate a stratification in the ionization structure of the broad-line region and are consistent with the virial relationship suggested by the analysis of similar active galaxies.

Description: This research focuses on the kinematics and evolution of the hot phase of the interstellar medium in the Galaxy. The plan was to measure the UV spectra of all hot stars observed with IUE, in order to identify and measure the main component and any high velocity components to the interstellar lines. Collection of data from higher resolution instruments on HST has been proposed for some of the interesting lines of sight. IUE spectra of 240 stars up to 8 kpc in 2 quadrants of the galactic plane have been examined to (1) estimate the total column density per kpc as a function of direction and distance, and (2) to obtain a lower limit to the number of high velocity components to the interstellar lines, thus giving an approximation of the number of conductive interfaces encountered per line of sight. By determining an approximation to the number of components per unit distance we aim to derive statistics on interfaces between hot and cold gas in the Galaxy. We find that 20% of the stars in this sample show at least one high velocity component in the C IV interstellar line. Two successful FUSE programs address this research and collected data for several of the lines of sight identified as locations of hot, expanding gas with the IUE data. One FUSE program is complete for the Vela SNR region. Data from another FUSE program to investigate the Cygnus superbubble region are being analyzed.

Description: The work completed includes the analysis of observations obtained during Cycle 6 of the Rossi X-ray Timing Explorer (RXTE). The project is part of a longer-term, continuing program to study the X-ray emission process in blazars in collaboration with Dr. Ian McHardy (U. of Southampton, UK). The goals of the program are to study the X-ray emission mechanism in blazars and the relation of the X-ray emission to changes in the relativistic jet. The program includes contemporaneous brightness and linear polarization monitoring at radio and optical wavelengths, total and polarized intensity imaging at 43 GHz with a resolution of 0.1 milliarcseconds with the VLBA, and well-sampled X-ray light curves obtained from a series of approved RXTE programs. During RXTE cycle 6, the project was awarded RXTE time to monitor three quasars: PKS 1510-089 two times per week, and 3C 273 and 3C 279 three times per week each. The X-ray data, including those from earlier cycles, were compared with radio measurements obtained 2-3 times per week in the centimeter-wave band with the Michigan radio telescope, monthly imaging observations with the VLBA at 43 GHz, and optical observations obtained at several telescopes around the world. Since it is the long-term X-ray light curve that is of greatest value, this report includes discussion of all the data collected through August 2003 rather than just the Cycle 6 data.

Description: The present report covers the first year of a grant which represents a direct continuation of NASA NAG5-4050, with the same title as before. It is dedicated as before to the discovery and characterization of new astrophysical molecules. This year, has been extremely productive, yielding many new discoveries of astronomical interest at both radio and optical wavelengths, and the publication or submission of the 15 papers listed below. Nearly all of these articles have or will soon appear in the leading refereed journals of astrophysics, chemical physics, physics, or molecular spectroscopy. One is a major invited review for Molecular Physics. One of our other invited reviews published in Spectrochimica Acta in 2001 was recently awarded the Sir Harold Thompson Memorial Award, annually given to the best paper in that journal. During the past year significant advances have been made by our group in the laboratory study of exotic silicon and carbon molecules of astronomical interest. The most exciting discoveries include the pure silicon cluster Si3, several novel silicon hydrides, and the detection of phenyl radical, C6H5, a fundamental reactive organic ring. In addition, the rotational spectra of many carbon chains terminated with Si, N, O, and other heteroatoms have also been detected for the first time. The laboratory astrophysics of the whole set is complete in the sense that the entire radio spectrum of each species has now been measured or can be calculated to very high accuracy. Nearly all of these newly found molecules are plausible candidates for the detection by radio astronomers in the interstellar gas or in circumstellar sources because they are similar in structure and composition to known astronomical species, and because most are calculated to possess large permanent dipole moments.

Description: When we last reported, our new microcalorimeter system was being prepared for delivery and permanent installation at the NIST EBIT. This occurred in June 2003 and check-out with the internal calibration source and EBIT plasma x-rays took place over the next several months during which time we modified several component parts to improve the performance. These changes included: 1) A redesign of the x-ray calibration source from a direct electron impact source to one that irradiates the microcalorimeter with fluorescent x-rays. The resulting calibration lines are free of bremsstrahlung background; 2) The microcalorimeter electronic circuit has been significantly improved to ensure long-term stability for the lengthy upcoming runs of the EBIT. Both the preamplifier feedback resistors were changed and the first stage of the preamplifier redesigned. Several photos of the new system are shown in slides 3 and 4. This microcalorimeter spectrometer only requires helium refills every three days (as opposed to every 24 hours in our earlier system) and it will hold a temperature of 65 mK for up to 48 hours (as opposed to 8 hours). Consequently, the efficiency of data acquisition will improve dramatically. The first x-ray spectra of the new calibration source made with the 4-element detector array is shown. An example of the temperature control capabilities of the ADR for a 23 hour interval is shown. The horizontal line shows the temperature stability (about +/- 3 micro kelvin). There are a few short-lived heating excursions caused by technical staff working on the EBIT machine simultaneously. During actual experimental runs these are absent. This temporal profile was interrupted to test additional components of the system; otherwise, the temperature controlling would have continued for another 24 hours.

Description: We have completed a first draft of a paper on the galaxy group ESO3060170, the hottest known fossil group. We have submitted a first draft of the paper but the final completion is delayed due to several issues mentioned by the referee that we wish to revisit and discuss in more detail. The XMM data was combined with Chandra data which allowed a rich set of projects. The paper discusses the north-south elongation which is similar to that of the central dominant galaxy as well as the galaxy distribution. We detect an X-ray 'finger' or small tail emanating from the central galaxy to the north, suggesting motion of the galaxy within the elongated gravitational potential. The overall agreement between XMM and Chandra data are excellent (although the XMM data extend to larger radii). Both data sets show a cool core centered on the dominant galaxy. Surprisingly, - the temperature maps and detailed spectra indicate that the finger of gas is NOT cool, but has the same temperature as the ambient gas. We extracted surface brightness profiles, deprojected gas density profiles, cooling time profiles, and entropy profiles. There is a sharp discontinuity in gas temperature where the surface brightness profile starts to rise rapidly at 10 kpc. This produces a decrease in the cooling time and the gas entropy within 10 kpc. The central cooling time (within 10 kpc) is less than 109 years and falls to almost half that value in the inner 5 kpc. Despite the very short cooling time, we find no evidence (even with the excellent statistics from XMM-Newton) for multi-phased gas, i.e., a cooling flow. We find two 'edges' associated with the gas distribution (common in peaked X-ray groups and galaxies). On large scales, the temperature profile is flat and disagrees with the profile predicted by Loken et al. (2003) from detailed numerical simulations. We studied the galaxy distribution within one virial radius. The galaxy concentration associated with the group is detectable only within 0.3 of the virial radius (450 kpc) given the available depth of the optical galaxy catalogs at present. We have derived total mass and gas mass distributions (from the X-ray data) and find the gas fraction approaches a constant 8% (for H0 = 70).

Description: During the past year we have been preparing our new microcalorimeter system for permanent delivery to the NIST EBIT. Unfortunately, there have been delays due to technical difficulties in the fabrication of the two-stage adiabatic demagnetization refrigerator and in the life expectancy of the thin windows used for internal thermal baffling of the infrared radiation. These problems have been solved and we are completing tests of the entire system and it will be set up at NIST during the first week of May. Several photos of the new system are shown in Figures 1A and 1B. This microcalorimeter spectrometer only requires helium refills every three days (as opposed to every 24 hours) and it will hold a temperature! of 65 mK for up to 48 hours (as opposed to 8 hours). Consequently, the efficiency of data acquisition will improve dramatically. In parallel we have published a paper that reviews our previous work (Takacs et al. 2003), especially on Fe XVII, in the context of recent measurements by other groups. This paper is included. We highlight a recent measurement of a broad band spectrum of Fe in Figure 2 that simultaneously includes L and K radiation. It is compared with the simulated spectrum of the Perseus Cluster that one could expect to obtain with a microcalorimeter in the focus of a grazing incidence telescope such as the one being designed for Constellation X. Both the charge state distributions and the relative intensity ratios of the emission lines within the particular charge state are very similar in the two spectra. This further demonstrates the importance and relevance of the laboratory measurements in predicting the components of cosmic spectra.

Description: The subject grant is for analysis of XMM data from the supernova remnant CTA1. Our investigation centered on the study of the compact source RX 50007.0+7302 that, based on our previous observations, appears to be a neutron star powering a wind nebula in the remnant interior. This compact source has also been suggested as the counterpart of the EGRET source 2EG J0008+7307. The analysis of the data from the compact source is complete. We find that the spectrum of the source is well described by a power law with the addition of a soft thermal component that may correspond to emission from hot polar cap regions or to cooling emission from a light element atmosphere over the entire star. There is evidence of extended emission on small spatial scales which may correspond to structure in the underlying synchrotron nebula. Extrapolation of the nonthermal emission component to gamma-ray energies yields a flux that is consistent with that of 2EG J0008+7307, thus strengthening the proposition that there is a gamma-ray emitting pulsar at the center of CTA 1. Our timing studies with the EPIC pn data revealed no evidence for pulsations, however; we set an upper limit of 61% on the pulsed fraction from this source. The results from this study were presented.

Description: The target galactic nucleus is ostensibly "normal," but the presence of water maser emission indicates that it may be an obscured AGN. Our primary goal has been to test this hypothesis through the detection hard X-ray emission and an X-ray spectrum characteristic of heavy absorption. The data for this program became available in May 2003. Light curves show that most of the data are good and that the background was well behaved. We have constructed images and found the X-ray emission in this target is a hard and may be extended. The presence of an AGN has been confirmed through observation with XMM. This is one of many edge-on normal galaxies that may contain AGN. The orientation of the extension will be compared to the structure of molecular gas observed at high angular resolution with VLBI. These radio observations were completed in August 2003. Correlation was completed in late 2003. The data have yet to be delivered by the ATNF. Since the start of the program, we have also obtained six radio spectra, on top of the two obtained at or close to the epoch of discovery for the maser. Early indications were that the maser emission exhibited unusually broad velocity extent and that it "jumped in velocity on time scales of weeks. This behavior has been speculated to be the signature of emission excited by jet-ISM interactions. After more extensive and more sensitive radio study, we find that (1) the emission does not in fact change velocity centroid substantially and (2) it comprises narrower, more "normal" emission features. The joint interpretation of radio and XMM data will be possible once the VLBI images and position-resolved X-ray spectra have been made, probably in summer 2004. If the radio emission arises in an accretion disk, rather than a jet, then we may predict that the VLBI images will exhibit velocity-position gradients orthogonal to the extension of X-ray emission, which presumably traces an optical narrow line or outflow region.

Description: We have completed a first draft of a paper on the galaxy group ES03060170, the hottest known fossil group. We have submitted a first draft of the paper but the final completion is delayed due to several issues mentioned by the referee that we wish to revisit and discuss in more detail. The XMM data was combined with Chandra data which allowed a rich set of projects. The paper discusses the north-south elongation which is similar to that of the central dominant galaxy as well as the galaxy distribution. We detect an X-ray "finger" or small tail emanating from the central galaxy to the north, suggesting motion of the galaxy within the elongated gravitational potential. The overall agreement between XMM and Chandra data are excellent (although the XMM data extend to larger radii). Both data sets show a cool core centered on the dominant galaxy. Surprisingly, the temperature maps and detailed spectra indicate that the finger of gas is NOT cool, but has the same temperature as the ambient gas. We extracted surface brightness profiles, deprojected gas density profiles, cooling time profiles, and entropy profiles. There is a sharp discontinuity in gas temperature where the surface brightness profile starts to rise rapidly at 10 kpc. This produces a decrease in the cooling time and the gas entropy within 10 kpc. The central cooling time (within 10 kpc) is less than l0(exp 9) years and falls to almost half that value in the inner 5 kpc. Despite the very short cooling time, we find no evidence (even with the excellent statistics from XMM-Newton) for multi-phased gas, i.e., a cooling flow. We find two "edges" associated with the gas distribution (common in peaked X-ray groups and galaxies). On large scales, the temperature profile is flat and disagrees with the profile predicted by Loken et al. (2003) from detailed numerical simulations. We studied the galaxy distribution within one virial radius. The galaxy concentration associated with the group is detectable only within 0.3 of the virial radius (450 kpc) given the available depth of the optical galaxy catalogs at present. We have derived total mass and gas mass distributions (from the X- ray data) and find the gas fraction approaches a constant 8% (for Ho = 70).

Description: This XMM project consisted of three observations of the nearby, hot galaxy cluster Triangulum Australis, one of the cluster center and two offsets. The goal was to measure the radial gas temperature profile out to large radii and derive the total gravitating mass within the radius of average mass overdensity 500. The central pointing also provides data for a detailed two-dimensional gas temperature map of this interesting cluster. We have analyzed all three observations. The derivation of the temperature map using the central pointing is complete, and the paper is soon to be submitted. During the course of this study and of the analysis of archival XMM cluster observations, it became apparent that the commonly used XMM background flare screening techniques are often not accurate enough for studies of the cluster outer regions. The information on the cluster's total masses is contained at large off-center distances, and it is precisely the temperatures for those low-brightness regions that are most affected by the detector background anomalies. In particular, our two offset observations of the Triangulum have been contaminated by the background flares ("bad cosmic weather") to a degree where they could not be used for accurate spectral analysis. This forced us to expand the scope of our project. We needed to devise a more accurate method of screening and modeling the background flares, and to evaluate the uncertainty of the XMM background modeling. To do this, we have analyzed a large number of archival EPIC blank-field and closed-cover observations. As a result, we have derived stricter background screening criteria. It also turned out that mild flares affecting EPIC-pn can be modeled with an adequate accuracy. Such modeling has been used to derive our Triangulum temperature map. The results of our XMM background analysis, including the modeling recipes, are presented in a paper which is in final preparation and will be submitted soon. It will be useful not only for our future analysis but for other XMM cluster observations as well.

Description: This grant funds a research program to use infrared extinction measurements to probe the detailed structure of dark molecular clouds and investigate the physical conditions which give rise to star and planet formation. The goals of the this program are to: 1) acquire deep infrared and molecular-line observations of a carefully selected sample of nearby dark clouds, 2) reduce and analyze the data obtained in order to produce detailed extinction maps of the clouds, 3) prepare results, where appropriate, for publication. A description of how these goals were met are included.

Description: In L-dwarfs and T-dwarfs the resonance lines of sodium and potassium are so profoundly pressure-broadened that their wings extend several hundred nanometers from line center. With accurate knowledge of the line profiles as a function of temperature and pressure: such lines can prove to be valuable diagnostics of the atmospheres of such objects. We have initiated a joint program of theoretical and experimental research to study the line-broadening of alkali atom resonance lines due to collisions with species such as helium and molecular hydrogen. Although potassium and sodium are the alkali species of most interest in the atmospheres of cool brown dwarfs and extrasolar giant planets, some of our theoretical focus this year has involved the calculation of pressure-broadening of lithium resonance lines by He, as a test of a newly developed suite of computer codes. In addition, theoretical calculations have been carried out to determine the leading long range van der Waals coefficients for the interactions of ground and excited alkali metal atoms with helium atoms, to within a probable error of 2%. Such data is important in determining the behavior of the resonance line profiles in the far wings. Important progress has been made on the experimental aspects of the program since the arrival of a postdoctoral fellow in September. A new absorption cell has been designed, which incorporates a number of technical improvements over the previous cell, including a larger cell diameter to enhance the signal, and fittings which allow for easier cleaning, thereby significantly reducing the instrument down-time.

Description: The main focus was the study of star-forming regions through high spectral- and spatial resolution observations of mm-wavelength lines, and through models of the observations. The main accomplishments were a) demonstration that more than 15 starless cores show substantial evidence of extended inward motion at about half the sound speed; b) observations of infall asymmetry in several cores, in lines of N2H(+) and DCO(+), low- depletion tracers of the "inner core"; c) observation of "infall asymmetry" of spectral lines over approx. 0.5 pc in the NGC1333 cluster-forming region; d) observations indicating that cores are nearly at rest with respect to their envelopes; and e) development of analytic, power-series solutions to the equations of motions for condensing 1-D systems (layers, cylinders and spheres).

Description: Here we discuss the Pioneer anomaly and present the next steps towards the understanding of its origin.

Description: This theory grant was awarded to study the curious nature, origin and evolution of hot gas in elliptical galaxies and their surrounding groups. Understanding the properties of this X-ray emitting gas has profound implications over the broad landscape of modern astrophysics: cosmology, galaxy formation, star formation, cosmic metal enrichment, galactic structure and dynamics, and the physics of hot gases containing dust and magnetic fields. One of our principal specific objectives was to interpret the marvelous new observations from the XMM and Chandru satellite X-ray telescopes.

Description: This report covers work carried out at the Smithsonian Astrophysical Observatory during the period 1 December 2003 to 30 November 2004 to support efforts to prepare the Kepler Input Catalog. The Catalog will be used to select the targets observed for planetary transits by Kepler.

Description: During this grant period, the physics of gamma-ray bursts was investigated. A number of new results have emerged. The importance of pair formation in high compactness burst spectra may help explain x-ray flashes; a universal jet shape is a likely explanation for the distribution of jet break times; gravitational waves may be copiously produced both in short bursts from compact mergers and in long bursts arising from collapsars; x-ray iron lines are likely to be due to interaction with the stellar atmosphere of the progenitor; prompt optical flashes from reverse shocks will give diagnostics on the Lorentz factor and the environment; GeV and TeV emission from bursts may be expected in the external shock; etc. The group working with the PI included postdocs Dr. Bing Zhang (now assistant professor at University of Nevada); Dr. Shiho Kobayashi; graduate student Lijun Gou; collaborators Drs. Tim Kallman and Martin Rees. Meszaros shared with Rees and Dr. Bohan Paczynsky the AAS Rossi Prize in 2000 for their work on the theory of gamma ray bursts. The refereed publications and conference proceedings resulting from this research are summarized below. The PI gave a number of invited talks at major conferences, also listed.

Description: We have made substantial progress in our program to understand the ionization of the Local Cloud. Our steady state photoionization models are remarkably successful in accounting for the ionization observed both by line of sight (column density) observations towards nearby stars and in situ measurements. By lifting some of the simplifying assumptions in our earlier modeling we are making further progress in this understanding. In future work on this project we will continue our exploration of photoionization models, focusing on making sense of the differences between the LIC and the BC.

Description: Citations and abstracts are provided for three publications supported by this grant.

Description: The on-orbit performance of the Chandra X-ray Observatory over its first five years of operation is reviewed. The Observatory is running smoothly and the scientific return continues to be outstanding.

Description: We present the long-term X-ray light curves and detailed spectral and timing analyses of XTE J1908+094 using the Rossi X-Ray Timing Explorer Proportional Counter Array observations covering two outbursts in 2002 and early 2003. At the onset of the first outburst, the source was found in a spectrally low/hard state lasting for approx.40 days, followed by a 3 day long transition to the high/soft state. The source flux (in 2- 10 keV) reached approx.100 mcrab on 2002 April 6, then decayed rapidly. In power spectra, we detect strong band-limited noise and varying low- frequency quasi-periodic oscillations that evolved from approx.0.5 to approx.5 Hz during the initial low/hard state of the source. We find that the second outburst closely resembled the spectral evolution of the first. The X-ray transient s overall outburst characteristics led us to classify XTE J1908+094 as a black hole candidate. Here we also derive precise X-ray position of the source using Chandra observations that were performed during the decay phase of the first outburst and following the second outburst.

Description: This program used the Far Ultraviolet Spectroscopic Explorer (FUSE) to measure the surface abundances of three helium-rich hot subdwarfs in the Galactic field, to test our hypothesis that such stars form via a late helium flash while descending the white dwarf cooling curve. If these stars form via a late He flash, the stellar envelope should be mixed with the interior, enhancing the surface carbon abundance enormously while depleting the surface hydrogen. The first observation of PG1544+488 was incomplete because detector 2 was turned off, but the repeat observations were fully successful. The observations of LB1766 had no problems and were fully successful. The observations of PG1127+019 became infeasible when FUSE lost some of its pointing control, so it was replaced with JL87, which was also observed successfully.

Description: The goal of this proposal is to detect or place sensitive upper limits on pulsations from the quiescent neutron-star X-ray binary Cen X-4. The observations were obtained on 1 March 2003. Unfortunately, due to problems with the standard pipeline processing of this observation, the data were not delivered until more than one year later, on 21 April 2004. No pulsations have yet been found in the data set. However, new optical spectra have been obtained which will lead to better constraints on the orbital period of Cen X-4 and which may enable detection of pulsations.

Description: 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.

Description: Spectral modeling and variability studies are reported in this paper, including fitting with relativistically blurred ionized-disk models designed to reproduce the radiation spectrum from the inner accretion disk around a black hole. PG 1404+226 is highly variable in flux and has a steep X-ray spectrum that varies over the course of the observation. Its strong observed soft X-ray excess can be well fitted with an ionized-disk model. The parameters derived from this model indicate that the central black hole is rotating, since the inner radius of the accretion disk lies within the last stable orbit for a non-rotating black hole. The harder spectrum of PG 1244426 also shows a soft X-ray excess, but it is better represented by a more conventional two-component phenomenonological model of a power-law and thermal emission.

Description: The broad-band X-ray continuum of NAB 0205424 is well constrained due to the excellent photon statistics obtained (about 97,700 counts), and its impressive soft X-ray excess is clearly apparent. The hard X-ray power law has become notably steeper than when NAB 0205424 was observed with ASCA, attesting to the presence of significant X-ray spectral variability. A strong and broad emission feature is detected from about 5 to 6.4 keV, and we have modeled this as a relativistic line emitted close to the black hole from a narrow annulus of the accretion disk. Furthermore, a strong X-ray flare is detected with a hard X-ray spectrum; this flare may be responsible for illuminating the inner line-emitting part of the accretion disk. The combined observational results can be broadly interpreted in terms of the "thundercloud model proposed by Merloni & Fabian (2001).

Description: The Rossi X-ray Timing Explorer PCA scans of the Galactic bulge (galactic longitude plus or minus 11 degrees) have detected 8 recent transients which have peak intensities of 10 to 400 mCrab. Some of the transient events have a fast rise and slow decay typical of accretion disk instabilities. It is common for these decays to be oscillatory, rather than steady, as if there are waves within the disk. There are also outbursts with symmetric light curves. In particular, the source in Terzan 2 which had a very long (decade) doubling of intensity peaking near the beginning of 1997, in 2004 has had two 30 day brightenings by a factor of 5 only 100 days apart. During each of these a burst was observed in snapshot observations near the peak. The source SLX 1735-269, also a burster, though not in our observations, has had irregularly repeated occurrences of fast swings between close to zero and 2-4 times normal. Some examples, such as the increase, drop, and slow recovery of GS 1826-238 suggest a change in the accretion disk such as emptying and refilling or a peculiar alignment. Follow up observations have provided deeper information about these transient sources and possible explanations for their behavior will be addressed.

Description: We have modeled the limb-brightened cavity seen in the new Spitzer Space Telescope IR images of the southwest lobe of HH 46/47 as the bow shock driven by an outflow from a young, low-mass star. We present models in which the outflow is a perfectly collimated, straight jet, models in which the jet precesses, and finally a model in which the outflow takes the form of a latitude-dependent wind. We study cases in which the outflow moves into a constant-density cloud and into a stratified cloud. We find that the best agreement with the observed cavity is obtained for the precessing jet in a stratified cloud. However, the straight jet (traveling in a stratified cloud) also gives cavity shapes close to the observed one. The latitude-dependent wind model that we have computed gives cavity shapes that are substantially wider than the observed cavity. We therefore conclude that the cavity seen in the Spitzer observations of the southwest lobe of the HH 46/47 outflow do not seem to imply the presence of a latitude-dependent wind, as it can be modeled successfully with a perfectly collimated jet model.

Description: We present the first rest-frame of approximately 4 microns detection of a Lyman break galaxy. The data were obtained using the 16 microns imaging capability of the Spitzer Infrared Spectrograph. The target object, J134026.44+634433.2, is an extremely luminous Lyman break galaxy at z=2.79, first identified in Sloan Digital Sky Survey (SDSS) spectra (as reported by Bentz et al.). The source is strongly detected with a flux of 0.94 +/- 0.02 mJy. Combining Spitzer and SDSS photometry with supporting ground-based J- and K-band data, we show that the spectral energy distribution is consistent with an actively star-forming galaxy. We also detect other objects in the Spitzer field of view, including a very red mid-infrared source. We find no evidence of a strong lens among the mid-infrared sources.

Description: We investigate extremely red objects (EROs) using near- and mid-infrared observations in five passbands (3.6 to 24 microns) obtained from the Spitzer Space Telescope, and deep ground-based R and K imaging. The great sensitivity of the Infrared Array Camera (IRAC) camera allows us to detect 64 EROs (a surface density of 2.90 +/- 0.36 arcmin(exp -2); [3.6](sub AB) is less than 23.7) in only 12 minutes of IRAC exposure time, by means of an R - [3.6] color cut (analogous to the traditional red R - K cut). A pure infrared K - [3.6] red cut detects a somewhat different population and may be more effective at selecting z greater than 1.3 EROs. We find approximately 17% of all galaxies detected by IRAC at 3.6 or 4.5 microns to be EROs. These percentages rise to about 40% at 5.8 microns, and about 60% at 8.0 microns. We utilize the spectral bump at 1.6 microns to divide the EROs into broad redshift slices using only near-infrared colors (2.2/3.6/4.5 microns). We conclude that two-thirds of all EROs lie at redshift z greater than 1.3. Detections at 24 microns imply that at least 11% of 0.6 less than z and less than 1.3 EROs and at least 22% of z greater than 1.3 EROs are dusty star-forming galaxies.

Description: We summarize the Differential Phase (DP) technique as well as the planned implementation at the Keck Interferometer. Multicolor phase measurements are potentially a powerful astrophysical probe - and can allow ground-based direct detection of extrasolar planets. Better than 0.1 mrad phase measurements in the infrared can allow the Keck Interferometer to detect radiation from the so-called hot-Jupiter or 'Roaster' class of planets. At JPL, we are presently developing and testing instrumentation that will enable these extremely sensitive measurements. First on-sky observations are expected to start in mid-2004. In this article we describe DP and other related techniques, provide an outline of the instrument and present results from preliminary laboratory experiments.

Description: Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Description: We employed density functional theory (DFT) calculations to model the photoionization behavior of benzene and small polycyclic aromatic hydrocarbons when they are embedded in a matrix of water ice in order to investigate issues raised by recent experimental work by Gudipati and Allamandola. The ionization energies of benzene, naphthalene, anthracene, and pyrene were found to be lowered by 1.5-2.1 eV in water ice. Low-lying vertical electronic excitation energies were computed with time-dependent DFT for both neutral and ionized species and are found in both cases to be remarkably unaffected by the ice matrix. Chemical behavior in ultraviolet-photoprocessed ices is also discussed, with a focus on electron recombination and pathways leading to phenol and analogous products.

Description: Two recent experimental studies have demonstrated that amino acids or amino acid precursors are generated when astrophysical ice analogs are subjected to ultraviolet (UV) irradiation at cryogenic temperatures. Understanding the complete phenomenology of photoprocessing is critical to elucidating chemical reaction mechanisms that can function within an ice matrix under very cold conditions. Pushing beyond the much better characterized study of photolytic dissociation of chemical bonds through electronic excitation, this work explored the ability of UV radiation present in the interstellar medium to ionize small molecules embedded in ices. Quantum chemical calculations, including bulk solvation effects, were used to study the ionization of hydrogen (H2), water, and methanol (CH3OH) bound in small clusters of water. Ionization potentials were found to be much smaller in the condensed phase than in the gas phase; even a small cluster can account for large changes in the ionization potentials in ice, as well as the known formation of an OH--H3O+ pair in the case of H2O photoionization. To gauge the impact of photoionization on subsequent grain chemistry, the reaction between OH and CO in the presence of H3O+ was studied and compared with the potential energy surface without hydronium present, which is relevant to chemistry following photolysis. The differences indicate that the reaction is somewhat more likely to proceed to products (H + CO2) in the case of photoionization.

Description: This paper discusses the Laser Astrometric Test of Relativity (LATOR) mission. LATOR is a Michelson-Morley-type experiment designed to test the pure tensor metric nature of gravitation the fundamental postulate of Einstein's theory of general relativity. With its focus on gravity's action on light propagation it complements other tests which rely on the gravitational dynamics of bodies.

Description: Conclusions: SIM is on track for a launch in 2010. Project 'Phase C' start in October 2005. SIM is a flexibly-pointed instrument Astrometric capability (wide angle) a) 4 pas parallax precision on targets down to V = 20 = Stellar and Galactic astrophysics; and b) 1 pas single measurement accuracy. Sensitive to terrestrial-mass planets around the nearest stars New proposal opportunities are: a) AO-2 for new (large) Key projects - 2005; and b) General Observer program, for smaller proposals - 2006.