Our goal in this paper is to derive a carbon-star luminosity function that will eventually be used to determine distances to galaxies at $50$-$60$ Mpc and hence yield a value of the Hubble constant. Cool N-type carbon stars exhibit redder near-infrared colours than oxygen-rich stars. Using Two Micron All Sky Survey near-infrared photometry and the Gaia Data Release 2, we identify carbon stars in the Magellanic Clouds (MC) and the Milky Way (MW). Carbon stars in the MC appear as a distinct horizontal feature in the near-infrared ($(J-K_s)_0$, $M_J$) colour-magnitude diagram. We build a colour selection ($1.4 < (J-K_s)_0 < 2$) and derive the luminosity function of the colour-selected carbon stars. We find the median absolute magnitude and the dispersion, in the J band, for the Large Magellanic Cloud and the Small Magellanic Cloud to be, respectively, ($\bar{M_J} = -6.284 \pm 0.004$, $\sigma = 0.352 \pm 0.005$) and ($\bar{M_J} = -6.160 \pm 0.015$, $\sigma = 0.365 \pm 0.014$). The difference between the MC may be explained by the lower metallicity of the Small Magellanic Cloud, but in any case it provides limits on the type of galaxy whose distance can be determined with this technique. To account for metallicity effects, we developed a composite magnitude, named C, for which the error-weighted mean C magnitude of both the MC are equal. Thanks to the next generation of telescopes (JWST, ELT, TMT), carbon stars could be detected in MC-type galaxies at distances out to $50$-$60$ Mpc. The final goal is to eventually try and improve the measurement of the Hubble constant while exploring the current tensions related to its value.

X-ray telescopes are powerful tools for the study of neutron stars and black holes. Several methods have already been developed (e.g., principal component analysis, nonlinear optimal filtering or high-rate processing method) to analyze the pulses that result from X-rays absorbed in superconducting transition-edge sensors. Our goal is to develop a lightweight, linear filter that will maximize energy and time resolution when X-ray photons are detected by transition-edge sensors. Furthermore, we find the minimal sampling rate that will not degrade the energy and time-resolution of these techniques. Our method is designed for the widest range of photon energies (from $0.1$ keV to $30$ keV). Transition-edge sensors exhibit a non-linear response that becomes more pronounced with increasing photon energy; therefore, we need to treat high-energy photons differently from low-energy photons. In order to retrieve the energy and the arrival time of the photon, we fit simulations of the evolution of the current including the typical noise sources in a sensor with simulated theoretical models. The curve-fitting parameters are interpolated to extract the energy and time resolution. For energies from $0.1$ keV to $30$ keV and with a sampling rate of $195$ kHz, we successfully obtain a $2\sigma$-energy resolution between $1.67$ eV and $6.43$ eV. Those results hold if the sampling rate decreases by a factor two. About time resolution, with a sampling rate of $195$ kHz we get a $2\sigma$-time resolution between $94$ ns and $0.55$ ns for a sensor with the physical parameters as those used in the HOLMES experiment. In order to make this method useful on a larger scale, it will be essential to get a more general description of the noise in a TES, and it will be necessary to develop a robust way to identify pile-up events.

Since the birth of X-ray Astronomy, spectacular advances have been seen in the imaging, spectroscopic and timing studies of the hot and violent X-ray Universe, and further leaps forward are expected in the future. On the other hand, polarimetry is very much lagging behind: after the measurements of the Crab Nebula and Scorpius X-1, obtained by OSO-8 in the 70s, no more observations have been performed in the classical X-ray band, even if some interesting results have been obtained in hard X-rays and in soft gamma-rays. The NASA/ASI mission IXPE, scheduled for the launch in 2021, is going to provide for the first time imaging X-ray polarimetry in the 2-8 keV band thanks to its photoelectric polarimeter, coupled with ~25'' angular resolution X-ray mirrors. Its orders of magnitude improvement in sensitivity with respect to the OSO-8 Bragg polarimeter implies scientifically meaningful polarimetric measurements for at least the brightest specimens of most classes of X-ray sources. In 2027, the Chinese-led mission eXTP should also be launched. In addition to timing and spectroscopic instruments, eXTP will have on board photoelectric polarimeters very similar to those of IXPE, but with a total effective area 2-3 times larger. Building on IXPE results, eXTP will increase the number of sources for which significant polarimetric measurements could be obtained. However, further progresses, such as exploring a broader energy range, considering a larger effective area, improving the angular resolution, and performing wide-field polarization measurements, are needed to reach a mature phase for X-ray polarimetry. In the first part of this White Paper we will discuss a few scientific cases in which a next generation X-ray Polarimetry mission can provide significant advances. In the second part, a possible concept for a medium-class Next Generation X-ray Polarimetry (NGXP) mission will be sketched.

In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional observing locations can be competitive, under certain circumstances, with larger and more expensive facilities located at poorer sites. In addition to research, these facilities enable many different types of educational experiences for wide range of people, from high school students to undergraduates to graduate students to postdocs. Canada should remain committed to partnering with large, international observatories such as CFHT, Gemini, and TMT, but it should also negotiate international agreements and commit funding to expand the use of small and moderate research observatories at domestic and international sites through coordination with the NRC, the Tri-Council, and the Canadian Foundation for Innovation. Both capital and operational costs (with site rental costs allowed) need to be included in support possibilities. CASCA should establish and maintain a small to moderate telescope expression of interest database that would help to facilitate Canadian institutions in organizing consortia, particularly for smaller institutions. The astronomical community should work with the NRC to make existing facilities more accessible to the astronomical community for research. This could involve, for example, automating the Plaskett and/or providing travel funds for supporting classical observing modes. Finally, a small to moderate aperture facility in the Arctic would be a world-class observatory and should be advanced over the next decade.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We have searched the Gaia DR2 catalogue for previously unknown hot white dwarfs in the direction of young open star clusters. The aim of this experiment was to try and extend the initial-final mass relation (IFMR) to somewhat higher masses, potentially providing a tension with the Chandrasekhar limit currently thought to be around 1.38 M$_{\odot}$. We discovered a particularly interesting white dwarf in the direction of the young $\sim$150 Myr old cluster Messier 47 (NGC 2422). All Gaia indicators (proper motion, parallax, location in the Gaia colour-magnitude diagram) suggest that it is a cluster member. Its spectrum, obtained from Gemini South, yields a number of anomalies: it is a DB (helium-rich atmosphere) white dwarf, it has a large magnetic field (2.5 MG), is of high mass ($\sim$1.06 M$_\odot$) and its colours are very peculiar --- particularly the redder ones ($r$, $i$, $z$ and $y$), which suggest that it has a late-type companion. This is the only magnetized, detached binary white dwarf with a non-degenerate companion of any spectral type known in or out of a star cluster. If the white dwarf is a cluster member, as all indicators suggest, its progenitor had a mass just over 6 M$_\odot$. It may, however, be telling an even more interesting story than the one related to the IFMR, one about the origin of stellar magnetic fields, Type I supernovae and gravitational waves from low mass stellar systems.

In this paper, we analyze stellar proper motions in the core of the globular cluster 47 Tucanae to explore the possibility of an intermediate-mass black hole (IMBH) influence on the stellar dynamics. Our use of short-wavelength photometry affords us an exceedingly clear view of stellar motions into the very center of the crowded core, yielding proper motions for $>$50,000 stars in the central 2'. We model the velocity dispersion profile of the cluster using an isotropic Jeans model. The density distribution is taken as a central IMBH point mass added to a combination of King templates. We individually model the general low-mass cluster objects (main sequence/giant stars), as well as the concentrated populations of heavy binary systems and dark stellar remnants. Using unbinned likelihood model fitting, we find that the inclusion of the concentrated populations in our model plays a crucial role in fitting for an IMBH mass. The concentrated binaries and stellar-mass black holes (BHs) produce a sufficient velocity dispersion signal in the core so as to make an IMBH unnecessary to fit the observations. We additionally determine that a stellar-mass BH retention fraction of $\gtrsim 8.5\%$ becomes incompatible with our observed velocities in the core.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

The energy loss of a rotationally powered pulsar is primarily carried away as electromagnetic radiation and a particle wind. Considering that the magnetic field strength of pulsars ranges from about $10^8$ to $10^{15}$ G, one could expect quantum electrodynamics (QED) to play a role in their spin-down, especially for strongly magnetized ones (magnetars). In fact several authors have argued that QED corrections will dominate the spin-down for slowly rotating stars. They called this effect quantum vacuum friction (QVF). However, QVF was originally derived using a problematic self-torque technique, which leads to a dramatic overestimation of this spin-down effect. Here, instead of using QVF, we explicitly calculate the energy loss from rotating neutron stars using the Poynting vector and a model for a particle wind, and we include the QED one-loop corrections. We express the excess emission as QED one-loop corrections to the radiative magnetic moment of a neutron star. We do find a small component of the spin-down luminosity that originates from the vacuum polarization. However, it never exceeds one percent of the classical magnetic dipole radiation in neutron stars for all physically interesting field strengths. Therefore, we find that the radiative corrections of QED are irrelevant in the energetics of neutron-star spin-down.

The detection of brown dwarfs in globular star clusters will allow us to break the degeneracies in age, mass and composition that affect our current models, and therefore to constrain the physics of their atmospheres and interiors. Furthermore, detecting brown dwarfs will help us constrain the properties of the clusters themselves, as they carry information about the clusters' age and dynamics. Great advancements in this field are to be expected in the next ten years, thanks to the extraordinary sensitivity in the infrared of upcoming telescopes like JWST and the ELTs.

The X-ray emission of neutron stars and black holes presents a rich phenomenology that can lead us to a better understanding of their nature and to address more general physics questions: Does general relativity apply in the strong gravity regime? Is spacetime around black holes described by the Kerr metric? This white paper considers how we can investigate these questions by studying reverberation mapping and quasi-periodic oscillations in accreting systems with a combination of high-spectral and high-timing resolution. In the near future, we will be able to study compact objects in the X-rays in a new way: advancements in transition-edge sensors (TES) technology will allow for electron-volt-resolution spectroscopy combined with nanoseconds-precision timing.

The advent of moderately high-resolution X-ray spectroscopy with Chandra and XMM promised to usher in a new age in the study of neutron stars: we thought we would study neutron stars like stars, with resolved absorption spectra revealing their surface chemical composition and physical conditions (e.g. surface gravity, pressure, temperature). Nature, however, did not cooperate in this endeavor, as observations of neutron stars have not revealed verified atomic absorption lines yet. In the near future, advancements in transition-edge sensors (TES) technology will allow for electron-volt-resolution spectroscopy combined with nanoseconds-precision timing. Combining these detectors with collector optics will also us to study neutron stars in much greater detail by achieving high-energy resolution with much larger collecting areas to uncover even weak spectral features over a wide range of the photon energies. Perhaps we will finally be able to study neutron stars like stars.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

Using parallaxes from Gaia DR2, we estimate the distance to the globular clusters 47 Tuc and NGC 362, taking advantage of the background stars in the Small Magellanic Cloud and quasars to account for various parallax systematics. We found the parallax to be dependent on the Gaia DR2 G apparent magnitude for stars with $13<\textrm{G}<18$ , where brighter stars have a lower parallax zero point than fainter stars. The distance to 47 Tuc was found to be $4.45\pm0.01\pm0.12$ kpc, and for NGC 362 $8.54\pm0.20\pm0.44$ kpc with random and systematic errors listed respectively. This is the first time a precise distance measurement directly using parallaxes has been done for either of these two globular clusters.

Radiative corrections of quantum electrodynamics cause a vacuum threaded by magnetic field to be birefringent. This means that radiation of different polarizations travels at different speeds. Even in the strong magnetic fields of astrophysical sources the difference in speed is small; however, it has profound consequences for the extent of polarization expected from strongly magnetized sources. We demonstrate how the birefringence arises from first principles, show how birefringence affects the polarization state of radiation and present recent calculations for the expected polarization from magnetars and X-ray pulsars.

The effect of vacuum birefringence is one of the first predictions of quantum electrodynamics (QED): the presence of a charged Dirac field makes the vacuum birefringent when threaded by magnetic fields. This effect, extremely weak for terrestrial magnetic fields, becomes important for highly magnetized astrophysical objects, such as accreting black holes. In the X-ray regime, the polarization of photons traveling in the magnetosphere of a black hole is not frozen at emission but is changed by the local magnetic field. We show that, for photons traveling along the plane of the disk, where the field is expected to be partially organized, this results in a depolarization of the X-ray radiation. Because the amount of depolarization depends on the strength of the magnetic field, this~effect can provide a way to probe the magnetic field in black-hole accretion disks and to study the role of magnetic fields in astrophysical accretion in general.

In the next decade, x-ray polarimetry will open a new window on the high-energy Universe, as several missions that include an x-ray polarimeter are currently under development. Observations of the polarization of x-rays coming from the accretion disks of stellar-mass and supermassive black holes are among the new polarimeters' major objectives. In this paper, we show that these observations can be affected by the quantum electrodynamic (QED) effect of vacuum birefringence: after an x-ray photon is emitted from the accretion disk, its polarization changes as the photon travels through the accretion disk's magnetosphere, as a result of the vacuum becoming birefringent in presence of a magnetic field. We show that this effect can be important for black holes in the energy band of the upcoming polarimeters, and has to be taken into account in a complete model of the x-ray polarization that we expect to detect from black-hole accretion disks, both for stellar mass and for supermassive black holes. We find that, for a chaotic magnetic field in the disk, QED can significantly decrease the linear polarization fraction of edge-on photons, depending on the spin of the hole and on the strength of the magnetic field. This effect can provide, for the first time, a direct way to probe the magnetic field strength close to the innermost stable orbit of black-hole accretion disks and to study the role of magnetic fields in astrophysical accretion in general.

Radio pulsars found in binary systems with short orbital periods are usually fast spinning as a consequence of recycling via mass transfer from their companion stars; this process is also thought to decrease the magnetic field of the neutron star being recycled. Here, we report on timing observations of the recently discovered binary PSR J1755$-$2550 and find that this pulsar is an exception: with a characteristic age of 2.1 Myr, it is relatively young; furthermore, with a spin period of 315 ms and a surface magnetic field strength at its poles of 0.88$\times$10$^{12}$ G the pulsar shows no sign of having been recycled. Based on its timing and orbital characteristics, the pulsar either has a massive white dwarf (WD) or a neutron star (NS) companion. To distinguish between these two cases, we searched radio observations for a potential recycled pulsar companion and analysed archival optical data for a potential WD companion. Neither work returned conclusive detections. We apply population synthesis modelling and find that both solutions are roughly equally probable. Our population synthesis also predicts a minimum mass of 0.90 M$_{\odot}$ for the companion star to PSR J1755$-$2550 and we simulate the systemic runaway velocities for the resulting WDNS systems which may merge and possibly produce Ca-rich supernovae. Whether PSR J1755$-$2550 hosts a WD or a NS companion star, it is certainly a member of a rare subpopulation of binary radio pulsars.

Cosmology is intrinsically intertwined with questions in fundamental physics. The existence of non-baryonic dark matter requires new physics beyond the Standard Model of elemenatary-particle interactions and Einstein's general relativity, as does the accelerating expansion of the universe. Current tensions between various cosmological measurements may be harbingers of yet more new physics. Progress on understanding dark matter and cosmic acceleration requires long term, high-precision measurements and excellent control of systematics, demanding observational programs that are often outside the discovery/characterization mode that drives many areas of astronomy. We outline potential programs through which the Hubble Space Telescope (HST) could have a major impact on issues in fundamental physics in the coming years. To realize this impact, we suggest the introduction of a "HST Fundamental Physics" observational program that would be subject to a modified proposal and review process.

Multi-epoch observations with ACS and WFC3 on HST provide a unique and comprehensive probe of stellar dynamics within 47 Tucanae. We confront analytic models of the globular cluster with the observed stellar proper motions that probe along the main sequence from just above 0.8 to 0.1M$_\odot$ as well as white dwarfs younger than one gigayear. One field lies just beyond the half-light radius where dynamical models (\eg lowered Maxwellian distributions) make robust predictions for the stellar proper motions. The observed proper motions in this outer field show evidence for anisotropy in the velocity distribution as well as skewness; the latter is evidence of rotation. The measured velocity dispersions and surface brightness distributions agree in detail with a rotating, anisotropic model of the stellar distribution function with mild dependence of the proper-motion dispersion on mass. However, the best fitting models under-predict the rotation and skewness of the stellar velocities. In the second field, centered on the core of the cluster, the mass segregation in proper motion is much stronger. Nevertheless the model developed in the outer field can be extended inward by taking this mass segregation into account in a heuristic fashion. The proper motions of the main-sequence stars yield a mass estimate of the cluster of $1.31 \pm 0.02 \times 10^6 \mathrm{M}_\odot$ at a distance of 4.7 kpc. By comparing the proper motions of a sample of giant and sub-giant stars with the observed radial velocities we estimate the distance to the cluster kinematically to be $4.29 \pm 0.47$ kpc.

Axions and axion-like particles are a leading model for the dark matter in the Universe; therefore, dark matter halos may be boson stars in the process of collapsing. We examine a class of static boson stars with a non-minimal coupling to gravity. We modify the gravitational density of the boson field to be proportional to an arbitrary power of the modulus of the field, introducing a non-standard coupling. We find a class of solutions very similar to Newtonian polytropic stars that we denote "quantum polytropes." These quantum polytropes are supported by a non-local quantum pressure and follow an equation very similar to the Lane-Emden equation for classical polytropes. Furthermore, we derive a simple condition on the exponent of the non-linear gravitational coupling, $\alpha>8/3$, beyond which the equilibrium solutions are unstable.

Using images from the Hubble Space Telescope Advanced Camera for Surveys, we measure the rate of cooling of white dwarfs in the globular cluster 47 Tucanae and compare it to modelled cooling curves. We examine the effects of the outer convective envelope reaching the nearly isothermal degenerate core and the release of latent heat during core crystallisation on the white dwarf cooling rates. For white dwarfs typical of 47 Tuc, the onset of these effects occur at similar times. The latent heat released during crystallisation is a small heat source. In contrast, the heat reservoir of the degenerate core is substantially larger. When the convective envelope reaches the nearly isothermal interior of the white dwarf, the star becomes brighter than it would be in the absence of this effect. Our modelled cooling curves that include this convective coupling closely match the observed luminosity function of the white dwarfs in 47 Tuc.

Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violation, focusing on galileon gravity as an example. For galaxies in-falling towards galaxy clusters, the supermassive black hole can be offset from the galaxy center away from the direction of the cluster. Hence, well resolved images of galaxies around nearby clusters can be used to identify the displaced black hole via the star cluster bound to it. We show that this signal is accessible with imaging surveys, both ongoing ones such as the Dark Energy Survey, and future ground and space based surveys. Already, the observation of the central black hole in M~87 places new constraints on the galileon parameters, which we present here. $\mathcal{O}(1)$ matter couplings are disfavored for a large region of the parameter space. We also find a novel phenomenon whereby the black hole can escape the galaxy completely in less than one billion years.

High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investigate the properties of the required scatterers, and whether the mechanism is consistent with the degree of pulsed emission observed in the tail of the giant flare. We conclude that the mechanism is consistent with current data, although the origin of the scattering population remains an open question. We propose an alternative picture in which the emission is closer to that star and is dominated by synchrotron radiation. The $RHESSI$ observations of the December 2004 flare modestly favor this latter picture. We assess the prospects for the Fermi Gamma-Ray Space Telescope to detect and characterize a similar high energy component in a future giant flare. Such a detection should help to resolve some of the outstanding issues.

We present a model to account for the observed debris disks around young white dwarfs and the presence of metal-lines in their spectra. Stellar evolution models predict that the mass-loss on the AGB will be pulsed; furthermore, observations indicate that the bulk of the mass-loss occurs on the AGB. In this case, if the progenitors of the white dwarfs had remnants of planetary formation like the Sun's Oort cloud or the Kuiper Belt and a planet lying within that cloud or nearby, we find that up to 2% of the planetesimals will fall either into planet-crossing orbits or into chaotic regions after the mass-loss, depending on the location and mass of the planet (from Mars to Neptune). This yields a sufficient mass of comets that can be scattered toward the star, form a debris disk and pollute the atmosphere.

Globular clusters are the oldest conglomerates of stars in our Galaxy and can be useful laboratories to test theories from stellar evolution to cosmology. In this paper, we present a new method to estimate the absolute age of a globular cluster from observations of its brown dwarfs. The transition region between the end of the main sequence and the brown dwarf regime is characterized by a dearth of objects as function of magnitude. The brightest of the cooling brown dwarfs is easily identified by an increase in density in the color-magnitude diagram as you go fainter in magnitudes, and these brightest brown dwarfs get fainter with age. By identifying the brightest brown dwarfs, it is thus possible to determine the age of a globular cluster within a 1 Gyr precision with four-sigma confidence. This new method, which is independent of current methods of age estimation and which does not rely on the knowledge of the cluster's distance from Earth, will become feasible thanks to the high spatial resolution and incredible infrared sensitivity of the James Webb Space Telescope.

We present a new and efficient algorithm for finding point sources in the photon event data stream from the Fermi Gamma-Ray Space Telescope, FermiFAST. The key advantage of FermiFAST is that it constructs a catalogue of potential sources very fast by arranging the photon data in a hierarchical data structure. Using this structure FermiFAST rapidly finds the photons that could have originated from a potential gamma-ray source. It calculates a likehihood ratio for the contribution of the potential source using the angular distribution of the photons within the region of interest. It can find within a few minutes the most significant half of the Fermi Third Point Source catalogue (3FGL) with nearly 80\% purity from the four years of data used to construct the catalogue. If a higher purity sample is desirable, one can achieve a sample that includes the most significant third of the Fermi 3FGL with only five percent of the sources unassociated with Fermi sources. Outside the galaxy plane, all but eight of the 580 FermiFAST detections are associated with 3FGL sources. And of these eight, six yield significant detections of greater than five sigma when a further binned likelihood analysis is performed. This software allows for rapid exploration of the Fermi data, simulation of the source detection to calculate the selection function of various sources and the errors in the obtained parameters of the sources detected.

The recent discovery of gravitational radiation from merging black holes poses a challenge of how to organize the electromagnetic follow-up of gravitational-wave events as well as observed bursts of neutrinos. We propose a technique to select the galaxies that are most likely to host the event given some assumptions of whether the particular event is associated with recent star formation, low metallicity stars or simply proportional to the total stellar mass in the galaxy. We combine data from the 2-MASS Photometric Redshift Galaxy Catalogue with results from galaxy formation simulations to develop observing strategies that potentially reduce the area of sky to search by up to a factor of two relative to an unweighted search of galaxies, and a factor twenty to a search over the entire LIGO localization region.

By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $\eta_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for this technique, the high concentration of stars near its centre substantially weakens the constraints. On the other hand, if a central black hole is found in NGC 205 or one of the other satellite ellipticals of M\,31, substantially better constraints could be obtained. In all cases the constraints could improve dramatically with better astrometry.

Blue stragglers (BSS) are stars whose position in the Color-Magnitude Diagram (CMD) places them above the main sequence turn-off (TO) point of a star cluster. Using data from the core of 47 Tuc in the ultraviolet (UV), we have identified various stellar populations in the CMD, and used their radial distributions to study the evolution and origin of BSS, and obtain a dynamical estimate of the mass of BSS systems. When we separate the BSS into two samples by their magnitude, we find that the bright BSS show a much more centrally concentrated radial distribution and thus higher mass estimate (over twice the TO mass for these BSS systems), suggesting an origin involving triple or multiple stellar systems. In contrast, the faint BSS are less concentrated, with a radial distribution similar to the main sequence (MS) binaries, pointing to the MS binaries as the likely progenitors of these BSS. Putting our data together with available photometric data in the visible and using MESA evolutionary models, we calculate the expected number of stars in each evolutionary stage for the normal evolution of stars and the number of stars coming from the evolution of BSS. The results indicate that BSS have a post-MS evolution comparable to that of a normal star of the same mass and a MS BSS lifetime of about 200-300 Myr. We also find that the excess population of asymptotic giant branch (AGB) stars in 47 Tuc is due to evolved BSS.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

The initial discovery of LIGO on 14 September 2015 was the inspiral merger and ring-down of the black hole binary at a distance of about 500~Mpc or a redshift of about 0.1. The search for electromagnetic counterparts for the inspiral of binary black holes is impeded by coarse initial source localisations and a lack of a compelling model for the counterpart; therefore, rapid electromagnetic follow-up is required to understand the astrophysical context of these sources. Because astrophysical sources of gravitational radiation are likely to reside in galaxies, it would make sense to search first in regions where the LIGO-Virgo probability is large and where the density of galaxies is large as well. Under the assumption that the probability of a gravitational-wave event from a given region of space is proportional to the density of galaxies within the probed volume, one can calculate an improved localisation of the position of the source simply by multiplying the LIGO-Virgo skymap by the density of galaxies in the range of redshifts. We propose using the 2-MASS Photometric Redshift Galaxy Catalogue for this purpose and demonstrate that using it can dramatically reduce the search region for electromagnetic counterparts.

We use the effects of mass segregation on the radial distribution of different stellar populations in the core of 47 Tucanae to find estimates for the masses of stars at different post main sequence evolutionary stages. We take samples of main sequence (MS) stars from the core of 47 Tucanae, at different magnitudes (i.e. different masses), and use the effects of this dynamical process to develop a relation between the radial distance (RD) at which the cumulative distribution reaches the 20th and 50th percentile, and stellar mass. From these relations we estimate the masses of different post MS populations. We find that mass remains constant for stars going through the evolutionary stages between the upper MS up to the horizontal branch (HB). By comparing RDs of the HB stars with stars of lower masses, we can exclude a mass loss greater than 0.09M during the red giant branch (RGB) stage at nearly the 3{\sigma} level. The slightly higher mass estimates for the asymptotic giant branch (AGB) are consistent with the AGB having evolved from somewhat more massive stars. The AGB also exhibits evidence of contamination by more massive stars, possibly blue stragglers (BSS), going through the RGB phase. We do not include the BSS in this paper due to the complexity of these objects, instead, the complete analysis of this population is left for a companion paper. The process to estimate the masses described in this paper are exclusive to the core of 47 Tuc.

We present a robust statistical analysis of the white dwarf cooling sequence in 47 Tucanae. We combine HST UV and optical data in the core of the cluster, Modules for Experiments in Stellar Evolution (MESA) white dwarf cooling models, white dwarf atmosphere models, artificial star tests, and a Markov Chain Monte Carlo (MCMC) sampling method to fit white dwarf cooling models to our data directly. We use a technique known as the unbinned maximum likelihood to fit these models to our data without binning. We use these data to constrain neutrino production and the thickness of the hydrogen layer in these white dwarfs. The data prefer thicker hydrogen layers $(q_\mathrm{H}=3.2\e{-5})$ and we can strongly rule out thin layers $(q_\mathrm{H}=10^{-6})$. The neutrino rates currently in the models are consistent with the data. This analysis does not provide a constraint on the number of neutrino species.

Giant $\gamma$-ray flares comprise the most extreme radiation events observed from magnetars. Developing on (sub)millisecond timescales and generating vast amounts of energy within a fraction of a second, the initial phase of these extraordinary bursts present a significant challenge for candidate trigger mechanisms. Here we assess and critically analyse the linear growth of the relativistic tearing instability in a globally twisted magnetosphere as the trigger mechanism for giant $\gamma$-ray flares. Our main constraints are given by the observed emission timescales, the energy output of the giant flare spike, and inferred dipolar magnetic field strengths. We find that the minimum growth time of the linear mode is comparable to the $e$-folding rise time, i.e. $\sim10^{-1}$ ms. With this result we constrain basic geometric parameters of the current sheet. We also discuss the validity of the presumption that the $e$-folding emission timescale may be equated with the growth time of an MHD instability.

By examining the diffusion of young white dwarfs through the core of the globular cluster 47 Tucanae, we estimate the time when the progenitor star lost the bulk of its mass to become a white dwarf. According to stellar evolution models of the white-dwarf progenitors in 47 Tucanae, we find this epoch to coincide approximately with the star ascending the asymptotic giant branch ($3.0 \pm 8.1$ Myr before the tip of the AGB) and more than ninety million years after the helium flash (with ninety-percent confidence). From the diffusion of the young white dwarfs we can exclude the hypothesis that the bulk of the mass loss occurs on the red-giant branch at the four-sigma level. Furthermore, we find that the radial distribution of horizontal branch stars is consistent with that of the red-giant stars and upper-main-sequence stars and inconsistent with the loss of more than 0.2 solar masses on the red-giant branch at the six-sigma level.

We present Hubble Space Telescope observations of the upper part (T_eff> 10 000 K) of the white dwarf cooling sequence in the globular cluster 47 Tucanae and measure a luminosity function of hot white dwarfs. Comparison with previous determinations from large scale field surveys indicates that the previously determined plateau at high effective temperatures is likely a selection effect, as no such feature is seen in this sample. Comparison with theoretical models suggests that the current estimates of white dwarf neutrino emission (primarily by the plasmon channel) are accurate, and variations are restricted to no more than a factor of two globally, at 95% confidence. We use these constraints to place limits on various proposed exotic emission mechanisms, including a non-zero neutrino magnetic moment, formation of axions, and emission of Kaluza-Klein modes into extra dimensions.

Using images from the Hubble Space Telescope Wide-Field Camera 3, we measure the rate of diffusion of stars through the core of the globular cluster 47 Tucanae using a sample of young white dwarfs identified in these observations. This is the first direct measurement of diffusion due to gravitational relaxation. We find that the diffusion rate $\kappa\approx 10-13$ arcsecond$^2$ Myr$^{-1}$ is consistent with theoretical estimates of the relaxation time in the core of 47 Tucanae of about 70 Myr.

The energy levels of the first few low-lying states of helium and lithium atoms in intense magnetic fields up to $\approx 10^8-10^9$~T are calculated in this study. A pseudospectral method is employed for the computational procedure. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these two- and three-electron systems in a self-consistent manner. The method exploits the natural symmetries of the problem without assumptions of any basis functions for expressing the wave functions of the electrons or the commonly employed adiabatic approximation. It is seen that the results obtained here for a few of the most tightly bound states of each of the atoms, helium and lithium, are in good agreement with findings elsewhere. In this regard, we report new data for two new states of lithium that have not been studied thus far in the literature. It is also seen that the pseudospectral method employed here is considerably more economical, from a computational point of view, than previously employed methods such as a finite-element based approach. The key enabling advantage of the method described here is the short computational times which are on the order of seconds for obtaining accurate results for heliumlike systems.

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci - a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

We present core radii for 54 Milky Way globular clusters determined by fitting King-Michie models to cumulative projected star count distributions. We find that fitting star counts rather than surface brightness profiles produces results that differ significantly due to the presence of mass segregation. The sample in each cluster is further broken down into various mass groups, each of which is fit independently, allowing us to determine how the concentration of each cluster varies with mass. The majority of the clusters in our sample show general agreement with the standard picture that more massive stars will be more centrally concentrated. We find that core radius vs. stellar mass can be fit with a two parameter power-law. The slope of this power-law is a value that describes the amount of mass segregation present in the cluster, and is measured independently of our distance from the cluster. This value correlates strongly with the core relaxation time and physical size of each cluster. Supplementary figures are also included showing the best fits and likelihood contours of fit parameters for all 54 clusters.

The cooling theory of neutron stars is corroborated by its comparison with observations of thermally emitting isolated neutron stars and accreting neutron stars in binary systems. An important ingredient for such an analysis is the age of the object, which, typically, is obtained from the spin-down history. This age is highly uncertain if the object's magnetic field varies appreciably over time. Other age estimators, such as supernova remnant ages and kinematic ages, only apply to few handful of neutron stars. We conduct a population synthesis study of the nearby isolated thermal emitters and obtain their ages statistically from the observed luminosity function of these objects. We argue that a more sensitive blind scan of the galactic disk with the upcoming space telescopes can help to constrain the ages to higher accuracy.

Based on the width of its main sequence, and an actual observed split when viewed through particular filters, it is widely accepted that 47 Tucanae contains multiple stellar populations. In this contribution, we divide the main-sequence of 47 Tuc into four color groups, which presumably represent stars of various chemical compositions. The kinematic properties of each of these groups is explored via proper-motions, and a strong signal emerges of differing proper-motion anisotropies with differing main-sequence color; the bluest main-sequence stars exhibit the largest proper-motion anisotropy which becomes undetectable for the reddest stars. In addition, the bluest stars are also the most centrally concentrated. A simiilar analysis for SMC stars, which are located in the background of 47 Tuc on our frames, yields none of the anisotropy exhibited by the 47 Tuc stars. We discuss implications of these results for possible formation scenarios of the various populations.

This paper presents new and efficient algorithms for matching stellar catalogues where the transformation between the coordinate systems of the two catalagoues is unknown and may include shearing. Finding a given object whether a star or asterism from the first catalogue in the second is logarithmic in time rather than polynomial, yielding a dramatic speed up relative to a naive implementation. Both acceleration of the matching algorithm and the ability to solve for arbitrary affine transformations not only will allow the registration of stellar catalogues and images that are now impossible to use but also will find applications in machine vision and other imaging applications.

We examine an external trigger mechanism that gives rise to the intense soft gamma-ray repeater (SGR) giant flares. Out of the three giant flares, two showcased the existence of a precursor, which we show to have had initiated the main flare. We develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. We show that the thinning time is on the order of the pre-flare quiescent time.

Diffraction is important when nearby substellar objects gravitationally lens distant stars. If the wavelength of the observation is comparable to the Schwarzschild radius of lensing object, diffraction leaves an observable imprint on the lensing signature. The SKA may have sufficient sensitivity to detect the typical sources, giant stars in the bulge. The diffractive signatures in a lensing event break the degeneracies between the mass of the lens, its distance and proper motion.

Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just beyond the half-light radius where standard models of globular cluster dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predictions for the stellar proper motions as a function of mass. The observed proper motions show no evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail with a straightforward model of the stellar distribution function. We do not find any evidence that the young white dwarfs have received a natal kick in contradiction with earlier results. Using the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the distance to NGC 6397 of $2.2^{+0.5}_{-0.7}$ kpc and a mass of the cluster of $1.1 \pm 0.1 \times 10^5 \mathrm{M}_\odot$ at the photometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory that will escape the cluster, yielding an estimate of the evaporation timescale, over which the number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, providing the first direct constraint on the relaxation time of the stars in a globular cluster of greater than or about 0.7 Gyr.

We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models are then used to determine the rate at which objects are leaving the main sequence, which must be the same as the rate at which objects are arriving on the white dwarf sequence in our field. The result is an empirically derived relation between temperature ($T_{eff}$) and time ($t$) on the white dwarf cooling sequence. Comparing this result to theoretical cooling models, we find general agreement with the expected slopes between 20,000K and 30,000K and between 6,000K and 20,000K, but the transition to the Mestel cooling rate of $T_{eff} \propto t^{-0.4}$ is found to occur at hotter temperatures, and more abruptly than is predicted by any of these models.

The mechanical properties of a neutron star crust, such as breaking strain and shear modulus, have implications for the detection of gravitational waves from a neutron star as well as bursts from Soft Gamma-ray Repeaters (SGRs). These properties are calculated here for three different crustal compositions for a non-accreting neutron star that results from three different cooling histories, as well as for a pure iron crust. A simple shear is simulated using molecular dynamics to the crustal compositions by deforming the simulation box. The breaking strain and shear modulus are found to be similar in the four cases, with a breaking strain of ~0.1 and a shear modulus of ~10^{30} dyne cm^{-2} at a density of \rho = 10^{14} g cm^{-3} for simulations with an initially perfect BCC lattice. With these crustal properties and the observed properties of {PSR J2124-3358} the predicted strain amplitude of gravitational waves for a maximally deformed crust is found to be greater than the observational upper limits from LIGO. This suggests that the neutron star crust in this case may not be maximally deformed or it may not have a perfect BCC lattice structure. The implications of the calculated crustal properties of bursts from SGRs are also explored. The mechanical properties found for a perfect BCC lattice structure find that crustal events alone can not be ruled out for triggering the energy in SGR bursts.

We investigate the possibility that a magnetic field may be present in the star $o-$Ceti (hereafter, Mira) and that the field plays a role in the star's mass loss. The model presented here is an application of an earlier derived theory that has been successfully employed for intermediate and high-mass evolved stars, and is now extended to the low-mass end. The modelling shows that it is possible to obtain a hybrid magnetohydrodynamic-dust-driven wind scenario for Mira, in which the role of a magnetic field in the equatorial plane of the star is dynamically important for producing a stellar wind. The wind velocity and the temperatures obtained from the model appear consistent with findings elsewhere.

We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, the proper motion dispersion in the core is found to be 179 +/- 17 microarcsec/yr. We find no evidence of anisotropy in the motions and no radial variation in the proper motions with respect to distance from the cluster center. We also set an upper limit on any central black hole to be ~150 Msun at 90% confidence level.

We present calculations for a magnetised hybrid wind model for Betelegeuse ($\alpha - $Orionis). The model is a direct application of our previously derived theory, combining a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star \citep[see][]{Thirumalai2010}. The resulting hybrid picture provides a mechanism for solving the problem of lifting stellar material up from the photosphere \citep[e.g.][] {Harper2009,Guandalini2006,Jura1984} and into the circumstellar envelope. It also predicts wind velocities in agreement with current estimates. Our approach reveals that magnetic fields in supergiant stars like Betelgeuse \citep[see][]{Auriere2010}, may play a vital role in determining the nature of the stellar outflow and consequently, opens a new avenue of investigation in the field of hybrid stellar winds.

We present calculations for a magnetised hybrid wind model for Asymptotic Giant Branch (AGB) stars. The model incorporates a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star. The resulting hybrid picture preserves traits of both types of winds. It is seen that this combination requires that the dust-parameter ($\Gamma_{d}$) be less than unity in order to achieve an outflow. The emergence of critical points in the wind changes the nature of the dust-driven outflow, simultaneously, the presence of a dust condensation radius changes the morphology of the magnetohydrodynamic (MHD) solutions for the wind. In this context, we additionally investigate the effect of having magnetic-cold spots on the equator of an AGB star and its implications for dust formation; which are seen to be consistent with previous findings.

The energy levels of hydrogen and helium atoms in strong magnetic fields are calculated in this study. The current work contains estimates of the binding energies of the first few low-lying states of these systems that are improvements upon previous estimates. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these one- and two-electron systems in a self-consistent manner. The method described herein is applicable to calculations of atomic structure in magnetic fields of arbitrary strength as it exploits the natural symmetries of the problem without assumptions of any basis functions for expressing the wave functions of the electrons or the commonly employed adiabatic approximation. The method is found to be readily extendable to systems with more than two electrons.

We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These images cover more than 60 square arcmins and extend over a radial range of 5-13.7 arcmin (6.5-17.9 pc) within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance modulus of (m - M)o=13.36+/-0.02+/-0.06 corresponding to a distance of 4.70+/-0.04+/-0.13 kpc, where the first error is random and the second is systematic. We also search the white dwarf photometry for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find no convincing cases within our sample.

The light pseudoscalar particle, dubbed the axion, borne out of the Peccei-Quinn solution to the strong CP problem in QCD remains elusive. One avenue of inferring its existence is through its coupling to electromagnetic radiation. So far, laboratory experiments have dedicated all efforts to detect the axion in the mass range $10^{-6} < m_a < 10^{-3}$ eV with a photon-axion coupling strength $g_{a\gamma\gamma} < 10^{-10}GeV^{-1}$, where the limits are derived from astrophysical considerations. In this study, we present a novel way of constraining $g_{a\gamma\gamma}$ by looking at the level of linear polarization in the radiation emerging from magnetic white dwarfs (mWDs). We find that photon-axion oscillations in WD magnetospheres can enhance the degree of linear polarization. Observing that most mWDs show only 5% linear polarization, we derive upper limits on $g_{a\gamma\gamma}$ for different axion masses.

We compute electromagnetic wave propagation through the magnetosphere of a magnetar. The magnetosphere is modeled as the QED vacuum and a cold, strongly magnetized plasma. The background field and electromagnetic waves are treated nonperturbatively and can be arbitrarily strong. This technique is particularly useful for examining non-linear effects in propagating waves. Waves travelling through such a medium typically form shocks; on the other hand we focus on the possible existence of waves that travel without evolving. Therefore, in order to examine the nonlinear effects, we make a travelling wave ansatz and numerically explore the resulting wave equations. We discover a class of solutions in a homogeneous plasma which are stabilized against forming shocks by exciting nonorthogonal components which exhibit strong nonlinear behaviour. These waves may be an important part of the energy transmission processes near pulsars and magnetars.

Microlensing is generally studied in the geometric optics limit. However, diffraction may be important when nearby substellar objects lens occult distant stars. In particular the effects of diffraction become more important as the wavelength of the observation increases. Typically if the wavelength of the observation is comparable to the Schwarzschild radius of lensing object, diffraction leaves an observable imprint on the lensing signature. The commissioning of the Square Kilometre Array (SKA) over the next decade begs the question of whether it will become possible to follow up lensing events with radio observations because the SKA may have sufficient sensitivity to detect the typical sources, giant stars in the bulge. The detection of diffractive lensing in a lensing event would place unique constraints on the mass of the lens and its distance. In particular it would distinguish rapidly moving stellar mass lenses (e.g. neutron stars) from slowly moving substellar objects such freely floating planets. An analysis of the sensitivity of the SKA along with new simple closed-form estimates of the expected signal applied to local exemplars for stellar radio emission reveals that this effect can nearly be detected with the SKA. If the radio emission from bulge giants is stronger than expected, the SKA could detect the diffractive microlensing signature from Earth-like interstellar planets in the solar neighborhood.

Gravitational lensing is generally treated in the geometric optics limit; however, when the wavelength of the radiation approaches or exceeds the Schwarzschild radius of the lens, diffraction becomes important. Although the magnification generated by diffractive gravitational lensing is well understood, the astrometric signatures of diffractive microlensing are first derived in this paper along with a simple closed-form bound for the astrometric shift. This simple bound yields the maximal shifts for substellar lenses in solar neighbourhood observed at 20~GHz, accessible to high sensitivity, high angular resolution radio telescopes such as the proposed Square Kilometre Array (SKA).

We examine two trigger mechanisms, one internal and the other external to the neutron star, that give rise to the intense soft gamma-ray repeater (SGR) giant flares. So far, three giant flares have been observed from the three out of the seven confirmed SGRs on March 5, 1979, August 27, 1998, and December 27, 2004. The last two events were found to be much more powerful than the first, and both showcased the existence of a precursor, that we show to have had initiated the main flare. In the internal mechanism, we propose that the strongly wound up poloidal magnetic field develops tangential discontinuities and dissipates its torsional energy in heating the crust. The timescale for the instability to develop coincides with the duration of the quiescent state that followed the precursor. Alternatively, we develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. Again, we show that the thinning time is on the order of the preflare quiescent time. This model naturally explains the origin of the observed nonthermal radiation during the flares, as well as the post flare radio afterglows.

We compare the expected abundance of cosmic-ray electrons and positrons from pulsars and magnetars. We assume that the distribution of infant pulsars and magnetars follows that of high-mass stars in the Milky Way and that the production rate of cosmic rays is proportional to the spin-down and magnetic-decay power of pulsars and magnetars, respectively. In combination with primary and secondary cosmic-ray leptons from other sources (especially supernova remnants), we find that both magnetars and pulsars can easily account for the observed cosmic-ray spectrum, in particular the dip seen by HESS at several TeV and the increase in positron fraction found by PAMELA.

We present a new pseudospectral algorithm for the calculation of the structure of atoms in strong magnetic fields. We have verified this technique for one, two and three-electron atoms in zero magnetic fields against laboratory results and find typically better than one-percent accuracy. We further verify this technique against the state-of-the-art calculations of hydrogen, helium and lithium in strong magnetic fields (up to about $2\times 10^{6}$ T) and find a similar level of agreement. The key enabling advantages of the algorithm are its simplicity (about 130 lines of commented code) and its speed (about $10^2-10^5$ times faster than finite-element methods to achieve similar accuracy).

We show that kicks generated by topological currents may be responsible for the large velocities seen in a number of pulsars. The majority of the kick builds up within the first second of the star's birth and generates a force about two orders of magnitude larger than a neutrino kick in the same temperature and magnetic field regime. Because of the nature of the topological currents the star's cooling is not affected until it reaches 10^9 K; thereafter the current replaces neutrino emission as the dominant cooling process. A requirement for the kick to occur is a suitably thin crust on the star; this leads us to speculate that pulsars with large kicks are quark stars and those with small kicks are neutron stars. If true this would be an elegant way to distinguish quark stars from neutron stars.

Microlensing and occultation are generally studied in the geometric optics limit. However, diffraction may be important when recently discovered Kuiper-Belt objects (KBOs) occult distant stars. In particular the effects of diffraction become more important as the wavelength of the observation and the distance to the KBO increase. For sufficiently distant and massive KBOs or Oort cloud objects not only is diffraction important but so is gravitational lensing. For an object similar to Eris but located in the Oort cloud, the signature of gravitational lensing would be detected easily during an occultation and would give constraints on the mass and radius of the object.

We review the polarization properties of X-ray emission from highly magnetized neutron stars, focusing on emission from the stellar surfaces. We discuss how x-ray polarization can be used to constrain neutron star magnetic field and emission geometry, and to probe strong-field quantum electrodynamics and possibly constrain the properties of axions.

We have investigated the crustal properties of neutron stars without fallback accretion. We have calculated the chemical evolution of the neutron star crust in three different cases (a modified Urca process without the thermal influence of a crust, a thick crust, and a direct Urca process with a thin crust) in order to determine the detailed composition of the envelope and atmosphere as the nuclear reactions freeze out. Using a nuclear reaction network up to technetium, we calculate the distribution of nuclei at various depths of the neutron star. The nuclear reactions quench when the cooling timescale is shorter than the inverse of the reaction rate. Trace light elements among the calculated isotopes may have enough time to float to the surface before the layer crystallizes and form the atmosphere or envelope of the neutron star. The composition of the neutron-star envelope determines the total photon flux from the surface, and the composition of the atmosphere determines the emergent spectrum. Our calculations using each of the three cooling models indicate that without accretion of fallback the neutron star atmospheres are dependent on the assumed cooling process of the neutron star. Each of the cooling methods have different elements composing the atmosphere: for the modified Urca process the atmosphere is $^{28}$Si, the thick crust has an atmosphere of $^{50}$Cr, and the thin crust has an atmosphere of $^{40}$Ca. In all three cases the atmospheres are composed of elements which are lighter then iron.

A fully relativistic treatment of Bernstein waves in an electron-positron pair plasma has remained too formidable a task owing to the very complex nature of the problem. In this article, we perform contour integration of the dielectric response function and numerically compute the dispersion curves for a uniform, magnetized, relativistic electron-positron pair plasma. The behavior of the dispersion solution for several cases with different plasma temperatures is highlighted. In particular, we find two wave modes that exist only for large wavelengths and frequencies similar to the cyclotron frequency in a moderately relativistic pair plasma. The results presented here have important implications for the study of those objects where a hot magnetized electron-positron plasma plays a fundamental role in generating the observed radiation.

Observations of white dwarfs in the globular clusters NGC 6397 and Omega Centauri indicate that these stars may get a velocity kick during their time as giants. If the mass loss while on the asymptotic giant branch is slightly asymmetric, the resulting white dwarf could be born with such a velocity kick. These energetic white dwarfs will impart their excess energy on other stars as they travel through the cluster. A Monte-Carlo simulation of the white-dwarfs kicks combined with estimate of the phase-space diffusion of the white dwarfs reveals that as the white dwarfs equilibrate, they lose most of their energy in the central region of the cluster. They could possibly augment the effect of binaries, delaying core collapse or increasing the size of globular cluster cores.

Recent observations of white dwarfs in globular clusters indicate that these stars may get a velocity kick during their time as giants. This velocity kick could originate naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. The kicks may be large enough to dramatically change the radial distribution of young white dwarfs, giving them larger energies than other stars in the cluster. As these energetic white dwarfs travel through the cluster they can impart their excess energy on the other stars in the cluster. This new heat source for globular clusters is expected to be largest during the clusters' youth.

We are investigating mass fractions on the crust of a neutron star which would remain after one year of cooling. We use cooling curves corresponding with various densities, or depths, of the neutron star just after its formation. We assume the modified Urca process dominates the energy budget of the outer layers of the star in order to calculate the temperature of the neutron star as a function of time. Using a nuclear reaction network up to technetium, we calculate how the distribution of nuclei quenches at various depths of the neutron star crust. The initial results indicate that $^{28}$Si is the lightest isotope to be optically thick on the surface after one year of cooling.

If the winds of an asymptotic-giant-branch stars are sufficiently strong are slightly asymmetric, they can alter the star's trajectory through a globular cluster; therefore, if these winds are asymmetric, one would expect young white dwarfs to be less radially concentrated than either their progenitors or older white dwarfs in globular clusters. This latter effect has recently been observed. Additionally the young white dwarfs should have larger typical velocities than their progenitors. After phase mixing this latter effect is vastly diminished relative to the changes in the spatial distribution of young white dwarfs with kicks, so it is more difficult to detect than the change in the spatial distribution. The most powerful kinematic signature is the change in the eccentricity of the orbits that is revealed through the distribution of the position angles of proper motion.

Recent observations of white dwarfs in globular clusters indicate that these stars may get a velocity kick during their time as giants. This velocity kick could originate naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. If white dwarfs get a kick comparable to the orbital velocity of the binary, the initial Runge-Lenz vector (eccentricity vector) of the orbit is damped to be replaced by a component pointing toward the cross product of the initial angular momentum and the force. The final eccentricity may be of order unity and if the kick is sufficiently large, the system may be disrupted. These results may have important ramifications for the evolution of binary stars and planetary systems

The wind of an asymptotic-giant-branch stars is sufficiently strong that if it is slightly asymmetric, it can propel the star outside of the open cluster of its birth or significantly alter its trajectory through a globular cluster; therefore, if these stellar winds are asymmetric, one would expect a deficit of white dwarfs of all ages in open clusters and for young white dwarfs to be less radially concentrated than either their progenitors or older white dwarfs in globular clusters. This latter effect has recently been observed. Hence, detailed studies of the radial distribution of young white dwarfs in globular clusters could provide a unique probe of mass loss on the asymptotic giant branch and during the formation of planetary nebulae both as a function of metallicity and a limited range of stellar mass.

Magnetars, neutron stars with ultra strong magnetic fields ($B\sim 10^{14} - 10^{15}$G), manifest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray pulsars. This study estimates the birthrate of magnetars to be $\sim$ 0.22 per century with a galactic population comprising of $\sim$17 objects. A population synthesis was carried out based on the five anomalous X-ray pulsars detected in the ROSAT All-Sky Survey by comparing their number to that of massive OB stars in a well defined volume. Additionally, the group of seven X-ray dim isolated neutron stars detected in the same survey were found to have a birthrate of $\sim$ 2 per century with a galactic population of $\sim$ 22,000 objects.

We propose that the presence of additional planets in extrasolar planetary systems can be detected by long-term transit timing studies. If a transiting planet is on an eccentric orbit then the presence of another planet causes a secular advance of the transiting planet's pericenter over and above the effect of general relativity. Although this secular effect is impractical to detect over a small number of orbits, it causes long-term differences in when future transits occur, much like the long-term decay observed in pulsars. Measuring this transit-timing delay would thus allow the detection of either one or more additional planets in the system or the first measurements of non-zero oblateness ($J_2$) of the central stars.

Recent observations and theoretical investigations of neutron stars indicate that their atmospheres consist not of hydrogen or iron but possibly other elements such as helium. We calculate the ionization and dissociation equilibrium of helium in the conditions found in the atmospheres of magnetized neutron stars. For the first time this investigation includes the internal degrees of freedom of the helium molecule. We found that at the temperatures and densities of neutron star atmospheres the rotovibrational excitations of helium molecules are populated. Including these excitations increases the expected abundance of molecules by up to two orders of magnitude relative to calculations that ignore the internal states of the molecule; therefore, if the atmospheres of neutron stars indeed consist of helium, helium molecules and possibly polymers will make the bulk of the atmosphere and leave signatures on the observed spectra from neutron stars. We applied our calculation to nearby radio-quiet neutron stars with B_{dipole}~10^{13}-10^{14} G. If helium comprises their atmospheres, our study indicates that isolated neutron stars with T_{BB}~10^6 K such as RXJ0720.4-3125 and RXJ1605.3+3249 will have He^+ ions predominantly, while isolated neutron stars with lower temperature (T_{BB}~5x10^5 K) such as RXJ1856.5-3754 and RXJ0420.0-5022 will have some fraction of helium molecules. If helium molecules are abundant, their spectroscopic signatures may be detected in the optical, UV and X-ray band. (Abridged)

Recent experiments suggest that polarized photons may couple significantly to pseudoscalar particles such as axions. We study the possible observational signatures of axion-photon coupling for radiation from magnetic stars, with particular focus on neutron stars. We present general methods for calculating the axion-photon conversion probability during propagation through a varying magnetized vacuum as well as across an inhomogeneous atmosphere. Partial axion-photon conversion may take place in the vacuum region outside the neutron star. Strong axion-photon mixing occurs due to a resonance in the atmosphere, and depending on the axion coupling strength and other parameters, significant axion-photon conversion can take place at the resonance. Such conversions may produce observable effects on the radiation spectra and polarization signals from the star. We also apply our results to axion-photon propagation in the Sun and in magnetic white dwarfs. We find that there is no appreciable conversion of solar axions to photons during the propagation.

The modes of vibration of hanging and partially supported strings provide useful analogies to scalar fields travelling through spacetimes that admit conformally flat spatial sections. This wide class of spacetimes includes static, spherically symmetric spacetimes. The modes of a spacetime where the scale factor depends as a power-law on one of the coordinates provide a useful starting point and yield a new classification of these spacetimes on the basis of the shape of the string analogue. The family of corresponding strings follow a family of curves related to the cycloid, denoted here as hypercycloids (for reasons that will become apparent). Like the spacetimes that they emulate these strings exhibit horizons, typically at their bottommost points where the string tension vanishes; therefore, hanging strings may provide a new avenue for the exploration of the quantum mechanics of horizons.

Owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamical fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially observable, fraction of the energy ($\sim 10^{33}$ ergs) can generate a non-thermal distribution of electrons and positrons far from the star. This paper examines the cooling and radiative output of these particles. Small-scale waves may produce only the non-thermal emission. The properties of this non-thermal emission in the absence of a fireball match those of the quiescent, non-thermal radiation recently observed non-thermal emission from several anomalous X-ray pulsars and soft-gamma repeaters. Initial estimates of the emission as a function of angle indicate that the non-thermal emission should be beamed and therefore one would expect this emission to be pulsed as well. According to this model the pulsation of the non-thermal emission should be between 90 and 180 degrees out of phase from the thermal emission from the stellar surface.

We study the propagation of polarization light through the magnetosphere of neutron stars. At intermediate frequencies (the optical through the infrared), both the birefringence induced by the plasma and by quantumelectrodynamics influence the observation polarization of radiation from the surface of the neutron star. Because these two processes compete in this regime, we find that polarization observations can constrain the properties of the neutron-star magnetosphere, specifically the total charge density. We calculate both the phase-resolved and the phase-averaged polarization signatures induced by magnetospheric birefringence.

The fact that we apparently live in an accelerating universe places limitations on where humans might visit. If the current energy density of the universe is dominated by a cosmological constant, a rocket could reach a galaxy observed today at a redshift of 1.7 on a one-way journey or merely 0.65 on a round trip. Unfortunately these maximal trips are impractical as they require an infinite proper time to traverse. However, calculating the rocket trajectory in detail shows that a rocketeer could nearly reach such galaxies within a lifetime (a long lifetime admittedly -- about 100 years). For less negative values of $w$ the maximal redshift increases becoming infinite for $w\geq -1/3$.

Ultramagnetized neutron stars or magnetars are magnetically powered neutron stars. Their strong magnetic fields dominate the physical processes in their crusts and their surroundings. The past few years have seen several advances in our theoretical and observational understanding of these objects. In spite of a surfeit of observations, their spectra are still poorly understood. I will discuss the emission from strongly magnetized condensed matter surfaces of neutron stars, recent advances in our expectations of the surface composition of magnetars and a model for the non-thermal emission from these objects.

Rossby waves (or r-modes) on the surface of a neutron star have become a leading model for the oscillations observed during the tail of Type-I X-ray bursts. Their frequency evolution matches well with the observed frequency drifts of the oscillations in the bursts, and the burning appears to excite these waves quite naturally. This paper addresses the detailed shape of the expected flux profiles from r-modes on neutron stars as a function of energy. R-modes naturally account for both the small amplitude of the observed oscillations and their lack of harmonic content. However, the model predicts that the oscillation at higher energies leads the lower energy variation. The observed oscillations have the opposite trend which possibly indicates that the higher energy photons are upscattered in the plasma surrounding the neutron star and therefore delayed.

In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect couples the direction of the polarization of photons leaving the NS surface, to the direction of the magnetic field along the ray's path. We analyze the consequences that this effect has on aligning the polarization vectors to generate large net polarizations, while considering thermal radiation originating from a thermal hydrogen atmosphere. Counter to previous predictions, we show that the thermal radiation should be highly polarized even in the optical. When detected, this polarization will be the first demonstration of vacuum birefringence. It could be used as a tool to prove the high magnetic field nature of AXPs and it could also be used to constrain physical NS parameters, such as $R/M$, to which the net polarization is sensitive.

Previously, we showed that, owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamic fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially observable, fraction of the energy ($\sim 10^{33}$ ergs) can generate a non-thermal distribution of electrons and positrons far from the star. In this paper, we examine the cooling and radiative output of these particles. We also investigate the properties of non-thermal emission in the absence of a fireball to understand the breakdown of fast modes that do not yield an optically thick pair plasma. This quiescent, non-thermal radiation associated with fast mode breakdown may account for the recently observed non-thermal emission from several anomalous X-ray pulsars and soft-gamma repeaters.

We propose a model to account for the bursts from soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) in which quantum electrodynamics (QED) plays a vital role. In our theory, that we term ``fast-mode breakdown,'' magnetohydrodynamic (MHD) waves that are generated near the surface of a neutron star and propagate outward through the magnetosphere will be modified by the polarization of the vacuum. For neutron star magnetic fields $B_\rmscr{NS} \gtrsim B_\rmscr{QED} \approx 4.4 \times 10^{13}$ G, the interaction of the wave fields with the vacuum produces non-linearities in fast MHD waves that can steepen in a manner akin to the growth of hydrodynamic shocks. Under certain conditions, fast modes can develop field discontinuities on scales comparable to an electron Compton wavelength, at which point the wave energy will be dissipated through electron-positron pair production. We show that this process operates if the magnetic field of the neutron star is sufficiently strong and the ratio of the wavelength of the fast mode to its amplitude is sufficiently small, in which case the wave energy will be efficiently converted into an extended pair-plasma fireball. The radiative output from this fireball will consist of hard X-rays and soft $\gamma$-rays, with a spectrum similar to those seen in bursts from SGRs and AXPs (abridged).

(abridged) An outburst of more than 80 individual bursts, similar to those seen from Soft Gamma Repeaters (SGRs), was detected from the Anomalous X-ray Pulsar (AXP) 1E 2259+586 in 2002 June. Coincident with this burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse profile and spin down of the underlying X-ray source. We present RXTE and XMM-Newton observations of 1E 2259+586 that show the evolution of the aforementioned source parameters during and following this episode. Specifically, we observe an X-ray flux increase by more than an order of magnitude having two distinct components. The first component is linked to the burst activity and decays within ~2 days during which the energy spectrum is considerably harder than during the quiescent state of the source. The second component decays over the year following the glitch according to a power law in time with an exponent -0.22 +/- 0.01. The pulsed fraction decreased initially to ~15% RMS, and the pulse profile changed significantly during the outburst. A glitch was observed in 1E 2259+586 that preceded the observed burst activity. A fraction of the glitch (~19%) recovered, although the recovery was not purely exponential. An exponential rise of ~20% of the frequency jump with a time scale of ~14 days results in a significantly better fit to the data, however, contamination from a systematic drift in the phase of the pulse profile cannot be excluded. The long-term post-glitch spin-down rate decreased in magnitude relative to the pre-glitch value. A comparison with SGR outburst properties, a physical interpretation of the results, and implications on the number of magnetar candidates in our Galaxy are discussed.

During a Type-I burst, the turbulent deflagation front may excite waves in the neutron star ocean and upper atmosphere with frequencies, $\omega \sim 1$ Hz. These waves may be observed as highly coherent flux oscillations during the burst. The frequencies of these waves changes as the upper layers of the neutron star cool which accounts for the small variation in the observed QPO frequencies. In principle several modes could be excited but the fundamental buoyant $r-$mode exhibits significantly larger variability for a given excitation than all of the other modes. An analysis of modes in the burning layers themselves and the underlying ocean shows that it is unlikely these modes can account for the observed burst oscillations. On the other hand, photospheric modes which reside in a cooler portion of the neutron star atmosphere may provide an excellent explanation for the observed oscillations.

We present a global linear stability analysis of nuclear fuel accumulating on the surface of an accreting neutron star and we identify the conditions under which thermonuclear bursts are triggered. The analysis reproduces all the recognized regimes of hydrogen and helium bursts, and in addition shows that at high accretion rates, near the limit of stable burning, there is a regime of ``delayed mixed bursts'' which is distinct from the more usual ``prompt mixed bursts.'' In delayed mixed bursts, a large fraction of the fuel is burned stably before the burst is triggered. Bursts thus have longer recurrence times, but at the same time have somewhat smaller fluences. Therefore, the parameter alpha, which measures the ratio of the energy released via accretion to that generated through nuclear reactions in the burst, is up to an order of magnitude larger than for prompt bursts. This increase in alpha near the threshold of stable burning has been seen in observations. We explore a wide range of mass accretion rates, neutron star radii and core temperatures, and calculate a variety of burst properties. From a preliminary comparison with data, we suggest that bursting neutron stars may have hot cores, with T_{core} >~ 10^{7.5} K, consistent with interior cooling via the modified URCA or similar low-efficiency process, rather than T_{core} ~ 10^7 K, as expected for the direct URCA process. There is also an indication that neutron star radii are somewhat small <~ 10 km. Both of these conclusions need to be confirmed by comparing more careful calculations with better data.

We study the rotational evolution of a protoneutron star with hyperons and nucleons or solely nucleons in its core due to the escape of the trapped neutrinos. It is found that at the early stage of its evolution, the stellar crust contracts significantly, consequently the star spins up. At the late stage, a the protoneutron star with hyperons, it keeps shrinking and spinning up till all the trapped neutrinos escape. Consequently, the distribution of the stellar initial spin periods is skewed toward shorter periods. For a protoneutron star with only nucleons, the expansion of its core dominates, and the stellar rotation slows down. After the neutrinos escape, the range of the spin periods is narrower than the initial one, but the distribution is still nearly uniform. If the hyperonic star is metastable, its rotational frequency accelerates distinguishedly before it collapses to a black hole.

The advent of large panoramic photometric surveys of the sky offers the possibly of exploring the association of gamma-ray bursts (GRBs) with supernovae. To date, a few gamma-ray bursts have been connected possibly with supernovae: GRB 980425 -- SN 1998bw, GRB 011121 -- SN 2001ke, GRB 970228 and GRB 980326. A combination of a large detection rate of GRBs and rapid coverage of a large portion of the sky to faint magnitude limits offers the possibility of detecting a supernova preceding an associated GRB or at least placing limits on the rate of association between these two phenomena and the time delay between them. This would provide important constraints on theoretical models for gamma-ray bursts.

General relativity limits the compactness of static stars. If the pressure of the fluid is positive and the density decreases with distance from the center, the value of the circumferential radius of the star must be greater than (9/4) G M/c^2, or equivalently the redshift of the surface must be less than two. If constraints on the equation of state of the material are relaxed, general relativity alone does not restrict the redshift of a static stellar surface. However, because black hole candidates in the universe generally accrete material from their environs, the process of accretion provides upper limits on the redshift of a astrophysical black-hole candidates.

In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect decouples the polarization modes of photons leaving the NS surface. Both the total intensity and the intensity in each of the two modes is preserved along a ray's path through the neutron-star magnetosphere. We analyze the consequences that this effect has on aligning the observed polarization vectors across the image of the stellar surface to generate large net polarizations. Counter to previous predictions, we show that the thermal radiation of NSs should be highly polarized even in the optical. When detected, this polarization will be the first demonstration of vacuum birefringence. It could be used as a tool to prove the high magnetic field nature of AXPs and it could also be used to constrain physical NS parameters, such as $R/M$, to which the net polarization is sensitive.

We present a broadband modeling of the afterglow of GRB 021004. The optical transient of this burst has been detected very early and densely sampled in several bands. Its light curve shows significant deviations from a simple power law. We use the data from the X-ray to the $I-$band gathered in the first month of observations, and examine three models. Two involve variations in the energy of the shock. The first (energy injection) allows only increases to the shock energy, while the second (patchy shell) allows the energy to increase or decrease. In the final model (clumpy medium) the energy of the shock is constant while the density varies. While all three models reproduce well the optical bands, the variable density model can best account for the X-ray data, and the energy-injection model has the poorest fit. None of the models can account for the modest color variations observed during the first few days of the burst.

Type I X-ray bursts are very common in neutron star X-ray binaries, but no Type I burst has been seen in the dozen or so binaries in which the accreting compact star is too massive to be a neutron star and therefore is identified as a black hole candidate. We have carried out a global linear stability analysis of the accumulating fuel on the surface of a compact star to identify the conditions under which thermonuclear bursts are triggered. Our analysis, which improves on previous calculations, reproduces the gross observational trends of bursts in neutron star systems. It further shows that, if black hole candidates have surfaces, they would very likely exhibit instabilities similar to those that lead to Type I bursts on neutron stars. The lack of bursts in black hole candidates is thus significant, and indicates that these objects have event horizons. We discuss possible caveats to this conclusion.

We report polarization measurements and photometry for the optical afterglow of the gamma-ray burst GRB 020405. We measured a highly significant 9.9% polarization (in V band) 1.3 days after the burst and argue that it is intrinsic to the GRB. The light curve decay is well fitted by a $t^{-1.72}$ power-law; we do not see any evidence for a break between 1.24 and 4.3 days after the burst. We discuss these measurements in the light of several models of GRB afterglows.

Andersson et al. and Bildsten proposed that the spin of accreting neutron stars is limited by removal of angular momentum by gravitational radiation which increases dramatically with the spin frequency of the star. Both Bildsten and Andersson et al. argued that the r-modes of the neutron star for sufficiently quickly rotating and hot neutron stars will grow due to the emission of gravitational radiation, thereby accounting for a time varying quadrupole component to the neutron star's mass distribution. However, Levin later argued that the equilibrium between spin-up due to accretion and spin-down due to gravitational radiation is unstable, because the growth rate of the r-modes and consequently the rate of gravitational wave emission is an increasing function of the core temperature of the star. The system executes a limit cycle, spinning up for several million years and spinning down in less than a year. However, the duration of the spin-down portion of the limit cycle depends sensitively on the amplitude at which the nonlinear coupling between different r-modes becomes important. As the duration of the spin-down portion increases the fraction of accreting neutron stars which may be emitting gravitational radiation increases while the peak flux in gravitational radiation decreases. Depending on the distribution of quickly rotating neutron stars in the Galaxy and beyond, the number of gravitational emitters detectable with LIGO may be large.

The thermal emission of strongly magnetized neutron-star atmospheres is thought to be highly polarized. However, because of the different orientations of the magnetic field over the surface of the neutron star (NS), it is commonly assumed that the net observed polarization will be significantly reduced as the polarization from different regions will cancel each other. We show that the birefringence of the magnetized QED vacuum decouples the polarization modes in the magnetosphere; therefore, the direction of the polarization follows the direction of the magnetic field up to a large distance from the stellar surface. At this distance, the rays that leave the surface and are destined for our detectors pass through only a small solid angle; consequently, the polarization direction of the emission originating in different regions will tend to align together. The net observed polarization of the thermal radiation of NSs should therefore be very large. Measurement of this polarization will be the first direct evidence of the birefringence of the magnetized vacuum due to QED and a direct probe of behavior of the vacuum at magnetic fields of order of and above the critical QED field of 4.4 x 10 13 G. The large observable polarization will also help us learn more about the atmospheric properties of NSs.

Roughly a dozen X-ray binaries are presently known in which the compact accreting primary stars are too massive to be neutron stars. These primaries are identified as black holes, though there is as yet no definite proof that any of the candidate black holes actually possesses an event horizon. We discuss how Type I X-ray bursts may be used to verify the presence of the event horizon in these objects. Type I bursts are caused by thermonuclear explosions when gas accretes onto a compact star. The bursts are commonly seen in many neutron star X-ray binaries, but they have never been seen in any black hole X-ray binary. Our model calculations indicate that black hole candidates ought to burst frequently if they have surfaces. Based on this, we argue that the lack of bursts constitutes strong evidence for the presence of event horizons in these objects.

If the potential of a scalar field, Phi, which currently provides the ``dark energy'' of the universe, has a negative minimum -M_0^4, then quantum-mechanical fluctuations could nucleate a bubble of Phi at a negative value of the potential. The bubble would expand at the speed of light. Given that no such bubble enveloped us in the past, we find that any minimum in V(Phi) must be separated from the current Phi value by more than min{1.5M_0, 0.21M_Pl}, where M_Pl is the Planck mass. We also show that vacuum decay renders a cyclic or ekpyrotic universe with M_0^4 > 10^{-10}M_Pl^4, untenable.

Recent developments in the spectropolarimetric study of compact objects, specifically black holes (stellar and massive) and neutron stars are reviewed. The lectures are organized around five topics: disks, jets, outflows, neutron stars and black holes. They emphasize physical mechanisms and are intended to bridge the gap between the fundamentals of polarimetry and the phenomenology of observed cosmic sources of polarized radiation, as covered by the other lecturers. There has been considerable recent progress in spectropolarimetry from radio through optical frequencies and this is producing some unique diagnostics of the physical conditions around compact objects. It is argued that there is a great need to develop a correspondingly sensitive polarimetric capability at ultraviolet through gamma-ray energies.

We present a simple, unified model which accounts for properties of the X-ray emission from the three middle-aged pulsars PSR 1055-52, PSR 0656+14 and PSR 0630+18 (Geminga). The X-ray radiation from these objects is pulsed more strongly at energies above a transition point around 0.5 keV. In addition, the phase of the pulses shifts by about 80-100 degrees around the same point. Geminga also has the peculiarity that its pulsed fraction decreases in the 0.3-0.5 keV energy range, attaining a minimum near 0.5 keV. We show that a two-component hydrogen atmosphere is able to account for these disparate features. In our model, the hotter component is powered by particle bombardment and is restricted to the polar regions, while the cooler one covers the entire stellar surface. The two components also differ in their emission patterns, with the hard and soft contributions coming from areas radiating into fan and pencil beams, respectively.

The energy source powering the X-ray emission from anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) is still uncertain. In one scenario, the presence of an ultramagnetized neutron star, or ``magnetar'', with B on the order of 10^{14} - 10^{15} G is invoked. To investigate this hypothesis, we have analyzed archival ASCA data for several known AXPs and SGRs, and fitted them with a model in which all or part of the X-ray flux originates as thermal emission from a magnetar. Our magnetar spectral model includes the effects of the anisotropy of the heat flow through an ultramagnetized neutron star envelope, reprocessing by a light element atmosphere, and general relativistic corrections to the observed spectrum. We obtain good fits to the data with radii for the emitting areas which are generally consistent with those expected for neutron stars, in contrast to blackbody (BB) fits, which imply much smaller radii. Furthermore, the inclusion of atmospheric effects results in inferred temperatures which are lower than those implied by BB fits, but however still too high to be accounted by thermal cooling alone. An extra source of heating (possibly due to magnetic field decay) is needed. Despite the harder tail in the spectrum produced by reprocessing of the outgoing flux through the atmosphere, spectral fits still require a considerable fraction of the flux to be in a power-law component.

The curved spacetime surrounding a rotating black hole dramatically alters the structure of nearby electromagnetic fields. The Wald field which is an asymptotically uniform magnetic field aligned with the angular momentum of the hole provides a convenient starting point to analyze the effects of radiative corrections on electrodynamics in curved spacetime. Since the curvature of the spacetime is small on the scale of the electron's Compton wavelength, the tools of quantum field theory in flat spacetime are reliable and show that a rotating black hole immersed in a magnetic field approaching the quantum critical value of $B_k=m^2 c^3/(e\hbar) \approx 4.4 \times 10^{13}$~G $\approx 1.3\times10^{-11}$ cm$^{-1}$ is unstable. Specifically, a maximally rotating three-solar-mass black hole immersed in a magnetic field of $2.3 \times 10^{12}$~G would be a copious producer of electron-positron pairs with a luminosity of $3 \times 10^{52}$ erg s$^{-1}$.

The spectacular first light observation by the Chandra X-Ray Observatory revealed an X-ray point source near the center of the 300 yr old Cas A supernova remnant. We present an analysis of the public X-ray spectral and timing data. No coherent pulsations were detected in the Chandra/HRC data. The 3-sigma upper limit on the pulsed fraction is <35% for P>20 ms. The Chandra/ACIS spectrum of the point source may be fit with an ideal blackbody (kT=0.5 keV), or with BB models modified by the presence of a NS atmosphere (kT=0.25-0.35 keV), but the temperature is higher and the inferred emitting area lower than expected for a 300 yr old NS according to standard cooling models. The spectrum may also be fit with a power law model (photon index 2.8-3.6). Both the spectral properties and the timing limits of the point source are inconsistent with a young Crab-like pulsar, but are quite similar to the properties of the anomalous X-ray pulsars. The spectral parameters are also very similar to those of the other radio-quiet X-ray point sources in the supernova remnants Pup A, RCW 103, and PKS 1209-52. Current limits on an optical counterpart for the Cas A point source rule out models that invoke fallback accretion onto a compact object if fallback disk properties are similar to those in quiescent low-mass X-ray binaries. However, the optical limits are marginally consistent with plausible alternative assumptions for a fallback disk. In this case, accreting NS models can explain the X-ray data, but an accreting BH model is not promising.

We examine properties of thermal radiation emitted by strongly magnetized neutron stars (NSs). In particular, we show that the pulsation amplitudes of the energy-integrated flux are an increasing function of the absorption column density to the source. This is especially important for the interpretation of the Anomalous $X$-ray Pulsars (AXPs) as cooling neutron stars with high magnetic fields. We show that the high-pulsation amplitudes observed in these objects are consistent with cooling models, if the large amount of absorption to these sources is taken into account. We also show that cooling models imply inferred radii of the emitting regions on the order of $\sim 5-6$ times smaller than the actual NS radii, again in agreement with observations.

The increase in spin frequency as the burning atmospheres of TypeI X-ray bursts cool provides a strong constraint on the radius of the underlying neutron star. If the change in spin frequency is due to a change in the thickness of the atmosphere, the radius of the star must exceed $3 G M/c^2$ for any equation of state and approximately $3.5 G M/c^2$ for most physically reasonable equations of state. This constraint arises because the direction of the Coriolis force for radial motion reverses for $R < G / c^2 (M + I / R^2)$. Furthermore, the marked change in the magnitude of the Coriolis force near compact stars provides a straightforward explanation for why the frequency of the quickly rotating bursters shifts by the same amount as the slow rotators; they are slightly more massive, 1.6M$_\odot$ versus 1.4M$_\odot$. \

The magnetic fields of white dwarfs distort their shape generating an anisotropic moment of inertia. A magnetized white dwarf which rotates obliquely relative to the symmetry axis has a mass quadrupole moment which varies in time, so it will emit gravitational radiation. LISA may be able to detect the gravitational waves from two nearby, quickly rotating white dwarfs.

Recent observations of 1E~161348-5055, the neutron-star candidate at the center of the supernova remnant RCW 103, show that a component of its emission varies sinusoidally with a period of approximately six hours. We argue that this period is what one would expect for a freely precessing neutron star with a spin period of about one second. We produce light curves for a freely precessing neutron star with a hotspot. By a suitable choice of parameters, we obtain light curves which are constant with rotational phase when the flux from the star reaches a maximum. At other phases of the precession, the flux varies as the star rotates but the total flux decreases by a factor of several. These models can explain the behavior observed from 1E~161348-5055 and predict that the spin period should be detectable at minimum flux from sufficiently sensitive measurements.

Extremely strong magnetic fields change the vacuum index of refraction. Although this polarization dependent effect is small for typical neutron stars, it is large enough to decouple the polarization states of photons traveling within the field. The photon states evolve adiabatically and follow the changing magnetic field direction. The combination of a rotating magnetosphere and a frequency dependent state decoupling predicts polarization phase lags between different wave bands, if the emission process takes place well within the light cylinder. This QED effect may allow observations to distinguish between different pulsar emission mechanisms and to reconstruct the structure of the magnetosphere.

We examine the timing history of several anomalous X-ray pulsars (AXPs), and find that they exhibit spin down irregularities statistically similar to those of radio pulsars. We propose that these irregularities are simply glitches like those of radio pulsars and fit glitching spin down solutions to the data available for 1E~2259+586, 1E~1048.1-5937 and 4U~0142+61. The inferred magnitude of the glitches ($\Delta \Omega$) is comparable to that exhibited in radio pulsar glitching. With these results we argue that the three AXPs that we studied are isolated neutron stars spinning down by magnetic dipole radiation and powered by a combination of cooling, magnetic field decay and magnetospheric particle bombardment.

Extremely strong magnetic fields change the vacuum index of refraction. This induces a lensing effect that is not unlike the lensing phenomenon in strong gravitational fields. The main difference between the two is the polarization dependency of the magnetic lensing, a behaviour that induces a handful of interesting effects. The main prediction is that the thermal emission of neutron stars with extremely strong magnetic fields is polarized - up to a few percent for the largest fields known. This potentially allows a direct method for measuring their magnetic fields.

We examine wave propagation and the formation of shocks in strongly magnetized plasmas by applying a variational technique and the method of characteristics to the coupled magnetohydrodynamic (MHD) and quantum-electrodynamic (QED) equations of motion. In sufficiently strong magnetic fields such as those found near neutron stars, not only is the plasma extremely relativistic but the effects of QED must be included to understand processes in the magnetosphere. As Thompson & Blaes [1] find, the fundamental modes in the extreme relativistic limit of MHD coupled with QED are two oppositely directed Alfv\'{e}n modes and the fast mode. QED introduces nonlinear couplings which affect the propagation of the fast mode such that waves traveling in the fast mode evolve as vacuum electromagnetic ones do in the presence of an external magnetic field [2] (Heyl & Hernquist 1998). The propagation of a single Alfv\'{e}n mode is unaffected but QED does alter the coupling between the Alfv\'{e}n modes.

Recently launched x-ray telescopes have discovered several candidate isolated neutron stars. The thermal radiation from these objects may potentially constrain our understanding of nuclear physics in a realm inaccessible to terrestrial experiments. To translate the observed fluxes from neutron stars into constraints, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We describe models of the thermal structure of the envelopes of neutron stars with magnetic fields up to 10^{14} G. Unlike earlier work, we infer the properties of envelope models in two dimensions and precisely account for the quantization of the electron phase space. Both dipole and uniformly magnetized envelopes are considered.

Ultramagnetized neutron stars or magnetars have been invoked to explain several astrophysical phenomena. We examine how the magnetic field of a magnetar will decay over time and how this decay affects the cooling of the object. We find that for sufficiently strong nascent fields, field decay alters the cooling evolution significantly relative to similarly magnetized neutron stars with constant fields. As a result, old magnetars can be expected to be bright in the soft X-ray band. The soft X-ray source RXJ~0720.4$-$3125 may well be the nearest such old magnetar.

We calculate the atomic structure of hydrogen and helium, atoms and molecules in an intense magnetic field, analytically and numerically with a judiciously chosen basis.

We examine the propagation of electromagnetic radiation through a strong magnetic field using the method of characteristics. Owing to nonlinear effects associated with vacuum polarization, such waves can develop discontinuities analogous to hydrodynamical shocks. We derive shock jump conditions and discuss the physical nature of these non-linear waves

Recent ROSAT measurements show that the x-ray emission from isolated neutron stars is modulated at the stellar rotation period. To interpret these measurements, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We present nearly analytic models of the thermal structure of the envelopes of ultramagnetized neutron stars. Specifically, we examine the limit in which only the ground Landau level is filled. We use the models to estimate the amplitude of modulation expected from non-uniformities in the surface temperatures of strongly magnetized neutron stars. In addition, we estimate cooling rates for stars with fields $B \sim 10^{15}-10^{16}$ G which are relevant to models that invoke ``magnetars'' to account for soft $\gamma$-ray emission from some repeating sources.

Recent observations of the compact source embedded within the supernova remnant RCW 103 rekindle interest in the origin of this object's emission. We contrast several models in which neutron-star cooling powers RCW 103. Specifically, either the presence of an accreted envelope or a sufficiently intense magnetic field can account for the X-ray emission from this object.

RX J0720.4-3125 has recently been identified as a pulsating soft X-ray source in the ROSAT all-sky survey with a period of 8.391 s. Its spectrum is well characterized by a black-body with a temperature of $8 \times 10^5$ K. We propose that the radiation from this object is thermal emission from a cooling neutron star. For this black-body temperature we can obtain a robust estimate of the object's age of $\sim 3 \times 10^5$ yr, yielding a polar field $\sim 10^{14}$ G for magnetic-dipole spin down and a value of ${\dot P}$ compatible with current observations.

Using recently calculated analytic and numerical models for the thermal structure of ultramagnetized neutron stars, we estimate the effects that ultrastrong magnetic fields $B \ge 10^{14}$ G have on the thermal evolution of a neutron star. Understanding this evolution is necessary to interpret models that invoke ``magnetars'' to account for soft $\gamma$-ray emission from some repeating sources.

Using recently calculated analytic models for the thermal structure of ultramagnetized neutron stars, we estimate the thermal fluxes from young ($t\sim 1000$ yr) ultramagnetized ($B \sim 10^{15}$ G) cooling neutron stars. We find that the pulsed X-ray emission from objects such as 1E 1841-045 and 1E 2259+586 as well as many soft-gamma repeaters can be explained by photon cooling if the neutron star possesses a thin insulating envelope of matter of low atomic weight at densities $\rho < 10^{7}-10^{8}$ g/cm$^3$. The total mass of this insulating layer is $M \sim 10^{-11}-10^{-8} M_\odot$.

We consider the field equations of a static magnetic field including one-loop QED corrections, and calculate the corrections to the field of a magnetic dipole. PACS: 12.20.Ds, 97.60.Jd, 97.60.Gb

We use an analytic form for the Heisenberg-Euler Lagrangian to calculate the birefringent and dichroic properties of the vacuum for arbitrarily strong wrenchless fields. PACS : 12.20.Ds, 42.25.Lc 97.60.Jd, 98.70.Rz

We determine the evolution of the galaxy luminosity function (LF) as a function of spectral type using the Autofib redshift survey, a compendium of over 1700 galaxy redshifts in various magnitude-limited samples spanning b_J=11.5-24.0. To carry out this analysis we have developed a cross-correlation technique which classifies faint galaxy spectra into one of six types based on local galaxy templates. Tests and simulations show that this technique yields classifications correct to within one type for more than 90% of the galaxies in our sample. We have also developed extensions of the step-wise maximum likelihood method and the STY parametric method for estimating LFs which are applicable to recovering an evolving LF from multiple samples. We find that: (i) The spectra and LF of E/S0 galaxies show no appreciable evolution out to at least z ~ 0.5. (ii) Early-type spirals show modest evolution, characterised by a gradual steepening of the faint end of their LF with redshift. (iii) Out to z ~ 0.5, the overall evolution of the galaxy population is dominated by changes seen in late-type spirals. The characteristic luminosity (L^*) of these galaxies appears to brighten with redshift and there are signs of strong density evolution (a rapid increase in \phi^*). These effects appear to be luminosity dependent so that the LF steepens at higher redshift. These trends are accompanied by a steep increase in the median [OII] equivalent width, implying a rapid increase in the star-formation rate with redshift at fixed luminosity---a given star-formation rate is found at higher redshift in galaxies of higher luminosity. We find that these conclusions are robust with respect to spectral classification errors and the luminosity function estimator.

We calculate the reaction cross-sections for the fusion of hydrogen and deuterium in strong magnetic fields as are believed to exist in the atmospheres of neutron stars. We find that in the presence of a strong magnetic field ($B \gsim 10^{12}$G), the reaction rates are many orders of magnitude higher than in the unmagnetized case. The fusion of both protons and deuterons are important over a neutron star's lifetime for ultrastrong magnetic fields ($B \sim 10^{16}$G). The enhancement may have dramatic effects on thermonuclear runaways and bursts on the surfaces of neutron stars.

We derive an analytic form for the Heisenberg-Euler Lagrangian in the limit where the component of the electric field parallel to the magnetic field is small. We expand these analytic functions to all orders in the field strength ($F_{\mu\nu}F^{\mu\nu}$) in the limits of weak and strong fields, and use these functions to estimate the pair-production rate in arbitrarily strong electric fields and the photon-splitting rate in arbitrarily strong magnetic fields.

We present a detailed determination of the restframe B-band galaxy luminosity function (LF) as a function of redshift and star formation activity from z=0 to z=0.75. The dataset used for this purpose is a combined sample of over 1700 redshifts spanning a wide range in apparent magnitude, 11.5<B<24.0, which we term the Autofib Redshift Survey. The sample includes various earlier magnitude-limited surveys constructed by our team as well as a new survey of 1026 redshifts measured for galaxies at intermediate magnitudes. The large range in apparent magnitude sampled allows us to investigate both the nature of the LF at low redshift (z<0.1) and possible evolution in its shape to z=0.75. We find that earlier bright surveys have underestimated the absolute normalisation of the LF. Because the shape of the local LF does not change with the survey apparent magnitude limit, it seems unlikely that the local deficiency arises from an underestimated population of low luminosity galaxies. Furthermore, surface brightness losses cannot be significant unless they conspire to retain the LF shape over a variety of detection thresholds. Our data directly demonstrates that the B-band LF evolves with redshift. This evolution is best represented as a steepening of the faint-end slope of the LF, from alpha=-1.1 at low redshift to alpha=-1.5 at z=0.5. Using [OII] emission as an indicator of star formation activity, we show that the LF of quiescent galaxies has remained largely unchanged since z=0.5, whereas the luminosity density of star-forming galaxies has declined by nearly a factor of 2. The steepening of the overall LF with lookback time is a direct consequence of the increasing space density of blue star-forming galaxies at moderate redshifts.

We predict the observable properties of the galaxy population in popular hierarchical models of galaxy formation. We employ a detailed semianalytic procedure which incorporates the formation and merging of dark matter halos, the shock heating and radiative cooling of gas, self-regulated star formation, the merging of galaxies within dark matter halos, and the spectral evolution of the stellar populations. We contrast the standard CDM cosmogony with variants of the CDM model having either a low value of H_0, or a low value of Omega with or without a cosmological constant. In addition, we compare galaxy formation in these CDM universes with a CHDM model. We find that although the models have some success in remedying the shortcomings of the standard CDM cosmogony, none of these new models produce broad agreement with the whole range of observations. Although the low-Omega and Omega+Lambda=1 CDM models reduce the discrepancy between the predicted and observed Tully-Fisher relations (the main weakness of galaxy formation in standard CDM), these models predict an inverted colour-magnitude relation and do not produce an exponential cut-off at the bright end of the galaxy luminosity function. All of our models predict recent star formation and exhibit galaxy colours bluer than observed, but this problem is far more severe in the CHDM model which produces colours about two magnitudes too blue in B-K. Unlike in the variants of the CDM model in the CHDM case this result is not dependent on our model of stellar feedback, but is instead directly caused by the late epoch of structure formation in this model.