We present our analysis of the LyC emission and escape fraction of 111 spectroscopically verified galaxies with and without AGN from $2.26<z<4.3$. We extended our ERS sample from Smith et al. (2018; arXiv:1602.01555) with 64 galaxies in the GOODS North and South fields using WFC3/UVIS F225W, F275W, and F336W mosaics we independently drizzled using the HDUV, CANDELS, and UVUDF data. Among the 17 AGN from the 111 galaxies, one provided a LyC detection in F275W at $m_{AB}=23.19$ mag (S/N $\simeq$ 133) and $GALEX$ NUV at $m_{AB}=23.77$ mag (S/N $\simeq$ 13). We simultaneously fit $SDSS$ and $Chandra$ spectra of this AGN to an accretion disk and Comptonization model and find $f_{esc}$ values of $f_{esc}^{F275W}\simeq 28^{+20}_{-4}$% and $f_{esc}^{NUV}\simeq 30^{+22}_{-5}$%. For the remaining 110 galaxies, we stack image cutouts that capture their LyC emission using the F225W, F275W, and F336W data of the GOODS and ERS samples, and both combined, as well as subsamples of galaxies with and without AGN, and $all$ galaxies. We find the stack of 17 AGN dominate the LyC production from $\langle z\rangle\simeq 2.3-4.3$ by a factor of $\sim$10 compared to all 94 galaxies without AGN. While the IGM of the early universe may have been reionized mostly by massive stars, there is evidence that a significant portion of the ionizing energy came from AGN.

We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $\beta$ are in a range of $-3.4\le \beta \le -2.6$. Unlike previous, tentative detections of $\beta \simeq -3$ in photometrically selected, low-luminosity galaxies, our LAEs are spectroscopically confirmed and luminous ($M_{\rm UV}<-20$ mag). We model their broadband spectral energy distributions (SEDs), and find that two $\beta\simeq-2.6\pm0.2$ galaxies can be well fitted with young and dust-free stellar populations. However, it becomes increasingly difficult to fit bluer galaxies. We explore further interpretations by including non-zero LyC escape fraction $f_{\rm esc}$, very low metallicities, and/or AGN contributions. Assuming $f_{\rm esc}\simeq0.2$, we achieve the bluest slopes $\beta\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $\beta\simeq-2.8$ and --2.9 ($\sigma_{\beta}=0.15$). Larger $f_{\rm esc}$ values and very low metallicities are not favored by the strong nebular line emission (evidenced by the IRAC flux) or the observed (IRAC 1 - IRAC 2) color. Finally, we find that the $\beta\simeq-2.9$ galaxy can potentially be well explained by the combination of a very young population with a high $f_{\rm esc}$ ($\ge0.5$) and an old, dusty population. We are not able to produce two $\beta \simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies.

We investigate the accuracy of 4000\AA/Balmer-break based redshifts by combining Hubble Space Telescope ({\it HST}) grism data with photometry. The grism spectra are from the Probing Evolution And Reionization Spectroscopically (PEARS) survey with {\it HST} using the G800L grism on the Advanced Camera for Surveys (ACS). The photometric data come from a compilation by the 3D-HST collaboration of imaging from multiple surveys (notably CANDELS and 3D-HST). We show evidence that spectrophotometric redshifts (SPZs) typically improve the accuracy of photometric redshifts by $\sim$17--60\%. Our SPZ method is a template fitting based routine which accounts for correlated data between neighboring points within grism spectra via the covariance matrix formalism and also accounts for galaxy morphology along the dispersion direction. We show that the robustness of the SPZ is directly related to the fidelity of the D4000 measurement. We also estimate the accuracy of continuum-based redshifts, i.e., for galaxies that do not contain strong emission lines, based on the grism data alone ($\sigma^{\rm NMAD}_{\Delta z/(1+z)}{\lesssim}0.06$). Given that future space-based observatories like WFIRST and Euclid will spend a significant fraction of time on slitless spectroscopic observations, we estimate number densities for objects with $\ | \mathrm{\Delta z/(1+z_{s})} \ | \leq 0.02$. We predict $\sim$700--4400 galaxies/degree$^2$ for galaxies with D4000$>$1.1 and $\ | \mathrm{\Delta z/(1+z_{s})} \ | \leq 0.02$ to a limiting depth of $i_{AB}$=24 mag. This is \emph{especially} important in the absence of an accompanying rich photometric dataset like the existing one for the CANDELS fields, where redshift accuracy from future surveys will rely only on the presence of a feature like the 4000\AA/Balmer breaks or the presence of emission lines within the grism spectra.

We present and analyze spatially-resolved maps for the observed $V$- and $g$-band to 3.6$\mu$m flux ratios and the inferred dust extinction values, $A_V$, for a sample of 257 nearby NGC and IC galaxies. Flux ratio maps are constructed using PSF-matched mosaics of SDSS $g$- and $r$-band images and Spitzer/IRAC 3.6$\mu$m mosaics, with all pixels contaminated by foreground stars or background objects masked out. By applying the $\beta_V$ method (Tamura et al. 2009, 2010), which was recently calibrated as a function of redshift and morphological type by Kim, Jansen, & Windhorst (2017), dust extinction maps were created for each galaxy. The typical 1-$\sigma$ scatter in $\beta_V$ around the average, both within a galaxy and in each morphological type bin, is $\sim$20%. Combined, these result in a $\sim$0.4 mag scatter in $A_V$. $\beta_V$ becomes insensitive to small-scale variations in stellar populations once resolution elements subtend an angle larger than that of a typical giant molecular cloud ($\sim$200pc). We find noticeably redder $V$$-$3.6$\mu$m colors in the center of star-forming galaxies and galaxies with a weak AGN. The derived intrinsic $V$$-$3.6$\mu$m colors for each Hubble type are generally consistent with the model predictions of Kim et al. (2017). Finally, we discuss the applicability of the $\beta_V$ dust-correction method to more distant galaxies, for which well-matched HST rest-frame visible and JWST rest-frame $\sim$3.5$\mu$m images will become available in the near-future.

Research over the past decade has shown diminishing empirical evidence for major galaxy mergers being a dominating or even important mechanism for the growth of supermassive black holes in galaxies and the triggering of optically or X-ray selected active galactic nuclei (AGN). We here for the first time test whether such a connection exists at least in the most plausible part of parameter space for this mechanism: the highest specific accretion rate broad-line AGNs at the peak epoch of black hole activity around z = 2. To that end we examine 21 galaxies hosting a high accreting black hole (L/Ledd > 0.7) observed with HST/WFC3 and 92 stellar mass- and redshift- matched inactive galaxies taken from the CANDELS survey. We removed the AGN point sources from their host galaxies and avoided bias in visual classification by adding and then subtracting mock point sources to and from the comparison galaxies, producing matched residual structures for both sets. The resulting samples were joined, randomized, and subsequently visually ranked with respect to perceived strength of structural distortions by 10 experts. The ensuing individual rankings were combined into a consensus sequence and from this we derived merger fractions for both samples. With the merger fractions f$_{m,agn}$ = 0.24 $\pm$ 0.09 for the AGN host galaxy sample and f$_{m,ina}$ = 0.19 $\pm$ 0.04 for the inactive galaxies, we find no significant difference between the AGN host galaxies and inactive galaxies. Also, both samples display comparable fractions of disk-dominated galaxies. These findings are consistent with previous studies for different AGN populations, and we conclude that even black hole growth at the highest specific accretion rates and at the peak of cosmic AGN activity is not predominantly caused by major mergers. (abriged)

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

We present Hubble Space Telescope WFC3-IR imaging in the fields of six apparently bright dusty star-forming galaxies (DSFGs) at $z$ = 2-4 identified by their rest-frame far-infrared colors using the Planck and Herschel space facilities. We detect near-infrared counterparts for all six submillimeter sources, allowing us to undertake strong-lensing analyses. One field in particular stands out for its prominent giant arcs, PLCK G165.7+67.0 (G165). After combining the color and morphological information, we identify 11 sets of image multiplicities in this one field. We construct a strong-lensing model constrained by this lensing evidence, which uncovers a bimodal spatial mass distribution, and from which we measure a mass of $(2.6 \pm 0.11)$ $\times$ $10^{14}$ $M_{\odot}$ within $\sim$250 kpc. The bright ($S_{350}$ $\approx$ 750 mJy) DSFG appears as two images: a giant arc with a spatial extent of 4.5" that is merging with the critical curve, and a lower-magnification counterimage that is detected in our new longer-wavelength ground- and space-based imaging data. Using our ground-based spectroscopy, we calculate a dynamical mass of $1.3^{+0.04}_{-0.70} \times 10^{15}$ $M_{\odot}$ to the same fixed radius, although this value may be inflated relative to the true value if the velocity distribution is enhanced in the line-of-sight direction. We suggest that the bimodal mass taken in combination with the weak X-ray flux and low SZ decrement may be explained as a pre-merger for which the intracluster gas is diluted along the line of sight, while the integrated surface mass density is supercritical to strong-lensing effects.

Direct observational searches for Population III (Pop III) stars at high redshift are faced with the question of how to select the most promising targets for spectroscopic follow-up. To help answer this, we use a large-scale cosmological simulation, augmented with a new subgrid model that tracks the fraction of pristine gas, to follow the evolution of high-redshift galaxies and the Pop III stars they contain. We generate rest-frame ultraviolet (UV) luminosity functions for our galaxies and find that they are consistent with current $z \ge 7 $ observations. Throughout the redshift range $7 \le z \le 15$ we identify "Pop III-bright" galaxies as those with at least 75% of their flux coming from Pop III stars. While less than 1% of galaxies brighter than $m_{\rm UV, AB} = 31.4$ mag are Pop III--bright in the range $7\leq z \leq8$, roughly 17% of such galaxies are Pop III--bright at $z=9$, immediately before reionization occurs in our simulation. Moving to $z=10$, $m_{\rm UV, AB} = 31.4$ mag corresponds to larger, more luminous galaxies and the Pop III-bright fraction falls off to 5%. Finally, at the highest redshifts, a large fraction (29% at $z=14$ and 41% at $z=15)$ of all galaxies are Pop III-bright regardless of magnitude. While $m_{\rm UV, AB} = 31.4$ mag galaxies are extremely rare during this epoch, we find that 13% of galaxies at $z = 14$ are Pop III-bright with $m_{\rm UV, AB} \le 33$ mag, an intrisic magnitude within reach of the James Webb Space Telescope using lensing. Thus, we predict that the best redshift to search for luminous Pop III--bright galaxies is just before reionization, while lensing surveys for fainter galaxies should push to the highest redshifts possible.

We present observations of escaping Lyman Continuum (LyC) radiation from 34 massive star-forming galaxies and 12 weak AGN with reliably measured spectroscopic redshifts at $z$$\simeq$2.3-4.1. We analyzed Hubble Space Telescope Wide Field Camera 3 (WFC3) mosaics of the Early Release Science field in three UVIS filters to sample the rest-frame LyC over this redshift range. With our best current assessment of the WFC3 systematics, we provide 1$\sigma$ upper limits for the average LyC emission of galaxies at $z$=2.35, 2.75, and 3.60 to $\sim$28.5, 28.1, and 30.7 mag in image stacks of 11-15 galaxies in the WFC3/UVIS F225W, F275W, and F336W, respectively. The LyC flux of weak AGN at $z$=2.62 and 3.32 are detected at 28.3 and 27.4 mag with SNRs of $\sim$2.7 and 2.5 in F275W and F336W for stacks of 7 and 3 AGN, respectively, while AGN at $z$=2.37 are constrained to $\gtrsim$27.9 mag at 1$\sigma$ in a stack of 2 AGN. The stacked AGN LyC light profiles are flatter than their corresponding non-ionizing UV continuum profiles out to radii of r$\lesssim$0."9, which may indicate a radial dependence of porosity in the ISM. With synthetic stellar SEDs fit to UV continuum measurements longwards of Ly$\alpha$ and IGM transmission models, we constrain the absolute LyC escape fractions to $f_{\rm esc}^{\rm abs}$$\simeq$$22^{+44}_{-22}$% at $z$=2.35 and $\lesssim$55% at $z$=2.75 and 3.60, respectively. All available data for galaxies, including published work, suggests a more sudden increase of $f_{\rm esc}$ with redshift at $z$$\simeq$2. Dust accumulating in (massive) galaxies over cosmic time correlates with increased HI column density, which may lead to reducing $f_{\rm esc}$ more suddenly at $z$$\lesssim$2. This may suggest that star-forming galaxies collectively contributed to maintaining cosmic reionization at redshifts $z$$\gtrsim$2-4, while AGN likely dominated reionization at $z$$\lesssim$2.

We present a Monte Carlo simulation designed to predict the vertical velocity dispersion of brown dwarfs in the Milky Way. We show that since these stars are constantly cooling, the velocity dispersion has a noticeable trend with spectral type. With realistic assumptions for the initial-mass function, star-formation history, and the cooling models, we show that the velocity dispersion is roughly consistent with what is observed for M dwarfs, decreases to cooler spectral types, and increases again for the coolest types in our study ($\sim$T9). We predict a minimum in the velocity dispersions for L/T transition objects, however the detailed properties of the minimum predominately depend on the star-formation history. Since this trend is due to brown dwarf cooling, we expect the velocity dispersion as a function of spectral type should deviate from constancy around the hydrogen-burning limit. We convert from velocity dispersion to vertical scale height using standard disk models, and present similar trends in disk thickness as a function of spectral type. We suggest that future, wide-field photometric and/or spectroscopic missions may collect sizable samples of distant ($\sim\!1$ kpc) of dwarfs that span the hydrogen-burning limit. As such, we speculate that such observations may provide a unique way of constraining the average spectral type of hydrogen-burning.

Recent work suggests that strong emission line, star-forming galaxies may be significant Lyman Continuum leakers. We combine archival HST broadband ultraviolet and optical imaging (F275W and F606W, respectively) with emission line catalogs derived from WFC3 IR G141 grism spectroscopy to search for escaping Lyman Continuum (LyC) emission from homogeneously selected $z\simeq$2.5 SFGs. We detect no escaping Lyman Continuum from SFGs selected on [OII] nebular emission (N=208) and, within a narrow redshift range, on [OIII]/[OII]. We measure 1$\sigma$ upper limits to the LyC escape fraction relative to the non-ionizing UV continuum from [OII] emitters, $f_{esc}<$5.6%, and strong [OIII]/[OII]$>$5 ELGs, $f_{esc}<$14.0%. Our observations are not deep enough to detect $f_{esc}\lesssim$10% typical of the low redshift Lyman continuum emitters. However, we find that this population represents a small fraction of the star-forming galaxy population at $z\simeq$2. Thus, unless the number of extreme emission line galaxies grows substantially to z$>$6, such galaxies may be insufficient for reionization. Deeper survey data in the rest-frame ionizing UV will be necessary to determine whether strong line ratios could be useful for pre-selecting LyC leakers at high redshift.

We present a study of the trade-off between depth and resolution using a large number of U-band imaging observations in the GOODS-North field (Giavalisco et al. 2004) from the Large Binocular Camera (LBC) on the Large Binocular Telescope (LBT). Having acquired over 30 hours of data (315 images with 5-6 mins exposures), we generated multiple image mosaics, starting with the best atmospheric seeing images (FWHM $\lesssim$0.8"), which constitute $\sim$10% of the total data set. For subsequent mosaics, we added in data with larger seeing values until the final, deepest mosaic included all images with FWHM $\lesssim$1.8" ($\sim$94% of the total data set). From the mosaics, we made object catalogs to compare the optimal-resolution, yet shallower image to the lower-resolution but deeper image. We show that the number counts for both images are $\sim$90% complete to $U_{AB}$ $\lesssim26$. Fainter than $U_{AB}$$\sim$ 27, the object counts from the optimal-resolution image start to drop-off dramatically (90% between $U_{AB}$ = 27 and 28 mag), while the deepest image with better surface-brightness sensitivity ($\mu^{AB}_{U}$$\lesssim$ 32 mag arcsec$^{-2}$) show a more gradual drop (10% between $U_{AB}$ $\simeq$ 27 and 28 mag). For the brightest galaxies within the GOODS-N field, structure and clumpy features within the galaxies are more prominent in the optimal-resolution image compared to the deeper mosaics. Finally, we find - for 220 brighter galaxies with $U_{AB}$$\lesssim$ 24 mag - only marginal differences in total flux between the optimal-resolution and lower-resolution light-profiles to $\mu^{AB}_{U}$$\lesssim$ 32 mag arcsec$^{-2}$. In only 10% of the cases are the total-flux differences larger than 0.5 mag. This helps constrain how much flux can be missed from galaxy outskirts, which is important for studies of the Extragalactic Background Light.

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

We analyze the relationship between star formation (SF), substructure, and supercluster environment in a sample of 107 nearby galaxy clusters using data from the Sloan Digital Sky Survey. Previous works have investigated the relationships between SF and cluster substructure, and cluster substructure and supercluster environment, but definitive conclusions relating all three of these variables has remained elusive. We find an inverse relationship between cluster SF fraction (f_SF) and supercluster environment density, calculated using the galaxy luminosity density field at a smoothing length of 8 h^-1 Mpc (D8). The slope of f_SF vs. D8 is -0.008 +/- 0.002. The f_SF of clusters located in low-density large-scale environments, 0.244 +/- 0.011, is higher than for clusters located in high-density supercluster cores, 0.202 +/- 0.014. We also divide superclusters, according to their morphology, into filament- and spider-type systems. The inverse relationship between cluster f_SF and large-scale density is dominated by filament- rather than spider-type superclusters. In high-density cores of superclusters, we find a higher f_SF in spider-type superclusters, 0.229 +/- 0.016, than in filament-type superclusters, 0.166 +/- 0.019. Using principal component analysis, we confirm these results and the direct correlation between cluster substructure and SF. These results indicate that cluster SF is affected by both the dynamical age of the cluster (younger systems exhibit higher amounts of SF); the large-scale density of the supercluster environment (high-density core regions exhibit lower amounts of SF); and supercluster morphology (spider-type superclusters exhibit higher amounts of SF at high densities).

We directly measure the thermal energy of the gas surrounding galaxies through the thermal Sunyaev-Zel'dovich (tSZ) effect. We perform a stacking analysis of microwave background images from the Atacama Cosmology Telescope, around 1179 massive quiescent elliptical galaxies at 0.5 <= z <= 1.0 ('low-z') and 3274 galaxies at 1.0 <= z <= 1.5 ('high-z'), selected using data from the Wide-Field Infrared Survey Explorer All-Sky Survey and the Sloan Digital Sky Survey (SDSS) within the SDSS Stripe-82 field. The gas surrounding these galaxies is expected to contain energy from past episodes of active galactic nucleus (AGN) feedback, and after using modeling to subtract undetected contaminants, we detect a tSZ signal at a significance of 0.9-sigma for our low-z galaxies and 1.8-sigma for our high-z galaxies. We then include data from the high-frequency Planck bands for a subset of 227 low-z galaxies and 529 high-z galaxies and find low-z and high-z tSZ detections of 1.0-sigma and 1.5-sigma, respectively. These results indicate an average thermal heating around these galaxies of 5.6(+5.9/-5.6) x 10^60 erg for our low-z galaxies and 7.0(+4.7/-4.4) x 10^60 erg for our high-z galaxies. Based on simple heating models, these results are consistent with gravitational heating without additional heating due to AGN feedback.

The most frequently proposed model for the origin of quasars holds that the high accretion rates seen in luminous active galactic nuclei are primarily triggered during major mergers between gas-rich galaxies. While plausible for decades, this model has only begun to be tested with statistical rigor in the past few years. Here we report on a Hubble Space Telescope study to test this hypothesis for $z=2$ quasars with high super-massive black hole masses ($M_\mathrm{BH}=10^9-10^{10}~M_\odot{}$), which dominate cosmic black hole growth at this redshift. We compare Wide Field Camera 3 $F160W$ (rest-frame $V$-band) imaging of 19 point source-subtracted quasar hosts to a matched sample of 84 inactive galaxies, testing whether the quasar hosts have greater evidence for strong gravitational interactions. Using an expert ranking procedure, we find that the quasar hosts are uniformly distributed within the merger sequence of inactive galaxies, with no preference for quasars in high-distortion hosts. Using a merger/non-merger cutoff approach, we recover distortion fractions of $f_\mathrm{m,qso}=0.39\pm{}0.11$ for quasar hosts and $f_\mathrm{m,gal}=0.30\pm{}0.05$ for inactive galaxies (distribution modes, 68% confidence intervals), with both measurements subjected to the same observational conditions and limitations. The slight enhancement in distorted fraction for quasar hosts over inactive galaxies is not significant, with a probability that the quasar fraction is higher of $P(f_\mathrm{m,qso}>f_\mathrm{m,gal}) = 0.78$ ($0.78\,\sigma{}$), in line with results for lower mass and lower $z$ AGN. We find no evidence that major mergers are the primary triggering mechanism for the massive quasars that dominate accretion at the peak of cosmic quasar activity.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

Energetic feedback due to active galactic nuclei (AGN) is likely to play an important role in the observed anti-hierarchical trend in the evolution of galaxies, and yet the energy injected into the circumgalactic medium by this process is largely unknown. One promising approach to constrain this feedback is through measurements of CMB spectral distortions due to the thermal Sunyaev-Zel'dovich (tSZ) effect, whose magnitude is directly proportional to the energy input by AGN. Here we co-add South Pole Telescope SZ (SPT-SZ) survey data around a large set of massive quiescent elliptical galaxies at z >= 0.5. We use data from the Blanco Cosmology Survey and VISTA Hemisphere Survey to create a large catalog of galaxies split up into two redshift bins, with 3394 galaxies at 0.5 <= z <= 1.0 and 924 galaxies at 1.0 <= z <= 1.5, with typical stellar masses of 1.5 x 10^11 M_Sun. We then co-add the emission around these galaxies, resulting in a measured tSZ signal at 2.2 sigma significance for the lower redshift bin and a contaminating signal at 1.1 sigma for the higher redshift bin. To remove contamination due to dust emission, we use SPT-SZ source counts to model a contaminant source population in both the SPT-SZ bands and Planck high-frequency bands for a subset of 937 low-redshift galaxies and 240 high-redshift galaxies. This increases our detection to 3.6 sigma for low redshifts and 0.9 sigma for high redshifts. We find the mean angularly-integrated Compton-y values to be 2.2 (-0.7+0.9) x 10^-7 Mpc^2 for low redshifts and 1.7 (-1.8+2.2) x 10^-7 Mpc^2 for high redshifts, corresponding to total thermal energies of 7.6 (-2.3+3.0) x 10^60 ergs and 6.0 (-6.3+7.7) x 10^60 ergs, respectively. These numbers are higher than expected from simple theoretical models that do not include AGN feedback, and serve as constraints that can be applied to current simulations of massive galaxy formation. (abridged)

We combine wide and deep galaxy number-count data from GAMA, COSMOS/G10, HST ERS, HST UVUDF and various near-, mid- and far- IR datasets from ESO, Spitzer and Herschel. The combined data range from the far-UV (0.15microns) to far-IR (500microns), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive the IGL and the extrapolated-IGL in all bands. We argue undetected low surface brightness galaxies and intra-cluster/group light is modest, and that our extrapolated-IGL measurements are an accurate representation of the extra-galactic background light. Our data agree with most earlier IGL estimates and with direct measurements in the far-IR, but disagree strongly with direct estimates in the optical. Close agreement between our results and recent very high-energy experiments (H.E.S.S. and MAGIC), suggest that there may be an additional foreground affecting the direct estimates. The most likely culprit could be the adopted Zodiacal light model. Finally we use a modified version of the two-component model to integrate the EBL and obtain measurements of the Cosmic Optical Background (COB) and Cosmic Infrared Background (CIB) of: $24^{+4}_{-4}$nWm$^{-2}$sr$^{-1}$ and $26^{+5}_{-5}$nWm$^{-2}$sr$^{-1}$ respectively (48:52%). Over the next decade, upcoming space missions such as Euclid and WFIRST, have the capacity to reduce the COB error to $<1\%$, at which point comparisons to the very high energy data could have the potential to provide a direct detection and measurement of the reionisation field.

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$\alpha$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$\alpha$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$\alpha$ emission are indistinguishable in terms of age, stellar mass, and SFR.

We have performed a search for muon neutrinos from dark matter annihilation in the center of the Sun with the 79-string configuration of the IceCube neutrino telescope. For the first time, the DeepCore sub-array is included in the analysis, lowering the energy threshold and extending the search to the austral summer. The 317 days of data collected between June 2010 and May 2011 are consistent with the expected background from atmospheric muons and neutrinos. Upper limits are therefore set on the dark matter annihilation rate, with conversions to limits on spin-dependent and spin-independent WIMP-proton cross-sections for WIMP masses in the range 20 - 5000 GeV. These are the most stringent spin-dependent WIMP-proton cross-sections limits to date above 35 GeV.

We investigate the relationship between star formation (SF) and level of relaxation in a sample of 379 galaxy clusters at z < 0.2. We use data from the Sloan Digital Sky Survey to measure cluster membership and level of relaxation, and to select star-forming galaxies based on mid-infrared emission detected with the Wide-Field Infrared Survey Explorer. For galaxies with absolute magnitudes M_r < -19.5, we find an inverse correlation between SF fraction and cluster relaxation: as a cluster becomes less relaxed, its SF fraction increases. Furthermore, in general, the subtracted SF fraction in all unrelaxed clusters (0.117 +/- 0.003) is higher than that in all relaxed clusters (0.097 +/- 0.005). We verify the validity of our SF calculation methods and membership criteria through analysis of previous work. Our results agree with previous findings that a weak correlation exists between cluster SF and dynamical state, possibly because unrelaxed clusters are less evolved relative to relaxed clusters.

We present an analysis of the stellar populations of 102 visually-selected early-type galaxies (ETGs) with spectroscopic redshifts (0.3<z<1.5) from observations in the Early Release Science program with the Wide Field Camera 3 (WFC3) on {\it Hubble Space Telescope} (HST). We fit one- and two-component synthetic stellar models to the ETGs UV-optical-near-IR spectral energy distributions and find a large fraction (~40%) are likely to have experienced a minor (f$\lesssim$10% of stellar mass) burst of recent (t$\lesssim$1 Gyr) star-formation. The measured ages and mass fraction of the young stellar populations do not strongly trend with measurements of galaxy morphology. We note that massive (log(M[$M_{\odot}$])>10.5) recently star-forming ETGs appear to have larger sizes. Furthermore, high-mass, quiescent ETGs identified with likely companions populate a distinct region in the size-mass parameter space, in comparison with the distribution of massive ETGs with evidence of RSF. We conclude that both mechanisms of the quenching of star-formation in disk-like ETGs and (gas-rich, minor) merger activity contribute to the formation of young stars and the size-mass evolution of intermediate redshift ETGs. The number of ETGs for which we have both HST WFC3 panchromatic (especially UV) imaging and spectroscopically-confirmed redshifts is relatively small, therefore a conclusion on the relative roles of both of these mechanisms remains an open question.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

We study of the role of major mergers (mass ratios >1:4) in driving size growth in high-redshift (1<z<2) spheroidal galaxies (SGs) with stellar masses between 10^9.5 MSun and 10^10.7 MSun. This is a largely unexplored mass range at this epoch, containing the progenitors of more massive SGs on which the bulk of the size-evolution literature is based. We visually split our SGs into systems that are relaxed and those that exhibit tidal features indicative of a recent merger. Numerical simulations indicate that, given the depth of our images, only tidal features due to major mergers will be detectable at these epochs (features from minor mergers being too faint), making the disturbed SGs a useful route for estimating major-merger-driven size growth. The disturbed SGs are offset in size from their relaxed counterparts, lying close to the upper envelope of the local size -- mass relation. The mean size ratio of the disturbed SGs to their relaxed counterparts is ~2. Combining this observed size growth with empirical major-merger histories from the literature suggests that the size evolution of a significant fraction (around two-thirds) of SGs in this mass range could be driven by major mergers. If, as is likely, our galaxies are progenitors of more massive M > 10^10.7 MSun SGs at z<1, then major mergers are also likely to play an important role in the size growth of at least some massive SGs in this mass range.

The IceCube Neutrino Observatory is a large Cherenkov detector instrumenting $1\,\mathrm{km}^3$ of Antarctic ice. The detector can be used to search for signatures of particle physics beyond the Standard Model. Here, we describe the search for non-relativistic, magnetic monopoles as remnants of the GUT (Grand Unified Theory) era shortly after the Big Bang. These monopoles may catalyze the decay of nucleons via the Rubakov-Callan effect with a cross section suggested to be in the range of $10^{-27}\,\mathrm{cm^2}$ to $10^{-21}\,\mathrm{cm^2}$. In IceCube, the Cherenkov light from nucleon decays along the monopole trajectory would produce a characteristic hit pattern. This paper presents the results of an analysis of first data taken from May 2011 until May 2012 with a dedicated slow-particle trigger for DeepCore, a subdetector of IceCube. A second analysis provides better sensitivity for the brightest non-relativistic monopoles using data taken from May 2009 until May 2010. In both analyses no monopole signal was observed. For catalysis cross sections of $10^{-22}\,(10^{-24})\,\mathrm{cm^2}$ the flux of non-relativistic GUT monopoles is constrained up to a level of $\Phi_{90} \le 10^{-18}\,(10^{-17})\,\mathrm{cm^{-2}s^{-1}sr^{-1}}$ at a 90% confidence level, which is three orders of magnitude below the Parker bound. The limits assume a dominant decay of the proton into a positron and a neutral pion. These results improve the current best experimental limits by one to two orders of magnitude, for a wide range of assumed speeds and catalysis cross sections.

[abridged] New near-infrared surveys, using the HST, offer an unprecedented opportunity to study rest-frame optical galaxy morphologies at z>1 and to calibrate automated morphological parameters that will play a key role in classifying future massive datasets like EUCLID or LSST. We study automated parameters (e.g. CAS, Gini, M20) of massive galaxies at 1<z<3, measure their dependence on wavelength and evolution with redshift and quantify the reliability of these parameters in discriminating between visually-determined morphologies, using machine learning algorithms. We find that the relative trends between morphological types observed in the low-redshift literature are preserved at z>1: bulge-dominated systems have systematically higher concentration and Gini coefficients and are less asymmetric and rounder than disk-dominated galaxies. However, at z>1, galaxies are, on average, 50% more asymmetric and have Gini and M20 values that are 10% higher and 20% lower respectively. In bulge-dominated galaxies, morphological parameters derived from the rest-frame UV and optical wavelengths are well correlated; however late-type galaxies exhibit higher asymmetry and clumpiness when measured in the rest-frame UV. We find that broad morphological classes (e.g. bulge vs. disk dominated) can be distinguished using parameters with high (80%) purity and completeness of 80%. In a similar vein, irregular disks and mergers can also be distinguished from bulges and regular disks with a contamination lower than 20%. However, mergers cannot be differentiated from the irregular morphological class using these parameters, due to increasingly asymmetry of non-interacting late-type galaxies at z>1. Our automated procedure is applied to the CANDELS GOODS-S field and compared with the visual classification recently released on the same area getting similar results.

We present the results of a first search for self-annihilating dark matter in nearby galaxies and galaxy clusters using a sample of high-energy neutrinos acquired in 339.8 days of live time during 2009/10 with the IceCube neutrino observatory in its 59-string configuration. The targets of interest include the Virgo and Coma galaxy clusters, the Andromeda galaxy, and several dwarf galaxies. We obtain upper limits on the cross section as a function of the weakly interacting massive particle mass between 300 GeV and 100 TeV for the annihilation into b bbar, W+W-, \tau+\tau-, \mu+\mu-, and \nu \nu bar. A limit derived for the Virgo cluster, when assuming a large effect from subhalos, challenges the weak interacting massive particle interpretation of a recently observed GeV positron excess in cosmic rays.

Accurate measurement of neutrino energies is essential to many of the scientific goals of large-volume neutrino telescopes. The fundamental observable in such detectors is the Cherenkov light produced by the transit through a medium of charged particles created in neutrino interactions. The amount of light emitted is proportional to the deposited energy, which is approximately equal to the neutrino energy for $\nu_e$ and $\nu_\mu$ charged-current interactions and can be used to set a lower bound on neutrino energies and to measure neutrino spectra statistically in other channels. Here we describe methods and performance of reconstructing charged-particle energies and topologies from the observed Cherenkov light yield, including techniques to measure the energies of uncontained muon tracks, achieving average uncertainties in electromagnetic-equivalent deposited energy of $\sim 15\%$ above 10 TeV.

We investigate the relationship between star formation (SF) and substructure in a sample of 107 nearby galaxy clusters using data from the Sloan Digital Sky Survey (SDSS). Several past studies of individual galaxy clusters have suggested that cluster mergers enhance cluster SF, while others find no such relationship. The SF fraction in multi-component clusters (0.228 +/- 0.007) is higher than that in single-component clusters (0.175 +/- 0.016) for galaxies with M^0.1_r < -20.5. In both single- and multi-component clusters, the fraction of star-forming galaxies increases with clustercentric distance and decreases with local galaxy number density, and multi-component clusters show a higher SF fraction than single-component clusters at almost all clustercentric distances and local densities. Comparing the SF fraction in individual clusters to several statistical measures of substructure, we find weak, but in most cases significant at greater than 2 sigma, correlations between substructure and SF fraction. These results could indicate that cluster mergers may cause weak but significant SF enhancement in clusters, or unrelaxed clusters exhibit slightly stronger SF due to their less evolved states relative to relaxed clusters.

Papers on atmospheric and diffuse UHE neutrino searches of all flavors submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

In this paper, we present a novel artificial intelligence (AI) program that identifies pulsars from recent surveys using image pattern recognition with deep neural nets---the PICS (Pulsar Image-based Classification System) AI. The AI mimics human experts and distinguishes pulsars from noise and interferences by looking for patterns from candidate. The information from each pulsar candidate is synthesized in four diagnostic plots, which consist of up to thousands pixel of image data. The AI takes these data from each candidate as its input and uses thousands of such candidates to train its ~9000 neurons. Different from other pulsar selection programs which use pre-designed patterns, the PICS AI teaches itself the salient features of different pulsars from a set of human-labeled candidates through machine learning. The deep neural networks in this AI system grant it superior ability in recognizing various types of pulsars as well as their harmonic signals. The trained AI's performance has been validated with a large set of candidates different from the training set. In this completely independent test, PICS ranked 264 out of 277 pulsar-related candidates, including all 56 previously known pulsars, to the top 961 (1%) of 90008 test candidates, missing only 13 harmonics. The first non-pulsar candidate appears at rank 187, following 45 pulsars and 141 harmonics. In other words, 100% of the pulsars were ranked in the top 1% of all candidates, while 80% were ranked higher than any noise or interference. The performance of this system can be improved over time as more training data are accumulated. This AI system has been integrated into the PALFA survey pipeline and has discovered six new pulsars to date.

We report on results of an all-sky search for high-energy neutrino events interacting within the IceCube neutrino detector conducted between May 2010 and May 2012. The search follows up on the previous detection of two PeV neutrino events, with improved sensitivity and extended energy coverage down to approximately 30 TeV. Twenty-six additional events were observed, substantially more than expected from atmospheric backgrounds. Combined, both searches reject a purely atmospheric origin for the twenty-eight events at the $4\sigma$ level. These twenty-eight events, which include the highest energy neutrinos ever observed, have flavors, directions, and energies inconsistent with those expected from the atmospheric muon and neutrino backgrounds. These properties are, however, consistent with generic predictions for an additional component of extraterrestrial origin.

We present HST/WFC3 narrow-band imaging of the starburst galaxy M83 targeting the hydrogen recombination lines (H$\beta$, H$\alpha$ and Pa$\beta$), which we use to investigate the dust extinction in the HII regions. We derive extinction maps with 6 parsec spatial resolution from two combinations of hydrogen lines (H$\alpha$/H$\beta$ and H$\alpha$/Pa$\beta$), and show that the longer wavelengths probe larger optical depths, with $A_V$ values larger by $\gtrsim$1 mag than those derived from the shorter wavelengths. This difference leads to a factor $\gtrsim$2 discrepancy in the extinction-corrected H$\alpha$ luminosity, a significant effect when studying extragalactic HII regions. By comparing these observations to a series of simple models, we conclude that a large diversity of absorber/emitter geometric configurations can account for the data, implying a more complex physical structure than the classical foreground "dust screen" assumption. However, most data points are bracketed by the foreground screen and a model where dust and emitters are uniformly mixed. When averaged over large ($\gtrsim$100--200 pc) scales, the extinction becomes consistent with a "dust screen", suggesting that other geometries tend to be restricted to more local scales. Moreover, the extinction in any region can be described by a combination of the foreground screen and the uniform mixture model with weights of 1/3 and 2/3 in the center ($\lesssim$2 kpc), respectively, and 2/3 and 1/3 for the rest of the disk. This simple prescription significantly improves the accuracy of the dust extinction corrections and can be especially useful for pixel-based analyses of galaxies similar to M83.

We report on the search for neutrino-induced particle-showers, so-called cascades, in the IceCube-40 detector. The data for this search was collected between April 2008 and May 2009 when the first 40 IceCube strings were deployed and operational. Three complementary searches were performed, each optimized for different energy regimes. The analysis with the lowest energy threshold (2 TeV) targeted atmospheric neutrinos. A total of 67 events were found, consistent with the expectation of 41 atmospheric muons and 30 atmospheric neutrino events. The two other analyses targeted a harder, astrophysical neutrino flux. The analysis with an intermediate threshold of 25 TeV lead to the observation of 14 cascade-like events, again consistent with the prediction of 3.0 atmospheric neutrino and 7.7 atmospheric muon events. We hence set an upper limit of $E^2 \Phi_{lim} \leq 7.46\times10^{-8}\,\mathrm{GeV sr^{-1} s^{-1} cm^{-2}}$ (90% C.L.) on the diffuse flux from astrophysical neutrinos of all neutrino flavors, applicable to the energy range 25 TeV to 5 PeV, assuming an $E_{\nu}^{-2}$ spectrum and a neutrino flavor ratio of 1:1:1 at the Earth. The third analysis utilized a larger and optimized sample of atmospheric muon background simulation, leading to a higher energy threshold of 100 TeV. Three events were found over a background prediction of 0.04 atmospheric muon events and 0.21 events from the flux of conventional and prompt atmospheric neutrinos. Including systematic errors this corresponds to a $2.7\sigma$ excess with respect to the background-only hypothesis. Our observation of neutrino event candidates above 100 TeV complements IceCube's recently observed evidence for high-energy astrophysical neutrinos.

We have searched for extremely high energy neutrinos using data taken with the IceCube detector between May 2010 and May 2012. Two neutrino induced particle shower events with energies around 1 PeV were observed, as reported previously. In this work, we investigate whether these events could originate from cosmogenic neutrinos produced in the interactions of ultra-high energy cosmic-rays with ambient photons while propagating through intergalactic space. Exploiting IceCube's large exposure for extremely high energy neutrinos and the lack of observed events above 100 PeV, we can rule out the corresponding models at more than 90% confidence level. The model independent quasi-differential 90% CL upper limit, which amounts to $E^2 \phi_{\nu_e + \nu_\mu + \nu_\tau} = 1.2 \times 10^{-7}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ at 1 EeV, provides the most stringent constraint in the energy range from 10 PeV to 10 EeV. Our observation disfavors strong cosmological evolution of the highest energy cosmic ray sources such as the Fanaroff-Riley type II class of radio galaxies.

A search for high-energy neutrinos was performed using data collected by the IceCube Neutrino Observatory from May 2009 to May 2010, when the array was running in its 59-string configuration. The data sample was optimized to contain muon neutrino induced events with a background contamination of atmospheric muons of less than 1%. These data, which are dominated by atmospheric neutrinos, are analyzed with a global likelihood fit to search for possible contributions of prompt atmospheric and astrophysical neutrinos, neither of which have yet been identified. Such signals are expected to follow a harder energy spectrum than conventional atmospheric neutrinos. In addition, the zenith angle distribution differs for astrophysical and atmospheric signals. A global fit of the reconstructed energies and directions of observed events is performed, including possible neutrino flux contributions for an astrophysical signal and atmospheric backgrounds as well as systematic uncertainties of the experiment and theoretical predictions. The best fit yields an astrophysical signal flux for $\nu_\mu + \bar\nu_\mu $ of $E^2 \cdot \Phi (E) = 0.25 \cdot 10^{-8} \textrm{GeV} \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{sr}^{-1}$, and a zero prompt component. Although the sensitivity of this analysis for astrophysical neutrinos surpasses the Waxman and Bahcall upper bound, the experimental limit at 90% confidence level is a factor of 1.5 above at a flux of $E^2 \cdot \Phi (E) = 1.44 \cdot 10^{-8} \textrm{GeV} \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{sr}^{-1}$.

The IceCube project has transformed one cubic kilometer of deep natural Antarctic ice into a Cherenkov detector. Muon neutrinos are detected and their direction inferred by mapping the light produced by the secondary muon track inside the volume instrumented with photomultipliers. Reconstructing the muon track from the observed light is challenging due to noise, light scattering in the ice medium, and the possibility of simultaneously having multiple muons inside the detector, resulting from the large flux of cosmic ray muons. This manuscript describes work on two problems: (1) the track reconstruction problem, in which, given a set of observations, the goal is to recover the track of a muon; and (2) the coincident event problem, which is to determine how many muons are active in the detector during a time window. Rather than solving these problems by developing more complex physical models that are applied at later stages of the analysis, our approach is to augment the detectors early reconstruction with data filters and robust statistical techniques. These can be implemented at the level of on-line reconstruction and, therefore, improve all subsequent reconstructions. Using the metric of median angular resolution, a standard metric for track reconstruction, we improve the accuracy in the initial reconstruction direction by 13%. We also present improvements in measuring the number of muons in coincident events: we can accurately determine the number of muons 98% of the time.

Papers on searches for Dark Matter and exotic particles submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on point source searches submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on ice properties, reconstruction and future developments submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on neutrino oscillation and supernova searches submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on cosmic rays submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

We present the results of a search for neutrino point sources using the IceCube data collected between April 2008 and May 2011 with three partially completed configurations of the detector: the 40-, 59- and 79-string configurations. The live-time of this data set are 1,040 days. An unbinned maximum likelihood ratio test was used to search for an excess of neutrinos above the atmospheric background at any given direction in the sky. By adding two more years of data with improved event selection and reconstruction techniques, the sensitivity was improved by a factor 3.5 or more with respect to the previously published results obtained with the 40-string configuration of IceCube. We performed an all-sky survey and a dedicated search using a catalog of \textit{a priori} selected objects observed by other telescopes. In both searches, the data are compatible with the background-only hypothesis. In the absence of evidence for a signal, we set upper limits on the flux of muon neutrinos. For an E$^{-2}$ neutrino spectrum, the observed limits are between 0.9 and $23.2\times 10^{-12}$ TeV$^{-1}$ cm$^{-2}$s$^{-1}$. We also report upper limits for neutrino emission from groups of sources which were selected according to theoretical models or observational parameters and analysed with a stacking approach.

Composite dark matter candidates, which can arise from new strongly-coupled sectors, are well-motivated and phenomenologically interesting, particularly in the context of asymmetric generation of the relic density. In this work, we employ lattice calculations to study the electromagnetic form factors of electroweak-neutral dark-matter baryons for a three-color, QCD-like theory with Nf = 2 and 6 degenerate fermions in the fundamental representation. We calculate the (connected) charge radius and anomalous magnetic moment, both of which can play a significant role for direct detection of composite dark matter. We find minimal Nf dependence in these quantities. We generate mass-dependent cross-sections for dark matter-nucleon interactions and use them in conjunction with experimental results from XENON100, excluding dark matter candidates of this type with masses below 10 TeV.

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

We present a full analysis of the Probing Evolution And Reionization Spectroscopically (PEARS) slitess grism spectroscopic data obtained with the Advanced Camera for Surveys on HST. PEARS covers fields within both the Great Observatories Origins Deep Survey (GOODS) North and South fields, making it ideal as a random survey of galaxies, as well as the availability of a wide variety of ancillary observations to support the spectroscopic results. Using the PEARS data we are able to identify star forming galaxies within the redshift volume 0< z<1.5. Star forming regions in the PEARS survey are pinpointed independently of the host galaxy. This method allows us to detect the presence of multiple emission line regions (ELRs) within a single galaxy. 1162 Ha, [OIII] and/or [OII] emission lines have been identified in the PEARS sample of ~906 galaxies down to a limiting flux of ~1e-18 erg/s/cm^2. The ELRs have also been compared to the properties of the host galaxy, including morphology, luminosity, and mass. From this analysis we find three key results: 1) The computed line luminosities show evidence of a flattening in the luminosity function with increasing redshift; 2) The star forming systems show evidence of disturbed morphologies, with star formation occurring predominantly within one effective (half-light) radius. However, the morphologies show no correlation with host stellar mass; and 3) The number density of star forming galaxies with M_* > 1e9} M_sun decreases by an order of magnitude at z<0.5 relative to the number at 0.5<z<0.9 in support of the argument for galaxy downsizing.

We report on the observation of two neutrino-induced events which have an estimated deposited energy in the IceCube detector of 1.04 $\pm$ 0.16 and 1.14 $\pm$ 0.17 PeV, respectively, the highest neutrino energies observed so far. These events are consistent with fully contained particle showers induced by neutral-current $\nu_{e,\mu,\tau}$ ($\bar\nu_{e,\mu,\tau}$) or charged-current $\nu_{e}$ ($\bar\nu_{e}$) interactions within the IceCube detector. The events were discovered in a search for ultra-high energy neutrinos using data corresponding to 615.9 days effective livetime. The expected number of atmospheric background is $0.082 \pm 0.004 \text{(stat)}^{+0.041}_{-0.057} \text{(syst)}$. The probability to observe two or more candidate events under the atmospheric background-only hypothesis is $2.9\times10^{-3}$ ($2.8\sigma$) taking into account the uncertainty on the expected number of background events. These two events could be a first indication of an astrophysical neutrino flux, the moderate significance, however, does not permit a definitive conclusion at this time.

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 \AA. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 \AA) that could be the dominant contribution to the total ionizing flux at z >= 6.

We report on the observation of a significant deficit of cosmic rays from the direction of the Moon with the IceCube detector. The study of this "Moon shadow" is used to characterize the angular resolution and absolute pointing capabilities of the detector. The detection is based on data taken in two periods before the completion of the detector: between April 2008 and May 2009, when IceCube operated in a partial configuration with 40 detector strings deployed in the South Pole ice, and between May 2009 and May 2010 when the detector operated with 59 strings. Using two independent analysis methods, the Moon shadow has been observed to high significance (> 6 sigma) in both detector configurations. The observed location of the shadow center is within 0.2 degrees of its expected position when geomagnetic deflection effects are taken into account. This measurement validates the directional reconstruction capabilities of IceCube.

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

Modern radio pulsar surveys produce a large volume of prospective candidates, the majority of which are polluted by human-created radio frequency interference or other forms of noise. Typically, large numbers of candidates need to be visually inspected in order to determine if they are real pulsars. This process can be labor intensive. In this paper, we introduce an algorithm called PEACE (Pulsar Evaluation Algorithm for Candidate Extraction) which improves the efficiency of identifying pulsar signals. The algorithm ranks the candidates based on a score function. Unlike popular machine-learning based algorithms, no prior training data sets are required. This algorithm has been applied to data from several large-scale radio pulsar surveys. Using the human-based ranking results generated by students in the Arecibo Remote Command enter programme, the statistical performance of PEACE was evaluated. It was found that PEACE ranked 68% of the student-identified pulsars within the top 0.17% of sorted candidates, 95% within the top 0.34%, and 100% within the top 3.7%. This clearly demonstrates that PEACE significantly increases the pulsar identification rate by a factor of about 50 to 1000. To date, PEACE has been directly responsible for the discovery of 47 new pulsars, 5 of which are millisecond pulsars that may be useful for pulsar timing based gravitational-wave detection projects.

Gamma-Ray Bursts (GRBs) have been proposed as a leading candidate for acceleration of ultra high-energy cosmic rays, which would be accompanied by emission of TeV neutrinos produced in proton-photon interactions during acceleration in the GRB fireball. Two analyses using data from two years of the IceCube detector produced no evidence for this neutrino emission, placing strong constraints on models of neutrino and cosmic-ray production in these sources.

A search for muon neutrinos from dark matter annihilations in the Galactic Center region has been performed with the 40-string configuration of the IceCube Neutrino Observatory using data collected in 367 days of live-time starting in April 2008. The observed fluxes were consistent with the atmospheric background expectations. Upper limits on the self-annihilation cross-section are obtained for dark matter particle masses ranging from 100 GeV to 10 TeV. In the case of decaying dark matter, lower limits on the lifetime have been determined for masses between 200 GeV and 20 TeV.

We report the first measurement of the atmospheric electron neutrino flux in the energy range between approximately 80 GeV and 6 TeV, using data recorded during the first year of operation of IceCube's DeepCore low energy extension. Techniques to identify neutrinos interacting within the DeepCore volume and veto muons originating outside the detector are demonstrated. A sample of 1029 events is observed in 281 days of data, of which 496 $\pm$ 66(stat.) $\pm$ 88(syst.) are estimated to be cascade events, including both electron neutrino and neutral current events. The rest of the sample includes residual backgrounds due to atmospheric muons and charged current interactions of atmospheric muon neutrinos. The flux of the atmospheric electron neutrinos is consistent with models of atmospheric neutrinos in this energy range. This constitutes the first observation of electron neutrinos and neutral current interactions in a very large volume neutrino telescope optimized for the TeV energy range.

The IceCube Neutrino Observatory, approximately 1 km^3 in size, is now complete with 86 strings deployed in the Antarctic ice. IceCube detects the Cherenkov radiation emitted by charged particles passing through or created in the ice. To realize the full potential of the detector, the properties of light propagation in the ice in and around the detector must be well understood. This report presents a new method of fitting the model of light propagation in the ice to a data set of in-situ light source events collected with IceCube. The resulting set of derived parameters, namely the measured values of scattering and absorption coefficients vs. depth, is presented and a comparison of IceCube data with simulations based on the new model is shown.

We analyze the spectral energy distributions (SEDs) of Lyman break galaxies (LBGs) at z=1-3 selected using the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) UVIS channel filters. These HST/WFC3 observations cover about 50 sq. arcmin in the GOODS-South field as a part of the WFC3 Early Release Science program. These LBGs at z=1-3 are selected using dropout selection criteria similar to high redshift LBGs. The deep multi-band photometry in this field is used to identify best-fit SED models, from which we infer the following results: (1) the photometric redshift estimate of these dropout selected LBGs is accurate to within few percent; (2) the UV spectral slope (beta) is redder than at high redshift (z>3), where LBGs are less dusty; (3) on average, LBGs at z=1-3 are massive, dustier and more highly star-forming, compared to LBGs at higher redshifts with similar luminosities (0.1L*<~L<~2.5L*), though their median values are similar within 1-sigma uncertainties. This could imply that identical dropout selection technique, at all redshifts, find physically similar galaxies; and (4) stellar masses of these LBGs are directly proportional to their UV luminosities with a logarithmic slope of ~0.46, and star-formation rates are proportional to their stellar masses with a logarithmic slope of ~0.90. These relations hold true --- within luminosities probed in this study --- for LBGs from z~1.5 to 5. The star-forming galaxies selected using other color-based techniques show similar correlations at z~2, but to avoid any selection biases, and for direct comparison with LBGs at z>3, a true Lyman break selection at z~2 is essential. The future HST UV surveys, both wider and deeper, covering a large luminosity range are important to better understand LBG properties, and their evolution.

We report on searches for neutrino sources at energies above 200 GeV in the Northern sky of the galactic plane, using the data collected by the South Pole neutrino telescopes IceCube and AMANDA. The galactic region considered here includes the Local Arm towards the Cygnus region and our closest approach to the Perseus Arm. The data have been collected between 2007 and 2009 when AMANDA was an integrated part of IceCube, which was still under construction and operated with 22-strings (2007-8) and 40-strings (2008-9) of optical modules deployed in the ice. By combining the larger IceCube detector with the lower energy threshold of the more compact AMANDA detector, we obtain an improved sensitivity at energies below $\sim$10 TeV with respect to previous searches. The analyses presented here are: a scan for point sources within the galactic plane; a search optimized for multiple and extended sources in the Cygnus region, which might be below the sensitivity of the point source scan; and studies of seven pre-selected neutrino source candidates. For one of them, Cygnus X-3, a time-dependent search for neutrinos in coincidence with observed radio and X-ray flares has been performed. No evidence of a signal is found, and upper limits are reported for each of the searches. We investigate neutrino spectra proportional to E$^{-2}$ and E$^{-3}$ to cover the entire range of possible spectra. The soft E$^{-3}$ spectrum results in an energy distribution similar to a source with cut-off below $\sim$50 TeV. For the considered region of the galactic plane, the 90% confidence level muon neutrino flux upper limits are in the range E$^3$dN/dE$\sim 5.4 - 19.5 \times 10^{-11} \rm{TeV^{2} cm^{-2} s^{-1}}$ for point-like neutrino sources in the energy region [180.0 GeV - 20.5 TeV]. These represent the most stringent upper limits for soft-spectra neutrino sources within the Galaxy reported to date.