We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

We present far-infrared (FIR) properties of an extremely luminous infrared galaxy (ELIRG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1 (WISE1013+6112). This ELIRG is selected as an IR-bright dust-obscured galaxy (DOG) based on the photometry from the Sloan digital sky survey (SDSS) and wide-field infrared survey explorer (WISE). In order to derive its accurate IR luminosity, we perform follow-up observations at 89 and 154 $\mu$m using the high-resolution airborne wideband camera-plus (HAWC+) on board the 2.7-m stratospheric observatory for infrared astronomy (SOFIA) telescope. We conduct spectral energy distribution (SED) fitting with CIGALE using 15 photometric data (0.4-1300 $\mu$m). We successfully pin down FIR SED of WISE1013+6112 and its IR luminosity is estimated to be $L_{\rm IR}$ = (1.62 $\pm$ 0.08) $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the universe. We determine the dust temperature of WISE1013+6112 is $T_{\rm dust}$ = 89 $\pm$ 3 K, which is significantly higher than that of other populations such as SMGs and FIR-selected galaxies at similar IR luminosities. The resultant dust mass is $M_{\rm dust}$ = (2.2 $\pm$ 0.1) $\times 10^{8}$ $M_{\odot}$. This indicates that WISE1013+6112 has a significant active galactic nucleus (AGN) and star-forming activity behind a large amount of dust.

Our current knowledge of cosmic star-formation history during the first two billion years (corresponding to redshift z >3) is mainly based on galaxies identified in rest-frame ultraviolet light. However, this population of galaxies is known to under-represent the most massive galaxies, which have rich dust content and/or old stellar populations. This raises the questions of the true abundance of massive galaxies and the star-formation-rate density in the early universe. Although several massive galaxies that are invisible in the ultraviolet have recently been confirmed at early epochs, most of them are extreme starbursts with star-formation rates exceeding 1000 solar masses per year, suggesting that they are unlikely to represent the bulk population of massive galaxies. Here we report submillimeter (wavelength 870um) detections of 39 massive star-forming galaxies at z > 3, which are unseen in the spectral region from the deepest ultraviolet to the near-infrared. With a space density of about $2 \times 10^{-5}$ per cubic megaparsec (two orders of magnitudes higher than extreme starbursts) and star-formation rates of 200 solar masses per year, these galaxies represent the bulk population of massive galaxies that have been missed from previous surveys. They contribute a total star-formation-rate density ten times larger than that of equivalently massive ultraviolet-bright galaxies at z >3. Residing in the most massive dark matter halos at their redshifts, they are probably the progenitors of the largest present-day galaxies in massive groups and clusters. Such a high abundance of massive and dusty galaxies in the early universe challenges our understanding of massive-galaxy formation.

Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$\beta$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution.

We present structural parameters and morphological properties of faint 450-um selected submillimeter galaxies (SMGs) from the JCMT Large Program, STUDIES, in the COSMOS-CANDELS region. Their properties are compared to an 850um selected and a matched star-forming samples. We investigate stellar structures of 169 faint 450-um sources (S450=2.8-29.6mJy; S/N>4) at z<3 using HST near-infrared observations. Based on our spectral energy distribution fitting, half of such faint SMGs (LIR=10^11.65+-0.98Lsun) lie above the star-formation rate (SFR)/stellar mass plane. The size-mass relation shows that these SMGs are generally similar to less-luminous star-forming galaxies selected by NUV-r vs. r-J colors. Because of the intrinsic luminosity of the sample, their rest-frame optical emission is less extended than the 850um sources (S850>2mJy), and more extended than the star-forming galaxies in the same redshift range. For the stellar mass and SFR matched sample at z~=1 and z~=2, the size differences are marginal between faint SMGs and the matched galaxies. Moreover, faint SMGs have similar Sersic indices and projected axis ratios as star-forming galaxies with the same stellar mass and SFR. Both SMGs and the matched galaxies show high fractions (~70%) of disturbed features at z~=2, and the fractions depend on the SFRs. These suggest that their star formation activity is related to galaxy merging, and the stellar structures of SMGs are similar to those of star-forming galaxies. We show that the depths of submillimeter surveys are approaching the lower luminosity end of star-forming galaxies, allowing us to detect galaxies on the main sequence.

This paper reports the first results of a direct dark matter search with the DEAP-3600 single-phase liquid argon (LAr) detector. The experiment was performed 2 km underground at SNOLAB (Sudbury, Canada) utilizing a large target mass, with the LAr target contained in a spherical acrylic vessel of 3600 kg capacity. The LAr is viewed by an array of PMTs, which would register scintillation light produced by rare nuclear recoil signals induced by dark matter particle scattering. An analysis of 4.44 live days (fiducial exposure of 9.87 tonne-days) of data taken with the nearly full detector during the initial filling phase demonstrates the detector performance and the best electronic recoil rejection using pulse-shape discrimination in argon, with leakage $<1.2\times 10^{-7}$ (90% C.L.) between 16 and 33 keV$_{ee}$. No candidate signal events are observed, which results in the leading limit on WIMP-nucleon spin-independent cross section on argon, $<1.2\times 10^{-44}$ cm$^2$ for a 100 GeV/c$^2$ WIMP mass (90% C.L.).

The Dark matter Experiment using Argon Pulse-shape discrimination (DEAP) has been designed for a direct detection search for particle dark matter using a single-phase liquid argon target. The projected cross section sensitivity for DEAP-3600 to the spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons is $10^{-46}~\rm{cm}^{2}$ for a 100 GeV/$c^2$ WIMP mass with a fiducial exposure of 3 tonne-years. This paper describes the physical properties and construction of the DEAP-3600 detector.

We present the discovery of an extremely-luminous dust-obscured galaxy (DOG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1. This DOG is selected as a candidate of extremely-luminous infrared (IR) galaxies based on the photometry from the Sloan Digital Sky Survey and Wide-field Infrared Survey Explorer. In order to derive its accurate IR luminosity, we perform follow-up observations at 450 and 850 $\mu$m using the Submillimetre Common User Bolometer Array 2 on the James Clerk Maxwell Telescope, and at 870 and 1300 $\mu$m using the Submillimeter Array, which enable us to pin down its IR Spectral Energy Distribution (SED). We perform SED fitting using 14 photometric data (0.4 - 1300 $\mu$m) and estimate its IR luminosity, $L_{\rm IR}$ (8-1000 $\mu$m), to be $2.2^{+1.5}_{-1.0}$ $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the Universe. The energy contribution from an active galactic nucleus (AGN) to the IR luminosity is $94^{+6}_{-20}$%, which indicates it is an AGN-dominated DOG. On the other hand, its stellar mass ($M_*$) and star formation rate (SFR) are $\log \,(M_\ast/M_{\odot})$ = $11.2^{+0.6}_{-0.2}$ and $\log \,({\rm SFR}/M_{\odot}\,{\rm yr}^{-1}$) = $3.1^{+0.2}_{-0.1}$, respectively, which means that this DOG can be considered as a starburst galaxy in $M_*$--SFR plane. This extremely-luminous DOG shows significant AGN and star forming activity that provides us an important laboratory to probe the maximum phase of the co-evolution of galaxies and supermassive black holes.

We present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available.

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

The identification of high-redshift massive galaxies with old stellar populations may pose challenges to some models of galaxy formation. However, to securely classify a galaxy as quiescent, it is necessary to exclude significant ongoing star formation, something that can be challenging to achieve at high redshift. In this letter, we analyse deep ALMA/870um and SCUBA-2/450um imaging of the claimed "post-starburst" galaxy ZF-20115 at z=3.717 that exhibits a strong Balmer break and absorption lines. The far-infrared imaging reveals a luminous starburst located 0.4+/-0.1 arcsec (~3kpc in projection) from the position of the rest-frame ultra-violet/optical emission, with an obscured star-formation rate of 100 Mo/yr. This star-forming component is undetected in the rest-frame ultraviolet but contributes significantly to the lower angular resolution photometry at restframe wavelengths >3500A, significantly complicating the determination of a reliable stellar mass. Importantly, in the presence of dust obscuration, strong Balmer features are not a unique signature of a post-starburst galaxy and are indeed frequently observed in infrared-luminous galaxies. We conclude that the ZF20015 system does not pose a challenge to current models of galaxy formation and that deep sub-/millimeter observations are a prerequisite for any claims of quiescence. The multi-wavelength observations of ZF20115 unveil a complex system with an intricate and spatially-varying star-formation history. ZF20115 demonstrates that understanding high-redshift obscured starbursts will only be possible with multi-wavelength studies that include high-resolution observations, available with the JWST, at mid-infrared wavelengths.

In this paper, we carefully calculated the tensor-to-scalar ratio, the running spectral index, and the running of running spectrum for (extra) natural inflation in order to compare with recent BICEP2 data, PLANCK satellite data and future 21 cm data. We discovered that the prediction for running spectral index and the running of running spectrum in natural inflation is different from that in the case of extra natural inflation. Near future observation for the running spectral index may achieve enough accuracy to allow us distinguishing between extra natural inflation from natural inflation. Distinguishing the models by using the running of running spectrum is not impossible but would be more challenging for future experiments.

In this paper, we consider high scale (100 TeV) supersymmetry (SUSY) breaking and realize the idea of hilltop supernatural inflation in concrete particle physics models based on flipped-SU(5) and Pati-Salam models in the framework of supersymmetric grand unified theories (SUSY GUTs). The inflaton can be a flat direction including right-handed sneutrino and the waterfall field is a GUT Higgs. The spectral index is $n_s=0.96$ which fits very well with recent data by PLANCK satellite. There is no both thermal and non-thermal gravitino problems. Non-thermal leptogenesis can be resulted from the decay of right-handed sneutrino which plays (part of) the role of inflaton.

A family of spin-lattice models are derived as convergent finite dimensional approximations to the rest frame kinetic energy of a barotropic fluid coupled to a massive rotating sphere. In not fixing the angular momentum of the fluid component, there is no Hamiltonian equations of motion of the fluid component of the coupled system. This family is used to formulate a statistical equilibrium model for the energy - relative enstrophy theory of the coupled barotropic fluid - rotating sphere system, known as the spherical model, which because of its microcanonical constraint on relative enstrophy, does not have the low temperature defect of the classical energy - enstrophy theory. This approach differs from previous works and through the quantum - classical mapping between quantum field theory in spatial dimension $d$ and classical statistical mechanics in dimension $d+1,$ provides a new example of Feynman's generalization of the Least Action Principle to problems that do not have a standard Lagrangian or Hamiltonian. A simple mean field theory for this statistical equlibrium model is formulated and solved, providing precise conditions on the planetary spin and relative enstrophy in order for phase transitions to occur at positive and negative critical temperatures, $T_{+}$ and $T_{-}.$

We propose a two-zero-texture general Zee model, compatible with the large mixing angle Mikheyev-Smirnov-Wolfenstein solution. The washing out of the baryon number does not occur in this model for an adequate parameter range. We check the consistency of a model with the constraints coming from flavor changing neutral current processes, the recent cosmic microwave background observation, and the Z-burst scenario.