Search

The upcoming Legacy Survey of Space and Time (LSST) at the Vera Rubin Observatory is expected to detect a few million transients per night, which will generate a live alert stream during the entire 10 years of the survey. This will be distributed via community brokers whose task is to select subsets of the stream and direct them to scientific communities. Given the volume and complexity of data, machine learning (ML) algorithms will be paramount for this task. We present the infrastructure tests and classification methods developed within the {\sc Fink} broker in preparation for LSST. This work aims to provide detailed information regarding the underlying assumptions, and methods, behind each classifier, enabling users to make informed follow-up decisions from {\sc Fink} photometric classifications. Using simulated data from the Extended LSST Astronomical Time-series Classification Challenge (ELAsTiCC), we showcase the performance of binary and multi-class ML classifiers available in {\sc Fink}. These include tree-based classifiers coupled with tailored feature extraction strategies, as well as deep learning algorithms. We introduce the CBPF Alert Transient Search (CATS), a deep learning architecture specifically designed for this task. Results show that {\sc Fink} classifiers are able to handle the extra complexity which is expected from LSST data. CATS achieved $97\%$ accuracy on a multi-class classification while our best performing binary classifier achieve $99\%$ when classifying the Periodic class. ELAsTiCC was an important milestone in preparing {\sc Fink} infrastructure to deal with LSST-like data. Our results demonstrate that {\sc Fink} classifiers are well prepared for the arrival of the new stream; this experience also highlights that transitioning from current infrastructures to Rubin will require significant adaptation of currently available tools., Comment: 18 pages, 29 figures. Submitted to A&A

Large surveys provide numerous non-targeted observations of small bodies (SSOs). The upcoming LSST of the Rubin observatory will be the largest source of SSO photometry in the next decade. With non-coordinated epochs of observation, colors, and therefore taxonomy and composition, can only be computed by comparing absolute magnitudes obtained in each filter by solving the phase function (evolution of brightness of the small body against the solar phase angle). Current models in use in the community (HG, HG12* , HG1G2) however fail to reproduce the long-term photometry of many targets due to the change in aspect angle between apparitions. We aim at deriving a generic yet simple phase function model accounting for the variable geometry of the SSOs over multiple apparitions. We propose the sHG1G2 phase function model in which we introduce a term describing the brightness changes due to spin orientation and polar oblateness. We apply this new model to 13,245,908 observations of 122,675 SSOs. These observations were acquired in the g and r filters with the Zwicky Transient Facility. We retrieve them and implement the new sHG1G2 model in Fink, a broker of alerts designed for the LSST. The sHG1G2 model leads to smaller residuals than other phase function models, providing a better description of the photometry of asteroids. We determine the absolute magnitude H and phase function coefficients (G1, G2) in each filter, the spin orientation (RA_0,DEC_0), and the polar-to-equatorial oblateness R for 95,593 Solar System Objects (SSOs), which constitutes about a tenfold increase in the number of characterised objects compared to current census. The application of the sHG1G2 model on ZTF alert data using the FINK broker shows that the model is appropriate to extract physical properties of asteroids from multi-band and sparse photometry, such as the forthcoming LSST survey.

We present a campaign designed to train the GRANDMA network and its infrastructure to follow up on transient alerts and detect their early afterglows. In preparation for O4 II campaign, we focused on GRB alerts as they are expected to be an electromagnetic counterpart of gravitational-wave events. Our goal was to improve our response to the alerts and start prompt observations as soon as possible to better prepare the GRANDMA network for the fourth observational run of LIGO-Virgo-Kagra (which started at the end of May 2023), and future missions such as SM. To receive, manage and send out observational plans to our partner telescopes we set up dedicated infrastructure and a rota of follow-up adcates were organized to guarantee round-the-clock assistance to our telescope teams. To ensure a great number of observations, we focused on Swift GRBs whose localization errors were generally smaller than the GRANDMA telescopes' field of view. This allowed us to bypass the transient identification process and focus on the reaction time and efficiency of the network. During 'Ready for O4 II', 11 Swift/INTEGRAL GRB triggers were selected, nine fields had been observed, and three afterglows were detected (GRB 220403B, GRB 220427A, GRB 220514A), with 17 GRANDMA telescopes and 17 amateur astronomers from the citizen science project Kilonova-Catcher. Here we highlight the GRB 220427A analysis where our long-term follow-up of the host galaxy allowed us to obtain a photometric redshift of $z=0.82\pm0.09$, its lightcurve elution, fit the decay slope of the afterglows, and study the properties of the host galaxy.

GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of observatories and from observational partners. Adding complementary data from the literature, we then derive essential physical parameters associated with the ejecta and external properties (i.e. the geometry and environment) of the GRB and compare with other analyses of this event. We spectroscopically confirm the presence of an associated supernova, SN2023pel, and we derive a photospheric expansion velocity of v $\sim$ 17$\times10^3$ km s$^{-1}$. We analyze the photometric data first using empirical fits of the flux and then with full Bayesian Inference. We again strongly establish the presence of a supernova in the data, with a maximum (pseudo-)bolometric luminosity of $5.75 \times 10^{42}$ erg/s, at $15.76^{+0.81}_{-1.21}$ days (in the observer frame) after the trigger, with a half-max time width of 22.0 days. We compare these values with those of SN1998bw, SN2006aj, and SN2013dx. Our best-fit model favours a very low density environment ($\log_{10}({n_{\rm ISM}/{\rm cm}^{-3}}) = -2.38^{+1.45}_{-1.60}$) and small values for the jet's core angle $\theta_{\rm core} = 1.54^{+1.02}_{-0.81} \ \rm{deg}$ and viewing angle $\theta_{\rm obs} = 0.76^{+1.29}_{-0.76} \ \rm{deg}$. GRB 230812B is thus one of the best observed afterglows with a distinctive supernova bump.

With the advent of large-scale astronomical surveys such as the Zwicky Transient Facility (ZTF), the number of alerts generated by transient, variable and moving astronomical objects is growing rapidly, reaching millions per night. Concerning solar system minor planets, their identification requires linking the alerts of many observations over a potentially large time, leading to a very large combinatorial number. This work aims to identify new candidates for solar system objects from massive alert data streams produced by large-scale surveys, such as the ZTF and the Vera C. Rubin Observatory's Legacy Survey of Space and Time. Our analysis used the Fink alert broker capabilities to reduce the 111,275,131 processed alerts from ZTF between November 2019 and December 2022 to only 389,530 new solar system alert candidates over the same period. We then implemented a linking algorithm, Fink-FAT, to create real-time trajectory candidates from alert data and extract orbital parameters. The analysis was validated on ZTF alert packets linked to confirmed solar system objects from the Minor Planet Center database. Finally, the results were confronted against follow-up observations. Between November 2019 and December 2022, Fink-FAT extracted 327 new orbits from solar system object candidates at the time of the observations, over which 65 were still unreported in the MPC database as of March 2023. After two late follow-up observation campaigns of six orbit candidates, four were associated with known solar system minor planets, and two remain unknown. Fink-FAT is deployed in the Fink broker and successfully analyzes in real time the alert data from the ZTF survey by regularly extracting new candidates for solar system objects. Our scalability tests also show that Fink-FAT can handle the even larger volume of alert data that the Rubin Observatory will send., Comment: Published in A&A

GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRANDMA+HXMT data, augmented with data from the literature up to 60 days. We then model numerically, using a Bayesian approach, the GRANDMA and HXMT-LE afterglow observations, that we augment with Swift-XRT and additional optical/NIR observations reported in the literature. We find that the GRB afterglow, extinguished by a large dust column, is most likely behind a combination of a large Milky-Way dust column combined with moderate low-metallicity dust in the host galaxy. Using the GRANDMA+HXMT-LE+XRT dataset, we find that the simplest model, where the observed afterglow is produced by synchrotron radiation at the forward external shock during the deceleration of a top-hat relativistic jet by a uniform medium, fits the multi-wavelength observations only moderately well, with a tension between the observed temporal and spectral evolution. This tension is confirmed when using the extended dataset. We find that the consideration of a jet structure (Gaussian or power-law), the inclusion of synchrotron self-Compton emission, or the presence of an underlying supernova do not improve the predictions, showing that the modelling of GRB22109A will require going beyond the most standard GRB afterglow model. Placed in the global context of GRB optical afterglows, we find the afterglow of GRB 221009A is luminous but not extraordinarily so, highlighting that some aspects of this GRB do not deviate from the global known sample despite its extreme energetics and the peculiar afterglow evolution., Comment: Accepted to ApJL for the special issue, 37 pages, 23 pages main text, 6 tables, 13 figures

We describe the fast transient classification algorithm in the center of the kilonova (KN) science module currently implemented in the Fink broker and report classification results based on simulated catalogs and real data from the ZTF alert stream. We used noiseless, homogeneously sampled simulations to construct a basis of principal components (PCs). All light curves from a more realistic ZTF simulation were written as a linear combination of this basis. The corresponding coefficients were used as features in training a random forest classifier. The same method was applied to long (>30 days) and medium (<30 days) light curves. The latter aimed to simulate the data situation found within the ZTF alert stream. Classification based on long light curves achieved 73.87% precision and 82.19% recall. Medium baseline analysis resulted in 69.30% precision and 69.74% recall, thus confirming the robustness of precision results when limited to 30 days of observations. In both cases, dwarf flares and point Type Ia supernovae were the most frequent contaminants. The final trained model was integrated into the Fink broker and has been distributing fast transients, tagged as KN_candidates, to the astronomical community, especially through the GRANDMA collaboration. We showed that features specifically designed to grasp different light curve behaviors provide enough information to separate fast (KN-like) from slow (non-KN-like) evolving events. This module represents one crucial link in an intricate chain of infrastructure elements for multi-messenger astronomy which is currently being put in place by the Fink broker team in preparation for the arrival of data from the Vera Rubin Observatory Legacy Survey of Space and Time., Comment: 8 Pages, 12 Figures, submitted to Astronomy and Astrophysics

We present on the status of POLARBEAR-2 A (PB2-A) focal plane fabrication. The PB2-A is the first of three telescopes in the Simon Array (SA), which is an array of three cosmic microwave background (CMB) polarization sensitive telescopes located at the POLARBEAR (PB) site in Northern Chile. As the successor to the PB experiment, each telescope and receiver combination is named as PB2-A, PB2-B, and PB2-C. PB2-A and -B will have nearly identical receivers operating at 90 and 150 GHz while PB2-C will house a receiver operating at 220 and 270 GHz. Each receiver contains a focal plane consisting of seven close-hex packed lenslet coupled sinuous antenna transition edge sensor bolometer arrays. Each array contains 271 di-chroic optical pixels each of which have four TES bolometers for a total of 7588 detectors per receiver. We have produced a set of two types of candidate arrays for PB2-A. The first we call Version 11 (V11) and uses a silicon oxide (SiOx) for the transmission lines and cross-over process for orthogonal polarizations. The second we call Version 13 (V13) and uses silicon nitride (SiNx) for the transmission lines and cross-under process for orthogonal polarizations. We have produced enough of each type of array to fully populate the focal plane of the PB2-A receiver. The average wirebond yield for V11 and V13 arrays is 93.2% and 95.6% respectively. The V11 arrays had a superconducting transition temperature (Tc) of 452 +/- 15 mK, a normal resistance (Rn) of 1.25 +/- 0.20 Ohms, and saturations powers of 5.2 +/- 1.0 pW and 13 +/- 1.2 pW for the 90 and 150 GHz bands respectively. The V13 arrays had a superconducting transition temperature (Tc) of 456 +/-6 mK, a normal resistance (Rn) of 1.1 +/- 0.2 Ohms, and saturations powers of 10.8 +/- 1.8 pW and 22.9 +/- 2.6 pW for the 90 and 150 GHz bands respectively.

GRANDMA is a world-wide collaboration with the primary scientific goal of studying gravitational-wave sources, discovering their electromagnetic counterparts and characterizing their emission. GRANDMA involves astronomers, astrophysicists, gravitational-wave physicists, and theorists. GRANDMA is now a truly global network of telescopes, with (so far) 30 telescopes in both hemispheres. It incorporates a citizen science programme (Kilonova-Catcher) which constitutes an opportunity to spread the interest in time-domain astronomy. The telescope network is an heterogeneous set of already-existing observing facilities that operate coordinated as a single observatory. Within the network there are wide-field imagers that can observe large areas of the sky to search for optical counterparts, narrow-field instruments that do targeted searches within a predefined list of host-galaxy candidates, and larger telescopes that are devoted to characterization and follow-up of the identified counterparts. Here we present an overview of GRANDMA after the third observing run of the LIGO/VIRGO gravitational-wave observatories in $2019-2020$ and its ongoing preparation for the forthcoming fourth observational campaign (O4). Additionally, we review the potential of GRANDMA for the discovery and follow-up of other types of astronomical transients., Comment: Submitted to the Proceedings of the SPIE, Astronomical Telescopes and Instrumentation 2022

We present our follow-up observations with GRANDMA of transient sources revealed by the Zwicky Transient Facility (ZTF). Over a period of six months, all ZTF triggers were examined in real time by a dedicated science module implemented in the Fink broker, which will be used for the data processing of the Vera C. Rubin Observatory. In this article, we present three selection methods to identify kilonova candidates. Out of more than 35 million candidates, a hundred sources have passed our selection criteria. Six were then followed-up by GRANDMA (by both professional and amateur astronomers). The majority were finally classified either as asteroids or as supernovae events. We mobilized 37 telescopes, bringing together a large sample of images, taken under various conditions and quality. To complement the orphan kilonova candidates (those without associated gamma-ray bursts, which were all), we included three additional supernovae alerts to conduct further observations of during summer 2021. We demonstrate the importance of the amateur astronomer community that contributed images for scientific analyzes of new sources discovered in a magnitude range r'=17-19 mag. We based our rapid kilonova classification on the decay rate of the optical source that should exceed 0.3 mag/day. GRANDMA's follow-up determined the fading rate within 1.5+/-1.2 days post-discovery, without waiting for further observations from ZTF. No confirmed kilonovae were discovered during our observing campaign. This work will be continued in the coming months in the view of preparing for kilonova searches in the next gravitational-wave observing run O4., Comment: 22 pages, 9 figures

LiteBIRD is a JAXA-led Strategic Large-Class mission designed to search for the existence of the primordial gravitational waves produced during the inflationary phase of the Universe, through the measurements of their imprint onto the polarization of the cosmic microwave background (CMB). These measurements, requiring unprecedented sensitivity, will be performed over the full sky, at large angular scales, and over 15 frequency bands from 34GHz to 448GHz. The LiteBIRD instruments consist of three telescopes, namely the Low-, Medium- and High-Frequency Telescope (respectively LFT, MFT and HFT). We present in this paper an overview of the design of the Medium-Frequency Telescope (89-224GHz) and the High-Frequency Telescope (166-448GHz), the so-called MHFT, under European responsibility, which are two cryogenic refractive telescopes cooled down to 5K. They include a continuous rotating half-wave plate as the first optical element, two high-density polyethylene (HDPE) lenses and more than three thousand transition-edge sensor (TES) detectors cooled to 100mK. We provide an overview of the concept design and the remaining specific challenges that we have to face in order to achieve the scientific goals of LiteBIRD., Comment: SPIE Conference

LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. JAXA selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with its expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD plans to map the cosmic microwave background (CMB) polarization over the full sky with unprecedented precision. Its main scientific objective is to carry out a definitive search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with an insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. To this end, LiteBIRD will perform full-sky surveys for three years at the Sun-Earth Lagrangian point L2 for 15 frequency bands between 34 and 448 GHz with three telescopes, to achieve a total sensitivity of 2.16 micro K-arcmin with a typical angular resolution of 0.5 deg. at 100GHz. We provide an overview of the LiteBIRD project, including scientific objectives, mission requirements, top-level system requirements, operation concept, and expected scientific outcomes., Comment: 20 pages, 9 figures

LiteBIRD has been selected as JAXA's strategic large mission in the 2020s, to observe the cosmic microwave background (CMB) $B$-mode polarization over the full sky at large angular scales. The challenges of LiteBIRD are the wide field-of-view (FoV) and broadband capabilities of millimeter-wave polarization measurements, which are derived from the system requirements. The possible paths of stray light increase with a wider FoV and the far sidelobe knowledge of $-56$ dB is a challenging optical requirement. A crossed-Dragone configuration was chosen for the low frequency telescope (LFT : 34--161 GHz), one of LiteBIRD's onboard telescopes. It has a wide field-of-view ($18^\circ \times 9^\circ$) with an aperture of 400 mm in diameter, corresponding to an angular resolution of about 30 arcminutes around 100 GHz. The focal ratio f/3.0 and the crossing angle of the optical axes of 90$^\circ$ are chosen after an extensive study of the stray light. The primary and secondary reflectors have rectangular shapes with serrations to reduce the diffraction pattern from the edges of the mirrors. The reflectors and structure are made of aluminum to proportionally contract from warm down to the operating temperature at $5\,$K. A 1/4 scaled model of the LFT has been developed to validate the wide field-of-view design and to demonstrate the reduced far sidelobes. A polarization modulation unit (PMU), realized with a half-wave plate (HWP) is placed in front of the aperture stop, the entrance pupil of this system. A large focal plane with approximately 1000 AlMn TES detectors and frequency multiplexing SQUID amplifiers is cooled to 100 mK. The lens and sinuous antennas have broadband capability. Performance specifications of the LFT and an outline of the proposed verification plan are presented., Comment: 21 pages, 14 figures

Counting pairs of galaxies or stars according to their distance is at the core of real-space correlation analyzes performed in astrophysics and cosmology. Upcoming galaxy surveys (LSST, Euclid) will measure properties of billions of galaxies challenging our ability to perform such counting in a minute-scale time relevant for the usage of simulations. The problem is only limited by efficient access to the data, hence belongs to the big data category. We use the popular Apache Spark framework to address it and design an efficient high-throughput algorithm to deal with hundreds of millions to billions of input data. To optimize it, we revisit the question of nonhierarchical sphere pixelization based on cube symmetries and develop a new one dubbed the "Similar Radius Sphere Pixelization" (SARSPix) with very close to square pixels. It provides the most adapted indexing over the sphere for all distance-related computations. Using LSST-like fast simulations, we compute autocorrelation functions on tomographic bins containing between a hundred million to one billion data points. In each case we achieve the construction of a standard pair-distance histogram in about 2 minutes, using a simple algorithm that is shown to scale, over a moderate number of nodes (16 to 64). This illustrates the potential of this new techniques in the field of astronomy where data access is becoming the main bottleneck. They can be easily adapted to other use-cases as nearest-neighbors search, catalog cross-match or cluster finding. The software is publicly available from https://github.com/astrolabsoftware/SparkCorr., Comment: published version

Using only cosmic microwave background polarization data from the POLARBEAR experiment, we measure $B$-mode polarization delensing on subdegree scales at more than $5\sigma$ significance. We achieve a 14% $B$-mode power variance reduction, the highest to date for internal delensing, and improve this result to 2% by applying for the first time an iterative maximum a posteriori delensing method. Our analysis demonstrates the capability of internal delensing as a means of improving constraints on inflationary models, paving the way for the optimal analysis of next-generation primordial $B$-mode experiments., Comment: Matches version published in Physical Review Letters

Apache Spark is a Big Data framework for working on large distributed datasets. Although widely used in the industry, it remains rather limited in the academic community or often restricted to software engineers. The goal of this paper is to show with practical uses-cases that the technology is mature enough to be used without excessive programming skills by astronomers or cosmologists in order to perform standard analyses over large datasets, as those originating from future galaxy surveys. To demonstrate it, we start from a realistic simulation corresponding to 10 years of LSST data taking (6 billions of galaxies). Then, we design, optimize and benchmark a set of Spark python algorithms in order to perform standard operations as adding photometric redshift errors, measuring the selection function or computing power spectra over tomographic bins. Most of the commands execute on the full 110 GB dataset within tens of seconds and can therefore be performed interactively in order to design full-scale cosmological analyses. A jupyter notebook summarizing the analysis is available at https://github.com/astrolabsoftware/1807.03078.

We report an improved measurement of the cosmic microwave background (CMB) $B$-mode polarization power spectrum with the POLARBEAR experiment at 150 GHz. By adding new data collected during the second season of observations (2013-2014) to re-analyzed data from the first season (2012-2013), we have reduced twofold the band-power uncertainties. The band powers are reported over angular multipoles $500 \leq \ell \leq 2100$, where the dominant $B$-mode signal is expected to be due to the gravitational lensing of $E$-modes. We reject the null hypothesis of no $B$-mode polarization at a confidence of 3.1$\sigma$ including both statistical and systematic uncertainties. We test the consistency of the measured $B$-modes with the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) framework by fitting for a single lensing amplitude parameter $A_L$ relative to the Planck best-fit model prediction. We obtain $A_L = 0.60 ^{+0.26} _{-0.24} ({\rm stat}) ^{+0.00} _{-0.04}({\rm inst}) \pm 0.14 ({\rm foreground}) \pm 0.04 ({\rm multi})$, where $A_{L}=1$ is the fiducial $\Lambda$CDM value, and the details of the reported uncertainties are explained later in the manuscript., Comment: 16 pages, 10 figures. Minor changes to match the published version. For data and figures, see http://bolo.berkeley.edu/polarbear/data/polarbear_BB_2017/

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017., Comment: 9pages,8figures

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $\sigma(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV., Comment: Accepted to Journal of Low Temperature Physics LTD16 Special Issue, Low Temperature Detector 16 Conference Proceedings, 5 pages, 1 figure

In this work we employ the pure-pseudo formalism devised to minimise the effects of the leakage on the variance of power spectrum estimates and discuss the limits on the tensor-to-scalar ratio, $r$, that could be realistically set by current and forthcoming measurements of the $B$-mode angular power spectrum. We compare those with the results obtained using other approaches: na\"{i}ve mode-counting, minimum-variance quadratic estimators, and re-visit the question of optimizing the sky coverage of small-scale, suborbital experiments in order to maximize the statistical significance of the detection of $r$. We show that the optimized sky coverage is largely insensitive to the adopted approach at least for reasonably compact sky patches. We find, however, that the mode-counting overestimates the detection significance by a factor $\sim1.17$ as compared to the lossless maximum variance approach and by a factor $\sim1.25$ as compared to the lossy pure pseudo-spectrum estimator. In a second time, we consider more realistic experimental configurations. With a pure pseudospectrum reconstruction of $B$-modes and considering only statistical uncertainties, we find that a detection of $r\sim0.11$, $r\sim0.0051$ and $r\sim0.0026$ at 99$\%$ of confidence level is within the reach of current sub-orbital experiments, future arrays of ground-based telescopes and a satellite mission, respectively. This means that an array of telescopes could be sufficient to discriminate between large- and small-field models of inflation, even if the $E$-to-$B$ leakage is consistently included but accounted for in the analysis. However, a satellite mission will be required to distinguish between different small-field models depending on the number of e-folds., Comment: 20 pages, 10 figures. Minor changes to match the published version. Published in Physical Review D

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations., Comment: 20 pages, 16 figures

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array., Comment: Presented at SPIE Astronomical Telescopes and Instrumentation 2014: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VII. Published in Proceedings of SPIE Volume 9153

We report a measurement of the B-mode polarization power spectrum in the cosmic microwave background (CMB) using the POLARBEAR experiment in Chile. The faint B-mode polarization signature carries information about the Universe's entire history of gravitational structure formation, and the cosmic inflation that may have occurred in the very early Universe. Our measurement covers the angular multipole range 500 < l < 2100 and is based on observations of an effective sky area of 25 square degrees with 3.5 arcmin resolution at 150 GHz. On these angular scales, gravitational lensing of the CMB by intervening structure in the Universe is expected to be the dominant source of B-mode polarization. Including both systematic and statistical uncertainties, the hypothesis of no B-mode polarization power from gravitational lensing is rejected at 97.1% confidence. The band powers are consistent with the standard cosmological model. Fitting a single lensing amplitude parameter A_BB to the measured band powers, A_BB = 1.12 +/- 0.61 (stat) +0.04/-0.12 (sys) +/- 0.07 (multi), where A_BB = 1 is the fiducial WMAP-9 LCDM value. In this expression, "stat" refers to the statistical uncertainty, "sys" to the systematic uncertainty associated with possible biases from the instrument and astrophysical foregrounds, and "multi" to the calibration uncertainties that have a multiplicative effect on the measured amplitude A_BB., Comment: 22 pages, 12 figures. v3 is updated to reflect the erratum published in ApJ 2017 848:73, which changed the rejection of the hypothesis of no B-mode polarization power from gravitational lensing from 97.2% to 97.1%

Gravitational lensing due to the large-scale distribution of matter in the cosmos distorts the primordial Cosmic Microwave Background (CMB) and thereby induces new, small-scale $B$-mode polarization. This signal carries detailed information about the distribution of all the gravitating matter between the observer and CMB last scattering surface. We report the first direct evidence for polarization lensing based on purely CMB information, from using the four-point correlations of even- and odd-parity $E$- and $B$-mode polarization mapped over $\sim30$ square degrees of the sky measured by the POLARBEAR experiment. These data were analyzed using a blind analysis framework and checked for spurious systematic contamination using null tests and simulations. Evidence for the signal of polarization lensing and lensing $B$-modes is found at 4.2$\sigma$ (stat.+sys.) significance. The amplitude of matter fluctuations is measured with a precision of $27\%$, and is found to be consistent with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model. This measurement demonstrates a new technique, capable of mapping all gravitating matter in the Universe, sensitive to the sum of neutrino masses, and essential for cleaning the lensing $B$-mode signal in searches for primordial gravitational waves., Comment: 7 pages, 3 figures. Matches Physical Review Letters accepted version. Submitted on 11 March 2014; accepted on 25 April 2014. The companion paper (arXiv:1312.6645) describes a measurement of polarization lensing in cross-correlation

We reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$\sigma$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$\sigma$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics., Comment: 6 pages, 2 figures. v2: replaced with version accepted for publication in Physical Review Letters. The companion paper (arXiv:1312.6646) describes a measurement of the polarization lensing power spectrum