The extreme ultraviolet (EUV) spectra of distant star-forming (SF) regions cannot be probed directly using either ground- or space-based telescopes due to the high cross-section for interaction of EUV photons with the interstellar medium. The mm/submm recombination lines of H and He, which can be observed from the ground can serve as a reliable probe of the EUV. These lines are weekly impacted by dust absorption, their fluxes are independent of metallicity, and depend linearly on the number of EUV continuum photons. We present a study based on ALMA observations of three Galactic ultra-compact HII regions and the starburst region Sgr B2(M), in which we reconstruct the key parameters of the EUV spectra of these sources using mm recombination lines of HI, HeI and HeII. We find that the EUV spectra between 13.6 and 54.4eV of the objects studied here have similar frequency dependence: L_nu~nu^{-4.5 \pm 0.4}. We compare the inferred values of the EUV spectral slopes with the values expected for a single stellar evolution model (Starburst99) and the Binary Population and Spectral Synthesis code (BPASS v2.1). We find that the observed spectral slope gamma~4.5 \pm 0.4 differs from the model predictions. For an instantaneous starburst model with age <5Myrs the observed spectral slopes are consistently steeper than predictions of either of the models. For an instantaneous starburst with age >5Myrs the observed spectral slope is in-between the values predicted by the two models. This may imply that the fraction of interacting binaries in HII regions is substantially lower than assumed in BPASS v2.1. The technique demonstrated here allows one to deduce the EUV spectra of SF regions providing critical insight into photon production rates at lambda>912A and can serve as calibration to starburst synthesis models, improving our understanding of SF and the properties of ionizing flux during reionization.

We address the claim of Weiss et al. (2018) that the radii of adjacent planets in Kepler multi-planet systems are correlated. We explore two simple toy models---in the first the radii of the planets are chosen at random from a single universal distribution, and in the second we postulate several types of system with distinct radius distributions. We show that an apparent correlation between the radii of adjacent planets similar to the one reported by Weiss et al. (2018) can arise in both models. In addition the second model fits all other aspects of the radius distribution, including the signal-to-noise distribution of the observed planets. We also comment on the validity of a commonly used correction that is used to estimate intrinsic planet occurrence rates, based on weighting planets by the inverse of their detectability.

Carefully accounting for neutrino transport is an essential component of many astrophysical studies. Solving the full transport equation is too expensive for most realistic applications, especially those involving multiple spatial dimensions. For such cases, resorting to approximations is often the only viable option for obtaining solutions. One such approximation, which recently became popular, is the M1 method. It utilizes the system of the lowest two moments of the transport equation and closes the system with an ad hoc closure relation. The accuracy of the M1 solution depends on the quality of the closure. Several closures have been proposed in the literature and have been used in various studies. We carry out an extensive study of these closures by comparing the results of M1 calculations with precise Monte Carlo calculations of the radiation field around spherically-symmetric protoneutron star models. We find that no closure performs consistently better or worse than others in all cases. The level of accuracy a given closure yields depends on the matter configuration, neutrino type, and neutrino energy. Given this limitation, the maximum entropy closure by Minerbo (1978) on average yields relatively accurate results in the broadest set of cases considered in this work.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We present 90 mas (37 pc) resolution ALMA imaging of Arp 220 in the CO (1-0) line and continuum at $\lambda = 2.6$ mm. The internal gas distribution and kinematics of both galactic nuclei are well-resolved for the first time. In the West nucleus, the major gas and dust emission extends out to 0.2\arcsec radius (74 pc); the central resolution element shows a strong peak in the dust emission but a factor 3 dip in the CO line emission. In this nucleus, the dust is apparently optically thick ($\tau_{\rm 2.6mm} \sim1$) at $\lambda = 2.6$ mm with a dust brightness temperature $\sim147$ K. The column of ISM at this nucleus is $\rm N_{H2} \geq 2\times10^{26}$ cm$^{-2}$, corresponding to $\sim$900 gr cm$^{-2}$. The East nucleus is more elongated with radial extent 0.3\arcsec or $\sim111$ pc. The derived kinematics of the nuclear disks provide a good fit to the line profiles, yielding the emissivity distributions, the rotation curves and velocity dispersions. In the West nucleus, there is evidence of a central Keplerian component requiring a central mass of $8\times10^8$ \msun. The intrinsic widths of the emission lines are $\Delta \rm v (FWHM)$ = 250 (West) and 120 (East) \kms. Given the very short dissipation timescales for turbulence ($\lesssim10^5$ yrs), we suggest that the line widths may be due to semi-coherent motions within the nuclear disks. The symmetry of the nuclear disk structures is impressive -- implying the merger timescale is significantly longer than the rotation period of the disks.

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

We examine the use of submm recombination lines of H, He and He$^+$ to probe the extreme ultraviolet (EUV) luminosity of starbursts (SB) and AGN. We find that the submm recombination lines of H, He and He$^+$ are in fact extremely reliable and quantitative probes of the EUV continuum at 13.6 eV to above 54.6 eV. At submm wavelengths, the recombination lines originate from low energy levels (n = 20 -- 50). The maser amplification, which poses significant problems for quantitative interpretation of the higher n, radio frequency recombination lines, is insignificant. Lastly, at submm wavelengths the dust extinction is minimal. The submm line luminosities are therefore directly proportional to the emission measures ($EM_{ION} = n_e \times n_{ion} \times \rm{volume}$) of their ionized regions. We also find that the expected line fluxes are detectable with ALMA and can be imaged at $\sim0.1$\arcsec ~resolution in low redshift ULIRGs. Imaging of the HI lines will provide accurate spatial and kinematic mapping of the star formation distribution in low-z IR-luminous galaxies. And the relative fluxes of the HI and HeII recombination lines will strongly constrain the relative contributions of starbursts and AGN to the luminosity. The HI lines should also provide an avenue to constraining the submm dust extinction curve.