To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed. In this study we present the first result of the DESTINYS survey (Disk Evolution Study Through Imaging of Nearby Young Stars). DESTINYS is an ESO/SPHERE large program that aims at studying disk evolution in scattered light, mainly focusing on a sample of low-mass stars (<1$M_\odot$) in nearby (~200 pc) star-forming regions. In this particular study we present the observations of the ET Cha (RECX 15) system, a nearby 'old' classical T Tauri star (5-8 Myr, ~100 pc), which is still strongly accreting. We use SPHERE/IRDIS in H-band polarimetric imaging mode to obtain high contrast images of the ET Cha system to search for scattered light from the circumstellar disk as well as thermal emission from close companions. We additionally employ VLT/NACO total intensity archival data taken in 2003. We report here the discovery of a low-mass (sub)stellar companion with SPHERE/IRDIS to ET Cha. We are estimating the mass of this new companion based on photometry. Depending on the system age it is a 5 Myr, 50 $M_{Jup}$ brown dwarf or an 8 Myr, 0.10 $M_\odot$ M-type pre-main-sequence star. We explore possible orbital solutions and discuss the recent dynamic history of the system. Independent of the precise companion mass we find that the presence of the companion likely explains the small size of the disk around ET Cha. The small separation of the binary pair indicates that the disk around the primary component is likely clearing from the outside in, explaining the high accretion rate of the system.

While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolution of $\sim0.3''$ ($\sim50$ au). The CO isotopologue emission traces four major components of the RU Lup system: a compact Keplerian disk with a radius of $\sim120$ au, a non-Keplerian ``envelope-like'' structure surrounding the disk and extending to $\sim260$ au from the star, at least five blueshifted spiral arms stretching up to 1000 au, and clumps outside the spiral arms located up to 1500 au in projection from RU Lup. We comment on potential explanations for RU Lup's peculiar gas morphology, including gravitational instability, accretion of material onto the disk, or perturbation by another star. RU Lup's extended non-Keplerian CO emission, elevated stellar accretion rate, and unusual photometric variability suggest that it could be a scaled-down Class II analog of the outbursting FU Ori systems.

The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems.

There is mounting evidence that the composition and structure of planetary systems are intimately linked to their birth environments. During the past decade, several spectral surveys probed the chemistry of the earliest stages of star formation and of late planet-forming disks. However, very little is known about the chemistry of intermediate protostellar stages, i.e. Class I Young Stellar Objects (YSOs), where planet formation may have already begun. We present here the first results of a 3mm spectral survey performed with the IRAM-30m telescope to investigate the chemistry of a sample of seven Class I YSOs located in the Taurus star-forming region. These sources were selected to embrace the wide diversity identified for low-mass protostellar envelope and disk systems. We present detections and upper limits of thirteen small ($N_{\rm atoms}\leq3$) C, N, O, and S carriers - namely CO, HCO$^+$, HCN, HNC, CN, N$_2$H$^+$, C$_2$H, CS, SO, HCS$^+$, C$_2$S, SO$_2$, OCS - and some of their D, $^{13}$C, $^{15}$N, $^{18}$O, $^{17}$O, and $^{34}$S isotopologues. Together, these species provide constraints on gas-phase C/N/O ratios, D- and $^{15}$N-fractionation, source temperature and UV exposure, as well as the overall S-chemistry. We find substantial evidence of chemical differentiation among our source sample, some of which can be traced back to Class I physical parameters, such as the disk-to-envelope mass ratio (proxy for Class I evolutionary stage), the source luminosity, and the UV-field strength. Overall, these first results allow us to start investigating the astrochemistry of Class I objects, however, interferometric observations are needed to differentiate envelope versus disk chemistry.

The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.

We report new HI observations of four z$\sim$0.05 star-forming galaxies undertaken during the commissioning phase of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). FAST is the largest single-dish telescope with a 500 meter aperture and a 19-Beam receiver. Exploiting the unprecedented sensitivity provided by FAST, we aim to study the atomic gas, via the HI 21cm emission line, in low-$z$ star-forming galaxies taken from the Valpara\'iso ALMA/APEX Line Emission Survey (VALES) project. Together with previous ALMA CO($J=1-0$) observations, the HI data provides crucial information to measure the gas mass and dynamics. As a pilot HI survey, we targeted four local star-forming galaxies at $z\sim0.05$. In particular, one of them has already been detected in HI by the Arecibo Legacy Fast ALFA survey (ALFALFA), allowing a careful comparison. We use an ON-OFF observing approach that allowed us to reach an rms of 0.7mJy/beam at a 1.7km/s velocity resolution within only 20 minutes ON-target integration time. We demonstrate the great capabilities of the FAST 19-beam receiver for pushing the detectability of the HI emission line of extra-galactic sources. The HI emission line detected by FAST shows good consistency with the previous ALFALFA results. Our observations are put in context with previous multi-wavelength data to reveal the physical properties of these low-$z$ galaxies. We find that the CO($J=1-0$) and HI emission line profiles are similar. The dynamical mass estimated from the HI data is an order of magnitude higher than the baryon mass and the dynamical mass derived from the CO observations, implying that the mass probed by dynamics of HI is dominated by the dark matter halo. In one case, a target shows an excess of CO($J=1-0$) in the line centre, which can be explained by an enhanced CO($J=1-0$) emission induced by a nuclear starburst showing high velocity dispersion.

We present the results of our systematic search for the reddest far-infrared (FIR) and submillimeter (sub-mm) galaxies using the data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and the SCUBA2 Cosmological Legacy Survey (S2CLS). The red FIR galaxies are "500~$\mu$m risers," whose spectral energy distributions increase with wavelength across the three FIR passbands of the Spectral and Photometric Imaging REceiver (SPIRE) of Herschel. Within 106.5 deg$^2$ of the HerMES fields, we have selected 629 500 $\mu$m risers. The red sub-mm galaxies are "SPIRE dropouts," which are prominent detections in the S2CLS 850 $\mu$m data but are extremely weak or invisible in the SPIRE bands. Within the 2.98 deg$^2$ common area of HerMES and S2CLS, we have selected 95 such objects. These very red sources could be dusty starbursts at high redshifts ($z\gtrsim 4$-6) because the peak of their cold-dust emission heated by star formation is shifted to the reddest FIR/sub-mm bands. The surface density of 500 $\mu$m risers is $\sim$8.2 deg$^{-2}$ at the $\geq 20$ mJy level in 500 $\mu$m, while that of SPIRE dropouts is $\sim$19.3 deg$^{-2}$ at the $\geq 5$ mJy level in 850 $\mu$m. Both type of objects could span a wide range of redshifts, however. Using deep radio data in these fields to further select the ones likely at the highest redshifts, we find that the surface density of $z>6$ candidates is 5.5 deg$^{-2}$ among 500 $\mu$m risers and is 0.8-13.6 deg$^{-2}$ among SPIRE dropouts. If this is correct, the dust-embedded star formation processes in such objects could contribute comparably as Lyman-break galaxies to the global SFR density at $z>6$.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $\mu$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $\mu$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2''$ to $1.5''$ ($1'' \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the S\'ersic index $\simeq 0.95$ and the S\'ersic radius of $\sim0.18'' (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $\mu$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius.

The slow solar wind is typically characterized as having low Alfv\'enicity. However, Parker Solar Probe (PSP) observed predominately Alfv\'enic slow solar wind during several of its initial encounters. From its first encounter observations, about 55.3\% of the slow solar wind inside 0.25 au is highly Alfv\'enic ($|\sigma_C| > 0.7$) at current solar minimum, which is much higher than the fraction of quiet-Sun-associated highly Alfv\'enic slow wind observed at solar maximum at 1 au. Intervals of slow solar wind with different Alfv\'enicities seem to show similar plasma characteristics and temperature anisotropy distributions. Some low Alfv\'enicity slow wind intervals even show high temperature anisotropies, because the slow wind may experience perpendicular heating as fast wind does when close to the Sun. This signature is confirmed by Wind spacecraft measurements as we track PSP observations to 1 au. Further, with nearly 15 years of Wind measurements, we find that the distributions of plasma characteristics, temperature anisotropy and helium abundance ratio ($N_\alpha/N_p$) are similar in slow winds with different Alfv\'enicities, but the distributions are different from those in the fast solar wind. Highly Alfv\'enic slow solar wind contains both helium-rich ($N_\alpha/N_p\sim0.045$) and helium-poor ($N_\alpha/N_p\sim0.015$) populations, implying it may originate from multiple source regions. These results suggest that highly Alfv\'enic slow solar wind shares similar temperature anisotropy and helium abundance properties with regular slow solar winds, and they thus should have multiple origins.

Theory and observation for the search for life on exoplanets via atmospheric "biosignature gases" is accelerating, motivated by the capabilities of the next generation of space- and ground-based telescopes. The most observationally accessible rocky planet atmospheres are those dominated by molecular hydrogen gas, because the low density of H$_2$-gas leads to an expansive atmosphere. The capability of life to withstand such exotic environments, however, has not been tested in this context. We demonstrate that single-celled microorganisms ($\textit{E. coli}$ and yeast) that normally do not inhabit H$_2$-dominated environments can survive and grow in a 100% H$_2$ atmosphere. We also describe the astonishing diversity of dozens of different gases produced by $\textit{E. coli}$, including many already proposed as potential biosignature gases (e.g., nitrous oxide, ammonia, methanethiol, dimethylsulfide, carbonyl sulfide, and isoprene). This work demonstrates the utility of lab experiments to better identify which kinds of alien environments can host some form of possibly detectable life.

We investigate whether the hard X-ray photon index (${\Gamma}$) versus accretion rate correlation for super-Eddington accreting quasars is different from that for sub-Eddington accreting quasars. We construct a sample of 113 bright quasars from the Sloan Digital Sky Survey Data Release 14 quasar catalog, including 38 quasars as the super-Eddington subsample and 75 quasars as the sub-Eddington subsample. We derive black-hole masses using a simple-epoch virial mass formula based on the ${\rm H\beta}$ lines, and we use the standard thin disk model to derive the dimensionless accretion rates ($\dot{\mathscr{M}}$) for our sample. The X-ray data for these quasars are collected from the Chandra and XMM-Newton archives. We fit the hard X-ray spectra using a single power-law model to obtain ${\Gamma}$ values. We find a statistically significant ($R_{\rm S}=0.43$, $p=7.75\times{10}^{-3}$) correlation between ${\Gamma}$ and $\dot{\mathscr{M}}$ for the super-Eddington subsample. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlation for the sub-Eddington subsample is also significant, but weaker ($R_{\rm S}=0.30$, $p=9.98\times{10}^{-3}$). Linear regression analysis shows that ${\rm \Gamma}=(0.34\pm0.11){\rm log}{\dot{\mathscr{M}}}+(1.71\pm0.17)$ and ${\rm \Gamma}=(0.09\pm0.04){\rm log}{\dot{\mathscr{M}}}+(1.93\pm0.04)$ for the super- and sub-Eddington subsamples, respectively. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlations of the two subsamples are different, suggesting different disk-corona connections in these two types of systems. We propose one qualitative explanation of the steeper ${\Gamma}$-$\dot{\mathscr{M}}$ correlation in the super-Eddington regime that involves larger seed photon fluxes received by the compact coronae from the thick disks in super-Eddington accreting quasars.

The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^17 eV, complementing the reach of IceCube. An ARIANNA observation of a neutrino source would provide strong insight into the enigmatic sources of cosmic rays. ARIANNA observes the radio emission from high energy neutrino interactions in the Antarctic ice. Among radio based concepts under current investigation, ARIANNA-200 would uniquely survey the vast majority of the southern sky at any instant in time, and an important region of the northern sky, by virtue of its location on the surface of the Ross Ice Shelf in Antarctica. The broad sky coverage is specific to the Moore's Bay site, and makes ARIANNA-200 ideally suited to contribute to the multi-messenger thrust by the US National Science Foundation, Windows on the Universe - Multi-Messenger Astrophysics, providing capabilities to observe explosive sources from unknown directions. The ARIANNA architecture is designed to measure the angular direction to within 3 degrees for every neutrino candidate, which too plays an important role in the pursuit of multi-messenger observations of astrophysical sources.

A comprehensive analysis of plasma and composition characteristics inside magnetic clouds (MCs) observed by the Advanced Composition Explorer (ACE) spacecraft from 1998 February to 2011 August is presented. The results show that MCs have specific interior structures, and MCs of different speeds show differences in composition and structure. Compared with the slow MCs, fast MCs have enhanced mean charge states of iron, oxygen, silicon, magnesium, $\mathrm{O^{7+}/O^{6+}}$, $\mathrm{C^{6+}/C^{5+}}$, $\mathrm{C^{6+}/C^{4+}}$ and $\mathrm{Fe^{\geq16+}/Fe_{total}}$ values. For ionic species in fast MCs, a higher atomic number represents a greater enhancement of mean charge state than slow MCs. We also find that both the fast and slow MCs display bimodal structure distribution in the mean iron charge state ($\mathrm{\langle Q\rangle Fe}$), which suggests that the existence of flux rope prior to the eruption is common. Furthermore, the $\mathrm{\langle Q\rangle Fe} $, $\mathrm{Fe^{\geq16+}/Fe_{total}}$, and $\mathrm{O^{7+}/O^{6+}}$ ratio distribution inside fast MCs have the feature that the posterior peak is higher than the anterior one. This result agrees with the "standard model" for CME/flares, by which magnetic reconnection occurs beneath the flux rope, thereby ionizing the ions of the posterior part of flux rope sufficiently by high-energy electron collisions or by direct heating in the reconnection region.

H$_2$CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an ALMA survey of H$_2$CO towards 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H$_2$CO is detected towards 13 disks and tentatively detected towards a 14th. We find both centrally-peaked and centrally-depressed emission morphologies, and half of the disks show ring-like structures at or beyond expected CO snowline locations. Together these morphologies suggest that H$_2$CO in disks is commonly produced through both gas-phase and CO-ice-regulated grain-surface chemistry. We extract disk-averaged and azimuthally-averaged H$_2$CO excitation temperatures and column densities for four disks with multiple H$_2$CO line detections. The temperatures are between 20-50K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H$_2$CO emission in disks is generally emerging from the warm molecular layer, with some contributions from the colder midplane. Applying the same H$_2$CO excitation temperatures to all disks in the survey, we find that H$_2$CO column densities span almost three orders of magnitude ($\sim 5 \times 10^{11} - 5 \times 10^{14} \mathrm{cm}^{-2}$). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H$_2$CO compared to T Tauri disks, possibly because of less CO freeze-out. More H$_2$CO observations towards Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H$_2$CO and other oxygen-bearing organics efficiently form during planet formation.

In this paper we present our imaging observations on the CO(7-6) line and its underlying continuum emission of the young submillimeter galaxy LESS 073 at redshift 4.755, using the Atacama Large Millimeter/submillimeter Array (ALMA). At the achieved resolution of $\sim$$1^{\prime\prime}.2\times0^{\prime\prime}.9$ ($8\times6$~kpc$^2$), the CO(7-6) emission is largely unresolved (with a deconvolved size of $1^{\prime\prime}.1(\pm0^{\prime\prime}.5) \times 0^{\prime\prime}.9(\pm0^{\prime\prime}.8)$), and the continuum emission is totally unresolved. The CO(7-6) line emission has an integrated flux of $0.86\pm0.08$~Jy km/s, and a line width of $343\pm40$ km/s. The continuum emission has a flux density of 0.51 mJy. By fitting the observed far-infrared (FIR) spectral energy distribution of LESS 073 with a single-temperature modified blackbody function, we obtained a dust temperature $T_{\rm dust}=57.6\pm3.5$ K, 60-to-100 $\mu$m flux density ratio $f_{60}/f_{100}=0.86\pm0.08$, and total infrared luminosity $L_{\rm IR}=(5.8\pm0.9) \times 10^{12}~L_\odot$. The SED-fit-based $f_{60}/f_{100}$ is consistent with those estimated from various line ratios as advocated by our earlier work, indicating that those proposed line-ratio-based method can be used to practically derive $f_{60}/f_{100}$ for high-$z$ sources. The total molecular gas mass of LESS 073 is $(3.3\pm1.7) \times10^{10}~M_\odot$, and the inferred gas depletion time is about 43 Myr.

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $\alpha$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

We determined the rotational periods of 3122 Florence, 3830 Trelleborg, and (131077) 2000 YH105 with the Harvard Clay Telescope and KeplerCam at the Fred L. Whipple Observatory. We found the rotational periods to be 2.3580 $\pm$ 0.0015 h, 17.059 $\pm$ 0.017 h, and 1.813 $\pm$ 0.00003 h, respectively. Our measurement of 3122 Florence's period agrees with Warner (2016), who reported 2.3580 $\pm$ 0.0002 h.

We present evidence for localised deviations from Keplerian rotation, i.e., velocity "kinks", in 8 of 18 circumstellar disks observed by the DSHARP program: DoAr 25, Elias 2-27, GW Lup, HD 143006, HD 163296, IM Lup, Sz 129 and WaOph 6. Most of the kinks are detected over a small range in both radial extent and velocity, suggesting a planetary origin, but for some of them foreground contamination prevents us from measuring their spatial and velocity extent. Because of the DSHARP limited spectral resolution and signal-to-noise in the 12CO J=2-1 line, as well as cloud contamination, the kinks are usually detected in only one spectral channel, and will require confirmation. The strongest circumstantial evidence for protoplanets in the absence of higher spectral resolution data and additional tracers is that, upon deprojection, we find that all of the candidate planets lie within a gap and/or at the end of a spiral detected in dust continuum emission. This suggests that a significant fraction of the dust gaps and spirals observed by ALMA in disks are caused by embedded protoplanets.

According to \emph{Wind} observations between June 2004 and May 2019, this Letter investigates the proton and alpha particle temperatures in the space of ($\theta_d$, $V_d/V_A$) for the first time, where $\theta_d$ and $V_d$ are the radial angle and magnitude of alpha$-$proton differential flow vector respectively, $V_A$ is the local Alfv\'en speed. Results show that the temperatures significantly depend on $\theta_d$ as well as $V_d/V_A$. In case of low proton parallel beta ($\beta_{p{\parallel}} < 1$), it is found that the proton perpendicular temperature is clearly enhanced when $\theta_d$ is small ($\lesssim 45^\circ$) and $V_d/V_A \gtrsim 0.5$. On the contrary, the perpendicular temperature of alpha particles is considerably enhanced when $\theta_d$ is large ($\gtrsim 90^\circ$) or $V_d/V_A$ is sufficiently small. The maximum of proton parallel temperature takes place at $\theta_d \sim 90^\circ$ accompanied by higher $\beta_{p{\parallel}}$ and by larger turbulence amplitude of magnetic fluctuations in inertial range. This study should present strong evidence for cyclotron resonance heating of protons and alpha particles in the solar wind. Other mechanisms including Landau resonance and stochastic heating are also proposed, which tend to have different ($\theta_d$, $V_d/V_A$) spaces than cyclotron resonance heating.

Most of the massive star-forming galaxies are found to have `inside-out' stellar mass growth modes, which means the inner parts of the galaxies mainly consist of the older stellar population, while the star forming in the outskirt of the galaxy is still ongoing. The high-resolution HST images from Hubble Deep UV Legacy Survey (HDUV) and Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) projects with the unprecedented depth in both F275W and F160W bands are the perfect data sets to study the forming and formed stellar distribution directly. We selected the low redshift ($0.05 < z_{\rm spec} < 0.3$) galaxy sample from the GOODS-North field where the HST F275W and F160W images are available. Then we measured the half light radius in F275W and F160W bands, which are the indicators of the star formation and stellar mass. By comparing the F275W and F160W half light radius, we find the massive galaxies are mainly follow the `inside-out' growth which is consistent with the previous results. Moreover, the HST F275W and F160W images reveal that some of the low-mass galaxies ($<10^8M_\odot$) have the `outside-in' growth mode: their images show a compact UV morphology, implying an ongoing star formation in the galaxy centre, the stars in the outskirts of the galaxies are already formed. The two modes transit smoothly at stellar mass range about $10^{8-9}M_\odot$ with a large scatter. We also try to identify the possible neighbour massive galaxies from the SDSS data, which represent the massive galaxy sample. We find that all of the spec-z selected galaxies have no massive galaxy nearby. Thus the `outside-in' mode we find in the low-mass galaxies are not likely originated from the environment.

We study axion strings of hyperlight axions coupled to photons. Hyperlight axions -- axions lighter than Hubble at recombination -- are a generic prediction of the string axiverse. These axions strings produce a distinct quantized polarization rotation of CMB photons which is $\mathcal{O}(\alpha_{\rm em})$. As the CMB light passes many strings, this polarization rotation converts E-modes to B-modes and adds up like a random walk. Using numerical simulations we show that the expected size of the final result is well within the reach of current and future CMB experiments through the measurement of correlations of CMB B-modes with E- and T-modes. The quantized polarization rotation angle is topological in nature and can be seen as a geometric phase. Its value depends only on the anomaly coefficient and is independent of other details such as the axion decay constant. Measurement of the anomaly coefficient by measuring this rotation will provide information about the UV theory, such as the quantization of electric charge and the value of the fundamental unit of charge. The presence of axion strings in the universe relies only on a phase transition in the early universe after inflation, after which the string network rapidly approaches an attractor scaling solution. If there are additional stable topological objects such as domain walls, axions as heavy as $10^{-15}$ eV would be accessible. The existence of these strings could also be probed by measuring the relative polarization rotation angle between different images in gravitationally lensed quasar systems.

Stochastic heating is a non-linear heating mechanism driven by the violation of magnetic moment invariance due to large-amplitude turbulent fluctuations producing diffusion of ions towards higher kinetic energies in the direction perpendicular to the magnetic field. It is frequently invoked as a mechanism responsible for the heating of ions in the solar wind. Here, we quantify for the first time the proton stochastic heating rate $Q_\perp$ at radial distances from the Sun as close as $0.16$ au, using measurements from the first two Parker Solar Probe encounters. Our results for both the amplitude and radial trend of the heating rate, $Q_\perp \propto r^{-2.5}$, agree with previous results based on the Helios data set at heliocentric distances from 0.3 to 0.9 au. Also in agreement with previous results, $Q_\perp$ is significantly larger in the fast solar wind than in the slow solar wind. We identify the tendency in fast solar wind for cuts of the core proton velocity distribution transverse to the magnetic field to exhibit a flat-top shape. The observed distribution agrees with previous theoretical predictions for fast solar wind where stochastic heating is the dominant heating mechanism.

Understanding the physical processes in the solar wind and corona which actively contribute to heating, acceleration, and dissipation is a primary objective of NASA's Parker Solar Probe (PSP) mission. Observations of coherent electromagnetic waves at ion scales suggests that linear cyclotron resonance and non-linear processes are dynamically relevant in the inner heliosphere. A wavelet-based statistical study of coherent waves in the first perihelion encounter of PSP demonstrates the presence of transverse electromagnetic waves at ion resonant scales which are observed in 30-50\% of radial field intervals. Average wave amplitudes of approximately 4 nT are measured, while the mean duration of wave events is of order 20 seconds; however long duration wave events can exist without interruption on hour-long timescales. Though ion scale waves are preferentially observed during intervals with a radial mean magnetic field, we show that measurement constraints, associated with single spacecraft sampling of quasi-parallel waves superposed with anisotropic turbulence, render the measured quasi-parallel ion-wave spectrum unobservable when the mean magnetic field is oblique to the solar wind flow; these results imply that the occurrence of coherent ion-scale waves is not limited to a radial field configuration. The lack of strong radial scaling of characteristic wave amplitudes and duration suggests that the waves are generated {\em{in-situ}} through plasma instabilities. Additionally, observations of proton distribution functions indicate that temperature anisotropy may drive the observed ion-scale waves.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

We report the detection of spiral substructure in both the gas velocity and temperature structure of the disk around TW~Hya, suggestive of planet-disk interactions with an unseen planet. Perturbations from Keplerian rotation tracing out a spiral pattern are observed in the SE of the disk, while significant azimuthal perturbations in the gas temperature are seen in the outer disk, outside 90~au, extending the full azimuth of the disk. The deviation in velocity is either $\Delta v_{\phi} \, / \, v_{\rm kep} \sim 0.1$ or $\Delta v_{z} \, / \, v_{\rm kep} \sim 0.01$ depending on whether the perturbation is in the rotational or vertical direction, while radial perturbations can be ruled out. Deviations in the gas temperature are $\pm 4$ K about the azimuthally averaged profile, equivalent to deviations of $\Delta T_{\rm gas} \, / \, T_{\rm gas} \sim 0.05$. Assuming all three structures can be described by an Archimedean spiral, measurements of the pitch angles of both velocity and temperature spirals show a radially decreasing trend for all three, ranging from 9$^{\circ}$ at 70 au, dropping to 3$^{\circ}$ at 200 au. Such low pitch-angled spirals are not readily explained through the wake of an embedded planet in the location of previously reported at 94 au, but rather require a launching mechanism which results in much more tightly wound spirals. Molecular emission tracing distinct heights in the disk is required to accurately distinguish between spiral launching mechanisms.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

The ubiquitous solar jets or jet-like activities are generally regarded as an important source of energy and mass input to the upper solar atmosphere and the solar wind. However, questions about their triggering and driving mechanisms are not completely understood. By taking advantage of high temporal and high spatial resolution stereoscopic observations taken by the Solar Dynamic Observatory (SDO) and the Solar Terrestrial Relations Observatory (STEREO), we report an intriguing two-sided-loop jet occurred on 2013 June 02, which was dynamically associated with the eruption of a mini-filament below an overlying large filament, and two distinct reconnection processes are identified during the formation stage. The SDO observations reveals that the two-sided-loop jet showed a concave shape with a projection speed of about 80 - 136. From the other view angle, the STEREO ahead observations clearly showed that the trajectory of the two arms of the two-sided-loop were along the cavity magnetic field lines hosting the large filament. Contrary to the well-accepted theoretical model, the present observation sheds new light on our understanding of the formation mechanism of two-sided-loop jets. Moreover, the eruption of the two-sided-loop jet not only supplied mass to the overlying large filament, but also provided a rare opportunity to diagnose the magnetic structure of the overlying large filament via the method of three-dimensional reconstruction.

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.

We show that the leading coupling between a shift symmetric inflaton and the Standard Model fermions leads to an induced electroweak symmetry breaking due to particle production during inflation, and as a result, a unique oscillating feature in non-Gaussianities. In this one parameter model, the enhanced production of Standard Model fermions dynamically generates a new electroweak symmetry breaking minimum, where the Higgs field classically rolls into. The production of fermions stops when the Higgs expectation value and hence the fermion masses become too large, suppressing fermion production. The balance between the above-mentioned effects gives the Standard Model fermions masses that are uniquely determined by their couplings to the inflaton. In particular, the heaviest Standard Model fermion, the top quark, can produce a distinct cosmological collider physics signature characterised by a one-to-one relation between amplitude and frequency of the oscillating signal, which is observable at future 21-cm surveys.

The detection of an oscillating pattern in the bispectrum of density perturbations could suggest the existence of a high-energy second minimum in the Higgs potential. If the Higgs field resided in this new minimum during inflation and was brought back to the electroweak vacuum by thermal corrections during reheating, the coupling of Standard Model particles to the inflaton would leave its imprint on the bispectrum. We focus on the fermions, whose dispersion relation can be modified by the coupling to the inflaton, leading to an enhanced particle production during inflation even if their mass during inflation is larger than the Hubble scale. This results in a large non-analytic contribution to non-Gaussianities, with an amplitude $f_{\rm NL}$ as large as $100$ in the squeezed limit, potentially detectable in future 21-cm surveys. Measuring the contributions from two fermions would allow us to compute the ratio of their masses, and to ascribe the origin of the signal to a new Higgs minimum. Such a discovery would be a tremendous step towards understanding the vacuum instability of the Higgs potential, and could have fascinating implications for anthropic considerations.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.

The extreme infrared (IR) luminosity of local luminous and ultra-luminous IR galaxies (U/LIRGs; 11 < log LIR /Lsun < 12 and log LIR /Lsun > 12, respectively) is mainly powered by star-formation processes triggered by mergers or interactions. While U/LIRGs are rare locally, at z > 1, they become more common, they dominate the star-formation rate (SFR) density, and a fraction of them are found to be normal disk galaxies. Therefore, there must be an evolution of the mechanism triggering these intense starbursts with redshift. To investigate this evolution, we present new optical SWIFT integral field spectroscopic H{\alpha}+[NII] observations of a sample of 9 intermediate-z (0.2 < z < 0.4) U/LIRG systems selected from Herschel 250{\mu}m observations. The main results are the following: (a) the ratios between the velocity dispersion and the rotation curve amplitude indicate that 10-25% (1-2 out of 8) might be compatible with being isolated disks while the remaining objects are interacting/merging systems; (b) the ratio between un-obscured and obscured SFR traced by H{\alpha} and LIR, respectively, is similar in both local and these intermediate-z U/LIRGs; and (c) the ratio between 250{\mu}m and the total IR luminosities of these intermediate-z U/LIRGs is higher than that of local U/LIRGs with the same LIR . This indicates a reduced dust temperature in these intermediate-z U/LIRGs. This, together with their already measured enhanced molecular gas content, suggests that the interstellar medium conditions are different in our sample of intermediate-z galaxies when compared to local U/LIRGs.

We present observations of the C-band $1_{10}-1_{11}$ (4.8 GHz) and Ku-band $2_{11}-2_{12}$ (14.5 GHz) K-doublet lines of H$_2$CO and the C-band $1_{10}-1_{11}$ (4.6 GHz) line of H$_2$$^{13}$CO toward a large sample of Galactic molecular clouds, through the Shanghai Tianma 65-m radio telescope (TMRT). Our sample with 112 sources includes strong H$_2$CO sources from the TMRT molecular line survey at C-band and other known H$_2$CO sources. All three lines are detected toward 38 objects (43 radial velocity components) yielding a detection rate of 34\%. Complementary observations of their continuum emission at both C- and Ku-bands were performed. Combining spectral line parameters and continuum data, we calculate the column densities, the optical depths and the isotope ratio H$_2$$^{12}$CO/H$_2$$^{13}$CO for each source. To evaluate photon trapping caused by sometimes significant opacities in the main isotopologue's rotational mm-wave lines connecting our measured K-doublets, and to obtain $^{12}$C/$^{13}$C abundance ratios, we used the RADEX non-LTE model accounting for radiative transfer effects. This implied the use of the new collision rates from \citet{Wiesenfeld2013}. Also implementing distance values from trigonometric parallax measurements for our sources, we obtain a linear fit of $^{12}$C/$^{13}$C = (5.08$\pm$1.10)D$_{GC}$ + (11.86$\pm$6.60), with a correlation coefficient of 0.58. D$_{GC}$ refers to Galactocentric distances. Our $^{12}$C/$^{13}$C ratios agree very well with the ones deduced from CN and C$^{18}$O but are lower than those previously reported on the basis of H$_2$CO, tending to suggest that the bulk of the H$_2$CO in our sources was formed on dust grain mantles and not in the gas phase.

We discuss models in which vacua other than our own can be directly observed in the present universe. Models with density-dependent vacuum structure can give rise to `non-lethal'-vacua: vacua with lower energy-density than our vacuum, but only in regions with finite Standard Model densities. These models provide an explicit example of a bubble which is confined to a finite region of space and produces potentially detectable signatures, unlike standard Coleman tunneling events where bubbles expand at the speed of light and are never directly observable. We study the expansion and contraction of a confined bubble created after a core-collapse supernova, focusing on energy deposition that may be observable in the vicinity of a supernova remnant due to the formation and evolution of a confined bubble.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

The gas pixel detector (GPD) is designed and developed for high-sensitivity astronomical X-ray polarimetry, which is a new window about to open in a few years. Due to the small mass, low power, and compact geometry of the GPD, we propose a CubeSat mission Polarimeter Light (PolarLight) to demonstrate and test the technology directly in space. There is no optics but a collimator to constrain the field of view to 2.3 degrees. Filled with pure dimethyl ether (DME) at 0.8 atm and sealed by a beryllium window of 100 micron thick, with a sensitive area of about 1.4 mm by 1.4 mm, PolarLight allows us to observe the brightest X-ray sources on the sky, with a count rate of, e.g., ~0.2 counts/s from the Crab nebula. The PolarLight is 1U in size and mounted in a 6U CubeSat, which was launched into a low Earth Sun-synchronous orbit on October 29, 2018, and is currently under test. More launches with improved designs are planned in 2019. These tests will help increase the technology readiness for future missions such as the enhanced X-ray Timing and Polarimetry (eXTP), better understand the orbital background, and may help constrain the physics with observations of the brightest objects.

The observation of gravitational waves from a binary neutron star merger by LIGO/VIRGO and the associated electromagnetic counterpart provides a high precision test of orbital dynamics, and therefore a new and sensitive probe of extra forces and new radiative degrees of freedom. Axions are one particularly well-motivated class of extensions to the Standard Model leading to new forces and sources of radiation, which we focus on in this paper. Using an effective field theory (EFT) approach, we calculate the first post-Newtonian corrections to the orbital dynamics, radiated power, and gravitational waveform for binary neutron star mergers in the presence of an axion. This result is applicable to many theories which add an extra massive scalar degree of freedom to General Relativity. We then perform a detailed forecast of the potential for Advanced LIGO to constrain the free parameters of the EFT, and map these to the mass $m_a$ and decay constant $f_a$ of the axion. At design sensitivity, we find that Advanced LIGO can potentially exclude axions with $m_a \lesssim 10^{-11} \ {\rm eV}$ and $f_a \sim (10^{14} - 10^{17}) \ {\rm GeV}$. There are a variety of complementary observational probes over this region of parameter space, including the orbital decay of binary pulsars, black hole superradiance, and laboratory searches. We comment on the synergies between these various observables.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics established by $^{12}$CO $J=2-1$ observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at $R\sim69$ au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large ($R>20$ au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks.

Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

To characterize the substructures induced in protoplanetary disks by the interaction between stars in multiple systems, we study the $1.25\,$mm continuum and the $^{12}$CO$(J=2-1)$ spectral line emission of the triple systems HT Lup and AS 205, at scales of $\approx 5\,$au, as part of the 'Disk Substructures at High Angular Resolution Project' (DSHARP). In the continuum emission, we find two symmetric spiral arms in the disk around AS 205 N, with pitch angle of $14^\circ$, while the southern component AS 205 S, itself a spectroscopic binary, is surrounded by a compact inner disk and a bright ring at a radius of $34\,$au. The $^{12}$CO line exhibits clear signatures of tidal interactions, with spiral arms, extended arc-like emission, and high velocity gas, possible evidence of a recent close encounter between the disks in the AS 205 system, as these features are predicted by hydrodynamic simulations of fly-by encounters. In the HT Lup system, we detect continuum emission from all three components. The primary disk, HT Lup A, also shows two-armed symmetric spiral structure with a pitch angle of $4^\circ$, while HT Lup B and C, located at $25$ and $434\,$au in projected separation from HT Lup A, are barely resolved with $\sim5$ and $\sim10\,$au in diameter, respectively. The gas kinematics for the closest pair indicates a different sense of rotation for each disk, which could be explained by either a counter rotation of the two disks in different, close to parallel, planes, or by a projection effect of these disks with a close to $90^\circ$ misalignment between them.

We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting the two innermost rings. Although our $^{12}$CO observations at similar spatial resolution do not show obvious substructure, they reveal an inner disk depleted of CO emission. From the continuum emission and the CO velocity field we find that the innermost ring has a higher inclination than the outermost rings and the arc. This is evidence for either a small ($\sim8^{\circ}$) or moderate ($\sim41^{\circ}$) misalignment between the inner and outer disk, depending on the specific orientation of the near/far sides of the inner/outer disk. We compare the observed substructures in the ALMA observations with recent scattered light data from VLT/SPHERE of this object. In particular, the location of narrow shadow lanes in the SPHERE image combined with pressure scale height estimates, favor a large misalignment of about $41^{\circ}$. We discuss our findings in the context of a dust-trapping vortex, planet-carved gaps, and a misaligned inner disk due to the presence of an inclined companion to HD 143006.

The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_\nu$ to dust surface density $\Sigma_d$.

A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low alpha_turb << 5e-4 and very large grains a_grain >> 0.1 cm can be excluded by the data for all the rings studied in this paper.

The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques.

We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained with concentric Gaussian rings. In particular, the dust emission consists of a series of narrow and closely spaced rings in the inner $\sim$60 au, two well-separated bright rings in the outer disk, centered at 74 and 120 au, and at least two fainter emission features at 90 and 130 au. We model the visibilities with a parametric representation of the radial surface brightness profile, consisting of a central core and 7 concentric Gaussian rings. Recent hydro-dynamical simulations of low viscosity disks show that super-Earth planets can produce the multiple gaps seen in AS 209 millimeter continuum emission. The $^{12}$CO line emission is centrally peaked and extends out to $\sim$300 au, much farther than the millimeter dust emission. We find axisymmetric, localized deficits of CO emission around four distinct radii, near 45, 75, 120 and 210 au. The outermost gap is located well beyond the edge of the millimeter dust emission, and therefore cannot be due to dust opacity and must be caused by a genuine CO surface density reduction, due either to chemical effects or depletion of the overall gas content.

ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and $^{12}$CO J=2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 67 and 100 au. They also reveal new small scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission inform about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate to the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet unseen planets.

We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 au FWHM) observations of their 240 GHz (1.25 mm) continuum emission. These data provide a first homogeneous look at the small-scale features in disks that are directly relevant to the planet formation process, quantifying their prevalence, morphologies, spatial scales, spacings, symmetry, and amplitudes, for targets with a variety of disk and stellar host properties. We find that these substructures are ubiquitous in this sample of large, bright disks. They are most frequently manifested as concentric, narrow emission rings and depleted gaps, although large-scale spiral patterns and small arc-shaped azimuthal asymmetries are also present in some cases. These substructures are found at a wide range of disk radii (from a few au to more than 100 au), are usually compact ($<$10 au), and show a wide range of amplitudes (brightness contrasts). Here we discuss the motivation for the project, describe the survey design and the sample properties, detail the observations and data calibration, highlight some basic results, and provide a general overview of the key conclusions that are presented in more detail in a series of accompanying articles. The DSHARP data -- including visibilities, images, calibration scripts, and more -- are released for community use at https://almascience.org/alma-data/lp/DSHARP.

The HIgh-Resolution Mid-infrarEd Spectrometer (HIRMES) is the 3rd Generation Instrument for the Stratospheric Observatory For Infrared Astronomy (SOFIA), currently in development at the NASA Goddard Space Flight Center (GSFC), and due for commissioning in 2019. By combining direct-detection Transition Edge Sensor (TES) bolometer arrays, grating-dispersive spectroscopy, and a host of Fabry-Perot tunable filters, HIRMES will provide the ability for High Resolution (R~100,000), Mid-Resolution (R~10,000), and Low-Resolution (R~600) slit-spectroscopy, and 2D Spectral Imaging (R~2000 at selected wavelengths) over the 25 - 122 {\mu}m mid-far infrared waveband. The driving science application is the evolution of proto-planetary systems via measurements of water-vapor, water-ice, deuterated hydrogen (HD), and neutral oxygen lines. However, HIRMES has been designed to be as flexible as possible to cover a wide range of science cases that fall within its phase-space, all whilst reaching sensitivities and observing powers not yet seen thus far on SOFIA, providing unique observing capabilities which will remain unmatched for decades.

Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.

[Abridged] We combine new CO(1-0) line observations of 24 intermediate redshift galaxies (0.03 < z < 0.28) along with literature data of galaxies at 0<z<4 to explore scaling relations between the dust and gas content using PAH 6.2 $\mu$m ($L_{6.2}$), CO ($L'_{\rm CO}$), and infrared ($L_{\rm IR}$) luminosities for a wide range of redshifts and physical environments. Our analysis confirms the existence of a universal $L_{6.2}-L_{\rm CO}$ correlation followed by normal star-forming galaxies (SFGs) and starbursts (SBs) at all redshifts. This relation is also followed by local ULIRGs that appear as outliers in the $L_{6.2}-L_{\rm IR}$ and $L_{\rm IR}-L'_{\rm CO}$ relations from the sequence defined by normal SFGs. The emerging tight ($\sigma \approx 0.26$ dex) and linear ($\alpha = 1.03$) relation between $L_{6.2}$ and $L_{\rm CO}$ indicates a $L_{6.2}$ to molecular gas ($M_{\rm H_2}$) conversion factor of $\alpha_{6.2} = M_{\rm H2}/L_{6.2} = (2.7\pm1.3) \times \alpha_{\rm CO}$, where $\alpha_{\rm CO}$ is the $L'_{\rm CO}$ to $M_{\rm H_2}$ conversion factor. We also find that on galaxy integrated scales, PAH emission is better correlated with cold rather than with warm dust emission, suggesting that PAHs are associated with the diffuse cold dust, which is another proxy for $M_{\rm H_2}$. Focusing on normal SFGs among our sample, we employ the dust continuum emission to derive $M_{\rm H_2}$ estimates and find a constant $M_{\rm H_2}/L_{6.2}$ ratio of $\alpha_{6.2} = 12.3 \ M_{\rm H_2}/{\rm L}_{\odot}$ ($\sigma\approx 0.3$ dex). We propose that the presented $L_{6.2}-L'_{\rm CO}$ and $L_{6.2}-M_{\rm H_2}$ relations will serve as useful tools for the determination of the physical properties of high-$z$ SFGs, for which PAH emission will be routinely detected by the James Webb Space Telescope.

We present new constraints on gas-phase C, N, and O abundances in the molecular layer of the IM Lup protoplanetary disk. Building on previous physical and chemical modeling of this disk, we use new ALMA observations of C$_2$H to constrain the C/O ratio in the molecular layer to be $\sim0.8$, i.e., higher than the solar value of $\sim0.54$. We use archival ALMA observations of HCN and H$^{13}$CN to show that no depletion of N is required (assuming an interstellar abundance of $7.5\times10^{-5}$ per H). These results suggest that an appreciable fraction of O is sequestered in water ice in large grains settled to the disk mid-plane. Similarly, a fraction of the available C is locked up in less volatile molecules. By contrast, N remains largely unprocessed, likely as N$_2$. This pattern of depletion suggests the presence of true abundance variations in this disk, and not a simple overall depletion of gas mass. If these results hold more generally, then combined CO, C$_2$H, and HCN observations of disks may provide a promising path for constraining gas-phase C/O and N/O during planet-formation. Together, these tracers offer the opportunity to link the volatile compositions of disks to the atmospheres of planets formed from them.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

Ly$\alpha$ Emitters (LAEs) may represent an important galaxy population in the low mass regime. We present our deep narrowband imaging surveys in the COSMOS and ECDF-S fields and study the properties of LAEs at $z=2.23\pm0.03$. The narrowband surveys conducted at Magellan II telescope allow us to obtain a sample of 452 LAEs reaching a $5\sigma$ limiting magnitude of $\sim26$ mag. Our Ly$\alpha$ luminosity functions extend to $10^{41.8}$ erg s$^{-1}$ with steep faint-end slope. Using multi-wavelength ancillary data, especially the deep Spitzer/IRAC 3.6$\mu$m and 4.5$\mu$m photometric data, we obtained reliable stellar mass estimates for 130 IRAC-detected LAEs, spanning a range of $8 < {\rm log}(M_\star/M_\odot)< 11.5$. For the remaining IRAC-undetected LAEs, the median-stacked spectral energy distribution yields a stellar mass of ${\rm log}(M_\star/M_\odot)=7.97^{+0.05}_{-0.07}$ and the rest-frame ultraviolet emission indicates a median star formation rate of ${\rm log} (SFR/M_\odot$ yr$^{-1})=-0.14\pm0.35$. There are six LAEs detected by the Spitzer/MIPS 24$\mu$m or even Herschel far-infrared observations. Taking into account the six MIR/FIR detected LAEs, our LAEs cover a wide range in the star formation rate (${\rm 1<SFR<2000}$ M$_\odot$ yr$^{-1}$). Although LAEs as a population are diverse in their stellar properties, they are mostly low-mass star-forming galaxies and follow the star formation main sequence relations or their extrapolations to the low-mass end, implying a normal star-forming nature of LAEs. The clustering analysis indicates that our LAEs reside in dark matter halos with ${\rm <\log(M_{h}/M_{\odot})> =10.8^{+0.56}_{-1.1}}$, suggesting that they are progenitors of local Large Magellanic Cloud-like galaxies.

This paper describes new deep 3.6 and 4.5 micron imaging of three UltraVISTA near-infrared survey stripes within the COSMOS field. The observations were carried out with Spitzer's Infrared Array Camera (IRAC) for the Spitzer Matching Survey of the Ultra-VISTA Deep Stripes (SMUVS). In this work we present our data reduction techniques, and document the resulting mosaics, coverage maps, and catalogs in both IRAC passbands for the three easternmost UltraVISTA survey stripes, covering a combined area of about 0.66 square degrees, of which 0.45 square degrees have at least 20 hr integration time. SMUVS reaches point-source sensitivities of about 25.0 AB mag at both 3.6 and 4.5 microns with a significance of 4-sigma accounting for both survey sensitivity and source confusion. To this limit the SMUVS catalogs contain a total of about 350,000 sources, each of which is detected significantly in at least one IRAC band. Because of its uniform and high sensitivity, relatively large area coverage, and the wide array of ancillary data available in COSMOS, the SMUVS survey will be useful for a large number of cosmological investigations. We will make all images and catalogues described herein publicly available via the Spitzer Science Center.

We present the results from our Atacama Large Millimeter/submillimeter Array (ALMA) imaging observations of the CO (7-6), [CI] 370 um (hereafter [CI]) and [NII] 205 um (hereafter [NI]I) lines and their underlying continuum emission of BRI 1335-0417, an infrared bright quasar at z = 4.407. At the achieved resolutions of 1.1" to 1.2" (or 7.5 to 8.2 kpc), the continuum at 205 and 372 um (rest-frame), the CO (7-6), and the [CI] emissions are at best barely resolved whereas the [NII] emission is well resolved with an ALMA beam de-convolved major axis of 1.3" (+/- 0.3") or 9 (+/-2) kpc. As a warm dense gas tracer, the CO (7-6) emission shows a more compact spatial distribution and a significantly higher peak velocity dispersion than the other two lines that probe lower density gas, a picture favoring a merger-triggered star formation (SF) scenario over an orderly rotating SF disk. The CO (7-6) data also indicate a possible QSO-driven gas outflow that reaches a maximum line-of-sight velocity of 500 to 600 km/s. The far-infrared (FIR) dust temperature (T_dust) of 41.5 K from a gray-body fit to the continuum agrees well with the average T_dust inferred from various line luminosity ratios. The resulting L_CO(7-6)/L_FIR luminosity ratio is consistent with that of local luminous infrared galaxies powered predominantly by SF. The CO(7-6) luminosity-inferred SF rate is 5.1 (+/-1.5) x 10^3 M_solar/yr . The system has an effective star-forming region of 1.7 (+1.7/-0.8) kpc in diameter and a molecular gas reservoir of ~5 x 10^{11} M_solar.

We analyze sparsely-sampled near-infrared (JHKs) light curves of a sample of 1781 Mira variable candidates in M33, originally discovered using I-band time-series observations. We extend our single-band semi-parametric Gaussian process modeling of Mira light curves to a multi-band version and obtain improved period determinations. We use our previous results on near-infrared properties of candidate Miras in the LMC to classify the majority of the M33 sample into Oxygen- or Carbon-rich subsets. We derive Period-Luminosity relations for O-rich Miras and determine a distance modulus for M33 of 24.80 +- 0.06 mag.

The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky.

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

Modern radio interferometers enable observations of spectral lines with unprecedented spatial resolution and sensitivity. In spite of these technical advances, many lines of interest are still at best weakly detected and therefore necessitate detection and analysis techniques specialized for the low signal-to-noise ratio (SNR) regime. Matched filters can leverage knowledge of the source structure and kinematics to increase sensitivity of spectral line observations. Application of the filter in the native Fourier domain improves SNR while simultaneously avoiding the computational cost and ambiguities associated with imaging, making matched filtering a fast and robust method for weak spectral line detection. We demonstrate how an approximate matched filter can be constructed from a previously observed line or from a model of the source, and we show how this filter can be used to robustly infer a detection significance for weak spectral lines. When applied to ALMA Cycle 2 observations of CH3OH in the protoplanetary disk around TW Hya, the technique yields a ~53% SNR boost over aperture-based spectral extraction methods, and we show that an even higher boost will be achieved for observations at higher spatial resolution. A Python-based open-source implementation of this technique is available under the MIT license at https://github.com/AstroChem/VISIBLE.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

With growing data from ongoing and future supernova surveys it is possible to empirically quantify the shapes of SNIa light curves in more detail, and to quantitatively relate the shape parameters with the intrinsic properties of SNIa. Building such relationship is critical in controlling systematic errors associated with supernova cosmology. Based on a collection of well-observed SNIa samples accumulated in the past years, we construct an empirical SNIa light curve model using a statistical method called the functional principal component analysis (FPCA) for sparse and irregularly sampled functional data. Using this method, the entire light curve of an SNIa is represented by a linear combination of principal component functions, and the SNIa is represented by a few numbers called principal component scores. These scores are used to establish relations between light curve shapes and physical quantities such as intrinsic color, interstellar dust reddening, spectral line strength, and spectral classes. These relations allow for descriptions of some critical physical quantities based purely on light curve shape parameters. Our study shows that some important spectral feature information is being encoded in the broad band light curves, for instance, we find that the light curve shapes are correlated with the velocity and velocity gradient of the Si II $\lambda$6355 line. This is important for supernova surveys, e.g., LSST and WFIRST. Moreover, the FPCA light curve model is used to construct the entire light curve shape, which in turn is used in a functional linear form to adjust intrinsic luminosity when fitting distance models.

Solar radio type III bursts are believed to be the most sensitive signature of near-relativistic electron beam propagation in the corona. A solar radio type IIIb-III pair burst with fine frequency structures, observed by the Low Frequency Array (LOFAR) with high temporal ($\sim10$ ms) and spectral (12.5 kHz) resolutions at 30 - 80 MHz, is presented. The observations show that the type III burst consists of many striae, which have a frequency scale of about 0.1 MHz in both the fundamental (plasma) and the harmonic (double plasma) emission. We investigate the effects of background density fluctuations based on the observation of striae structure to estimate the density perturbation in solar corona. It is found that the spectral index of the density fluctuation spectrum is about $-1.7$, and the characteristic spatial scale of the density perturbation is around $700$ km. This spectral index is very close to a Kolmogorov turbulence spectral index of $-5/3$, consistent with a turbulent cascade. This fact indicates that the coronal turbulence may play the important role of modulating the time structures of solar radio type III bursts, and the fine structure of radio type III bursts could provide a useful and unique tool to diagnose the turbulence in the solar corona.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We analyze high angular resolution ALMA observations of the TW Hya disk to place constraints on the CO and dust properties. We present new, sensitive observations of the $^{12}$CO $J = 3-2$ line at a spatial resolution of 8 AU (0\farcs14). The CO emission exhibits a bright inner core, a shoulder at $r\approx70$ AU, and a prominent break in slope at $r\approx90$ AU. Radiative transfer modeling is used to demonstrate that the emission morphology can be reasonably reproduced with a $^{12}$CO column density profile featuring a steep decrease at $r\approx15$ AU and a secondary bump peaking at $r\approx70$ AU. Similar features have been identified in observations of rarer CO isotopologues, which trace heights closer to the midplane. Substructure in the underlying gas distribution or radially varying CO depletion that affects much of the disk's vertical extent may explain the shared emission features of the main CO isotopologues. We also combine archival 1.3 mm and 870 $\mu$m continuum observations to produce a spectral index map at a spatial resolution of 2 AU. The spectral index rises sharply at the continuum emission gaps at radii of 25, 41, and 47 AU. This behavior suggests that the grains within the gaps are no larger than a few millimeters. Outside the continuum gaps, the low spectral index values of $\alpha\approx 2$ indicate either that grains up to centimeter size are present, or that the bright continuum rings are marginally optically thick at millimeter wavelengths.

In this work we present new APEX/SEPIA Band-5 observations targeting the CO ($J=2\text{-}1$) emission line of 24 Herschel-detected galaxies at $z=0.1-0.2$. Combining this sample {with} our recent new Valpara\'iso ALMA Line Emission Survey (VALES), we investigate the star formation efficiencies (SFEs = SFR/$M_{\rm H_{2}}$) of galaxies at low redshift. We find the SFE of our sample bridges the gap between normal star-forming galaxies and Ultra-Luminous Infrared Galaxies (ULIRGs), which are thought to be triggered by different star formation modes. Considering the $\rm SFE'$ as the SFR and the $L'_{\rm CO}$ ratio, our data show a continuous and smooth increment as a function of infrared luminosity (or star formation rate) with a scatter about 0.5 dex, instead of a steep jump with a bimodal behaviour. This result is due to the use of a sample with a much larger range of sSFR/sSFR$_{\rm ms}$ using LIRGs, with luminosities covering the range between normal and ULIRGs. We conclude that the main parameters controlling the scatter of the SFE in star-forming galaxies are the systematic uncertainty of the $\alpha_{\rm CO}$ conversion factor, the gas fraction and physical size.

The recent direct observation of gravitational waves (GW) from merging black holes opens up the possibility of exploring the theory of gravity in the strong regime at an unprecedented level. It is therefore interesting to explore which extensions to General Relativity (GR) could be detected. We construct an Effective Field Theory (EFT) satisfying the following requirements. It is testable with GW observations; it is consistent with other experiments, including short distance tests of GR; it agrees with widely accepted principles of physics, such as locality, causality and unitarity; and it does not involve new light degrees of freedom. The most general theory satisfying these requirements corresponds to adding to the GR Lagrangian operators constructed out of powers of the Riemann tensor, suppressed by a scale comparable to the curvature of the observed merging binaries. The presence of these operators modifies the gravitational potential between the compact objects, as well as their effective mass and current quadrupoles, ultimately correcting the waveform of the emitted GW.

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

Previous studies revealed that electromagnetic cyclotron waves (ECWs) near the proton cyclotron frequency exist widely in the solar wind, and the majority of ECWs are left-handed (LH) polarized waves. Using the magnetic field data from the STEREO mission, this Letter carries out a survey of ECWs over a long period of 7 years, and calculates the occurrence rates of ECWs with different polarization senses. Results show that the occurrence rate is nearly a constant for the ECWs with right-handed polarization, but it varies significantly for the ECWs with LH polarization. Further investigation of plasma conditions reveals that the LH ECWs take place preferentially in a plasma characterized by higher temperature, lower density, and larger velocity. Some considerable correlations between the occurrence rate of LH ECWs and the properties of ambient plasmas are discussed. The present research may provide evidence for effect of alpha particles on generation of ECWs.

We study the near-infrared properties of 690 Mira candidates in the central region of the Large Magellanic Cloud, based on time-series observations at JHKs. We use densely-sampled I-band observations from the OGLE project to generate template light curves in the near infrared and derive robust mean magnitudes at those wavelengths. We obtain near-infrared Period-Luminosity relations for Oxygen-rich Miras with a scatter as low as 0.12 mag at Ks. We study the Period-Luminosity-Color relations and the color excesses of Carbon-rich Miras, which show evidence for a substantially different reddening law.

We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field <U> with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years.

Milli-magnetically charged particles generically appear in scenarios with kinetic mixing. We present model independent bounds on these particles coming from magnetars. Schwinger pair production discharges the magnetic field of the magnetar. Thus the existence of large magnetic fields at magnetars place strong bounds on the milli-magnetic charge to be smaller than $10^{-18}$ over a large mass range.

H2CO is one of the most readily detected organic molecules in protoplanetary disks. Yet its distribution and dominant formation pathway(s) remain largely unconstrained. To address these issues, we present ALMA observations of two H2CO lines (3_{12}-2_{11} and 5_{15}-4_{14}) at 0".5 (~30 au) spatial resolution toward the disk around the nearby T Tauri star TW Hya. Emission from both lines is spatially resolved, showing a central depression, a peak at 0".4 radius, and a radial decline at larger radii with a bump at ~1", near the millimeter continuum edge. We adopt a physical model for the disk and use toy models to explore the radial and vertical H2CO abundance structure. We find that the observed emission implies the presence of at least two distinct H2CO gas reservoirs: (1) a warm and unresolved inner component (<10 au), and (2) an outer component that extends from ~15 au to beyond the millimeter continuum edge. The outer component is further constrained by the line ratio to arise in a more elevated disk layer at larger radii. The inferred H2CO abundance structure agrees well with disk chemistry models, which predict efficient H2CO gas-phase formation close to the star, and cold H2CO grain surface formation, through H additions to condensed CO, followed by non-thermal desorption in the outer disk. The implied presence of active grain surface chemistry in the TW Hya disk is consistent with the recent detection of CH3OH emission, and suggests that more complex organic molecules are formed in disks, as well.

We present the discovery of 1847 Mira candidates in the Local Group galaxy M33 using a novel semi-parametric periodogram technique coupled with a Random Forest classifier. The algorithms were applied to ~2.4x10^5 I-band light curves previously obtained by the M33 Synoptic Stellar Survey. We derive preliminary Period-Luminosity relations at optical, near- & mid-infrared wavelengths and compare them to the corresponding relations in the Large Magellanic Cloud.

The deuterium enrichment of molecules is sensitive to their formation environment. Constraining patterns of deuterium chemistry in protoplanetary disks is therefore useful for probing how material is inherited or reprocessed throughout the stages of star and planet formation. We present ALMA observations at $\sim0.6"$ resolution of DCO$^+$, H$^{13}$CO$^+$, DCN, and H$^{13}$CN in the full disks around T Tauri stars AS 209 and IM Lup, the transition disks around T Tauri stars V4046 Sgr and LkCa 15, and the full disks around Herbig Ae stars MWC 480 and HD 163296. We also present ALMA observations of HCN in the IM Lup disk. DCN, DCO$^+$, and H$^{13}$CO$^+$ are detected in all disks, and H$^{13}$CN in all but the IM Lup disk. We find efficient deuterium fractionation for the sample, with estimates of disk-averaged DCO$^+$/HCO$^+$ and DCN/HCN abundance ratios ranging from $\sim0.02-0.06$ and $\sim0.005-0.08$, respectively, which is comparable to values reported for other ISM environments. The relative distributions of DCN and DCO$^+$ vary between disks, suggesting that multiple formation pathways may be needed to explain the diverse emission morphologies. In addition, gaps and rings observed in both H$^{13}$CO$^+$ and DCO$^+$ emission provide new evidence that DCO$^+$ bears a complex relationship with the location of the midplane CO snowline.

We explore the dependence of the incidence of moderate-luminosity ($L_{X} = 10^{41.9-43.7}$ erg s$^{-1}$) AGNs and the distribution of their accretion rates on host color at 0.5 < z < 2.5, using deep X-ray data in GOODS fields. We use extinction-corrected rest-frame U-V colors to divide both AGN hosts and non-AGN galaxies into red sequence (quiescent), green valley (transition), and blue cloud (star-forming) populations. We find that both the AGN fraction at fixed stellar mass and its evolution with redshift are dependent on host colors. Most notably, red galaxies have the lowest AGN fraction (~5\%) at z~1 yet with most rapid redshift evolution, increasing by a factor of 5 (~24\%) at z~2. Green galaxies exhibit the highest AGN fraction across all redshifts, which is most pronounced at z~2 with more than half of them hosting an AGN at $M_{*} > 10^{10.6} M_{\odot}$. Together with the high AGN fraction in red galaxies at z~2, this indicates that X-ray AGNs could be important in both transforming blue galaxies into red ones and subsequently maintaining their quiescence at high redshift. Furthermore, consistent with low-redshift studies, we find that the probability of hosting an AGN in the total galaxy population can be characterized by a universal Eddington ratio ($p(\lambda_{Edd}) \sim \lambda_{Edd}^{-0.4}$) and a moderate redshift evolution. Yet consistent with their different AGN fractions, different populations appear to also have different $p(\lambda_{Edd})$ with red galaxies exhibiting more rapid redshift evolution than green and blue ones. Evidence for a steeper power law of $p(\lambda_{Edd})$ in red galaxies is also presented, though larger samples are needed to confirm. These results suggest that the AGN accretion or the growth of supermassive black holes is related to their host properties, and may also influence their hosts in a different mode dependent on the host color.

Nitrogen fractionation is commonly used to assess the thermal history of Solar System volatiles. With ALMA it is for the first time possible to directly measure 14N/15N ratios in common molecules during the assembly of planetary systems. We present ALMA observations of the H13CN and HC15N J=3-2 lines at 0".5 angular resolution, toward a sample of six protoplanetary disks, selected to span a range of stellar and disk structure properties. Adopting a typical 12C/13C ratio of 70, we find comet-like 14N/15N ratios of 80-160 in 5/6 of the disks (3 T Tauri and 2 Herbig Ae disks) and lack constraints for one of the T Tauri disks (IM Lup). There are no systematic differences between T Tauri and Herbig Ae disks, or between full and transition disks within the sample. In addition, no correlation is observed between disk-averaged D/H and 14N/15N ratios in the sample. One of the disks, V4046 Sgr, presents unusually bright HCN isotopologue emission, enabling us to model the radial profiles of H13CN and HC15N. We find tentative evidence of an increasing 14N/15N ratio with radius, indicating that selective photodissociation in the inner disk is important in setting the 14N/15N ratio during planet formation.

The Ly$\alpha$ escape fraction is a key measure to constrain the neutral state of the intergalactic medium and then to understand how the universe was fully reionized. We combine deep narrowband imaging data from the custom-made filter NB393 and the $H_{2}S$1 filter centered at 2.14 $\mu$m to examine the Ly$\alpha$ emitters and H$\alpha$ emitters at the same redshift $z=2.24$. The combination of these two populations allows us to determine the Ly$\alpha$ escape fraction at $z=2.24$. Over an area of 383 arcmin$^{2}$ in the Extended Chandra Deep Field South (ECDFS), 124 Ly$\alpha$ emitters are detected down to NB393 = 26.4 mag at the 5$\sigma$ level, and 56 H$\alpha$ emitters come from An14. Of these, four have both Ly$\alpha$ and H$\alpha$ emissions (LAHAEs). We measure the individual/volumetric Ly$\alpha$ escape fraction by comparing the observed Ly$\alpha$ luminosity/luminosity density to the extinction-corrected H$\alpha$ luminosity/luminosity density. We revisit the extinction correction for H$\alpha$ emitters using the Galactic extinction law with the color excess for nebular emission. We also adopt the Calzetti extinction law together with an identical color excess for stellar and nebular regions to explore how the uncertainties in extinction correction affect our results. In both cases, an anti-correlation between the Ly$\alpha$ escape fraction and dust attenuation is found among the LAHAEs, suggesting that dust absorption is responsible for the suppression of the escaping Ly$\alpha$ photons. However, the estimated Ly$\alpha$ escape fraction of individual LAHAEs varies up to ~3 percentage points between the two methods of extinction correction. We find the global Ly$\alpha$ escape fraction at $z=2.24$ to be ($3.7\pm1.4$)% in the ECDFS. The variation in the color excess of the extinction causes a discrepancy of ~1 percentage point in the global Ly$\alpha$ escape fraction.

We develop a non-linear semi-parametric Gaussian process model to estimate periods of Miras with sparsely-sampled light curves. The model uses a sinusoidal basis for the periodic variation and a Gaussian process for the stochastic changes. We use maximum likelihood to estimate the period and the parameters of the Gaussian process, while integrating out the effects of other nuisance parameters in the model with respect to a suitable prior distribution obtained from earlier studies. Since the likelihood is highly multimodal for period, we implement a hybrid method that applies the quasi-Newton algorithm for Gaussian process parameters and search the period/frequency parameter over a dense grid. A large-scale, high-fidelity simulation is conducted to mimic the sampling quality of Mira light curves obtained by the M33 Synoptic Stellar Survey. The simulated data set is publicly available and can serve as a testbed for future evaluation of different period estimation methods. The semi-parametric model outperforms an existing algorithm on this simulated test data set as measured by period recovery rate and quality of the resulting Period-Luminosity relations.

Solar flares are the most powerful explosions occurring in the solar system, which may lead to disastrous space weather events and impact various aspects of our Earth. So far, it is still a big challenge in modern astrophysics to understand the origin of solar flares and predict their onset. Based on the analysis of soft X-ray emission observed by the Geostationary Operational Environmental Satellite (GOES), this work reported a new discovery of very long-periodic pulsations occurred in the preflare phase before the onset of solar flares (preflare-VLPs). These pulsations are typically with period of 8 - 30 min and last for about 1 - 2 hours. They are possibly generated from LRC oscillations of plasma loops where electric current dominates the physical process during magnetic energy accumulation in the source region. The preflare-VLP provides an essential information for understanding the triggering mechanism and origin of solar flares, and may help us to response to solar explosions and the corresponding disastrous space weather events as a convenient precursory indicator.

The spatial distribution of gas and solids in protoplanetary disks determines the composition and formation efficiency of planetary systems. A number of disks show starkly different distributions for the gas and small grains compared to millimeter-centimeter sized dust. We present new Atacama Large Millimeter/Submillimeter Array (ALMA) observations of the dust continuum, CO, $^{13}$CO, and C$^{18}$O in the IM Lup protoplanetary disk, one of the first systems where this dust-gas dichotomy was clearly seen. The $^{12}$CO is detected out to a radius of 970 AU, while the millimeter continuum emission is truncated at just 313 AU. Based upon this data, we have built a comprehensive physical and chemical model for the disk structure, which takes into account the complex, coupled nature of the gas and dust and the interplay between the local and external environment. We constrain the distributions of gas and dust, the gas temperatures, the CO abundances, the CO optical depths, and the incident external radiation field. We find that the reduction/removal of dust from the outer disk exposes this region to higher stellar and external radiation and decreases the rate of freeze-out, allowing CO to remain in the gas out to large radial distances. We estimate a gas-phase CO abundance of $5\%$ of the ISM value and a low external radiation field ($G_0\lesssim4$). The latter is consistent with that expected from the local stellar population. We additionally find tentative evidence for ring-like continuum substructure, suggestions of isotope-selective photodissociation, and a diffuse gas halo.

We give the three-dimensional dynamical autonomous systems for most of the popular scalar field dark energy models including (phantom) quintessence, (phantom) tachyon, k-essence and general non-canonical scalar field models, change the dynamical variables from variables $(x, y, \lambda)$ to observable related variables $(w_{\phi}, \Omega_{\phi}, \lambda)$, and show the intimate relationships between those scalar fields that the three-dimensional system of k-essence can reduce to (phantom) tachyon, general non-canonical scalar field can reduce to (phantom) quintessence and k-essence can also reduce to (phantom) quintessence for some special cases. For the applications of the three-dimensional dynamical systems, we investigate several special cases and give the exactly dynamical solutions in detail. In the end of this paper, we argue that, it is more convenient and also has more physical meaning to express the differential equations of dynamical systems in $(w_{\phi}, \Omega_{\phi}, \lambda)$ instead of variables $(x, y, \lambda)$ and to investigate the dynamical system in 3-Dimension instead of 2-Dimension. We also raise a question about the possibility of the chaotic behavior in the spatially flat single scalar field FRW cosmological models in the presence of ordinary matter.

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models.

This work reports a peculiar and interesting train of microwave type III pair bursts in the impulsive rising phase of a solar flare on 2011 September 26. The observations include radio spectrometers at frequency of 0.80 - 2.00 GHz, hard X-ray (RHESSI and FERMI), EUV images of SWAP/PROBA-2 and magnetogram of HMI/SDO. By using a recently developed method (Tan et al. 2016a), we diagnosed the plasma density, temperature, plasma beta, magnetic field near the source region, the energy of energetic electrons and the distance between the acceleration region and the emission start sites of type III bursts. From the diagnostics, we find that: (1) The plasma density, temperature, magnetic field, and the distance between the acceleration region and the emission start sites almost have no obvious variations during the period of type III pair trains, while the energy of electrons has an obvious peak value which is consistent to the hard X-ray emission. (2) The plasma beta is much higher than an unity showing a highly dynamic process near the emission start site of type III bursts. (3) Although the reversed-slope type III branches drift slower at one order of magnitude than that of the normal type III branches, the related downgoing and upgoing electrons still could have same order of magnitude of energy. These facts indicate that both of the upgoing and downgoing electrons are possibly accelerated by similar mechanism and in a small source region. This diagnostics can help us to understand the microphysics in the source region of solar bursts.

Millimeter observations of CO isotopologues are often used to make inferences about protoplanetary disk gas density and temperature structures. The accuracy of these estimates depends on our understanding of CO freezeout and desorption from dust grains. Most models of these processes indicate that CO column density decreases monotonically with distance from the central star due to a decrease in gas density and freezeout beyond the CO snowline. We present ALMA Cycle 2 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=2-1$ emission that instead suggest CO enhancement in the outer disk of T Tauri star AS 209. Most notably, the C$^{18}$O emission consists of a central peak and a ring at a radius of $\sim1''$ (120 AU), well outside the expected CO snowline. We propose that the ring arises from the onset of CO desorption near the edge of the millimeter dust disk. CO desorption exterior to a CO snowline may occur via non-thermal processes involving cosmic rays or high-energy photons, or via a radial thermal inversion arising from dust migration.

Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector.

We report on the identification of the old stellar population galaxy candidates at z > 5. We developed a new infrared color selection scheme to isolate galaxies with the strong Balmer breaks at z > 5, and applied it to the ultra-deep and wide infrared survey data from the Spitzer Extended Deep Survey (SEDS) and the UKIRT Infrared Deep Sky Survey. The eight objects satisfying K - [3.6] > 1.3 and K - [3.6] > 2.4 ([3.6] - [4.5]) + 0.6 are selected in the 0.34 deg^2 SEDS Ultra Deep Survey field. Rich multi-wavelength imaging data from optical to far-infrared are also used to reject blending sources and strong nebular line emitters, and we finally obtained the three most likely evolved galaxies at z > 5. Their stacked spectral energy distribution is fitted well with the old stellar population template with M_{*} = (7.5+-1.5) x 10^{10} Msun, star formation rate = 0.9 +- 0.2 Msun yr^{-1}, dust A_V < 1, and age = 0.7+-0.4 Gyr at z = 5.7+-0.6, where the dusty star-forming galaxies at z ~ 2.8 are disfavored because of the faintness in the 24um. The stellar mass density of these evolved galaxy candidates, (6+-4) x 10^4 Msun Mpc^{-3}, is much lower than that of star-forming galaxies, but the non-zero fraction suggests that initial star-formation and quenching have been completed by z ~ 6.

We present a novel automatic adaptive aperture photometry algorithm for measuring the total magnitudes of merging galaxies with irregular shapes. First, we use a morphological pattern recognition routine for identifying the shape of an irregular source in a background-subtracted image. Then, we extend the shape of the source by using the Dilation image operation to obtain an aperture that is quasi-homomorphic to the shape of the irregular source. The magnitude measured from the homomorphic aperture would thus have minimal contamination from the nearby background. As a test of our algorithm, we applied our technique to the merging galaxies observed by the Sloan Digital Sky Survey (SDSS) and the Canada-France-Hawaii Telescope (CFHT). Our results suggest that the adaptive homomorphic aperture algorithm can be very useful for investigating extended sources with irregular shapes and sources in crowded regions.

The diffuse interstellar band (DIB) at 5797.1 \AA\ is simulated based on three premises: (1) The carrier of the DIB is polar as concluded by T. Oka et al. from the anomalous spectrum toward Herschel 36. (2) The sharp central feature observed by P. J. Sarre's group is the $Q$-branch of a parallel band of a prolate top. (3) The radiative temperature of the environment is $T_r$=2.73 K. A $^2 \Pi \leftarrow$ $^2\Pi$ transition of a linear radical is simulated. Results depend on 10 parameters, with the rotational constant $B$ being the most critical. Comparisons of calculated spectra with observed data constrain $B$, which in turn constrains the number of heavy atoms to 5 $\leq n \leq 7$. Upper limits based on the Her 36 spectrum and lower limits based on stability against photodissociation are also discussed. The latter is based on the assumption that the DIB molecules are produced top-down from the breakdown of dust rather than bottom-up by chemical reactions. A difficulty with this limit is that J. P. Maier's laboratory has observed many molecules within this size range, some of which have been tested astronomically. Candidates are discussed in light of many prolate tops observed by the P. Thaddeus and M. C. McCarthy microwave group.

A new air-shower core-detector array (YAC : Yangbajing Air-shower Core-detector array) has been developed to measure the primary cosmic-ray composition at the "knee" energies in Tibet, China, focusing mainly on the light components. The prototype experiment (YAC-I) consisting of 16 detectors has been constructed and operated at Yangbajing (4300 m a.s.l.) in Tibet since May 2009. YAC-I is installed in the Tibet-III AS array and operates together. In this paper, we performed a Monte Carlo simulation to check the sensitivity of YAC-I+Tibet-III array to the cosmic-ray light component of cosmic rays around the knee energies, taking account of the observation conditions of actual YAC-I+Tibet-III array. The selection of light component from others was made by use of an artificial neural network (ANN). The simulation shows that the light-component spectrum estimated by our methods can well reproduce the input ones within 10\% error, and there will be about 30\% systematic errors mostly induced by the primary and interaction models used. It is found that the full-scale YAC and the Tibet-III array is powerful to study the cosmic-ray composition, in particular, to obtain the energy spectra of protons and helium nuclei around the knee energies.

We present Herschel PACS $\&$ SPIRE far-infrared (FIR) and sub-mm imaging observations for a large K-band selected sample of 88 close major-merger pairs of galaxies (H-KPAIRs) in 6 photometric bands (70, 100, 160, 250, 350, and 500 $\mu m$). Among 132 spiral galaxies in the 44 spiral-spiral (S$+$S) pairs and 44 spiral-elliptical (S$+$E) pairs, 113 are detected in at least one Herschel band. Star formation rate (SFR) and dust mass ($M_{\rm dust}$) are derived from the IR SED fitting. Mass of total gas ($M_{\rm gas}$) is estimated by assuming a constant dust-to-gas mass ratio of 0.01. Star forming spiral galaxies (SFGs) in S$+$S pairs show significant enhancements in both specific star formation rate (sSFR) and star formation efficiency (SFE), while having nearly the same gas mass, compared to control galaxies. On the other hand, for SFGs in S$+$E pairs, there is no significant sSFR enhancement and the mean SFE enhancement is significantly lower than that of SFGs in S$+$S pairs. This suggests an important role for the disc-disc collision in the interaction induced star formation. The $M_{\rm gas}$ of SFGs in S$+$E pairs is marginally lower than that of their counterparts in both S$+$S pairs and the control sample. Paired galaxies with and without interaction signs do not differ significantly in their mean sSFR and SFE. As found in previous works, this much larger sample confirms the primary and secondary spirals in S+S pairs follow a Holmberg effect correlation on sSFR.

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

Boulders are ubiquitously found on the surfaces of small rocky bodies in the inner solar system and their spatial and size distributions give insight into the geological evolution and collisional history of the parent bodies. Using images acquired by the Chang'e-2 spacecraft, more than 200 boulders have been identified over the imaged area of the near-Earth asteroid Toutatis. The cumulative boulder size frequency distribution (SFD) shows a steep slope of -4.4 $\pm$ 0.1, which is indicative of a high degree of fragmentation. Similar to Itokawa, Toutatis probably has a rubble-pile structure, as most boulders on its surface cannot solely be explained by impact cratering. The significantly steeper slope for Toutatis' boulder SFD compared to Itokawa may imply a different preservation state or diverse formation scenarios. In addition, the cumulative crater SFD has been used to estimate a surface crater retention age of approximately 1.6 $\pm$ 0.3 Gyr.

Recently, AMS-02 reported their observed results of cosmic rays(CRs). In addition to the AMS-02 data, we add HESS data to estimate the spectra of CR electrons and the diffuse gamma rays above TeV. In the conventional diffusion model, a global analysis is performed on the spectral features of CR electrons and the diffuse gamma rays by GALRPOP package. The results show that the spectrum structure of the primary component of CR electrons can not be fully reproduced by a simple power law and the relevant break is around hundred GeV. At 99\% C.L., the injection indices above the break decrease from 2.54 to 2.35, but the ones below the break are only in the range 2.746 - 2.751. The spectrum of CR electrons does not need to add TeV cutoff to match the features of HESS data too. Based on the difference between the fluxes of CR electrons and the primary component of them, the predicted excess of CR positrons is consistent with the interpretations as pulsar or dark matter. In the analysis of the Galactic diffuse gamma rays with the indirect constraint of AMS-02 and HESS data, it is found that the fluxes of Galactic diffuse gamma rays are consistent with GeV data of Fermi-LAT in the high latitude regions. The results indicate that the inverse Compton scattering(IC) is the dominant component in the range of the hundred GeV to tens of TeV respectively from the high latitude regions to the low ones, and in the all regions of Galaxy the flux of diffuse gamma rays is less than CR electrons at the energy scale of 20 TeV.

We have discovered an optically rich galaxy cluster at z=1.7089 with star formation occurring in close proximity to the central galaxy. The system, SpARCS104922.6+564032.5, was detected within the Spitzer Adaptation of the red-sequence Cluster Survey, (SpARCS), and confirmed through Keck-MOSFIRE spectroscopy. The rest-frame optical richness of Ngal(500kpc) = 30+/-8 implies a total halo mass, within 500kpc, of ~3.8+/-1.2 x 10^14 Msun, comparable to other clusters at or above this redshift. There is a wealth of ancillary data available, including Canada-France-Hawaii Telescope optical, UKIRT-K, Spitzer-IRAC/MIPS, and Herschel-SPIRE. This work adds submillimeter imaging with the SCUBA2 camera on the James Clerk Maxwell Telescope and near-infrared imaging with the Hubble Space Telescope (HST). The mid/far-infrared (M/FIR) data detect an Ultra-luminous Infrared Galaxy spatially coincident with the central galaxy, with LIR = 6.2+/-0.9 x 10^12 Lsun. The detection of polycyclic aromatic hydrocarbons (PAHs) at z=1.7 in a Spitzer-IRS spectrum of the source implies the FIR luminosity is dominated by star formation (an Active Galactic Nucleus contribution of 20%) with a rate of ~860+/-30 Msun/yr. The optical source corresponding to the IR emission is likely a chain of of > 10 individual clumps arranged as "beads on a string" over a linear scale of 66 kpc. Its morphology and proximity to the Brightest Cluster Galaxy imply a gas-rich interaction at the center of the cluster triggered the star formation. This system indicates that wet mergers may be an important process in forming the stellar mass of BCGs at early times.

Observations of deuterium fractionation in the solar system, and in interstellar and circumstellar material, are commonly used to constrain the formation environment of volatiles. Toward protoplanetary disks, this approach has been limited by the small number of detected deuterated molecules, i.e. DCO$^+$ and DCN. Based on ALMA Cycle 2 observations toward the disk around the T Tauri star AS 209, we report the first detection of N$_2$D$^+$ (J=3-2) in a protoplanetary disk. These data are used together with previous Submillimeter Array observations of N$_2$H$^+$ (J=3-2) to estimate a disk-averaged D/H ratio of 0.3--0.5, an order of magnitude higher than disk-averaged ratios previously derived for DCN/HCN and DCO$^+$/HCO$^+$ around other young stars. The high fractionation in N$_2$H$^+$ is consistent with model predictions. The presence of abundant N$_2$D$^+$ toward AS 209 also suggests that N$_2$D$^+$ and the N$_2$D$^+$/N$_2$H$^+$ ratio can be developed into effective probes of deuterium chemistry, kinematics, and ionization processes outside the CO snowline of disks.

The Far-InfraRed Spectroscopic Explorer (FIRSPEX) is a candidate mission in response to a bi-lateral Small-mission call issued by the European Space Agency (ESA) and the Chinese Academy of Sciences (CAS). FIRSPEX is a small satellite (~1m telescope) operating from Low Earth Orbit (LEO). It consists of a number of heterodyne detection bands targeting key molecular and atomic transitions in the terahertz (THz) and Supra-Terahertz (>1 THz) frequency range. The FIRSPEX bands are: [CII] 158 microns (1.9 THz), [NII] 205 microns (1.46 THz), [CI] 370 microns (0.89 THz), CO(6-5) 433 microns (0.69 THz). The primary goal of FIRSPEX is to perform an unbiased all sky spectroscopic survey in four far-infrared lines delivering the first 3D-maps (high spectral resolution) of the Galaxy. The spectroscopic surveys will build on the heritage of Herschel and complement the broad-band all-sky surveys carried out by the IRAS and AKARI observatories. In addition FIRSPEX will enable targeted observations of nearby and distant galaxies allowing for an in-depth study of the ISM components.

Aiming at the observation of cosmic-ray chemical composition at the "knee" energy region, we have been developinga new type air-shower core detector (YAC, Yangbajing Air shower Core detector array) to be set up at Yangbajing (90.522$^\circ$ E, 30.102$^\circ$ N, 4300 m above sea level, atmospheric depth: 606 g/m$^2$) in Tibet, China. YAC works together with the Tibet air-shower array (Tibet-III) and an underground water cherenkov muon detector array (MD) as a hybrid experiment. Each YAC detector unit consists of lead plates of 3.5 cm thick and a scintillation counter which detects the burst size induced by high energy particles in the air-shower cores. The burst size can be measured from 1 MIP (Minimum Ionization Particle) to $10^{6}$ MIPs. The first phase of this experiment, named "YAC-I", consists of 16 YAC detectors each having the size 40 cm $\times$ 50 cm and distributing in a grid with an effective area of 10 m$^{2}$. YAC-I is used to check hadronic interaction models. The second phase of the experiment, called "YAC-II", consists of 124 YAC detectors with coverage about 500 m$^2$. The inner 100 detectors of 80 cm $\times $ 50 cm each are deployed in a 10 $\times$ 10 matrix from with a 1.9 m separation and the outer 24 detectors of 100 cm $\times$ 50 cm each are distributed around them to reject non-core events whose shower cores are far from the YAC-II array. YAC-II is used to study the primary cosmic-ray composition, in particular, to obtain the energy spectra of proton, helium and iron nuclei between 5$\times$$10^{13}$ eV and $10^{16}$ eV covering the "knee" and also being connected with direct observations at energies around 100 TeV. We present the design and performance of YAC-II in this paper.

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.

In this work, we investigate the rotational dynamics of the ginger-shaped near-Earth asteroid 4179 Toutatis, which was closely observed by Chang'e-2 at a distance of $770\pm120~$ meters from the asteroid's surface during the outbound flyby \citep{Huang2013} on 13 December 2012. A sequence of high-resolution images was acquired during the flyby mission. In combination with ground-based radar observations collected over the last two decades, we analyze these flyby images and determine the orientation of the asteroid at the flyby epoch. The 3-1-3 Euler angles of the conversion matrix from the J2000 ecliptic coordinate system to the body-fixed frame are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$, respectively. The least-squares method is utilized to determine the rotational parameters and spin state of Toutatis. The characteristics of the spin-state parameters and angular momentum variations are extensively studied using numerical simulations, which confirm those reported by \citet{Takahashi2013}. The large amplitude of Toutatis' precession is assumed to be responsible for its tumbling attitude as observed from Earth. Toutatis' angular momentum orientation is determined to be described by $\lambda_{H}=180.2^{+0.2^\circ}_{-0.3^\circ}$ and $\beta_{H}=-54.75^{+0.15^\circ}_{-0.10^\circ}$, implying that it has remained nearly unchanged for two decades. Furthermore, using Fourier analysis to explore the change in the orientation of Toutatis' axes, we reveal that the two rotational periods are 5.38 and 7.40 days, respectively, consistent with the results of the former investigation. Hence, our investigation provides a clear understanding of the state of the rotational dynamics of Toutatis.

We present the results of deep \chandra\ imaging of the central region of the Extended Groth Strip, the AEGIS-X Deep (AEGIS-XD) survey. When combined with previous \chandra\ observations of a wider area of the strip, AEGIS-X Wide (AEGIS-XW; Laird et~al. 2009), these provide data to a nominal exposure depth of 800ks in the three central ACIS-I fields, a region of approximately $0.29$~deg$^{2}$. This is currently the third deepest X-ray survey in existence, a factor $\sim 2-3$ shallower than the Chandra Deep Fields (CDFs) but over an area $\sim 3$ times greater than each CDF. We present a catalogue of 937 point sources detected in the deep \chandra\ observations. We present identifications of our X-ray sources from deep ground-based, Spitzer, GALEX and HST imaging. Using a likelihood ratio analysis, we associate multi band counterparts for 929/937 of our X-ray sources, with an estimated 95~\% reliability, making the identification completeness approximately 94~\% in a statistical sense. Reliable spectroscopic redshifts for 353 of our X-ray sources are provided predominantly from Keck (DEEP2/3) and MMT Hectospec, so the current spectroscopic completeness is $\sim 38$~per cent. For the remainder of the X-ray sources, we compute photometric redshifts based on multi-band photometry in up to 35 bands from the UV to mid-IR. Particular attention is given to the fact that the vast majority the X-ray sources are AGN and require hybrid templates. Our photometric redshifts have mean accuracy of $\sigma=0.04$ and an outlier fraction of approximately 5\%, reaching $\sigma=0.03$ with less than 4\% outliers in the area covered by CANDELS . The X-ray, multi-wavelength photometry and redshift catalogues are made publicly available.

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.

A good observation of preflare activities is important for us to understand the origin and triggering mechanism of solar flares, and to predict the occurrence of solar flares. This work presents the characteristics of microwave spectral fine structures as preflare activities of four solar flares observed by Ond\v{r}ejov radio spectrograph in the frequency range of 0.8--2.0 GHz. We found that these microwave bursts which occurred 1--4 minutes before the onset of flares have spectral fine structures with relatively weak intensities and very short timescales. They include microwave quasi-periodic pulsations (QPP) with very short period of 0.1-0.3 s and dot bursts with millisecond timescales and narrow frequency bandwidths. Accompanying these microwave bursts, there are filament motions, plasma ejection or loop brightening on the EUV imaging observations and non-thermal hard X-ray emission enhancements observed by RHESSI. These facts may reveal certain independent non-thermal energy releasing processes and particle acceleration before the onset of solar flares. They may be conducive to understand the nature of solar flares and predict their occurrence.

We present Herschel far-IR photometry and spectroscopy as well as ground based CO observations of an intermediate redshift (0.21 < z < 0.88) sample of Herschel-selected (ultra)-luminous infrared galaxies (L_IR > 10^11.5L_sun). With these measurements we trace the dust continuum, far-IR atomic line emission, in particular [CII]\,157.7microns, as well as the molecular gas of z~0.3 (U)LIRGs and perform a detailed investigation of the interstellar medium of the population. We find that the majority of Herschel-selected intermediate redshift (U)LIRGs have L_CII/L_FIR ratios that are a factor of about 10 higher than that of local ULIRGs and comparable to that of local normal and high-$z$ star forming galaxies. Using our sample to bridge local and high-z [CII] observations, we find that the majority of galaxies at all redshifts and all luminosities follow a L_CII-L_FIR relation with a slope of unity, from which local ULIRGs and high-z AGN dominated sources are clear outliers. We also confirm that the strong anti-correlation between the L_CII/L_FIR ratio and the far-IR color L_60/L_100 observed in the local Universe holds over a broad range of redshifts and luminosities, in the sense that warmer sources exhibit lower L_CII/L_FIR at any epoch. Intermediate redshift ULIRGs are also characterised by large molecular gas reservoirs and by lower star formation efficiencies compared to that of local ULIRGs. The high L_CII/L_FIR ratios, the moderate star formation efficiencies (L_LIR/L_CO or L_IR/M_gas) and the relatively low dust temperatures of our sample (which are also common characteristics of high-z star forming galaxies with ULIRG-like luminosities) indicate that the evolution of the physical properties of (U)LIRGs between the present day and z > 1 is already significant by z ~ 0.3.

Microwave emission with spectral zebra pattern structures (ZPs) is observed frequently in solar flares and the Crab pulsar. The previous observations show that ZP is only a structure overlapped on the underlying broadband continuum with slight increments and decrements. This work reports an extremely unusual strong ZP burst occurring just at the beginning of a solar flare observed simultaneously by two radio telescopes located in China and Czech Republic and by the extreme ultraviolet (EUV) telescope on board NASA's satellite Solar Dynamics Observatory on 2013 April 11. It is a very short and super strong explosion whose intensity exceeds several times that of the underlying flaring broadband continuum emission, lasting for just 18 s. EUV images show that the flare starts from several small flare bursting points (FBPs). There is a sudden EUV flash with extra enhancement in one of these FBPs during the ZP burst. Analysis indicates that the ZP burst accompanying EUV flash is an unusual explosion revealing a strong coherent process with rapid particle acceleration, violent energy release, and fast plasma heating simultaneously in a small region with short duration just at the beginning of the flare.

A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $\alpha_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.).

We provide a catalog of 391 mid-infrared-selected (MIR, 24$\mu$m) broad-emission-line (BEL, type 1) quasars in the 22 deg$^2$ SWIRE Lockman Hole field. This quasar sample is selected in the MIR from Spitzer MIPS with $S_{\rm 24} > 400\mu$Jy, jointly with an optical magnitude limit of r (AB) $<$ 22.5 for broad line identification. The catalog is based on MMT and SDSS spectroscopy to select BEL quasars, extends the SDSS coverage to fainter magnitudes and lower redshifts, and recovers a more complete quasar population. The MIR-selected quasar sample peaks at $z\sim$1.4, and recovers a significant and constant (20\%) fraction of extended objects with SDSS photometry across magnitudes, which was not included in the SDSS quasar survey dominated by point sources. This sample also recovers a significant population of $z < 3$ quasars at $i > 19.1$. We then investigate the continuum luminosity and line profiles of these MIR quasars, and estimate their virial black hole masses and the Eddington ratios. The SMBH mass shows evidence of downsizing, though the Eddington ratios remain constant at $1 < z < 4$. Compared to point sources in the same redshift range, extended sources at $z < 1$ show systematically lower Eddington ratios. The catalog and spectra are publicly available online.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

We give a full investigation on the dynamics of power-law kinetic quintessence $L(X, \phi)=V(\phi)(-X+X^2)$ by considering the potential related parameter $\Gamma$($=\frac{V V''}{V'^2}$) as a function of another potential parameter $\lambda$($=\frac{V'}{\kappa V^{3/2}}$), which correspondingly extends the analysis of the dynamical system of our universe from two-dimension to three-dimension. Beside the critical points found in previous papers, we find a new de-Sitter-like dominant attractor(\textit{cp$6$}) and give its stable condition using the center manifold theorem. For the dark energy dominant solution(\textit{cp$6$} and \textit{cp$7$}), it could be distinguished from canonical quintessence and tachyon models since the sound speed $c_s^2=0$ or $c_s^2\ll 1$. For the scaling solution (\textit{cp$8$}), it is very interesting that the sound speed $c_s^2=1/5$ while it behaves as ordinary matter. We therefore point out that the power-law kinetic quintessence should have different signatures on cold dark matter power spectrum and cosmic microwave background both at early time when this scalar field is an early dark energy with $\Omega_\phi$ being non-negligible at high redshift and at late time when it drives the accelerating expansion. We even do not know whether there are any degeneracies of the impacts between these two epoches. They are expected to be investigated in future.

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

Using deep narrow-band $H_2S1$ and $K_{s}$-band imaging data obtained with CFHT/WIRCam, we identify a sample of 56 H$\alpha$ emission-line galaxies (ELGs) at $z=2.24$ with the 5$\sigma$ depths of $H_2S1=22.8$ and $K_{s}=24.8$ (AB) over 383 arcmin$^{2}$ area in the ECDFS. A detailed analysis is carried out with existing multi-wavelength data in this field. Three of the 56 H$\alpha$ ELGs are detected in Chandra 4 Ms X-ray observation and two of them are classified as AGNs. The rest-frame UV and optical morphologies revealed by HST/ACS and WFC3 deep images show that nearly half of the H$\alpha$ ELGs are either merging systems or with a close companion, indicating that the merging/interacting processes play a key role in regulating star formation at cosmic epoch z=2-3; About 14% are too faint to be resolved in the rest-frame UV morphology due to high dust extinction. We estimate dust extinction from SEDs. We find that dust extinction is generally correlated with H$\alpha$ luminosity and stellar mass (SM). Our results suggest that H$\alpha$ ELGs are representative of star-forming galaxies (SFGs). Applying extinction correction for individual objects, we examine the intrinsic H$\alpha$ luminosity function (LF) at $z=2.24$, obtaining a best-fit Schechter function characterized by a faint-end slope of $\alpha=-1.3$. This is shallower than the typical slope of $\alpha \sim -1.6$ in previous works based on constant extinction correction. We demonstrate that this difference is mainly due to the different extinction corrections. The proper extinction correction is thus key to recovering the intrinsic LF as the extinction globally increases with H$\alpha$ luminosity. Moreover, we find that our H$\alpha$ LF mirrors the SM function of SFGs at the same cosmic epoch. This finding indeed reflects the tight correlation between SFR and SM for the SFGs, i.e., the so-called main sequence.

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $\lambda_{\rm rest} <$ 3.5$ \mu$m; a hot dust component peaking at $\lambda_{rest} \sim$ 5$\,\mu$m; and cold dust component which becomes significant for $\lambda_{\rm rest} >$ 50$\,\mu$m. Six objects in the sample are detected at 250 and 350$ \mu$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ \mu$m detections corresponding to a rest-frame 4.5$ \mu$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered.

Starting from a sample of 24 \micron\ sources in the Extended Groth Strip, we use 3.6 to 8 \micron\ color criteria to select ultraluminous infrared galaxies (ULIRGs) at $z\sim2$. Spectroscopy from 20-38 \micron\ of 14 objects verifies their nature and gives their redshifts. Multi-wavelength data for these objects imply stellar masses ${>}10^{11}$ \Msun\ and star formation rates $\ge$410 \Msun yr$^{-1}$. Four objects of this sample observed at 1.6 \micron\ (rest-frame visible) with {\it HST}/WFC3 show diverse morphologies, suggesting that multiple formation processes create ULIRGs. Four of the 14 objects show signs of active galactic nuclei, but the luminosity appears to be dominated by star formation in all cases.

We report the first results from a spectroscopic survey of the [CII] 158um line from a sample of intermediate redshift (0.2<z<0.8) (ultra)-luminous infrared galaxies, (U)LIRGs (LIR>10^11.5 Lsun), using the SPIRE-Fourier Transform Spectrometer (FTS) on board the Herschel Space Observatory. This is the first survey of [CII] emission, an important tracer of star-formation, at a redshift range where the star-formation rate density of the Universe increases rapidly. We detect strong [CII] 158um line emission from over 80% of the sample. We find that the [CII] line is luminous, in the range (0.8-4)x10^(-3) of the far-infrared continuum luminosity of our sources, and appears to arise from photodissociation regions on the surface of molecular clouds. The L[CII]/LIR ratio in our intermediate redshift (U)LIRGs is on average ~10 times larger than that of local ULIRGs. Furthermore, we find that the L[CII]/LIR and L[CII]/LCO(1-0) ratios in our sample are similar to those of local normal galaxies and high-z star-forming galaxies. ULIRGs at z~0.5 show many similarities to the properties of local normal and high-z star forming galaxies. Our findings strongly suggest that rapid evolution in the properties of the star forming regions of luminous infrared galaxies is likely to have occurred in the last 5 billion years.

[abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5-6.2sigma in the very deep Spitzer/IRAC 4.5 micron data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z=10.2+-0.4) is robustly detected also at 3.6 micron (6.9sigma), revealing a flat UV spectral energy distribution with a slope beta=-2.0+-0.2, consistent with demonstrated trends with luminosity at high redshift. The abundance of such luminous candidates suggests that the luminosity function evolves more significantly in phi_* than in L_* at z>~8 with a higher number density of bright sources than previously expected. Despite the discovery of these luminous candidates, the cosmic star formation rate density for galaxies with SFR >0.7 M_sun/yr shows an order-of-magnitude increase in only 170 Myr from z ~ 10 to z ~ 8, consistent with previous results. Based on the IRAC detections, we derive galaxy stellar masses at z~10, finding that these luminous objects are typically 10^9 M_sun. The cosmic stellar mass density at z~10 is log10 rho_* = 4.7^+0.5_-0.8 M_sun Mpc^-3 for galaxies brighter than M_UV~-18. The remarkable brightness, and hence luminosity, of these z~9-10 candidates highlights the opportunity for deep spectroscopy to determine their redshift and nature, demonstrates the value of additional search fields covering a wider area to understand star-formation in the very early universe, and highlights the opportunities for JWST to map the buildup of galaxies at redshifts much earlier than z~10.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

We present the results of a combined galaxy population analysis for the host galaxies of active galactic nuclei (AGN) identified at 0 < z < 1.4 within the SDSS, Bootes and DEEP2 surveys. We identified AGN in a uniform and unbiased manner at X-ray, infrared and radio wavelengths. Supermassive black holes undergoing radiatively-efficient accretion (detected as X-ray and/or infrared AGN) appear to be hosted in a separate and distinct galaxy population than AGN undergoing powerful mechanically dominated accretion (radio AGN). Consistent with some previous studies, radiatively efficient AGN appear to be preferentially hosted in modest star-forming galaxies, with little dependence on AGN or galaxy luminosity. AGN exhibiting radio-emitting jets due to mechanically-dominated accretion are almost exclusively observed in massive, passive galaxies. Crucially, we now provide strong evidence that the observed host-galaxy trends are independent of redshift. In particular, these different accretion-mode AGN have remained as separate galaxy populations throughout the last 9 Gyr. Furthermore, it appears that galaxies hosting AGN have evolved along the same path as galaxies that are not hosting AGN with little evidence for distinctly separate evolution.

Out of several dozen z > 7 candidate galaxies observed spectroscopically, only five have been confirmed via Lyman-alpha emission, at z=7.008, 7.045, 7.109, 7.213 and 7.215. The small fraction of confirmed galaxies may indicate that the neutral fraction in the intergalactic medium (IGM) rises quickly at z > 6.5, as Lyman-alpha is resonantly scattered by neutral gas. However, the small samples and limited depth of previous observations makes these conclusions tentative. Here we report the results of a deep near-infrared spectroscopic survey of 43 z > 6.5 galaxies. We detect only a single galaxy, confirming that some process is making Lyman-alpha difficult to detect. The detected emission line at 1.0343 um is likely to be Lyman-alpha emission, placing this galaxy at a redshift z = 7.51, an epoch 700 million years after the Big Bang. This galaxy's colors are consistent with significant metal content, implying that galaxies become enriched rapidly. We measure a surprisingly high star formation rate of 330 Msol/yr, more than a factor of 100 greater than seen in the Milky Way. Such a galaxy is unexpected in a survey of our size, suggesting that the early universe may harbor more intense sites of star-formation than expected.

The properties of K-band selected galaxies (K_AB<24) in the z = 3.09 SSA22 protocluster field are studied. 430 galaxies at 2.6 < z_phot < 3.6 are selected as potential protocluster members in a 112 arcmin^2 area based on their photometric redshifts. We find that \approx 20% of the massive galaxies with stellar masses >10^11 M_sun at z_phot \sim 3.1 have colors consistent with those of quiescent galaxies with ages > 0.5 Gyr. This fraction increases to \approx 50% after correcting for unrelated foreground/background objects. We also find that 30% of the massive galaxies are heavily reddened dusty star-forming galaxies. Few such quiescent galaxies at similar redshifts are seen in typical survey fields. An excess surface density of 24\mu m sources at z_phot \sim 3.1 is also observed, implying the presence of dusty star-formation activity in the protocluster. Cross-correlation with the X-ray data indicates that the fraction of K-band selected protocluster galaxies hosting active galactic nuclei (AGN) is also high compared with the field. The sky distribution of the quiescent galaxies, the 24\mu m sources, and the X-ray AGNs show clustering around a density peak of z=3.1 Ly\alpha emitters (LAEs). A significant fraction of the massive galaxies have already become quiescent, while the dusty star-formation is still active in the SSA22 protocluster. These findings indicate that we are witnessing the formation epoch of massive early-type galaxies at the center of predecessors to present-day rich galaxy clusters.

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

(Abridged) In tandem with observational datasets, we utilize realistic mock catalogs, based on a semi-analytic galaxy formation model, constructed specifically for Pan-STARRS1 Medium Deep Surveys in order to assess the performance of the Probability Friends-of-Friends (PFOF, Liu et al.) group finder, and aim to develop a grouping optimization method applicable to surveys like Pan-STARRS1. Producing mock PFOF group catalogs under a variety of photometric redshift accuracies ({\sigma}{\Delta}z/(1+zs)), we find that catalog purities and completenesses from ``good' {\sigma}{\Delta}z/(1+zs)) ~ 0.01) to ``poor' {\sigma}{\Delta}z/(1+zs)) ~ 0.07) photo-zs gradually degrade respectively from 77% and 70% to 52% and 47%. To avoid model dependency of the mock for use on observational data we apply a ``subset optimization' approach, using spectroscopic-redshift group data from the target field to train the group finder for application to that field, as an alternative method for the grouping optimization. We demonstrate this approach using these spectroscopically identified groups as the training set, i.e. zCOSMOS groups for PFOF searches within PS1 Medium Deep Field04 (PS1MD04) and DEEP2 EGS groups for searches in PS1MD07. We ultimately apply PFOF to four datasets spanning the photo-z uncertainty range from 0.01 to 0.06 in order to quantify the dependence of group recovery performance on photo-z accuracy. We find purities and completenesses calculated from observational datasets broadly agree with their mock analogues. Further tests of the PFOF algorithm are performed via matches to X-ray clusters identified within the PS1MD04 and COSMOS footprints. Across over a decade in group mass, we find PFOF groups match ~85% of X-ray clusters in COSMOS and PS1MD04, but at a lower statistical significance in the latter.

A new EAS hybrid experiment has been designed by constructing a YAC (Yangbajing Air shower Core) detector array inside the existing Tibet-III air shower array. The first step of YAC, called "YAC-I", consists of 16 plastic scintillator units (4 rows times 4 columns) each with an area of 40 cm * 50 cm which is used to check hadronic interaction models used in AS simulations. A Monte Carlo study shows that YAC-I can record high energy electromagnetic component in the core region of air showers induced by primary particles of several tens TeV energies where the primary composition is directly measured by space experiments. It may provide a direct check of the hadronic interaction models currently used in the air shower simulations in the corresponding energy region. In present paper, the method of the observation and the sensitivity of the characteristics of the observed events to the different interaction models are discussed.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

After carefully cross-identifying a previously discovered GALEX-selected Lyman Break Galaxy (LBG) candidates one-to-one with their optical counterparts in the field of the CDF-S, we re-estimate their photometric redshifts using multi-wavelength data from UV, optical to NIR. We refine a new updated sample of 383 LBGs at $0.7\la z \la 1.4$. Most LBGs are classified as starburst and irregular types. Ages spread from several Myr to 1.5Gyr. Their dust-corrected star formation rates (SFRs) and stellar masses ($M_*$) are from $4\my$ to $220\my$ and from $2.3\times 10^8 \msun$ to $4 \times 10^{11} \msun$. The rest-frame FUV luminosity function of LBGs are presented. LBGs of irregular types mainly distribute along the "main sequence" of star forming galaxies while most LBGs of starburst types locate in the starburst region. A "downsizing" effect is clearly found and LBGs distribute in the "blue" cloud. HST images in F606W ($V$ band) and F850LP ($z$ band) are taken from the GEMS and GOODS-S surveys. SExtractor and GALFIT are applied to get their morphological parameters. A morphological sample of 142 LBGs with reliable results of $\sersic$ and sizes in both bands is defined. We find that LBGs at $z\sim 1 $ are dominated by disk-like galaxies. Correlations between photometric and morphological properties of LBGs are investigated. Strong correlations between their half-light radii and $M_*$, i.e., size-stellar mass relations, are found in both bands. Physical connections between correlations and the "downsizing" effect are discussed.

We present the analysis of Spitzer/IRAC and NIR imaging observation of AFGL 5157, an active star forming region. In the IRAC images, this region shows strong PAH emissions in channel 4 and H$_2$ emissions in channel 2. Many of the H$_2$ features are aligned to form jet-like structures. Three bipolar jets in the NH$_3$ core region and a couple of jets northwest of the core have been identified. We identify the possible driving agents of the bipolar jets and show them to be very young. An embedded cluster has been detected in the NH$_3$ core; many members in the cluster show their SEDs increasing from JHK bands toward longer wavelength, indicative of their early evolutionary stages. Millimeter and sub-millimeter continuum emissions are found to coincide spatially with these presumable Class 0/I sources, in the NH$_3$ core and the NW subregion. The existence of H$_2$ bipolar jets and very young stellar objects suggests that star formation is still going on at present epoch in these subregions. Combining the information from previous studies, we propose a sequential star formation scenario in the whole AFGL 5157 region.

We present a multi-wavelength study of a 3.6 $\mu$m-selected galaxy sample in the Extended Groth strip. The sample is complete for galaxies with stellar mass $>10^{9.5}$ \Msun and redshift $0.4<z<1.2$. In this redshift range, the IRAC 3.6 $\mu$m band measures the rest-frame near-infrared band, permitting nearly unbiased selection with respect to both quiescent and star-forming galaxies. The numerous spectroscopic redshifts available in the EGS are used to train an Artificial Neural Network to estimate photometric redshifts. The distribution of photometric redshift errors is Gaussian with standard deviation ${\sim}0.025(1+z)$, and the fraction of redshift failures (${>}3\sigma$ errors) is about 3.5%. A new method of validation based on pair statistics confirms the estimate of standard deviation even for galaxies lacking spectroscopic redshifts. Basic galaxy properties measured include rest-frame $U-B$ colors, $B$- and $K$-band absolute magnitudes, and stellar masses. We divide the sample into quiescent and star-forming galaxies according to their rest-frame $U-B$ colors and 24 to 3.6 \micron\ flux density ratios and derive rest $K$-band luminosity functions and stellar mass functions for quiescent, star forming, and all galaxies. The results show that massive, quiescent galaxies were in place by $z\approx1$, but lower mass galaxies generally ceased their star formation at later epochs.

The shutdown of star formation in galaxies is generally termed `quenching'. Although quenching may occur through a variety of processes, the exact mechanism(s) that is in fact responsible for quenching is still in question. This paper addresses quenching by searching for traces of possible quenching processes through their effects on galaxy structural parameters such as surface stellar mass density and Sersic index (n). We analyze the rest-frame U-B color correlations versus these structural parameters using a sample of galaxies in the redshift range 0.5< z<0.8 from the DEEP2/AEGIS survey. We find that Sersic index (n) has the smallest overlap region among all tested parameters and resembles a step-function with a threshold value of n=2.3. There exists, however, a significant population of outliers with blue colors yet high n values that seem to contradict this behavior. We hypothesize that their Sersic values may be distorted by bursts of star formation, AGNs, and/or poor fits, leading us to consider central surface stellar mass density as an alternative to Sersic index. Not only does it correct the outliers, it also forms a tight relationship with color, suggesting that the innermost structure of galaxies is most physically linked with quenching. Furthermore, at z~0.65, the majority of the blue cloud galaxies cannot simply fade onto the red sequence since their GIM2D bulge masses are only half as large on average as the bulge masses of similar red sequence galaxies, thus demonstrating that stellar mass must absolutely increase at the centers of galaxies as they quench. We discuss a two-stage model for quenching in which galaxy star formation rates are controlled by their dark halos while they are still in the blue cloud and a second quenching process sets in later, associated with the central stellar mass build-up.

We propose a model for Dirac leptogenesis based on a non-anomalous $U(1)^{\prime}$ gauged family symmetry. The anomaly cancellation conditions are satisfied with no new chiral fermions other than the three right-handed neutrinos, giving rise to stringent constraints among the charges. Realistic masses and mixing angles are obtained for all fermions. The model predicts neutrinos of the Dirac type with naturally suppressed masses. Dirac leptogenesis is achieved through the decay of the flavon fields. The cascade decays of the vector-like heavy fermions in the Froggatt-Nielsen mechanism play a crucial role in the separation of the primodial lepton numbers. We find that a large region of parameter space of the model gives rise to a sufficient cosmological baryon number asymmetry through Dirac leptogenesis.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

Infrared 3.6 to 8 micron images of the Extended Groth Strip yield plausible counterpart identifications for all but one of 510 radio sources in the AEGIS20 S(1.4 GHz) > 50 micro-Jy sample. This is the first such deep sample that has been effectively 100% identified. Achieving the same identification rate at R-band would require observations reaching R_AB > 27. Spectroscopic redshifts are available for 46% of the sample and photometric redshifts for an additional 47%. Almost all of the sources with 3.6 micron AB magnitudes brighter than 19 have spectroscopic redshifts z < 1.1, while fainter objects predominantly have photometric redshifts with 1 \lapprox z \lapprox 3. Unlike more powerful radio sources that are hosted by galaxies having large stellar masses within a relatively narrow range, the AEGIS20 counterparts have stellar masses spanning more than a factor of 10 at z \sim 1. The sources are roughly 10--15% starbursts at z \lapprox 0.5 and 20--25% AGNs mostly at z > 1 with the remainder of uncertain nature.

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

We present a study on Spectral Energy Distributions, Morphologies, and star formation for an IRAC-selected extremely red object sample in the GOODS Chandra Deep Field-South. This work was enabled by new HST/WFC3 near-IR imaging from the CANDELS survey as well as the deepest available X-ray data from Chandra 4 Ms observations. This sample consists of 133 objects with the 3.6um limiting magnitude of [3.6] = 21.5, and is approximately complete for galaxies with M >10^{11}M_sun at 1.5 < z < 2.5. We classify this sample into two types, quiescent and star-forming galaxies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further morphological study of this sample show a consistent result with the observed color classification. The classified quiescent galaxies are bulge dominated and star-forming galaxies in the sample have disk or irregular morphologies. Our observed infrared color classification is also consistent with the rest-frame color (U-V vs V-J) classification. We also found that quiescent and star-forming galaxies are well separated in the nonparametric morphology parameter (Gini vs M_{20}) diagram measuring their concentration and clumpiness: quiescent galaxies have Gini coefficient higher than 0.58 and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star formation quenching process must lead to or be accompanied by the increasing galaxy concentration. One prominent morphological feature of this sample is that disks are commonly seen in this massive galaxy sample at 1.5 < z < 2.5: 30% of quiescent galaxies and 70% of star forming galaxies with M >10^{11}M_sun have disks in their rest-frame optical morphologies. The prevalence of these extended, relatively undisturbed disks challenges the merging scenario as the main mode of massive galaxy formation.

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

We report detection of four IRAC sources in the GOODS-South field with an extremely red color of H$-$[3.6]$>$4.5. The four sources are not detected in the deep HST WFC3 H-band image with H$_{limit}$=28.3 mag. We find that only 3 types of SED templates can produce such a red H$-$[3.6] color: a very dusty SED with the Calzetti extinction of A$_V$=16 mag at z=0.8; a very dusty SED with the SMC extinction of A$_V$=8 mag at z=2.0$\sim$2.2; and an 1Gyr SSP with A$_V \sim$0.8 at z=5.7. We argue that these sources are unlikely dusty galaxies at z$\leq$2.2 based on absent strong MIPS 24$\mu$m emission. The old stellar population model at z$>$4.5 remains a possible solution for the 4 sources. At z$>$4.5, these sources have stellar masses of Log(M$_*$/M$_{\odot}$)=10.6$\sim$11.2. One source, ERS-1, is also a type-II X-ray QSO with L$_{2-8keV}$=1.6$\times 10^{44}$ erg s$^{-1}$. One of the four sources is an X-ray QSO and another one is a HyperLIRG, suggesting a galaxy-merging scenario for the formation of these massive galaxies at high redshifts.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present broad band photometry and photometric redshifts for 187611 sources located in ~0.5deg^2 in the Lockman Hole area. The catalog includes 389 X-ray detected sources identified with the very deep XMM-Newton observations available for an area of 0.2 deg^2. The source detection was performed on the Rc, z' and B band images and the available photometry is spanning from the far ultraviolet to the mid infrared, reaching in the best case scenario 21 bands. Astrometry corrections and photometric cross-calibrations over the entire dataset allowed the computation of accurate photometric redshifts. Special treatment is undertaken for the X-ray sources, the majority of which is active galactic nuclei. Comparing the photometric redshifts to the available spectroscopic redshifts we achieve for normal galaxies an accuracy of \sigma_{\Delta z/(1+z)}=0.036, with 12.7% outliers, while for the X-ray detected sources the accuracy is \sigma_{\Delta z/(1+z)}=0.069, with 18.3% outliers, where the outliers are defined as sources with |z_{phot}-z_{spec}|>0.15 (1+z_{spec})}. These results are a significant improvement over the previously available photometric redshifts for normal galaxies in the Lockman Hole, while it is the first time that photometric redshifts are computed and made public for AGN for this field.

Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy. However, the overall spectrum of the EBL between 0.1 and 1000 microns has never been determined directly neither from observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncertainties in our modelling based directly on the data are quantified, and their consequences for attenuation of very-high-energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level roughly matching galaxy count data. Independent efforts from IR and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

We study the self-gravitational effect on the equation of state (EoS) of Bose and Fermi gases in thermal equilibrium at the end of reheating, the period after quark-hadron transition and before Big Bang Nucleosynthesis (BBN). After introducing new grand canonical partition functions based on the work of Uhlenbeck and Gropper, we notice some interesting features of the newly developed EoSs with distinct behaviors of relativistic and non-relativistic gases under self-gravity. The usual negligence of the self-gravitational effect when solving the background expansion of the early universe is justified with numerical results, showing the magnitude of the self-gravitational modification of the state constant to be less than $O(10^{-78})$. This helps us to clarify the background thermal evolution of the primordial patch. Such clarification is crucial in testing gravity theories, evaluating inflation models and determining element abundances in BBN.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy, but the overall spectrum of the EBL between 0.1-1000 microns has never been determined directly from galaxy spectral energy distribution (SED) observations over a wide redshift range. The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. This is achieved from the observed evolution of the rest-frame K-band galaxy luminosity function up to redshift 4 (Cirasuolo et al. 2010), combined with a determination of galaxy SED-type fractions. These are based on fitting SWIRE templates to a multiwavelength sample of about 6000 galaxies in the redshift range from 0.2 to 1 from the All-wavelength Extended Groth Strip International Survey (AEGIS). The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and AGN galaxies in that redshift range are estimated, and two alternative extrapolations of SED-types to higher redshifts are considered. This allows calculation of the evolution of the luminosity densities from the UV to the IR, the evolving star formation rate density of the universe, the evolving contribution to the bolometric EBL from the different galaxy populations including AGN galaxies and the buildup of the EBL. Our EBL calculations are compared with those from a semi-analytic model, from another observationally-based model and observational data. The EBL uncertainties in our modeling based directly on the data are quantified, and their consequences for attenuation of very high energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL is well constrained from the UV to the mid-IR, but independent efforts from infrared and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

We study the evolution of the scaling relations between maximum circular velocity, stellar mass and optical half-light radius of star-forming disk-dominated galaxies in the context of LCDM-based galaxy formation models. Using data from the literature combined with new data from the DEEP2 and AEGIS surveys we show that there is a consistent observational and theoretical picture for the evolution of these scaling relations from z\sim 2 to z=0. The evolution of the observed stellar scaling relations is weaker than that of the virial scaling relations of dark matter haloes, which can be reproduced, both qualitatively and quantitatively, with a simple, cosmologically-motivated model for disk evolution inside growing NFW dark matter haloes. In this model optical half-light radii are smaller, both at fixed stellar mass and maximum circular velocity, at higher redshifts. This model also predicts that the scaling relations between baryonic quantities evolve even more weakly than the corresponding stellar relations. We emphasize, though, that this weak evolution does not imply that individual galaxies evolve weakly. On the contrary, individual galaxies grow strongly in mass, size and velocity, but in such a way that they move largely along the scaling relations. Finally, recent observations have claimed surprisingly large sizes for a number of star-forming disk galaxies at z \sim 2, which has caused some authors to suggest that high redshift disk galaxies have abnormally high spin parameters. However, we argue that the disk scale lengths in question have been systematically overestimated by a factor \sim 2, and that there is an offset of a factor \sim 1.4 between H\alpha sizes and optical sizes. Taking these effects into account, there is no indication that star forming galaxies at high redshifts (z\sim 2) have abnormally high spin parameters.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We present first insights into the far-IR properties for a sample of IRAC and MIPS-24um detected Lyman Break Galaxies (LBGs) at z ~ 3, as derived from observations in the northern field of the Great Observatories Origins Survey (GOODS-N) carried out with the PACS instrument on board the Herschel Space Observatory. Although none of our galaxies are detected by Herschel, we employ a stacking technique to construct, for the first time, the average spectral energy distribution of infrared luminous LBGs from UV to radio wavelengths. We derive a median IR luminosity of L_{IR} = 1.6 x 10^12 Lo, placing the population in the class of ultra luminous infrared galaxies (ULIRGs). Complementing our study with existing multi-wavelength data, we put constraints on the dust temperature of the population and find that for their L_{IR}, MIPS-LBGs are warmer than submm-luminous galaxies while they fall in the locus of the L_{IR}-T_{d} relation of the local ULIRGs. This, along with estimates based on the average SED, explains the marginal detection of LBGs in current sub-mm surveys and suggests that these latter studies introduce a bias towards the detection of colder ULIRGs in the high-z universe, while missing high-z ULIRGS with warmer dust.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We report the detection and successful modeling of the unusual 9.7\mum Si--O stretching silicate emission feature in the type 1 (i.e. face-on) LINER nucleus of M81. Using the Infrared Spectrograph (IRS) instrument on Spitzer, we determine the feature in the central 230 pc of M81 to be in strong emission, with a peak at ~10.5\mum. This feature is strikingly different in character from the absorption feature of the galactic interstellar medium, and from the silicate absorption or weak emission features typical of galaxies with active star formation. We successfully model the high signal-to-noise ratio IRS spectra with porous silicate dust using laboratory-acquired mineral spectra. We find that the most probable fit uses micron-sized, porous grains of amorphous silicate and graphite. In addition to silicate dust, there is weak PAH emission present (particularly at 11.3\mum, arising from the C--H out-of-plane bending vibration of relatively large PAHs of ~500--1000 C atoms) whose character reflects the low-excitation AGN environment, with some evidence that small PAHs of ~100--200 C atoms (responsible for the 7.7\mum C--C stretching band) in the immediate vicinity of the nucleus have been preferentially destroyed. (abstract continues)

We present high-resolution 345 GHz interferometric observations of two extreme luminous (L_{IR}>10^{13} L_sun), submillimetre-selected galaxies (SMGs) in the COSMOS field with the Submillimeter Array (SMA). Both targets were previously detected as unresolved point-sources by the SMA in its compact configuration, also at 345 GHz. These new data, which provide a factor of ~3 improvement in resolution, allow us to measure the physical scale of the far-infrared in the submillimetre directly. The visibility functions of both targets show significant evidence for structure on 0.5-1 arcsec scales, which at z=1.5 translates into a physical scale of 5-8 kpc. Our results are consistent with the angular and physical scales of two comparably luminous objects with high-resolution SMA followup, as well as radio continuum and CO sizes. These relatively compact sizes (<5-10 kpc) argue strongly for merger-driven starbursts, rather than extended gas-rich disks, as the preferred channel for forming SMGs. For the most luminous objects, the derived sizes may also have important physical consequences; under a series of simplifying assumptions, we find that these two objects in particular are forming stars close to or at the Eddington limit for a starburst.

We present Spitzer observations for a sample of close major-merger galaxy pairs (KPAIR sample) selected from 2MASS/SDSS-DR3 cross-matches. The goals are to study the star formation activity in these galaxies and to set a local bench mark for the cosmic evolution of close major mergers. The Spitzer KPAIR sample (27 pairs, 54 galaxies) includes all spectroscopically confirmed S+S and S+E pairs in a parent sample that is complete for primaries brighter than K=12.5 mag, projected separations of 5< s < 20 kpc/h, and mass ratios<2.5. The Spitzer data consist of images in 7 bands (3.6, 4.5, 5.8, 8, 24, 70, 160 um). Compared to single spiral galaxies in a control sample, only spiral galaxies in S+S pairs show significantly enhanced specific star formation rate (sSFR=SFR/M), whereas spiral galaxies in S+E pairs do not. Furthermore, the SFR enhancement of spiral galaxies in S+S pairs is highly mass-dependent. Only those with $\rm M \gsim 10^{10.5} M_\sun$ show significant enhancement. Relatively low mass ($\rm M \sim 10^{10} M_\sun$) spirals in S+S pairs have about the same SFR/M compared to their counterparts in the control sample. There is evidence for a correlation between the global star formation activities (but not the nuclear activities) of the component galaxies in massive S+S major-merger pairs (the "Holmberg effect"). There is no significant difference in the SFR/M between the primaries and the secondaries, nor between spirals of SEP<1 and those of SEP.1. The contribution of KPAIR galaxies to the cosmic SFR density in the local universe is only 1.7%.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

Based on the joint-observations of the radio broadband spectral emissions of solar eclipse on August 1, 2008 at Jiuquan (total eclipse) and Huairou (partial eclipse) at the frequencies of 2.00 -- 5.60 GHz (Jiuquan), 2.60 -- 3.80 GHZ (Chinese solar broadband radiospectrometer, SBRS/Huairou), and 5.20 -- 7.60 GHz (SBRS/Huairou), the authors assemble a successive series of broadband spectrum with a frequency of 2.60 -- 7.60 GHz to observe the solar eclipse synchronously. This is the first attempt to analyze the solar eclipse radio emission under the two telescopes located at different places with broadband frequencies in the periods of total and partial eclipse. With these analyses, the authors made a new semiempirical model of the coronal plasma density of the quiet Sun and made a comparison with the classic models.

We performed an analysis of the mid-infared properties of the 12micron Seyfert sample, a complete unbiased 12micron flux limited sample of local Seyfert galaxies selected from the IRAS Faint Source Catalog, based on low resolution spectra obtained with the Infrared Spectrograph (IRS) on-board Spitzer Space Telescope. A detailed presentation of this analysis is dicussed in Wu et al. (2009). We find that on average, the 15-30micron slope of the continuum is -0.85+/-0.61 for Seyfert 1s and -1.53+/-0.84 for Seyfert 2s, and there is substantial scatter in each type. Moreover, nearly 32% of Seyfert 1s, and 9% of Seyfert 2s, display a peak in the mid-infrared spectrum at 20micron, which is attributed to an additional hot dust component. The Polycyclic Aromatic Hydrocarbon (PAH) equivalent width decreases with increasing dust temperature, as indicated by the global infrared color of the host galaxies. However, no statistical difference in PAH equivalent width is detected between the two Seyfert types, 1 and 2, of the same bolometric luminosity. Finally, we propose a new method to estimate the AGN contribution to the integrated 12micron galaxy emission, by subtracting the "star formation" component in the Seyfert galaxies, making use of the tight correlation between PAH 11.2micron luminosity and 12micron luminosity for star forming galaxies.

We have undertaken a deep (sigma~1.1 mJy) 1.1-mm survey of the z=0.54 cluster MS 0451.6-0305 using the AzTEC camera on the James Clerk Maxwell Telescope. We detect 36 sources with S/N>3.5 in the central 0.10 deg^2 and present the AzTEC map, catalogue and number counts. We identify counterparts to 18 sources (50%) using radio, mid-infrared, Spitzer IRAC and Submillimeter Array data. Optical, near- and mid-infrared spectral energy distributions are compiled for the 14 of these galaxies with detectable counterparts, which are expected to contain all likely cluster members. We then use photometric redshifts and colour selection to separate background galaxies from potential cluster members and test the reliability of this technique using archival observations of submillimetre galaxies. We find two potential MS 0451-03 members, which, if they are both cluster galaxies have a total star-formation rate (SFR) of ~100 solar masses per year -- a significant fraction of the combined SFR of all the other galaxies in MS 0451-03. We also examine the stacked rest-frame mid-infrared, millimetre and radio emission of cluster members below our AzTEC detection limit and find that the SFRs of mid-IR selected galaxies in the cluster and redshift-matched field populations are comparable. In contrast, the average SFR of the morphologically classified late-type cluster population is ~3 times less than the corresponding redshift-matched field galaxies. This suggests that these galaxies may be in the process of being transformed on the red-sequence by the cluster environment. Our survey demonstrates that although the environment of MS 0451-03 appears to suppress star-formation in late-type galaxies, it can support active, dust-obscured mid-IR galaxies and potentially millimetre-detected LIRGs.

We present new observational determinations of the evolution of the 2-10keV X-ray luminosity function (XLF) of AGN. We utilise data from a number of surveys including both the 2Ms Chandra Deep Fields and the AEGIS-X 200ks survey, enabling accurate measurements of the evolution of the faint end of the XLF. We combine direct, hard X-ray selection and spectroscopic follow-up or photometric redshift estimates at z<1.2 with a rest-frame UV colour pre-selection approach at higher redshifts to avoid biases associated with catastrophic failure of the photometric redshifts. Only robust optical counterparts to X-ray sources are considered using a likelihood ratio matching technique. A Bayesian methodology is developed that considers redshift probability distributions, incorporates selection functions for our high redshift samples, and allows robust comparison of different evolutionary models. We find that the XLF retains the same shape at all redshifts, but undergoes strong luminosity evolution out to z~1, and an overall negative density evolution with increasing redshift, which thus dominates the evolution at earlier times. We do not find evidence that a Luminosity-Dependent Density Evolution, and the associated flattening of the faint-end slope, is required to describe the evolution of the XLF. We find significantly higher space densities of low-luminosity, high-redshift AGN than in prior studies, and a smaller shift in the peak of the number density to lower redshifts with decreasing luminosity. The total luminosity density of AGN peaks at z=1.2+/-0.1, but there is a mild decline to higher redshifts. We find >50% of black hole growth takes place at z>1, with around half in Lx<10^44 erg/s AGN.

[abridged] We present constraints on the nature of the first galaxies selected at 350 microns. The sample includes galaxies discovered in the deepest blank-field survey at 350 microns (in the Bootes Deep Field) and also later serendipitous detections in the Lockman Hole. Spectral energy distribution templates are fit to identified counterparts, and the sample is found to comprise IR-luminous galaxies at 1<z<3 predominantly powered by star formation. The first spectrum of a 350-micron-selected galaxy provides an additional confirmation, showing prominent dust grain features typically associated with star-forming galaxies. Compared to submillimeter galaxies selected at 850 and 1100 microns, galaxies selected at 350 microns have a similar range of far-infrared color temperatures. However, no 350-micron-selected sources are reliably detected at 850 or 1100 microns. Galaxies in our sample with redshifts 1<z<2 show a tight correlation between the far- and mid-infrared flux densities, but galaxies at higher redshifts show a large dispersion in their mid- to far-infrared colors. The 350-micron flux densities (15<S(350)<40 mJy) place these objects near the Herschel/SPIRE 350-micron confusion threshold, with the lower limit on the star formation rate density suggesting the bulk of the 350-micron contribution will come from less luminous infrared sources and normal galaxies.

We present results of a large survey of the mid--IR properties of 248 Lyman Break Galaxies with confirmed spectroscopic redshift using deep Spitzer/IRAC observations in six cosmological fields. We model the Spectral Energy Distributions (SEDs) employing a revised version of the Bruzual and Charlot synthesis population code that incorporates a new treatment of the TP--AGB phase (CB07). Our primary aim is to investigate the impact of the AGB phase in the stellar masses of the LBGs, and compare our new results with previous stellar mass estimates. Based on the new CB07 code we find that the stellar masses of LBGs are smaller on average by a factor of ~1.4 compared to previous estimates. LBGs with 8um and/or 24um detections show higher masses (M~10^11 Mo) than LBGs faint in the IRAC bands (M~10^9 Mo). The ages of these massive LBGs are considerably higher than the rest of the population, indicating that they have been star-forming for at least ~1 Gyr. We also show how the addition of the IRAC bands, improves the accuracy of the estimated stellar masses and reduced the scatter on the derived M/L ratios. In particular, we present a tight correlation between the 8um IRAC band (rest-frame K for galaxies at z~3) and the stellar mass. We calculate the number density of massive (M > 10^11 Mo) LBGs and find it to be $\Phi$= (1.12 $\pm$ 0.4) x 10^(-5) Mpc^(-3), ~1.5 times lower than that found by previous studies. Finally, based on UV-corrected SFRs we investigate the SFR-stellar mass correlation at z~3, find it similar to the one observed at other redshifts and show that our data place the peak of the evolution of the specific star formation rate at z~3.

We present results from a continuing interferometric survey of high-redshift submillimeter galaxies with the Submillimeter Array, including high-resolution (beam size ~2 arcsec) imaging of eight additional AzTEC 1.1mm selected sources in the COSMOS Field, for which we obtain six reliable (peak S/N>5 or peak S/N>4 with multiwavelength counterparts within the beam) and two moderate significance (peak S/N>4) detections. When combined with previous detections, this yields an unbiased sample of millimeter-selected SMGs with complete interferometric followup. With this sample in hand, we (1) empirically confirm the radio-submillimeter association, (2) examine the submillimeter morphology - including the nature of submillimeter galaxies with multiple radio counterparts and constraints on the physical scale of the far infrared - of the sample, and (3) find additional evidence for a population of extremely luminous, radio-dim submillimeter galaxies that peaks at higher redshift than previous, radio-selected samples. In particular, the presence of such a population of high-redshift sources has important consequences for models of galaxy formation - which struggle to account for such objects even under liberal assumptions - and dust production models given the limited time since the Big Bang.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

We present new mid-infrared and submillimetre observations for a sample of eight high redshift type-2 QSOs located in the Chandra Deep Field South. The sources are X-ray absorbed with luminosities in excess of 10^44 erg/s. Two of the targets have robust detections, S/N > 4, while a further three targets are marginally detected with S/N > =2.5. All sources are detected in multiple mid-infrared bands with the Spitzer Space Telescope. The multiwavelength spectral energy distributions (SEDs) of the type-2 QSOs are compared to those of two local ultraluminous galaxies (Arp220 and IR22491) in order to assess contributions from a star-forming component in various parts of the SED. We suggest that their submillimetre emission is possibly due to a starburst while a large fraction of the mid-infrared energy is likely to originate in the obscured central quasar. Using the mid-infrared and submm observations we derive infrared luminosities which are found to be in excess of L>10^12Lsun. The submillimetre (850micron) to X-ray (2 keV) spectral indices (alpha_SX) span a wide range. About half of the type-2 QSOs have values typical for a Compton-thick AGN with only 1 per cent of the nuclear emission seen through scattering and, the remaining with values typical of submm-bright galaxies. Combining the available observational evidence we outline a possible scenario for the early stages of evolution of these sources.

We present the Spitzer Space Telescope Infrared Array Camera (IRAC) observations for a sample of local elliptical galaxies to study later stages of AGN activities. A sample of 36 elliptical galaxies is selected from the Palomar spectroscopic survey. We detect nuclear non-stellar infrared emission in 9 of them. There is unambiguous evidence of circumnuclear dust in these 9 galaxies in their optical images. We also find a remarkable correlation between the infrared excess emission and the nuclear radio/X-ray emission, suggesting that infrared excess emission is tightly related with nuclear activity. Possible origin of infrared excess emission from hot dust heated by the central AGN is supported by spectral indices of IR excess emission.

We present low-resolution 5.5-35 um spectra for 103 galaxies from the 12 um Seyfert sample, a complete unbiased 12 um flux limited sample of local Seyfert galaxies selected from the IRAS Faint Source Catalog, obtained with the Infrared Spectrograph (IRS) on-board Spitzer Space Telescope. For 70 of the sources observed in the IRS mapping mode, uniformly extracted nuclear spectra are presented for the first time. We performed an analysis of the continuum emission, the strength of the Polycyclic Aromatic Hydrocarbon (PAH) and astronomical silicate features of the sources. We find that on average, the 15-30 um slope of the continuum is alpha_{15-30}=-0.85+-0.61 for Seyfert 1s and -1.53+-0.84 for Seyfert 2s, and there is substantial scatter in each type. Moreover, nearly 32% of Seyfert 1s, and 9% of Seyfert 2s, display a peak in the mid-infrared spectrum at 20 um, which is attributed to an additional hot dust component. The PAH equivalent width decreases with increasing dust temperature, asindicated by the global infrared color of the host galaxies. However, no statistical difference in PAH equivalent width is detected between the two Seyfert types, 1 and 2, of the same bolometric luminosity. The silicate features at 9.7 and 18um in Seyfert 1 galaxies are rather weak, while Seyfert 2s are more likely to display strong silicate absorption. Those Seyfert 2s with the highest silicate absorption also have high infrared luminosity and high absorption (hydrogen column density N_H>10^23 cm^-2 as measured from the X-rays. Finally, we propose a new method to estimate the AGN contribution to the integrated 12 um galaxy emission, by subtracting the "star formation" component in the Seyfert galaxies, making use of the tight correlation between PAH 11.2 um luminosity and 12 um luminosity for star forming galaxies.

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we present low-resolution (64 < lambda / dlambda < 124), mid-infrared (20-38 micron) spectra of 23 high-redshift ULIRGs detected in the Bootes field of the NOAO Deep Wide-Field Survey. All of the sources were selected to have 1) fnu(24 micron) > 0.5 mJy; 2) R-[24] > 14 Vega mag; and 3) a prominent rest-frame 1.6 micron stellar photospheric feature redshifted into Spitzer's 3-8 micron IRAC bands. Of these, 20 show emission from polycyclic aromatic hydrocarbons (PAHs), usually interpreted as signatures of star formation. The PAH features indicate redshifts in the range 1.5 < z < 3.0, with a mean of <z>=1.96 and a dispersion of 0.30. Based on local templates, these sources have extremely large infrared luminosities, comparable to that of submillimeter galaxies. Our results confirm previous indications that the rest-frame 1.6 micron stellar bump can be efficiently used to select highly obscured starforming galaxies at z~2, and that the fraction of starburst-dominated ULIRGs increases to faint 24 micron flux densities. Using local templates, we find that the observed narrow redshift distribution is due to the fact that the 24 micron detectability of PAH-rich sources peaks sharply at z = 1.9. We can analogously explain the broader redshift distribution of Spitzer-detected AGN-dominated ULIRGs based on the shapes of their SEDs. Finally, we conclude that z~2 sources with a detectable 1.6 micron stellar opacity feature lack sufficient AGN emission to veil the 7.7 micron PAH band.

We present detailed clustering analysis of a large K-band selected local galaxy sample, which is constructed from the 2MASS and the SDSS and consists of $82,486$ galaxies with $10 < K < 13.5$ and $0.01 < z < 0.1$. The two-point correlation function of the magnitude-limited sample in real space at small scales is well described by a power law $\xi(r)=(r/6.44\pm0.23)^{-1.81\pm0.02}$. The pairwise velocity dispersion is derived from the anisotropic two-point correlation function and find the dispersion $\sigma_{12}=685\pm 17\kms$ if its scale invariance is assumed, which is larger than values measured in optical bands selected galaxy samples. We further investigate the dependence of the two-point correlation function and the $\sigma_{12}$ on the $g-r$ color and the $K$-band luminosity, obtain similar results to previous works in optical bands. Comparing a mock galaxy sample with our real data indicates that the semi-analytical model can not mimic the $\sigma_{12}$ in observation albeit it can approximate the two-point correlation function within measurement uncertainties.

We present $u'g'R$ optical images taken with the MMT/Megacam and the Subaru/Suprime of the Extended Groth Strip survey. The total survey covers an area of about $\sim 1$ degree$^2$, including four sub-fields and is optimized for the study of galaxies at $z\sim3$. Our methods for photometric calibration in AB magnitudes, the limiting magnitude and the galaxy number count are described. A sample of 1642 photometrically selected candidate LBGs to an apparent $R_{AB}$ magnitude limit of 25.0 is present. The average sky surface density of our LBGs sample is $\sim$ 1.0 arcmin$^{-2}$, slightly higher than the previous finding.

We analyze a sample of galaxies chosen to have F24>0.5 mJy and satisfy a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and PAH luminosities, to which we add multi-wavelength broadband photometry. Stellar population modeling and IRS spectra together demonstrate that the double criteria used to select this sample have efficiently isolated massive star-forming galaxies at z~1.9. This is the first starburst-dominated ULIRG sample at high redshift with total infrared luminosity measured directly from FIR and millimeter photometry, and as such gives us the first accurate view of broadband SEDs for starburst galaxies at extremely high luminosity and at all wavelengths. Similar broadband data are assembled for three other galaxy samples -- local starburst galaxies, local AGN/ULIRGS, and a second 24mu-luminous z~1.9 sample dominated by AGN. L(PAH)/L(IR) for the new z~1.9 starburst sample is the highest ever seen, some three times higher than in local starbursts, whereas in AGNs this ratio is depressed below the starburst trend, often severely. Several pieces of evidence imply that AGNs exist in this starburst dominated sample, except two of which even host very strong AGN, while they still have very strong PAH emission. The ACS images show most objects have very extended morphologies in the rest-frame UV band, thus extended distribution of PAH molecules. Such an extended distribution prevents further destruction PAH molecules by central AGNs. We conclude that objects in this sample are ULIRGs powered mainly by starburst; and the total infrared luminosity density contributed by this type of objects is 0.9-2.6x 10^7 Lsol/Mpc^3.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

We present millimetre observations of a sample of 12 high redshift ultraluminous infrared galaxies (ULIRGs) in the Extended Growth Strip (EGS). These objects were initially selected on the basis of their observed mid--IR colours (0.0 < [3.6]-[4.5] < 0.4 and -0.7 < [3.6]-[8.0] < 0.5) to lie at high redshift 1.5 < z < 3, and subsequent 20-38 micron mid-IR spectroscopy confirms that they lie in a narrow redshift window centered on z=2. We detect 9/12 of the objects in our sample at high significance (>3 sigma) with a mean 1200\micron flux of <F_1200> = 1.6+/-0.1 mJy. Our millimetre photometry, combined with existing far-IR photometry from the Far-IR Deep Extragalactic Legacy (FIDEL) Survey and accurate spectroscopic redshifts, places constraints both sides of the thermal dust peak. This allows us to estimate the dust properties, including the far--IR luminosity, dust temperature, and dust mass. We find that our sample is similar to other high-z and intermediate-z ULIRGs, and local systems, but has a different dust selection function than submillimeter-selected galaxies. Finally, we use existing 20cm radio continuum imaging to test the far-IR/radio correlation at high redshift. We find that our sample is consistent with the local relation, implying little evolution. Furthermore, this suggests that our sample selection method is efficient at identifying ultraluminous, starburst--dominated systems within a very narrow redshift range centered at z~2.

We present the results of a Spitzer IRAC and MIPS 24 micron study of extended Lyman-alpha clouds (or Lyman-alpha Blobs, LABs) within the SSA22 filamentary structure at z = 3.09. We detect 6/26 LABs in all IRAC filters, four of which are also detected at 24 micron, and find good correspondence with the 850 micron measurements of Geach et al. 2005. An analysis of the rest-frame ultraviolet, optical, near- and mid-infrared colors reveals that these six systems exhibit signs of nuclear activity (AGN)and/or extreme star formation. Notably, they have properties that bridge galaxies dominated by star formation (Lyman-break galaxies; LBGs) and those with AGNs (LBGs classified as QSOs). The LAB systems not detected in all four IRAC bands, on the other hand, are, as a group, consistent with pure star forming systems, similar to the majority of the LBGs within the filament. These results indicate that the galaxies within LABs do not comprise a homogeneous population, though they are also consistent with scenarios in which the gas halos are ionized through a common mechanism such as galaxy-scale winds driven by the galaxies within them, or gravitational heating of the collapsing cloud itself.

We present high resolution submillimeter interferometric imaging of two of the brightest high-redshift submillimeter galaxies known: GN20 and AzTEC1 at 0.8 and 0.3 arcsec resolution respectively. Our data - the highest resolution submillimeter imaging of high redshift sources accomplished to date - was collected in three different array configurations: compact, extended, and very extended. We derive angular sizes of 0.6 and 1.0 arcsec for GN20 and 0.3 and 0.4 arcsec for AzTEC1 from modeling their visibility functions as a Gaussian and elliptical disk respectively. Because both sources are B-band dropouts, they likely lie within a relatively narrow redshift window around z~4, which indicates their angular extent corresponds to physical scales of 4-8 and 1.5-3 kpc respectively for the starburst region. By way of a series of simple assumptions, we find preliminary evidence that these hyperluminous starbursts - with star formation rates >1000 $M_\odot$ yr$^{-1}$ - are radiating at or close to their Eddington limit. Should future high resolution observations indicate that these two objects are typical of a population of high redshift Eddington-limited starbursts, this could have important consequences for models of star formation and feedback in extreme environments.

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the IRAC colour-colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (<100 Myr old), extreme starbursts at z>2 also occupy the same general area as AGNs in the IRAC colour-colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (>1 mag bluer in the observed [4.5]-[8.0] colour). Only about 20% of the objects overlap in the colour-colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red IR colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of IR colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Spitzer MIPS 24 micron data.

While both isothermal and NFW-based mass models for galaxy clusters are widely adopted in strong-lensing studies, they cannot easily be distinguished based solely on observed positions of arcs and arclets. We compare the magnifications predicted for giant arcs obtained from isothermal and NFW profiles, taking axially-symmetric and asymmetric mass distributions into account. We find that arc magnifications can differ strongly between the two types of density profiles even if the image morphology is well reproduced. Magnifications by lenses with NFW density profiles are usually larger than those for lenses with singular or nearly singular isothermal density profiles, unless the latter have large cores. Asymmetries play an important role. We illustrate our results with the two well-studied clusters MS 2137 and A~370. We confirm earlier results showing that both isothermal and NFW mass models can very well reproduce the observed arcs, radial arcs and other arclets. While the mass model for MS 2137 is not very well constrained, the two types of mass models produce strongly differing critical curves and caustics for A~370. We find that the NFW mass profile is preferred for A~370. We identify new arclet candidates in the field of A~370. Redshift estimates allowed by the lens model are consistent with results in the literature, where available. Three newly found counter-images are suggested to arise from an active, dust-enshrouded star-forming galaxy at $z\approx1.1$.

The Extended Groth Strip (EGS) is one of the premier fields for extragalactic deep surveys. Deep observations of the EGS with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope cover an area of 0.38 square degrees to a 50% completeness limit of 1.5 uJy at 3.6 um. The catalog comprises 57434 objects detected at 3.6 um, with 84%, 28%, and 24% also detected at 4.5, 5.8, and 8.0 um. Number counts are consistent with results from other Spitzer surveys. Color distributions show that the EGS IRAC sources comprise a mixture of populations: low-redshift star-forming galaxies, quiescent galaxies dominated by stellar emission at a range of redshifts, and high redshift galaxies and AGN.

We investigate whether the mean star formation activity of star-forming galaxies from z=0 to z=0.7 in the GOODS-S field can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in galaxies. We selected star-forming galaxies at z=0 and at z=0.7 in IR and in UV to have access to all the recent star formation. We focused on galaxies with a stellar mass ranging between 10^{10} and 10^{11} M_sun for which the results are not biased by the selections. We compared the data to chemical evolution models developed for spiral galaxies and originally built to reproduce the main characteristics of the Milky Way and nearby spirals without fine-tuning them for the present analysis. We find a shallow decrease in the specific star formation rate (SSFR) when the stellar mass increases. The evolution of the SSFR characterizing both UV and IR selected galaxies from z=0 to z=0.7 is consistent with the models built to reproduce the present spiral galaxies. There is no need to strongly modify of the physical conditions in galaxies to explain the average evolution of their star formation from z=0 to z=0.7. We use the models to predict the evolution of the star formation rate and the metallicity on a wider range of redshift and we compare these predictions with the results of semi-analytical models.

We present photometric analysis of deep mid-infrared observations obtained by Spitzer/IRAC covering the fields Q1422+2309, Q2233+1341, DSF2237a,b, HDFN, SSA22a,b and B20902+34, giving the number counts and the depths for each field. In a sample of 751 LBGs lying in those fields, 443, 448, 137 and 152 are identified at 3.6microns, 4.5microns, 5.8microns, 8.0microns IRAC bands respectively, expanding their spectral energy distribution to rest-near-infrared and revealing that LBGs display a variety of colours. Their rest-near-infrared properties are rather inhomogeneous, ranging from those that are bright in IRAC bands and exhibit [R]-[3.6] > 1.5 colours to those that are faint or not detected at all in IRAC bands with [R]-[3.6] < 1.5 colours and these two groups of LBGs are investigated. We compare the mid-IR colours of the LBGs with the colours of star-forming galaxies and we find that LBGs have colours consistent with star-foming galaxies at z~3. The properties of the LBGs detected in the 8microns IRAC band (rest frame K-band) are examined separately, showing that they exhibit redder [R]-[3.6] colours than the rest of the population and that IRAC 8microns band can be used as a diagnostic tool, to separate AGN dominated objects from normal star-forming galaxies at z~3

We examine the infrared properties of 43 high redshift (0.1<z<1.2), infrared-luminous galaxies in the Extended Groth Strip (EGS), selected by a deep 70 micron survey with the Multiband Imaging Photometer on Spitzer (MIPS). In addition and with reference to starburst-type Spectral Energy Distributions (SEDs), we derive a set of equations for estimating the total infrared luminosity (L_IR) in the range 8-1000 microns using photometry from at least one MIPS band. 42 out of 43 of our sources' optical/infrared SEDs (lambda_observed < 160 microns) are starburst-type, with only one object displaying a prominent power-law near-infrared continuum. For a quantitative analysis, models of radiation transfer in dusty media are fit onto the infrared photometry, revealing that the majority of galaxies are represented by high extinction, A_v>35 and for a large fraction (~50 per cent) the SED turns over into the Rayleigh-Jeans regime at wavelengths longward of 90 microns. For comparison, we also fit semi-empirical templates based on local galaxy data, however, these underestimate the far-infrared SED shape by a factor of at least 2 and in extreme cases up to 10 for the majority (~70 per cent) of the sources. Further investigation of SED characteristics reveals that the mid-infrared (70/24 microns) continuum slope is decoupled from various galaxy properties such as the total infrared luminosity and far-infrared peak, quantified by the L_160/L_70 ratio. In view of these results, we propose that these high-redshift galaxies have different properties to their local counterparts, in the sense that large amounts of dust cause heavy obscuration and are responsible for an additional cold emissive component, appearing as a far-infrared excess in their SEDs.

We investigate the environment of infrared luminous galaxies (L$_{IR}$[8-1000 $\mu $m$] >10^{11}$L$_{\sun}$). We focus on the redshift range 0.7 $\leq$ z $\leq$ 1, where these galaxies dominate the star formation activity and play a significant role in galaxy evolution. We employ MIPS 24$\mu$m data to identify infrared galaxies in the Extended Groth Strip (EGS). We use a local density indicator to probe the environment on few Mpc scales and a group member catalog, both of which make use of the DEEP2 spectroscopic redshift catalog, to quantify the environment of these galaxies. We find that the local environment of LIRGs and ULIRGs is intermediate between that of blue and red galaxies. LIRGs and ULIRGs avoid underdense environments and inhabit local environments that are more dense on average than those of other DEEP2 galaxies at similar redshifts. However, when the comparison sample of the non-IR DEEP2 galaxies is restricted to have the same range of stellar mass, color, or luminosity as the IR--galaxies, there is no longer any significant difference in environment; the IR-galaxies follow the same trends in the color-environment and luminosity-environment relations observed at z$\sim$1. We also find that about 30% of the LIRGs and ULIRGs belong to groups, associated with a minimum dark matter halo of 6$\times10^{12}$M$_{\odot}$h$^{-1}$. The group members constitute 20 % of the sources responsible for the IR star formation rate density and comoving energy density at z$\sim$1.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

We present the {\it Spitzer} Space Telescope InfraRed Array Camera (IRAC) and Multiband Imaging Photometer (MIPS) observations of the elliptical galaxy NGC 315. After removal of the host galaxy's stellar emission, we detected for the first time an infrared-red nucleus in NGC 315. We measured the spectral energy distribution (SED) for this active nucleus with wavelength range covering from radio to X-ray, and obtained the bolometric luminosity of $\rm L_{bol} \approx 1.9 \times 10^{43} ergs s^{-1}$, corresponding to an extremely low Eddington ratio (L/L$_{\rm Edd}$) of 4.97 $\times$ 10$^{-4}$. Our results confirm that the physical nature of the nucleus of NGC 315 is a low-luminosity AGN, consistent with the recent optical and {\it Chandra} X-ray observations.

We investigate the stellar populations of a sample of 162 Lyman-alpha emitting galaxies (LAEs) at z = 3.1 in the Extended Chandra Deep Field South, using deep Spitzer IRAC data available from the GOODS and SIMPLE surveys to derive reliable stellar population estimates. We divide the LAEs according to their rest-frame near-IR luminosities into IRAC-detected and IRAC-undetected samples. About 70% of the LAEs are undetected in 3.6 micron down to [3.6] = 25.2 AB. Stacking analysis reveals that the average stellar population of the IRAC-undetected sample has an age of ~ 200 Myr and a mass of ~ 3x10^8 solar masses, consistent with the expectation that LAEs are mostly young and low-mass galaxies. On the other hand, the IRAC-detected LAEs are on average significantly older and more massive, with an average age > 1 Gyr and mass ~ 10^10 solar masses. Comparing the IRAC colors and magnitudes of the LAEs to z ~ 3 Lyman break galaxies (LBGs) shows that the IRAC-detected LAEs lie at the faint blue end of the LBG color-magnitude distribution, suggesting that IRAC-detected LAEs may be the low mass extension of the LBG population. We also present tentative evidence for a small fraction (~ 5%) of obscured AGN within the LAE sample. Our results suggest that LAEs posses a wide range of ages and masses. Additionally, the presence of evolved stellar populations inside LAEs suggests that the Lyman-alpha luminous phase of galaxies may either be a long-lasting or recurring phenomenon.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host AGN. The LAEs exhibit a moderate spatial correlation length of r_0=3.6+0.8-1.0 Mpc, corresponding to a bias factor b=1.7+0.3-0.4, which implies median dark matter halo masses of log10(M_med) = 10.9+0.5-0.9 M_sun. Comparing the number density of LAEs, (1.5+-0.3)E-3/Mpc^3, with the number density of these halos finds a mean halo occupation ~1-10%. The evolution of galaxy bias with redshift implies that most z=3.1 LAEs evolve into present-day galaxies with L<2.5L*, whereas other z>3 galaxy populations typically evolve into more massive galaxies. Halo merger trees show that z=0 descendants occupy halos with a wide range of masses, with a median descendant mass close to that of L*. Only 30% of LAEs have sufficient stellar mass (>~3E9 M_sun) to yield detections in deep Spitzer-IRAC imaging. A two-population SED fit to the stacked UBVRIzJK+[3.6,4.5,5.6,8.0]micron fluxes of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0+0.6-0.4 E9 M_sun), moderate star formation rate (2+-1 M_sun/yr), a young component age of 20+30-10 Myr, and little dust (A_V<0.2). The best fit model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed.

We present the quantitative rest-frame B morphological evolution and galaxy merger fractions at 0.2 < z < 1.2 as observed by the All-wavelength Extended Groth Strip International Survey (AEGIS). We use the Gini coefficent and M_20 to identify major mergers and classify galaxy morphology for a volume-limited sample of 3009 galaxies brighter than 0.4 L_B^*, assuming pure luminosity evolution of 1.3 M_B per unit redshift. We find that the merger fraction remains roughly constant at 10 +/- 2% for 0.2 < z < 1.2. The fraction of E/S0/Sa increases from 21+/- 3% at z ~ 1.1 to 44 +/- 9% at z ~ 0.3, while the fraction of Sb-Ir decreases from 64 +/- 6% at z ~ 1.1 to 47 +/- 9% at z ~ 0.3. The majority of z < 1.2 Spitzer MIPS 24 micron sources with L(IR) > 10^11 L_sun are disk galaxies, and only ~ 15% are classified as major merger candidates. Edge-on and dusty disk galaxies (Sb-Ir) are almost a third of the red sequence at z ~ 1.1, while E/S0/Sa makeup over 90% of the red sequence at z ~ 0.3. Approximately 2% of our full sample are red mergers. We conclude (1) the galaxy merger rate does not evolve strongly between 0.2 < z < 1.2; (2) the decrease in the volume-averaged star-formation rate density since z ~ 1 is a result of declining star-formation in disk galaxies rather than a disappearing population of major mergers; (3) the build-up of the red sequence at z < 1 can be explained by a doubling in the number of spheroidal galaxies since z ~ 1.2.

Ultraviolet (UV) galaxies have been selected from GALEX. The presence of a FUV-dropout in their spectral energy distributions proved to be a very complete (83.3%) but not very efficient (21.4%) tool for identifying Lyman Break Galaxies (LBGs) at z~1. We divide the LBG sample into two sub-classes: red LBGs (RLBGs) detected at 24 micron which are mainly Luminous IR Galaxies (LIRGs) and blue LBGs (BLBGs) undetected at 24 microns down to 83 microJy. Two of the RLBGs are also detected at 70 micron. The median SED of the RLBGs is similar (above lambda~1 micron) to the dusty starburst HR10. However, unlike local (U)LIRGs, RLBGs are UV bright objects. We suggest that these objects contain a large amount of dust but that some bare stellar populations are also directly visible. The median SED of the BLBGs is consistent with their containing the same stellar population as the RLBGs but with a lower dust content. The luminosity function of our LBG sample at z~1 is similar to the luminosity function of NUV-selected galaxies at the same redshift. The integrated luminosity densities of z~1 LBGs and NUV-selected galaxies are very consistent. We show that star formation rates (SFRs) estimated from UV measurements and corrected using the IRX-beta method provide average total SFR_TOT in agreement with SFR_UV + SFR_dust. However, IRX-beta-based SFR_TOT shows a large dispersion. Summing up the detected UV (1150A rest-frame) and IR-based star formation rates of the detected objects, we find that only one third of the total (i.e. UV + dust) LBG SFR resides in BLBGs and two thirds in RLBGs, even though most LBGs at z~1 are BLBGs. On the other hand, the total SFR of LBGs accounts for only 11% of the total SFR at z~1. Finally, we observe a regular decrease of L_TIR / L_FUV from z=0 to z~2 for UV-selected samples.

We present the results of a study on the properties and evolution of massive (M_* > 10^11 M_0) galaxies at z~0.4 - 2 utilising Keck spectroscopy, near-Infrared Palomar imaging, and Hubble, Chandra, and Spitzer data covering fields targeted by the DEEP2 galaxy spectroscopic survey. Our sample is K band selected based on wide-area NIR imaging from the Palomar Observatory Wide-Field Infrared Survey, which covers 1.53 deg^2 to K_s,vega~20.5. Our major findings include: (i) statistically the mass and number densities of M_* > 10^11 M_0 galaxies show little evolution between z = 0 - 1, and from z ~ 0 - 2 for M_* > 10^11.5 M_0 galaxies. (ii) Using Hubble ACS imaging, we find that M_* > 10^11 selected galaxies show a nearly constant elliptical fraction of ~70-90% at all redshifts. The remaining objects are peculiars possibly undergoing mergers at z > 0.8, while spirals dominate the remainder at lower redshifts. (iii) We find that only a fraction (~60%) of massive galaxies with M_* > 10^11 M_0 are on the red-sequence at z~1.4, while nearly 100% evolve onto it by z~0.4. (iv) By utilising Spitzer MIPS imaging and [OII] line fluxes we argue that M_* > 10^11.5 galaxies have a steeply declining star formation rate density ~(1+z)^6. By examining the contribution of star formation to the evolution of the mass function, as well as the merger history through the CAS parameters, we determine that M_* >10^11 M_0 galaxies undergo on average 0.9^+0.7_-0.5 major mergers at 0.4 < z < 1.4. (v) A high (5%) fraction of all M_* > 10^11 M_0 galaxies are X-ray emitters. Roughly half of these are morphologically distorted ellipticals or peculiars. We compare our results with the Millennium simulation, finding that the number and mass densities of M_* > 10^11.5 M_0 galaxies are under predicted by a factor of > 100.

We have used the Submillimeter Array to image a flux limited sample of seven submillimeter galaxies, selected by the AzTEC camera on the JCMT at 1.1 mm, in the COSMOS field at 890um with 2" resolution. All of the sources - two radio-bright and five radio-dim - are detected as single point-sources at high significance (> 6\sigma), with positions accurate to 0.2" that enable counterpart identification at other wavelengths observed with similarly high angular resolution. All seven have IRAC counterparts, but only two have secure counterparts in deep HST/ACS imaging. As compared to the two radio-bright sources in the sample, and those in previous studies, the five radio-dim sources in the sample (1) have systematically higher submillimeter-to-radio flux ratios, (2) have lower IRAC 3.6-8.0um fluxes, and (3) are not detected at 24um. These properties, combined with size constraints at 890um (\theta < 1.2"), suggest that the radio-dim submillimeter galaxies represent a population of very dusty starbursts, with physical scales similar to local ultraluminous infrared galaxies, and an average redshift higher than radio-bright sources.

We present the rest-frame ultraviolet through near infrared spectral energy distribution for an interacting Lyman break galaxy at a redshift z=4.42, the highest redshift merging system known with clearly resolved tidal features. The two objects in this system - HDF-G4 and its previously unidentified companion - are both B_{435} band dropouts, have similar V_{606}-i_{775} and i_{775}-z_{850} colors, and are separated by 1", which at z=4.42 corresponds to 7 kpc projected nuclear separation; all indicative of an interacting system. Fits to stellar population models indicate a stellar mass of M_\star = 2.6\times 10^{10} M_\odot, age of \tau_\star = 720 My, and exponential star formation history with an e-folding time \tau_0 = 440 My. Using these derived stellar populations as constraints, we model the HDF-G4 system using hydrodynamical simulations, and find that it will likely evolve into a quasar by z\sim3.5, and a quiescent, compact spheroid by z\sim 2.5 similar to those observed at z > 2. And, the existence of such an object supports galaxy formation models in which major mergers drive the high redshift buildup of spheroids and black holes.

We study the infrared (IR) properties of high-redshift galaxies using deep Spitzer 24, 70, and 160 micron data. Our primary interest is to improve the constraints on the total IR luminosities, L(IR), of these galaxies. We combine the Spitzer data in the southern Extended Chandra Deep Field with a K-band-selected galaxy sample and photometric redshifts from the Multiwavelength Survey by Yale-Chile. We used a stacking analysis to measure the average 70 and 160 micron flux densities of 1.5 < z < 2.5 galaxies as a function of 24 micron flux density, X-ray activity, and rest-frame near-IR color. Galaxies with 1.5 < z < 2.5 and S(24)=53-250 micro-Jy have L(IR) derived from their average 24-160 micron flux densities within factors of 2-3 of those derived from the 24 micron flux densities only. However, L(IR) derived from the average 24-160 micron flux densities for galaxies with S(24) > 250 micro-Jy and 1.5 < z < 2.5 are lower than those derived using only the 24 micron flux density by factors of 2-10. Galaxies with S(24) > 250 micro-Jy have S(70)/S(24) flux ratios comparable to sources with X-ray detections or red rest-frame IR colors, suggesting that warm dust possibly heated by AGN may contribute to the high 24 micron emission. Based on the average 24-160 micron flux densities, nearly all 24 micron-selected galaxies at 1.5 < z < 2.5 have L(IR) < 6 x 10^12 solar luminosities, which if attributed to star formation corresponds to < 1000 solar masses per year. This suggests that high redshift galaxies may have similar star formation efficiencies and feedback processes as local analogs. Objects with L(IR) > 6 x 10^12 solar luminosities are quite rare, with a surface density ~ 30 +/- 10 per sq. deg, corresponding to ~ 2 +/- 1 x 10^-6 Mpc^-3 over 1.5 < z < 2.5.

A large position angle misalignment between the stellar bar and the distribution of dust in the late-type barred spiral NGC 3488 was discovered, using mid-infrared images from the Spitzer Space Telescope and optical images from the Sloan Digital Sky Survey (SDSS). The angle between the bar and dust patterns was measured to be 25+-2deg, larger than most of the misalignments found previously in barred systems based on Ha or HI/CO observations. The stellar bar is bright at optical and 3.6um, while the dust pattern is more prominent in the 8um band but also shows up in the SDSS u and g-band images, suggesting a rich interstellar medium environment harboring ongoing star formation. This angular misalignment is unlikely to have been caused by spontaneous bar formation. We suggest that the stellar bar and the dust pattern may have different formation histories, and that the large misalignment was triggered by a tidal interaction with a small companion. A statistical analysis of a large sample of nearby galaxies with archival Spitzer data indicates that bar structure such as that seen in NGC 3488 is quite rare in the local Universe.

A SST survey in the NOAO Deep-Wide Field in Bo\"otes provides a complete, 8-micron-selected sample of galaxies to a limiting (Vega) magnitude of 13.5. In the 6.88 deg$^2$ field sampled, 79% of the 4867 galaxies have spectroscopic redshifts, allowing an accurate determination of the local (z<0.3) galaxy luminosity function. Stellar and dust emission can be separated on the basis of observed galaxy colors. Dust emission (mostly PAH) accounts for 80% of the 8 micron luminosity, stellar photospheres account for 19%, and AGN emission accounts for roughly 1 %. A sub-sample of the 8 micron-selected galaxies have blue, early-type colors, but even most of these have significant PAH emission. The luminosity functions for the total 8 micron luminosity and for the dust emission alone are both well fit by Schechter functions. For the 8 micron luminosity function, the characteristic luminosity is \nu L_{\nu}^*(8.0 \micron) = 1.8 \times 10^{10}$ \Lsun while for the dust emission alone it is 1.6 x 10^{10}$ \Lsun\null. The average 8 \micron luminosity density at z<0.3 is 3.1 x 10^7 \Lsun Mpc^{-3}, and the average luminosity density from dust alone is 2.5 x 10^7 \Lsun Mpc^{-3}. This luminos ity arises predominantly from galaxies with 8 \micron luminosities ($\nu L_{\nu}$) between $2\times 10^9$ and $2 x 10^{10}$ \Lsun, i.e., normal galaxies, not LIRGs or ULIRGs.

Using data from the DEEP2 Galaxy Redshift Survey and HST/ACS imaging in the Extended Groth Strip, we select nearly 100 interacting galaxy systems including kinematic close pairs and morphologically identified merging galaxies. Spitzer MIPS 24 micron fluxes of these systems reflect the current dusty star formation activity, and at a fixed stellar mass (M_{*}) the median infrared luminosity (L_{IR}) among merging galaxies and close pairs of blue galaxies is twice (1.9 +/- 0.4) that of control pairs drawn from isolated blue galaxies. Enhancement declines with galaxy separation, being strongest in close pairs and mergers and weaker in wide pairs compared to the control sample. At z ~ 0.9, 7.1% +/- 4.3% of massive interacting galaxies (M_{*} > 2*10^{10} M_{solar}) are found to be ULIRGs, compared to 2.6% +/- 0.7% in the control sample. The large spread of IR luminosity to stellar mass ratio among interacting galaxies suggests that this enhancement may depend on the merger stage as well as other as yet unidentified factors (e.g., galaxy structure, mass ratio, orbital characteristics, presence of AGN or bar). The contribution of interacting systems to the total IR luminosity density is moderate (<= 36 %).

We analyze star formation (SF) as a function of stellar mass (M*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS). For 2905 field galaxies, complete to 10^10(10^10.8) Msun at z<0.7(1), with Keck spectroscopic redshifts out to z=1.1, we compile SF rates (SFR) from emission lines, GALEX, and Spitzer MIPS 24 micron photometry, optical-NIR M* measurements, and HST morphologies. Galaxies with reliable signs of SF form a distinct "main sequence (MS)", with a limited range of SFR at a given M* and z (1 sigma < +-0.3 dex), and log(SFR) approximately proportional to log(M*). The range of log(SFR) remains constant to z>1, while the MS as a whole moves to higher SFR as z increases. The range of SFR along the MS constrains the amplitude of episodic variations of SF, and the effect of mergers on SFR. Typical galaxies spend ~67(95)% of their lifetime since z=1 within a factor of <~ 2(4) of their average SFR at a given M* and z. The dominant mode of the evolution of SF since z~1 is apparently a gradual decline of the average SFR in most individual galaxies, not a decreasing frequency of starburst episodes, or a decreasing factor by which SFR are enhanced in starbursts. LIRGs at z~1 seem to mostly reflect the high SFR typical for massive galaxies at that epoch. The smooth MS may reflect that the same set of few physical processes governs star formation prior to additional quenching processes. A gradual process like gas exhaustion may play a dominant role.

Based on the catalog of merging galaxies in the Early Data Release (EDR) of the Sloan Digital Sky Survey (SDSS), the differences of star formation history between merging galaxies and field galaxies are studied statistically by means of three spectroscopic indicators the 4000-\r{A} break strength, the Balmer absorption-line index, and the specific star formation rate. It is found that for early-type merging galaxies the interactions will not induce significant enhancement of the star-formation activity because of its stability and lack of cool gas. On the other hand, late-type merging galaxies always in general display more active star formation than field galaxies on different timescales within about 1Gyr. We also conclude that the mean stellar ages of late-type merging galaxies are younger than those of late-type field galaxies.

We simulate the bulge formation in very late-type dwarf galaxies from circumnuclear super star clusters (SSCs) moving in a responding cuspy dark matter halo (DMH). The simulations show that (1) the response of DMH to sinking of SSCs is detectable only in the region interior to about 200 pc. The mean logarithmic slope of the responding DM density profile over that area displays two different phases: the very early descent followed by ascent till approaching to 1.2 at the age of 2 Gyrs. (2) the detectable feedbacks of the DMH response on the bulge formation turned out to be very small, in the sense that the formed bulges and their paired nuclear cusps in the fixed and the responding DMH are basically the same, both are consistent with $HST$ observations. (3) the yielded mass correlation of bulges to their nuclear (stellar) cusps and the time evolution of cusps' mass are accordance with recent findings on relevant relations. In combination with the consistent effective radii of nuclear cusps with observed quantities of nuclear clusters, we believe that the bulge formation scenario that we proposed could be a very promising mechanism to form nuclear clusters.

We use deep GALEX images of CDFS in UV to define the first large sample of 420 Lyman Break Galaxies at z~1. We use a PSF fitting to estimate UV magnitudes on these deep crowded images. Deep Spitzer IRAC and MIPS provide the first detection of a large sample of Lyman Break Galaxies in the mid- to far-infrared range. We are therefore able to study and compare the UV and TIR emission of Lyman Break Galaxies. We find that about 15% of the LBG sample are strong emitters at 24 microns (Red LBGs). Most of them are Luminous IR Galaxies (LIRGs) while the rest (Blue LBGs) are undetected at the 83 microJy level of MIPS GTO image. We find that Blue LBGs have a Spectral Energy Distribution similar to high redshift Lyman Break Galaxies. Finally, the dust-to-FUV ratio of this sample is compared with similar ratios at z=1 and z~2. This work suggests an evolution (decrease) of the dust-to-FUV ratio with the redshift.

We have measured the contribution of submillimeter and mid-infrared sources to the extragalactic background radiation at 70 and 160um. Specifically, we have stacked flux in 70 and 160um Spitzer Space Telescope (Spitzer) observations of the Canada-UK Deep Sub-millimeter Survey 14h field at the positions of 850um sources detected by SCUBA and also 8 and 24um sources detected by Spitzer. We find that per source, the SCUBA galaxies are the strongest and the 8um sources the weakest contributors to the background flux at both 70 and 160um. Our estimate of the contribution of the SCUBA sources is higher than previous estimates. However, expressed as a total contribution, the full 8um source catalogue accounts for twice the total 24um source contribution and ~10 times the total SCUBA source contribution. The 8um sources account for the majority of the background radiation at 160um with a flux of 0.87+/-0.16 MJy/sr and at least a third at 70um with a flux of 0.103+/-0.019 MJy/sr. These measurements are consistent with current lower limits on the background at 70 and 160um. Finally, we have investigated the 70 and 160um emission from the 8 and 24um sources as a function of redshift. We find that the average 70um flux per 24um source and the average 160um flux per 8 and 24um source is constant over all redshifts, up to z~4. In contrast, the low-redshift half (z<1) of the of 8um sample contributes approximately four times the total 70um flux of the high-redshift half. These trends can be explained by a single non-evolving SED.

We present 0.5 -160 micron Spectral Energy Distributions (SEDs) of galaxies, detected at 70microns with the Multiband Imaging Photometer for Spitzer (MIPS), using broadband imaging data from Spitzer and ground-based telescopes. Spectroscopic redshifts, in the range 0.2<z<1.5, have been measured as part of the Deep Extragalactic Evolutionary Probe2 (DEEP2) project. Based on the SEDs we explore the nature and physical properties of the sources. Using the optical spectra we derive Hbeta and [OII]-based Star Formation Rates (SFR) which are 10-100 times lower than SFR estimates based on IR and radio. The median offset in SFR between optical and IR is reduced by a factor of ~3 when we apply a typical extinction corrections. We investigate mid-to-far infrared correlations for low redshift (>0.5) and high redshift (0.5<z<1.2) bins. Using this unique ``far-infrared'' selected sample we derive an empirical mid to far-infrared relationship that can be used to estimate the infrared energy budget of galaxies in the high-redshift universe. Our sample can be used as a template to translate far-infrared luminosities into bolometric luminosities for high redshift objects.

We present a study of three Lyman-alpha emitting galaxies (LAEs), selected via a narrow-band survey in the GOODS northern field, and spectroscopically confirmed to have redshifts of z ~ 5.65. Using HST ACS and Spitzer IRAC data, we constrain the rest-frame UV-to-optical spectral energy distributions (SEDs) of the galaxies. Fitting stellar population synthesis models to the observed SEDs, we find best-fit stellar populations with masses between ~ 10^9 - 10^10 M_sun and ages between ~ 5 - 100 Myr, assuming a simple starburst star formation history. However, stellar populations as old as 700 Myr are admissible if a constant star formation rate model is considered. Very deep near-IR observations may help to narrow the range of allowed models by providing extra constraints on the rest-frame UV spectral slope. Our narrow-band selected objects and other IRAC-detected z ~ 6 i'-dropout galaxies have similar 3.6 um magnitudes and z' - [3.6] colors, suggesting that they posses stellar populations of similar masses and ages. This similarity may be the result of a selection bias, since the IRAC-detected LAEs and i'-dropouts probably only sample the bright end of the luminosity function. On the other hand, our LAEs have blue i' - z' colors compared to the i'-dropouts, and would have been missed by the i'-dropout selection criterion. A better understanding of the overlap between the LAE and the i'-dropout populations is necessary in order to constrain the properties of the overall high-redshift galaxy population, such as the total stellar mass density at z ~ 6.

We describe AEGIS20 - a radio survey of the Extended Groth Strip (EGS) conducted with the Very Large Array (VLA) at 1.4GHz. The resulting catalog contains 1,123 emitters and is sensitive to ultraluminous starbursts to z<=1.3, well matched to the redshift range of the DEEP2 spectroscopic survey in this region. We use stacking techniques to explore the uJy-level emission from a variety of galaxy populations selected via conventional criteria - Lyman-break galaxies (LBGs), distant red galaxies (DRGs), UV-selected galaxies and extremely red objects (EROs) - determining their properties as a function of color, magnitude and redshift and their extinction-free contributions to the history of star formation. We confirm the familiar pattern that the star-formation-rate (SFR) density, SFRD, rises by at least ~5x from z=0-1, although we note highly discrepant UV- and radio-based SFR estimates. Our radio-based SFRs become more difficult to interpret at z>1 where correcting for contamination by radio-loud active galactic nuclei (AGN) comes at the price of rejecting luminous starbursts. Whilst stacking radio images is a useful technique, accurate radio-based SFRs for z>>1 galaxies require precise redshifts and extraordinarily high-fidelity radio data to identify and remove accretion-related emission.

We analyze very deep HST, VLT and Spitzer photometry of galaxies at 2<z<3.5 in the Hubble Deep Field South. The sample is selected from the deepest public K-band imaging currently available. We show that the rest-frame U-V vs V-J color-color diagram is a powerful diagnostic of the stellar populations of distant galaxies. Galaxies with red rest-frame U-V colors are generally red in rest-frame V-J as well. However, at a given U-V color a range in V-J colors exists, and we show that this allows us to distinguish young, dusty galaxies from old, passively evolving galaxies. We quantify the effects of IRAC photometry on estimates of masses, ages, and the dust content of z>2 galaxies. The estimated distributions of these properties do not change significantly when adding IRAC data to the UBVIJHK photometry. However, for individual galaxies the addition of IRAC can improve the constraints on the stellar populations, especially for red galaxies: uncertainties in stellar mass decrease by a factor of 2.7 for red (U-V > 1) galaxies, but only by a factor of 1.3 for blue (U-V < 1) galaxies. We find a similar color-dependence of the improvement for estimates of age and dust extinction. In addition, the improvement from adding IRAC depends on the availability of full near-infrared JHK coverage; if only K-band were available, the mass uncertainties of blue galaxies would decrease by a more substantial factor 1.9. Finally, we find that a trend of galaxy color with stellar mass is already present at z>2. The most massive galaxies at high redshift have red rest-frame U-V colors compared to lower mass galaxies even when allowing for complex star formation histories.

As part of the "All Wavelength Extended Groth Strip International Survey" (AEGIS), we describe the panchromatic characterization of an X-ray luminous active galactic nucleus (AGN) in a merging galaxy at z=1.15. This object is detected at infrared (8mic, 24mic, 70mic, 160mic), submillimeter (850mic) and radio wavelengths, from which we derive a bolometric luminosity L_bol ~ 9x10^12 Lsol. We find that the AGN clearly dominates the hot dust emission below 40mic but its total energetic power inferred from the hard X-rays is substantially less than the bolometric output of the system. About 50% of the infrared luminosity is indeed produced by a cold dust component that probably originates from enshrouded star formation in the host galaxy. In the context of a coeval growth of stellar bulges and massive black holes, this source might represent a ``transition'' object sharing properties with both quasars and luminous starbursts. Study of such composite galaxies will help address how the star formation and disk-accretion phenomena may have regulated each other at high redshift and how this coordination may have participated to the build-up of the relationship observed locally between the masses of black holes and stellar spheroids.

The All-wavelength Extended Groth Strip International Survey (AEGIS) team presents broad-band spectral energy distributions (SEDs), from X-ray to radio wavelengths, for 71 galaxies spanning the redshift range 0.55-1.16 (<z>~0.7). Galaxies with secure redshifts were selected from a small (22 arcminute-square) sub-section of the Keck/DEIMOS galaxy redshift survey in the Extended Groth Strip field that has also been targeted for deep panchromatic imaging by Chandra (X-ray), GALEX (ultraviolet), Canada-France-Hawaii Telescope (optical), Hubble Space Telescope (optical/near infrared), Palomar Observatory (near infrared), Spitzer (mid/far infrared), and the Very Large Array (radio.) The absolute magnitude of the typical galaxy in our sample is M_B=-19.82. The ultraviolet to mid-infrared portion of their spectral energy distributions (SEDs) are found to be bracketed by two stellar-only model SEDs: an early burst followed by passive evolution and a constant star-formation rate since early times; this suggests that few of these galaxies are undergoing major starbursts. Approximately half the galaxies show a mid- to far-infrared excess relative to the model SEDs, consistent with thermal emission from interstellar dust. Two objects have power-law SEDs, indicating that they are dominated by active galactic nuclei; both are detected in X-rays. The galaxies are grouped by rest-frame color,quantitative optical morphology, and [OII] emission line strength (possible indicator of star formation). On average, the panchromatic SEDs of the galaxies, from the ultraviolet to the infrared, follow expected trends: redder SEDs are associated with red U-B, early-type morphology, and low [OII] emission, and vice versa for blue SEDs.

We report the detection of rest--frame 6.2 and 7.7 \micron emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01. This is currently the highest redshift galaxy where these PAH emission features have been detected. The total infrared luminosity inferred from the MIPS 24 \micron and radio flux density is 2$\times10^{13}$ L$_{\odot}$, which qualifies this object as a so--called hyperluminous infrared galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG has very strong 6.2 and 7.7 \micron PAH emission. We argue that intense star formation dominates the infrared emission of this source, although we cannot rule out the presence of a deeply obscured AGN. This LBG appears to be a distorted system in the HST ACS F606W and F814W images, possibly indicating that a significant merger or interaction is driving the large IR luminosity.

We study 87 Extremely Red Objects (EROs), selected both to have color redder than R-[3.6] = 4.0, and to have confirmed spectroscopic redshifts. Together, these two constraints result in this sample populating a fairly narrow redshift range at 0.76 < z < 1.42. The key new ingredient included here is deep Spitzer Space Telescope InfraRed Array Camera (IRAC) data. Based on [3.6]-[8.0] color, we demonstrate that it is possible to classify EROs into early-type, dusty starburst, or power-law (AGN) types. We present ultraviolet to mid-infrared spectral energy distributions (SEDs) and Advanced Camera for Surveys (ACS) images, both of which support our simple IRAC color classification.

A new 6-cm survey of almost 0.6 square degrees to a limit of 0.55-mJy/beam (10-sigma) finds 37 isolated radio sources and 7 radio source pairs (not necessarily physical companions). IRAC counterparts are identified for at least 92% of the radio sources within the area of deep IRAC coverage, which includes 31 isolated sources and 6 pairs. This contrasts with an identification rate of <74% to R<23.95 in visible light. Eight of the IRAC galaxies have power law spectral energy distributions, implying that the mid-infrared emission comes from a powerful AGN. The remaining 26 IRAC galaxies show stellar emission in the mid-infrared, probably in most of these galaxies because the stellar emission is bright enough to outshine an underlying AGN. The infrared colors suggest that the majority of these galaxies are bulge-dominated and have redshifts between approximately 0.5 and 1.0. Visible spectra from the DEEP2 redshift survey, available for 11 galaxies, are consistent with this suggestion. The IRAC galaxies fall into two distinct groups in a color-magnitude diagram, one group (the "stripe") includes all the AGN. The other group (the "blue clump") has blue 3.6 to 8 micron colors and a small range of 8 micron magnitudes. This separation should be useful in classifying galaxies found in other radio surveys.

We determine the importance of redshift-dependent systematic effects in the determination of stellar masses from broad band spectral energy distributions (SEDs), using high quality kinematic and photometric data of early-type galaxies at z~1 and z~0. We find that photometric masses of z~1 galaxies can be systematically different, by up to a factor of 2, from photometric masses of z~0 galaxies with the same dynamical mass. The magnitude of this bias depends on the choice of stellar population synthesis model and the rest-frame wavelength range used in the fits. The best result, i.e., without significant bias, is obtained when rest-frame optical SEDs are fitted with models from Bruzual&Charlot (2003). When the SEDs are extended to the rest-frame near-IR, a bias is introduced: photometric masses of the z~1 galaxies increase by a factor of 2 relative to the photometric masses of the z~0 galaxies. When we use the Maraston (2005) models, the photometric masses of the z~1 galaxies are low relative to the photometric masses of the z~0 galaxies by a factor of ~1.8. This offset occurs both for fits based on rest-frame optical SEDs, and fits based on rest-frame optical+near-IR SEDs. The results indicate that model uncertainties produce uncertainties as high as a factor of 2.5 in mass estimates from rest-frame near-IR photometry, independent of uncertainties due to unknown star formation histories.

Accurate photometric redshifts are calculated for nearly 200,000 galaxies to a 4.5 micron flux limit of ~13 uJy in the 8.5 deg^2 Spitzer/IRAC Shallow survey. Using a hybrid photometric redshift algorithm incorporating both neural-net and template-fitting techniques, calibrated with over 15,000 spectroscopic redshifts, a redshift accuracy of \sigma = 0.06(1+z) is achieved for 95% of galaxies at 0<z<1.5. The accuracy is \sigma = 0.12(1+z) for 95% of AGN at 0<z<3. Redshift probability functions, central to several ongoing studies of the galaxy population, are computed for the full sample. We demonstrate that these functions accurately represent the true redshift probability density, allowing the calculation of valid confidence intervals for all objects. These probability functions have already been used to successfully identify a population of Spitzer-selected high redshift (z>1) galaxy clusters. We present one such spectroscopically confirmed cluster at <z>=1.24, ISCS J1434.5+3427. Finally, we present a measurement of the 4.5 micron-selected galaxy redshift distribution.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

We describe the spectral energy distributions (SEDs) of Lyman Break Galaxies (LBGs) at z~3 using deep mid-infrared and optical observations of the Extended Groth Strip, obtained with IRAC and MIPS on board Spitzer and from the ground, respectively. We focus on LBGs with detections at all four IRAC bands, in particular the 26 galaxies with IRAC 8 micron band (rest--frame K-band) detections. We use stellar population synthesis models and probe the stellar content of these galaxies. Based on best--fit continuous star-formation models we derive estimates of the stellar mass for these LBGs. As in previous studies, we find that a fraction of LBGs have very red colors and large estimated stellar masses (M_stellar > 5 x 10 ^{10} M_solar): the present Spitzer data allow us, for the first time, to study these massive LBGs in detail. We discuss the link between these LBGs and submm-luminous galaxies.We find that the number density of these massive LBGs at high redshift is higher than predicted by current semi-analytic models of galaxy evolution.

We identify 17 possible 8.0 micron-selected counterparts to the submillimeter galaxies in the CUDSS 14-hour field, derived from deep imaging carried out with the IRAC and MIPS instruments aboard the Spitzer Space Telescope. Ten of the 17 counterparts are not the same as those previously identified at shorter wavelengths. We argue that 8.0 micron selection offers a better means for identifying counterparts to submillimeter galaxies than near-infrared or optical selection. Based on the panchromatic SEDs, most counterparts appear to be powered by ongoing star formation. Power-law fits to the SEDs suggest that five objects in the 8.0 micron-selected counterpart sample harbor dominant AGNs; a sixth object is identified as a possible AGN. The 3.6 to 8.0 micron colors of the infrared-selected counterparts are significantly redder than the general IRAC galaxy population in the CUDSS 14-hour field.

We have studied the infrared properties of the brightest cluster galaxies (BCGs) located in the cores of X-ray-luminous clusters at 0.15 < z < 0.35. The majority of the BCGs are not particularly infrared-luminous compared with other massive early-type galaxies, suggesting that the cluster environment has little influence on the infrared luminosities of the BCGs. The exceptions, however, are the BCGs in the three X-ray-brightest clusters in the sample, A1835, Z3146, and A2390. These BCGs have a prominent far-infrared peak in their spectral energy distributions (SEDs), and two of them (those in A1835 and Z3146) can be classified as luminous infrared galaxies (LIRGs: L_{IR} > 10^{11} L_{sun}). Although radio AGNs are found to be prevalent among the BCGs, the infrared luminosities of these three BCGs, judged from the infrared SED signatures, are likely to be powered by star formation. Considering the overall trend that clusters with shorter radiative gas cooling times harbor more infrared-luminous BCGs, the enhanced star formation may be caused by the cooling cluster gas accreting onto the BCGs.

We are currently undertaking a Spitzer GTO program to image ~30 massive lensing clusters at moderate redshift with both IRAC and MIPS. By taking advantage of the gravitatinoal lensing power of these clusters, we will study the population of faint galaxies that are below the nominal Spitzer detection limits. Here, we present a few examples of our science programs.

Mid-infrared observations of Active Galactic Nuclei (AGN) are important for understanding of the physical conditions around the central accretion engines. Chandra and XMM-Newton X-ray observations of a 300 arcmin^2 region in the Extended Groth Strip are used to select a sample of~150 AGN. The Spitzer instruments IRAC and MIPS detect 68-80% of these sources, which show a wide range of mid-infrared properties. About 40% of the sources have red power-law spectral energy distributions (f_nu ~ nu^alpha, alpha<0) in the 3.6-8 um IRAC bands. In these sources the central engine dominates the emission at both X-ray and IR wavelengths. Another 40% of the sources have blue mid-IR spectral energy distributions (alpha>0) with their infrared emission dominated by the host galaxy; the remaining 20% are not well-fit by a power law. Published IRAC color criteria for AGN select most of the red sources, but only some of the blue sources. As with all other known methods, selecting AGN with mid-IR colors will not produce a sample that is simultaneously complete and reliable. The IRAC SED type does not directly correspond to X-ray spectral type (hard/soft). The mid-IR properties of X-ray-detected Lyman-break, radio, submillimeter, and optically-faint sources vary widely and, for the most part, are not distinct from those of the general X-ray/infrared source population. X-ray sources emit 6-11% of the integrated mid--IR light, making them significant contributors to the cosmic infrared background.

We investigate the redshift dependence of the contribution to the extragalactic far-infrared/sub-millimeter background from galaxies detected by the Spitzer Space Telescope at 8um and 24um. Using seven-band optical to mid-infrared photometry, we estimate photometric redshifts for the Spitzer sources which appear to be mostly L* galaxies at a median redshift of z=1.0. These sources, extracted from deep 8um and 24um mosaics of the CUDSS 14-hour field with 5sigma limits of 5.8uJy and 70uJy respectively, exhibit significant 850um and 450um emission as observed by SCUBA. At 850um, after removing >=4sigma sources and those securely identified in our previous cross-matching paper, we measure stacked flux at the significance level of 4.4sigma and 2.9sigma from the full 8um and 24um galaxy catalogue respectively. At 450um, flux is detected from all 8um galaxies at the level of 3.5sigma, while there is no significant emission from the 24um galaxies. We find that the 850um flux is emitted almost exclusively at z>~1.3 from the Spitzer sources with 0.44mJy (4.7sigma) per 8um source and 0.51mJy (2.8sigma) per 24um source. This corresponds to a contribution of (16+/-3)% toward the 850um extra-galactic background from the 8um sources and (5.0+/-1.8)% from the 24um sources. At 450um, only the 8um sources within the redshift interval 1<z<2 exhibit significant emission with an average flux per source of 3.35mJy (3.0sigma). This is a contribution of (37+/-12)% to the 450um background.