Wind is the process that connects Mars' climate system. Measurements of Mars atmospheric winds from orbit would dramatically advance our understanding of Mars and help prepare for human exploration of the Red Planet. Multiple instrument candidates are in development and will be ready for flight in the next decade. We urge the Decadal Survey to make these measurements a priority for 2023-2032.

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength.

The Atmospheric Chemistry Suite (ACS) instrument onboard the ExoMars Trace Gas Orbiter (TGO) ESA-Roscosmos mission began science operations in March 2018. ACS Mid InfraRed (MIR) channel notably provides solar occultation observations of the martian atmosphere in the 2.3 - 4.2 $\mu$m spectral range. Here we use these observations to characterize water ice clouds before and during the MY 34 Global Dust Storm (GDS). We developed a method to detect water ice clouds with mean particle size $\leq$ 2 $\mu$m, and applied it to observations gathered between $L_s=165^\circ$ and $L_s=243^\circ$. We observe a shift in water ice clouds maximum altitudes from about 60 km before the GDS to above 90 km during the storm. These very high altitude, small-sized ($r_\mathrm{eff} \leq 0.3$ $\mu$m) water ice clouds are more frequent during MY34 compared to non-GDS years at the same season. Particle size frequently decreases with altitude, both locally within a given profile and globally in the whole dataset. We observe that the maximum altitude at which a given size is observed can increase during the GDS by several tens of km for certain sizes. We notably notice some large water ice particles ($r_\mathrm{eff}\geq1.5$ $\mu$m) at surprisingly high altitudes during the GDS (50 - 70 km). These results suggest that GDS can significantly impact the formation and properties of high altitude water ice clouds as compared to the usual perihelion dust activity.

Aims: Spectral and temporal analysis of the NuSTAR observation Galactic Be-XRB Swift J1845.7-0037. during its recent outburst. Methods: For the spectral analysis we use both phenomenological and physics-based models. We employ an often used empirical model to identify the main characteristics of the spectral shape in relation to nominal spectral characteristics of X-ray pulsars. Additionally, we used the latest version of Bulk \& Thermal comptonization model (BW), to assess the validity of the spectral components required by the empirical model and to investigate the origin of the hard X-ray emission. We also analyzed the source light-curve, studying the pulse shape at different energy ranges and tracking the spectral evolution with pulse phase by using the model independent hardness ratio (HR). Results: We find that while both the empirical and physical (BW) spectral models can produce good spectral fits, the BW model returns physically plausible best-fit values for the source parameters and does not require any additional spectral components to the non-thermal, accretion column emission. The BW model also yielded an estimation of the neutron star magnetic field placing it in the 10^12G range. Conclusions: Our results, show that the spectral and temporal characteristics of the source emission are consistent with the scattering processes expected for radiation dominated shocks within the accretion column of highly magnetized accreting neutron stars. We further indicate that physically-derived spectral models such as BW, can be used to tentatively infer fundamental source parameters, in the absence of more direct observational signatures.

We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.

We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.

The Neutron Star Interior Composition Explorer (NICER) has observed seven thermonuclear X-ray bursts from the Low Mass X-ray Binary (LMXB) neutron star 4U 1728-34 from the start of the mission's operations until February of 2019. Three of these bursts show oscillations in their decaying tail with frequencies that are within 1 Hz of the previously detected burst oscillations from this source. Two of these burst oscillations have unusual properties: They have large fractional rms amplitudes of $ 48 \pm 9 \%$ and $ 46 \pm 9 \%$, and they are detected only at photon energies above 6 keV. By contrast, the third detected burst oscillation is compatible with previous observations of this source, with a fractional rms amplitude of $7.7 \pm 1.5\%$ rms in the 0.3 to 6.2 keV energy band. We discuss the implications of these large-amplitude burst oscillations, finding they are difficult to explain with the current theoretical models for X-ray burst tail oscillations.

Thermally emitting X-ray isolated neutron stars represent excellent targets for testing cooling surface emission and atmosphere models, which are used to infer physical parameters of the neutron star. Among the seven known members of this class, RX J1605.3+3249 is the only one that still lacks confirmation of its spin period. Here we analyze NICER and XMM-Newton observations of RX J1605.3+3249, in order to address its timing and spectral behavior. Contrary to a previous tentative detection, but in agreement with the recent work by Pires et al. (2019), we find no significant pulsation with pulsed fraction higher than 1.3% (3{\sigma}) for periods above 150 ms. We also find a limit of 2.6% for periods above 2 ms, despite searches in different energy bands. The X-ray spectrum can be fit by either a double-blackbody model or by a single-temperature magnetized atmosphere model, both modified by a Gaussian absorption line at ~0.44 keV. The origin of the absorption feature as a proton cyclotron line or as an atomic transition in the neutron star atmosphere is discussed. The predictions of the best-fit X-ray models extended to IR, optical and UV bands are compared with archival data. Our results are interpreted in the framework of a fallback disk scenario.

We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231$-$1411 using NICER. The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s$^{-1}$ in the 0.35-1.5 keV band. We analyzed Fermi LAT data to update the pulse timing model covering 10 years of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and $\gamma$ rays provides insight into the geometry of the system.

We present evidence for three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). Our timing analysis reveals three sudden spin-down events of magnitudes $\Delta\nu = -23, -30,$ and $-43 \,\mu$Hz (fractional amplitudes $\Delta\nu/\nu = -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfluid components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse.

We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of 7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3-10 keV flux increased by 5% while the 30-60 keV flux decreased significantly by 35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical NPEX continuum with an added broad Gaussian emission component around the spectral peak at 20 keV. Taken together with the observed Pdot value obtained from Fermi/GBM, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase of the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field.

Studying the physical processes occurring in the region just above the magnetic poles of strongly magnetized, accreting binary neutron stars is essential to our understanding of stellar and binary system evolution. Perhaps more importantly, it provides us with a natural laboratory for studying the physics of high temperature and high density plasmas exposed to extreme radiation, gravitational, and magnetic fields. Observations over the past decade have shed new light on the manner in which plasma falling at velocities near the speed of light onto a neutron star surface is halted. Recent advances in modeling these processes have resulted in direct measurement of the magnetic fields and plasma properties. On the other hand, numerous physical processes have been identified that challenge our current picture of how the accretion process onto neutron stars works. Observation and theory are our essential tools in this regime because the extreme conditions cannot be duplicated on Earth. This white paper gives an overview of the current theory, the outstanding theoretical and observational challenges, and the importance of addressing them in contemporary astrophysics research.

Interstellar dust plays a central role in shaping the detailed structure of the interstellar medium, thus strongly influencing star formation and galaxy evolution. Dust extinction provides one of the main pillars of our understanding of interstellar dust while also often being one of the limiting factors when interpreting observations of distant objects, including resolved and unresolved galaxies. The ultraviolet (UV) and mid-infrared (MIR) wavelength regimes exhibit features of the main components of dust, carbonaceous and silicate materials, and therefore provide the most fruitful avenue for detailed extinction curve studies. Our current picture of extinction curves is strongly biased to nearby regions in the Milky Way. The small number of UV extinction curves measured in the Local Group (mainly Magellanic Clouds) clearly indicates that the range of dust properties is significantly broader than those inferred from the UV extinction characteristics of local regions of the Milky Way. Obtaining statistically significant samples of UV and MIR extinction measurements for all the dusty Local Group galaxies will provide, for the first time, a basis for understanding dust grains over a wide range of environments. Obtaining such observations requires sensitive medium-band UV, blue-optical, and mid-IR imaging and followup R ~ 1000 spectroscopy of thousands of sources. Such a census will revolutionize our understanding of the dependence of dust properties on local environment providing both an empirical description of the effects of dust on observations as well as strong constraints on dust grain and evolution models.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

We have performed a comprehensive spectral and timing analysis of the first NuSTAR observation of the high-mass X-ray binary 4U 1538-522. The observation covers the X-ray eclipse of the source, plus the eclipse ingress and egress. We use the new measurement of the mid-eclipse time to update the orbital parameters of the system and find marginally-significant evolution in the orbital period, with $\dot{P}_{\rm orb}/P_{\rm orb} = \left(-0.95 \pm 0.37\right) \times 10^{-6}$ yr$^{-1}$. The cyclotron line energy is found approximately 1.2 keV higher than RXTE measurements from 1997--2003, in line with the increased energy observed by Suzaku in 2012 and strengthening the case for secular evolution of 4U 1538-522's CRSF. We additionally characterize the behavior of the iron fluorescence and emission lines and line-of-sight absorption as the source moves into and out of eclipse.

We report on the discovery by the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) of the accreting millisecond X-ray pulsar IGR J17591-2342, detecting coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation. This implies an orbital period of ~8.8 hours and a projected semi-major axis of ~1.23 lt-s. From the binary mass function, we estimate a minimum companion mass of 0.42 solar masses, obtained assuming a neutron star mass of 1.4 solar masses and an inclination angle lower than 60 degrees, as suggested by the absence of eclipses or dips in the light-curve of the source. The broad-band energy spectrum is dominated by Comptonisation of soft thermal seed photons with a temperature of ~0.7 keV by electrons heated to 21 keV. We also detect black-body-like thermal direct emission compatible with an emission region of a few kilometers and temperature compatible with the seed source of Comptonisation. A weak Gaussian line centered on the iron K-alpha; complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, measured during the outburst fading.

Accretion disks around neutron stars regularly undergo sudden strong irradiation by Type I X-ray bursts powered by unstable thermonuclear burning on the stellar surface. We investigate the impact on the disk during one of the first X-ray burst observations with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. The burst is seen from Aql X-1 during the hard spectral state. In addition to thermal emission from the neutron star, the burst spectrum exhibits an excess of soft X-ray photons below 1 keV, where NICER's sensitivity peaks. We interpret the excess as a combination of reprocessing by the strongly photoionized disk and enhancement of the pre-burst persistent flux, possibly due to Poynting Robertson drag or coronal reprocessing. This is the first such detection for a short sub-Eddington burst. As these bursts are observed frequently, NICER will be able to study how X-ray bursts affect the disk and corona for a range of accreting neutron star systems and disk states.

We present results of recent Neutron Star Interior Composition Explorer observations of the accreting millisecond X-ray pulsar IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38 minute orbital period. NICER observed the source for approximately 26 ksec over a 5.3 day span in 2017 August, and again for 14 and 11 ksec in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hz in Rossi X-ray Timing Explorer data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z^2 method finds a best-fitting circular orbit with a period of 2278.21 s (37.97 min), a projected semi-major axis of 0.00390 lt-sec, and a barycentric pulsar frequency of 163.6561105 Hz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 e-8} solar masses, presently the smallest known for a stellar binary. The minimum donor mass ranges from about 0.005 - 0.007 solar masses, for a neutron star mass from 1.2 - 2 solar masses. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175 - 0.0155 solar masses and 19 deg < i < 27.5 deg, where the lower and upper bounds correspond to 1.4 and 2 solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 +- 0.11 % (0.3 - 3.2 keV).

Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin-frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation (QPO) near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L~1E+38 erg/s, the pulse profiles transition from single-peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turn-over to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ~1E+13 G.

The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron Star Interior Composition Explorer (NICER) mission, which is scheduled to launch in 2017 and will be hosted as an externally attached payload on the International Space Station (ISS). During NICER's 18-month baseline science mission to understand ultra-dense matter through observations of neutron stars in the soft X-ray band, SEXTANT will, for the first-time, demonstrate real-time, on-board X-ray pulsar navigation. Using NICER/SEXTANT as an example, we describe the factors that determine the measurement errors on pulse times of arrival, including source and background count rates, and pulse profile shapes. We then describe properties of the SEXTANT navigation pulsar catalog and prospects for growing it once NICER launches. Finally, we describe the factors affecting the prediction of pulse arrival times in advance, including variable interstellar propagation effect and red timing noise. Together, all of these factors determine how well a particular realization of an X-ray pulsar-based navigation system will perform.

We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO2 absorptions and the known CO2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising inter-annual increase in water-ice cloud opacities.

Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 10^12 G) close to the surface of an accreting X-ray pulsar. The line profiles of CRSFs cannot be described by an analytic expression. Numerical methods such as Monte Carlo (MC) simulations of the scattering processes are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time consuming simulations possible on convenient time scales. Furthermore, it can generate responses to mono-energetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo .

We present a concept for using a polarization sensitive multispectral lidar to map the seasonal distribution and exchange of volatiles among the reservoirs of the Martian surface and atmosphere. The LIDAR instrument will be a multi-wavelength, altitude-resolved, active near-infrared (NIR, with 10 bands around 1.6 microns) instrument to measure the reflected intensity and polarization of backscattered radiation from planetary surfaces and atmospheres. The proposed instrument would be ideally suited for a mission to Mars to comprehensively investigate the nature and seasonal distributions of volatiles and aerosols. The investigation would include the abundance of atmospheric dust and condensed volatiles, surface and cloud/aerosol grain sizes and shapes, ice and dust particle microphysics and also variations in atmospheric chemistry during multiple overflight local times throughout polar night and day. Such an instrument would be ideal for mapping and detection of recently detected CO2 frost phenomena and H2O and CO2 precipitation events in the polar regions of Mars. Herein we discuss the applicability of this instrument to detect and map sublimation/deposition 'mode flips' recently discovered by Brown et al. (2016) using the CRISM passive infrared sensor on Mars Reconnaissance Orbiter.

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01<B/B_{crit}<=0.12, where B_{crit}=4.413x10^{13} G and 3keV<=kT<15keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo.

We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of the Her X-1 35-day accretion disk precession period. This model allows us to estimate the accretion rate, the Comptonizing temperature of the radiating plasma, the radius of the magnetic polar cap, and the average scattering opacity parameters in the accretion column. This is in contrast to previous phenomenological models that characterized the shape of the X-ray spectrum but could not determine the physical parameters of the accretion flow. We describe the spectral fitting details and discuss the interpretation of the accretion flow physical parameters.

The past decade has seen significant progress on the direct detection and characterization of young, self-luminous giant planets at wide orbital separations from their host stars. Some of these planets show evidence for disequilibrium processes like transport-induced quenching in their atmospheres; photochemistry may also be important, despite the large orbital distances. These disequilibrium chemical processes can alter the expected composition, spectral behavior, thermal structure, and cooling history of the planets, and can potentially confuse determinations of bulk elemental ratios, which provide important insights into planet-formation mechanisms. Using a thermo/photochemical kinetics and transport model, we investigate the extent to which disequilibrium chemistry affects the composition and spectra of directly imaged giant exoplanets. Results for specific "young Jupiters" such as HR 8799 b and 51 Eri b are presented, as are general trends as a function of planetary effective temperature, surface gravity, incident ultraviolet flux, and strength of deep atmospheric convection. We find that quenching is very important on young Jupiters, leading to CO/CH4 and N2/NH3 ratios much greater than, and H2O mixing ratios a factor of a few less than, chemical-equilibrium predictions. Photochemistry can also be important on such planets, with CO2 and HCN being key photochemical products. Carbon dioxide becomes a major constituent when stratospheric temperatures are low and recycling of water via the H2 + OH reaction becomes kinetically stifled. Young Jupiters with effective temperatures <~ 700 K are in a particularly interesting photochemical regime that differs from both transiting hot Jupiters and our own solar-system giant planets.

We present an attempt to extract information about the comet 9P/Tempel 1 nucleus from the characteristics of the ejecta cloud produced by the impactor of the Deep Impact mission. For this purpose we use two techniques. We first study the shadow cast on the nucleus surface by the ejecta cloud and investigate how areas of different brightness are related to the varying optical thickness or albedo of the ejecta cloud. The shadow was seen during the first 2.0 seconds after the impact (afterward it became obscured by the ejecta cloud). We have found that all brightness variations in the shadow are the result of the surface inhomogeneities, indicating that during first 2.0 seconds the ejecta cloud was homogeneous within the MRI spatial resolution. Our second technique is to study the obscuration of the nucleus limb by the ejecta. This study covers the period 0.76- 68.8 seconds after impact and is based on comparison of the ejecta cloud brightness on the limb and just beyond the limb. At this stage we do see inhomogeneities in the ejecta cloud that relate to the albedo and optical thickness variations in the ejected dust. Specifically, we have found two distinct bands of low optical thickness and one band of a high optical thickness. Based on crater formation ideas we estimate the depth of excavation of the ejected material for the found inhomogeneities and, thus, define a potential layering structure for the comet nucleus, Our estimates suggest that the low-optical thickness material was excavated from a depth of 15-18 and 30-32 meters in the case the porous nucleus material and 37-46 and 87-93 meters in the case of a non-porous nucleus material, and a layer of high optical thickness originated from the depth 9-11 m for porous material or 20-23 m for non-porous material. Based on the crater diameter estimates, we expect that the real depth of the layers is between these two cases.

Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura linearization method and the Nakayama instability boundary conditions, we discover that there are regions of parameter space where disks/shocks with outflows can be stable or unstable. In regions of instability, we find that preshock deceleration is always unstable to the zeroth mode with zero frequency of oscillation, but is always stable to the fundamental mode and overtones. Furthermore, we also find that preshock acceleration is always unstable to the zeroth mode and that the fundamental mode and overtones become increasingly less stable as the shock location moves away from the horizon when the disk half-height expands above $\sim 12$ gravitational radii at the shock radius. In regions of stability, we demonstrate the zeroth mode to be stable for the velocity profiles that exhibit preshock acceleration and deceleration. Moreover, for models that are linearly unstable, our model suggests the possible existence of quasi-periodic oscillations (QPOs) with ratios 2:3 and 3:5. These ratios are believed to occur in stellar and supermassive black hole candidates, for example, in GRS 1915+105 and Sgr A*, respectively. We expect that similar QPO ratios also exist in regions of stable shocks.

We report on radio timing and multiwavelength observations of the 4.66 ms redback pulsar J1048+2339, which was discovered in an Arecibo search targeting the Fermi-LAT source 3FGLJ1048.6+2338. Two years of timing allowed us to derive precise astrometric and orbital parameters for the pulsar. PSR J1048+2339 is in a 6-hour binary, and exhibits radio eclipses over half the orbital period and rapid orbital period variations. The companion has a minimum mass of 0.3 solar masses, and we have identified a $V \sim 20$ variable optical counterpart in data from several surveys. The phasing of its $\sim 1$~mag modulation at the orbital period suggests highly efficient and asymmetric heating by the pulsar wind, which may be due to an intrabinary shock that is distorted near the companion, or to the companion's magnetic field channeling the pulsar wind to specific locations on its surface. We also present gamma-ray spectral analysis of the source and preliminary results from searches for gamma-ray pulsations using the radio ephemeris.

New low-resolution UV spectra of a sample of reddened OB stars in M31 were obtained with HST/STIS to study the wavelength dependence of interstellar extinction and the nature of the underlying dust grain populations. Extinction curves were constructed for four reddened sightlines in M31 paired with closely matching stellar atmosphere models. The new curves have a much higher S/N than previous studies. Direct measurements of N(H I) were made using the Ly$\alpha$ absorption lines enabling gas-to-dust ratios to be calculated. The sightlines have a range in galactocentric distance of 5 to 14 kpc and represent dust from regions of different metallicities and gas-to-dust ratios. The metallicities sampled range from Solar to 1.5 Solar. The measured curves show similarity to those seen in the Milky Way and the Large Magellanic Cloud. The Maximum Entropy Method was used to investigate the dust composition and size distribution for the sightlines observed in this program finding that the extinction curves can be produced with the available carbon and silicon abundances if the metallicity is super-Solar.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

In this document, we summarize the main capabilities of the James Webb Space Telescope (JWST) for performing observations of Mars. The distinctive vantage point of JWST at the Sun-Earth Lagrange point (L2) will allow sampling the full observable disk, permitting the study of short-term phenomena, diurnal processes (across the East-West axis) and latitudinal processes between the hemispheres (including seasonal effects) with excellent spatial resolutions (0.07 arcsec at 2 {\mu}m). Spectroscopic observations will be achievable in the 0.7-5 {\mu}m spectral region with NIRSpec at a maximum resolving power of 2700, and with 8000 in the 1-1.25 {\mu}m range. Imaging will be attainable with NIRCam at 4.3 {\mu}m and with two narrow filters near 2 {\mu}m, while the nightside will be accessible with several filters in the 0.5 to 2 {\mu}m. Such a powerful suite of instruments will be a major asset for the exploration and characterization of Mars. Some science cases include the mapping of the water D/H ratio, investigations of the Martian mesosphere via the characterization of the non-LTE CO$_2$ emission at 4.3 {\mu}m, studies of chemical transport via observations of the O$_2$ nightglow at 1.27 {\mu}m, high cadence mapping of the variability dust and water ice clouds, and sensitive searches for trace species and hydrated features on the Martian surface. In-flight characterization of the instruments may allow for additional science opportunities.

In a search with the Parkes radio telescope of 56 unidentified Fermi-LAT gamma-ray sources, we have detected 11 millisecond pulsars (MSPs), 10 of them discoveries, of which five were reported in Kerr et al. (2012). We did not detect radio pulsations from another six pulsars now known in these sources. We describe the completed survey, which included multiple observations of many targets done to minimize the impact of interstellar scintillation, acceleration effects in binary systems, and eclipses. We consider that 23 of the 39 remaining sources may still be viable pulsar candidates. We present timing solutions and polarimetry for five of the MSPs, and gamma-ray pulsations for PSR J1903-7051 (pulsations for five others were reported in the second Fermi-LAT catalog of gamma-ray pulsars). Two of the new MSPs are isolated and five are in >1 d circular orbits with 0.2-0.3 Msun presumed white dwarf companions. PSR J0955-6150, in a 24 d orbit with a ~0.25 Msun companion but eccentricity of 0.11, belongs to a recently identified class of eccentric MSPs. PSR J1036-8317 is in an 8 hr binary with a >0.14 Msun companion that is probably a white dwarf. PSR J1946-5403 is in a 3 hr orbit with a >0.02 Msun companion with no evidence of radio eclipses.

Context. The diffusive shock acceleration mechanism has been widely accepted as the acceleration mechanism for galactic cosmic rays. While self-consistent hybrid simulations have shown how power-law spectra are produced, detailed information on the interplay of diffusive particle motion and the turbulent electromagnetic fields responsible for repeated shock crossings are still elusive. Aims. The framework of test-particle theory is applied to investigate the effect of diffusive shock acceleration by inspecting the obtained cosmic-ray energy spectra. The resulting energy spectra can be obtained this way from the particle motion and, depending on the prescribed turbulence model, the influence of stochastic acceleration through plasma waves can be studied. Methods. A numerical Monte-Carlo simulation code is extended to include collisionless shock waves. This allows one to trace the trajectories of test particle while they are being accelerated. In addition, the diffusion coefficients can be obtained directly from the particle motion, which allows for a detailed understanding of the acceleration process. Results. The classic result of an energy spectrum with $E^{-2}$ is only reproduced for parallel shocks, while, for all other cases, the energy spectral index is reduced depending on the shock obliqueness. Qualitatively, this can be explained in terms of the diffusion coefficients in the directions that are parallel and perpendicular to the shock front.

This is an overview of pulsar accretion modeling. The physics of pulsar accretion, i.e., the process of plasma flow onto the neutron star surface, can be constrained from the spectral properties of the X-ray source. We discuss a new implementation of the physical continuum model developed by Becker and Wolff (2007, ApJ 654, 435). The model incorporates Comptonized blackbody, bremsstrahlung, and cyclotron emission. We discuss preliminary results of applying the new tool to the test cases of Suzaku data of Cen X-3 and XTE J1946+274. Cen X-3 is a persistent accreting pulsar with an O-star companion observed during a bright period. XTE J1946+274 is a transient accreting pulsar with a Be companion observed during a dim period. Both sources show spectra that are well described with an empirical Fermi Dirac cutoff power law model. We extend the spectral analysis by making the first steps towards a physical description of Cen X-3 and XTE J1946+274.

We present new Chandra and XMM-Newton observations of a sample of eight radio-quiet Gamma-ray pulsars detected by the Fermi Large Area Telescope. For all eight pulsars we identify the X-ray counterpart, based on the X-ray source localization and the best position obtained from Gamma-ray pulsar timing. For PSR J2030+4415 we found evidence for an about 10 arcsec-long pulsar wind nebula. Our new results consolidate the work from Marelli et al. 2011 and confirm that, on average, the Gamma-ray--to--X-ray flux ratios (Fgamma/Fx) of radio-quiet pulsars are higher than for the radio-loud ones. Furthermore, while the Fgamma/Fx distribution features a single peak for the radio-quiet pulsars, the distribution is more dispersed for the radio-loud ones, possibly showing two peaks. We discuss possible implications of these different distributions based on current models for pulsar X-ray emission.

We report on the interannual variability of the atmospheric ice/dust cycle in the Martian polar regions for Mars Years 28-30. We used CRISM emission phase function measurements to derive atmospheric dust optical depths and data from the MARCI instrument to derive atmospheric water ice optical depths. We have used autocorrelation and cross correlation functions in order to quantify the degree to which dust and ice are correlated throughout both polar regions during Mars Years 28-29. We find that in the south polar region, dust has the tendency to "self clear", demonstrated by negative autocorrelation around the central peak. This does not occur in the north polar region. In the south polar region, dust and ice are temporally and spatially anti correlated. In the north polar region, this relationship is reversed, however temporal correlation of northern dust and ice clouds is weak - 6 times weaker than the anticorrelation in the south polar region. Our latitudinal autocorrelation functions allow us to put average spatial sizes of event cores and halos. Dust events in the south are largest, affecting almost the entire pole, whereas dust storms are smaller in the north. Ice clouds in north are similar in latitudinal extent to those in the south (both have halos < 10{\deg}). Using cross-correlation functions of water ice and dust, we find that dust events temporally lag ice events by 35-80 degrees of solar longitude in the north and south poles, which is likely due to seasonality of dust and ice events.

XSS J12270-4859 is an X-ray binary associated with the Fermi LAT gamma-ray source 1FGL J1227.9-4852. In 2012 December, this source underwent a transition where the X-ray and optical luminosity dropped and the spectral signatures of an accretion disc disappeared. We report the discovery of a 1.69 millisecond pulsar (MSP), PSR J1227-4853, at a dispersion measure of 43.4 pc cm$^{-3}$ associated with this source, using the GMRT at 607 MHz. This demonstrates that, post-transition, the system hosts an active radio MSP. This is the third system after PSR J1023+0038 and PSR J1824-2452I showing evidence of state switching between radio MSP and low-mass X-ray binary (LMXB) states. We report timing observations of PSR J1227-4853 with the GMRT and Parkes, which give a precise determination of the rotational and orbital parameters of the system. The companion mass measurement of 0.17 to 0.46 M$_{sun}$ suggests that this is a redback system. PSR J1227-4853 is eclipsed for about 40% of its orbit at 607 MHz; with additional short-duration eclipses at all orbital phases. We also find that the pulsar is very energetic, with a spin-down luminosity of ~ 10$^{35}$ erg s$^{-1}$. We report simultaneous imaging and timing observations with the GMRT, which suggests that eclipses are caused by absorption, rather than dispersion smearing or scattering.

We have produced a multiannual climatology of airborne dust from Martian year 24 to 31 using multiple datasets of retrieved or estimated column optical depths. The datasets are based on observations of the Martian atmosphere from April 1999 to July 2013 made by different orbiting instruments: the Thermal Emission Spectrometer (TES) aboard Mars Global Surveyor, the Thermal Emission Imaging System (THEMIS) aboard Mars Odyssey, and the Mars Climate Sounder (MCS) aboard Mars Reconnaissance Orbiter (MRO). The procedure we have adopted consists of gridding the available retrievals of column dust optical depth (CDOD) from TES and THEMIS nadir observations, as well as the estimates of this quantity from MCS limb observations. Our gridding method calculates averages and uncertainties on a regularly spaced, but possibly incomplete, spatio-temporal grid, using an iterative procedure weighted in space, time, and retrieval uncertainty. In order to evaluate strengths and weaknesses of the resulting gridded maps, we validate them with independent observations of CDOD. We have statistically analyzed the irregularly gridded maps to provide an overview of the dust climatology on Mars over eight years, specifically in relation to its interseasonal and interannual variability. Finally, we have produced multiannual, regular daily maps of CDOD by spatially interpolating the irregularly gridded maps using a kriging method. These synoptic maps are used as dust scenarios in the Mars Climate Database version 5, and are useful in many modelling applications in addition to forming a basis for instrument intercomparisons. The derived dust maps for the eight available Martian years are publicly available and distributed with open access.

X-ray observations of quiescent X-ray binaries have the potential to provide insight into the structure and the composition of neutron stars. EXO 0748-676 had been actively accreting for over 24 yr before its outburst ceased in late 2008. Subsequent X-ray monitoring revealed a gradual decay of the quiescent thermal emission that can be attributed to cooling of the accretion-heated neutron star crust. In this work, we report on new Chandra and Swift observations that extend the quiescent monitoring to ~5 yr post-outburst. We find that the neutron star temperature remained at ~117 eV between 2009 and 2011, but had decreased to ~110 eV in 2013. This suggests that the crust has not fully cooled yet, which is supported by the lower temperature of ~95 eV that was measured ~4 yr prior to the accretion phase in 1980. Comparing the data to thermal evolution simulations reveals that the apparent lack of cooling between 2009 and 2011 could possibly be a signature of convection driven by phase separation of light and heavy nuclei in the outer layers of the neutron star.

We analyzed Deep Impact High Resolution Instrument (HRI) images acquired within the first seconds after collision of the Deep Impact impactor with the nucleus of comet 9P/Tempel 1. These images reveal an optically thick ejecta plume that casts a shadow on the surface of the nucleus. Using the 3D radiative transfer code HYPERION we simulated light scattering by the ejecta plume, taking into account multiple scattering of light from the ejecta, the surrounding nuclear surface and the actual observational geometry (including an updated plume orientation geometry that accounts for the latest 9P/Tempel 1 shape model). Our primary dust model parameters were the number density of particles, their size distribution and composition. We defined the composition through the density of an individual particle and the ratio of its material constituents, which we considered to be refractories, ice and voids. The results of our modeling indicate a dust/ice mass ratio for the ejecta particles of at least 1. To further constrain the parameters of the model, we checked for consistency between the ejecta mass resulting from our modeling with the ejecta mass estimated by the crater formation modeling. Constraining the particle size distribution by results of other studies of the Deep Impact ejecta, we find the number density of ejecta particles equal to ~10^4 particles/cm^3 at the base of the plume.

We present the scientific case to build a multiple-wavelength, active, near-infrared (NIR) instrument to measure the reflected intensity and polarization characteristics of backscattered radiation from planetary surfaces and atmospheres. We focus on the ability of such an instrument to enhance, perhaps revolutionize, our understanding of climate, volatiles and astrobiological potential of modern-day Mars.

Context: Accretion-powered X-ray pulsars show highly energy-dependent and complex pulse-profile morphologies. Significant deviations from the average pulse profile can appear, in particular close to the cyclotron line energies. These deviations can be described as energy-dependent phase lags, that is, as energy-dependent shifts of main features in the pulse profile. Aims: Using a numerical study we explore the effect of cyclotron resonant scattering on observable, energy-resolved pulse profiles. Methods: We generated the observable emission as a function of spin phase, using Monte Carlo simulations for cyclotron resonant scattering and a numerical ray-tracing routine accounting for general relativistic light-bending effects on the intrinsic emission from the accretion columns. Results: We find strong changes in the pulse profile coincident with the cyclotron line energies. Features in the pulse profile vary strongly with respect to the average pulse profile with the observing geometry and shift and smear out in energy additionally when assuming a non-static plasma. Conclusions: We demonstrate how phase lags at the cyclotron energies arise as a consequence of the effects of angular redistribution of X-rays by cyclotron resonance scattering in a strong magnetic field combined with relativistic effects. We also show that phase lags are strongly dependent on the accretion geometry. These intrinsic effects will in principle allow us to constrain a system's accretion geometry.

Dust aerosol plays a fundamental role in the behavior and evolution of the Martian atmosphere. The first five Mars years of Mars Exploration Rover data provide an unprecedented record of the dust load at two sites. This record is useful for characterization of the atmosphere at the sites and as ground truth for orbital observations. Atmospheric extinction optical depths have been derived from solar images after calibration and correction for time-varying dust that has accumulated on the camera windows. The record includes local, regional, and globally extensive dust storms. Comparison with contemporaneous thermal infrared data suggests significant variation in the size of the dust aerosols, with a 1 {\mu}m effective radius during northern summer and a 2 {\mu}m effective radius at the onset of a dust lifting event. The solar longitude (LS) 20-136{\deg} period is also characterized by the presence of cirriform clouds at the Opportunity site, especially near LS=50 and 115{\deg}. In addition to water ice clouds, a water ice haze may also be present, and carbon dioxide clouds may be present early in the season. Variations in dust opacity are important to the energy balance of each site, and work with seasonal variations in insolation to control dust devil frequency at the Spirit site.

We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

Massive young stellar objects (YSOs), like low-mass YSOs, are thought to be surrounded by optically thick envelopes and/or discs and are observed to have associated regions that produce polarized light at near-infrared wavelengths. These polarized regions are thought to be lower-density outflows along the polar axes of the YSO envelopes. Using the 0.2 arcsec spatial resolution of the Near-Infrared Camera and Multi-Object Spectrometer on the Hubble Space Telescope we are examining the structure of the envelopes and outflow regions of massive YSOs in star-forming regions within a few kpc of the Sun. Here we report on 2 micron polarimetry of Mon R2-IRS3, S140-IRS1, and AFGL 2591. All three sources contain YSOs with highly-polarized monopolar outflows, with Mon R2-IRS3 containing at least two YSOs in a small cluster. The central stars of all four YSOs are also polarized, with position angles perpendicular to the directions of the outflows. We infer that this polarization is due to scattering and absorption by aligned grains. We have modelled our observations of S140-IRS1 and AFGL 2591 as light scattered and absorbed both by spherical grains and by elongated grains that are aligned by magnetic fields. Models that best reproduce the observations have a substantial toroidal component to the magnetic field in the equatorial plane. Moreover, the toroidal magnetic field in the model that best fits AFGL 2591 extends a large fraction of the height of the model cavity, which is 10^5 au. We conclude that the massive YSOs in this study all show evidence of the presence of a substantial toroidal magnetic field.

We present the synthesis of CRISM EPF and MARCI data to examine the evolution of atmospheric water ice and dust opacity at both poles for MY 28-30.

We have updated our publicly available dust radiative transfer code (HOCHUNK3D) to include new emission processes and various 3-D geometries appropriate for forming stars. The 3-D geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2-D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium. We include the option for simple power-law envelope geometry, which combined with fractal clumping, and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from PAHs and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: 1) Outflow cavities may be more clumpy than infalling envelopes. 2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband SED slope; and related to this, 3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as YSOs. 4) Our hydrostatic equilibrium models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.

Using the Giant Metrewave Radio Telescope (GMRT) we performed deep observations to search for radio pulsations in the directions of unidentified Fermi Large Area Telescope (LAT) gamma-ray sources. We report the discovery of an eclipsing black-widow millisecond pulsar, PSR J1544+4937, identified with the un-cataloged gamma-ray source Fermi J1544.2+4941. This 2.16 ms pulsar is in a 2.9 hours compact circular orbit with a very low-mass companion (Mc > 0.017 Msun). At 322 MHz this pulsar is found to be eclipsing for 13% of its orbit, whereas at 607 MHz the pulsar is detected throughout the low-frequency eclipse phase. Variations in the eclipse ingress phase are observed, indicating a clumpy and variable eclipsing medium. Moreover, additional short-duration absorption events are observed around the eclipse boundaries. Using the radio timing ephemeris we were able to detect gamma-ray pulsations from this pulsar, confirming it as the source powering the gamma-ray emission.

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3--50\,keV band. In agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7keV. Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10keV, the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the continuum.

Accreting X-ray pulsars are among the best observed objects of X-ray astronomy with a rich data set of observational phenomena in the spectral and timing domain. While the general picture for these sources is well established, the detailed physics behind the observed phenomena are often subject of debate. We present recent observational, theoretical and modeling results for the structure and dynamics of the accretion column in these systems. Our results indicate the presence of different accretion regimes and possible explanations for observed variations of spectral features with luminosity.

A model for the formation and distribution of sedimentary rocks on Mars is proposed. The rate-limiting step is supply of liquid water from seasonal melting of snow or ice. The model is run for a O(10^2) mbar pure CO2 atmosphere, dusty snow, and solar luminosity reduced by 23%. For these conditions snow only melts near the equator, and only when obliquity >40 degrees, eccentricity >0.12, and perihelion occurs near equinox. These requirements for melting are satisfied by 0.01-20% of the probability distribution of Mars' past spin-orbit parameters. Total melt production is sufficient to account for aqueous alteration of the sedimentary rocks. The pattern of seasonal snowmelt is integrated over all spin-orbit parameters and compared to the observed distribution of sedimentary rocks. The global distribution of snowmelt has maxima in Valles Marineris, Meridiani Planum and Gale Crater. These correspond to maxima in the sedimentary-rock distribution. Higher pressures and especially higher temperatures lead to melting over a broader range of spin-orbit parameters. The pattern of sedimentary rocks on Mars is most consistent with a Mars paleoclimate that only rarely produced enough meltwater to precipitate aqueous cements and indurate sediment. The results suggest intermittency of snowmelt and long globally-dry intervals, unfavorable for past life on Mars. This model makes testable predictions for the Mars Science Laboratory rover at Gale Crater. Gale Crater is predicted to be a hemispheric maximum for snowmelt on Mars.

Accretion-powered X-ray pulsars exhibit significant variability of the Cyclotron Resonance Scattering Feature (CRSF) centroid energy on pulse-to-pulse timescales, and also on much longer timescales. Two types of spectral variability are observed. For sources in group 1, the CRSF energy is negatively correlated with the variable source luminosity, and for sources in group 2, the opposite behavior is observed. The physical basis for this bimodal behavior is currently not understood. We explore the hypothesis that the accretion dynamics in the group 1 sources is dominated by radiation pressure near the stellar surface, and that Coulomb interactions decelerate the gas to rest in the group 2 sources. We derive a new expression for the critical luminosity such that radiation pressure decelerates the matter to rest in the supercritical sources. The formula for the critical luminosity is evaluated for 5 sources, using the maximum value of the CRSF centroid energy to estimate the surface magnetic field strength. The results confirm that the group 1 sources are supercritical and the group 2 sources are subcritical, although the situation is less clear for those highly variable sources that cross over the critical line. We also explain the variation of the CRSF energy with luminosity as a consequence of the variation of the characteristic emission height. The sign of the height variation is opposite in the supercritical and subcritical cases, hence creating the observed bimodal behavior.

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J1803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6e35 erg/s and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J0622+3749, J1620-4927, J1746-3239, J2028+3332, J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| > 10 deg). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2e11 G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3e33 erg/s) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars.

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of {\gamma}-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its {\gamma}-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of {\gamma}-ray emission in the off-pulse phase interval at the ~ 6{\sigma} level. The flux from this emission in the energy range E \geq 100 MeV is F_100 = (1.73\pm0.40)\times10^(-8) photons/cm^2/s and is best fitted by a power law with a photon index of {\Gamma} = 2.54\pm0.14. The pulsed {\gamma}-ray flux in the same energy range is F_100 = (3.95\pm0.07)\times10^(-7) photons/cm^2/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of {\Gamma} = 1.41 \pm 0.23 and a cutoff energy E_c = 4.04 \pm 0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the {\gamma}-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

We have used the Spitzer Space Telescope Infrared Spectrograph (IRS) to observe the 5-37 micron thermal emission of comet 73P/Schwassmann-Wachmann 3 (SW3), components B and C. We obtained low spectral resolution (R ~ 100) data over the entire wavelength interval, along with images at 16 and 22 micron. These observations provided an unprecedented opportunity to study nearly pristine material from the surface and what was until recently the interior of an ecliptic comet - cometary surface having experienced only two prior perihelion passages, and including material that was totally fresh. The spectra were modeled using a variety of mineral types including both amorphous and crystalline components. We find that the degree of silicate crystallinity, ~ 35%, is somewhat lower than most other comets with strong emission features, while its abundance of amorphous carbon is higher. Both suggest that SW3 is among the most chemically primitive solar system objects yet studied in detail, and that it formed earlier or farther from the sun than the bulk of the comets studied so far. The similar dust compositions of the two fragments suggests that these are not mineralogically heterogeneous, but rather uniform throughout their volumes. Atomic abundances derived from the spectral models indicates a depletion of O compared to solar photospheric values, despite the inclusion of water ice and gas in the models. Atomic C may be solar or slightly sub-solar, but its abundance is complicated by the potential contribution of spectrally featureless mineral species to the portion of the spectra most sensitive to the derication of the C abundance. We find a relatively high bolometric albedo, ~ 0.13 for the dust, considering the large amount of dark carbonaceous material, but consistent with the presence of abundant small particles and strong emission features.

We present new Spitzer/IRS spectra of two hot R Coronae Borealis (RCB) stars, one in the Galaxy,V348 Sgr, and one lying in the LMC, HV 2671. These two objects may constitute a link between the RCB stars and the late Wolf-Rayet ([WCL]) class of central stars of planetary nebula (CSPNe) such as CPD -56 8032 that has little or no hydrogen in their atmospheres. HV 2671 and V348 Sgr are members of a rare subclass that has significantly higher effective temperatures than most RCB stars, but sharing the traits of hydrogen deficiency and dust formation that define the cooler RCB stars. The [WC] CSPNe star, CPD -56 8032, displays evidence for dual-dust chemistry showing both PAHs and crystalline silicates in its mid-IR spectrum. HV 2671 shows strong PAH emission but shows no sign of having crystalline silicates. The spectrum of V348 Sgr is very different from those of CPD -56 8032 and HV 2671. The PAH emission seen strongly in the other two stars is not present. Instead, the spectrum is dominated by a broad emission centered at about 8.2 micron. The mid-IR spectrum of CPD -56 8032 shows emission features that may be associated with C60. The other two stars do not show evidence for C60. HV 2671 has also been detected by Herschel/PACS and SPIRE. V348 Sgr and CPD -56 8032 have been detected by AKARI/FIS. These data were combined with Spitzer, IRAS, 2MASS and other photometry to produce their spectral energy distributions from the visible to the far-IR. Monte Carlo radiative transfer modeling was used to study the circumstellar dust around these stars. HV 2671 and CPD -56 8032 require both a flared inner disk with warm dust and an extended diffuse envelope with cold dust to to fit their SEDs. The SED of V348 Sgr can be fit with a much smaller disk and envelope.

The time variations of spectral properties of dark martian surface features are investigated using the OMEGA near-IR dataset. The analyzed period covers two Mars years, spanning from early 2004 to early 2008 (includes the 2007 global dust event). Radiative transfer modeling indicates that the apparent albedo variations of low to mid-latitude dark regions are consistent with those produced by the varying optical depth of atmospheric dust as measured simultaneously from the ground by the Mars Exploration Rovers. We observe only a few significant albedo changes that can be attributed to surface phenomena. They are small-scaled and located at the boundaries between bright and dark regions. We then investigate the variations of the mean particle size of aerosols using the evolution of the observed dark region spectra between 1 and 2.5 {\mu}m. Overall, we find that the observed changes in the spectral slope are consistent with a mean particle size of aerosols varying with time between 1 and 2 {\mu}m. Observations with different solar zenith angles make it possible to characterize the aerosol layer at different altitudes, revealing a decrease of the particle size of aerosols as altitude increases.

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nan\c{c}ay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

In late 2008, the quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 started a transition from outburst to quiescence, after it had been actively accreting for more than 24 years. In a previous work, we discussed Chandra and Swift observations obtained during the first five months after this transition. Here, we report on further X-ray observations of EXO 0748-676, extending the quiescent monitoring to 1.6 years. Chandra and XMM-Newton data reveal quiescent X-ray spectra composed of a soft, thermal component that is well-fitted by a neutron star atmosphere model. An additional hard powerlaw tail is detected that changes non-monotonically over time, contributing between 4 and 20 percent to the total unabsorbed 0.5-10 keV flux. The combined set of Chandra, XMM-Newton and Swift data reveals that the thermal bolometric luminosity fades from ~1E34 to 6E33 (D/7.4 kpc)^2 erg/s, whereas the inferred neutron star effective temperature decreases from ~124 to 109 eV. We interpret the observed decay as cooling of the neutron star crust and show that the fractional quiescent temperature change of EXO 0748-676 is markedly smaller than observed for three other neutron star X-ray binaries that underwent prolonged accretion outbursts.

The Advanced X-ray Timing Array (AXTAR) is a mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument, the Large Area Timing Array (LATA) is a collimated instrument with 2-50 keV coverage and over 3 square meters effective area. The LATA is made up of an array of supermodules that house 2-mm thick silicon pixel detectors. AXTAR will provide a significant improvement in effective area (a factor of 7 at 4 keV and a factor of 36 at 30 keV) over the RXTE PCA. AXTAR will also carry a sensitive Sky Monitor (SM) that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and technical concept for AXTAR and present results from a preliminary mission design study.

We report the discovery of eight gamma-ray pulsars in blind frequency searches using the LAT, onboard the Fermi Gamma-ray Space Telescope. Five of the eight pulsars are young (tau_c<100 kyr), energetic (Edot>10^36 erg/s), and located within the Galactic plane (|b|<3 deg). The remaining three are older, less energetic, and located off the plane. Five pulsars are associated with sources included in the LAT bright gamma-ray source list, but only one, PSR J1413-6205, is clearly associated with an EGRET source. PSR J1023-5746 has the smallest characteristic age (tau_c=4.6 kyr) and is the most energetic (Edot=1.1E37 erg/s) of all gamma-ray pulsars discovered so far in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages (tau_c~1 Myr) and are the least energetic (Edot~5E33 erg/s) of the newly-discovered pulsars. We present the timing models, light curves, and detailed spectral parameters of the new pulsars. We used recent XMM observations to identify the counterpart of PSR J2055+25 as XMMU J205549.4+253959. In addition, publicly available archival Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023-5746 as a faint, highly absorbed source, CXOU J102302.8-574606. The large X-ray absorption indicates that this could be among the most distant gamma-ray pulsars detected so far. PSR J1023-5746 is positionally coincident with the TeV source HESS J1023-575, located near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident with a 4.3 sigma excess reported by Milagro at a median energy of 35 TeV. Deep radio follow-up observations of the eight pulsars resulted in no detections of pulsations and upper limits comparable to the faintest known radio pulsars, indicating that these can be included among the growing population of radio-quiet pulsars in our Galaxy being uncovered by the LAT, and currently numbering more than 20.

We report our complete database of X-ray eclipse timings of the low mass X-ray binary EXO0748-676 observed by the Rossi X-Ray Timing Explorer (RXTE) satellite. As of this writing we have accumulated 443 full X-ray eclipses, 392 of which have been observed with the Proportional Counter Array on RXTE. These include both observations where an eclipse was specifically targeted and those eclipses found in the RXTE data archive. Eclipse cycle count has been maintained since the discovery of the EXO0748-676 system in February 1985. We describe our observing and analysis techniques for each eclipse and describe improvements we have made since the last compilation by Wolff et al. (2002). The principal result of this paper is the database containing the timing results from a seven-parameter fit to the X-ray light curve for each observed eclipse along with the associated errors in the fitted parameters. Based on the standard O-C analysis, EXO0748-676 has undergone four distinct orbital period epochs since its discovery. In addition, EXO0748-676 shows small-scale events in the O-C curve that are likely due to short-lived changes in the secondary star.

Million degree plasmas are ubiquitous in the Universe, and examples include the atmospheres of white dwarfs; accretion phenomena in young stars, cataclysmic variables and active galactic nuclei; the coronae of stars; and the interstellar medium of our own galaxy and of others. The bulk of radiation from million degree plasmas is emitted at extreme ultraviolet (EUV) wavelengths, which includes critical spectral features containing diagnostic information often not available at other wavelengths. With underpinning by a mature instrument technology, there is great opportunity here for exciting discoveries.

The observations of Mars by the CRISM and OMEGA hyperspectral imaging spectrometers require correction for photometric, atmospheric and thermal effects prior to the interpretation of possible mineralogical features in the spectra. Here, we report on a simple, yet non-trivial, adaptation to the commonly-used volcano-scan correction technique for atmospheric CO_2, which allows for the improved detection of minerals with intrinsic absorption bands at wavelengths between 1.9-2.1 $\mu$m. This volcano-scan technique removes the absorption bands of CO_2 by ensuring that the Lambert albedo is the same at two wavelengths: 1.890 $\mu$m and 2.011 $\mu$m, with the first wavelength outside the CO_2 gas bands and the second wavelength deep inside the CO_2 gas bands. Our adaptation to the volcano-scan technique moves the first wavelength from 1.890 $\mu$m to be instead within the gas bands at 1.980 $\mu$m, and for CRISM data, our adaptation shifts the second wavelength slightly, to 2.007 $\mu$m. We also report on our efforts to account for a slight ~0.001 $\mu$m shift in wavelengths due to thermal effects in the CRISM instrument.

We discuss the DISORT-based radiative transfer pipeline ('CRISM_LambertAlb') for atmospheric and thermal correction of MRO/CRISM data acquired in multispectral mapping mode (~200 m/pixel, 72 spectral channels). Currently, in this phase-one version of the system, we use aerosol optical depths, surface temperatures, and lower-atmospheric temperatures, all from climatology derived from Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) data, and surface altimetry derived from MGS Mars Orbiter Laser Altimeter (MOLA). The DISORT-based model takes as input the dust and ice aerosol optical depths (scaled to the CRISM wavelength range), the surface pressures (computed from MOLA altimetry, MGS-TES lower-atmospheric thermometry, and Viking-based pressure climatology), the surface temperatures, the reconstructed instrumental photometric angles, and the measured I/F spectrum, and then outputs a Lambertian albedo spectrum. The Lambertian albedo spectrum is valuable geologically since it allows the mineralogical composition to be estimated. Here, I/F is defined as the ratio of the radiance measured by CRISM to the solar irradiance at Mars divided by $\pi$. After discussing the capabilities and limitations of the pipeline software system, we demonstrate its application on several multispectral data cubes: the outer northern ice cap of Mars, Tyrrhena Terra, and near the landing site for the Phoenix mission. For the icy spectra near the northern polar cap, aerosols need to be included in order to properly correct for the CO_2 absorption in the H_{2}O ice bands at wavelengths near 2.0 $\mu$m. In future phases of software development, we intend to use CRISM data directly in order to retrieve the spatiotemporal maps of aerosol optical depths, surface pressure and surface temperature.

The quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 recently started the transition to quiescence following an accretion outburst that lasted more than 24 years. We report on two Chandra and twelve Swift observations performed within five months after the end of the outburst. The Chandra spectrum is composed of a soft, thermal component that fits to a neutron star atmosphere model with kT^inf~0.12 keV, joined by a hard powerlaw tail that contributes ~20% of the total 0.5-10 keV unabsorbed flux. The combined Chandra/Swift data set reveals a relatively hot and luminous quiescent system with a temperature of kT^inf~0.11-0.13 keV and a bolometric thermal luminosity of ~8.1E33-1.6E34 (d/7.4 kpc)^2 erg/s. We discuss our results in the context of cooling neutron star models.

We report evidence of magnetic activity associated with the secondary star in the EXO 0748-676 low mass X-ray binary system. An analysis of a sequence of five consecutive X-ray eclipses observed during December 2003 with the RXTE satellite brings out a feature occurring during ingress we interpret as the X-ray photoelectric absorption shadow, as seen by an observer at Earth, of a plasma structure suspended above the surface of the secondary star. The light curve feature consists of an initial drop in count rate to near zero (the absorption shadow) with a very short rebound to a significant fraction of the pre-ingress count rate and then a final plunge to totality over a total time scale of ~25 s. The ingress feature persists for at least 5 consecutive orbital periods (a total of ~19 hr), and possibly up to 5 days in our data. Our data also show significant post-egress dipping during this eclipse sequence, unusual for this source, indicating possible secondary star mass ejection during this episode.

Accretion-powered X-ray pulsars are among the most luminous X-ray sources in the Galaxy. However, despite decades of theoretical and observational work since their discovery, no satisfactory model for the formation of the observed X-ray spectra has emerged. In this paper, we report on a self-consistent calculation of the spectrum emerging from a pulsar accretion column that includes an explicit treatment of the bulk and thermal Comptonization occurring in the radiation-dominated shocks that form in the accretion flows. Using a rigorous eigenfunction expansion method, we obtain a closed-form expression for the Green's function describing the upscattering of monochromatic radiation injected into the column. The Green's function is convolved with bremsstrahlung, cyclotron, and blackbody source terms to calculate the emergent photon spectrum. We show that energization of photons in the shock naturally produces an X-ray spectrum with a relatively flat continuum and a high-energy exponential cutoff. Finally, we demonstrate that our model yields good agreement with the spectra of the bright pulsar Her X-1 and the low luminosity pulsar X Per.

We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.

We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron images at two different spatial resolutions, and with radio continuum maps. We determine an accurate median ratio of IRAC 8.0-micron/MSX\8.3-micron fluxes, of 1.55+/-0.15. From robust spectral energy distributions of these regions we conclude that the present 8.0-micron diffuse calibration of the SST is 36% too high compared with the MSX validated calibration, perhaps due to scattered light inside the camera. This is an independent confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science Center (SSC). From regression analyses we find that 843-MHz radio fluxes of HII regions and mid-infrared (MIR) fluxes are linearly related for MSX at 8.3-microns and Spitzer at 8.0 microns, confirming the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse continuum fluxes is 600 times smaller in nonthermal than thermal regions, making it a sharp discriminant. The ratios are largely independent of morphology up to a size of ~24 arcsec. We provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse haloes, perhaps infrared counter- parts to weakly ionized radio haloes found around compact HII regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse HII regions.

We develop a new theoretical model for the spectral formation process in accretion-powered X-ray pulsars based on a detailed treatment of the bulk and thermal Comptonization occurring in the accreting, shocked gas. A rigorous eigenfunction expansion method is employed to obtain the analytical solution for the Green's function describing the scattering of radiation injected into the column from a monochromatic source located at an arbitrary height above the stellar surface. The emergent spectrum is calculated by convolving the Green's function with source terms corresponding to bremsstrahlung, cyclotron, and blackbody emission. The energization of the photons in the shock, combined with cyclotron absorption, naturally produces an X-ray spectrum with a relatively flat continuum shape and a high-energy quasi-exponential cutoff. We demonstrate that the new theory successfully reproduces the phase-averaged spectra of the bright pulsars Her X-1, LMC X-4, and Cen X-3. In these luminous sources, it is shown that the emergent spectra are dominated by Comptonized bremsstrahlung emission.

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation.

We have observed an unusually strong X-ray burst as a part of our regular eclipse timing observations of the low mass binary system EXO0748-676. The burst peak flux was 5.2x10^-8 ergs cm^-2 s^-1, approximately five times the normal peak X-ray burst flux observed from this source by RXTE. Spectral fits to the data strongly suggest that photospheric radius expansion occurred during the burst. In this Letter we examine the properties of this X-ray burst, which is the first example of a radius expansion burst from EXO0748-676 observed by RXTE. We find no evidence for coherent burst oscillations. Assuming that the peak burst luminosity is the Eddington luminosity for a 1.4 solar mass neutron star we derive a distance to EXO0748-676 of 7.7 kpc for a helium-dominated burst photosphere and 5.9 kpc for a hydrogen-dominated burst photosphere.

We present an extinction analysis of 9 paths through the LMC and SMC based on FUSE observations. To date, just two LMC sight lines have probed dust grain composition and size distributions in the Clouds using spectra including wavelengths as short as 950 A. We supplement these with results from 4 regions distinguished by their IR through UV extinction curves and grouped as LMCAvg, LMC2, SMC bar and SMC wing. Despite the distinct characters of Milky Way and Magellanic Cloud extinction, our results are generally analogous to those found for Galactic curves in that the FUSE portions of each extinction curve are described reasonably well by FM curves fitted only to longer wavelength data and lack any dramatic new extinction features, and any deviations from the CCM formalism continue into FUV wavelengths. An MEM analysis of these curves suggests that LMCAvg and SMC wing sight lines require more silicon and/or carbon in dust than current abundance measurements would allow, while the requirements for LMC2 and SMC bar sight lines do not fully tax the available reservoirs. An intermediate product of this analysis is the measurement of new H_2 abundances in the Magellanic Clouds. Collectively considering Cloud sight lines that possess significant H_2 column densities, E(B-V)/N(HI) ratios are reduced by significant factors relative to the Galactic mean, whereas the corresponding E(B-V)/N(H_2) values more closely resemble their Galactic counterpart. These trends reflect the fact that among these sight lines f(H_2)-values are lower than those common in the Milky Way for paths with similar degrees of reddening.

Accretion-powered X-ray pulsars are among the most luminous X-ray sources in the Galaxy. However, despite decades of theoretical and observational work since their discovery, no satisfactory model for the formation of the observed X-ray spectra has emerged. In particular, the previously available theories are unable to reproduce the power-law variation observed at high energies in many sources. In this paper, we present the first self-consistent calculation of the spectrum emerging from a pulsar accretion column that includes an explicit treatment of the energization occurring in the shock. Using a rigorous eigenfunction expansion method based on the exact dynamical solution for the velocity profile in the column, we obtain a closed-form expression for the Green's function describing the upscattering of radiation injected into the column from a monochromatic source located at the top of the thermal mound, near the base of the flow. The Green's function is convolved with a Planck distribution to calculate the radiation spectrum resulting from the reprocessing of blackbody photons emitted by the thermal mound. We demonstrate that the energization of the photons in the shock naturally produces an X-ray spectrum with a power-law shape at high energies and a blackbody shape at low energies, in agreement with many observations of accreting X-ray pulsars.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emission as thermal, and the probable presence of polycylic aromatic hydrocarbons in and around the source. These lead to the conclusion that G313.3+00.3 is a planetary nebula. This object is of particular interest because it was discovered solely through radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts as a proof of concept for the discovery of many more highly extinguished planetary nebulae. G313.3+00.3 is well-resolved by both the instruments with which it was identified, and suffers extreme reddening due to its location in the Scutum-Crux spiral arm.

We present the first self-consistent model for the dynamics and the radiative transfer occurring in bright X-ray pulsar accretion columns, with a special focus on the role of the shock in energizing the emerging X-rays. The pressure inside the accretion column of a luminous X-ray pulsar is dominated by the photons, and consequently the equations describing the coupled radiative-dynamical structure must be solved simultaneously. Spectral formation in these sources is therefore a complex, nonlinear phenomenon. We obtain the analytical solution for the Green's function describing the upscattering of monochromatic radiation injected into the column from the thermal mound located near the base of the flow. The Green's function is convolved with a Planck distribution to model the X-ray spectrum resulting from the reprocessing of blackbody photons produced in the thermal mound. These photons diffuse through the infalling gas and eventually escape out the walls of the column, forming the observed X-ray spectrum. We show that the resulting column-integrated, phase-averaged spectrum has a power-law shape at high energies and a blackbody shape at low energies, in agreement with the observational data for many X-ray pulsars.

We present extinction curves that include data down to far ultraviolet wavelengths (FUV; 1050 - 1200 A) for nine Galactic sight lines. The FUV extinction was measured using data from the Far Ultraviolet Spectroscopic Explorer. The sight lines were chosen for their unusual extinction properties in the infrared through the ultraviolet; that they probe a wide range of dust environments is evidenced by the large spread in their measured ratios of total-to-selective extinction, R_V = 2.43 - 3.81. We find that extrapolation of the Fitzpatrick & Massa relationship from the ultraviolet appears to be a good predictor of the FUV extinction behavior. We find that predictions of the FUV extinction based upon the Cardelli, Clayton & Mathis (CCM) dependence on R_V give mixed results. For the seven extinction curves well represented by CCM in the infrared through ultraviolet, the FUV extinction is well predicted in three sight lines, over-predicted in two sight lines, and under-predicted in 2 sight lines. A Maximum Entropy Method analysis using a simple three component grain model shows that seven of the nine sight lines in the study require a larger fraction of grain materials to be in dust when FUV extinction is included in the models. Most of the added grain material is in the form of small (radii < 200 A) grains.

We critique the method of constructing extinction curves of lensing galaxies using multiply imaged QSOs. If one of the two QSO images is lightly reddened or if the dust along both sightlines has the same properties then the method works well and produces an extinction curve for the lensing galaxy. These cases are likely rare and hard to confirm. However, if the dust along each sightline has different properties then the resulting curve is no longer a measurement of extinction. Instead, it is a measurement of the difference between two extinction curves. This "lens difference curve'' does contain information about the dust properties, but extracting a meaningful extinction curve is not possible without additional, currently unknown information. As a quantitative example, we show that the combination of two Cardelli, Clayton, & Mathis (CCM) type extinction curves having different values of R(V) will produce a CCM extinction curve with a value of R(V) which is dependent on the individual R(V) values and the ratio of V band extinctions. The resulting lens difference curve is not an average of the dust along the two sightlines. We find that lens difference curves with any value of R(V), even negative values, can be produced by a combination of two reddened sightlines with different CCM extinction curves with R(V) values consistent with Milky Way dust (2.1 < R(V) < 5.6). This may explain extreme values of R(V) inferred by this method in previous studies. But lens difference curves with more normal values of R(V) are just as likely to be composed of two dust extinction curves with R(V) values different than that of the lens difference curve. While it is not possible to determine the individual extinction curves making up a lens difference curve, there is information about a galaxy's dust contained in the lens difference curves.

We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

Orbital period changes are an important diagnostic for understanding low mass X-ray binary (LMXB) accretion-induced angular momentum exchange and overall system evolution. We present our most recent results for the eclipse timing of the LMXB EXO0748-676. Since its discovery in 1985 it has apparently undergone three distinct orbital period "epochs", each characterized by a different orbital period than the previous epoch. We outline the orbital period behavior for EXO0748-676 over the past 18 years and discuss the implications of this behavior in light of current theoretical ideas for LMXB evolution.

The wavelength dependence of interstellar polarization due to dust in M31 has been observed along four sightlines. Only one sightline had been measured previously.The globular clusters, S78, S150, S233 and Baade 327 were used as point sources to probe the interstellar dust in M31. The Serkowski law produces good fits for all the sightlines although the relationship between K and lambda(max) may be different from that found in the Galaxy. The results of this study imply that the slope K/lambda(max) may be significantly larger in M31. The Serkowski curves are significantly narrower than those of the same lambda(max) in the Galaxy and may require extreme modifications to the size distributions of silicate particles. The fits for the four sightlines reveal values of lambda(max) ranging from 4800 to 5500 A. These are consistent with average values of lambda(max) measured in the Galaxy and the Magellanic Clouds. The range measured for M31 corresponds to R(V) values of 2.7 to 3.1. The range in R(V) seen in the Galaxy is 2.5 to 5.5 implying, for this small sample, that the average size ofinterstellar grains in M31 is typically smaller than that seen for Galactic grains if the nature of the grains is the same. Also, the polarization efficiency for these sightlines is large although some bias is expected since sightlines known to have significant interstellar polarization were selected for the sample.

Narayan and Heyl (2002) have developed a theoretical framework to convert suitable upper limits on type I X-ray bursts from accreting black hole candidates (BHCs) into evidence for an event horizon. However, no appropriate observational limit exists in the literature. In this paper we survey 2101.2 ks of data from the Unconventional Stellar Aspect (USA) X-ray timing experiment and 5142 ks of data from the Rossi X-ray Timing Explorer (RXTE) experiment to obtain a formal constraint of this type. 1122 ks of neutron star data yield a population averaged mean burst rate of 1.69 x 10-5 bursts s-1 while 6081 ks of BHC data yield a 95% confidence level upper limit of 4.9 x 10-7 bursts s-1. This is the first published limit of this type for Black Hole Candidates. Applying the theoretical framework of Narayan and Heyl (2002) we calculate regions of unstable luminosity where the neutron stars are expected to burst and the BHCs would be expected to burst if they had a surface. In this unstable luminosity region 464 ks of neutron star data yield an averaged mean burst rate of 4.1 x 10-5 bursts s-1 and 1512 ks of BHC data yield a 95% confidence level upper limit of 2.0 x 10-6 bursts s-1, and a limit of > 10\sigma that BHCs do not burst with a rate similar to the rate of neutron stars in these unstable regions. This gives further evidence that BHCs do not have surfaces unless there is some new physics occurring on their surface.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

We present an exhaustive, quantitative comparison of all of the known extinction curves in the Small and Large Magellanic Clouds (SMC and LMC) with our understanding of the general behavior of Milky Way extinction curves. The R_V dependent CCM relationship and the sample of extinction curves used to derive this relationship is used to describe the general behavior of Milky Way extinction curves. The ultraviolet portion of the SMC and LMC extinction curves are derived from archival IUE data, except for one new SMC extinction curve which was measured using HST/STIS observations. The optical extinction curves are derived from new (for the SMC) and literature UBVRI photometry (for the LMC). The near-infrared extinction curves are calculated mainly from 2MASS photometry supplemented with DENIS and new JHK photometry. For each extinction curve, we give R_V = A(V)/E(B-V) and N(HI) values which probe the same dust column as the extinction curve. We compare the properties of the SMC and LMC extinction curves with the CCM relationship three different ways: each curve by itself, the behavior of extinction at different wavelengths with R_V, and behavior of the extinction curve FM fit parameters with R_V. As has been found previously, we find that a small number of LMC extinction curves are consistent with the CCM relationship, but majority of the LMC and all of the SMC curves do not follow the CCM relationship. For the first time, we find that the CCM relationship seems to form a bound on the properties of all of the LMC and SMC extinction curves. This result strengthens the picture of dust extinction curves exhibit a continuum of properties between those found in the Milky Way and the SMC Bar. (abridged)

We have obtained new STIS/HST spectra to search for structure in the ultraviolet interstellar extinction curve, with particular emphasis on a search for absorption features produced by polycyclic aromatic hydrocarbons (PAHs). The presence of these molecules in the interstellar medium has been postulated to explain the infrared emission features seen in the 3-13 $\mu$m spectra of numerous sources. UV spectra are uniquely capable of identifying specific PAH molecules. We obtained high S/N UV spectra of stars which are significantly more reddened than those observed in previous studies. These data put limits on the role of small (30-50 carbon atoms) PAHs in UV extinction and call for further observations to probe the role of larger PAHs. PAHs are of importance because of their ubiquity and high abundance inferred from the infrared data and also because they may link the molecular and dust phases of the interstellar medium. A presence or absence of ultraviolet absorption bands due to PAHs could be a definitive test of this hypothesis. We should be able to detect a 20 \AA wide feature down to a 3$\sigma$ limit of $\sim$0.02 A$_V$. No such absorption features are seen other than the well-known 2175 \AA bump.

We present 2-D radiation transfer models of Class I Protostars and show the effect of including more realistic geometries on the resulting spectral energy distributions and images. We begin with a rotationally flattened infalling envelope as our comparison model, and add a flared disk and bipolar cavity. The disk affects the spectral energy distribution most strongly at edge-on inclinations, causing a broad dip at about 10 um (independent of the silicate feature) due to high extinction and low scattering albedo in this wavelength region. The bipolar cavities allow more direct stellar+disk radiation to emerge into polar directions, and more scattering radiation to emerge into all directions. The wavelength-integrated flux, often interpreted as luminosity, varies with viewing angle, with pole-on viewing angles seeing 2-4 times as much flux as edge-on, depending on geometry. Thus, observational estimates of luminosity should take into account the inclination of a source. The envelopes with cavities are significantly bluer in near-IR and mid-IR color-color plots than those without cavities. Using 1-D models to interpret Class I sources with bipolar cavities would lead to an underestimate of envelope mass and an overestimate of the implied evolutionary state. We compute images at near-, mid-, and far-IR wavelengths. We find that the mid-IR colors and images are sensitive to scattering albedo, and that the flared disk shadows the midplane on large size scales at all wavelengths plotted. Finally, our models produce polarization spectra which can be used to diagnose dust properties, such as albedo variations due to grain growth. Our results of polarization across the 3.1 um ice feature agree well with observations for ice mantles covering 5% of the radius of the grains.

We present an analysis of long term X-ray monitoring observations of Circinus X-1 (Cir X-1) made with four different instruments: Vela 5B, Ariel V ASM, Ginga ASM, and RXTE ASM, over the course of more than 30 years. We use Lomb-Scargle periodograms to search for the ~16.5 day orbital period of Cir X-1 in each of these data sets and from this derive a new orbital ephemeris based solely on X-ray measurements, which we compare to the previous ephemerides obtained from radio observations. We also use the Phase Dispersion Minimization (PDM) technique, as well as FFT analysis, to verify the periods obtained from periodograms. we obtain dynamic periodograms (both Lomb-Scargle and PDM) of Cir X-1 during the RXTE era, showing the period evolution of Cir X-1, and also displaying some unexplained discrete jumps in the location of the peak power.

Employing the Maximum Entropy Method algorithm, we fit interstellar extinction measurements which span the wavelength range 0.125-3 micron. We present a uniform set of MEM model fits, all using the same grain materials, optical constants and abundance constraints. In addition, we are taking advantage of improved UV and IR data and better estimates of the gas-to-dust ratio. The model fits cover the entire range of extinction properties that have been seen in the Galaxy and the Magellanic Clouds. The grain models employed for this presentation are the simplistic homogeneous spheres models (i.e., Mathis, Rumpl, & Nordsieck 1977) with two (graphite, silicate) or three (graphite, silicate, amorphous carbon) components. Though such usage is only a first step, the results do provide interesting insight into the use of grain size as a diagnostic of dust environment. We find that the SMC Bar extinction curve cannot be fit using carbon grains alone. This is a challenge to the recent observational result indicating little silicon depletion in the SMC.

We have analyzed observations of the Acidalia hemisphere of Mars taken by the Hubble Space Telescope's Near-Infrared Camera Multi-Object Spectrograph (HST/NICMOS) during July of 1997 (Ls = 152 degrees, northern Martian summer). The data consist of images at ~60 km/pixel resolution, using both narrow- and medium-band filters specifically selected to allow us to study the hydration state of the Martian surface. Calibration was performed by comparison to Phobos-2 ISM observations of overlapping regions, and atmospheric gas correction was performed by modeling the atmosphere for each pixel using a line-by-line radiative transfer code coupled with the MOLA altimetry data. Our results indicate the presence of at least three spectrally different large-scale (>1000 km diameter) terrains corresponding to the dark regions of northern Acidalia, the southern hemisphere classical dark terrain, and the classical intermediate terrain adjacent to southern Acidalia. We also identified two other spectrally unique terrains, corresponding to the northern polar ice cap, and to the southern winter polar hood. Comparisons with mineral spectra indicate the possibility of different H2O- or OH-bearing (i.e., hydroxides and/or hydrates) minerals existing both in northern Acidalia and in the nearby intermediate albedo terrain. Hydrated minerals do not appear to be spectrally important components of the southern hemisphere dark terrains imaged by HST in 1997.

Space Telescope Imaging Spectrograph (STIS) and archival International Ultraviolet Explorer (IUE) data were used to search for structure in the ultraviolet interstellar extinction curve. Between 1250 and 3200 A, in addition to the well-known 2175 A bump, there is only one fairly small (10%) absorption feature at 2700 A. The feature is very broad (>400 A). It coincides extremely well in position and shape with laboratory measurements of presolar diamond grains found in meteorites. As this feature appears in both the new STIS and archival IUE data, it is unlikely to be an instrumental effect. This feature appears to be ubiquitous in the ISM as it is seen along sightlines toward various clouds with reddenings in the range of E(B-V) = 0.35 to 1.73 mag. The identification with diamonds remains tentative since it is based on a single, weak feature. If the feature is due to interstellar diamonds then they require an order of magnitude more absorption strength per gram/cm^2 than the presolar meteoritic diamonds measured in the lab.

We present 7 eclipse timings of the low mass X-ray binary EXO0748-676 obtained with the USA experiment during 1999-2000 as well as 122 eclipse timings obtained with RXTE during 1996-2000. According to our analysis, the mean orbital period has increased by ~8 ms between the pre-RXTE era (1985-1990) and the RXTE/USA era (1996-2000). This corresponds to an orbital period derivative of P(orb)/(dP(orb)/dt)~2x10^7 years. However, neither a constant orbital period derivative nor any other simple ephemeris provides an acceptable fit to the data: individual timings of eclipse centers have residuals of up to 15 or more seconds away from our derived smooth ephemerides. When we consider all published eclipse timing data including those presented here, a model that includes observational measurement error, cumulative period jitter, and underlying period evolution is found to be consistent with the timing data. We discuss several physical mechanisms for LMXB orbital evolution in an effort to account for the change in orbital period and the observed intrinsic jitter in the mid-eclipse times.

This paper reports X-ray spectral observations of a relatively nearby z = 0.048 BL Lacertae (BL Lac) object 1ES1959+650, which is a potential TeV emitter. The observations include 31 short pointings made by the Unconventional Stellar Aspect (USA) Experiment on board the Advanced Research and Global Observation Satellite (ARGOS), and 17 pointings by the PCA on board the Rossi X-ray Timing Explorer (RXTE). Most of these observations were spaced by less than 1 day. \es was detected by the ARGOS USA detector in the range 1-16 keV, and by the PCA in the 2-16 keV range but at different times. During the closely spaced RXTE observations beginning on 2000 July 28, an ending of one flare and a start of another are visible, associated with spectral changes, where the photon index Gamma ranges between ~ 1.4 and 1.7, and the spectrum is harder when the source is brighter. This implies that 1ES1959 is an XBL-type blazar, with the X-ray emission likely to originate via the synchrotron process. The USA observations reveal another flare that peaked on 2000 November 14 and doubled the flux within a few days, again associated with spectral changes of the same form. The spectral variability correlated with the flux and timing characteristics of this object that are similar to those of other nearby BL Lacs, and suggest relativistic beaming with a Doppler factor delta > 1.6 and magnetic fields of the order of a few mG. We also suggest that the steady component of the X-ray emission -- present in this object as well as in other XBLs -- may be due to the large-scale relativistic jet (such as measured by Chandra in many radio-loud AGN), but pointing very closely to our line of sight.

We investigate the observational signatures of a straightforward evolutionary scenario for protoplanetary disks in which the disk mass of small (50 micron) particles decreases homologously with time, but the disk structure and stellar parameters do not change. Our goal is to identify optimal infrared spectral indicators of the existence of disks, their structure, and mass evolution that may be tested with the upcoming SIRTF mission. We present simulated spectral energy distributions and colors over a wide range of masses. The SED is most sensitive to disk mass in the far-IR and longer wavelengths, which is already known from millimeter and radio observations. As the disk mass decreases, the excess emission of the disk over the stellar photosphere diminishes more rapidly at the longest rather than at short wavelengths. At near-infrared wavelengths, the disk remains optically thick to stellar radiation over a wide range of disk mass, resulting in a slower decline in the SED in this spectral regime. Therefore, near-IR excesses (K-L) provide a robust means of detecting disks in star clusters down to 1E-7 solar masses, while the far-IR excess probes the disk mass. Reducing the disk mass results in a clear progression in color-color diagrams with low mass disks displaying the bluest colors. We interpret color-color diagrams for Taurus-Auriga sources in the context of decreasing disk mass.

We present spectral energy distribution (SED) models for the edge-on classical T Tauri star HH30 IRS that indicate dust grains have grown to larger than 50 microns within its circumstellar disk. The disk geometry and inclination are known from previous modeling of multiwavelength Hubble Space Telescope images and we use the SED to constrain the dust size distribution. Model spectra are shown for different circumstellar dust models: a standard ISM mixture and larger grain models. As compared to ISM grains, the larger dust grain models have a shallower wavelength dependent opacity. Models with the larger dust grains provide a good match to the currently available data, but mid and far-IR observations are required to more tightly constrain the dust size distribution. The accretion luminosity in our models is L_acc<0.2 L_star corresponding to an accretion rate of 4E-9M_sun/yr. Dust size distributions that are simple power-law extensions (i.e., no exponential cutoff) yield acceptable fits to the optical/near-IR but too much emission at mm wavelengths and require larger disk masses. Such a simple size distribution would not be expected in an environment such as the disk of HH30 IRS, particularly over such a large range in grain sizes. However, its ability to adequately characterize the grain populations may be determined from more complete observational sampling of the SED in the mid to far-IR.

We present a linear diffusion model for the evolution of the double-peaked outburst in the transient source XTEJ1118+480. The model treats the two outbursts as episodic mass deposition at the outer radius of the disk followed by evolution of disk structure according to a diffusion process. We demonstrate that light curves with fast-rise, exponential decay profile are a general consequence of the diffusion process. Deconvolution of the light curve proves to be feasible and gives an input function specifying mass deposition at the outer disk edge as well as the total mass of the disk, both as functions of time. The derived evolution of total disk mass can be correlated with the observed evolution of the ~0.1 Hz QPO in the source reported in Wood et al. (2000).

We report on timing and spectral observations of the 2000 outburst of XTE J1550-564 made by the Unconventional Stellar Aspect (USA) Experiment on board the Advanced Research and Global Observation Satellite (ARGOS). We observe a low-frequency quasi-periodic oscillation (LFQPO) with a centroid frequency that tends to increase with increasing flux and a fractional rms amplitude which is correlated with the hardness ratio. The evolution of the hardness ratio (4--16 keV/1--4 keV) with time and source flux is examined. The hardness-intensity diagram (HID) shows a cyclical movement in the clockwise direction and possibly indicates the presence of two independent accretion flows. We observe a relationship between the USA 4--16 keV count rate and radio observations and discuss this in the context of previously observed correlations between X-ray, radio, optical and IR data. We examine our results in the context of models invoking two accretion flows: a thin disk and a hot sub-Keplerian flow.

Ground based imaging, imaging polarimetry, and recent Hubble Space Telescope WFPC2 and NICMOS images of protostars have revealed very complex scattered light patterns that cannot be entirely explained by 2-D radiation transfer models. We present here for the first time radiation transfer models of T Tau and IRAS 04016+2610 that are fully 3-D, with the aim of investigating the effects on image morphology of multiple illuminating sources and infalling envelopes that have been shaped by multiple outflows. For T Tau we have constructed scattered light models where the illumination of the surrounding envelope is by a binary with each source surrounded by its own small circumstellar disk or envelope. We find that the asymmetries in the WFPC2 image of T Tau can be reproduced if the disks in the binary system are misaligned, consistent with a recently discovered bipolar outflow believed to originate from the secondary. For IRAS 04016+2610 we find that the observed scattered light pattern can be reproduced by scattering in an envelope with cavities carved by two sets of bipolar outflows, suggestive of an embedded binary system.

We present Hubble Space Telescope (HST) Near-Infrared Camera and Multi-object Spectrometer (NICMOS) observations of the reflection nebulosity associated with the T Tauri star HH 30. The images show the scattered light pattern characteristic of a highly inclined, optically thick disk with a prominent dustlane whose width decreases with increasing wavelength. The reflected nebulosity exhibits a lateral asymmetry in the upper lobe on the opposite side to that reported in previously published Wide Field Planetary Camera 2 (WFPC2) images. The radiation transfer model which most closely reproduces the data has a flared accretion disk with dust grains larger than standard interstellar medium grains by a factor of approximately 2.1. A single hotspot on the stellar surface provides the necessary asymmetry to fit the images and is consistent with previous modeling of the light curve and images. Photometric analysis results in an estimated extinction of Av>~80; however, since the photometry measures only scattered light rather than direct stellar flux, this a lower limit. The radiative transfer models require an extinction of Av = 7,900.

In this paper and a companion paper we present the DIRTY model, a Monte Carlo radiative transfer code, self-consistently including dust heating and emission, and accounting for the effects of the transient heating of small grains. The code is completely general; the density structure of the dust, the number and type of heating sources, and their geometric configurations can be specified arbitrarily within the model space. Source photons are tracked through the scattering and absorbing medium using Monte Carlo techniques and the effects of multiple scattering are included. The dust scattering, absorbing, and emitting properties are calculated from realistic dust models derived by fitting observed extinction curves in Local Group galaxies including the Magellanic Clouds and the Milky Way. The dust temperature and the emitted dust spectrum are calculated self consistently from the absorbed energy including the effects of temperature fluctuations in small grains. Dust self-absorption is also accounted for, allowing the treatment of high optical depths, by treating photons emitted by the dust as an additional heating source and adopting an iterative radiative transfer scheme. As an illustrative case, we apply the DIRTY radiative transfer code to starburst galaxies wherein the heating sources are derived from stellar evolutionary synthesis models. Within the context of the starburst model, we examine the dependence of the UV to FIR SED, dust temperatures, and dust masses predicted by DIRTY on variations of the input parameters.

The USA experiment on ARGOS and RXTE have exensively observed the X-ray transient XTE J1118+480 during its recent outburst in 2000 April--June. We present detailed monitoring of the evolution of a low frequency QPO which drifts from 0.07 Hz to 0.15 Hz during the outburst. We examine possible correlations of the QPO frequency with the flux and spectral characteristics of the source, and compare this QPO to low frequency QPOs observed in other black hole candidates.

We report results of VRI photometric monitoring of the T Tauri star plus disk system HH30. We find that HH30 is highly variable over timescales of a few days with Delta V~1.5 mag, Delta I~1.1 mag. Furthermore we find hints of periodicity with periodograms indicating possible periods of 11.6 and 19.8 days. The VRI photometry is available through the anonymous ftp service. We model the variability with Monte Carlo radiation transfer simulations for a spotted star plus disk system and find that the large variability is best reproduced with a single hot spot and circumstellar grains that are larger than typical interstellar grains. The apparent existence of a single hot spot and the need for larger, more forward throwing grains is consistent with previous modeling of HST imagery.

We have studied the UV extinction properties along 30 Galactic sightlines using data from the International Ultraviolet Explorer (IUE) archive that have never been previously examined. These distant (d > 1 kpc) sightlines were selected to investigate the distribution and physical conditions of gas located in low density regions of the Galactic disk and halo. The average densities along these sightlines are extremely low. It is likely that they are dominated by the warm intercloud medium and have little contribution from the cold cloud medium. We find that a subsample of these sightlines has extinction curves with weak bumps and very steep far-UV extinction reminiscent of the Magellanic clouds. These sightlines all lie in the region bounded by 325 < l < 0 and -5 > b > -11. The gas along these sightlines shows forbidden velocities which may indicate that the dust has been subject to shocks. This type of low density sightline may mimic the environments found in the Magellanic Clouds. Large values of N(Ca II)/N(Na I) indicating low depletion are associated with steep far-UV extinction. A possible correlation exists between decreasing bump strength and increasing far-UV steepness for extinction curves in the Galaxy and the Magellanic Clouds.

The USA Experiment is a new X-ray timing experiment with large collecting area and microsecond time resolution capable of conducting a broad program of studies of galactic X-ray binaries. USA is one of nine experiments aboard the Advanced Research and Global Observation Satellite which was launched February 23, 1999. USA is a collimated proportional counter X-ray telescope with about 1000 cm$^2$ of effective area per detector with two detectors sensitive to photons in the energy range 1-15 keV. A unique feature of USA is that photon events are time tagged by reference to an onboard GPS receiver allowing precise absolute time and location determination. We will present an overview of the USA instrument, capabilities, and scientific observing plan as well as the current status of the instrument.

We describe planetesimal accretion calculations for the dusty ring observed in the nearby A0 star HR 4796A. Models with initial masses of 10-20 times the minimum mass solar nebula produce a ring of width 7-15 AU and height 0.3-0.6 AU at 70 AU in roughly 10 Myr. The ring has a radial optical depth of 1. These results agree with limits derived from infrared images and from the excess infrared luminosity.

We present new images of the reflection nebula surrounding the R Coronae Borealis Star, UW Cen. This nebula, first detected in 1990, has changed its appearance significantly. At the estimated distance of UW Cen, this nebula is approximately 0.6 ly in radius so the nebula cannot have physically altered in only 8 years. Instead, the morphology of the nebula appears to change as different parts are illuminated by light from the central star modulated by shifting thick dust clouds near its surface. These dust clouds form and dissipate at irregular intervals causing the well-known declines in the R Coronae Borealis (RCB) stars. In this way, the central star acts like a lighthouse shining through holes in the dust clouds and lighting up different portions of the nebula. The existence of this nebula provides clues to the evolutionary history of RCB stars possibly linking them to the Planetary Nebulae and the final helium shell flash stars.