The Advanced Molybdenum-based Rare process Experiment (AMoRE) searches for neutrino-less double-beta (0{\nu}\b{eta}\b{eta}) decay of 100Mo in enriched molybdate crystals. The AMoRE crystals must have low levels of radioactive contamination to achieve low background signals with energies near the Q-value of the 100Mo 0{\nu}\b{eta}\b{eta} decay. To produce low-activity crystals, radioactive contaminants in the raw materials used to form the crystals must be controlled and quantified. 100EnrMoO3 powder, which is enriched in the 100Mo isotope, is of particular interest as it is the source of 100Mo in the crystals. A high-purity germanium detector having 100% relative efficiency, named CC1, is being operated in the Yangyang underground laboratory. Using CC1, we collected a gamma spectrum from a 1.6-kg 100EnrMoO3 powder sample enriched to 96.4% in 100Mo. Activities were analyzed for the isotopes 228Ac, 228Th, 226Ra, and 40K. They are long-lived naturally occurring isotopes that can produce background signals in the region of interest for AMoRE. Activities of both 228Ac and 228Th were < 1.0 mBq/kg at 90% confidence level (C.L.). The activity of 226Ra was measured to be 5.1 \pm 0.4 (stat) \pm 2.2 (syst) mBq/kg. The 40K activity was found as < 16.4 mBq/kg at 90% C.L.

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

We present $Spitzer$ 4.5$\mu$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s).

The Center for Underground Physics (CUP) of the Institute for Basic Science (IBS) is searching for the neutrinoless double-beta decay (0{\nu}\b{eta}\b{eta}) of 100Mo in the molybdate crystals of the AMoRE experiment. The experiment requires pure scintillation crystals to minimize internal backgrounds that can affect the 0{\nu}\b{eta}\b{eta} signal. For the last few years, we have been growing and studying Li2MoO4 crystals in a clean-environment facility to minimize external contamination during the crystal growth. Before growing Li2100MoO4 crystal, we have studied Li2natMoO4 crystal growth by a conventional Czochralski (CZ) grower. We grew a few different kinds of Li2natMO4 crystals using different raw materials in a campaign to minimize impurities. We prepared the fused Al2O3 refractories for the growth of ingots. Purities of the grown crystals were measured with high purity germanium detectors and by inductively coupled plasma mass spectrometry. The results show that the Li2MoO4 crystal has purity levels suitable for rare-event experiments. In this study, we present the growth of Li2MoO4 crystals at CUP and their purities.

This paper describes the liquid scintillator veto system for the COSINE-100 dark matter experiment and its performance. The COSINE-100 detector consists of 8 NaI(Tl) crystals immersed in 2,200~L of linear alkylbenzene based liquid scintillator. The liquid scintillator tags between 65 and 75\% of the internal $^{40}$K background in the 2--6~keV energy region. We also describe the background model for the liquid scintillator, which is primarily used to assess its energy calibration and threshold.

Protostars and young stars are strongly spatially "clustered" or "correlated" within their natal giant molecular clouds (GMCs). We demonstrate that such clustering leads to the conclusion that the incident bolometric radiative flux upon a random young star/disc is enhanced (relative to volume-averaged fluxes) by a factor which increases with the total stellar mass of the complex. Because the Galactic cloud mass function is top-heavy, the typical star in our Galaxy experienced a much stronger radiative environment than those forming in well-observed nearby (but relatively small) clouds, exceeding fluxes in the Orion Nebular Cluster by factors of $\gtrsim$30. Heating of the circumstellar disc around a median young star is dominated by this external radiation beyond $\sim 50\,$AU. And if discs are not well-shielded by ambient dust, external UV irradiation can dominate over the host star down to sub-AU scales. Another consequence of stellar clustering is an extremely broad Galaxy-wide distribution of incident flux (spanning $>10$ decades), with half the Galactic star formation in a substantial "tail" towards even more intense background radiation. We also show that the strength of external irradiation is amplified super-linearly in high-density environments such as the Galactic centre, starbursts, or high-redshift galaxies.

The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (tSZ) effect. The third measure determines the lowest order relativistic correction to the tSZ signal, which is seeing increased observational relevance. Our procedure allows us to accurately model the relativistic SZ (rSZ) contribution and we show that a $\gtrsim 10\%-40\%$ underestimation of this rSZ cluster temperature is expected when applying standard X-ray relations. The correction also exhibits significant mass and redshift evolution, as we demonstrate here. We present the mass dependence of each temperature measure alongside their profiles and a short analysis of the temperature dispersion as derived from the aforementioned simulations. We also discuss a new relation connecting the temperature and Compton-$y$ parameter, which can be directly used for rSZ modelling. Simple fits to the obtained scaling relations and profiles are provided. These should be useful for future studies of the rSZ effect and its relevance to cluster cosmology.

Super-Earths and mini-Neptunes exhibit great diversity in their compositional and orbital properties. Their bulk densities span a large range, from those dense enough to be purely rocky to those needing a substantial contribution from volatiles to their volumes. Their orbital configurations range from compact, circular multi-transiting systems like Kepler-11 to systems like our Solar System's terrestrial planets with wider spacings and modest but significant eccentricities and mutual inclinations. Here we investigate whether a continuum of formation conditions resulting from variation in the amount of solids available in the inner disk can account for the diversity of orbital and compositional properties observed for super Earths, including the apparent dichotomy between single transiting and multiple transiting system. We simulate in situ formation of super-Earths via giant impacts and compare to the observed Kepler sample. We find that intrinsic variations among disks in the amount of solids available for in situ formation can account for the orbital and compositional diversity observed among Kepler's transiting planets. Our simulations can account for the planets' distributions of orbital period ratios, transit duration ratios, and transit multiplicity; higher eccentricities for single than multi transiting planets; smaller eccentricities for larger planets; scatter in the mass-radius relation, including lower densities for planets with masses measured with TTVs than RVs; and similarity in planets' sizes and spacings within each system. Our findings support the theory that variation among super-Earth and mini-Neptune properties is primarily locked in by different in situ formation conditions, rather than arising stochastically through subsequent evolution.

We present ground-based infrared transit observations for four dynamically interacting \textit{Kepler} planets, including Kepler-29b, Kepler-36c, KOI-1783.01, and Kepler-177c, obtained using the Wide-field Infrared Camera on the Hale 200" telescope at Palomar Observatory. By utilizing an engineered diffuser and custom guiding software, we mitigate time-correlated telluric and instrumental noise sources in these observations. We achieve an infrared photometric precision comparable to or better than that of space-based observatories such as the \textit{Spitzer Space Telescope}, and detect transits with greater than 3$\sigma$ significance for all planets. For Kepler-177c ($J=13.9$) our measurement uncertainties are only $1.2\times$ the photon noise limit and 1.9 times better than the predicted photometric precision for \textit{Spitzer} IRAC photometry of this same target. We find that a single transit observation obtained $4-5$ years after the end of the original \textit{Kepler} mission can reduce dynamical mass uncertainties by as much as a factor of three for these systems. Additionally, we combine our new observations of KOI-1783.01 with information from the literature to confirm the planetary nature of this system. We discuss the implications of our new mass and radius constraints in the context of known exoplanets with low incident fluxes, and we note that Kepler-177c may be a more massive analog to the currently known super-puffs given its core mass (3.8$\pm0.9M_\Earth$) and large gas-to-core ratio (2.8$\pm0.7$). Our demonstrated infrared photometric performance opens up new avenues for ground-based observations of transiting exoplanets previously thought to be restricted to space-based investigation.

There is a long standing debate about whether or not the annual modulation signal reported by the DAMA/LIBRA collaboration is induced by Weakly Interacting Massive Particles~(WIMP) in the galaxy's dark matter halo scattering from nuclides in their NaI(Tl) crystal target/detector. This is because regions of WIMP-mass vs. WIMP-nucleon cross-section parameter space that can accommodate the DAMA/LIBRA-phase1 modulation signal in the context of the standard WIMP dark matter galactic halo and isospin-conserving~(canonical), spin-independent~(SI) WIMP-nucleon interactions have been excluded by many of other dark matter search experiments including COSINE-100, which uses the same NaI(Tl) target/detector material. Moreover, the recently released DAMA/LIBRA-phase2 results are inconsistent with an interpretation as WIMP-nuclide scattering via the canonical SI interaction and prefer, instead, isospin-violating or spin-dependent interactions. Dark matter interpretations of the DAMA/LIBRA signal are sensitive to the NaI(Tl) scintillation efficiency for nuclear recoils, which is characterized by so-called quenching factors~(QF), and the QF values used in previous studies differ significantly from recently reported measurements, which may have led to incorrect interpretations of the DAMA/LIBRA signal. In this article, the compatibility of the DAMA/LIBRA and COSINE-100 results, in light of the new QF measurements is examined for different possible types of WIMP-nucleon interactions. The resulting allowed parameter space regions associated with the DAMA/LIBRA signal are explicitly compared with 90\% confidence level upper limits from the initial 59.5~day COSINE-100 exposure. With the newly measured QF values, the allowed 3$\sigma$ regions from the DAMA/LIBRA data are still generally excluded by the COSINE-100 data.

Compton scattering between electrons and photons plays a crucial role in astrophysical plasmas. Many important aspects of this process can be captured by using the so-called Compton scattering kernel. For isotropic media, exact analytic expressions (valid at all electron and photon energies) do exist but are hampered by numerical issues and often are presented in complicated ways. In this paper, we summarize, simplify and improve existing analytic expressions for the Compton scattering kernel, with an eye on clarity and physical understanding. We provide a detailed overview of important properties of the kernel covering a wide range of energies and highlighting aspects that have not been appreciated as much previously. We discuss analytic expressions for the moments of the kernel, comparing various approximations and demonstrating their precision. We also illustrate the properties of the scattering kernel for thermal electrons at various temperatures and photon energies, introducing new analytic approximations valid to high temperatures. The obtained improved formulae for the kernel and its moments should prove useful in many astrophysical computations, one of them being the evolution of spectral distortions of the cosmic microwave background in the early Universe. A novel code, CSpack, for efficient computations of the Compton scattering kernel and its properties (in the future also including anisotropies in the initial electron and photon distributions) is being developed in a series of papers and will be available within one month.

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

COSINE-100 is a direct detection dark matter search experiment that uses a 106 kg array of eight NaI(Tl) crystals that are kept underground at the Yangyang Underground Laboratory to avoid cosmogenic activation of radioisotopes by cosmic rays. Even though the cosmogenic activity is declining with time, there are still significant background rates from the remnant nuclides. In this paper, we report measurements of cosmogenic isotope contaminations with less than one year half-lives that are based on extrapolations of the time dependent activities of their characteristic energy peaks to activity rates at the time the crystals were deployed underground. For longer-lived $^{109}$Cd ($T_{1/2}=1.6$ y) and $^{22}$Na ($T_{1/2}=2.6$ y), we investigate time correlations of characteristic $\gamma$/X-ray peaks. The inferred sea-level production rates are compared with caluclations based on the ACTIVIA and MENDL-2 model calculations and experimental data. For $^{3}$H, which has a long, 12.3 year half-life, we evaluated the activity levels from the exposure times and determined a cosmogenic activation rate that is consistent with other measurements.

Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $\delta$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted flux tube from the solar interior into a non-magnetized stratified atmosphere. Results. We emerge two $\Omega$-shaped loops forming a quadrupolar region. The emerging flux forms two initially separated bipoles, that later come in contact creating a $\delta$-spot central region. Above the two bipoles, two magnetic lobes expand and interact through a series of current sheets at the interface between them. Two recurrent confined eruptions are produced. In both cases, reconnection between sheared low-lying field lines forms a flux rope. Reconnection between the two lobes high in the atmosphere forms retracting field lines that push against the flux rope, creating a current sheet between them. It also forms a third magnetic lobe between the two emerged ones, that later acts as a strapping field. The flux rope eruptions are triggered when the reconnection between the flux ropes and the field above them becomes efficient enough to remove the tension of the overlying field. These reconnection events occur internally in the system. The first erupting flux rope almost fully reconnects with the overlying field. The second eruption is confined by the overlying field. During the confined eruption, the flux rope is enhanced in size, flux and twist, similar to confined-flare-to-flux-rope observations. Proxies of the emission reveal the two erupting filaments channels. A flare arcade is formed only in the second eruption due to the longer-lasting and more efficient reconnection below the flux rope.

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 a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and our new Y-band images from the Very Large Telescope SPHERE instrument. We use the knowledge of the disc's radial extent inferred from ALMA observations and the grain size distribution found by SED fitting to generate semi-dynamical dust models of the disc. We compare the models to scattered light and thermal emission data and find that a disc with a maximum of the surface density at 110~au and shallow edges can describe both thermal emission and scattered light observations. This suggests that grains close to the blow-out limit and large grains stem from the same planetesimal population and are mainly influenced by radiation pressure. The influence of inwards transport processes could not be analysed in this study.

Sub-Saturns straddle the boundary between gas-rich Jupiters and gas-poor super-Earths/sub-Neptunes. Their large radii (4--8$R_\oplus$) suggest that their gas-to-core mass ratios range $\sim$0.1--1.0. With their envelopes as massive as their cores, sub-Saturns are just on the verge of runaway gas accretion; they are expected to be significantly less populous than gas giants. Yet, the observed occurrence rates of sub-Saturns and Jupiters are comparable within $\sim$100 days. We show that in these inner regions of planetary systems, the growth of sub-Saturns/Jupiters is ultimately limited by local and global hydrodynamic flows---runaway accretion terminates and the formation of gas giants is suppressed. We derive a simple analytic formula for the local hydrodynamic accretion rate---an expression that has been previously reported only as an empirical fit to numerical simulations. Evolving simultaneously the background disk gas and the gas accretion onto planetary cores, we find that both the ubiquity of super-Earths/sub-Neptunes and the rarity of gas-rich planets are best explained if an underlying core mass distribution is peaked at $\sim$4.3$M_\oplus$. Within a finite disk lifetime $\sim$10 Myrs, massive cores ($\gtrsim 10M_\oplus$) can become either gas-poor or gas-rich depending on when they assemble but smaller cores ($\lesssim 10M_\oplus$) can only become gas-poor. This wider range of possible outcomes afforded by more massive cores may explain why metal-rich stars harbor a more diverse set of planets.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

We recommend an intensive effort to survey and understand the obliquity distribution of small close-in extrasolar planets over the coming decade. The orbital obliquities of exoplanets--i.e., the relative orientation between the planetary orbit and the stellar rotation--is a key tracer of how planets form and migrate. While the orbital obliquities of smaller planets are poorly explored today, a new generation of facilities coming online over the next decade will make such observations possible en masse. Transit spectroscopic observations with the extremely large telescopes will enable us to measure the orbital obliquities of planets as small as $\sim2R_{\oplus}$ around a wide variety of stars, opening a window into the orbital properties of the most common types of planets. This effort will directly contribute to understanding the formation and evolution of planetary systems, a key objective of the National Academy of Sciences' Exoplanet Science Strategies report.

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 use radial velocity observations to search for long-period gas giant companions in systems hosting inner super-Earth (1-4 R_Earth, 1-10 M_Earth) planets to constrain formation and migration scenarios for this population. We consistently re-fit published RV datasets for 65 stars and find 9 systems with statistically significant trends indicating the presence of an outer companion. We combine these RV data with AO images to constrain the masses and semi-major axes of these companions. We quantify our sensitivity to the presence of long-period companions by fitting the sample with a power law distribution and find an occurrence rate of 39+/-7% for companions 0.5-20 M_Jup and 1-20 AU. Half of our systems were discovered by the transit method and half were discovered by the RV method. While differences in RV baselines and number of data points between the two samples lead to different sensitivities to distant companions, we find that occurrence rates of gas giant companions in each sample are consistent at the 0.5$\sigma$ level. We compare the frequency of Jupiter analogs in these systems to the equivalent rate from field star surveys and find that Jupiter analogs are more common around stars hosting super-Earths. We conclude that the presence of outer gas giants does not suppress the formation of inner super-Earths, and that these two populations of planets instead appear to be correlated. We also find that the stellar metallicities of systems with gas giant companions are higher than those without companions, in agreement with the well-established metallicity correlation from RV surveys of field stars.

The plasma emission process, which is the mechanism for solar type II and type III radio bursts phenomena, is studied by means of particle-in-cell and weak turbulence simulation methods. By plasma emission, it is meant as a loose description of a series of processes, starting from the solar flare associated electron beam exciting Langmuir and ion-acoustic turbulence, and subsequent partial conversion of beam energy into the radiation energy by nonlinear processes. Particle-in-cell (PIC) simulation is rigorous but the method is computationally intense, and it is difficult to diagnose the results. Numerical solution of equations of weak turbulence (WT) theory, termed WT simulation, on the other hand, is efficient and naturally lends itself to diagnostics since various terms in the equation can be turned on or off. Nevertheless, WT theory is based upon a number of assumptions. It is, therefore, desirable to compare the two methods, which is carried out for the first time in the present paper with numerical solutions of the complete set of equations of the WT theory and with two-dimensional electromagnetic PIC simulation. Upon making quantitative comparisons it is found that WT theory is largely valid, although some discrepancies are also found. The present study also indicates that it requires large computational resources in order to accurately simulate the radiation emission processes, especially for low electron beam speeds. Findings from the present paper thus imply that both methods may be useful for the study of solar radio emissions as they are complementary.

We present the most sensitive direct imaging and radial velocity (RV) exploration of $\epsilon$ Eridani to date. $\epsilon$ Eridani is an adolescent planetary system, reminiscent of the early Solar system. It is surrounded by a prominent and complex debris disk which is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 years ago, but has met with scrutiny due to possible confusion with stellar noise. We confirm the planet with a new compilation and analysis of precise RV data spanning 30 years, and combine it with upper limits from our direct imaging search, the most sensitive ever performed. The deep images were taken in the Ms band (4.7$\mu$m) with the vortex coronagraph recently installed in W.M. Keck Observatory's infrared camera NIRC2, which opens a sensitive window for planet searches around nearby adolescent systems. The RV data and direct imaging upper limit maps were combined in an innovative joint Bayesian analysis, providing new constraints on the mass and orbital parameters of the elusive planet. $\epsilon$ Eridani b has a mass of $0.78^{+0.38}_{-0.12}$ $M_{Jup}$ and is orbiting $\epsilon$ Eridani at about $3.48\pm 0.02$ AU with a period of $7.37 \pm 0.07$ years. The eccentricity of $\epsilon$ Eridani b's orbit is $0.07^{+0.06}_{-0.05}$, an order of magnitude smaller than early estimates and consistent with a circular orbit. We discuss our findings from the standpoint of planet-disk interactions and prospects for future detection and characterization with the James Webb Space Telescope.

A solar active region (AR) that produces at least one M- or X-class major flare tends to produce multiple flares during its passage across the solar disk. It will be interesting if we can estimate how flare-productive a given major flaring AR is for a time interval of several days, by investigating time series of its photospheric magnetic field properties. For this, we studied 93 major flaring ARs observed from 2010 to 2016 by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). More specifically, for each AR under study, the mean and fluctuation were calculated from an 8-day time series of each of 18 photospheric magnetic parameters extracted from the Space-weather HMI Active Region Patch (SHARP) vector magnetogram products at 12-min cadence. We then compared these with the AR's 8-day flare index, which is defined as the sum of soft X-ray peak fluxes of flares produced in the AR during the same interval of the 8-day SHARP parameter time series. As a result, it is found that the 8-day flare index is well correlated with the mean and/or fluctuation values of some magnetic parameters (with correlation coefficients of 0.6-0.7 in log-log space). Interestingly, the 8-day flare index shows a slightly better correlation with the fluctuation than the mean for the SHARP parameters associated with the surface integral of photospheric magnetic free energy density. We also discuss how the correlation varies if the 8-day flare index is compared with the mean or fluctuation calculated from an initial portion of the SHARP parameter time series.

Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs (<<1%) are "real" but correspond to early stages, the true "per-freefall" SFE where most stars actually form can be much larger. Very high (>>10%) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20x larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artifact. We also explore individual dense "clumps" within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ~2x larger than observed, possibly indicating physics other than feedback from massive (main sequence) stars is needed to regulate their collapse.

Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with the largest accretion cross-sections also have the fastest radial drift speeds and are the most easily drained out of discs. We present a global model of planet formation by pebble accretion that keeps track of the disc's mass budget. Cores, each initialized with a lunar mass, grow from discs whose finite stores of mm-cm sized pebbles drift inward across all radii in viscously accreting gas. For every 1 $M_\oplus$ netted by a core, at least 10 $M_\oplus$ and possibly much more are lost to radial drift. Core growth rates are typically exponentially sensitive to particle Stokes number, turbulent Mach number, and solid surface density. This exponential sensitivity, when combined with disc migration, tends to generate binary outcomes from 0.1-30 AU: either sub-Earth cores remain sub-Earth, or explode into Jupiters, with the latter migrating inward to varying degrees. When Jupiter-breeding cores assemble from mm-cm sized pebbles, they do so in discs where such particles drain out in $\sim$10$^5$ yr or less; such fast-draining discs do not fit mm-wave observations.

Recent distance ladder determinations of the Hubble constant $H_0$ disagree at about the $3.5\sigma$ level with the value determined from Planck measurements of the cosmic microwave background (CMB) assuming a $\Lambda$CDM cosmology. This discrepancy has prompted speculation that new physics might be required beyond that assumed in the $\Lambda$CDM model. In this paper, we apply the inverse distance ladder to fit a parametric form of $H(z)$ to baryon acoustic oscillation (BAO) and Type Ia supernova data together with priors on the sound horizon at the end of the radiation drag epoch, $r_d$. We apply priors on $r_d$, based on inferences from either Planck or the Wilkinson Microwave Anistropy Probe (WMAP), and demonstrate that these values are consistent with CMB-independent determinations of $r_d$ derived from measurements of the primordial deuterium abundance, BAO and supernova data assuming the $\Lambda$CDM cosmology. The $H(z)$ constraints that we derive are independent of detailed physics within the dark sector at low redshifts, relying only on the validity of the Friedmann-Robertson-Walker (FRW) metric of General Relativity. For each assumed prior on $r_d$, we find consistency with the inferred value of $H_0$ and the Planck $\Lambda$CDM value and corresponding tension with the distance ladder estimate.

The majority of gas giants (planets of masses $\gtrsim10^2 M_\oplus$) are found to reside at distances beyond $\sim1$ au from their host stars. Within 1 au, the planetary population is dominated by super-Earths of $2-20 M_\oplus$. We show that this dichotomy between inner super-Earths and outer gas giants can be naturally explained should they form in nearly inviscid disks. In laminar disks, a planet can more easily repel disk gas away from its orbit. The feedback torque from the pile-up of gas inside the planet's orbit slows down and eventually halts migration. A pressure bump outside the planet's orbit traps pebbles and solids, starving the core. Gas giants are born cold and stay cold: more massive cores are preferentially formed at larger distances, and they barely migrate under disk feedback. We demonstrate this using 2D hydrodynamical simulations of disk-planet interaction lasting up to $10^5$ years: we track planet migration and pebble accretion until both come to an end by disk feedback. Whether cores undergo runaway gas accretion to become gas giants or not is determined by computing 1D gas accretion models. Our simulations show that in an inviscid minimum mass solar nebula, gas giants do not form inside $\sim$0.5 au, nor can they migrate there while the disk is present. We also explore the dependence on disk mass, and find that gas giants form further out in less massive disks.

Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such environments can be optically thin to both incident starlight and internally generated thermal radiation. We construct time-dependent models of such atmospheres, showing that super-Earths/sub-Neptunes can accrete their $\sim$1%-by-mass gas envelopes, and super-puffs/sub-Saturns their $\sim$20%-by-mass envelopes, over a wide range of nebular depletion histories requiring no fine tuning. Although nascent atmospheres can exhibit stratospheric temperature inversions effected by atomic Fe and various oxides that absorb strongly at visible wavelengths, the rate of gas accretion remains controlled by the radiative-convective boundary (rcb) at much greater pressures. For dusty envelopes, the temperature at the rcb $T_{\rm rcb} \simeq 2500$ K is still set by ${\rm H}_2$ dissociation; for dust-depleted envelopes, $T_{\rm rcb}$ tracks the temperature of the visible or thermal photosphere, whichever is deeper, out to at least $\sim$5 AU. The rate of envelope growth remains largely unchanged between the old radiative diffusion models and the new optically thin models, reinforcing how robustly super-Earths form as part of the endgame chapter in disk evolution.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We present results from the first measurement of axial range components of fiducialized neutron induced nuclear recoil tracks using the DRIFT directional dark matter detector. Nuclear recoil events are fiducialized in the DRIFT experiment using temporal charge carrier separations between different species of anions in 30:10:1 Torr of CS$_2$:CF$_4$:O$_2$ gas mixture. For this measurement, neutron-induced nuclear recoil tracks were generated by exposing the detector to $^{252}$Cf source from different directions. Using these events, the sensitivity of the detector to the expected axial directional signatures were investigated as the neutron source was moved from one detector axis to another. Results obtained from these measurements show clear sensitivity of the DRIFT detector to the axial directional signatures in this fiducialization gas mode.

Sub-Neptunes around FGKM dwarfs are evenly distributed in log orbital period down to $\sim$10 days, but dwindle in number at shorter periods. Both the break at $\sim$10 days and the slope of the occurrence rate down to $\sim$1 day can be attributed to the truncation of protoplanetary disks by their host star magnetospheres at co-rotation. We demonstrate this by deriving planet occurrence rate profiles from empirical distributions of pre-main-sequence stellar rotation periods. Observed profiles are better reproduced when planets are distributed randomly in disks---as might be expected if planets formed in situ---rather than piled up near disk edges, as would be the case if they migrated in by disk torques. Planets can be brought from disk edges to ultra-short ($< 1$ day) periods by asynchronous equilibrium tides raised on their stars. Tidal migration can account for how ultra-short period planets (USPs) are more widely spaced than their longer period counterparts. Our picture provides a starting point for understanding why metal-rich stars tend to harbor more USPs, and why the sub-Neptune population drops at $\sim$10 days regardless of whether the host star is of type FGK or early M. We predict planet occurrence rates around A stars to also break at short periods, but at $\sim$1 day instead of $\sim$10 days because A stars rotate faster than lower mass stars (this prediction presumes that the planetesimal building blocks of planets can drift inside the dust sublimation radius).

We present results from a 54.7 live-day shielded run of the DRIFT-IId detector, the world's most sensitive, directional, dark matter detector. Several improvements were made relative to our previous work including a lower threshold for detection, a more robust analysis and a tenfold improvement in our gamma rejection factor. After analysis, no events remain in our fiducial region leading to an exclusion curve for spin-dependent WIMP-proton interactions which reaches 0.28 pb at 100 GeV/c^2 a fourfold improvement on our previous work. We also present results from a 45.4 live-day unshielded run of the DRIFT-IId detector during which 14 nuclear recoil-like events were observed. We demonstrate that the observed nuclear recoil rate of 0.31+/-0.08 events per day is consistent with detection of ambient, fast neutrons emanating from the walls of the Boulby Underground Science Facility.

Measuring precise planet masses, densities, and orbital dynamics in individual planetary systems is an important pathway toward understanding planet formation. The WASP-47 system has an unusual architecture that motivates a complex formation theory. The system includes a hot Jupiter ("b") neighbored by interior ("e") and exterior ("d") sub-Neptunes, and a long-period eccentric giant planet ("c"). We simultaneously modeled transit times from the Kepler K2 Mission and 118 radial velocities to determine precise masses, densities, and Keplerian orbital elements of the WASP-47 planets. Combining RVs and TTVs provides a better estimate of the mass of planet d (13.6\pm2.0~M_\odot) than obtained with only RVs ($12.75\pm2.70~M_\odot$) or TTVs ($16.1\pm3.8~M_\odot$). Planets e and d have high densities for their size, consistent with a history of photo-evaporation and/or formation in a volatile-poor environment. Through our RV and TTV analysis, we find that the planetary orbits have eccentricities similar to the solar system planets. The WASP-47 system has three similarities to our own solar system: (1) the planetary orbits are nearly circular and coplanar, (2) the planets are not trapped in mean motion resonances, and (3) the planets have diverse compositions. None of the current single-process exoplanet formation theories adequately reproduce these three characteristics of the WASP-47 system (or our solar system). We propose that WASP-47, like the solar system, formed in two stages: first, the giant planets formed in a gas-rich disk and migrated to their present locations, and second, the high-density sub-Neptunes formed in situ in a gas-poor environment.

This study presents a catalog of 8107 molecular clouds that covers the entire Galactic plane and includes 98% of the $^{12}$CO emission observed within $b\pm5^\circ$. The catalog was produced using a hierarchical cluster identification method applied to the result of a Gaussian decomposition of the Dame et al. data. The total H$_2$ mass in the catalog is $1.2\times10^9 M_\odot$, in agreement with previous estimates. We find that 30% of the sight lines intersect only a single cloud, with another 25% intersecting only two clouds. The most probable cloud size is $R\sim30$ pc. We find that $M\propto R^{2.2\pm0.2}$, with no correlation between the cloud surface density, $\Sigma$, and $R$. In contrast with the general idea, we find a rather large range of values of $\Sigma$, from 2 to $300 M_\odot$ pc$^{-2}$, and a systematic decrease with increasing Galactic radius, $R_{\rm gal}$. The cloud velocity dispersion and the normalization $\sigma_0=\sigma_v/R^{1/2}$ both decrease systematically with $R_{\rm gal}$. When studied over the whole Galactic disk, there is a large dispersion in the line width-size relation, and a significantly better correlation between $\sigma_v$ and $\Sigma\,R$. The normalization of this correlation is constant to better than a factor of two for $R_{\rm gal}<20$ kpc. This relation is used to disentangle the ambiguity between near and far kinematic distances. We report a strong variation of the turbulent energy injection rate. In the outer Galaxy it may be maintained by accretion through the disk and/or onto the clouds, but neither source can drive the 100 times higher cloud-averaged injection rate in the inner Galaxy.

The measurement of the direction of WIMP-induced nuclear recoils is a compelling but technologically challenging strategy to provide an unambiguous signature of the detection of Galactic dark matter. Most directional detectors aim to reconstruct the dark-matter-induced nuclear recoil tracks, either in gas or solid targets. The main challenge with directional detection is the need for high spatial resolution over large volumes, which puts strong requirements on the readout technologies. In this paper we review the various detector readout technologies used by directional detectors. In particular, we summarize the challenges, advantages and drawbacks of each approach, and discuss future prospects for these technologies.

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

Star formation on galactic scales is known to be a slow process, but whether it is slow on smaller scales is uncertain. We cross-correlate 5469 giant molecular clouds (GMCs) from a new all-sky catalog with 256 star forming complexes (SFCs) to build a sample of 191 SFC-GMC complexes---collections of multiple clouds each matched to 191 SFCs. The total mass in stars harbored by these clouds is inferred from WMAP free-free fluxes. We measure the GMC mass, the virial parameter, the star formation efficiency $\epsilon$ and the star formation rate per free-fall time $\epsilon_{\rm ff}$. Both $\epsilon$ and $\epsilon_{\rm ff}$ range over 3--4 orders of magnitude. We find that 68.3% of the clouds fall within $\sigma_{\log\epsilon}=0.79\pm0.22\,{\rm dex}$ and $\sigma_{\log\epsilon_{\rm ff}}=0.91\pm0.22\,{\rm dex}$ about the median. Compared to these observed scatters, a simple model with a time independent $\epsilon_{\rm ff}$ that depends on the host GMC properties predicts $\sigma_{\log\epsilon_{\rm ff}}=$0.12-0.24. Allowing for a time-variable $\epsilon_{\rm ff}$, we can recover the large dispersion in the rate of star formation. This strongly suggests that star formation in the Milky Way is a dynamic process on GMC scales. We also show that the surface star formation rate profile of the Milky Way correlates well with the molecular gas surface density profile.

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks. These interpretations suggest that the farthest reaches of planetary systems are perturbed by eccentric planets, possibly just a few Earth masses each.

We report the discovery and confirmation of two sub-Saturn planets orbiting a bright (V = 11.3), metal-rich ([Fe/H] = 0.42 $\pm$ 0.04 dex) G3 dwarf in the K2 Campaign 2 field. The planets are 5.68 $\pm$ 0.56 Earth-radii and 7.82 $\pm$ 0.72 Earth-radii and have orbital periods of 20.8851 $\pm$ 0.0003 d and 42.3633$\pm$0.0006 d, near to the 2:1 mean-motion resonance. We obtained 32 radial velocities (RVs) with Keck/HIRES and detected the reflex motion due to EPIC-203771098b and c. These planets have masses of 21.0 $\pm$ 5.4 Earth-masses and 27.0 $\pm$ 6.9 Earth-masses, respectively. With low densities of 0.63 $\pm$ 0.25 g/cc and 0.31 $\pm$ 0.12 g/cc, respectively, the planets require thick envelopes of H/He to explain their large sizes and low masses. Interior structure models predict that the planets have fairly massive cores of 17.6 $\pm$ 4.3 Earth-masses and 16.1 $\pm$ 4.2 Earth-masses, respectively. They may have formed exterior to their present locations, accreted their H/He envelopes at large orbital distances, and migrated in as a resonant pair. The proximity to resonance, large transit depths, and host star brightness offer rich opportunities for TTV follow-up. Finally, the low surface gravities of the EPIC-203771098 planets make them favorable targets for transmission spectroscopy by HST, Spitzer, and JWST.

The giant impact phase of terrestrial planet formation establishes connections between super-Earths' orbital properties (semimajor axis spacings, eccentricities, mutual inclinations) and interior compositions (the presence or absence of gaseous envelopes). Using N-body simulations and analytic arguments, we show that spacings derive not only from eccentricities, but also from inclinations. Flatter systems attain tighter spacings, a consequence of an eccentricity equilibrium between gravitational scatterings, which increase eccentricities, and mergers, which damp them. Dynamical friction by residual disk gas plays a critical role in regulating mergers and in damping inclinations and eccentricities. Systems with moderate gas damping and high solid surface density spawn gas-enveloped super-Earths with tight spacings, small eccentricities, and small inclinations. Systems in which super-Earths coagulate without as much ambient gas, in disks with low solid surface density, produce rocky planets with wider spacings, larger eccentricities, and larger mutual inclinations. A combination of both populations can reproduce the observed distributions of spacings, period ratios, transiting planet multiplicities, and transit duration ratios exhibited by Kepler super-Earths. The two populations, both formed in situ, also help to explain observed trends of eccentricity vs. planet size, and bulk density vs. method of mass measurement (radial velocities vs. transit timing variations). Simplifications made in this study --- including the limited timespan of the simulations, and the approximate treatments of gas dynamical friction and gas depletion history --- should be improved upon in future work to enable a detailed quantitative comparison to the observations.

The riddle posed by super-Earths (1-4$R_\oplus$, 2-20$M_\oplus$) is that they are not Jupiters: their core masses are large enough to trigger runaway gas accretion, yet somehow super-Earths accreted atmospheres that weigh only a few percent of their total mass. We show that this puzzle is solved if super-Earths formed late, as the last vestiges of their parent gas disks were about to clear. This scenario would seem to present fine-tuning problems, but we show that there are none. Ambient gas densities can span many (up to 9) orders of magnitude, and super-Earths can still robustly emerge after $\sim$0.1-1 Myr with percent-by-weight atmospheres. Super-Earth cores are naturally bred in gas-poor environments where gas dynamical friction has weakened sufficiently to allow constituent protocores to merge. So little gas is present at the time of core assembly that cores hardly migrate by disk torques: formation of super-Earths can be in situ. The picture --- that close-in super-Earths form in a gas-poor (but not gas-empty) inner disk, fed continuously by gas that bleeds inward from a more massive outer disk --- recalls the largely evacuated but still accreting inner cavities of transitional protoplanetary disks. We also address the inverse problem presented by super-puffs: an uncommon class of short-period planets seemingly too voluminous for their small masses (4-10$R_\oplus$, 2-6$M_\oplus$). Super-puffs easily acquire their thick atmospheres as dust-free, rapidly cooling worlds outside $\sim$1 AU where nebular gas is colder, less dense, and therefore less opaque. Unlike super-Earths which can form in situ, super-puffs migrated in to their current orbits; they are expected to form the outer links of mean-motion resonant chains, and to exhibit greater water content. We close by confronting observations and itemizing remaining questions.

Planets acquire atmospheres from their parent circumstellar disks. We derive a general analytic expression for how the atmospheric mass grows with time $t$, as a function of the underlying core mass $M_{\rm core}$ and nebular conditions, including the gas metallicity $Z$. Planets accrete as much gas as can cool: an atmosphere's doubling time is given by its Kelvin-Helmholtz time. Dusty atmospheres behave differently from atmospheres made dust-free by grain growth and sedimentation. The gas-to-core mass ratio (GCR) of a dusty atmosphere scales as GCR $\propto t^{0.4} M_{\rm core}^{1.7} Z^{-0.4} \mu_{\rm rcb}^{3.4}$, where $\mu_{\rm rcb} \propto 1/(1-Z)$ (for $Z$ not too close to 1) is the mean molecular weight at the innermost radiative-convective boundary. This scaling applies across all orbital distances and nebular conditions for dusty atmospheres; their radiative-convective boundaries, which regulate cooling, are not set by the external environment, but rather by the internal microphysics of dust sublimation, H$_2$ dissociation, and the formation of H$^-$. By contrast, dust-free atmospheres have their radiative boundaries at temperatures $T_{\rm rcb}$ close to nebular temperatures $T_{\rm out}$, and grow faster at larger orbital distances where cooler temperatures, and by extension lower opacities, prevail. At 0.1 AU in a gas-poor nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.9} M_{\rm core}^{1.6} Z^{-0.4} \mu_{\rm rcb}^{3.3}$, while beyond 1 AU in a gas-rich nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.5} M_{\rm core}^1 Z^{-0.4}\mu_{\rm rcb}^{2.2}$. We confirm our analytic scalings against detailed numerical models for objects ranging in mass from Mars (0.1 $M_\oplus$) to the most extreme super-Earths (10-20 $M_\oplus$), and explain why heating from planetesimal accretion cannot prevent the latter from undergoing runaway gas accretion.

Planets with sizes between those of Earth and Neptune divide into two populations: purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger gas-enveloped planets possessing voluminous and optically thick atmospheres. We show that whether a planet forms rocky or gas-enveloped depends on the solid surface density of its parent disk. Assembly times for rocky cores are sensitive to disk solid surface density. Lower surface densities spawn smaller planetary embryos; to assemble a core of given mass, smaller embryos require more mergers between bodies farther apart and therefore exponentially longer formation times. Gas accretion simulations yield a rule of thumb that a rocky core must be at least 2$M_\oplus$ before it can acquire a volumetrically significant atmosphere from its parent nebula. In disks of low solid surface density, cores of such mass appear only after the gas disk has dissipated, and so remain purely rocky. Higher surface density disks breed massive cores more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We test model predictions against observations, using planet radius as an observational proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface density. Theory can explain the observation that metal-rich stars host predominantly gas-enveloped planets.

The addition of O2 to gas mixtures in time projection chambers containing CS2 has recently been shown to produce multiple negative ions that travel at slightly different velocities. This allows a measurement of the absolute position of ionising events in the z (drift) direction. In this work, we apply the z-fiducialisation technique to a directional dark matter search. In particular, we present results from a 46.3 live-day source-free exposure of the DRIFT-IId detector run in this completely new mode. With full-volume fiducialisation, we have achieved the first background-free operation of a directional detector. The resulting exclusion curve for spin-dependent WIMP-proton interactions reaches 1.1 pb at 100 GeV/c2, a factor of 2 better than our previous work. We describe the automated analysis used here, and argue that detector upgrades, implemented after the acquisition of these data, will bring an additional factor of >3 improvement in the near future.

We present numerical evidence of dynamic star formation in which the accreted stellar mass grows superlinearly with time, roughly as $t^2$. We perform simulations of star formation in self-gravitating hydrodynamic and magneto-hydrodynamic turbulence that is continuously driven. By turning the self-gravity of the gas in the simulations on or off, we demonstrate that self-gravity is the dominant physical effect setting the mass accretion rate at early times before feedback effects take over, contrary to theories of turbulence-regulated star formation. We find that gravitational collapse steepens the density profile around stars, generating the power-law tail on what is otherwise a lognormal density probability distribution function. Furthermore, we find turbulent velocity profiles to flatten inside collapsing regions, altering the size-linewidth relation. This local flattening reflects enhancements of turbulent velocity on small scales, as verified by changes to the velocity power spectra. Our results indicate that gas self-gravity dynamically alters both density and velocity structures in clouds, giving rise to a time-varying star formation rate. We find that a substantial fraction of the gas that forms stars arrives via low density flows, as opposed to accreting through high density filaments.

We report on numerical experiment of the recurrent onset of helical "blowout" jets in an emerging flux region (EFR). We find that these jets are running with velocities of $\sim 100-250\,\mathrm{kms^{-1}}$ and they transfer a vast amount of heavy plasma into the outer solar atmosphere. During their emission, they undergo an untwisting motion as a result of reconnection between the twisted emerging and the non-twisted pre-existing magnetic field in the solar atmosphere. For the first time in the context of blowout jets, we provide a direct evidence that their untwisting motion is associated with the propagation of torsional Alfv\'en waves in the corona.

Close-in super-Earths having radii 1--4 $R_\oplus$ may possess hydrogen atmospheres comprising a few percent by mass of their rocky cores. We determine the conditions under which such atmospheres can be accreted by cores from their parent circumstellar disks. Accretion from the nebula is problematic because it is too efficient: we find that 10-$M_\oplus$ cores embedded in solar metallicity disks tend to undergo runaway gas accretion and explode into Jupiters, irrespective of orbital location. The threat of runaway is especially dire at $\sim$0.1 AU, where solids may coagulate on timescales orders of magnitude shorter than gas clearing times; thus nascent atmospheres on close-in orbits are unlikely to be supported against collapse by planetesimal accretion. The time to runaway accretion is well approximated by the cooling time of the atmosphere's innermost convective zone, whose extent is controlled by where H$_2$ dissociates. Insofar as the temperatures characterizing H$_2$ dissociation are universal, timescales for core instability tend not to vary with orbital distance --- and to be alarmingly short for 10-$M_\oplus$ cores. Nevertheless, in the thicket of parameter space, we identify two scenarios, not mutually exclusive, that can reproduce the preponderance of percent-by-mass atmospheres for super-Earths at $\sim$0.1 AU, while still ensuring the formation of Jupiters at $\gtrsim 1$ AU. Scenario (a): planets form in disks with dust-to-gas ratios that range from $\sim$20$\times$ solar at 0.1 AU to $\sim$2$\times$ solar at 5 AU. Scenario (b): the final assembly of super-Earth cores from mergers of proto-cores --- a process that completes quickly at $\sim$0.1 AU once begun --- is delayed by gas dynamical friction until just before disk gas dissipates completely.

Radon gas emanating from materials is of interest in environmental science and also a major concern in rare event non-accelerator particle physics experiments such as dark matter and double beta decay searches, where it is a major source of background. Notable for dark matter experiments is the production of radon progeny recoils (RPRs), the low energy (~100 keV) recoils of radon daughter isotopes, which can mimic the signal expected from WIMP interactions. Presented here are results of measurements of radon emanation from detector materials in the 1 metre cubed DRIFT-II directional dark matter gas time projection chamber experiment. Construction and operation of a radon emanation facility for this work is described, along with an analysis to continuously monitor DRIFT data for the presence of internal 222Rn and 218Po. Applying this analysis to historical DRIFT data, we show how systematic substitution of detector materials for alternatives, selected by this device for low radon emanation, has resulted in a factor of ~10 reduction in internal radon rates. Levels are found to be consistent with the sum from separate radon emanation measurements of the internal materials and also with direct measurement using an attached alpha spectrometer. The current DRIFT detector, DRIFT-IId, is found to have sensitivity to 222Rn of 2.5 {\mu}Bq/l with current analysis efficiency, potentially opening up DRIFT technology as a new tool for sensitive radon assay of materials.

Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial conditions and reconnection micro-physics is investigated with a parameter study. In particular, we aim to understand the stability of the coronal flux rope in the context of X-point collapse and the effects of boundary driving in both unstratified and stratified atmospheres. In the absence of driving, we assess the behavior of waves in the vicinity of the X-point. With boundary driving applied, we study the effects of reconnection micro-physics and atmospheric stratification on the eruption. We find that the Chen-Shibata equilibrium can be unstable to an X-point collapse even in the absence of driving due to wave accumulation at the X-point. However, the equilibrium can be stabilized by reducing the compressibility of the plasma, which allows small-amplitude waves to pass through the X-point without accumulation. Simulations with the photospheric boundary driving evaluate the impact of reconnection micro-physics and atmospheric stratification on the resulting dynamics: we show the evolution of the system to be determined primarily by the structure of the global magnetic fields with little sensitivity to the micro-physics of magnetic reconnection; and in a stratified atmosphere, we identify a novel mechanism for producing quasi-periodic behavior at the reconnection site behind a rising flux rope as a possible explanation of similar phenomena observed in solar and stellar flares.

We study the surface brightness profiles of a sample of brightest cluster galaxies (BCGs) with 0.3<z<0.9. The BCGs are selected from the first Red-sequence Cluster Survey and an X-ray cluster survey. The surface brightness profiles of the BCGs are measured using HST ACS images, and most of them can be well modeled by a Sersic profile with index ~ 6 and half-light radius ~ 30 kpc. Although the single Sersic model fits the profiles well, we argue that the systematics in the sky background measurement and the coupling between the model parameters make the comparison of the model parameters ambiguous. Direct comparison of the BCG profiles, on the other hand, has revealed an inside-out growth for these most massive galaxies: as the mass of a BCG increases, the central mass density of the galaxy increases slowly ($\rho_{1kpc} \propto M*^{0.39}$), while the slope of the outer profile grows continuously shallower ($\alpha \propto M*^{-2.5}$). Such a fashion of growth continues down to the less massive early-type galaxies (ETGs), without apparent distinction between BCGs and non-BCGs. For the very massive ETGs and BCGs, the slope of the Kormendy relation starts to trace the slope of the surface brightness profiles and becomes insensitive to subtle profile evolution. These results are generally consistent with dry mergers being the major driver of the mass growth for BCGs and massive ETGs. We also find strong correlations between the richness of clusters and the properties of BCGs: the more massive the clusters are, the more massive the BCGs ($M*_{bcg} \propto M_{clusters}^{0.6}$) and the shallower their surface brightness profiles. After taking into account the bias in the cluster samples, we find the masses of the BCGs have grown by at least a factor of 1.5 from z=0.5 to z=0, in contrast to the previous findings of no evolution. Such an evolution validates the expectation from the LCDM model.

The DRIFT-IId dark matter detector is a m$^3$-scale low-pressure TPC with directional sensitivity to WIMP-induced nuclear recoils. Its primary backgrounds were due to alpha decays from contamination on the central cathode. Efforts to reduce these backgrounds led to replacing the 20 \mu m wire central cathode with one constructed from 0.9 \mu m aluminized mylar, which is almost totally transparent to alpha particles. Detailed modeling of the nature and origin of the remaining backgrounds led to an in-situ, ppt-sensitive assay of alpha decay backgrounds from the central cathode. This led to further improvements in the thin-film cathode resulting in over 2 orders of magnitude reduction in backgrounds compared to the wire cathode. Finally, the addition of O$_2$ to CS$_2$ gas was found to produce multiple species of electronegative charge carriers, providing a method to determine the absolute position of nuclear recoils and reject all known remaining backgrounds while retaining a high efficiency for nuclear recoil detection.

Using spectroscopic data, we support the possibility of suprathermal distributions of coronal ions by fitting the equivalent kappa functions to their emission line profiles. We fit different kappa and Gaussian model functions to line profiles of the strong Fe xv line at 284.16 Angstroms, across two large-field spectroscopic rasters taken in a solar active region. Both single- and double-component Gaussian models are applied, as well as two kappa models, one with a free width parameter allowing for and the other with a constrained width that precludes "microturbulence". We then compare the goodness of fit of the computed best fits for each model. The kappa distribution is a generalization, or superset, of the Maxwellian, so they are able to fit line profiles more precisely than a Gaussian. In most of the data, the best-fit kappa model produces much lower residuals across the profile than any single Gaussian and sometimes double Gaussian. Most importantly, the distribution of estimated kappa values is found to lie mostly in the low-kappa range, implying ion populations far from Maxwellian. Even when the width is removed as a free parameter of fit, the kappa model is still able to fit the data credibly, again with low best-fit values of kappa. We find the shape of the Fe xv line, in the vast majority of the data analyzed, to be indicative of a highly suprathermal ion population.

The conversion of gas into stars is a fundamental process in astrophysics and cosmology. Stars are known to form from the gravitational collapse of dense clumps in interstellar molecular clouds, and it has been proposed that the resulting star formation rate is proportional to either the amount of mass above a threshold gas surface density, or the gas volume density. These star-formation prescriptions appear to hold in nearby molecular clouds in our Milky Way Galaxy's disk as well as in distant galaxies where the star formation rates are often much larger. The inner 500 pc of our Galaxy, the Central Molecular Zone (CMZ), contains the largest concentration of dense, high-surface density molecular gas in the Milky Way, providing an environment where the validity of star-formation prescriptions can be tested. Here we show that by several measures, the current star formation rate in the CMZ is an order-of-magnitude lower than the rates predicted by the currently accepted prescriptions. In particular, the region 1 deg < l < 3.5 deg, |b| < 0.5 deg contains ~10^7 Msun of dense molecular gas -- enough to form 1000 Orion-like clusters -- but the present-day star formation rate within this gas is only equivalent to that in Orion. In addition to density, another property of molecular clouds, such as the amplitude of turbulent motions, must be included in the star-formation prescription to predict the star formation rate in a given mass of molecular gas.

We analyze Spitzer GLIMPSE, MSX, and WMAP images of the Milky Way to identify 8 micron and free-free sources in the Galaxy. Seventy-two of the eighty-eight WMAP sources have coverage in the GLIMPSE and MSX surveys suitable for identifying massive star forming complexes (SFC). We measure the ionizing luminosity functions of the SFCs and study their role in the turbulent motion of the Galaxy's molecular gas. We find a total Galactic free-free flux f_{\nu} = 46177.6 Jy; the 72 WMAP sources with full 8 micron coverage account for 34263.5 Jy (~75%), with both measurements made at \nu=94GHz (W band). We find a total of 280 SFCs, of which 168 have unique kinematic distances and free-free luminosities. We use a simple model for the radial distribution of star formation to estimate the free-free and ionizing luminosity for the sources lacking distance determinations. The total dust-corrected ionizing luminosity is Q = 2.9 \pm 0.5 x 10^53 photons s^-1, which implies a galactic star formation rate of 1.2 \pm 0.2 M_{\sun} yr^-1. We present the (ionizing) luminosity function of the SFCs, and show that 24 sources emit half the ionizing luminosity of the Galaxy. The SFCs appear as bubbles in GLIMPSE or MSX images; the radial velocities associated with the bubble walls allow us to infer the expansion velocity of the bubbles. We calculate the kinetic luminosity of the bubble expansion and compare it to the turbulent luminosity of the inner molecular disk. SFCs emitting 80% of the total galactic free-free luminosity produce a kinetic luminosity equal to 65% of the turbulent luminosity in the inner molecular disk. This suggests that the expansion of the bubbles is a major driver of the turbulent motion of the inner Milky Way molecular gas.

We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the red and blue shifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor 5 or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.

The current status of the DRIFT (Directional Recoil Identification From Tracks) experiment at Boulby Mine is presented, including the latest limits on the WIMP spin-dependent cross-section from 1.5 kg days of running with a mixture of CS2 and CF4. Planned upgrades to DRIFT IId are detailed, along with ongoing work towards DRIFT III, which aims to be the world's first 10 m3-scale directional Dark Matter detector.

Data are presented from the DRIFT-IId detector housed in the Boulby mine in northeast England. A 0.8 m^3 fiducial volume, containing partial pressures of 30 Torr CS2 and 10 Torr CF4, was exposed for a duration of 47.4 live-time days with sufficient passive shielding to provide a neutron free environment within the detector. The nuclear recoil events seen are consistent with a remaining low level background from the decay of progeny of radon daughters attached to the central cathode of the detector. However, energy depositions from such events must drift across the entire width of the detector, and thus display large diffusion upon reaching the readout planes of the device. Exploiting this feature, it is shown to be possible to reject energy depositions from these radon decay progeny events while still retaining sensitivity to nuclear recoil events. The response of the detector is then interpreted, using the F nuclei content of the gas, in terms of sensitivity to proton spin-dependent WIMP-nucleon interactions, displaying a minimum in sensitivity cross section at 0.5 pb for a WIMP mass of 100 GeV/c^2.

SONYC - Substellar Objects in Nearby Young Clusters - is a survey program to investigate the frequency and properties of brown dwarfs down to masses below the Deuterium burning limit in nearby star forming regions. In this second paper, we present results on the ~1 Myr old cluster Rho Ophiuchi, combining our own deep optical and near-infrared imaging using Subaru with photometry from the 2-Micron All Sky Survey and the Spitzer Space Telescope. Of the candidates selected from iJKs photometry, we have confirmed three -- including a new brown dwarf with a mass close to the Deuterium limit -- as likely cluster members through low-resolution infrared spectroscopy. We also identify 27 sub-stellar candidates with mid-infrared excess consistent with disk emission, of which 16 are new and 11 are previously spectroscopically confirmed brown dwarfs. The high and variable extinction makes it difficult to obtain the complete sub-stellar population in this region. However, current data suggest that its ratio of low-mass stars to brown dwarfs in similar to those reported for several other clusters, though higher than what was found for NGC 1333 in Scholz et al. 2009.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

SONYC -- Substellar Objects in Nearby Young Clusters -- is a survey program to investigate the frequency and properties of substellar objects with masses down to a few times that of Jupiter in nearby star-forming regions. Here we present the first results from SONYC observations of NGC1333, a ~1Myr old cluster in the Perseus star-forming complex. We have carried out extremely deep optical and near-infrared imaging in four bands (i', z', J, K) using Suprime-Cam and MOIRCS instruments at the Subaru telescope. The survey covers 0.25sqdeg and reaches completeness limits of 24.7mag in the i'-band and 20.8mag in the J-band. We select 196 candidates with colors as expected for young, very low-mass objects. Follow-up multi-object spectroscopy with MOIRCS is presented for 53 objects. We confirm 19 objects as likely brown dwarfs in NGC1333, seven of them previously known. For 11 of them, we confirm the presence of disks based on Spitzer/IRAC photometry. The effective temperatures for the brown dwarf sample range from 2500K to 3000K, which translates to masses of ~0.015 to 0.1Ms. For comparison, the completeness limit of our survey translates to mass limits of 0.004Ms for Av<~5mag or 0.008Ms for Av<~ 10mag. Compared with other star-forming regions, NGC1333 shows an overabundance of brown dwarfs relative to low-mass stars, by a factor of 2-5. On the other hand, NGC1333 has a deficit of planetary-mass objects: Based on the surveys in SOrionis, the ONC and Cha I, the expected number of planetary-mass objects in NGC1333 is 8-10, but we find none. It is plausible that our survey has detected the minimum mass limit for star formation in this particular cluster, at around 0.012-0.02Ms. If confirmed, our findings point to significant regional/environmental differences in the number of brown dwarfs and the minimum mass of the IMF. (abridged)

The recent discovery of binary objects in the Kuiper-belt opens an invaluable window into past and present conditions in the trans-Neptunian part of the Solar System. For example, knowledge of how these objects formed can be used to impose constraints on planetary formation theories. We have recently proposed a binary-object formation model based on the notion of chaos-assisted capture. Here we present a more detailed analysis with calculations performed in the spatial (three-dimensional) three- and four-body Hill approximations. It is assumed that the potential binary partners are initially following heliocentric Keplerian orbits and that their relative motion becomes perturbed as these objects undergo close encounters. First, the mass, velocity, and orbital element distribu- tions which favour binary formation are identified in the circular and elliptical Hill limits. We then consider intruder scattering in the circular Hill four-body problem and find that the chaos-assisted capture mechanism is consistent with observed, apparently randomly distributed, binary mutual orbit inclinations. It also predicts asymmetric distributions of retrograde versus prograde orbits. The time-delay induced by chaos on particle transport through the Hill sphere is analogous to the formation of a resonance in a chemical reaction. Implications for binary formation rates are considered and the 'fine-tuning' problem recently identified by Noll et al. (2007) is also addressed.

The discovery that many trans-neptunian objects exist in pairs, or binaries, is proving invaluable for shedding light on the formation, evolution and structure of the outer Solar system. Based on recent systematic searches it has been estimated that up to 10% of Kuiper-belt objects might be binaries. However, all examples discovered to-date are unusual, as compared to near-Earth and main-belt asteroid binaries, for their mass ratios of order unity and their large, eccentric orbits. In this article we propose a common dynamical origin for these compositional and orbital properties based on four-body simulations in the Hill approximation. Our calculations suggest that binaries are produced through the following chain of events: initially, long-lived quasi-bound binaries form by two bodies getting entangled in thin layers of dynamical chaos produced by solar tides within the Hill sphere. Next, energy transfer through gravitational scattering with a low-mass intruder nudges the binary into a nearby non-chaotic, stable zone of phase space. Finally, the binary hardens (loses energy) through a series of relatively gentle gravitational scattering encounters with further intruders. This produces binary orbits that are well fitted by Kepler ellipses. Dynamically, the overall process is strongly favored if the original quasi-bound binary contains comparable masses. We propose a simplified model of chaotic scattering to explain these results. Our findings suggest that the observed preference for roughly equal mass ratio binaries is probably a real effect; that is, it is not primarily due to an observational bias for widely separated, comparably bright objects. Nevertheless, we predict that a sizeable population of very unequal mass Kuiper-belt binaries is likely awaiting discovery.

Single-armed, stationary density waves can propagate in very weakly self-gravitating gas disks dominated by a central mass. Examples include circumstellar disks of protostars and molecular disks in galactic nuclei. We explore the linear and nonlinear dynamics of such waves. Variational methods yield nonlinear versions of the dispersion relation, angular momentum flux, and propagation velocity in the tight-winding limit. The pitch angle increases with amplitude until the tight-winding approximation breaks down. We also find a series of nonlinear logarithmic spirals which is exact in the limit of small disk mass and which extends to large pitch angle.

This paper describes observations of a new phenomenon in evolved mass-losing AGB stars: the presence of two winds with different expansion velocities. CO(2-1) and CO(3-2) line emission was observed for 45 AGB stars at high velocity resolution and double winds found in 20% of the sample. Highly asymmetric lines were found in six other stars. The data tentatively suggest that double winds occur when the star undergoes a change (pulsational mode, chemical composition) and that the very narrow components represent the onset of a new phase of mass loss.