Magnetic confinement of stellar winds leads to the formation of magnetospheres, which can be sculpted into Centrifugal Magnetospheres (CMs) by rotational support of the corotating plasma. The conditions required for the CMs of magnetic early B-type stars to yield detectable emission in H$\alpha$ -- the principal diagnostic of these structures -- are poorly constrained. A key reason is that no detailed study of the magnetic and rotational evolution of this population has yet been performed. Using newly determined rotational periods, modern magnetic measurements, and atmospheric parameters determined via spectroscopic modelling, we have derived fundamental parameters, dipolar oblique rotator models, and magnetospheric parameters for 56 early B-type stars. Comparison to magnetic A- and O-type stars shows that the range of surface magnetic field strength is essentially constant with stellar mass, but that the unsigned surface magnetic flux increases with mass. Both the surface magnetic dipole strength and the total magnetic flux decrease with stellar age, with the rate of flux decay apparently increasing with stellar mass. We find tentative evidence that multipolar magnetic fields may decay more rapidly than dipoles. Rotational periods increase with stellar age, as expected for a magnetic braking scenario. Without exception, all stars with H$\alpha$ emission originating in a CM are 1) rapid rotators, 2) strongly magnetic, and 3) young, with the latter property consistent with the observation that magnetic fields and rotation both decrease over time.

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

Atmospheric parameters determined via spectral modelling are unavailable for many of the known magnetic early B-type stars. We utilized high-resolution spectra together with NLTE models to measure effective temperatures $T_{\rm eff}$ and surface gravities $\log{g}$ of stars for which these measurements are not yet available. We find good agreement between our $T_{\rm eff}$ measurements and previous results obtained both photometrically and spectroscopically. For $\log{g}$, our results are compatible with previous spectroscopic measurements; however, surface gravities of stars previously determined photometrically have been substantially revised. We furthermore find that $\log{g}$ measurements obtained with HARPSpol are typically about 0.1 dex lower than those from comparable instruments. Luminosities were determined using Gaia Data Release 2 parallaxes. We find Gaia parallaxes to be unreliable for bright stars ($V<6$ mag) and for binaries; in these cases we reverted to Hipparcos parallaxes. In general we find luminosities systematically lower than those previously reported. Comparison of $\log{g}$ and $\log{L}$ to available rotational and magnetic measurements shows no correlation between either parameter with magnetic data, but a clear slow-down in rotation with both decreasing $\log{g}$ and increasing $\log{L}$, a result compatible with the expectation that magnetic braking should lead to rapid magnetic spindown that accelerates with increasing mass-loss.

NU Ori is a massive spectroscopic and visual binary in the Orion Nebula Cluster, with 4 components: Aa, Ab, B, and C. The B0.5 primary (Aa) is one of the most massive B-type stars reported to host a magnetic field. We report the detection of a spectroscopic contribution from the C component in high-resolution ESPaDOnS spectra, which is also detected in a Very Large Telescope Interferometer (VLTI) dataset. Radial velocity (RV) measurements of the inner binary (designated Aab) yield an orbital period of 14.3027(7) d. The orbit of the third component (designated C) was constrained using both RVs and interferometry. We find C to be on a mildly eccentric 476(1) d orbit. Thanks to spectral disentangling of mean line profiles obtained via least-squares deconvolution we show that the Zeeman Stokes $V$ signature is clearly associated with C, rather than Aa as previously assumed. The physical parameters of the stars were constrained using both orbital and evolutionary models, yielding $M_{\rm Aa} = 14.9 \pm 0.5 M_\odot$, $M_{\rm Ab} = 3.9 \pm 0.7 M_\odot$, and $M_{\rm C} = 7.8 \pm 0.7 M_\odot$. The rotational period obtained from longitudinal magnetic field $\langle B_z \rangle$ measurements is $P_{\rm rot} = 1.09468(7)$ d, consistent with previous results. Modeling of $\langle B_z \rangle$ indicates a surface dipole magnetic field strength of $\sim 8$ kG. NU Ori C has a magnetic field strength, rotational velocity, and luminosity similar to many other stars exhibiting magnetospheric H$\alpha$ emission, and we find marginal evidence of emission at the expected level ($\sim$1% of the continuum).

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

The rotational and magnetic properties of many magnetic hot stars are poorly characterized, therefore the MiMeS and BinaMIcS collaborations have collected extensive high-dispersion spectropolarimetric datasets of these targets. We present longitudinal magnetic field measurements $\langle B_{\rm z}\rangle$ for 52 early B-type stars (B5 to B0), with which we attempt to determine their rotational periods $P_{\rm rot}$. Supplemented with high-resolution spectroscopy, low-resolution DAO circular spectropolarimetry, and archival Hipparcos photometry, we determined $P_{\rm rot}$ for 10 stars, leaving only 5 stars for which $P_{\rm rot}$ could not be determined. Rotational ephemerides for 14 stars were refined via comparison of new to historical magnetic measurements. The distribution of $P_{\rm rot}$ is very similar to that observed for the cooler Ap/Bp stars. We also measured $v\sin{i}$ and $v_{\rm mac}$ for all stars. Comparison to non-magnetic stars shows that $v\sin{i}$ is much lower for magnetic stars, an expected consequence of magnetic braking. We also find evidence that $v_{\rm mac}$ is lower for magnetic stars. LSD profiles extracted using single-element masks revealed widespread, systematic discrepancies in $\langle B_{\rm z}\rangle$ between different elements: this effect is apparent only for chemically peculiar stars, suggesting it is a consequence of chemical spots. Sinusoidal fits to H line $\langle B_{\rm z}\rangle$ measurements (which should be minimally affected by chemical spots), yielded evidence of surface magnetic fields more complex than simple dipoles in 6 stars for which this has not previously been reported; however, in all 6 cases the second- and third-order amplitudes are small relative to the first-order (dipolar) amplitudes.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

We present our analysis of HD~35502 based on high- and medium-resolution spectropolarimetric observations. Our results indicate that the magnetic B5IVsnp star is the primary component of a spectroscopic triple system and that it has an effective temperature of $18.4\pm0.6\,{\rm kK}$, a mass of $5.7\pm0.6\,M_\odot$, and a polar radius of $3.0^{+1.1}_{-0.5}\,R_\odot$. The two secondary components are found to be essentially identical A-type stars for which we derive effective temperatures ($8.9\pm0.3\,{\rm kK}$), masses ($2.1\pm0.2\,M_\odot$), and radii ($2.1\pm0.4\,R_\odot$). We infer a hierarchical orbital configuration for the system in which the secondary components form a tight binary with an orbital period of $5.66866(6)\,{\rm d}$ that orbits the primary component with a period of over $40\,{\rm yrs}$. Least-Squares Deconvolution (LSD) profiles reveal Zeeman signatures in Stokes $V$ indicative of a longitudinal magnetic field produced by the B star ranging from approximately $-4$ to $0\,{\rm kG}$ with a median uncertainty of $0.4\,{\rm kG}$. These measurements, along with the line variability produced by strong emission in H$\alpha$, are used to derive a rotational period of $0.853807(3)\,{\rm d}$. We find that the measured $v\sin{i}=75\pm5\,{\rm km\,s}^{-1}$ of the B star then implies an inclination angle of the star's rotation axis to the line of sight of $24^{+6}_{-10}\degree$. Assuming the Oblique Rotator Model, we derive the magnetic field strength of the B star's dipolar component ($14^{+9}_{-3}\,{\rm kG}$) and its obliquity ($63\pm13\degree$). Furthermore, we demonstrate that the calculated Alfv\'{e}n radius ($41^{+17}_{-6}\,R_\ast$) and Kepler radius ($2.1^{+0.4}_{-0.7}\,R_\ast$) place HD~35502's central B star well within the regime of centrifugal magnetosphere-hosting stars.

We present results from an analysis of all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season. This includes the first Keck Array observations at 95 GHz. The maps reach a depth of 50 nK deg in Stokes $Q$ and $U$ in the 150 GHz band and 127 nK deg in the 95 GHz band. We take auto- and cross-spectra between these maps and publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353 GHz. An excess over lensed-LCDM is detected at modest significance in the 95x150 $BB$ spectrum, and is consistent with the dust contribution expected from our previous work. No significant evidence for synchrotron emission is found in spectra such as 23x95, or for correlation between the dust and synchrotron sky patterns in spectra such as 23x353. We take the likelihood of all the spectra for a multi-component model including lensed-LCDM, dust, synchrotron and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio $r$), using priors on the frequency spectral behaviors of dust and synchrotron emission from previous analyses of WMAP and Planck data in other regions of the sky. This analysis yields an upper limit $r_{0.05}<0.09$ at 95% confidence, which is robust to variations explored in analysis and priors. Combining these $B$-mode results with the (more model-dependent) constraints from Planck analysis of CMB temperature plus BAO and other data, yields a combined limit $r_{0.05}<0.07$ at 95% confidence. These are the strongest constraints to date on inflationary gravitational waves.