The Simons Observatory (SO) is an upcoming cosmic microwave background (CMB) experiment located on Cerro Toco, Chile, that will map the microwave sky in temperature and polarization in six frequency bands spanning 27 to 285 GHz. SO will consist of one 6-meter Large Aperture Telescope (LAT) fielding $\sim$30,000 detectors and an array of three 0.42-meter Small Aperture Telescopes (SATs) fielding an additional 30,000 detectors. This synergy will allow for the extremely sensitive characterization of the CMB over angular scales ranging from an arcmin to tens of degrees, enabling a wide range of scientific output. Here we focus on the SATs targeting degree angular scales with successive dichroic instruments observing at Mid-Frequency (MF: 93/145 GHz), Ultra-High-Frequency (UHF: 225/285 GHz), and Low-Frequency (LF: 27/39 GHz). The three SATs will be able to map $\sim$10% of the sky to a noise level of 2 $\mu$K-arcmin when combining 93 and 145 GHz. The multiple frequency bands will allow the CMB to be separated from galactic foregrounds (primarily synchrotron and dust), with the primary science goal of characterizing the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma \left(r\right) \approx 0.003$.

We report the development of the achromatic half-wave plate (AHWP) at millimeter wave for cosmic microwave background polarization experiments. We fabricate an AHWP consisting of nine a-cut sapphire plates based on the Pancharatnam recipe to cover a wide frequency range. The modulation efficiency and the phase are measured in a frequency range of 33 to 260 GHz with incident angles up to 10 degrees. We find the measurements at room temperature are in good agreement with the predictions. This is the broadest demonstration of the AHWP at the millimeter wave.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

LiteBIRD is a proposed JAXA satellite mission to measure the CMB B-mode polarization with unprecedented sensitivity ($\sigma_r\sim 0.001$). To achieve this goal, $4676$ state-of-the-art TES bolometers will observe the whole sky for 3 years from L2. These detectors, as well as the SQUID readout, are extremely susceptible to EMI and other instrumental disturbances e.g. static magnetic field and vibration. As a result, careful analysis of the interference between the detector system and the rest of the telescope instruments is essential. This study is particularly important during the early phase of the project, in order to address potential problems before the final assembly of the whole instrument. We report our plan for the preparation of a cryogenic testbed to study the interaction between the detectors and other subsystems, especially a polarization modulator unit consisting of a magnetically-rotating half wave plate. We also present the requirements, current status and preliminary results.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

We present visible and near-infrared observations of a near-Earth object (NEO), 2012 $\mathrm{TC_4}$. The NEO 2012 $\mathrm{TC_4}$ approached close to the Earth at a distance of about 50,000 km in October 2017. This close approach provided a practical exercise for planetary defense. This apparition was also an appropriate opportunity to investigate 2012 $\mathrm{TC_4}$, which is a monolithic asteroid \citep{Polishook13}. We conducted the observation campaign of 2012 $\mathrm{TC_4}$ using six small- and medium-sized telescopes. The multiband photometry analysis showed that the taxonomic class of 2012 $\mathrm{TC_4}$ to be an X-type. In particular, we successfully obtained the high time resolution lightcurve of 2012 $\mathrm{TC_4}$ with the Tomo-e Gozen camera, which is the world's first wide-field CMOS camera, mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The shape and rotational motion models of 2012 $\mathrm{TC_4}$ were derived from the lightcurve. When 2012 $\mathrm{TC_4}$ was assumed to be a triaxial ellipsoid, the rotational and precession periods were 8.47 $\pm$ 0.01 min and 12.25 $\pm$ 0.01 min, respectively, with the long axis mode. This indicates that 2012 $\mathrm{TC_4}$ is a tumbling and monolithic asteroid. The shape models showed that the plausible axial lengths to be 6.2 $\times$ 8.0 $\times$ 14.9~m or 3.3 $\times$ 8.0 $\times$ 14.3~m. The flattened and elongated shape indicates that 2012 $\mathrm{TC_4}$ is a fragment produced by a impact event. We also estimated the excitation timescale, which implied that the impact event happened within $\sim$3 $\times$ 10$^{5}$ yr and 2012 $\mathrm{TC_4}$ has a fresh surface.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

We study the stellar dynamics of the first star clusters after intermediate-mass black holes (IMBHs) are formed via runaway stellar collisions. We use the outputs of cosmological simulations of Sakurai et al. (2017) to follow the star cluster evolution in a live dark matter (DM) halo. Mass segregation within a cluster promotes massive stars to be captured by the central IMBH occasionally, causing tidal disruption events (TDEs). We find that the TDE rate scales with the IMBH mass as $\dot{N}_{\rm TDE}\sim0.3\,{\rm Myr}^{-1}(M_{\rm IMBH}/1000\,{\rm M}_{\odot})^2$. The DM component affects the star cluster evolution by stripping stars from the outer part. When the DM density within the cluster increases, the velocity dispersion of the stars increases, and then the TDE rate decreases. By the TDEs, the central IMBHs grow to as massive as $700-2500\,{\rm M}_{\odot}$ in 15 million years. The IMBHs are possible seeds for the formation of supermassive BHs observed at $z\gtrsim 6-7$, if a large amount of gas is supplied through galaxy mergers and/or large-scale gas accretion, or they might remain as IMBHs from the early epochs to the present-day Universe.

We describe the design of a cryogenic rotation stage (CRS) for use with the cryogenic half-wave plate (CHWP) polarization modulator on the POLARBEAR-2b and POLARBEAR-2c (PB2b/c) cosmic microwave background (CMB) experiments, the second and third installments of the Simons Array. Rapid modulation of the CMB polarization signal using a CHWP suppresses 1/f contamination due to atmospheric turbulence and allows a single polarimeter to measure both polarization states, mitigating systematic effects that arise when differencing orthogonal detectors. To modulate the full detector array while avoiding excess photon loading due to thermal emission, the CHWP must have a clear-aperture diameter of > 450 mm and be cooled to < 100 K. We have designed a 454-mm-clear-aperture, < 65 K CRS using a superconducting magnetic bearing driven by a synchronous magnetic motor. We present the specifications for the CRS, its interfacing to the PB2b/c receiver cryostat, its performance in a stand-alone test, and plans for future work.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

We perform a set of cosmological simulations of early structure formation with incorporating baryonic streaming motions. We present a case where a significantly elongated gas cloud with $\sim\!10^4\,$solar masses (${\rm M_\odot}$) is formed in a pre-galactic ($\sim\!10^7\,{\rm M_\odot}$) dark halo. The gas streaming into the halo compresses and heats the massive filamentary cloud to a temperature of $\sim\!10,000\,$Kelvin. The gas cloud cools rapidly by atomic hydrogen cooling, and then by molecular hydrogen cooling down to $\sim\!400\,$Kelvin. The rapid decrease of the temperature and hence of the Jeans mass triggers fragmentation of the filament to yield multiple gas clumps with a few hundred solar masses. We estimate the mass of the primordial star formed in each fragment by adopting an analytic model based on a large set of radiation hydrodynamics simulations of protostellar evolution. The resulting stellar masses are in the range of $\sim\!50$-$120\,{\rm M_\odot}$. The massive stars gravitationally attract each other and form a compact star cluster. We follow the dynamics of the star cluster using a hybrid $N$-body simulation. We show that massive star binaries are formed in a few million years through multi-body interactions at the cluster center. The eventual formation of the remnant black holes will leave a massive black hole binary, which can be a progenitor of strong gravitational wave sources similar to those recently detected by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO).

The driving mechanism of massive outflows observed in high-mass star-forming regions is investigated using three-dimensional magnetohydrodynamics (MHD) and protostellar evolution calculations. In our previous paper, we showed that the mass outflow rate depends strongly on the mass accretion rate onto the circumstellar disk around a high-mass protostar, and massive outflows may be driven by the magnetic effect in high-mass star-forming cores. In the present study, in order to verify that the MHD disk wind is the primary driving mechanism of massive outflows, we quantitatively compare outflow properties obtained through simulations and observations. Since the outflows obtained through simulations are slightly younger than those obtained through observations, the time-integrated quantities of outflow mass, momentum, and kinetic energy are slightly smaller than those obtained through observations. On the other hand, time-derivative quantities of mass ejection rate, outflow momentum flux, and kinetic luminosity obtained through simulations are in very good agreement with those obtained through observations. This indicates that the MHD disk wind greatly contributes to the massive outflow driving from high-mass protostars, and the magnetic field might significantly control the high-mass star formation process.

We study the formation of massive black holes in the first star clusters. We first locate star-forming gas clouds in proto-galactic haloes of $\gtrsim \!10^7\,{\rm M}_{\odot}$ in cosmological hydrodynamics simulations and use them to generate the initial conditions for star clusters with masses of $\sim \!10^5\,{\rm M}_{\odot}$. We then perform a series of direct-tree hybrid $N$-body simulations to follow runaway stellar collisions in the dense star clusters. In all the cluster models except one, runaway collisions occur within a few million years, and the mass of the central, most massive star reaches $\sim \!400-1900\,{\rm M}_{\odot}$. Such very massive stars collapse to leave intermediate-mass black holes (IMBHs). The diversity of the final masses may be attributed to the differences in a few basic properties of the host haloes such as mass, central gas velocity dispersion, and mean gas density of the central core. Finally, we derive the IMBH mass to cluster mass ratios, and compare them with the observed black hole to bulge mass ratios in the present-day Universe.

The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by Hosokawa et al. The mass accretion rates were taken from numerical hydrodynamics models of Vorobyov & Basu computing the circumstellar disk evolution starting from the gravitational collapse of pre-stellar cloud cores of various mass and angular momentum. The resulting stellar evolution tracks were compared with the isochrones and isomasses calculated using non-accreting models. We find that mass accretion in the initial several Myr of protostellar evolution can have a strong effect on the subsequent evolution of young stars and brown dwarfs. The disagreement between accreting and non-accreting models in terms of the total stellar luminosity L_st, stellar radius R_st and effective temperature T_eff depends on the thermal efficiency of accretion, i.e., on the fraction of accretion energy absorbed by the central object. The largest mismatch is found for the cold accretion case, in which essentially all accretion energy is radiated away. The relative deviations in L_st and R_st in this case can reach 50% for 1.0-Myr-old objects and remain notable even for 10-Myr-old objects. In the hot and hybrid accretion cases, in which a constant fraction of accretion energy is absorbed, the disagreement between accreting and non-accreting models becomes less pronounced, but still remains notable for 1.0-Myr-old objects. These disagreements may lead to the wrong age estimate for objects of (sub-)solar mass when using the isochrones based on non-accreting models (abridged).

The relation between the mass accretion rate onto the circumstellar disc and the rate of mass ejection by magnetically driven winds is investigated using three-dimensional magnetohydrodynamics simulations. Using a spherical cloud core with a varying ratio of thermal to gravitational energy, which determines the mass accretion rate onto the disc, to define the initial conditions, the outflow propagation for approximately 10^4 yr after protostar formation is then calculated for several cloud cores. The mass ejection rate and accretion rate are comparable only when the magnetic energy of the initial cloud core is comparable to the gravitational energy. Consequently, in strongly magnetised clouds a higher mass accretion rate naturally produces both massive protostars and massive outflows. The simulated outflow mass, momentum, kinetic energy and momentum flux agree well with observations, indicating that massive stars form through the same mechanism as low-mass stars but require a significantly strong magnetic field to launch massive outflows.

We study formation and long-term evolution of a circumstellar disk in a collapsing molecular cloud core using a resistive magnetohydrodynamic simulation. While the formed circumstellar disk is initially small, it grows as accretion continues and its radius becomes as large as 200 AUs toward the end of the Class-I phase. A pair of grand-design spiral arms form due to gravitational instability in the disk, and they transfer angular momentum in the highly resistive disk. Although the spiral arms disappear in a few rotations as expected in a classical theory, new spiral arms form recurrently as the disk soon becomes unstable again by gas accretion. Such recurrent spiral arms persist throughout the Class-0 and I phase. We then perform synthetic observations and compare our model with a recent high-resolution observation of a young stellar object Elias 2-27, whose circumstellar disk has grand design spiral arms. We find good agreement between our theoretical model and the observation. Our model suggests that the grand design spiral arms around Elias 2-27 are consistent with material arms formed by gravitational instability. If such spiral arms commonly exist in young circumstellar disks, it implies that young circumstellar disks are considerably massive and gravitational instability is the key process of angular momentum transport.

We perform one-dimensional radiation hydrodynamical simulations to solve accretion flows onto massive black holes (BHs) with a very high rate. Assuming that photon trapping limits the luminosity emerging from the central region to $L\lesssim L_{\rm Edd}$, Inayoshi, Haiman & Ostriker (2016) have shown that an accretion flow settles to a "hyper-Eddington" solution, with a steady and isothermal ($T\simeq 8000$ K) Bondi profile reaching $\gtrsim 5000$ times the Eddington accretion rate $\dot{M}_{\rm Edd}\equiv L_{\rm Edd}/c^2$. Here we address the possibility that gas accreting with finite angular momentum forms a bright nuclear accretion disc, with a luminosity exceeding the Eddington limit ($1\lesssim L/L_{\rm Edd}\lesssim 100$). Combining our simulations with an analytic model, we find that a transition to steady hyper-Eddington accretion still occurs, as long as the luminosity remains below $L/L_{\rm Edd}\lesssim 35~(M_{\rm BH}/10^4~M_\odot)^{3/2}(n_\infty/10^5~{\rm cm^{-3}})(T_\infty/10^4~{\rm K})^{-3/2}(r_{\star}/10^{14}~{\rm cm})^{-1/2}$, where $n_\infty$ and $T_\infty$ are the density and temperature of the ambient gas, and $r_\star$ is the radius of the photosphere, at which radiation emerges. If the luminosity exceeds this value, accretion becomes episodic. Our results can be accurately recovered in a toy model of an optically thick spherical shell, driven by radiation force into a collapsing medium. When the central source is dimmer than the above critical value, the expansion of the shell is halted and reversed by ram pressure of the collapsing medium, and by shell's weight. Our results imply that rapid, unimpeded hyper-Eddington accretion is possible even if the luminosity of the central source far exceeds the Eddington limit, and can be either steady or strongly episodic.

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a self-gravitating circumstellar disk. We show that the stellar evolution with such episodic accretion differs qualitatively from that expected with a constant accretion rate, even if the mean accretion rate is the same. Unlike the case of constant mass accretion, whereby the star expands roughly following $R_* \simeq 2.6 \times 10^3 R_\odot (M_*/100~M_\odot)^{1/2}$, the protostar can substantially contract during the quiescent phases between accretion bursts. The stellar effective temperature and ionizing photon emissivity increase accordingly as the star contracts, which can cause strong ionizing feedback and halt the mass accretion onto the star. With a fixed duration of the quiescent phase $\Delta t_{\rm q}$, such contraction occurs in early evolutionary phases, i.e. for $M_* \lesssim 10^3~M_\odot$ with $\Delta t_{\rm q} \simeq 10^3$ yr. For later epochs and larger masses but the same $\Delta t_{\rm q}$, contraction is negligible even during quiescent phases. With larger quiescent times $\Delta t_{\rm q}$, however, the star continues to contract during quiescent phases even for the higher stellar masses. We show that such behavior is well understood by comparing the interval time and the thermal relaxation time for a bloated surface layer. We conclude that the UV radiative feedback becomes effective if the quiescent phase associated by the burst accretion is longer than $\sim 10^3$ yr, which is possible in an accretion disk forming in the direct collapse model.