The Crab pulsar has striking radio emission properties, with the two dominant pulse components -- the main pulse and the interpulse -- consisting entirely of giant pulses. The emission is scattered in both the Crab nebula and the interstellar medium, causing multi-path propagation and thus scintillation. We study the scintillation of the Crab's giant pulses using phased Westerbork data at 1668\,MHz. We find that giant pulse spectra correlate at only $\sim 2 \%$, much lower than the $1/3$ correlation expected from a randomized signal imparted with the same impulse response function. In addition, we find that the main pulse and the interpulse appear to scintillate differently, with the scintillation of the interpulse correlating to the main pulse at a lower amplitude over a longer time and frequency range. These lines of evidence suggest that the giant pulse emission regions are extended, and that the main pulse and interpulse arise in physically distinct regions which are resolved by the scattering screen. Assuming the scattering takes place in the nebular filaments, the emission regions are of order a light cylinder radius, as projected on the sky. With further VLBI and multi-frequency data, it may be possible to measure the distance to the scattering screens, the size of giant pulse emission regions, and the physical separation between the pulse components.

Binaries harbouring millisecond pulsars enable a unique path to determine neutron star masses: radio pulsations reveal the motion of the neutron star, while that of the companion can be characterised through studies in the optical range. PSR J1012+5307 is a millisecond pulsar in a 14.5-h orbit with a helium-core white dwarf companion. In this work we present the analysis of an optical spectroscopic campaign, where the companion star absorption features reveal one of the lightest known white dwarfs. We determine a white dwarf radial velocity semi-amplitude of K_2 = 218.9 +- 2.2 km/s, which combined with that of the pulsar derived from the precise radio timing, yields a mass ratio of q=10.44+- 0.11. We also attempt to infer the white dwarf mass from observational constraints using new binary evolution models for extremely low-mass white dwarfs, but find that they cannot reproduce all observed parameters simultaneously. In particular, we cannot reconcile the radius predicted from binary evolution with the measurement from the photometric analysis (R_WD=0.047+-0.003 Rsun). Our limited understanding of extremely low-mass white dwarf evolution, which results from binary interaction, therefore comes as the main factor limiting the precision with which we can measure the mass of the white dwarf in this system. Our conservative white dwarf mass estimate of M_WD = 0.165 +- 0.015 Msun, along with the mass ratio enables us to infer a pulsar mass of M_NS = 1.72 +- 0.16 Msun. This value is clearly above the canonical 1.4 Msun, therefore adding PSR J1012+5307 to the growing list of massive millisecond pulsars.

The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale.

The delay time distribution of Type Ia supernovae (the time-dependent rate of supernovae resulting from a burst of star formation) has been measured using different techniques and in different environments. Here, we study in detail the distribution for field galaxies, using the SDSS DR7 Stripe 82 supernova sample. We improve a technique we introduced earlier, which is based on galaxy color and luminosity, and is insensitive to details of the star formation history, to include the normalization. Assuming a power-law dependence of the supernova rate with time, ${\rm DTD}(t)=A(t/{\rm 1\,Gyr})^{s}$, we find a power-law index $s= -1.34 ^{+0.19} _{-0.17}$ and a normalization $\rm{log}\ A= -12.15 ^{+0.10} _{-0.13}\, {\rm dex}(M_\odot ^{-1}\, \rm{yr}^{-1})$, corresponding to a number of type Ia supernovae integrated over a Hubble time of $PE = 0.004^{+0.002} _{-0.001}\, M_\odot ^{-1}$. We also implement a method used by Maoz and collaborators, which is based on star formation history reconstruction, and find that this gives a consistent result for the slope, but a lower, marginally inconsistent normalization. With our normalization, the distribution for field galaxies is made consistent with that derived for cluster galaxies. Comparing the inferred distribution with predictions from different evolutionary scenarios for type Ia supernovae, we find that our results are intermediate between the various predictions and do not yet constraint the evolutionary path leading to SNe~Ia.

How an otherwise inert carbon-oxygen white dwarf can be made to explode as a Type Ia supernova remains unknown. A promising test of theoretical models is to constrain the distribution of material that is left unburned, in particular of carbon. So far, most investigations used line identification codes to detect carbon in the ejecta, a method that cannot be readily compared against model predictions because it requires assumed opacities and temperatures. Here, we instead use tomographic techniques to investigate the amount of carbon in the inner layers of SN~2011fe, starting from the previously published tomographic analysis of Mazzali et al. From the presence of the carbon feature in the optical at early epochs and its disappearance later on, we derive an average carbon mass fraction between 0.001 and 0.05 for velocities in the range $13500 \lesssim v \lesssim 16000\ \rm{km\ s^{-1}}$, and an upper limit of 0.005 inside that region. Based on our models and the assumed density profile, only small amounts of carbon should be in the neutral state, too little to be responsible for features seen in near-infrared spectra that were previously identified as due to neutral carbon; We discuss possible reasons for this discrepancy. We compare our results against a suite of explosion models, although uncertainties in both the models and our simulations make it difficult to draw definitive conclusions.

Many radio pulsars have stable pulse profiles, but some exhibit mode changing where the profile switches between two or more quasi-stable modes of emission. So far, these effects had only been seen in relatively slow pulsars, but we show here that the pulse profile of PSR B1957+20, a millisecond pulsar, switches between two modes, with a typical time between mode changes of only $1.7$ s (or $\sim\!1000$ rotations), the shortest observed so far. The two modes differ in both intensity and polarization, with relatively large differences in the interpulse and much more modest ones in the main pulse. We find that the changes in the interpulse precede those in the main pulse by $\sim\!25$ ms, placing an empirical constraint on the timescale over which mode changes occurs. We also find that the properties of the giant pulses emitted by PSR B1957+20 are correlated with the mode of the regular emission: their rate and the rotational phase at which they are emitted both depend on mode. Furthermore, the energy distribution of the giant pulses emitted near the main pulse depends on mode as well. We discuss the ramifications for our understanding of the radio emission mechanisms as well as for pulsar timing experiments.

In regions with strongly varying electron density, radio emission can be magnified significantly by plasma lensing. In the presence of magnetic fields, magnification in time and frequency will be different for two circular polarizations. We show how these effects can be used to measure or constrain the magnetic field parallel to the line of sight, $B_\parallel$, as well as its spatial structure, $\sigma_{B_\parallel}$, in the lensing region. In addition, we discuss how generalized Faraday rotation can constrain the strength of the perpendicular field, $B_\perp$. We attempt to make such measurements for the Black Widow pulsar, PSR~B1957+20, in which plasma lensing was recently discovered. For this system, pressure equilibrium suggests $B\gtrsim 20\,$G at the interface between the pulsar and companion winds, where the radio eclipse starts and ends, and where most lensing occurs. We find no evidence for large-scale magnetic fields, with, on average, $B_\parallel=0.02\pm0.09\,$G over the egress lensing region. From individual lensing events, we strongly constrain small scale magnetic structure to $\sigma_B<10\,$mG, thus excluding scenarios with a strong but rapidly varying field. Finally, from the lack of reduction of average circular polarization in the same region, we rule out a strong, quasi-transverse field. We cannot identify any plausible scenario in which a large magnetic field in this system is concealed, leaving the nature of the interface between the pulsar and companion winds an enigma. Our method can be applied to other sources showing plasma lensing, including other eclipsing pulsars and fast radio bursts, to study the local properties of the magnetic field.

PSR J0636+5128 is a millisecond pulsar in one of the most compact pulsar binaries known, with a 96\,min orbital period. The pulsar mass function suggests a very low-mass companion, similar to that seen in so-called "black widow" binaries. Unlike in most of those, however, no radio eclipses by material driven off from the companion were seen leading to the possibility that the companion was a degenerate remnant of a carbon-oxygen white dwarf. We report the discovery of the optical counterpart of its companion in images taken with the Gemini North and Keck~I telescopes. The companion varies between $r=25$ and $r=23$ on the 96\,min orbital period of the binary, caused by irradiation from the pulsar's energetic wind. We modeled the multi-color lightcurve using parallax constraints from pulsar timing and determine a companion mass of $(1.71\pm0.23)\times 10^{-2}\,M_\odot$, a radius of $(7.6\pm1.4)\times 10^{-2}\,R_\odot$, and a mean density of $54\pm26\,{\rm g\,cm}^{-3}$, all for an assumed neutron star mass of $1.4\,M_\odot$. This makes the companion to PSR J0636+5128 one of the densest of the "black widow" systems. Modeling suggests that the composition is not predominantly hydrogen, perhaps due to an origin in an ultra-compact X-ray binary.

Radio pulsars scintillate because their emission travels through the ionized interstellar medium via multiple paths, which interfere with each other. It has long been realized that the scattering screens responsible for the scintillation could be used as `interstellar lenses' to localize pulsar emission regions. Most scattering screens, however, only marginally resolve emission components, limiting results to statistical inferences and detections of small positional shifts. Since screens situated close to the source have better resolution, it should be easier to resolve emission regions of pulsars located in high density environments such as supernova remnants or binaries in which the pulsar's companion has an ionized outflow. Here, we report events of extreme plasma lensing in the `Black Widow' pulsar, PSR~B1957+20, near the phase in its 9.2 hour orbit in which its emission is eclipsed by its companion's outflow. During the lensing events, the flux is enhanced by factors of up to 70--80 at specific frequencies. The strongest events clearly resolve the emission regions: they affect the narrow main pulse and parts of the wider interpulse differently. We show that the events arise naturally from density fluctuations in the outer regions of the outflow, and infer a resolution of our lenses comparable to the pulsar's radius, about 10\,km. Furthermore, the distinct frequency structures imparted by the lensing are reminiscent of what is observed for the repeating fast radio burst FRB 121102, providing observational support for the idea that this source is observed through, and thus at times strongly magnified by, plasma lenses.

We present near infrared (NIR) spectroscopy of the nearby supernova 2014J obtained $\sim$450 d after explosion. We detect the [Ni II] 1.939 $\mu$m line in the spectra indicating the presence of stable $^{58}$Ni in the ejecta. The stable nickel is not centrally concentrated but rather distributed as the iron. The spectra are dominated by forbidden [Fe II] and [Co II] lines. We use lines, in the NIR spectra, arising from the same upper energy levels to place constraints on the extinction from host galaxy dust. We find that that our data are in agreement with the high $A_V$ and low $R_V$ found in earlier studies from data near maximum light. Using a $^{56}$Ni mass prior from near maximum light $\gamma$-ray observations, we find $\sim$0.05 M$_\odot$ of stable nickel to be present in the ejecta. We find that the iron group features are redshifted from the host galaxy rest frame by $\sim$600 km s$^{-1}$.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

Type Ia supernovae are generally agreed to arise from thermonuclear explosions of carbon-oxygen white dwarfs. The actual path to explosion, however, remains elusive, with numerous plausible parent systems and explosion mechanisms suggested. Observationally, type Ia supernovae have multiple subclasses, distinguished by their lightcurves and spectra. This raises the question whether these reflect that multiple mechanisms occur in nature, or instead that explosions have a large but continuous range of physical properties. We revisit the idea that normal and 91bg-like supernovae can be understood as part of a spectral sequence, in which changes in temperature dominate. Specifically, we find that a single ejecta structure is sufficient to provide reasonable fits of both the normal type Ia supernova SN~2011fe and the 91bg-like SN~2005bl, provided that the luminosity and thus temperature of the ejecta are adjusted appropriately. This suggests that the outer layers of the ejecta are similar, thus providing some support of a common explosion mechanism. Our spectral sequence also helps to shed light on the conditions under which carbon can be detected in pre-maximum SN~Ia spectra -- we find that emission from iron can "fill in" the carbon trough in cool SN~Ia. This may indicate that the outer layers of the ejecta of events in which carbon is detected are relatively metal poor compared to events where carbon is not detected.

The interstellar medium scatters radio waves which causes pulsars to scintillate. For intrinsically short bursts of emission, the observed signal should be a direct measurement of the impulse response function. We show that this is indeed the case for giant pulses from PSR B1957+20: from baseband observations at 327 MHz, we demonstrate that the observed voltages of a bright pulse allow one to coherently descatter nearby ones. We find that while the scattering timescale is $12.3\,\mu$s, the power in the descattered pulses is concentrated within a span almost two orders of magnitude shorter, of $\lesssim\!200\,$ns. This sets an upper limit to the intrinsic duration of the giant pulses. We verify that the response inferred from the giant pulses is consistent with the scintillation pattern obtained by folding the regular pulsed emission, and that it decorrelates on the same timescale, of~$84\,$s. In principle, with large sets of giant pulses, it should be possible to constrain the structure of the scattering screen much more directly than with other current techniques, such as holography on the dynamic spectrum and cyclic spectroscopy.

We observed a nearby millisecond pulsar J2124-3358 with the Hubble Space Telescope in broad far-UV (FUV) and optical filters. The pulsar is detected in both bands with fluxes F(1250-2000 A)= (2.5+/-0.3)x10^-16 erg/s/cm^2 and F(3800-6000 A)=(6.4+/-0.4)x10^-17 erg/s/cm^2, which correspond to luminosities of ~5.8x10^27 and 1.4x10^27 erg/s, for d=410 pc and E(B-V)=0.03. The optical-FUV spectrum can be described by a power-law model, f_nu~nu^alpha, with slope alpha=0.18-0.48 for a conservative range of color excess, E(B-V)=0.01-0.08. Since a spectral flux rising with frequency is unusual for pulsar magnetospheric emission in this frequency range, it is possible that the spectrum is predominantly magnetospheric (power law with alpha<0) in the optical while it is dominated by thermal emission from the neutron star surface in the FUV. For a neutron star radius of 12 km, the surface temperature would be between 0.5x10^5 and 2.1x10^5 K, for alpha ranging from -1 to 0, E(B-V)=0.01-0.08, and d=340-500 pc. In addition to the pulsar, the FUV images reveal extended emission spatially coincident with the known Halpha bow shock, making PSR J2124-3358 the second pulsar (after PSR J0437-4715) with a bow shock detected in FUV.

White dwarf (WD) binary mergers are possible progenitors to a number of unusual stars and transient phenomena, including type Ia supernovae. To date, simulations of mergers have not included magnetic fields, even though they are believed to play a significant role in the evolution of the merger remnant. We simulated a 0.625 - 0.65 $M_{\odot}$ carbon-oxygen WD binary merger in the magnetohydrodynamic moving mesh code Arepo. Each WD was given an initial dipole field with a surface value of $\sim10^3$ G. As in simulations of merging double neutron star binaries, we find exponential field growth within Kelvin-Helmholtz instability-generated vortices during the coalescence of the two stars. The final field has complex geometry, and a strength $>10^{10}$ G at the center of the merger remnant. Its energy is $\sim2\times10^{47}$ ergs, $\sim0.2$% of the remnant's total energy. The strong field likely influences further evolution of the merger remnant by providing a mechanism for angular momentum transfer and additional heating, potentially helping to ignite carbon fusion.

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

The hierarchical triple system PSR J0337+1715 offers an unprecedented laboratory to study secular evolution of interacting systems and to explore the complicated mass-transfer history that forms millisecond pulsars and helium-core white dwarfs. The latter in particular, however, requires knowledge of the properties of the individual components of the system. Here we present precise optical spectroscopy of the inner companion in the PSR J0337+1715 system. We confirm it as a hot, low-gravity DA white dwarf with Teff=15,800+/-100 K and log(g)=5.82+/-0.05. We also measure an inner mass ratio of 0.1364+/-0.0015, entirely consistent with that inferred from pulsar timing, and a systemic radial velocity of 29.7+/-0.3 km/s. Combined with the mass (0.19751 Msun) determined from pulsar timing, our measurement of the surface gravity implies a radius of 0.091+/-0.005 Rsun; combined further with the effective temperature and extinction, the photometry implies a distance of 1300+/-80 pc. The high temperature of the companion is somewhat puzzling: with current models, it likely requires a recent period of unstable hydrogen burning, and suggests a surprisingly short lifetime for objects at this phase in their evolution. We discuss the implications of these measurements in the context of understanding the PSR J0337+1715 system, as well as of low-mass white dwarfs in general.

Gravitationally bound three-body systems have been studied for hundreds of years and are common in our Galaxy. They show complex orbital interactions, which can constrain the compositions, masses, and interior structures of the bodies and test theories of gravity, if sufficiently precise measurements are available. A triple system containing a radio pulsar could provide such measurements, but the only previously known such system, B1620-26 (with a millisecond pulsar, a white dwarf, and a planetary-mass object in an orbit of several decades), shows only weak interactions. Here we report precision timing and multi-wavelength observations of PSR J0337+1715, a millisecond pulsar in a hierarchical triple system with two other stars. Strong gravitational interactions are apparent and provide the masses of the pulsar (1.4378(13) Msun, where Msun is the solar mass and the parentheses contain the uncertainty in the final decimal places) and the two white dwarf companions (0.19751(15) Msun and 0.4101(3) Msun), as well as the inclinations of the orbits (both approximately 39.2 degrees). The unexpectedly coplanar and nearly circular orbits indicate a complex and exotic evolutionary past that differs from those of known stellar systems. The gravitational field of the outer white dwarf strongly accelerates the inner binary containing the neutron star, and the system will thus provide an ideal laboratory in which to test the strong equivalence principle of general relativity.

The merger of two carbon-oxygen white dwarfs can lead either to a spectacular transient, stable nuclear burning or a massive, rapidly rotating white dwarf. Simulations of mergers have shown that the outcome strongly depends on whether the white dwarfs are similar or dissimilar in mass. In the similar-mass case, both white dwarfs merge fully and the remnant is hot throughout, while in the dissimilar case, the more massive, denser white dwarf remains cold and essentially intact, with the disrupted lower mass one wrapped around it in a hot envelope and disk. In order to determine what constitutes "similar in mass" and more generally how the properties of the merger remnant depend on the input masses, we simulated unsynchronized carbon-oxygen white dwarf mergers for a large range of masses using smoothed-particle hydrodynamics. We find that the structure of the merger remnant varies smoothly as a function of the ratio of the central densities of the two white dwarfs. A density ratio of 0.6 approximately separates similar and dissimilar mass mergers. Confirming previous work, we find that the temperatures of most merger remnants are not high enough to immediately ignite carbon fusion. During subsequent viscous evolution, however, the interior will likely be compressed and heated as the disk accretes and the remnant spins down. We find from simple estimates that this evolution can lead to ignition for many remnants. For similar-mass mergers, this would likely occur under sufficiently degenerate conditions that a thermonuclear runaway would ensue.

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

Most millisecond pulsars with low-mass companions are in systems with either helium-core white dwarfs or non-degenerate ("black widow" or "redback") stars. A candidate counterpart to PSR J1816+4510 was identified by Kaplan et al. (2012) whose properties were suggestive of both types of companions although identical to neither. We have assembled optical spectroscopy of the candidate companion and confirm that it is part of the binary system with a radial velocity amplitude of 343+/-7 km/s, implying a high pulsar mass, Mpsr*sin^3i=1.84+/-0.11 Msun, and a companion mass Mc*sin^3i=0.192+/-0.012 Msun, where i is the inclination of the orbit. The companion appears similar to proto-white dwarfs/sdB stars, with a gravity log(g)=4.9+/-0.3, and effective temperature 16000+/-500 K. The strongest lines in the spectrum are from hydrogen, but numerous lines from helium, calcium, silicon, and magnesium are present as well, with implied abundances of roughly ten times solar (relative to hydrogen). As such, while from the spectrum the companion to PSR J1816+4510 is superficially most similar to a low-mass white dwarf, it has much lower gravity, is substantially larger, and shows substantial metals. Furthermore, it is able to produce ionized gas eclipses, which had previously been seen only for low-mass, non-degenerate companions in redback or black widow systems. We discuss the companion in relation to other sources, but find we understand neither its nature nor its origins. Thus, the system is interesting for understanding unusual stellar products of binary evolution, as well as, independent of its nature, for determining neutron-star masses.

In the last few years, over 43 millisecond radio pulsars have been discovered by targeted searches of unidentified gamma-ray sources found by the Fermi Gamma-Ray Space Telescope. A large fraction of these millisecond pulsars are in compact binaries with low-mass companions. These systems often show eclipses of the pulsar signal and are commonly known as black widows and redbacks because the pulsar is gradually destroying its companion. In this paper, we report on the optical discovery of four strongly irradiated millisecond pulsar companions. All four sources show modulations of their color and luminosity at the known orbital periods from radio timing. Light curve modelling of our exploratory data shows that the equilibrium temperature reached on the companion's dayside with respect to their nightside is consistent with about 10-30% of the available spin-down energy from the pulsar being reprocessed to increase the companion's dayside temperature. This value compares well with the range observed in other irradiated pulsar binaries and offers insights about the energetics of the pulsar wind and the production of gamma-ray emission. In addition, this provides a simple way of estimating the brightness of irradiated pulsar companions given the pulsar spin-down luminosity. Our analysis also suggests that two of the four new irradiated pulsar companions are only partially filling their Roche lobe. Some of these sources are relatively bright and represent good targets for spectroscopic follow-up. These measurements could enable, among other things, mass determination of the neutron stars in these systems.

The nature of the progenitors of type Ia supernovae is still a mystery. While plausible candidates are known for both the single degenerate and double degenerate models, the observed numbers fall significantly short of what is required to reproduce the type Ia supernovae rate. Some of the most promising single-degenerate type Ia progenitors are recurrent novae and super-soft sources (SSS). White dwarfs with higher mass transfer rates can also be type Ia supernova progenitors. For these rapidly accreting white dwarfs (RAWD), more material than is needed for steady burning accretes on the white dwarf, and extends the white dwarf's photosphere. Unlike super-soft sources, such objects will likely not be detectable at soft X-ray energies, but will be bright at longer wavelengths, such as the far ultraviolet (UV). Possible examples include LMC N66 and the V Sagittae stars. We present a survey using multi-object spectrographs looking for RAWD in the central core of the SMC, from objects selected to be bright in the far UV and with blue far UV-V colors. While we find some unusual objects, and recover known planetary nebula and WR stars, we detect no candidate RAWD. The upper limits from this non-detection depending on our expectations of what a RAWD should look like, as well assumptions about the internal extinction of the SMC. Assuming they resemble LMC N66 or fainter versions of Wolf-Rayet stars we set an upper limit of 10-14 RAWD in the SMC. However, our survey is unlikely to detect objects like V Sge, and hence we cannot set meaningful upper limits if RAWD generally resemble these.

SDSS 1355+0856 was identified as a hot white dwarf (WD) with a binary companion from time-resolved SDSS spectroscopy as part of the ongoing SWARMS survey. Follow-up observations with the ARC 3.5m telescope and the MMT revealed weak emission lines in the central cores of the Balmer absorption lines during some phases of the orbit, but no line emission during other phases. This can be explained if SDSS 1355+0856 is a detached WD+M dwarf binary similar to GD 448, where one of the hemispheres of the low-mass companion is irradiated by the proximity of the hot white dwarf. Based on the available data, we derive a period of 0.11438 +- 0.00006 days, a primary mass of 0.46 +- 0.01 solar masses, a secondary mass between 0.083 and 0.097 solar masses, and an inclination larger than 57 degrees. This makes SDSS 1355+0856 one of the shortest period post-common envelope WD+M dwarf binaries known, and one of only a few where the primary is likely a He-core white dwarf, which has interesting implications for our understanding of common envelope evolution and the phenomenology of cataclysmic variables. The short cooling time of the WD (25 Myr) implies that the system emerged from the common envelope phase with a period very similar to what we observe today, and was born in the period gap of cataclysmic variables.

Thermonuclear supernovae result when interaction with a companion reignites nuclear fusion in a carbon-oxygen white dwarf, causing a thermonuclear runaway, a catastrophic gain in pressure, and the disintegration of the whole white dwarf. It is usually thought that fusion is reignited in near-pycnonuclear conditions when the white dwarf approaches the Chandrasekhar mass. I briefly describe two long-standing problems faced by this scenario, and our suggestion that these supernovae instead result from mergers of carbon-oxygen white dwarfs, including those that produce sub-Chandrasekhar mass remnants. I then turn to possible observational tests, in particular those that test the absence or presence of electron captures during the burning.

PSR J1738+0333 is one of the four millisecond pulsars known to be orbited by a white dwarf companion bright enough for optical spectroscopy. Of these, it has the shortest orbital period, making it especially interesting for a range of astrophysical and gravity related questions. We present a spectroscopic and photometric study of the white dwarf companion and infer its radial velocity curve, effective temperature, surface gravity and luminosity. We find that the white dwarf has properties consistent with those of low-mass white dwarfs with thick hydrogen envelopes, and use the corresponding mass-radius relation to infer its mass; M_WD = 0.181 +/- +0.007/-0.005 solar masses. Combined with the mass ratio q=8.1 +/- 0.2 inferred from the radial velocities and the precise pulsar timing ephemeris, the neutron star mass is constrained to M_PSR = 1.47 +/- +0.07/-0.06 solar masses. Contrary to expectations, the latter is only slightly above the Chandrasekhar limit. We find that, even if the birth mass of the neutron star was only 1.20 solar masses, more than 60% of the matter that left the surface of the white dwarf progenitor escaped the system. The accurate determination of the component masses transforms this system in a laboratory for fundamental physics by constraining the orbital decay predicted by general relativity. Currently, the agreement is within 1 sigma of the observed decay. Further radio timing observations will allow precise tests of white dwarf models, assuming the validity of general relativity.

We present optical spectroscopic measurements of the eclipsing High Mass X-ray Binary XMMU J013236.7+303228 in M33. Based on spectra taken at multiple epochs of the 1.73d binary orbital period we determine physical as well as orbital parameters for the donor star. We find the donor to be a B1.5IV sub-giant with effective temperature T=22,000-23,000 K. From the luminosity, temperature and known distance to M33 we derive a radius of R = 8.9 \pm 0.5 R_sun. From the radial--velocity measurements, we determine a velocity semi-amplitude of K_opt = 63 \pm 12 km/sec. Using the physical properties of the B-star determined from the optical spectrum, we estimate the star's mass to be M_opt = 11 \pm 1 M_sun. Based on the X-ray spectrum, the compact companion is likely a neutron star, although no pulsations have yet been detected. Using the spectroscopically derived B-star mass we find the neutron star companion mass to be M_X = 2.0 \pm 0.4 M_sun, consistent with the neutron star mass in the HMXB Vela X-1, but heavier than the canonical value of 1.4 M_sun found for many millisecond pulsars. We attempt to use as an additional constraint that the B star radius inferred from temperature, flux, and distance, should equate the Roche radius, since the system accretes by Roche lobe overflow. This leads to substantially larger masses, but from trying to apply the technique to known systems, we find that the masses are consistently overestimated. Attempting to account for that in our uncertainties, we derive M_X = 2.2^{+0.8}_{-0.6} M_sun and M_opt =13 \pm 4 M_sun. We conclude that precise constraints require detailed modeling of the shape of the Roche surface.

We present the results of a multiplicity survey of 212 T Tauri stars in the Chamaeleon I and Taurus-Auriga star-forming regions, based on high-resolution spectra from the Magellan Clay 6.5 m telescope. From these data, we achieved a typical radial velocity precision of ~80 m/s with slower rotators yielding better precision, in general. For 174 of these stars, we obtained multi-epoch data with sufficient time baselines to identify binaries based on radial velocity variations. We identified eight close binaries and four close triples, of which three and two, respectively, are new discoveries. The spectroscopic multiplicity fractions we find for Cha I (7%) and Tau-Aur (6%) are similar to each other, and to the results of field star surveys in the same mass and period regime. However, unlike the results from imaging surveys, the frequency of systems with close companions in our sample is not seen to depend on primary mass. Additionally, we do not find a strong correlation between accretion and close multiplicity. This implies that close companions are not likely the main source of the accretion shut down observed in weak-lined T Tauri stars. Our results also suggest that sufficient radial velocity precision can be achieved for at least a subset of slowly rotating young stars to search for hot Jupiter planets.

Middle-aged, cooling neutron stars are observed both as relatively rapidly spinning radio pulsars and as more slowly spinning, strongly magnetized isolated neutron stars (INSs), which stand out by their thermal X-ray spectra. The difference between the two classes may be that the INSs initially had much stronger magnetic fields, which decayed. To test this, we used the Chandra X-ray Observatory to observe 1RXS J072559.8-261229, a possible X-ray counterpart to PSR J0726-2612, which, with its 3.44s period and 3e13G inferred magnetic field strength, is the nearest and least extincted among the possible slowly-spinning, strong-field INS progenitors (it likely is in the Gould Belt, at ~1 kpc). We confirm the identification and find that the pulsar has a spectrum consistent with being purely thermal, with blackbody temperature kT=87+/-5 eV and radius R=5.7+2.6-1.3 km at a distance of 1 kpc. We detect sinusoidal pulsations at twice the radio period with a semi-amplitude of 27\pm5%. The properties of PSR J0726-2612 strongly resemble those of the INSs, except for its much shorter characteristic age of 200 kyr (instead of several Myr). We conclude that PSR J0726-2612 is indeed an example of a young INS, one that started with a magnetic field strength on the low end of those inferred for the INSs, and that, therefore, decayed by a relatively small amount. Our results suggest that the long-period, strong-field pulsars and the INSs are members of the same class, and open up new opportunities to understand the puzzling X-ray and optical emission of the INSs through radio observations of PSR J0726-2612.

We present extensive spectral and photometric observations of the recycled pulsar/white-dwarf binary containing PSR J0437-4715, which we analyzed together with archival X-ray and gamma-ray data, to obtain the complete mid-infrared to gamma-ray spectrum. We first fit each part of the spectrum separately, and then the whole multi-wavelength spectrum. We find that the optical-infrared part of the spectrum is well fit by a cool white dwarf atmosphere model with pure hydrogen composition. The model atmosphere (Teff = 3950pm150K, log g=6.98pm0.15, R_WD=(1.9pm0.2)e9 cm) fits our spectral data remarkably well for the known mass and distance (M=0.25pm0.02Msun, d=156.3pm1.3pc), yielding the white dwarf age (tau=6.0pm0.5Gyr). In the UV, we find a spectral shape consistent with thermal emission from the bulk of the neutron star surface, with surface temperature between 1.25e5 and 3.5e5K. The temperature of the thermal spectrum suggests that some heating mechanism operates throughout the life of the neutron star. The temperature distribution on the neutron star surface is non-uniform. In the X-rays, we confirm the presence of a high-energy tail which is consistent with a continuation of the cut-off power-law component (Gamma=1.56pm0.01, Ecut=1.1pm0.2GeV) that is seen in gamma-rays and perhaps even extends to the near-UV.

We present an analysis and interpretation of the Kepler binary system KOI 1224. This is the fourth binary found with Kepler that consists of a thermally bloated, hot white dwarf in a close orbit with a more or less normal star of spectral class A or F. As we show, KOI 1224 contains a white dwarf with Teff = 14400 +/- 1100 K, mass = 0.20 +/- 0.02 Msun, and radius = 0.103 +/- 0.004 Rsun, and an F-star companion of mass = 1.59 +/- 0.07 Msun that is somewhat beyond its terminal-age main sequence. The orbital period is quite short at 2.69802 days. The ingredients that are used in the analysis are the Kepler binary light curve, including the detection of the Doppler boosting effect; the NUV and FUV fluxes from the Galex images of this object; an estimate of the spectral type of the F-star companion; and evolutionary models of the companion designed to match its effective temperature and mean density. The light curve is modelled with a new code named Icarus which we describe in detail. Its features include the full treatment of orbital phase-resolved spectroscopy, Doppler boosting, irradiation effects and transits/eclipses, which are particularly suited to irradiated eclipsing binaries. We interpret the KOI 1224 system in terms of its likely evolutionary history. We infer that this type of system, containing a bloated hot white dwarf, is the direct descendant of an Algol-type binary. In spite of this basic understanding of the origin of KOI 1224, we discuss a number of problems associated with producing this type of system with this short of an short orbital period.

Missing from the usual considerations of nuclear burning white dwarfs as Type Ia supernova progenitors are systems with very higher mass transfer rates, where more material than is needed for steady burning accretes on the white dwarf. This will move out the photosphere of the white dwarf, causing it to emit at longer wavelengths. Thus, we propose the name ultra-soft source (USS) for these objects. We present a VLT/FLAMES survey looking for USS in the SMC, selected to be bright in the far UV and with blue far UV-V colors. While we find some unusual objects, and recover known planetary nebula and WR stars, we detect no objects with strong He II lines, which should be a signature of USS. This null result either puts an upper limit on the number of USS in the SMC, or shows that we do not understand what the optical spectra of such objects will look like. We also discuss the unusual LMC [WN] planetary nebula LMC N66 as a possible example of an USS. It has a luminosity consistent with that expected, and its spectra show incompletely CNO-processed material - strong helium lines, some hydrogen, enhanced nitrogen and depleted carbon. It also shows periodic outbursts. USS may resemble N66 in quiescence. However, it lacks a FUV excess, contrary to our predictions.

Recent studies have shown that for suitable initial conditions both super- and sub-Chandrasekhar mass carbon-oxygen white dwarf mergers produce explosions similar to observed SNe Ia. The question remains, however, how much fine tuning is necessary to produce these conditions. We performed a large set of SPH merger simulations, sweeping the possible parameter space. We find trends for merger remnant properties, and discuss how our results affect the viability of our recently proposed sub-Chandrasekhar merger channel for SNe Ia.

We present a phase-coherent timing solution for RX J0420.0-5022, the coolest (kT=45 eV) and fastest-spinning (P=3.45 s) of the seven so-called isolated neutron stars (INSs). Using 14 observations with the XMM-Newton spacecraft in 2010-2011, we were able to measure a spin-down rate nudot=(-2.3+/-0.2)e-15 Hz/s (Pdot=(2.8+/-0.3)e-14 s/s), from which we infer a dipolar magnetic field of 1.0e13 G. With reasonable confidence we were able to extend the timing solution back to archival XMM-Newton from 2002 and 2003, giving the same solution but with considerably more precision. This gives RX J0420 the lowest dipole magnetic field of the INSs. Our spectroscopy does not confirm the broad absorption feature at 0.3 keV hinted at in earlier observations, although difficulties in background subtraction near that energy make conclusions difficult. With this, all 6 of the INSs that have confirmed periodicities now have constrained spin-downs from coherent solutions. The evidence that the INSs are qualitatively different from rotation-powered pulsars now appears robust, with the likely conclusion that their characteristic ages are systematically older than their true ages, because their fields have decayed. The field decay probably also causes them to be systematically hotter than pulsars of the same (true) age.

We present Hubble Space Telescope optical and ultraviolet photometry for five nearby, thermally emitting neutron stars. With these measurements, all seven such objects have confirmed optical and ultraviolet counterparts. Combining our data with archival space-based photometry, we present spectral energy distributions for all sources and measure the "optical excess": the factor by which the measured photometry exceeds that extrapolated from X-ray spectra. We find that the majority have optical and ultraviolet fluxes that are inconsistent with that expected from thermal (Rayleigh-Jeans) emission, exhibiting more flux at longer wavelengths. We also find that most objects have optical excesses between 5 and 12, but that one object (RX J2143.0+0654) exceeds the X-ray extrapolation by a factor of more than 50 at 5000 A, and that this is robust to uncertainties in the X-ray spectra and absorption. We consider explanations for this ranging from atmospheric effects, magnetospheric emission, and resonant scattering, but find that none is satisfactory.

The most massive neutron stars constrain the behavior of ultra-dense matter, with larger masses possible only for increasingly stiff equations of state. Here, we present evidence that the black widow pulsar, PSR B1957+20, has a high mass. We took spectra of its strongly irradiated companion and found an observed radial-velocity amplitude of K_obs=324+/-3 km/s. Correcting this for the fact that, due to the irradiation, the center of light lies inward relative to the center of mass, we infer a true radial-velocity amplitude of K_2=353+/-4 km/s and a mass ratio q=M_PSR/M_2=69.2+/-0.8. Combined with the inclination i=65+/-2 deg inferred from models of the lightcurve, our best-fit pulsar mass is M_PSR=2.40+/-0.12 M_sun. We discuss possible systematic uncertainties, in particular in the lightcurve modeling. Taking an upper limit of i<85 deg based on the absence of radio eclipses at high frequency, combined with a conservative lower-limit to the motion of the center of mass, K_2>343 km/s (q>67.3), we infer a lower limit to the pulsar mass of M_PSR>1.66 M_sun.

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

Type Ia supernovae are generally thought to be due to the thermonuclear explosions of carbon-oxygen white dwarfs with masses near the Chandrasekhar mass. This scenario, however, has two long-standing problems. First, the explosions do not naturally produce the correct mix of elements, but have to be finely tuned to proceed from sub-sonic deflagration to super-sonic detonation. Second, population models and observations give formation rates of near-Chandrasekhar white dwarfs that are far too small. Here, we suggest that type Ia supernovae instead result from mergers of roughly equal-mass carbon-oxygen white dwarfs, including those that produce sub-Chandrasekhar mass remnants. Numerical studies of such mergers have shown that the remnants consist of rapidly rotating cores that contain most of the mass and are hottest in the center, surrounded by dense, small disks. We argue that the disks accrete quickly, and that the resulting compressional heating likely leads to central carbon ignition. This ignition occurs at densities for which pure detonations lead to events similar to type Ia supernovae. With this merger scenario, we can understand the type Ia rates, and have plausible reasons for the observed range in luminosity and for the bias of more luminous supernovae towards younger populations. We speculate that explosions of white dwarfs slowly brought to the Chandrasekhar limit---which should also occur---are responsible for some of the "atypical" type Ia supernovae.

We report the discovery of 1RXS J173006.4+033813, a polar cataclysmic variable with a period of 120.21 min. The white dwarf primary has a magnetic field of B = 42+6-5 MG, and the secondary is a M3 dwarf. The system shows highly symmetric double peaked photometric modulation in the active state as well as in quiescence. These arise from a combination of cyclotron beaming and ellipsoidal modulation. The projected orbital velocity of the secondary is K2 = 390+-4 km/s. We place an upper limit of 830+-65 pc on the distance.

Giant planets are usually thought to form within a few tens of AU of their host stars, and hence it came as a surprise when we found what appeared to be a planetary mass (~0.008 Msun) companion around the 5 Myr-old solar mass star 1RXS J160929.1-210524 in the Upper Scorpius association. At the time, we took the object's membership in Upper Scorpius -- established from near-infrared, H- and K-band spectroscopy -- and its proximity (2.2", or 330 AU) to the primary as strong evidence for companionship, but could not verify their common proper motion. Here, we present follow-up astrometric measurements that confirm that the companion is indeed co-moving with the primary star, which we interpret as evidence that it is a truly bound planetary mass companion. We also present new J-band spectroscopy and 3.0-3.8 microns photometry of the companion. Based on a comparison with model spectra, these new measurements are consistent with the previous estimate of the companion effective temperature of 1800+/-200 K. We present a new estimate of the companion mass based on evolution models and the calculated bolometric luminosity of the companion; we obtain a value of 0.008 (-0.002/+0.003) Msun, again consistent with our previous result. Finally, we present angular differential imaging observations of the system allowing us to rule out additional planets in the system more massive than 1, 2 and 8 Mjup at projected separations larger than 3" (~440 AU), 0.7" (~100 AU) and 0.35" (~50 AU), respectively. This companion is the least massive known to date at such a large orbital distance; it shows that objects in the planetary mass range exist at orbital separations of several hundred AU, posing a serious challenge for current formation models.

Among the initial results from Kepler were two striking lightcurves, for KOI 74 and KOI 81, in which the relative depths of the primary and secondary eclipses showed that the more compact, less luminous object was hotter than its stellar host. That result became particularly intriguing because a substellar mass had been derived for the secondary in KOI 74, which would make the high temperature challenging to explain; in KOI 81, the mass range for the companion was also reported to be consistent with a substellar object. We re-analyze the Kepler data and demonstrate that both companions are likely to be white dwarfs. We also find that the photometric data for KOI 74 show a modulation in brightness as the more luminous star orbits, due to Doppler boosting. The magnitude of the effect is sufficiently large that we can use it to infer a radial velocity amplitude accurate to 1 km/s. As far as we are aware, this is the first time a radial-velocity curve has been measured photometrically. Combining our velocity amplitude with the inclination and primary mass derived from the eclipses and primary spectral type, we infer a secondary mass of 0.22+/-0.03 Msun. We use our estimates to consider the likely evolutionary paths and mass-transfer episodes of these binary systems.

SDSS 1257+5428 is a white dwarf in a close orbit with a companion that has been suggested to be a neutron star. If so, it hosts the closest known neutron star, and its existence implies a great abundance of similar systems and a rate of white-dwarf neutron-star mergers similar to that of the type Ia supernova rate. Here, we present high signal-to-noise spectra of SDSS 1257+5428, which confirm an independent finding that the system is in fact composed of two white dwarfs, one relatively cool and with low mass, and the other hotter and more massive. With this, the demographics and merger rate are no longer puzzling (various factors combine to lower the latter by more than two orders of magnitude). We show that the spectra are fit well with a combination of two hydrogen model atmospheres, as long as the lines of the higher-gravity component are broadened significantly relative to what is expected from just pressure broadening. Interpreting this additional broadening as due to rotation, the inferred spin period is short, about 1 minute. Similarly rapid rotation is only seen in accreting white dwarfs that are magnetic; empirically, it appears that in non-magnetized white dwarfs, accreted angular momentum is lost by nova explosions before it can be transferred to the white dwarf. This suggests that the massive white dwarf in SDSS 1257+5428 is magnetic as well, with B~10^5 G. Alternatively, the broadening seen in the spectral lines could be due to a stronger magnetic field, of ~10^6 G. The two models could be distinguished by further observations.

We analyze the X-ray scattering halos around three Galactic Anomalous X-ray Pulsars (AXPs) in order to constrain the distance and the optical extinction of each source. We obtain surface brightness distributions from EPIC-pn data obtained with XMM-Newton, compare the profiles of different sources, and fit them with a model based on the standard theory of X-ray scattering by dust grains, both for a uniform distribution of dust along the line of sight, and for dust distributions constrained by previous measurements. Somewhat surprisingly, we find that for all three sources, the uniform distribution reproduces the observed surface brightness as well as or better than the distributions that are informed by previous constraints. Nevertheless, the inferred total dust columns are robust, and serve to confirm that previous measurements based on interstellar edges in high-resolution X-ray spectra and on modelling of broad-band X-ray spectra were reliable. Specifically, we find Av ~= 4, 6, and 8 mag for 4U 0142+61, 1E 1048.1-5937, and 1RXS J170849.0-400910, respectively. For 1E 1048.1-5937, this is well in excess of the extinction expected towards a HI bubble along the line of sight, thus casting further doubt on the suggested association with the source.

The nearby isolated neutron stars are a group of seven relatively slowly rotating neutron stars that show thermal X-ray spectra, most with broad absorption features. They are interesting both because they may allow one to determine fundamental neutron-star properties by modeling their spectra, and because they appear to be a large fraction of the overall neutron-star population. Here, we describe a series of XMM-Newton observations of the nearby isolated neutron star RX J0806.4-4123, taken as part of larger program of timing studies. From these, we limit the spin-down rate to dnu/dt=(-4.3+/-2.3)*10^{-16} Hz/s. This constrains the dipole magnetic field to be <3.7e13 G at 2sigma, significantly less than the field of 1e14 G implied by simple models for the X-ray absorption found at 0.45 keV. We confirm that the spectrum is thermal and stable (to within a few percent), but find that the 0.45 keV absorption feature is broader and more complex than previously thought. Considering the population of isolated neutron stars, we find that magnetic field decay from an initial field of 3e14 G accounts most naturally for their timing and spectral properties, both qualitatively and in the context of the models for field decay of Pons and collaborators.

The cores of neutron stars harbor the highest matter densities known to occur in nature, up to several times the densities in atomic nuclei. Similarly, magnetic field strengths can exceed the strongest fields generated in terrestrial laboratories by ten orders of magnitude. Hyperon-dominated matter, deconfined quark matter, superfluidity, even superconductivity are predicted in neutron stars. Similarly, quantum electrodynamics predicts that in strong magnetic fields the vacuum becomes birefringent. The properties of matter under such conditions is governed by Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED), and the close study of the properties of neutron stars offers the unique opportunity to test and explore the richness of QCD and QED in a regime that is utterly beyond the reach of terrestrial experiments. Experimentally, this is almost virgin territory.

We analyze the variability in accretion-related emission lines for 40 Classical T Tauri stars to probe the extent of accretion variations in young stellar objects. Our analysis is based on multi-epoch high-resolution spectra for young stars in Tau-Aur and Cha I. For all stars, we obtain typically four spectra, covering timescales from hours to months. As proxies for the accretion rate, we use the H-alpha 10% width and the CaII-8662 line flux. We find that while the two quantities are correlated, their variability amplitude is not. Converted to accretion rates, the CaII fluxes indicate typical accretion rate changes of 0.35 dex, with 32% exceeding 0.5 dex, while H-alpha 10% width suggests changes of 0.65 dex, with 66% exceeding 0.5 dex. We conclude that CaII fluxes are a more robust quantitative indicator of accretion than H-alpha 10% width, and that intrinsic accretion rate changes typically do not exceed 0.5 dex on timescales of days to months. The maximum extent of the variability is reached after a few days, suggesting that rotation is the dominant cause of variability. We see a decline of the inferred accretion rates towards later spectral types, reflecting the dM/dt vs. M relationship. There is a gap between accretors and non-accretors, pointing to a rapid shutdown of accretion. We conclude that the ~2 orders of magnitude scatter in the dM/dt vs. M relationship is dominated by object-to-object scatter instead of intrinsic source variability.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

We present a comprehensive study of rotation, disk and accretion signatures for 144 T Tauri stars in the young (~2 Myr old) Chamaeleon I and Taurus-Auriga star forming regions based on multi-epoch high-resolution optical spectra from the Magellan Clay 6.5 m telescope supplemented by mid-infared photometry from the Spitzer Space Telescope. In contrast to previous studies in the Orion Nebula Cluster and NGC 2264, we do not see a clear signature of disk braking in Tau-Aur and Cha I. We find that both accretors and non-accretors have similar distributions of v sin i. The rotational velocities in both regions show a clear mass dependence, with F--K stars rotating on average about twice as fast as M stars, consistent with results reported for other clusters of similar age. Similarly, we find the upper envelope of the observed values of specific angular momentum j varies as M^0.5 for our sample which spans a mass range of ~0.16 to ~3 M_sun. This power law complements previous studies in Orion which estimated j is proportional to M^0.25 for < ~2 Myr stars in the same mass regime, and a sharp decline in j with decreasing mass for older stars (~10 Myr) with M < 2 M_sun. For a subsample of 67 objects with mid-IR photometry, we examine the connection between accretion signatures and dusty disks: in the vast majority of cases (63/67), the two properties correlate well, which suggests that the timescale of gas accretion is similar to the lifetime of inner disks.

Magnetic field estimates for nearby isolated neutron stars (INS) help to constrain both the characteristics of the population and the nature of their peculiar X-ray spectra. From a series of XMM-Newton observations of RX J2143.0+0654, we measure a spin-down rate of -4.6e-16 +/- 2.0e-16 Hz/s. While this does not allow a definitive measurement of the dipole magnetic field strength, fields of >1e14 G such as those inferred from the presence of a spectral absorption feature at 0.75keV are excluded. Instead, the field is most likely around 2e13 G, very similar to those of other INS. We not only suggest that this similarity most likely reflects the influence of magnetic field decay on this population, but also discuss a more speculative possibility that it results from peculiar conditions on the neutron-star surface. We find no evidence for spectral variability above the ~2% level. We confirm the presence of the 0.75-keV feature found earlier, and find tentative evidence for an additional absorption feature at 0.4 keV.

We present Gemini near-infrared adaptive optics imaging and spectroscopy of a planetary mass candidate companion to 1RXS J160929.1-210524, a roughly solar-mass member of the 5 Myr-old Upper Scorpius association. The object, separated by 2.22" or 330 AU at ~150 pc, has infrared colors and spectra suggesting a temperature of 1800(-100/+200) K, and spectral type of L4(-2/+1). The H- and K-band spectra provide clear evidence of low surface gravity, and thus youth. Based on the widely used DUSTY models, we infer a mass of 8(-2/+4)Mjupiter. If gravitationally bound, this would be the lowest mass companion imaged around a normal star thus far, and its existence at such a large separation would pose a serious challenge to theories of star and planet formation.

We estimate cluster ages from lithium depletion in five pre-main-sequence groups found within 100 pc of the Sun: TW Hydrae Association, Eta Chamaeleontis Cluster, Beta Pictoris Moving Group, Tucanae-Horologium Association and AB Doradus Moving Group. We determine surface gravities, effective temperatures and lithium abundances for over 900 spectra through least squares fitting to model-atmosphere spectra. For each group, we compare the dependence of lithium abundance on temperature with isochrones from pre-main-sequence evolutionary tracks to obtain model dependent ages. We find that the Eta Chamaelontis Cluster and the TW Hydrae Association are the youngest, with ages of 12+/-6 Myr and 12+/-8 Myr, respectively, followed by the Beta Pictoris Moving Group at 21+/-9 Myr, the Tucanae-Horologium Association at 27+/-11 Myr, and the AB Doradus Moving Group at an age of at least 45 Myr (where we can only set a lower limit since the models -- unlike real stars -- do not show much lithium depletion beyond this age). Here, the ordering is robust, but the precise ages depend on our choice of both atmospheric and evolutionary models. As a result, while our ages are consistent with estimates based on Hertzsprung-Russell isochrone fitting and dynamical expansion, they are not yet more precise. Our observations do show that with improved models, much stronger constraints should be feasible: the intrinsic uncertainties, as measured from the scatter between measurements from different spectra of the same star, are very low: around 10 K in effective temperature, 0.05 dex in surface gravity, and 0.03 dex in lithium abundance.

The Guitar Nebula is an H-alpha nebula produced by the interaction of the relativistic wind of a very fast pulsar, PSR B2224+65, with the interstellar medium. It consists of a ram-pressure confined bow shock near its head and a series of semi-circular bubbles further behind, the two largest of which form the body of the Guitar. We present a scenario in which this peculiar morphology is due to instabilities in the back flow from the pulsar bow shock. From simulations, these back flows appear similar to jets and their kinetic energy is a large fraction of the total energy in the pulsar's relativistic wind. We suggest that, like jets, these flows become unstable some distance down-stream, leading to rapid dissipation of the kinetic energy into heat, and the formation of an expanding bubble. We show that in this scenario the sizes, velocities, and surface brightnesses of the bubbles depend mostly on observables, and that they match roughly what is seen for the Guitar. Similar instabilities may account for features seen in other bow shocks.

We present optical time series spectroscopy of the pulsating white dwarf star G29-38 taken at the Very Large Telescope (VLT). By measuring the variations in brightness, Doppler shift, and line shape of each spectrum, we explore the physics of pulsation and measure the spherical degree ($\ell$) of each stellar pulsation mode. We measure the physical motion of the g-modes correlated with the brightness variations for three of the eight pulsation modes in this data set. The varying line shape reveals the spherical degree of the pulsations, an important quantity for properly modeling the interior of the star with asteroseismology. Performing fits to the H$\beta$, H$\gamma$, and H$\delta$ lines, we quantify the changing shape of the line and compare them to models and previous time series spectroscopy of G~29-38. These VLT data confirm several $\ell$ identifications and add four new values, including an additional $\ell$=2 and a possible $\ell$=4. In total from both sets of spectroscopy of G29-38, eleven modes now have known spherical degrees.

We present the results of a multiplicity survey of 126 stars spanning ~0.1-3 solar masses in the ~2-Myr-old Chamaeleon I star-forming region, based on adaptive optics imaging with the ESO Very Large Telescope. Our observations have revealed 30 binaries and 6 triples, of which 19 and 4, respectively, are new discoveries. The overall multiplicity fraction we find for Cha I (~30%) is similar to those reported for other dispersed young associations, but significantly higher than seen in denser clusters and the field, for comparable samples. Both the frequency and the maximum separation of Cha I binaries decline with decreasing mass, while the mass ratios approach unity; conversely, tighter pairs are more likely to be equal mass. We confirm that brown dwarf companions to stars are rare, even at young ages at wide separations. Based on follow-up spectroscopy of two low-mass substellar companion candidates, we conclude that both are likely background stars. The overall multiplicity fraction in Cha I is in rough agreement with numerical simulations of cloud collapse and fragmentation, but its observed mass dependence is less steep than predicted. The paucity of higher-order multiples, in particular, provides a stringent constraint on the simulations, and seems to indicate a low level of turbulence in the prestellar cores in Cha I.

RX J1856.5-3754 is the X-ray brightest among the nearby isolated neutron stars. Its X-ray spectrum is thermal, and is reproduced remarkably well by a black-body, but its interpretation has remained puzzling. One reason is that the source did not exhibit pulsations, and hence a magnetic field strength--vital input to atmosphere models--could not be estimated. Recently, however, very weak pulsations were discovered. Here, we analyze these in detail, using all available data from the XMM-Newton and Chandra X-ray observatories. From frequency measurements, we set a 2-sigma upper limit to the frequency derivative of \dot\nu<1.3e-14 Hz/s. Trying possible phase-connected timing solutions, we find that one solution is far more likely than the others, and we infer a most probable value of \dot\nu=(-5.98+/-0.14)e-16 Hz/s. The inferred magnetic field strength is 1.5e13 G, comparable to what was found for similar neutron stars. From models, the field seems too strong to be consistent with the absence of spectral features for non-condensed atmospheres. It is sufficiently strong, however, that the surface could be condensed, but only if it is consists of heavy elements like iron. Our measurements imply a characteristic age of about 4 Myr. This is longer than the cooling and kinematic ages, as was found for similar objects, but at almost a factor ten, the discrepancy is more extreme. A puzzle raised by our measurement is that the implied rotational energy loss rate of about 3e30 erg/s is orders of magnitude smaller than what was inferred from the H-alpha nebula surrounding the source.

After our tentative detection of an optical counterpart to CXOU J010043.1-721134 from archival Hubble Space Telescope (HST) imaging, we have followed up with further images in four bands. Unfortunately, the source originally identified is not confirmed. We provide deep photometric limits in four bands and accurate photometry of field stars around the location of the magnetar.

Characterizing multiplicity in the very low mass (VLM) domain is a topic of much current interest and fundamental importance. Here we report on a near-IR AO imaging survey of 31 young brown dwarfs and VLM stars, 28 of which are in Chamaeleon I, using the ESO VLT. Our survey is sensitive enough to detect equal mass binaries down to separations of 0.04-0.07" (6-10 AU at 160 pc) and, typically, companions with mass ratios as low as 0.2 outside of 0.2" (~30 AU). We resolve the suspected 0.16" (~26 AU) binary Cha_Halpha 2 and present two new binaries, Hn 13 and CHXR 15, with separations of 0.13" (~20 AU) and 0.30" (~50 AU); the latter is one of the widest VLM systems discovered to date. We do not find companions around the majority of our targets, giving a binary frequency of 11% [+9,-6], thus confirming the trend for a lower binary frequency with decreasing mass. By combining our work with previous surveys, we arrive at the largest sample of young VLMOs (72) with high angular resolution imaging to date. Its multiplicity fraction is in statistical agreement with that for VLMOs in the field. In addition we note that many field stellar binaries with lower binding energies and/or wider cross sections have survived dynamical evolution and that statistical models suggest tidal disruption by passing stars is unlikely to affect the binary properties of our systems. Thus, we argue that there is no significant evolution of multiplicity with age among brown dwarfs and VLM stars in OB and T associations between a few Myr to several Gyr. Instead, the observations to date suggest that VLM objects are either less likely to be born in fragile multiple systems than solar mass stars or such systems are disrupted very early.

We present a comprehensive study of disks around 81 young low-mass stars and brown dwarfs in the nearby ~2-Myr-old Chamaeleon I star-forming region. We use mid-infrared photometry from the Spitzer Space Telescope, supplemented by findings from ground-based high-resolution optical spectroscopy and adaptive optics imaging. We derive disk fractions of 52 (+/-6) % and 58 (+6/-7) % based on 8-micron and 24-micron colour excesses, respectively, consistent with those reported for other clusters of similar age. Within the uncertainties, the disk frequency in our sample of K3-M8 objects in Cha I does not depend on stellar mass. Diskless and disk-bearing objects have similar spatial distributions. There are no obvious transition disks in our sample, implying a rapid timescale for the inner disk clearing process; however, we find two objects with weak excess at 3-8 microns and substantial excess at 24 microns, which may indicate grain growth and dust settling in the inner disk. For a sub-sample of 35 objects with high-resolution spectra, we investigate the connection between accretion signatures and dusty disks: in the vast majority of cases (29/35) the two are well correlated, suggesting that, on average, the timescale for gas dissipation is similar to that for clearing the inner dust disk. The exceptions are six objects for which dust disks appear to persist even though accretion has ceased or dropped below measurable levels. Adaptive optics images of 65 of our targets reveal that 17 have companions at (projected) separations of 10-80 AU. Of the five <20 AU binaries, four lack infrared excess, possibly indicating that a close companion leads to faster disk dispersal. The closest binary with excess is separated by ~20 AU, which sets an upper limit of ~8 AU for the outer disk radius. (abridged)

We present the results of a spectroscopic monitoring campaign of the OB-star companions to the eclipsing X-ray pulsars SMC X-1, LMC X-4 and Cen X-3. High-resolution optical spectra obtained with UVES on the ESO Very Large Telescope are used to determine the radial-velocity orbit of the OB (super)giants with high precision. The excellent quality of the spectra provides the opportunity to measure the radial-velocity curve based on individual lines, and to study the effect of possible distortions of the line profiles due to e.g. X-ray heating on the derived radial-velocity amplitude. Several spectral lines show intrinsic variations with orbital phase. The magnitude of these variations depends on line strength, and thus provides a criterion to select lines that do not suffer from distortions. The undistorted lines show a larger radial-velocity amplitude than the distorted lines, consistent with model predictions. Application of our line-selection criteria results in a mean radial-velocity amplitude K(Opt) of 20.2 +/- 1.1, 35.1 +/- 1.5, and 27.5 +/- 2.3 km/s (1 sigma errors), for the OB companion to SMC X-1, LMC X-4 and Cen X-3, respectively. Adding information on the projected rotational velocity of the OB companion (derived from our spectra), the duration of X-ray eclipse and orbital parameters of the X-ray pulsar (obtained from literature), we arrive at a neutron star mass of 1.06^{+0.11}_{-0.10}, 1.25^{+0.11}_{-0.10} and 1.34^{+0.16}_{-0.14} M{sun} for SMC X-1, LMC X-4 and Cen X-3, respectively. The mass of SMC X-1 is near the minimum mass (~1 M{sun}) expected for a neutron star produced in a supernova. We discuss the implications of the measured mass distribution on the neutron-star formation mechanism, in relation to the evolutionary history of the massive binaries.

We have observed the 5.54s anomalous X-ray pulsar XTE J1810-197 at radio, millimeter, and infrared (IR) wavelengths, with the aim of learning about its broad-band spectrum. At the IRAM 30m telescope, we have detected the magnetar at 88 and 144GHz, the highest radio-frequency emission ever seen from a pulsar. At 88GHz we detected numerous individual pulses, with typical widths ~2ms and peak flux densities up to 45Jy. Together with nearly contemporaneous observations with the Parkes, Nancay, and Green Bank telescopes, we find that in late 2006 July the spectral index of the pulsar was -0.5<alpha<0 over the range 1.4-144GHz. Nine dual-frequency Very Large Array and Australia Telescope Compact Array observations in 2006 May-September are consistent with this finding, while showing variability of alpha with time. We infer from the IRAM observations that XTE J1810-197 remains highly linearly polarized at millimeter wavelengths. Also, toward this pulsar, the transition frequency between strong and weak scattering in the interstellar medium may be near 50GHz. At Gemini, we detected the pulsar at 2.2um in 2006 September, at the faintest level yet observed, K_s=21.89+-0.15. We have also analyzed four archival IR Very Large Telescope observations (two unpublished), finding that the brightness fluctuated within a factor of 2-3 over a span of 3 years, unlike the monotonic decay of the X-ray flux. Thus, there is no correlation between IR and X-ray flux, and it remains uncertain whether there is any correlation between IR and radio flux.

We present results of the CFHT adaptive optics search for companions of a homogeneous group of contact binary stars, as a contribution to our attempts to prove a hypothesis that these binaries require a third star to become so close as observed. In addition to companions directly discovered at separations of >=1", we introduced a new method of AO image analysis utilizing distortions of the AO diffraction ring pattern at separations of 0.07"-1". Very close companions, with separations in the latter range were discovered in systems HV Aqr, OO Aql, CK Boo, XY Leo, BE Scl, and RZ Tau. More distant companions were detected in V402 Aur, AO Cam, V2082 Cyg. Our results provide a contribution to the mounting evidence that the presence of close companions is a very common phenomenon for very close binaries with orbital periods <1 day.

We have used a set of dedicated astrometric data from the Hubble Space Telescope to measure the parallax and proper motion of the nearby neutron star RX J0720.4-3125. At each of eight epochs over two years, we used the High Resolution Camera of the Advanced Camera for Surveys to measure the position of the B=26.6 target to a precision of ~2 mas (~0.07 pix) relative to 22 other stars. From these data we measure a parallax of pi=2.8+/-0.9 mas (for a distance of 360+170-90 pc) and a proper motion of mu=107.8+/-1.2 mas/yr. Exhaustive testing of every stage of our analysis suggests that it is robust, with a maximum systematic uncertainty on the parallax of 0.4 mas. The distance is compatible with earlier estimates made from scaling the optical emission of RX J0720.4-3125 relative to the even closer neutron star RX J1856.5-3754. The distance and proper motion imply a transverse velocity of 180+90-40 km/s, comparable to velocities observed for radio pulsars. The speed and direction suggest an origin for RX J0720.4-3125 in the Trumpler 10 OB association ~0.7 Myr ago, with a possible range of 0.5-1.0 Myr given by the uncertainty in the distance.

We have obtained observations of the magnetar XTE J1810-197 with the Very Long Baseline Array at two epochs separated by 106 days, at wavelengths of 6 cm and 3.6 cm. Comparison of the positions yields a proper motion value of 13.5+-1.0 mas/yr at an equatorial position angle of 209.4+-2.4 deg (east of north). This value is consistent with a lower-significance proper motion value derived from infrared observations of the source over the past three years, also reported here. Given its distance of 3.5+-0.5 kpc, the implied transverse velocity corrected to the local standard of rest is 212+-35 km/s (1 sigma). The measured velocity is slightly below the average for normal young neutron stars, indicating that the mechanism(s) of magnetar birth need not lead to high neutron star velocities. We also use Australia Telescope Compact Array, Very Large Array, and these VLBA observations to set limits on any diffuse emission associated with the source on a variety of spatial scales, concluding that the radio emission from XTE J1810-197 is >96% pulsed.

RX J0720.4-3125 is an isolated neutron star that, uniquely in its class, has shown changes in its thermal X-ray spectrum. We use new spectra taken with Chandra's Low Energy Transmission Grating Spectrometer, as well as archival observations, to try to understand the timescale and nature of these changes. We construct lightcurves, which show both small, slow variations on a timescale of years, and a larger event that occurred more quickly, within half a year. From timing, we find evidence for a `glitch' coincident with this larger event, with a fractional increase in spin frequency of 5x10^{-8}. We compare the `before' and `after' spectra with those from RX J1308.6+2127, an isolated neutron star with similar temperature and magnetic field strength, but with a much stronger absorption feature in its spectrum. We find that the `after' spectrum can be represented remarkably well by the superposition of the `before' spectrum, scaled by two thirds, and the spectrum of RX J1308.6+2127, thus suggesting that the event affected approximately one third of the surface. We speculate the event reflects a change in surface composition caused by, e.g., an accretion episode.

We have collected data spanning seven years of observations of the magnetar 4U 0142+61 in the infrared, optical and soft X-rays. These combine our own observations and analysis of archival data. We find that the source is variable in the optical, in contrast to what had been previously reported, that the K-band flux can vary by over a magnitude on the time-scale of days, and that the X-ray pulsed flux is not obviously correlated with either the total X-ray flux or infrared and optical fluxes. Furthermore, from multi-color photometry of the source within single nights, we conclude that there are two separate components to the infrared emission. The overall picture is unclear, and prompts the need for further, more frequent observations.

We present timing measurements, astrometry, and high-resolution spectra of a number of nearby, thermally emitting, isolated neutron stars. We use these to infer magnetic field strengths and distances, but also encounter a number of puzzles. We discuss three specific ones in detail: (i) For RX J0720.4-3125 and RX J1308.6+2127, the characteristic ages are in excess of 1~Myr, while their temperatures and kinematic ages indicate that they are much younger; (ii) For RX J1856.5-3754, the brightness temperature for the optical emission is in excess of that measured at X-ray wavelengths for reasonable neutron-star radii; (iii) For RX J0720.4-3125, the spectrum changed from an initially featureless state to one with an absorption feature, yet there was only a relatively small change in effective temperature. Furthermore, we attempt to see whether the spectra of all seven sourced, in six of which absorption features have now been found, can be understood in the context of strongly magnetised hydrogen atmospheres. We find that the energies of the absorption features can be reproduced, but that the featureless spectra of some sources, especially the Wien-like high-energy tails, remain puzzling.

The X-ray spectra of Anomalous X-ray Pulsars have long been fit by smooth, empirical models such as the sum of a black-body plus a power law. These reproduce the ~0.5 to 10 keV range well, but fail at lower and higher energies, grossly over-predicting the optical and under-predicting the hard X-ray emission. A poorly constrained source of uncertainty in determining the true, intrinsic spectra, in particular at lower energies, is the amount of interstellar extinction. In previous studies, extinction column densities with small statistical errors were derived as part of the fits of the spectra to simple continuum models. Different choices of model, however, each produced statistically acceptable fits, but a wide range of columns. Here, we attempt to measure the interstellar extinction in a model-independent way, using individual absorption edges of the elements O, Fe, Ne, Mg and Si in X-ray grating spectra taken with XMM-Newton. We find that our inferred equivalent hydrogen column density NH for 4U 0142+61 is a factor of 1.4 lower than the typically quoted value from black-body plus power-law fits, and is now consistent with estimates based on the dust scattering halo and visual extinction. For three other sources, we find column densities consistent with earlier estimates. We use our measurements to recover the intrinsic spectra of the AXPs empirically, without making assumptions on what the intrinsic spectral shapes ought to be. We find that the power-law components that dominate at higher energies do not extend below the thermal peak.

We present spectroscopic and photometric observations of the optical counterpart to PSR J1911-5958A, a millisecond pulsar located towards the globular cluster NGC 6752. We measure radial velocities from the spectra and determine the systemic radial velocity of the binary and the radial-velocity amplitude of the white-dwarf orbit. Combined with the pulsar orbit obtained from radio timing, we infer a mass ratio of Mpsr/Mwd=7.36+-0.25. The spectrum of the counterpart is that of a hydrogen atmosphere, showing Balmer absorption lines upto H12, and we identify the counterpart as a helium-core white dwarf of spectral type DA5. Comparison of the spectra with hydrogen atmosphere models yield a temperature Teff=10090+-150 K and a surface gravity log g=6.44+-0.20 cm s^-2. Using mass-radius relations appropriate for low-mass helium-core white dwarfs, we infer the white-dwarf mass Mwd=0.18+-0.02 Msun and radius Rwd=0.043+-0.009 Rsun. Combined with the mass ratio, this constrains the pulsar mass to Mpsr=1.40^+0.16_-0.10 Msun. If we instead use the white-dwarf spectrum and the distance of NGC 6752 to determine the white-dwarf radius, we find Rwd=0.058+-0.004 Rsun. For the observed temperature, the mass-radius relations predict a white-dwarf mass of Mwd=0.175+-0.010 Msun, constraining the pulsar mass to Mpsr=1.34+-0.08 Msun. We find that the white-dwarf radius determined from the spectrum and the systemic radial velocity of the binary are only marginally consistent with the values that are expected if PSR J1911-5958A is associated with NGC 6752. We discuss possible causes to explain this inconsistency, but conclude that our observations do not conclusively confirm nor disprove the assocation of the pulsar binary with the globular cluster.

We identify `red clump stars' - core helium-burning giants - among 2MASS stars and use them to measure the run of reddening with distance in the direction of each of the Galactic Anomalous X-ray Pulsars (AXP). We combine this with extinction estimates from X-ray spectroscopy to infer distances and find that the locations of all AXP are consistent with being in Galactic spiral arms. We also find that the 2-10 keV luminosities implied by our distances are remarkably similar for all AXPs, being all around 1.3e35 erg s-1. Furthermore, using our distances to estimate effective black-body emitting radii, we find that the radii are tightly anti-correlated with pulsed fraction, and somewhat less tightly anti-correlated with black-body temperature. We find no obvious relationship of any property with the dipole magnetic field strength inferred from the spin-down rate.

We report the findings of a comprehensive study of disk accretion and related phenomena in four of the nearest young stellar associations spanning 6-30 million years in age, an epoch that may coincide with the late stages of planet formation. We have obtained ~650 multi-epoch high-resolution optical spectra of 100 low-mass stars that are likely members of the eta Chamaeleontis (~6 Myr), TW Hydrae (~8 Myr), beta Pictoris (~12 Myr) and Tucanae-Horologium (~30 Myr) groups. Our data were collected over 12 nights between 2004 December - 2005 July on the Magellan Clay 6.5m telescope. Based on H$\alpha$ line profiles, along with a variety of other emission lines, we find clear evidence of on-going accretion in three out of 11 eta Cha stars and two out of 32 TW Hydrae members. None of the 57 beta Pic or Tuc-Hor members shows measurable signs of accretion. Together, these results imply significant evolution of the disk accretion process within the first several Myr of a low-mass star's life. While a few disks can continue to accrete for up to ~10 Myr, our findings suggest that disks accreting for beyond that timescale are rather rare. This result provides an indirect constraint on the timescale for gas dissipation in inner disks and, in turn, on gas giant planet formation. All accretors in our sample are slow rotators, whereas non-accretors cover a large range in rotational velocities. This may hint at rotational braking by disks at ages up to ~8 Myr. Our multi-epoch spectra confirm that emission-line variability is common even in somewhat older T Tauri stars, among which accretors tend to show particularly strong variations. Thus, our results indicate that accretion and wind activity undergo significant and sustained variations throughout the lifetime of accretion disks.

We have analyzed both archival and new infrared imaging observations of the field of the Anomalous X-ray Pulsar 1RXS J170849.0-400910, in search of the infrared counterpart. This field has been previously investigated, and one of the sources consistent with the position of the AXP suggested as the counterpart. We, however, find that this object is more likely a background star, while another object within the positional error circle has non-stellar colors and shows evidence for variability. These two pieces of evidence, along with a consistency argument for the X-ray-to-infrared flux ratio, point to the second source being the more likely infrared counterpart to the AXP.

We report the detection of probable optical counterparts for two Intermediate Mass Binary Pulsar (IMBP) systems, PSR J1528-3146 and PSR J1757-5322. Recent radio pulsar surveys have uncovered a handful of these systems with putative massive white dwarf companions, thought to have an evolutionary history different from that of the more numerous class of Low Mass Binary Pulsars (LMBPs) with He white dwarf companions. The study of IMBP companions via optical observations offers us several new diagnostics: the evolution of main sequence stars near the white-dwarf-neutron star boundary, the physics of white dwarfs close to the Chandrasekhar limit, and insights into the recycling process by which old pulsars are spun up to high rotation frequencies. We were unsuccessful in our attempt to detect optical counterparts of PSR J1141-6545, PSR J1157-5112, PSR J1435-6100, and PSR J1454-5846.

It is unclear how very close binary stars form, given that during the pre-main-sequence phase the component stars would have been inside each other. One hypothesis is that they formed further apart, but were brought in closer after formation by gravitational interaction with a third member of the system. If so, all close binaries should be members of triple (or higher-order) systems. As a test of this prediction, we present a search for the signature of third components in archival spectra of close binaries. In our sample of 75 objects, 23 show evidence for the presence of a third component, down to a detection limit of tertiary flux contributions of about 0.8% at 5200 \AA\ (considering only contact and semi-detached binaries, we find 20 out of 66). In a homogeneous subset of 58 contact binaries, we are fairly confident that the 15 tertiaries we have detected are all tertiaries present with mass ratios $0.28\lesssim M_3/M_{12}\lesssim0.75$ and implied outer periods $P\lesssim10^6{\rm d}$. We find that if the frequency of tertiaries was the same as that of binary companions to solar-type stars, we would have expected to detect about 12 tertiaries. In contrast, if all contact binaries were in triple systems, one would expect about 20. Thus, our results are not conclusive, but sufficiently suggestive to warrant further studies.

We present optical and near-infrared observations with Keck of the binary millisecond pulsar PSR J0751+1807. We detect a faint, red object - with R=25.08+-0.07, B-R=2.5+-0.3, and R-I=0.90+-0.10 - at the celestial position of the pulsar and argue that it is the white dwarf companion of the pulsar. The colours are the reddest among all known white dwarfs, and indicate a very low temperature, Teff~4000 K. This implies that the white dwarf cannot have the relatively thick hydrogen envelope that is expected on evolutionary grounds. Our observations pose two puzzles. First, while the atmosphere was expected to be pure hydrogen, the colours are inconsistent with this composition. Second, given the low temperature, irradiation by the pulsar should be important, but we see no evidence for it. We discuss possible solutions to these puzzles.

We present a phase-connected timing solution for the nearby isolated neutron star RX J1308.6+2127 (RBS 1223). From dedicated Chandra observations as well as archival Chandra and XMM-Newton data spanning a period of five years, we demonstrate that the 10.31-sec pulsations are slowing down steadily at a rate of Pdot=1.120(3)e-13 s/s. Under the assumption that this is due to magnetic dipole torques, we infer a characteristic age of 1.5 Myr and a magnetic field strength of 3.4e13 G. As with RX J0720.4-3125, the only other radio-quiet thermally emitting isolated neutron star for which a timing solution has been derived, the field strength is consistent with what was inferred earlier from the presence of a strong absorption feature in its X-ray spectrum. Furthermore, both sources share that the characteristic age is in excess of the cooling age inferred from standard cooling models. The sources differ, however, in their timing noise: while RX J0720.4-3125 showed considerable timing noise, RX J1308.6+2127 appears relatively stable.

We present high-speed, multi-colour optical photometry of the anomalous X-ray pulsar 4U 0142+61, obtained with ULTRACAM on the 4.2-m William Herschel Telescope. We detect 4U 0142+61 at magnitudes of i'=23.7+-0.1, g'=27.2+-0.2 and u'>25.8, consistent with the magnitudes found by Hulleman et al.(2004) and hence confirming their discovery of both a spectral break in the optical and a lack of long-term optical variability. We also confirm the discovery of Kern & Martin (2002) that 4U 0142+61 shows optical pulsations with an identical period (~8.7 s) to the X-ray pulsations. The rms pulsed fraction in our data is 29+-8%, 5-7 times greater than the 0.2-8 keV X-ray rms pulsed fraction. The optical and X-ray pulse profiles show similar morphologies and appear to be approximately in phase with each other, the former lagging the latter by only 0.04+-0.02 cycles. In conjunction with the constraints imposed by X-ray observations, the results presented here favour a magnetar interpretation for the anomalous X-ray pulsars.

We present the results of a dedicated effort to measure the spin-down rate of the nearby isolated neutron star RX J0720.4-3125. Comparing arrival times of the 8.39-sec pulsations for data from Chandra we derive an unambiguous timing solution for RX J0720.4-3125 that is accurate to <0.1 cycles over >5 years. Adding data from XMM and ROSAT, the final solution yields Pdot=(6.98+/-0.02)x10^(-14) s/s; for dipole spin-down, this implies a characteristic age of 2 Myr and a magnetic field strength of 2.4e13 G. The phase residuals are somewhat larger than those for purely regular spin-down, but do not show conclusive evidence for higher-order terms or a glitch. From our timing solution as well as recent X-ray spectroscopy, we concur with recent suggestions that RX J0720.4-3125 is most likely an off-beam radio pulsar with a moderately high magnetic field.

We present optical and infrared observations taken with the Very Large Telescope of the eclipsing binary pulsar system PSR B1718-19. The candidate companion of the pulsar, identified earlier in Hubble Space Telescope observations, has been detected in all three bands, R, I, and J. These detections allowed us to derive constraints on temperature, radius, and mass, pointing to a companion that has expanded to a radius between one of a main sequence star and one at the Roche-limit. We focus on the role of tidal circularisation in the system, which will have transformed the initially eccentric orbit expected from formation scenarios into the nearly circular orbit presently observed. Based on simple energy balance arguments, we are able to draw a picture of the companion's evolution resulting from the energy deposition in the star due to circularisation. In this picture, our measurement of the companion's parameters is consistent with the expected initial eccentricity. However, with the present understanding of tidal dissipation it remains difficult to account for the short time in which the system was circularised.

We report the precise optical and X-ray localization of the 3.2 ms accretion-powered X-ray pulsar XTE J1814-338 with data from the Chandra X-Ray Observatory as well as optical observations conducted during the 2003 June discovery outburst. Optical imaging of the field during the outburst of this soft X-ray transient reveals an R = 18 star at the X-ray position. This star is absent (R > 20) from an archival 1989 image of the field and brightened during the 2003 outburst, and we therefore identify it as the optical counterpart of XTE J1814-338. The best source position derived from optical astrometry is R.A. = 18h13m39.s04, Dec.= -33d46m22.3s (J2000). The featureless X-ray spectrum of the pulsar in outburst is best fit by an absorbed power-law (with photon index = 1.41 +- 0.06) plus blackbody (with kT = 0.95 +- 0.13 keV) model, where the blackbody component contributes approximately 10% of the source flux. The optical broad-band spectrum shows evidence for an excess of infrared emission with respect to an X-ray heated accretion disk model, suggesting a significant contribution from the secondary or from a synchrotron-emitting region. A follow-up observation performed when XTE J1814-338 was in quiescence reveals no counterpart to a limiting magnitude of R = 23.3. This suggests that the secondary is an M3 V or later-type star, and therefore very unlikely to be responsible for the soft excess, making synchroton emission a more reasonable candidate.

We present the results of a monitoring campaign of three eclipsing high-mass X-ray binaries (HMXBs: SMC X-1, LMC X-4 and Cen X-3). High-resolution VLT/UVES spectra are used to measure the radial velocities of these systems with high accuracy. We show that the subsequent mass determination of the neutron stars in these systems is significantly improved and discuss the implications of this result.

We present photometry of the Anomalous X-ray pulsar 1E1048.1-5937 in the infrared and optical, taken at Magellan and the VLT. The object is detected in the I, J and Ks bands under excellent conditions. We find that the source has varied greatly in its infrared brightness and present these new magnitudes. No correlation is found between the infrared flux and spin-down rate, but the infrared flux and X-ray flux may be anti-correlated. Assuming nominal reddening values, the resultant spectral energy distribution is found to be inconsistent with the only other AXP SED available (for 4U0142+61). We consider the effect of the uncertainty in the reddening to the source on its SED. We find that although both the X-ray and infrared fluxes have varied greatly for this source, the most recent flux ratio is remarkably consistent with what is is found for other AXPs. Finally, we discuss the implications of our findings in the context of the magnetar model.

We present the results of a near-infrared monitoring program of the Anomalous X-ray Pulsar 1E 2259+586, performed at the Gemini Observatory. This program began three days after the pulsar's 2002 June outburst, and spans ~1.5 years. We find that after an initial increase associated with the outburst, the near-infrared flux decreased continually and reached the pre-burst quiescent level after about one year. We compare both the near-infrared flux enhancement and its decay to those of the X-ray afterglow, and find them to be remarkably consistent. Fitting simple power laws to the RXTE pulsed flux and near-infrared data for t>1 day post-burst, we find the following decay indices: alpha=-0.21+/-0.01 (X-ray), alpha=-0.21+/-0.02 (near-infrared), where flux is a function of time such that F is proportional to t^alpha. This suggests that the enhanced infrared and X-ray fluxes have a physical link post-outburst, most likely from the neutron-star magnetosphere.

Optical observations of the companions of pulsars can help determine the properties of the binaries, as well as those of their components, and give clues to the preceding evolution. In this review, we first describe the different classes of binary pulsars, and present a table with a summary what is known about their optical counterparts. Next, we focus on the class of pulsars that have low-mass, helium-core white dwarf companions. We discuss attempts to determine the masses of both components using optical spectroscopy, and compare the pulsar spin-down ages with cooling ages of the white dwarfs. We confirm that for a given age, the lowest-mass white dwarfs are much hotter than the more massive ones, consistent with recent evolutionary models, although with one glaring exception. We discuss the case of PSR B0820+02, where the cooling age indicates a braking index less than 3, and conclude by describing how cooling ages can be used to test formation scenarios for PSR J1911-5958A, a pulsar binary in the outskirts of NGC 6752.

The pulsating white dwarf star PY Vul (G~185-32) exhibits pulsation modes with peculiar properties that set it apart from other variable stars in the ZZ Ceti (DAV) class. These peculiarities include a low total pulsation amplitude, a mode with bizarre amplitudes in the ultraviolet, and a mode harmonic that exceeds the amplitude of its fundamental. Here, we present optical, time series spectroscopy of PY Vul acquired with the Keck II LRIS spectrograph. Our analysis has revealed that the mode with unusual UV amplitudes also has distinguishing characteristics in the optical. Comparison of its line profile variations to models suggests that this mode has a spherical degree of four. We show that all the other peculiarities in this star are accounted for by a dominant pulsation mode of l=4, and propose this hypothesis as a solution to the mysteries of PY Vul.

I review attempts made to determine the properties of neutron stars. I focus on constraints on the maximum mass that a neutron star can have, and on attempts to measure neutron-star radii. So far, there appears to be only one neutron star for which there is strong evidence that its mass is above the canonical 1.4 Msun, viz., Vela X-1, for which a mass close to 1.9 Msun is found. Prospects for progress appear brightest for studies of systems in which the neutron star should have accreted substantial amounts of matter. While for individual systems the evidence that neutron stars can have high masses is weak, the ensemble appears to show that masses around 1.6 Msun are possible. For the radius determination, most attempts have focussed on neutron stars in low-mass X-ray binaries in which accretion has temporarily shut down. These neutron stars are easiest to model, since they should have pure Hydrogen atmospheres and low magnetic fields. To obtain accurate radii, however, requires precise distances and very high quality data.

We present X-ray spectra taken with XMM-Newton of RX J1605.3+3249, the third brightest in the class of nearby, thermally emitting neutron stars. In contrast to what is the case for the brightest object, RX J1856.5-3754, we find that the spectrum of RX J1605.3+3249 cannot be described well by a pure black body, but shows a broad absorption feature at 27\AA (0.45 keV). With this, it joins the handful of isolated neutron stars for which spectral features arising from the surface have been detected. We discuss possible mechanisms that might lead to the features, as well as the overall optical to X-ray spectral energy distribution, and compare the spectrum with what is observed for the other nearby, thermally emitting neutron stars. We conclude that we may be observing absorption due to the proton cyclotron line, as was suggested for the other sources, but weakened due to the strong-field quantum electrodynamics effect of vacuum resonance mode conversion.

We report the discovery of PSR J1909-3744, a 2.95 millisecond pulsar in a nearly circular 1.53 day orbit. Its narrow pulse width of 43 microseconds allows pulse arrival times to be determined with great accuracy. We have spectroscopically identified the companion as a moderately hot (T = 8500 K) white dwarf with strong absorption lines. Radial velocity measurements of the companion will yield the mass ratio of the system. Our timing data suggest the presence of Shapiro delay; we expect that further timing observations, combined with the mass ratio, will allow the first accurate determination of a millisecond pulsar mass. We have measured the timing parallax and proper motion for this pulsar which indicate a transverse velocity of 140 (+80/-40) km/s. This pulsar's stunningly narrow pulse profile makes it an excellent candidate for precision timing experiments that attempt to detect low frequency gravitational waves from coalescing supermassive black hole binaries.

I describe recent high-resolution X-ray spectroscopy of surface emission from nearby, thermally emitting neutron stars. I focus on RX J0720.4-3125, RX J1308.6+2127, and RX J1605.3+3249, all of which have similar temperature, but differ in the presence and strength of absorption features in their spectra. I discuss possible causes for the absorption we see in two sources, and conclude that it may be proton cyclotron line absorption, but weakened due to the strong-field quantum electrodynamics effect of vacuum resonance mode conversion.

We present time-resolved optical spectrophotometry of the pulsating hydrogen atmosphere (DA) white dwarf G 117-B15A. We find three periodicities in the pulsation spectrum (215s, 272s, and 304s) all of which have been found in earlier studies. By comparing the fractional wavelength dependence of the pulsation amplitudes (chromatic amplitudes) with models, we confirm a previous report that the strongest mode, at 215s, has l=1. The chromatic amplitude for the 272s mode is very puzzling, showing an increase in fractional amplitude with wavelength that cannot be reproduced by the models for any l at optical wavelengths. Based on archival HST data, we show that while the behaviour of the 215s mode at ultra-violet wavelengths is as expected from models, the weird behaviour of the 272s periodicity is not restricted to optical wavelengths in that it fails to show the expected increase in fractional amplitude towards shorter wavelengths. We discuss possible causes for the discrepancies found for the 272s variation, but find that all are lacking, and conclude that the nature of this periodicity remains unclear.

The detection of optical/infrared counterparts to Anomalous X-ray Pulsars (AXPs) has greatly increased our understanding of these systems. Models for the AXP phenomenon were based upon their X-ray emission, and all but the magnetar model made predictions for the optical/infrared that have now been falsified. With detections in hand, detailed studies of the optical/infrared to X-ray flux ratios, variability, and the spectral energy distributions have become possible. We present new data on two AXPs taken with Keck and Magellan, and compare the results with predictions made in the context of the magnetar model, in which the emission is due to ion currents flowing in the >10e14 G magnetosphere of young neutron stars.

We report on the identification of the optical counterpart of the binary millisecond pulsar PSR J1911-5958A, located in the outskirts of the globular cluster NGC 6752. At the position of the pulsar we find an object with V=22.08, B-V=0.38, U-B=-0.49. The object is blue with respect to the cluster main sequence by 0.8 magnitudes in B-V. We argue that the object is the white dwarf companion of the pulsar. Comparison with white dwarf cooling models shows that this magnitude and colors are consistent with a low-mass white dwarf at the distance of NGC 6752. If associated with NGC 6752, the white dwarf is relatively young, <2 Gyr, which sets constraints on the formation of the binary and its ejection from the core of the globular cluster.

We have re-analysed all available high-resolution ultraviolet IUE spectra of the high-mass X-ray binary HD153919/4U1700-37. The radial velocity semi-amplitude of 20.6 +/- 1.0 km/s and orbital eccentricity of 0.22 +/- 0.04 agree very well with the values obtained earlier from optical spectra. They disagree with earlier conclusions for the same data reduced by Heap & Corcoran (1992) and by Stickland & Lloyd (1993).

We have detected the optical counterpart to the nearby isolated neutron star RX J1605.3+3249 using observations from the Space Telescope Imaging Spectrometer aboard the Hubble Space Telescope. The counterpart, with m_{50CCD}=26.84+/-0.07 mag and very blue colors, lies close to the ROSAT HRI error circle and within the Chandra error circle. The spectrum is consistent with a Rayleigh-Jeans tail whose emission is a factor of ~14 above the extrapolation of the X-ray blackbody, and the source has an unabsorbed X-ray-to-optical flux ratio of log(fX/fopt)=4.4, similar to that of other isolated neutron stars. This confirms the classification of RX J1605.3+3249 as a neutron star.

We present the analysis of time-resolved, high resolution spectra of the cool white dwarf pulsator, G29-38. From measuring the Doppler shifts of the H-alpha core, we detect velocity changes as large as 16.5 km/s and conclude that they are due to the horizontal motions associated with the g-mode pulsations on the star. We detect seven pulsation modes from the velocity time-series and identify the same modes in the flux variations. We discuss the properties of these modes and use the advantage of having both velocity and flux measurements of the pulsations to test the convective driving theory proposed for DAV stars. Our data show limited agreement with the expected relationships between the amplitude and phases of the velocity and flux modes. Unexpectedly, the velocity curve shows evidence for harmonic distortion, in the form of a peak in the Fourier transform whose frequency is the exact sum of the two largest frequencies. Combination frequencies are a characteristic feature of the Fourier transforms of light curves of G29-38, but before now have not been detected in the velocities, nor does published theory predict that they should exist. We compare our velocity combination frequency to combination frequencies found in the analysis of light curves of G29-38, and discuss what might account for the existence of velocity combinations with the properties we observe. We also use our high-resolution spectra to determine if either rotation or pulsation can explain the truncated shape observed for the DAV star's line core. We are able to eliminate both mechanisms: the average spectrum does not fit the rotationally broadened model and the time-series of spectra provides proof that the pulsations do not significantly truncate the line.

We report on Keck optical BVRI images and spectroscopy of the companion of the binary millisecond pulsar PSR J0218+4232. A faint bluish (V=24.2, B-V=0.25) counterpart is observed at the pulsar location. Spectra of this counterpart reveal Balmer lines which confirm that the companion is a Helium-core white dwarf. We find that the white dwarf has a temperature of Teff=8060+-150 K. Unfortunately, the spectra are of insufficient quality to put a strong constraint on the surface gravity, although the best fit is for low log g and hence low mass (~0.2 Msun), as expected. We compare predicted white dwarf cooling ages with the characteristic age of the pulsar and find similar values for white dwarf masses of 0.19 to 0.3 Msun. These masses would imply a distance of 2.5 to 4 kpc to the system. The spectroscopic observations also enable us to estimate the mass ratio between the white dwarf and the pulsar. We find q=7.5+-2.4, which is consistent with the current knowledge of white dwarf companions to millisecond pulsars.

We present radio, optical, ultraviolet, and X-ray observations of the isolated, thermally-emitting neutron star RX J0720.4-3125 using the Parkes radio telescope, the Very Large Array, the Hubble Space Telescope, and the Chandra X-ray Observatory. From these data we show that the optical/UV spectrum of RX J0720.4-3125 is not well fit by a Rayleigh-Jeans tail as previously thought, but is instead best fit by either a single non-thermal power-law or a combination of a Rayleigh-Jeans tail and a non-thermal power-law. Taken together with the X-ray spectrum, we find the best model for RX J0720.4-3125 to be two blackbodies plus a power-law, with the cool blackbody implying a radius of 11-13 km at an assumed distance of 300 pc. This is similar to many middle aged (10^{5-6} yr) radio pulsars such as PSR B0656+14, evidence supporting the hypothesis that RX J0720.4-3125 is likely to be an off-beam radio pulsar. The radio data limit the flux at 1.4 GHz to be <0.24 mJy, or a luminosity limit of 4*pi*d^2*F < 3e25*d_300^2 ergs/s, and we see no sign of extended nebulosity, consistent with expectations for a pulsar like RX J0720.4-3125.

We report on our findings of the bright, pulsating, helium atmosphere white dwarf GD 358, based on time-resolved optical spectrophotometry. We identify 5 real pulsation modes and at least 6 combination modes at frequencies consistent with those found in previous observations. The measured Doppler shifts from our spectra show variations with amplitudes of up to 5.5 km/s at the frequencies inferred from the flux variations. We conclude that these are variations in the line-of-sight velocities associated with the pulsational motion. We use the observed flux and velocity amplitudes and phases to test theoretical predictions within the convective driving framework, and compare these with similar observations of the hydrogen atmosphere white dwarf pulsators (DAVs). The wavelength dependence of the fractional pulsation amplitudes (chromatic amplitudes) allows us to conclude that all five real modes share the same spherical degree, most likely, l=1. This is consistent with previous identifications based solely on photometry. We find that a high signal-to-noise mean spectrum on its own is not enough to determine the atmospheric parameters and that there are small but significant discrepancies between the observations and model atmospheres. The source of these remains to be identified. While we infer T_eff=24kK and log g~8.0 from the mean spectrum, the chromatic amplitudes, which are a measure of the derivative of the flux with respect to the temperature, unambiguously favour a higher effective temperature, 27kK, which is more in line with independent determinations from ultra-violet spectra.

We present a high-resolution echelle spectrum taken with the Very Large Telescope and analyse archival Hubble Space Telescope photometry of the recently identified companion of the eclipsing millisecond radio pulsar PSR J1740-5340 in the globular cluster NGC 6397. From the spectrum, we show that the companion is metal poor, as expected for a member of NGC 6397. Using synthetic photometry, and assuming a true distance modulus of 12.13 +/- 0.15 mag and a colour excess of E(B-V) = 0.179, we derive a radius of 1.60 +/- 0.17 solar radii and an effective temperature of 5410 +/- 50 K, implying a luminosity of 2.0 +/- 0.4 solar luminosities (1 sigma errors). These properties make it similar to the so-called `sub-subgiants' in the old open cluster M67 and `red stragglers' in 47 Tuc, which have luminosities comparable to those of turn-off stars, but cooler temperatures and larger radii. The light curve of the companion is well described by ellipsoidal variations, and despite the incomplete (~60%) phase coverage, we are able to derive good constraints on a number of the system parameters. In particular, for the inclination we find a 2 sigma lower limit of i > 48 deg (i > 46 deg at 3 sigma). Assuming a pulsar mass of 1.2 M_MSP < 2.4 solar masses, this implies a companion mass in the range 0.14 < M_comp < 0.38 solar masses. Combined with the photometric constraint, we find a best fit for i ~ 50 deg and M_comp ~ 0.3 solar masses. We infer a Roche lobe filling factor by radius of ~97%. [Abstract shortened]

Using a very deep observation with HST/STIS, we have searched for an optical counterpart to the nearby radio-quiet isolated neutron star RX J1308.6+2127 (RBS 1223). We have identified a single object in the 90% Chandra error circle that we believe to be the optical counterpart. This object has $m_{50CCD}=28.56\pm0.13$ mag, which translates approximately to an unabsorbed flux of $F_{\lambda}=(1.7 \pm 0.3)e-20$ ergs/s/cm^2/A at 5150 A or an X-ray-to-optical flux ratio of $log(f_X/f_opt)=4.9$. This flux is a factor of $\approx 5$ above the extrapolation of the black-body fit to the X-ray spectrum, consistent with the optical spectra of other isolated neutron stars. Without color information we cannot conclude that this source is indeed the counterpart of RX J1308.6+2127. If not, then the counterpart must have $m_{50CCD} > 29.6$ mag, corresponding to a flux that is barely consistent with the extrapolation of the black-body fit to the X-ray spectrum.

RX J0720.4-3125 is the third brightest neutron star in the soft X-ray sky and has been a source of mystery since its discovery, as its long 8-s period separates it from the population of typical radio pulsars. Three models were proposed for this source: a neutron star accreting from the interstellar medium, an off-beam radio pulsar, or an old, cooling magnetar. Using data from Chandra, ROSAT, and BeppoSAX we are able to place an upper limit to the period derivative, $|\dot{P}| < 3.6\times10^{-13}{s s}^{-1}$ (3-$\sigma$). While our upper limit on $\dot P$ allows for the accretion model, this model is increasingly untenable for another similar but better studied neutron star, RX J1856.5-3754, and we therefore consider the accretion model unlikely for RX J0720.4-3125. We constrain the initial magnetic field of RX J0720.4-3125 to be $\lsim 10^{14}$ G based on cooling models, suggesting that it is not and never was a magnetar, but is instead middle-aged neutron star. We propose that it is either a long-period high-magnetic field pulsar with $\dot P\sim 10^{-13}{s s}^{-1}$ similar to PSR J1814-1744, or a neutron star born with an initial period of $\approx 8.3$ s and $\dot P\sim 10^{-15}{s s}^{-1}$. The proximity of RX J0720.4-3125 is strongly suggestive of a large population of such objects; if so, radio pulsar surveys must have missed many of these sources.

The location of the red edge of the ZZ Ceti instability strip is defined observationally as being the lowest temperature for which a white dwarf with a H-rich atmosphere (DA) is known to exhibit periodic brightness variations. Whether this cut-off in flux variations is actually due to a cessation of pulsation or merely due to the attenuation of any variations by the convection zone, rendering them invisible, is not clear. The latter is a theoretical possibility because with decreasing effective temperature, the emergent flux variations become an ever smaller fraction of the amplitude of the flux variations in the interior. In contrast to the flux variations, the visibility of the velocity variations associated with the pulsations is not thought to be similarly affected. Thus, models imply that were it still pulsating, a white dwarf just below the observed red edge should show velocity variations. In order to test this possibility, we used time-resolved spectra of three DA white dwarfs that do not show photometric variability, but which have derived temperatures only slightly lower than the coolest ZZ Ceti variables. We find that none of our three targets show significant periodic velocity variations, and set 95% confidence limits on amplitudes of 3.0, 5.2, and 8.8km/s. Thus, for two out of our three objects, we can rule out velocity variations as large as 5.4km/s observed for the strongest mode in the cool white dwarf pulsator ZZ Psc. In order to verify our procedures, we also examined similar data of a known ZZ Ceti, HL Tau 76. Applying external information from the light curve, we detect significant velocity variations for this object with amplitudes of up to 4km/s, making it only the third DAV for which velocity variations have been measured. (abridged)

We present a detailed analysis of time-resolved optical spectra of the ZZ Ceti white dwarf, HS 0507+0434B. Using the wavelength dependence of observed mode amplitudes, we deduce the spherical degree, l, of the modes, most of which have l=1. The presence of a large number of combination frequencies (linear sums or differences of the real modes) enabled us not only to test theoretical predictions but also to indirectly infer spherical and azimuthal degrees of real modes that had no observed splittings. In addition to the above, we measure line-of-sight velocities from our spectra. We find only marginal evidence for periodic modulation associated with the pulsation modes: at the frequency of the strongest mode in the lightcurve, we measure an amplitude of 2.6+/-1.0 km/s, which has a probability of 2% of being due to chance; for the other modes, we find lower values. Our velocity amplitudes and upper limits are smaller by a factor of two compared to the amplitudes found in ZZ Psc. We find that this is consistent with expectations based on the position of HS 0507+0434B in the instability strip. Combining all the available information from data such as ours is a first step towards constraining atmospheric properties in a convectionally unstable environment from an observational perspective.

RX J1856.5-3754, a bright soft X-ray source believed to be the nearest thermally emitting neutron star, has commanded and continues to command intense interest from X-ray missions. One of the main goals is to determine the radius of this neutron star. An integral part of the determination is an accurate parallax. Walter (2001) analyzed Hubble Space Telescope (HST) data and derived a parallax, pi=16.5+/-2.3 mas. Combining this distance with the angular radius derived from blackbody fits to observations of RX J1856.5-3754 with ROSAT, EUVE, HST, Pons et al. (2001) derived an observed radius ("radiation radius"), R_{infty}= 7 km. This value is smaller than the radii calculated from all proposed equations-of-state (EOS) of dense baryonic matter (Haensel 2001). Here, we have analyzed the same HST data and find pi=7+/-2 mas. We have verified our result using a number of different, independent techniques, and find the result to be robust. Adopting our parallax, the radius of RX J1856.5-3754 is R_{infty}=15+/-6 km. This radius falls squarely in the range of radii, 12--16 km, expected from calculations of neutron star structure for different equations of state. With additional HST observations, the parallax estimate can be improved to the point by which the inferred radius can constrain the choice of EOS.

We report on an approximately twelve hour long X-ray flare from the low-mass X-ray binary KS 1731-260. The flare has a rise time of less than 13 min and declines exponentially with a decay time of 2.7 hours. The flare emission is well described by black-body radiation with peak temperature of 2.4 keV. The total energy release from the event is 10^{42} erg (for an assumed distance of 7 kpc). The flare has all the characteristics of thermo-nuclear X-ray bursts (so-called type I X-ray bursts), except for its very long duration and therefore large energy release (factor of 1500-4000 longer and 250-425 more energy than normal type I X-ray bursts from this source). The flare is preceded by a short and weak X-ray burst, possibly of type I. Days to weeks before the flare, type I X-ray bursts were seen at a rate of ~3 per day. However, after the flare type I X-ray bursting ceased for at least a month, suggesting that the superburst affected the type I bursting behaviour. The persistent emission is not significantly different during the non-bursting period. We compare the characteristics of this event with similar long X-ray flares, so-called superbursts, seen in other sources (4U 1735-44, 4U 1820-30, 4U 1636-53, Ser X-1, GX 3+1). The event seen from KS 1731-260 is the longest reported so far. We discuss two possible mechanisms that might cause these superbursts, unstable carbon burning (as proposed recently) and electron capture by protons with subsequent capture of the resulting neutrons by heavy nuclei.

After 3 years of quiescence, the globular cluster NGC 6440 exhibited a bright transient X-ray source turning on in August 2001, as noted with the RXTE All-Sky Monitor. We carried out a short target of opportunity observation with the Chandra X-ray Observatory and are able to associate the transient with the brightest of 24 X-ray sources detected during quiescence in July 2000 with Chandra. Furthermore, we securely identify the optical counterpart and determine that the 1998 X-ray outburst in NGC 6440 was from the same object. This is the first time that an optical counterpart to a transient in a globular cluster is securely identified. Since the transient is a type I X-ray burster, it is established that the compact accretor is a neutron star. Thus, this transient provides an ideal case to study the quiescent emission in the optical and X-ray of a transiently accreting neutron star while knowing the distance and reddening accurately. One model that fits the quiescent spectrum is an absorbed power law plus neutron star hydrogen atmosphere model. We find an intrinsic neutron star radius of 17_{-12}^{+31} km and an unabsorbed bolometric luminosity for the neutron star atmosphere of (2.1+/-0.8)E33 erg/s which is consistent with predictions for a cooling neutron star.

I discuss constraints on the equation of state of matter at supra-nuclear densities that can be derived from observations of neutron stars. I focus on recent work on Vela X-1, which may well be substantially more massive than the canonical 1.4 solar masses, and on the prospects offered by the `isolated' or `thermally-emitting' neutron stars. I present a re-analysis of the LETG observations of the brightest of the isolated neutron stars, RX J1856.5-3754, and find that the spectrum shows, unlike all other neutron stars observed so far, a significant absorption feature, near 39 Angstrom.

We present HST/STIS, Keck J- and K_s-band, and VLA 332-MHz and 1400-MHz images of the region around the Soft Gamma-ray Repeater SGR 1900+14. No non-stellar sources were detected at the position of SGR 1900+14, giving 3-sigma limits of S_{332 MHz} < 6 mJy, S_{1.4 GHz} < 1.5 mJy, m_{50CCD} > 29.0 mag, J > 22.8 mag, K_s > 20.8 mag (point sources), and S_{332 MHz} < 6.1 mJy arcmin^{-2}, S_{1.4 GHz} < 6.2 mJy arcmin^{-2} (extended emission). Given the very high extinction, the STIS and J-band non-detections are not constraining, but the K_s-band limit rules out basic accretion disk models for SGR 1900+14, and is almost comparable in depth to the K_s-band detection for the anomalous X-ray pulsar 4U 0142+61. Finally, we report the detection in this field of three new candidate supernova remnants, SNRs G043.5+00.6, G042.0-00.1, and G041.5+00.4.

We present near-infrared and optical observations of the field of the Anomalous X-ray Pulsar 1E 2259+58.6 taken with the Keck telescope. We derive a subarcsecond Chandra position and tie it to our optical reference frame using other stars in the field. We find a very faint source, Ks = 21.7\pm0.2 mag, with a position coincident with the Chandra position. We argue that this is the counterpart. In the J, I, and R bands, we derive (2 sigma) limits of 23.8, 25.6 and 26.4 mag, respectively. As with 4U 0142+61, for which a counterpart has previously been found, our results are inconsistent with models in which the source is powered by accretion from a disk, but may be consistent with the magnetar model.

We present spectroscopy and H alpha imaging of a faint nebula surrounding the X-ray bright, nearby neutron star RX J1856.5-3754. The nebula shows no strong lines other than the Hydrogen Balmer lines and has a cometary-like morphology, with the apex being approximately 1 arcsec ahead of the neutron star, and the tail extending up to at least 25 arcsec behind it. We find that the current observations can be satisfactorily accounted for by two different models. In the first, the nebula is similar to ``Balmer-dominated'' cometary nebulae seen around several radio pulsars, and is due to a bow shock in the ambient gas arising from the supersonic motion of a neutron star with a relativistic wind. In the second model, the nebula is an ionisation nebula, but of a type not observed before (though expected to exist), in which the ionisation and heating are very rapid compared to recombination and cooling. We use a detailed hydrodynamical simulation to confirm that a shape like that observed can be reproduced. Independent of the model, we find that RX J1856.5-3754 is not kept hot by accretion. If it is young and cooling, the lack of pulsations at X-ray wavelengths is puzzling. Using phenomenological arguments, we suggest that RX J1856.5-3754 may have a relatively weak, few 10^11 G, magnetic field. If so, it would be ironic that the two brightest nearby neutron stars, RX J1856.5-3754 and RX J0720.4-3125, may well represent the extreme ends of the neutron star magnetic field distribution, one a weak field neutron star and another a magnetar. [Abridged]

We present spectroscopy and imaging with the Very Large Telescope (VLT) of the neutron star RX J1856.5-2754. Little is known about the nature of this source other than that it is a nearby hot neutron star. Our VLT spectrum does not show any strong emission or absorption features. With considerable care to photometric calibration, we obtain photometric measurements over the optical and ultra-violet (UV) using our VLT observations and a detailed analysis of archival {\em Hubble Space Telescope} data. We find that the entire optical to UV spectral energy distribution is well described by a slightly reddened Rayleigh-Jeans tail ($f_\lambda=(2.96\pm0.06)\times10^{-19} (\lambda/5000\un{\AA})^{-4} 10^{-0.4(0.12\pm0.05)(A_\lambda/A_V-1.138)} \un{erg}\un{s^{-1}}\un{cm^{-2}}\un{\AA^{-1}}$, where $A_\lambda/A_V$ is the reddening curve; implied $V=25.58\pm0.02$). The reddening is consistent with the interstellar absorption inferred from X-ray spectroscopy. The simplest explanation for this Rayleigh-Jeans emission is that the optical-UV radiation arises from thermal emission from the surface of the neutron star. The high degree to which the data conform to the Rayleigh-Jeans tail significantly limits contributions from other sources of emission. In particular, our observations are inconsistent with the presence of an accretion disk and also strongly constrain the amount of magnetospheric emission from this enigmatic neutron star.

We measured the radial-velocity curve of HD77581, the B-supergiant companion of the X-ray pulsar Vela X-1, using 183 high-resolution optical spectra obtained in a nine-month campaign. We derive radial-velocity amplitudes for different lines and wavelength regions, and find all are consistent with each other, as well as with values found in previous analyses. We show that one apparent exception, an anomalously low value derived from ultra-violet spectra obtained with the IUE, was due to an error in the analysis procedures. We re-analyse all IUE spectra, and combine the resulting velocities with the ones derived from the new optical spectra presented here, as well as those derived from optical spectra published earlier. As in previous analyses, the radial velocities show strong deviations from those expected for a pure Keplerian orbit. The deviations likely are related to the pronounced line-profile variations seen in our spectra. It turns out that systematic deviations as a function of orbital phase are present. Our best estimate of the radial-velocity amplitude, Kopt = 21.7 +/- 1.6 km/s, has an uncertainty not much reduced to that found in previous analyses, in which the systematic deviations had not been taken into account. Combining our velocity amplitude with the accurate orbital elements of the X-ray pulsar, we infer M_ns sin^3i = 1.78 +/- 0.15 Msun.

We present K band photometry of the region including the radio error box of the soft X-ray transient 4U 1630-47 during its 1998 outburst. We detect a variable source at K=16.1 mag located inside the radio error circle which we identify as the counterpart to the X-ray source. We discuss the properties of the source, and conclude that it is most likely a black-hole X-ray binary similar to 4U 1543-47, GRO J1655-40 or SAX J1819.3-2525, containing a relatively early-type secondary.

Soft Gamma-ray Repeaters (SGRs) are among the most enigmatic sources known today. Exhibiting huge X- and Gamma-ray bursts and flares, as well as soft quiescent X-ray emission, their energy source remains a mystery. Just as mysterious are the Anomalous X-ray pulsars (AXPs), which share many of the same characteristics. Thanks to recent Chandra observations, SGR 0526-66, the first SGR, now appears to be a transition object bridging the two classes, and therefore observations of it have implications for both SGRs and AXPs. The two most popular current models for their persistent emission are accretion of a fossil disk or decay of an enormous (~10^15 G) magnetic field in a magnetar. We show how deep optical observations of SGR 0526-66, the only SGR with small enough optical extinction for meaningful observations, show no evidence of an optical counterpart. These observation place strong new constraints on both accretion disk and magnetar models, and suggest that the spectral energy distribution may peak in the hard-UV. Almost all accretion disks are excluded by the optical data, and a magnetar would require a ~10^15-10^16 G field.

The energy source of the anomalous X-ray pulsars is not well understood, hence their designation as anomalous. Unlike binary X-ray pulsars, no companions are seen, so the energy cannot be supplied by accretion of matter from a companion star. The loss of rotational energy, which powers radio pulsars, is insufficient to power AXPs. Two models are generally considered: accretion from a large disk left over from the birth process, or decay of a very strong magnetic field (10^15 G) associated with a 'magnetar'. The lack of counterparts at other wavelengths has hampered progress in our understanding of these objects. Here, we present deep optical observations of the field around 4U 0142+61, which is the brightest AXP in X-rays. We find an object with peculiar optical colours at the position of the X-ray source, and argue that it is the optical counterpart. The optical emission is too faint to admit the presence of a large accretion disk, but may be consistent with magnetospheric emission from a magnetar.

The BeppoSAX Wide Field Cameras detected a transient in NGC6440 on 1998 Aug 22. ROSAT did not detect this source on 1998 Sep 8, indicating that the X-ray flux had decreased by a factor 400 at least, and/or that the X-ray spectrum had become appreciably softer. Analysis of archival ROSAT HRI data reveals two sources associated with NGC6440; one of these may be the transient in quiescence. We have also obtained B and R images of NGC6440 on 1998 Aug 26 and again on 1999 July 15, when the transient had returned to quiescence. Subtraction of these images reveals one object in the core which was brighter in B, but not in R, during the X-ray outburst. We estimate B=22.7 and (B-R)o<0 on 1998 Aug 26, which makes it a viable candidate counterpart.

We present Keck R and I band images of the field of the anomalous X-ray pulsar 1E 2259+58.6. We derive an improved X-ray position from archival ROSAT HRI observations by correcting for systematic (boresight) errors. Within the corresponding error circle, no object is found on the Keck images, down to limiting magnitudes R = 25.7 and I = 24.3. We discuss the constraints imposed by these limits, and conclude that it is unlikely that 1E 2259+58.6 is powered by accretion from a disk, irrespective of whether it is in a binary or not, unless the binary is extremely compact.

I review attempts made to determine the properties of neutron stars. The focus is on the maximum mass that a neutron star can have, or, conversely, the minimum mass required for the formation of a black hole. There appears to be only one neutron star for which there is strong evidence that its mass is above the canonical 1.4 solar masses, viz., Vela X-1, for which a mass close to 1.9 solar masses is found. Prospects for progress appear brightest for studies of systems in which the neutron star should have accreted substantial amounts of matter.

We report on Keck and {\em Hubble Space Telescope} observations of the binary millisecond pulsar PSR B1855+09. We detect its white-dwarf companion and measure $\mv=25.90\pm0.12$ and $\mi=24.19\pm0.11$ (Vega system). From the reddening-corrected color, $(\mv-\mi)_0=1.06\pm0.21$, we infer a temperature $\Teff=4800\pm800 $K. The white-dwarf mass is known accurately from measurements of the Shapiro delay of the pulsar signal, $\Mcomp=0.258^{+0.028}_{-0.016} \Msun$. Hence, given a cooling model, one can use the measured temperature to determine the cooling age. The main uncertainty in the cooling models for such low-mass white dwarfs is the amount of residual nuclear burning, which is set by the thickness of the hydrogen layer surrounding the helium core. From the properties of similar systems, it has been inferred that helium white dwarfs form with thick hydrogen layers, with mass $\simgt3\times10^{-3} \Msun$, which leads to significant additional heating. This is consistent with expectations from simple evolutionary models of the preceding binary evolution. For PSR B1855+09, though, such models lead to a cooling age of $\sim10 $Gyr, which is twice the spin-down age of the pulsar. It could be that the spin-down age were incorrect, which would call the standard vacuum dipole braking model into question. For two other pulsar companions, however, ages well over 10 Gyr are inferred, indicating that the problem may lie with the cooling models. There is no age discrepancy for models in which the white dwarfs are formed with thinner hydrogen layers ($\simlt3\times10^{-4} \Msun$).

We present time-resolved spectrophotometry of the pulsating DA white dwarf G 29-38. As in previous broad-band photometry, the light curve shows the presence of a large number of periodicities. Many of these are combination frequencies, i.e., periodicities occurring at frequencies that are sums or differences of frequencies of stronger, real modes. We identify at least six real modes, and at least five combination frequencies. We measure line-of-sight velocities for our spectra and detect periodic variations at the frequencies of five of the six real modes, with amplitudes of up to 5 km/s. We argue that these variations reflect the horizontal surface motion associated with the g-mode pulsations. No velocity signals are detected at any of the combination frequencies, confirming that the flux variations at these frequencies do not reflect physical pulsation, but rather mixing of frequencies due to a non-linear transformation in the outer layers of the star. We discuss the amplitude ratios and phase differences found for the velocity and light variations, as well as those found for the real modes and their combination frequencies, both in a model-independent way and in the context of models based on the convective-driving mechanism. In a companion paper, we use the wavelength dependence of the amplitudes of the modes to infer their spherical degree.

We have used time-resolved spectroscopy to measure the colour dependence of pulsation amplitudes in the DAV white dwarf G 29-38. Model atmospheres predict that mode amplitudes should change with wavelength in a manner that depends on the spherical harmonic degree l of the mode. This dependence arises from the convolution of mode geometry with wavelength-dependent limb darkening. Our analysis of the six largest normal modes detected in Keck observations of G 29-38 reveals one mode with a colour dependence different from the other five, permitting us to identify the l value of all six modes and to test the model predictions. The Keck observations also show pulsation amplitudes that are unexpectedly asymmetric within absorption lines. We show that these asymmetries arise from surface motions associated with the non-radial pulsations (which are discussed in detail in a companion paper). By incorporating surface velocity fields into line profile calculations, we are able to produce models that more closely resemble the observations.

We report on Keck and Hubble Space Telescope optical observations of the eclipsing binary pulsar system PSR B1718-19, in the direction of the globular cluster NGC 6342. These reveal a faint star ($m_{\rm F702W}=25.21\pm0.07$; Vega system) within the pulsar's 0\farcs5 radius positional error circle. This may be the companion. If it is a main-sequence star in the cluster, it has radius $\rcomp\simeq0.3 \rsun$, temperature $\teff\simeq3600 $K, and mass $\mcomp\simeq0.3 \msun$. In many formation models, however, the pulsar (spun up by accretion or newly formed) and its companion are initially in an eccentric orbit. If so, for tidal circularization to have produced the present-day highly circular orbit, a large stellar radius is required, i.e., the star must be bloated. Using constraints on the radius and temperature from the Roche and Hayashi limits, we infer from our observations that $\rcomp\simlt0.44 \rsun$ and $\teff\simgt3300 $K. Even for the largest radii, the required efficiency of tidal dissipation is larger than expected for some prescriptions.

Pulsars in close, eccentric binary systems are usually assumed to have another neutron star as a companion. These double neutron star binaries have proven to be the best laboratories for experimental General Relativity and are the most secure candidates for gravitational wave interferometers. We present deep B, V, and R images of the field containing the eccentric binary pulsar system PSR B2303+46. We find a faint, blue object (B=26.60+/-0.09; (B-R)_0=-0.4+/-0.2) coincident with the timing position. We suggest this object is the optical counterpart to the PSR B2303+46 system. The counterpart is too bright to reflect emission from the pulsar or a neutron star companion. Most likely, the companion of PSR B2303+46 is not a neutron star but a massive white dwarf. We show that the observations are consistent with a hot white dwarf companion (T_eff>~5*10^4 K) with cooling age equal the characteristic age of the pulsar (t_cool~30 Myr) and mass within the range set by timing observations and the Chandrasekhar mass (1.2<M_C<1.4 M_sun). Given the eccentric orbit, the white dwarf must have formed before the neutron star, from what was originally the more massive star in the binary. Due to mass transfer, the originally less massive star could become sufficiently massive to end its life in a supernova explosion and form the radio pulsar. We constrain the mass of the pulsar to be in the range 1.24<M_psr<1.44 M_sun.

For the second time in 27 years a bright transient X-ray source has been detected coincident with the globular cluster NGC 6440. It was found to be active in August, 1998, with the Wide Field Camera and the narrow field instruments on the BeppoSAX spacecraft, and with the All-Sky Monitor and the Proportional Counter Array on the RossiXTE spacecraft. Four X-ray bursts were detected, at least one of which shows the characteristics of a thermonuclear flash on a neutron star, in analogy with some ~20 optically identified low-mass X-ray binaries. The broad-band spectrum is hard as is common among low-mass X-ray binaries of lower luminosity (>~10^37 erg/s) and can be explained by a Comptonized model. During the burst the >30 keV emission brightened, consistent with part of the burst emission being Compton up scattered within ~10^11 cm.

Enigmatic transitions between spin-up and spin-down have been observed in several X-ray pulsars accreting matter via an accretion disk. In these transitions, the torque changes sign but remains at nearly the same magnitude. It has been noted previously that alternating prograde and retrograde disk flows would explain many features of the torque reversals, although it has been unclear how a stable retrograde disk could be formed. We suggest that the reversals may be related to the disk at times being warped to such an extent that the inner region becomes tilted by more than 90 degrees. This region would thus become retrograde, leading to a negative torque. Accretion disk models can show such behavior, if account is taken of a warping instability due to irradiation. The resulting `flip-overs' of the inner parts of the disk can reproduce most characteristics of the observations, although it remains unclear what sets the timescale on which the phenomenon occurs. If this model were correct, it would have a number of ramifications, for instance that in the spin-down state the X-ray source would mostly be observed through the accretion disk.

RX J0720.4-3125 is an unidentified bright soft X-ray source which shows pulsations at a 8.39 s period and has a thermal spectrum. We present deep B and R band images of its X-ray localization. We find one possible counterpart in the X-ray error box, with magnitudes B=26.6+/-0.2 and R=26.9+/-0.3. The very high X-ray to optical flux ratio confirms that this object is an isolated neutron star. We discuss possible models and conclude that only two are consistent with the data and at the same time are able to draw from a large enough population to make finding one nearby likely. In our opinion the second criterion provides a stringent constraint but appears to have been ignored so far. The first model, suggested earlier, is that RX J0720.4-3125 is a weakly magnetized neutron star accreting from the interstellar medium. The second is that it is a relatively young, highly magnetized neutron star, a ``magnetar'', which is kept hot by magnetic field decay.

We present a spectrum of the cool (T_eff = 900 K) brown dwarf Gliese 229B. This spectrum, with a relatively high signal-to-noise ratio per spectral resolution element (> 30), spans the wavelength range from 0.837 microns to 5.0 microns. We identify a total of four different major methane absorption features, including the fundamental band at 3.3 microns, at least four steam bands, and two neutral cesium features. We confirm the recent detection of carbon monoxide (CO) in excess of what is predicted by thermochemical equilibrium calculations. Carbon is primarily involved in a chemical balance between methane and CO at the temperatures and pressures present in the outer parts of a brown dwarf. At lower temperatures, the balance favors methane, while in the deeper, hotter regions, the reaction reverses to convert methane into CO. The presence of CO in the observable part of the atmosphere is therefore a sensitive indicator of vertical flows. The high signal-to-noise ratio in the 1 to 2.5 microns region permits us to place constraints on the quantity of dust in the atmosphere of the brown dwarf. We are unable to reconcile the observed spectrum with synthetic spectra that include the presences of dust. The presence of CO but lack of dust may be a clue to the location of the boundaries of the outer convective region of the atmosphere: The lack of dust may mean that it is not being conveyed into the photosphere by convection, or that it exists in patchy clouds. If the dust is not in clouds, but rather sits below the outer convective region, we estimate that the boundary between outer convective and inner radiative layers is between 1250 K and 1600 K, in agreement with recent models.

We present infrared spectroscopy of the low-mass X-ray binary GX 1+4/V2116 Oph. This symbiotic binary consists of a 2-min accretion-powered pulsar and an M5 III red giant. A strong He I 1.083 micron emission line with a pronounced P Cygni profile was observed. From the blue edge of this feature, we infer an outflow velocity of 250(50) km/s. This is an order of magnitude faster than a typical red giant wind, and we suggest that radiation from the accretion disk or the neutron star may contribute to the acceleration of the outflow. We infer a wind mass loss rate of around 10^-6 Msun/yr. Accretion from such a strong stellar wind provides a plausible alternative to Roche lobe overflow for supplying the accretion disk which powers the X-ray source. The H I Paschen beta and He I 1.083 micron lines showed no evidence for the dramatic changes previously reported in some optical lines, and no evidence for pulsations at the 2-min pulsar period.

I introduce the two classes of pulsar, white-dwarf binaries, and describe for each what we have learned from a specific system, PSR J1012+5307 and PSR B0655+64, respectively, summarising what has been done, presenting new results, and discussing what the future may hold. Briefly, for the companion of PSR J1012+5307 we find a DA spectrum, and infer a mass of about 0.16Msun, the lowest among all spectroscopically identified white dwarfs. Combined with a radial-velocity orbit, a neutron-star mass between 1.5 and 3.2Msun (95% conf.) is derived. The companion of PSR B0655+64 shows strong Swan C2 bands, i.e., it is a DQ star. Unlike anything reported for other DQs, however, it shows variations in strength of the bands by a factor two. Most likely, the variations are periodic, with a period of about 9.7h. This is substantially shorter than the 1-day orbital period, which can likely be understood in terms of its past evolution.

I review what properties of the white-dwarf companions of recycled pulsars can be inferred from optical observations, and discuss how these can help us understand the characteristics and evolution of these binaries. I focus on spectroscopic observations, describing results obtained recently, and looking forward to what may come.

We report on spectroscopy of the white-dwarf companion of the millisecond radio pulsar PSR J1012+5307. We find strong Balmer absorption lines, as would be expected for a cool DA white dwarf. The profiles are much narrower than usual, however, and lines are seen up to H12, indicating that the companion has a low gravity and hence a low mass. This is consistent with the expectation---based on evolutionary considerations and on the mass function---that it is a low-mass white dwarf with a helium core. By comparing the spectra to model atmospheres, we derive an effective temperature $T_{\rm{}eff}=8550\pm25\,$K and a surface gravity $\log{}g=6.75\pm0.07$ (cgs units). Using the Hamada-Salpeter mass-radius relation for helium white dwarfs, with an approximate correction for finite-temperature effects, we infer a mass $\mwd=0.16\pm0.02\,\msun$. This is the lowest mass among all spectroscopically identified white dwarfs. We determine radial velocities from our spectra, and find a radial-velocity amplitude of $280\pm15\,\kms$. With the pulsar's radial-velocity amplitude, the mass ratio $\mpsr/\mwd=13.3\pm0.7$. From all constraints, we find that with 95\% confidence $1.5<\mpsr/\msun<3.2$.

We present orbital-phase resolved I and K-band spectroscopy of Cygnus X-3. All spectra show emission lines characteristic of Wolf-Rayet stars of the WN subclass. On time scales longer than about one day, the line strengths show large changes, both in flux and in equivalent width. In addition, the line ratios change, corresponding to a variation in spectral subtype of WN6/7 to WN4/5. We confirm the finding that at times when the emission lines are weak, they shift in wavelength as a function of orbital phase, with maximum blueshift coinciding with infrared and X-ray minimum, and maximum redshift half an orbit later. Furthermore, we confirm the prediction -- made on the basis of previous observations -- that at times when the emission lines are strong, no clear wavelength shifts are observed. We describe a simplified, but detailed model for the system, in which the companion of the X-ray source is a Wolf-Rayet star whose wind is at times ionised by the X-ray source, except for the part in the star's shadow. With this model, the observed spectral variations can be reproduced with only a small number of free parameters. We discuss and verify the ramifications of this model, and find that, in general, the observed properties can be understood. We conclude that Cyg X-3 is a Wolf-Rayet/X-ray binary.

We present spectra of the white-dwarf companions of the radio pulsars 0655+64 and 0820+02. For the latter, we find a spectrum showing strong lines of hydrogen, i.e., that of a DA star. From modelling these lines, the mass of a white dwarf can, in principle, be determined accurately, thus leading to constraints on the evolution of the binary and the mass of the neutron star. Our present spectrum is not of sufficient quality to set a strong limit, but it does indicate that the white dwarf most likely has a low mass. This is consistent with the star being a helium white dwarf, as would be expected from considerations of the mass function and the preceding evolution. From similar considerations, the companion of PSR 0655+64 is expected to be a more massive, carbon-oxygen white dwarf. This is confirmed by our spectra, which show the Swan bands of molecular carbon, making it a DQ star. Unlike what is observed in other DQ stars, the strength of the Swan C2 bands changes drastically, by a factor two in about two hours. We suggests this reflects large-scale surface inhomogeneities which are rotated in and out of the observed hemisphere. If so, this would imply that the white dwarf rotates supersynchronously.

We present new high-resolution, high signal-to-noise optical spectra of HD 77581, the optical counterpart of the X-ray source Vela X-1, and determine radial velocities from these spectra, as well as from high-resolution IUE spectra and from digitized photographic spectra. The measured velocities show strong deviations from a pure Keplerian radial-velocity curve, which are autocorrelated within one night, but not from one night to another. Since lines of different ions exhibit very similar changes in profile, these deviations most likely reflect large-scale motions of the stellar surface akin to non-radial pulsations. A possible cause could be that the varying tidal force exerted by the neutron star in its eccentric orbit excites high-order pulsation modes in the optical star which interfere constructively for short time intervals. The effect of such velocity excursions on the orbital solution is estimated by means of a Monte-Carlo simulation technique. We investigate sources of systematic error, due to, e.g., the tidal deformation of the star, and find, in particular, evidence for a systematic perturbation of the radial velocity near the time of velocity minimum. This possible distortion severely compromises the accuracy of the radial-velocity amplitude, leading to a 95\% confidence range of 18.0--28.2\,\kms. The corresponding 95\% confidence limits of the masses are given by $\MX=1.9_{-0.5}^{+0.7}\,\Msun$ and $\Mopt=23.5_{-1.5}^{+2.2}\,\Msun$.

We analyze the currently available observations of X-ray binaries in a consistent way, to re-determine the masses of the neutron stars in these systems. In particular, our attention is focussed on a realistic and consistent assessment of observational uncertainties and sources of systematic error. Confidence limits for these new mass estimates are generally less optimistic than previously assumed. The available observations, including data on six radio pulsars, do not firmly constrain the equation of state of neutron star matter. In particular, a firm upper mass limit cannot yet be established. An improvement of the accuracy of optical data holds the key to further progress.

We have obtained infrared spectra of the star suggested to be the counterpart of the soft gamma-ray repeater (SGR) 1806-20. We found strong emission lines similar to those seen in the spectra of the rare Luminous Blue Variables and B[e] stars. A He I absorption line is also seen, from which we infer a spectral type O9--B2. This classification, in combination with the minimum distance of $\simgt6$ kpc inferred from its extinction, makes the star one of the most luminous in the Galaxy. We infer that it is a companion to SGR 1806-20, and suggest that the presence of a companion is somehow related to the SGR phenomenon.

The soft gamma-ray repeater (SGR) 1806$-$20 is associated with the center-brightened non-thermal nebula G~10.0$-$0.3, thought to be a plerion. As in other plerions, a steady \Xray\ source, AX~1805.7$-$2025, has been detected coincident with the peak of the nebular radio emission. Vasisht et al.\ have shown that the radio peak has a core-jet appearance, and argue that the core marks the true position of the SGR. At optical wavelengths, we detect three objects in the vicinity of the radio core. Only for the star closest to the core, barely visible in the optical but bright in the infrared ($K=8.4\,$mag.), the reddening is consistent with the high extinction ($A_V\simeq30\,$mag.) that has been inferred for AX~1805.7$-$2025. From the absence of CO band absorption, we infer that the spectral type of this star is earlier than late~G/early~K. The large extinction probably arises in a molecular cloud located at a distance of 6$\,$kpc, which means that the star, just like AX~1805.7$-$2025, is in or behind this cloud. This implies that the star is a supergiant. Since supergiants are rare, a chance coincidence with the compact radio core is very unlikely. To our knowledge, there are only three other examples of luminous stars embedded in non-thermal radio nebulae, SS~433, \mbox{Cir X-1} and G~70.7+1.2. Given this and the low coincidence probability, we suggest that the bright star is physically associated with SGR~1806$-$20, making it the first stellar identification of a high-energy transient.