The population of Milky Way satellite galaxies is of great interest for cosmology, fundamental physics, and astrophysics. They represent the faint end of the galaxy luminosity function, are the most dark-matter dominated objects in the local Universe, and contain the oldest and most metal-poor stellar populations. Recent surveys have revealed around 60 satellites, but this could represent less than half of the total. Characterization of these systems remains a challenge due to their low luminosity. We consider the gravitational wave observatory LISA as a potential tool for studying these satellites through observations of their short-period double white dwarf populations. LISA will observe the entire sky without selection effects due to dust extinction, complementing optical surveys, and could potentially discover massive satellites hidden behind the disk of the galaxy.

The contour method is a new approach to calculating the non-adiabatic pulsation frequencies of stars. These frequencies can be found by solving for the complex roots of a characteristic equation constructed from the linear non-adiabatic stellar pulsation equations. A complex-root solver requires an initial trial frequency for each non-adiabatic root. A standard method for obtaining initial trial frequencies is to use a star's adiabatic pulsation frequencies, but this method can fail to converge to non-adiabatic roots, especially as the growth and/or damping rate of the pulsations becomes large. The contour method provides an alternative way for obtaining initial trial frequencies that robustly converges to non-adiabatic roots, even for stellar models with extremely non-adiabatic pulsations and thus large growth/damping rates. We describe the contour method implemented in the GYRE stellar pulsation code and use it to calculate the non-adiabatic pulsation frequencies of $10 \, \rm M_\odot$ and $20 \,\rm M_\odot$ $\beta$ Cephei models, and of a $0.9 \, \rm M_\odot$ extreme helium model.

The label switching problem arises in the Bayesian analysis of models containing multiple indistinguishable parameters with arbitrary ordering. Any permutation of these parameters is equivalent, therefore models with many such parameters have extremely multi-modal posterior distributions. It is difficult to sample efficiently from such posteriors. This paper discusses a solution to this problem which involves carefully mapping the input parameter space to a high dimensional hypertriangle. It is demonstrated that this solution is efficient even for large numbers of parameters and can be easily applied alongside any stochastic sampling algorithm. This method is illustrated using two example problems from the field of gravitational wave astronomy.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

The Tayler instability (TI) is a non-axisymmetric linear instability of an axisymmetric toroidal magnetic field in magneto-hydrostatic equilibrium (MHSE). Spruit (1999, 2002) has proposed that in a differentially rotating radiative region of a star, the TI drives a dynamo which generates magnetic fields that can efficiently transport angular momentum; a parameterized version of this dynamo has been implemented in stellar structure and evolution codes and shown to be important for determining interior spin. Numerical simulations, however, have yet to definitively demonstrate the operation of the dynamo. A criterion for the MHSE to develop the TI was derived using fully-compressible magneto-hydrodynamics, while numerical simulations of dynamical processes in stars frequently use an anelastic approximation. This motivates us to derive a new anelastic Tayler instability (anTI) criterion. We find that some MHSE configurations unstable in the fully-compressible case, become stable in the anelastic case. We find and characterize the unstable modes of a simple family of cylindrical MHSE configurations using numerical calculations, and discuss the implications for fully non-linear anelastic simulations.

Super-massive black hole binaries are expected to produce a GW signal in the nano-Hertz frequency band which may be detected by PTAs in the coming years. The signal is composed of both stochastic and individually resolvable components. Here we develop a generic Bayesian method for the analysis of resolvable sources based on the construction of `null-streams' which cancel the part of the signal held in common for each pulsar (the Earth-term). For an array of $N$ pulsars there are $N-2$ independent null-streams that cancel the GW signal from a particular sky location. This method is applied to the localisation of quasi-circular binaries undergoing adiabatic inspiral. We carry out a systematic investigation of the scaling of the localisation accuracy with signal strength and number of pulsars in the PTA. Additionally, we find that source sky localisation with the International PTA data release one is vastly superior than what is achieved by its constituent regional PTAs.

We propose a novel methodology to select host galaxy candidates of future pulsar timing array (PTA) detections of resolved gravitational waves (GWs) from massive black hole binaries (MBHBs). The method exploits the physical dependence of the GW amplitude on the MBHB chirp mass and distance to the observer, together with empirical MBH mass-host galaxy correlations, to rank potential host galaxies in the mass-redshift plane. This is coupled to a null-stream based likelihood evaluation of the GW amplitude and sky position in a Bayesian framework that assigns to each galaxy a probability of hosting the MBHB generating the GW signal. We test our algorithm on a set of realistic simulations coupling the likely properties of the first PTA resolved GW signal to synthetic all-sky galaxy maps. For a foreseeable PTA sky-localization precision of 100 squared degrees, we find that the GW source is hosted with 50% (90%) probability within a restricted number of <50 (<500) potential hosts. These figures are orders of magnitude smaller than the total number of galaxies within the PTA sky error-box, enabling extensive electromagnetic follow-up campaigns on a limited number of targets.

Context. While rotation has a major impact on stellar structure and evolution, its effects are not well understood. Thanks to high- quality and long timebase photometric observations obtained with recent space missions, we are now able to study stellar rotation more precisely. Aims. We aim to constrain radial differential rotation profiles in gamma Doradus (gamma Dor) stars, and to develop new theoretical seismic diagnosis for such stars with rapid and potentially non-uniform rotation. Methods. We derive a new asymptotic description which accounts for the impact of weak differential near-core rotation on gravity- mode period spacings. The theoretical predictions are illustrated from pulsation computations with the code GYRE and compared with observations of gamma Dor stars. When possible, we also derive the surface rotation rates in these stars by detecting and analysing signatures of rotational modulation, and compute the core-to-surface rotation ratios. Results. Stellar rotation has to be strongly differential before its effects on period spacing patterns can be detected, unless multiple period spacing patterns can be compared. Six stars in our sample exhibit a single unexplained period spacing pattern of retrograde modes. We hypothesise that these are Yanai modes. Finally, we find signatures of rotational spot modulation in the photometric data of eight targets. Conclusions. If only one period spacing pattern is detected and analysed for a star, it is difficult to detect differential rotation. A rigidly rotating model will often provide the best solution. Differential rotation can only be detected when multiple period spacing patterns have been found for a single star or its surface rotation rate is known as well. This is the case for eight stars in our sample, revealing surface-to-core rotation ratios between 0.95 and 1.05.

Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly-pulsating B (SPB) stars can significantly modify the stars' internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the GYRE oscillation code and the MESASTAR stellar evolution code. Focusing first on a $4.21$ $M_\odot$ model, we simulate 1,000 years of stellar evolution under the combined effects of the torque due to a single unstable prograde g mode (with an amplitude chosen on the basis of observational constraints), and diffusive angular momentum transport due to convection, overshooting, and rotational instabilities. We find that the g mode rapidly extracts angular momentum from the surface layers, depositing it deeper in the stellar interior. The angular momentum transport is so efficient that by the end of the simulation the initially non-rotating surface layers are spun in the retrograde direction to $\approx30\%$ of the critical rate. However, the additional inclusion of magnetic stresses in our simulations, almost completely inhibits this spin-up. Expanding our simulations to cover the whole instability strip, we show that the same general behavior is seen in all SPB stars. After providing some caveats to contextualize our results, we hypothesize that the observed slower surface rotation of SPB stars (as compared to other B-type stars) may be the direct consequence of the angular momentum transport that our simulations demonstrate.

White dwarfs (WDs) can be tidally disrupted only by massive black holes (MBHs) with masses less than $\sim10^5 M_\odot$. These tidal interactions feed material to the MBH well above its Eddington limit, with the potential to launch a relativistic jet. The corresponding beamed emission is a promising signpost to an otherwise quiescent MBH of relatively low mass. We show that the mass transfer history, and thus the lightcurve, are quite different when the disruptive orbit is parabolic, eccentric, or circular. The mass lost each orbit exponentiates in the eccentric-orbit case leading to the destruction of the WD after several tens of orbits. We examine the stellar dynamics of clusters surrounding MBHs to show that single-passage WD disruptions are substantially more common than repeating encounters. The $10^{49}$ erg s$^{-1}$ peak luminosity of these events makes them visible to cosmological distances. They may be detectible at rates of as many as tens per year by instruments like Swift. In fact, WD-disruption transients significantly outshine their main-sequence star counterparts, and are the most likely tidal interaction to be detected arising from MBHs with masses less than $10^5 M_\odot$. The detection or non-detection of such WD-disruption transients by Swift is, therefore, a powerful tool to constrain lower end of the MBH mass function. The emerging class of ultra-long gamma ray bursts all have peak luminosities and durations reminiscent of WD disruptions, offering a hint that WD-disruption transients may already be present in existing datasets.

We have studied cloudy taenite, metallographic cooling rates, and shock effects in 30 H3-6 chondrites to elucidate the thermal and early impact history of the H chondrite parent body. We focused on H chondrites with Ar-Ar ages greater than 4.4 Gyr and unshocked and mildly shocked H chondrites, as strongly shocked chondrites with such old ages are very rare. Cooling rates for most H chondrites at 500 C are 10-50 C/Myr and do not decrease systematically with increasing petrologic type as predicted by the onion-shell model in which types 3 to 5 are arranged in concentric layers around a type 6 core. Some type 4 chondrites cooled slower than some type 6 chondrites and type 3 chondrites did not cool faster than other types, contrary to the onion-shell model. Cloudy taenite particle sizes, which range from 40 to 120 nm, are inversely correlated with metallographic cooling rates and show that the latter were not compromised by shock heating. The three H4 chondrites that were used to develop the onion-shell model, Ste. Marguerite, Beaver Creek, and Forest Vale, cooled through 500 C at ~5000 C/Myr. Our thermal modeling shows that these rates are 50 times higher than could be achieved in a body that was heated by 26Al and cooled without disturbance by impact. Published Ar-Ar ages do not decrease systematically with increasing petrologic type but do correlate inversely with cloudy taenite particle size suggesting that impact mixing decreased during metamorphism. Metal and silicate compositions in regolith breccias show that impacts mixed material after metamorphism without causing significant heating. Impacts during metamorphism created Portales Valley and two other H6 chondrites with large metallic veins, excavated the fast-cooled H4 chondrites around 3-4 Myr after accretion, and mixed petrologic types.

We present initial results of the first panoramic search for high-amplitude near-infrared variability in the Galactic Plane. We analyse the widely separated two-epoch K-band photometry in the 5th and 7th data releases of the UKIDSS Galactic Plane Survey. We find 45 stars with \Delta K > 1 mag, including 2 previously known OH/IR stars and a Nova. Even though the mid-plane is not yet included in the dataset, we find the majority (66%) of our sample to be within known star forming regions (SFRs), with two large concentrations in the Serpens OB2 association (11 stars) and the Cygnus-X complex (12 stars). Sources in SFRs show spectral energy distributions (SEDs) that support classification as Young Stellar Objects (YSOs). This indicates that YSOs dominate the Galactic population of high amplitude infrared variable stars at low luminosities and therefore likely dominate the total high amplitude population. Spectroscopic follow up of the DR5 sample shows at least four stars with clear characteristics of eruptive pre-main-sequence variables, two of which are deeply embedded. Our results support the recent concept of eruptive variability comprising a continuum of outburst events with different timescales and luminosities, but triggered by a similar physical mechanism involving unsteady accretion. Also, we find what appears to be one of the most variable classical Be stars.

GRB 110918A is the brightest long GRB detected by Konus-WIND during its 19 years of continuous observations and the most luminous GRB ever observed since the beginning of the cosmological era in 1997. We report on the final IPN localization of this event and its detailed multiwavelength study with a number of space-based instruments. The prompt emission is characterized by a typical duration, a moderare $E_{peak}$ of the time-integrated spectrum, and strong hard-to-soft evolution. The high observed energy fluence yields, at z=0.984, a huge isotropic-equivalent energy release $E_{iso}=(2.1\pm0.1)\times10^{54}$ erg. The record-breaking energy flux observed at the peak of the short, bright, hard initial pulse results in an unprecedented isotropic-equivalent luminosity $L_{iso}=(4.7\pm0.2)\times10^{54}$erg s$^{-1}$. A tail of the soft gamma-ray emission was detected with temporal and spectral behavior typical of that predicted by the synchrotron forward-shock model. Swift/XRT and Swift/UVOT observed the bright afterglow from 1.2 to 48 days after the burst and revealed no evidence of a jet break. The post-break scenario for the afterglow is preferred from our analysis, with a hard underlying electron spectrum and ISM-like circumburst environment implied. We conclude that, among multiple reasons investigated, the tight collimation of the jet must have been a key ingredient to produce this unusually bright burst. The inferred jet opening angle of 1.7-3.4 deg results in reasonable values of the collimation-corrected radiated energy and the peak luminosity, which, however, are still at the top of their distributions for such tightly collimated events. We estimate a detection horizon for a similar ultraluminous GRB of $z\sim7.5$ for Konus-WIND, and $z\sim12$ for Swift/BAT, which stresses the importance of GRBs as probes of the early Universe.

In this paper, we present results from the complete set of cosmic microwave background (CMB) radiation temperature anisotropy observations made with the Arcminute Cosmology Bolometer Array Receiver (ACBAR) operating at 150 GHz. We include new data from the final 2005 observing season, expanding the number of detector-hours by 210% and the sky coverage by 490% over that used for the previous ACBAR release. As a result, the band-power uncertainties have been reduced by more than a factor of two on angular scales encompassing the third to fifth acoustic peaks as well as the damping tail of the CMB power spectrum. The calibration uncertainty has been reduced from 6% to 2.1% in temperature through a direct comparison of the CMB anisotropy measured by ACBAR with that of the dipole-calibrated WMAP5 experiment. The measured power spectrum is consistent with a spatially flat, LambdaCDM cosmological model. We include the effects of weak lensing in the power spectrum model computations and find that this significantly improves the fits of the models to the combined ACBAR+WMAP5 power spectrum. The preferred strength of the lensing is consistent with theoretical expectations. On fine angular scales, there is weak evidence (1.1 sigma) for excess power above the level expected from primary anisotropies. We expect any excess power to be dominated by the combination of emission from dusty protogalaxies and the Sunyaev-Zel'dovich effect (SZE). However, the excess observed by ACBAR is significantly smaller than the excess power at ell > 2000 reported by the CBI experiment operating at 30 GHz. Therefore, while it is unlikely that the CBI excess has a primordial origin; the combined ACBAR and CBI results are consistent with the source of the CBI excess being either the SZE or radio source contamination.

We present the first installment of HI sources extracted from the Arecibo Legacy Fast ALFA (ALFALFA) extragalactic survey, initiated in 2005. Sources have been extracted from 3-D spectral data cubes and then examined interactively to yield global HI parameters. A total of 730 HI detections are catalogued within the solid angle 11h44m < R.A.(J2000) < 14h00m and +12deg < Dec.(J2000) < +16deg, and redshift range -1600 \kms < cz < 18000 \kms. In comparison, the HI Parkes All-Sky Survey (HIPASS) detected 40 HI signals in the same region. Optical counterparts are assigned via examination of digital optical imaging databases. ALFALFA HI detections are reported for three distinct classes of signals: (a) detections, typically with S/N > 6.5; (b) high velocity clouds in the Milky Way or its periphery; and (c) signals of lower S/N (to ~ 4.5) which coincide spatially with an optical object of known similar redshift. Although this region of the sky has been heavily surveyed by previous targeted observations based on optical flux-- or size-- limited samples, 69% of the extracted sources are newly reported HI detections. The resultant positional accuracy of HI sources is 20" (median). The median redshift of the sample is ~7000 \kms and its distribution reflects the known local large scale structure including the Virgo cluster. Several extended HI features are found in the vicinity of the Virgo cluster. A small percentage (6%) of HI detections have no identifiable optical counterpart, more than half of which are high velocity clouds in the Milky Way vicinity; the remaining 17 objects do not appear connected to or associated with any known galaxy.

We report improved measurements of temperature anisotropies in the cosmic microwave background (CMB) radiation made with the Arcminute Cosmology Bolometer Array Receiver (ACBAR). In this paper, we use a new analysis technique and include 30% more data from the 2001 and 2002 observing seasons than the first release to derive a new set of band-power measurements with significantly smaller uncertainties. The planet-based calibration used previously has been replaced by comparing the flux of RCW38 as measured by ACBAR and BOOMERANG to transfer the WMAP-based BOOMERANG calibration to ACBAR. The resulting power spectrum is consistent with the theoretical predictions for a spatially flat, dark energy dominated LCDM cosmology including the effects of gravitational lensing. Despite the exponential damping on small angular scales, the primary CMB fluctuations are detected with a signal-to-noise ratio of greater than 4 up to multipoles of l=2000. This increase in the precision of the fine-scale CMB power spectrum leads to only a modest decrease in the uncertainties on the parameters of the standard cosmological model. At high angular resolution, secondary anisotropies are predicted to be a significant contribution to the measured anisotropy. A joint analysis of the ACBAR results at 150 GHz and the CBI results at 30 GHz in the multipole range 2000 < l < 3000 shows that the power, reported by CBI in excess of the predicted primary anisotropy, has a frequency spectrum consistent with the thermal Sunyaev-Zel'dovich effect and inconsistent with primary CMB. The results reported here are derived from a subset of the total ACBAR data set; the final ACBAR power spectrum at 150 GHz will include 3.7 times more effective integration time and 6.5 times more sky coverage than is used here.

We describe the Arcminute Cosmology Bolometer Array Receiver (ACBAR); a multifrequency millimeter-wave receiver designed for observations of the Cosmic Microwave Background (CMB) and the Sunyaev-Zel'dovich effect in clusters of galaxies. The ACBAR focal plane consists of a 16-pixel, background-limited, 240 mK bolometer array that can be configured to observe simultaneously at 150, 220, 280, and 350 GHz. With 4-5' FWHM Gaussian beam sizes and a 3 degree azimuth chop, ACBAR is sensitive to a wide range of angular scales. ACBAR was installed on the 2 m Viper telescope at the South Pole in January 2001. We describe the design of the instrument and its performance during the 2001 and 2002 observing seasons.

We present detections of the Sunyaev-Zeldovich Effect (SZE) at 150GHz and 275GHz for the X-ray luminous z=0.299 cluster 1E0657-67. These observations were obtained as part of an X-ray, weak lensing, and SZE survey of nearby X-ray bright clusters. The SZE maps were made with the ACBAR (150, 210, 275 GHz) bolometer array installed at the Viper telescope located at the South Pole. We also present preliminary results from a blind SZE cluster survey.

During 2001 and 2002, observations of several X-ray bright clusters of galaxies were conducted using the ACBAR Bolometer Array on the South Pole Viper telescope. A multi-frequency analysis of these clusters is currently underway. This multi-frequency analysis includes 150, 220 and 275 GHz data from ACBAR, X-ray imaging and spectroscopy from Chandra and XMM-Newton, and weak lensing data from the CTIO 4m Blanco telescope. We describe here how ACBAR can be used to create fully sampled cluster images and present such images for four of the clusters in our sample; Abell 3266, Abell 3827, Abell S1063 and 1E0657-56. In these images, the Sunyaev-Zel'dovich Effect is clearly detected.

We review the first science results from the Arcminute Cosmology Bolometer Array Receiver (ACBAR); a multifrequency millimeter-wave receiver optimized for observations of the Cosmic Microwave Background (CMB) and the Sunyaev-Zel'dovich (SZ) effect in clusters of galaxies. ACBAR was installed on the 2 m Viper telescope at the South Pole in January 2001 and the results presented here incorporate data through July 2002. We present the power spectrum of the CMB at 150 GHz over the range \ell = 150 - 3000 measured by ACBAR as well as estimates for the values of the cosmological parameters within the context of Lambda-CDM models. We find that the inclusion of Omega_Lambda greatly improves the fit to the power spectrum. We also observe a slight excess of small-scale anisotropy at 150 GHz; if interpreted as power from the SZ effect of unresolved clusters, the measured signal is consistent with CBI and BIMA within the context of the SZ power spectrum models tested.

As part of a combined Sunyaev-Zeldovich Effect (SZE), X-ray and weak lensing survey of low redshift (z<0.1) X-ray clusters, we present SZE images of the z=0.059 X-ray cluster Abell 3266 at three observing frequencies (150, 220, 275 GHz) and after the spectral subtraction of primary Cosmic Microwave Background (CMB) anisotropies. These images were generated using the ACBAR bolometer array operated on the Viper telescope at the South Pole. The multi-frequency data from ACBAR should allow us to overcome one of the main obstacles facing the analysis of SZE observations of nearby clusters, i.e. contamination from primary Cosmic Microwave Background (CMB) anisotropies.

We report an investigation of cosmological parameters based on the measurements of anisotropy in the cosmic microwave background radiation (CMB) made by ACBAR. We use the ACBAR data in concert with other recent CMB measurements to derive Bayesian estimates of parameters in inflation-motivated adiabatic cold dark matter models. We apply a series of additional cosmological constraints on the shape and amplitude of the density power spectrum, the Hubble parameter and from supernovae to further refine our parameter estimates. Previous estimates of parameters are confirmed, with sensitive measurements of the power spectrum now ranging from \ell \sim 3 to 2800. Comparing individual best model fits, we find that the addition of \Omega_\Lambda as a parameter dramatically improves the fits. We also use the high-\ell data of ACBAR, along with similar data from CBI and BIMA, to investigate potential secondary anisotropies from the Sunyaev-Zeldovich effect. We show that the results from the three experiments are consistent under this interpretation, and use the data, combined and individually, to estimate \sigma_8 from the Sunyaev-Zeldovich component.

We report the first measurements of anisotropy in the cosmic microwave background (CMB) radiation with the Arcminute Cosmology Bolometer Array Receiver (ACBAR). The instrument was installed on the 2.1m Viper telescope at the South Pole in January 2001; the data presented here are the product of observations up to and including July 2002. The two deep fields presented here, have had offsets removed by subtracting lead and trail observations and cover approximately 24 deg^2 of sky selected for low dust contrast. These results represent the highest signal to noise observations of CMB anisotropy to date; in the deepest 150GHz band map, we reached an RMS of 8.0\mu K per 5' beam. The 3 degree extent of the maps, and small beamsize of the experiment allow the measurement of the CMB anisotropy power spectrum over the range \ell = 150-3000 with resolution of \Delta \ell=150. The contributions of galactic dust and radio sources to the observed anisotropy are negligible and are removed in the analysis. The resulting power spectrum is found to be consistent with the primary anisotropy expected in a concordance \Lambda CDM Universe.