We present a search for anisotropic cosmic birefringence in 500 deg$^2$ of southern sky observed at 150 GHz with the SPTpol camera on the South Pole Telescope. We reconstruct a map of cosmic polarization rotation anisotropies using higher-order correlations between the observed cosmic microwave background (CMB) $E$ and $B$ fields. We then measure the angular power spectrum of this map, which is found to be consistent with zero. The non-detection is translated into an upper limit on the amplitude of the scale-invariant cosmic rotation power spectrum, $L(L+1)C_L^{\alpha\alpha}/2\pi < 0.10 \times 10^{-4}$ rad$^2$ (0.033 deg$^2$, 95% C.L.). This upper limit can be used to place constraints on the strength of primordial magnetic fields, $B_{1 \rm Mpc} < 17 {\rm nG} $ (95% C.L.), and on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma} < 4.0 \times 10^{-2}/H_I $ (95% C.L.), where $H_I$ is the inflationary Hubble scale. For the first time, we also cross-correlate the CMB temperature fluctuations with the reconstructed rotation angle map, a signal expected to be non-vanishing in certain theoretical scenarios, and find no detectable signal. We perform a suite of systematics and consistency checks and find no evidence for contamination.

We report a B-mode power spectrum measurement from the cosmic microwave background (CMB) polarization anisotropy observations made using the SPTpol instrument on the South Pole Telescope. This work uses 500 deg$^2$ of SPTpol data, a five-fold increase over the last SPTpol B-mode release. As a result, the bandpower uncertainties have been reduced by more than a factor of two, and the measurement extends to lower multipoles: $52 < \ell < 2301$. Data from both 95 and 150 GHz are used, allowing for three cross-spectra: 95 GHz x 95 GHz, 95 GHz x 150 GHz, and 150 GHz x 150 GHz. B-mode power is detected at very high significance; we find $P(BB < 0) = 5.8 \times 10^{-71}$, corresponding to a $18.1 \sigma$ detection of power. An upper limit is set on the tensor-to-scalar ratio, $r < 0.44$ at 95% confidence (the expected $1 \sigma$ constraint on $r$ given the measurement uncertainties is 0.22). We find the measured B-mode power is consistent with the Planck best-fit $\Lambda$CDM model predictions. Scaling the predicted lensing B-mode power in this model by a factor Alens, the data prefer Alens = $1.17 \pm 0.13$. These data are currently the most precise measurements of B-mode power at $\ell > 320$.

We report new measurements of millimeter-wave power spectra in the angular multipole range $2000 \le \ell \le 11,000$ (angular scales $5^\prime \gtrsim \theta \gtrsim 1^\prime$). By adding 95 and 150\,GHz data from the low-noise 500 deg$^2$ SPTpol survey to the SPT-SZ three-frequency 2540 deg$^2$ survey, we substantially reduce the uncertainties in these bands. These power spectra include contributions from the primary cosmic microwave background, cosmic infrared background, radio galaxies, and thermal and kinematic Sunyaev-Zel'dovich (SZ) effects. The data favor a thermal SZ (tSZ) power at 143\,GHz of $D^{\rm tSZ}_{3000} = 3.42 \pm 0.54~ \mu {\rm K}^2$ and a kinematic SZ (kSZ) power of $D^{\rm kSZ}_{3000} = 3.0 \pm 1.0~ \mu {\rm K}^2$. This is the first measurement of kSZ power at $\ge 3\,\sigma$. We study the implications of the measured kSZ power for the epoch of reionization, finding the duration of reionization to be $\Delta z_{re} = 1.0^{+1.6}_{-0.7}$ ($\Delta z_{re}< 4.1$ at 95% confidence), when combined with our previously published tSZ bispectrum measurement.

We present cosmological constraints based on the cosmic microwave background (CMB) lensing potential power spectrum measurement from the recent 500 deg$^2$ SPTpol survey, the most precise CMB lensing measurement from the ground to date. We fit a flat $\Lambda$CDM model to the reconstructed lensing power spectrum alone and in addition with other data sets: baryon acoustic oscillations (BAO) as well as primary CMB spectra from Planck and SPTpol. The cosmological constraints based on SPTpol and Planck lensing band powers are in good agreement when analysed alone and in combination with Planck full-sky primary CMB data. With weak priors on the baryon density and other parameters, the CMB lensing data alone provide a 4\% constraint on $\sigma_8\Omega_m^{0.25} = 0.0593 \pm 0.025$.. Jointly fitting with BAO data, we find $\sigma_8=0.779 \pm 0.023$, $\Omega_m = 0.368^{+0.032}_{-0.037}$, and $H_0 = 72.0^{+2.1}_{-2.5}\,\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1} $, up to $2\,\sigma$ away from the central values preferred by Planck lensing + BAO. However, we recover good agreement between SPTpol and Planck when restricting the analysis to similar scales. We also consider single-parameter extensions to the flat $\Lambda$CDM model. The SPTpol lensing spectrum constrains the spatial curvature to be $\Omega_K = -0.0007 \pm 0.0025$ and the sum of the neutrino masses to be $\sum m_{\nu} < 0.23$ eV at 95\% C.L. (with Planck primary CMB and BAO data), in good agreement with the Planck lensing results. With the differences in the $S/N$ of the lensing modes and the angular scales covered in the lensing spectra, this analysis represents an important independent check on the full-sky Planck lensing measurement.

We study the polarisation properties of extragalactic sources at 95 and 150 GHz in the SPTpol 500 deg$^2$ survey. We estimate the polarised power by stacking maps at known source positions, and correct for noise bias by subtracting the mean polarised power at random positions in the maps. We show that the method is unbiased using a set of simulated maps with similar noise properties to the real SPTpol maps. We find a flux-weighted mean-squared polarisation fraction $\langle p^2 \rangle= [8.9\pm1.1] \times 10^{-4}$ at 95 GHz and $[6.9\pm1.1] \times 10^{-4}$ at 150~GHz for the full sample. This is consistent with the values obtained for a sub-sample of active galactic nuclei. For dusty sources, we find 95 per cent upper limits of $\langle p^2 \rangle_{\rm 95}<16.9 \times 10^{-3}$ and $\langle p^2 \rangle_{\rm 150}<2.6 \times 10^{-3}$. We find no evidence that the polarisation fraction depends on the source flux or observing frequency. The 1-$\sigma$ upper limit on measured mean squared polarisation fraction at 150 GHz implies that extragalactic foregrounds will be subdominant to the CMB E and B mode polarisation power spectra out to at least $\ell\lesssim5700$ ($\ell\lesssim4700$) and $\ell\lesssim5300$ ($\ell\lesssim3600$), respectively at 95 (150) GHz.

We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel'dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confirmed 81 candidates as galaxy clusters. We use these follow-up images and archival images to estimate photometric redshifts for 66 galaxy clusters and spectroscopic observations to obtain redshifts for 13 systems. An additional 2 galaxy clusters are confirmed using the overdensity of near-infrared galaxies only, and are presented without redshifts. We find that 15 candidates (18% of the total sample) are at redshift of $z \geq 1.0$, with a maximum confirmed redshift of $z_{\rm{max}} = 1.38 \pm 0.10$. We expect this catalog to contain every galaxy cluster with $M_{500c} > 2.6 \times 10^{14} M_\odot h^{-1}_{70}$ and $z > 0.25$ in the survey area. The mass threshold is approximately constant above $z = 0.25$, and the complete catalog has a median mass of approximately $ M_{500c} = 2.7 \times 10^{14} M_\odot h^{-1}_{70}$. Compared to previous SPT works, the increased depth of the millimeter-wave data (11.2 and 6.5 $\mu$K-arcmin at 95 and 150 GHz, respectively) makes it possible to find more galaxy clusters at high redshift and lower mass.

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $\xi \ge 5$ and an additional 204 systems at $4<\xi<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $\xi>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($\lambda-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $\sigma$ level---with the relations intersecting at $\lambda=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

We present a measurement of the cosmic microwave background (CMB) lensing potential using 500 deg$^2$ of 150 GHz data from the SPTpol receiver on the South Pole Telescope. The lensing potential is reconstructed with signal-to-noise per mode greater than unity at lensing multipoles $L \lesssim 250$, using a quadratic estimator on a combination of CMB temperature and polarization maps. We report measurements of the lensing potential power spectrum in the multipole range of $100< L < 2000$ from sets of temperature-only, polarization-only, and minimum-variance estimators. We measure the lensing amplitude by taking the ratio of the measured spectrum to the expected spectrum from the best-fit $\Lambda$CDM model to the $\textit{Planck}$ 2015 TT+lowP+lensing dataset. For the minimum-variance estimator, we find $A_{\rm{MV}} = 0.944 \pm 0.058{\rm (Stat.)}\pm0.025{\rm (Sys.)}$; restricting to only polarization data, we find $A_{\rm{POL}} = 0.906 \pm 0.090 {\rm (Stat.)} \pm 0.040 {\rm (Sys.)}$. Considering statistical uncertainties alone, this is the most precise polarization-only lensing amplitude constraint to date (10.1 $\sigma$), and is more precise than our temperature-only constraint. We perform null tests and consistency checks and find no evidence for significant contamination.

We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $\lambda \ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8\sigma$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements.

We present measurements of the $E$-mode polarization angular auto-power spectrum ($EE$) and temperature-$E$-mode cross-power spectrum ($TE$) of the cosmic microwave background (CMB) using 150 GHz data from three seasons of SPTpol observations. We report the power spectra over the spherical harmonic multipole range $50 < \ell \leq 8000$, and detect nine acoustic peaks in the $EE$ spectrum with high signal-to-noise ratio. These measurements are the most sensitive to date of the $EE$ and $TE$ power spectra at $\ell > 1050$ and $\ell > 1475$, respectively. The observations cover 500 deg$^2$, a fivefold increase in area compared to previous SPTpol analyses, which increases our sensitivity to the photon diffusion damping tail of the CMB power spectra enabling tighter constraints on \LCDM model extensions. After masking all sources with unpolarized flux $>50$ mJy we place a 95% confidence upper limit on residual polarized point-source power of $D_\ell = \ell(\ell+1)C_\ell/2\pi <0.107\,\mu{\rm K}^2$ at $\ell=3000$, suggesting that the $EE$ damping tail dominates foregrounds to at least $\ell = 4050$ with modest source masking. We find that the SPTpol dataset is in mild tension with the $\Lambda CDM$ model ($2.1\,\sigma$), and different data splits prefer parameter values that differ at the $\sim 1\,\sigma$ level. When fitting SPTpol data at $\ell < 1000$ we find cosmological parameter constraints consistent with those for $Planck$ temperature. Including SPTpol data at $\ell > 1000$ results in a preference for a higher value of the expansion rate ($H_0 = 71.3 \pm 2.1\,\mbox{km}\,s^{-1}\mbox{Mpc}^{-1}$ ) and a lower value for present-day density fluctuations ($\sigma_8 = 0.77 \pm 0.02$).

A promising supermassive black hole seed formation channel is that of direct collapse from primordial gas clouds. We perform a suite of 3D hydrodynamics simulations of an isolated turbulent gas cloud to investigate conditions conducive to forming massive black hole seeds via direct collapse, probing the impact of cloud metallicity, gas temperature floor and cooling physics on cloud fragmentation. We find there is no threshold in metallicity which produces a sharp drop in fragmentation. When molecular cooling is not present, metallicity has little effect on fragmentation. When molecular cooling is present, fragmentation is suppressed by at most $\sim 25\%$, with the greatest suppression seen at metallicities below $2\%$ solar. A gas temperature floor $\sim 10^{4}$K produces the largest drop in fragmentation of any parameter choice, reducing fragmentation by $\sim 60\%$. At metallicities below $2\%$ solar or at temperatures $\sim 10^{3}$K we see a reduction in fragmentation $\sim 20-25 \%$. For a cloud of metallicity $2\%$ solar above and a temperature below $10^3$K, the detailed choices of temperature floor, metallicity, and cooling physics have little impact on fragmentation.

Models of modified gravity offer promising alternatives to the concordance $\Lambda$CDM cosmology to explain the late-time acceleration of the universe. A popular such model is $f(R)$ gravity, in which the Ricci scalar in the Einstein-Hilbert action is replaced by a general function of it. We study the $f(R)$ model of Hu & Sawicki (2007), which recovers standard General Relativity in high density regimes, while reproducing the desired late-time acceleration at cosmological scales. We run a suite of high resolution zoom simulations using the ECOSMOG code to examine the effect of $f(R)$ gravity on the properties of a halo that is analogous to the Virgo cluster. We show that the velocity dispersion profiles can potentially discriminate between $f(R)$ models and $\Lambda$CDM, and provide complementary analysis of lensing signal profiles to explore the possibility to further distinguish the different $f(R)$ models. Our results confirm the techniques explored by Cabre et al. (2012) to quantify the effect of environment in the behavior of $f(R)$ gravity, and we extend them to study halo satellites at various redshifts. We find that the modified gravity effects in our models are most observable at low redshifts, and that effects are generally stronger for satellites far from the center of the main halo. We show that the screening properties of halo satellites trace very well that of dark matter particles, which means that low-resolution simulations in which subhalos are not very well resolved can in principle be used to study satellite properties. We discuss observables, particularly for halo satellites, that can potentially be used to constrain the observational viability of $f(R)$ gravity.

Metal poor globular clusters (MPGCs) are a unique probe of the early universe, in particular the reionization era. Systems of globular clusters in galaxy clusters are particularly interesting as it is in the progenitors of galaxy clusters that the earliest reionizing sources first formed. Although the exact physical origin of globular clusters is still debated, it is generally admitted that globular clusters form in early, rare dark matter peaks (Moore et al. 2006; Boley et al. 2009). We provide a fully numerical analysis of the Virgo cluster globular cluster system by identifying the present day globular cluster system with exactly such early, rare dark matter peaks. A popular hypothesis is that that the observed truncation of blue metal poor globular cluster formation is due to reionization (Spitler et al. 2012; Boley et al. 2009; Brodie & Strader 2006); adopting this view, constraining the formation epoch of MPGCs provides a complementary constraint on the epoch of reionization. By analyzing both the line of sight velocity dispersion and the surface density distribution of the present day distribution we are able to constrain the redshift and mass of the dark matter peaks. We find and quantify a dependence on the chosen line of sight of these quantities, whose strength varies with redshift, and coupled with star formation efficiency arguments find a best fitting formation mass and redshift of $\simeq 5 \times 10^8 \rm{M}_\odot$ and $z\simeq 9$. We predict $\simeq 300$ intracluster MPGCs in the Virgo cluster. Our results confirm the techniques pioneered by Moore et al. (2006) when applied to the the Virgo cluster and extend and refine the analytic results of Spitler et al. (2012) numerically.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We determine the orbits of four double degenerate systems (DDs), composed of two white dwarfs, and of two white dwarf -- M dwarf binaries. The four DDs, WD1022+050, WD1428+373, WD1824+040, and WD2032+188, show orbital periods of 1.157155(5) d, 1.15674(2) d, 6.26602(6) d and 5.0846(3) d respectively. These periods combined with estimates for the masses of the brighter component, based on their effective temperatures, allow us to constrain the masses of the unseen companions. We estimate that the upper limit for the contribution of the unseen companions to the total luminosity in the four DDs ranges between 10 and 20 per cent. In the case of the two white dwarf - M dwarf binaries, WD1042-690 and WD2009+622, we calculate the orbital parameters by fitting simultaneously the absorption line from the white dwarf and the emission core from the M-dwarf. Their orbital periods are 0.337083(1) d and 0.741226(2) d respectively. We find signatures of irradiation on the inner face of WD2009+622's companion. We calculate the masses of both components from the gravitational redshift and the mass-radius relationship for white dwarfs and find masses of 0.75 -- 0.78 Msun and 0.61 -- 0.64 Msun for WD1042-690 and WD2009+622 respectively. This indicates that the stars probably reached the asymptotic giant branch in their evolution before entering a common envelope phase. These two white dwarf - M dwarf binaries will become cataclysmic variables, although not within a Hubble time, with orbital periods below the period gap.

Although a large number of orbital periods of cataclysmic variable stars (CVs) have been measured, comparison of period and luminosity distributions with evolutionary theory is affected by strong selection effects. A test has been discovered which is independent of these selection effects and is based upon the kinematics of CVs (Kolb & Stehle, 1996). If the standard models of evolution are correct then long-period (P_orb > 5 hrs) CVs should be typically less than 1.5 Gyr old, and their line-of-sight velocity dispersion ($\sigma_\gamma$) should be small. We present results from a pilot study which indicate that this postulate is indeed true. Four long-period dwarf novae (EM Cyg, V426 Oph, SS Cyg and AH Her) were observed over a complete orbit, in order that accurate radial velocities be obtained. We find values of -1.7, 5.4, 15.4 and 1.8 km/s with uncertainties of order 3 km/s, referred to the dynamical Local Standard of Rest (LSR), leading to a dispersion of ~ 8 km/s. Calculation of a 95 per cent confidence interval gives the result 4 < \sigma_\gamma < 28 km/s compared to a prediction of 15 km/s. We also have an improved determination of mass donor spectral type, K_2 and q for the four systems.

Short period double degenerates (DDs) are close white dwarf - white dwarf binary stars which are the result of the evolution of interacting binary stars. We present the first definitive measurements of the mass ratio for two DDs, WD0136+768 and WD1204+450, and an improved measurement of the mass ratio for WD0957-666. We compare the properties of the 6 known DDs with measured mass ratios to the predictions of various theoretical models. We confirm the result that standard models for the formation of DDs do not predict sufficient DDs with mass ratios near 1. We also show that the observed difference in cooling ages between white dwarfs in DDs is a useful constraint on the initial mass ratio of the binary. A more careful analysis of the properties of the white dwarf pair WD1704+481.2 leads us to conclude that the brighter white dwarf is older than its fainter companion. This is the opposite of the usual case for DDs and is caused by the more massive white dwarf being smaller and cooling faster. The mass ratio in the sense (mass of younger star)/(mass of older star) is then 1.43+-0.06 rather than the value 0.70+-0.03 given previously.

We present 594 radial velocity measurements for 71 white dwarfs obtained during our search for binary white dwarfs and not reported elsewhere. We identify three excellent candidate binaries, which require further observations to confirm our preliminary estimates for their orbital periods, and one other good candidate. We investigate whether our data support the existence of a population of single, low mass (<~0.5 solar masses) white dwarfs (LMWDs). These stars are difficult to explain in standard models of stellar evolution. We find that a model with a mixed single/binary population is at least ~20 times more likely to explain our data than a pure binary population. This result depends on assumed period distributions for binary LMWDs, assumed companion masses and several other factors. Therefore, the evidence in favour of the existence of a population of single LMWDs is not sufficient, in our opinion, to firmly establish the existence of such a population, but does suggest that extended observations of LMWDs to obtain a more convincing result would be worthwhile .

We have used the radial velocity variations of two sdB stars previously reported to be binaries to establish their orbital periods. They are PG0940+068, (P=8.33d) and PG1247+554 (P=0.599d). The minimum masses of the unseen companions, assuming a mass of 0.5 solar masses for the sdB stars, are 0.090 +/- 0.003 solar masses for PG1247+554 and 0.63 +/- 0.02 solar masses for PG0940+068. The nature of the companions is not constrained further by our data.

We present spectroscopy and photometry of GD 448, a detached white dwarf - M dwarf binary with a period of 2.47h. We find that the NaI 8200A feature is composed of narrow emission lines due to irradiation of the M dwarf by the white dwarf within broad absorption lines that are essentially unaffected by heating. Combined with an improved spectroscopic orbit and gravitational red shift measurement from spectra of the H-alpha line, we are able to derive masses for the white dwarf and M dwarf directly (0.41 +/- 0.01 solar masses and 0.096 +/- 0.004 solar masses, respectively). We use a simple model of the CaII emission lines to establish the radius of the M dwarf assuming the emission from its surface to be proportional to the incident flux per unit area from the white dwarf. The radius derived is 0.125 +/- 0.020 solar radii. The M dwarf appears to be a normal main-sequence star in terms of its mass and radius and is less than half the size of its Roche lobe. The thermal timescale of the M dwarf is much longer than the cooling age of the white dwarf so we conclude that the M dwarf was unaffected by the common-envelope phase. The anomalous width of the H-alpha emission from the M dwarf remains to be explained, but the strengh of the line may be due to X-ray heating of the M dwarf due to accretion onto the white dwarf from the M dwarf wind.

WD1704+481 is a visual binary in which both components are white dwarfs. We present spectra of the H-alpha line of both stars which show that one component (WD1704+481.2 = Sanduleak B = GR 577) is a close binary with two white dwarf components. Thus, WD1704+481 is the first known triple degenerate star. From radial velocity measurements of the close binary we find an orbital period of 0.1448d, a mass ratio, q=Mbright/Mfaint of q=0.70+-0.03 and a difference in the gravitational redshifts of 11.5+-2.3km/s. The masses of the close pair of white dwarfs predicted by the mass ratio and gravitational redshift difference combined with theoretical cooling curves are 0.39+-0.05 solar mass and 0.56+-0.07 solar masses. WD1704+481 is therefore also likely to be the first example of a double degenerate in which the less massive white dwarf is composed of helium and the other white dwarf is composed of carbon and oxygen.

We have discovered that the spectrum of the well-known dwarf nova EM Cyg is contaminated by light from a K2-5V star (in addition to the K-type mass donor star). The K2-5V star contributes approximately 16 per cent of the light from the system and if not taken into account has a considerable effect upon radial velocity measurements of the mass donor star. We obtain a new radial velocity amplitude for the mass donor star of K2 = 202 +/- 3 km/s, which compares with the value of K2 = 135 +/- 3 km/s obtained in Stover, Robinson & Nather's classic 1981 study of EM Cyg. The revised value of the amplitude combined with a measurement of rotational broadening of the mass donor vsini = 140 +/- 6 km/s, leads to a new mass ratio of q = M2/M1 = 0.88 +/- 0.05. This solves a long standing problem with EM Cyg because Stover et al.'s measurements indicated a mass ratio q > 1, a value which should have led to dynamically unstable mass transfer for the secondary mass deduced by Stover et al. The revised value of the mass ratio combined with the orbital inclination i = 67 +/- 2 degrees leads to masses of 0.99 +/- 0.12 Msun and 1.12 +/- 0.08 Msun for the mass donor and white dwarf respectively. The mass donor is evolved, since it has a later spectral type (K3) than its mass would imply. We discuss whether the K star could be physically associated with EM Cyg or not, and present the results of the spectroscopic study.

We have discovered that the white dwarf PG 2329+267 is magnetic, and assuming a centered dipole structure, has a dipole magnetic field strength of approximately 2.3MG. This makes it one of only approximately 4% of isolated white dwarfs with a detectable magnetic field. Linear Zeeman splitting as well as quadratic Zeeman shifts are evident in the hydrogen Balmer sequence and circular spectropolarimetry reveals 10% circular polarisation in the two displaced sigma components of Halpha. We suggest from comparison with spectra of white dwarfs of known mass that PG 2329+267 is more massive than typical isolated white dwarfs, in agreement with the hypothesis that magnetic white dwarfs evolve from magnetic chemically peculiar Ap and Bp type main sequence stars.

We have discovered that the detached double degenerate binary WD 0957-666 has an orbital period of 1.46 hours, rather than the 1.15 day orbital period reported earlier. This is the shortest period example of such a system yet discovered. We obtain a unique period, which fits both our and earlier data. At this period the emission of gravitational radiation will cause the binary to merge within approximately 2.0 x 10*8 years. This system represents a population of short orbital period binaries which will merge within a Hubble time, and so could account for type Ia supernovae, although due to the low mass of both stars (0.3 to 0.4 solar masses), it is unlikely to become a supernova itself. We have detected the companion star and have measured a mass ratio of q = 1.15. This is the third double degenerate for which q has been measured and all three have q close to 1, which is in conflict with the predicted mass ratio distribution which peaks at 0.7. This system is viewed close to edge on, and we estimate that the probability of this system undergoing eclipses is 15 %.