Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present ``total light" stacking techniques spanning 3.4-500{\mu}m aimed at revealing the total cluster IR emission, including low mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our stacking techniques to 232 well-studied massive (log M200/Msun~13.8) clusters across multiple z bins, recovering extended cluster emission at all wavelengths, typically at >5sigma. We measure the averaged near- to far-IR radial profiles and SEDs, quantifying the total stellar and dust content. The near-IR radial profiles are well described by an NFW model with a high (c~7) concentration parameter. Dust emission is similarly concentrated, albeit suppressed at small radii (r<0.2Mpc). The measured SEDs lack warm dust, consistent with the colder SEDs expected for low mass galaxies. We derive total stellar masses consistent with the theoretical Mhalo-M_star relation and specific-star formation rates that evolve strongly with redshift, echoing that of massive (log Mstar/Msun>10) cluster galaxies. Separating out the massive galaxy population reveals that the majority of cluster far-IR emission (~70-80%) is provided by the low mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low mass cluster galaxies.

In this work, we collect quasi-simultaneous infrared, optical, X-ray and $\gamma$-ray data of 60 $Fermi$-4LAC flat spectrum radio quasars (FSRQs). In the framework of the conventional one-zone leptonic model, we investigate the physical properties of $Fermi$-4LAC FSRQs' jets by modeling their quasi-simultaneous spectral energy distributions (SEDs). Our main results are summarized as follows. (1) There is a linear correlation between synchrotron peak frequency and curvature of the electron energy distribution. As suggested by previous works, the slope of the best linear fitting equation of this correlation is consistent with statistic acceleration which needs a fluctuation of fractional acceleration gain. (2) The gamma-ray dissipation regions are located at the range from 0.1 to 10 pc away from the super-massive black hole, and located outside the broad-line region (BLR) and within the dusty torus (DT). (3) A size relation $P_{\rm e}$ (the kinetic power carried in relativistic electrons) $\sim$ $P_{\rm B}$ (Poynting flux) $\leq$ $P_{\rm r}$ (the radiative power ) $<$ $P_{\rm p}$ (the kinetic power in cold protons) is found in our modeling. Among them, $P_{\rm e}\sim P_{\rm B}$ suggests that SEDs of almost all FSRQs with parameters are close to equipartition between the magnetic field and the relativistic electrons. The $P_{\rm e} < P_{\rm r}$ suggest that the most energy of the relativistic electrons are dissipated by EC radiation for FSRQs. (4) There is an anti-correlation between the peak energy of SEDs ($\gamma_{\rm peak}$) and the jet power ($P_{\rm jet}$), which is consistent with the blazar sequence.

Very high energy photons from cosmological gamma-ray bursts (GRBs) are expected to interact with extragalactic background light (EBL) and produce electron-positron pairs when they propagate through intergalactic medium (IGM). These relativistic pairs will then up-scatter cosmic microwave background (CMB) photons and emit secondary GeV emission. Meanwhile, the motion of these pairs are deflected by intergalactic magnetic field (IGMF), so the secondary GeV photons arrive later than the primary emission. It has been suggested that the properties of the secondary GeV emission can be used to constrain IGMF. Recently, TeV gamma-ray emission has been detected, for the first time, from a GRB (GRB 190114C) by the MAGIC telescope and its steep ${\rm \gamma-ray}$ spectrum shows a clear evidence of absorption by EBL. We then constrain the IGMF with the GeV flux limit obtained from the $Fermi$-LAT observations. We find a limit of $>10^{-19.5}$ G for the coherence length of $\lambda \leq 1$ Mpc. Although this limit is weaker than that obtained by using blazars, it represents the first limit from ${\rm \gamma-ray}$ observations of GRBs, which provides an independent constraint on IGMF. We also find that, for transient ${\rm \gamma-ray}$ sources, one can choose a favorable time window to search for the echo emission at a particular energy.

We continuously monitored the blazar S5 0716+714 in the optical $g$, $r$ and $i$ bands from Nov. 10, 2017 to Jun. 06, 2019. The total number of observations is 201 nights including 26973 data points. This is a very large quasi-simultaneous multicolor sample for the blazar. The average time spans and time resolutions are 3.4 hours and 2.9 minutes per night, respectively. During the period of observations, the target source in the $r$ band brightens from $14^{\rm m}.16$ to $12^{\rm m}.29$ together with five prominent sub-flares, and then first becomes fainter to $14^{\rm m}.76$ and again brightens to $12^{\rm m}.94$ with seven prominent sub-flares. For the long-term variations, we find a strong flatter when brighter (FWB) trend at a low flux state and then a weak FWB trend at a higher flux state. A weak FWB trend at a low flux state and then a strong FWB trend at a higher flux state are also reported. Most of sub-flares show the strong FWB trends, except for two flares with a weak FWB trend. The particle acceleration and cooling mechanisms together with the superposition of different FWB-slopes from sub-flares are likely to explain the optical color behaviours. A scenario of bent jet is discussed.

The low energy (LE) X-ray telescope is one of the three main instruments of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). It is equipped with Swept Charge Device (SCD) sensor arrays with a total geometrical area of 384 cm2 and an energy band from 0.7 keV to 13 keV. In order to evaluate the particle induced X-ray background and the cosmic X-ray background simultaneously, LE adopts collimators to define four types of Field Of Views (FOVs). LE is constituted of three detector boxes (LEDs) and an electric control box (LEB) and achieves a good energy resolution of 140 eV at 5.9 keV, an excellent time resolution of 0.98 ms, as well as an extremely low pileup (<1% at 18000 cts/s). Detailed performance tests and calibration on the ground have been performed, including energy-channel relation, energy response, detection efficiency and time response.

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $p\gamma$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

The origin of the multi-band activities (outbursts/flares) of blazars is still a heavily debated topic. Shock and magnetic reconnection have long been considered as possible triggers for the multi-band activities. In this paper, we present an exploration of the origin of multi-band activities for a high-redshift (z =1.8385) FSRQ PKS 1502+106. Utilizing multi-band data from radio to $\gamma$-ray and optical polarization observations, we investigate two dramatic activities in detail: a $\gamma$-ray dominated outburst in 2015 and an optical dominated outburst in 2017. Our main results are as follows. (I) A fast $\gamma$-ray flare with a flux-doubling time-scale as short as 1-hr in 2015 is discovered. Based on the variability time-scale, the physical parameters of the flaring region (e.g, minimum Doppler factor, emission region size, etc.) are constrained. At the peak of the flare, the $\gamma$-ray spectrum hardens to $\Gamma_{\gamma} = 1.82\pm0.04$ and exhibits an obvious curvature/break characteristic that is caused by the typical "cooling break". Modelings of multi-band SEDs reveal a very hard electronic energy spectrum with the electronic spectral index of $1.07\pm0.53$. This result suggests that this fast $\gamma$-ray flare may be triggered by magnetic reconnection. (II) During the outburst in 2017, the optical polarization degree and optical fluxes show a very tight correlation. By analyzing Stokes parameters of polarization observations, our results show that this outburst could be triggered by a transverse shock with a compression ratio of $\eta> 2.2$, and the magnetic field intensity of the shock emission region is about $0.032$ G.

A central compact object (CCO, e.g. a black hole) with an accretion disk has been suggested as the common central engine of various astrophysical phenomena, such as gamma-ray bursts (GRBs), tidal disruption events (TDEs) and active galactic nuclei (AGNs). A jet powered by such a system might precess due to the misalignment of the angular momenta of the CCO and accretion disk. Some quasi-periodic behaviors observed in the light curves of these phenomena can be well interpreted within the framework of a precessing jet model. In this paper, we study the emission polarization of precessing jets in the three kinds of phenomena. The polarization angle also shows a gradual change for the synchrotron emission in both the random and toroidal magnetic field configurations with the precessing jet, while it can only change abruptly by $90^\circ$ for the non-precessing top-hat jet. Polarization properties are periodic due to the assumptions made in our model. The polarization observations are crucial to confirm the precession nature of jets in GRBs, TDEs and AGNs.

Although many high-energy neutrinos detected by the IceCube telescope are believed to have anextraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedenteddiscovery of a high-energy muon neutrino event coincident with a multiwavelength flare from ablazar, TXS 0506+056, shed some light on the origin of the neutrinos. It is usually believed that ablazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH).Here we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclearinteractions between the accelerated cosmic-ray protons in the relativistic jet and the dense gasclouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in thejet, which could be much smaller than that required in the conventional hadronic model whichinstead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare fromthe optical to gamma-ray band can be well explained by invoking a second radiation zone in thejet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time thanthe multiwavelength flare so the neutrino event is not necessarily correlated with the flare but it is probably accompanied by a spectrum hardening above a few GeV.

Carbon monoxide (CO) observations show a luminosity$-$line-width correlation that evolves with redshift. We present a method to use CO measurements alone to infer the molecular gas fraction ($f_{\rm mol}$) and constrain the CO$-$H$_2$ conversion factor ($\alpha_{\rm CO}$). We compile from the literature spatially integrated low-$J$ CO observations of six galaxy populations, including a total of 449 galaxies between $0.01 \leq z \leq 3.26$. The CO data of each population provide an estimate of the $\alpha_{\rm CO}$-normalized mean molecular gas fraction ($f_{\rm mol}/\alpha_{\rm CO}$). The redshift evolution of the luminosity$-$line-width correlation thus indicates an evolution of $f_{\rm mol}/\alpha_{\rm CO}$. We use a Bayesian-based Monte-Carlo Markov Chain sampler to derive the posterior probability distribution functions of $f_{\rm mol}/\alpha_{\rm CO}$ for these galaxy populations, accounting for random inclination angles and measurement errors in the likelihood function. We find that the molecular gas fraction evolves rapidly with redshift, $f_{\rm mol} \propto (1+z)^\beta$ with $\beta \simeq 2$, for both normal star-forming and starburst galaxies. Furthermore, the evolution trend agrees well with that inferred from the Kennicutt-Schmidt relation and the star-forming main sequence. Finally, at $z < 0.1$ normal star-forming galaxies require a $\sim5\times$ larger $\alpha_{\rm CO}$ than starburst galaxies to match their molecular gas fractions, but at $z > 1$ both star-forming types exhibit sub-Galactic $\alpha_{\rm CO}$ values and normal star-forming galaxies appear more gas-rich than starbursts. Future applications of this method include calibrating Tully-Fisher relations without inclination correction and inferring the evolution of the atomic gas fraction with HI observations.

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies.

We study the requirement on the jet power in the conventional p-$\gamma$ models (photopion production and Bethe-Heitler pair production) for TeV BL Lac objects. We select a sample of TeV BL Lac objects whose SEDs are difficult to be explained by the one-zone leptonic model. Based on the relation between the p-$\gamma$ interaction efficiency and the opacity of $\gamma\gamma$ absorption, we find that the detection of TeV emission poses upper limits on the p-$\gamma$ interaction efficiencies in these sources and hence minimum jet powers can be derived accordingly. We find that the obtained minimum jet powers exceed the Eddington luminosity of the supermassive black holes. Implications for the accretion mode of the supermassive black hole in these BL Lac objects and the origin of their TeV emissions are discussed.

We investigate the spatial distribution of Lyman alpha emitting galaxies (LAEs) at $z \approx 2.67$, selected from the NOAO Deep Wide-Field Survey (NDWFS), using two-point statistics and topological diagnostics adopted from network science. We measure the clustering length, $r_0 \approx 4 h^{-1}$ Mpc, and the bias, $b_{LAE} = 2.2^{+0.2}_{-0.1}$. Fitting the clustering with halo occupation distribution (HOD) models results in two disparate possibilities: (1) where the fraction of central galaxies is $<$1% in halos of mass $>10^{12}$$M_\odot$; and (2) where the fraction is $\approx$20%. We refer to these two scenarios as the `Dusty Core Scenario' for Model#1 since most of central galaxies in massive halos are dead in Ly$\alpha$ emission, and the `Pristine Core Scenario' for Model#2 since the central galaxies are bright in Ly$\alpha$ emission. Traditional two-point statistics cannot distinguish between these disparate models given the current data sets. To overcome this degeneracy, we generate mock catalogs for each HOD model using a high resolution $N$-body simulation and adopt a network statistics approach, which provides excellent topological diagnostics for galaxy point distributions. We find three topological anomalies from the spatial distribution of observed LAEs, which are not reproduced by the HOD mocks. We find that Model#2 matches better all network statistics than Model#1, suggesting that the central galaxies in $> 10^{12} h^{-1} M_\odot$ halos at $z \approx 2.67$ need to be less dusty to be bright as LAEs, potentially implying some replenishing channels of pristine gas such as the cold mode accretion.

The brightest 500$\,\mu$m source in the XMM field, HXMM01, is a rare merger of luminous starburst galaxies at $z=2.3$ with a dust-obscured star-formation rate of 2,000$\,M_{\odot}\,{\rm yr}^{-1}$. Here we present high-resolution spectroscopic observations of HXMM01 with the Atacama Large Millimeter/submillimeter Array (ALMA). We detect line emission from ${\rm CO\,{\it J}=7\to6}$, [C I]$\,{{^3P_2}\to{^3P_1}}$, and p-${\rm H_{2}O}\,{2_{11}}\to{2_{02}}$ and continuum emission at $230\,$GHz. At a spatial resolution of 0.2" and a spectral resolution of 40$\,\rm km\,s^{-1}$, the source is resolved into three distinct components, which are spatially and dynamically associated within a projected radius of 20$\,$kpc and a radial velocity range of 2,000$\,\rm km\,s^{-1}$. For two major components, our Bayesian-based tilted-ring modeling of the ALMA spectral cubes shows almost flat rotation curves peaking at $\sim500\,\rm km\,s^{-1}$ at galactocentric distances between 2 and 5$\,$kpc. Each of them has a dynamical mass of $\sim10^{11}\,M_\odot$. The combination of the dynamical masses and the archival ${\rm CO\,{\it J}=1\to0}$ data places strong upper limits on the CO$\to$H$_2$ conversion factor of $\alpha_{\rm CO}\lesssim1.4-2.0\,{M_{\odot}}\,\rm (K\,km\,s^{-1}\,pc^{2})^{-1}$. These limits are significantly below the Galactic inner disk $\alpha_{\rm CO}$ value of $4.3\,{M_{\odot}}\,\rm (K\,km\,s^{-1}\,pc^{2})^{-1}$ but are consistent with those of local starbursts. Therefore, the previously estimated short gas depletion timescale of $\sim200\,$Myr remains unchanged.

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe.

We have identified 105 galaxy pairs at z ~ 0.04 with the MaNGA integral-field spectroscopic data. The pairs have projected separations between 1 kpc and 30 kpc, and are selected to have radial velocity offsets less than 600 km/s and stellar mass ratio between 0.1 and 1. The pair fraction increases with both the physical size of the integral-field unit and the stellar mass, consistent with theoretical expectations. We provide the best-fit analytical function of the pair fraction and find that ~3% of M* galaxies are in close pairs. For both isolated galaxies and paired galaxies, active galactic nuclei (AGN) are selected using emission-line ratios and H_alpha equivalent widths measured inside apertures at a fixed physical size. We find AGNs in ~24% of the paired galaxies and binary AGNs in ~13% of the pairs. To account for the selection biases in both the pair sample and the MaNGA sample, we compare the AGN comoving volume densities with those expected from the mass- and redshift-dependent AGN fractions. We find a strong (~5x) excess of binary AGNs over random pairing and a mild (~20%) deficit of single AGNs. The binary AGN excess increases from ~2x to ~6x as the projected separation decreases from 10-30 kpc to 1-10 kpc. Our results indicate that pairing of galaxies preserves the AGN duty cycle in individual galaxies but increases the population of binary AGNs through correlated activities. We suggest tidally-induced galactic-scale shocks and AGN cross-ionization as two plausible channels to produce low-luminosity narrow-line-selected binary AGNs.

We present a $^{13}\mathrm{CO} (J = 1 \rightarrow 0)$ mapping survey of 12 nearby galaxies from the CARMA STING sample. The line intensity ratio $\mathcal{R} \equiv I[^{12}\mathrm{CO} (J = 1 \rightarrow 0)]/I[^{13}\mathrm{CO} (J = 1 \rightarrow 0)]$ is derived to study the variations in molecular gas properties. For 11 galaxies where it can be measured with high significance, the spatially resolved $\mathcal{R}$ on (sub-)kiloparsec scales varies by up to a factor of 3--5 within a galaxy. Lower $\mathcal{R}$ values are usually found in regions with weaker $^{12}\rm CO$. We attribute this apparent trend to a bias against measuring large $\mathcal{R}$ values when $^{12}\rm CO$ is weak. Limiting our analysis to the $^{12}\rm CO$ bright regions that are less biased, we do not find $\mathcal{R}$ on (sub)kpc scales correlate with galactocentric distance, velocity dispersion or the star formation rate. The lack of correlation between SFR and $\mathcal{R}$ indicates that the CO optical depth is not sensitive to stellar energy input, or that any such sensitivity is easily masked by other factors. Extending the analysis to all regions with $\rm ^{12}CO$ emission by spectral stacking, we find that 5 out of 11 galaxies show higher stacked $\mathcal{R}$ for galactocentric radii of $\gtrsim 1$ kpc and $\Sigma_{\mathrm{SFR}} \lesssim 0.1 \rm \ M_{sun} yr^{-1} kpc^{-2}$, which could result from a greater contribution from diffuse gas. Moreover, significant galaxy-to-galaxy variations are found in $\mathcal{R}$, but the global $\mathcal{R}$ does not strongly depend on dust temperature, inclination, or metallicity of the galaxy.

We collect a sample of 78 Fermi detected blazars with broad line region (BLR) data and the quasi-simultaneous multi-wave band data. By analyzing the sample, we find: (1) For whole blazar sample, there exist significant correlations between radio, $\gamma$-ray luminosity and BLR luminosity. The slope of the best-fit linear regression equation is close to the theoretical value 1. Our results provided another support of jet-disk symbiosis. (2) For FSRQ sample, we find a significant correlation between radio luminosity and optical luminosity, and weaker correlations between X-ray luminosity and radio, optical luminosity. It indicates that the X-ray band radiation mechanism is slightly different from radio and optical bands. The X-ray emission is diluted by inverse Compton scattering (IC). There is a significant correlation between BLR luminosity and $\gamma$-ray luminosity, but no correlations between X-ray luminosity and BLR, $\gamma$-ray luminosity, which suggests that the seed soft photons are different for IC process. From these results we suggest that the X-ray band emission is diluted by synchrotron self-Compton (SSC) process, $\gamma$-ray band emission is produced by external Compton (EC) process.

Bolometric luminosity is a basic physical parameter that is widely used in the study of blazars. Due to the lack of simultaneous full wavelength data, several estimators of the bolometric luminosity are being used in practice. In this paper, we study and evaluate the reliability and significance of six estimators, the $5GHz$ luminosity, the 1keV luminosity, the $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity+the $\gamma$-ray luminosity and the 1keV luminosity+the $\gamma$-ray luminosity, by analyzing the linear correlations between the integrated bolometric luminosity and them. Our main results are as follows. (i) All the six estimators are reliable in the sense that they are all significant correlated with the bolometric luminosity. (ii) Ranking from the higher significance of the reliability to lower one the six estimators are the the $5GHz$ luminosity+the $\gamma$-ray luminosity, the 1keV luminosity+the $\gamma$-ray luminosity, $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity and the 1keV luminosity. (iii) We suggest that the bolometric luminosity can be well estimated by the $\gamma$-ray luminosity using the best linear equation that given in this paper for Fermi FSRQs. (iv) According to the linear regressions obtained in the analysis, we provide calibration for each estimator.

We report the detection of diffuse Ly$\alpha$ emission, or Ly$\alpha$ halos (LAHs), around star-forming galaxies at $z\approx3.78$ and $2.66$ in the NOAO Deep Wide-Field Survey Bo\"otes field. Our samples consist of a total of $\sim$1400 galaxies, within two separate regions containing spectroscopically confirmed galaxy overdensities. They provide a unique opportunity to investigate how the LAH characteristics vary with host galaxy large-scale environment and physical properties. We stack Ly$\alpha$ images of different samples defined by these properties and measure their median LAH sizes by decomposing the stacked Ly$\alpha$ radial profile into a compact galaxy-like and an extended halo-like component. We find that the exponential scale-length of LAHs depends on UV continuum and Ly$\alpha$ luminosities, but not on Ly$\alpha$ equivalent widths or galaxy overdensity parameters. The full samples, which are dominated by low UV-continuum luminosity Ly$\alpha$ emitters ($M_{\rm UV} \gtrsim -21$), exhibit LAH sizes of 5$\,-\,6\,$kpc. However, the most UV- or Ly$\alpha$-luminous galaxies have more extended halos with scale-lengths of 7$\,-\,9\,$kpc. The stacked Ly$\alpha$ radial profiles decline more steeply than recent theoretical predictions that include the contributions from gravitational cooling of infalling gas and from low-level star formation in satellites. On the other hand, the LAH extent matches what one would expect for photons produced in the galaxy and then resonantly scattered by gas in an outflowing envelope. The observed trends of LAH sizes with host galaxy properties suggest that the physical conditions of the circumgalactic medium (covering fraction, HI column density, and outflow velocity) change with halo mass and/or star-formation rates.

We fit the spectral energy distributions (SEDs) of members of a large sample of Fermi 2LAC blazars to synchrotron and inverse Compton (IC) models. Our main results are as follows. (i) As suggested by previous works, the correlation between peak frequency and curvature can be explained by statistical or stochastic particle acceleration mechanisms. For BL Lacs, we find a linear correlation between synchrotron peak frequency and its curvature. The slope of the correlation is consistent with the stochastic acceleration mechanisms and confirm previous studies. For FSRQs, we also find a linear correlation, but its slope cannot be explained by previous theoretical models. (ii) We find a significant correlation between IC luminosity and synchrotron luminosity. The slope of the correlation of FSRQs is consistent with the EC process. And the slope of the correlation of BL Lac is consistent with the SSC process. (iii) We find several significant correlations between IC curvature and several basic parameters of blazars (black hole mass, broad line luminosity, the Lorentz factor of jet). We also find significant correlations between bolometric luminosity and these basic parameters of blazars which suggest that the origin of jet is a mixture of the mechanisms proposed by Blandford $\&$ Znajek and by Blandford $\&$ Payne.

We report high spatial resolution observations of Giant Molecular Clouds (GMCs) in the nearby spiral galaxies NGC 6946, M101 and NGC 628 obtained with the CARMA telescope. We observed CO(1-0) over regions with active star formation, and higher resolution observations of CO(2-1) have allowed us to resolve some of the largest GMCs. Using a Bayesian fitting approach, we generate scaling relations for the sizes, line widths, and virial masses of the structures identified in this work. We do not find evidence for a tight power law relation between size and line width, although the limited dynamic range in cloud size remains a clear issue in our analysis. Additionally, we use a Bayesian approach to analyze the Kennicutt-Schmidt (K-S) relation for the identified structures. We find that the distribution of slopes are broadly distributed, mainly due to the limited dynamic range of our measured H2 mass surface density, and being most consistent with super-linear relations. On the other hand, when we use the Bayesian approach to analyze the K-S relation for a uniform grid, the distributions of slopes is consistent with sub-linear relations. On-arm regions tend to have higher star formation rates than inter-arm regions. As in NGC 6946, in M101 we find regions where the star formation efficiency shows marked peaks at specific galoctocentric radii. On the other hand, the distribution of SFE in NGC 628 is more contiguous. We hypothesize that differences in the distribution of SFE may be indicative of different processes driving the spiral structure.

We present a study of the vertical structure of the gaseous and stellar disks in a sample of edge-on galaxies (NGC 4157, 4565, and 5907) using BIMA/CARMA 12CO (J = 1 --> 0), VLA H I, and Spitzer 3.6 micron data. In order to take into account projection effects when we measure the disk thickness as a function of radius, we first obtain the inclination by modeling the radio data. Using the measurement of the disk thicknesses and the derived radial profiles of gas and stars, we estimate the corresponding volume densities and vertical velocity dispersions. Both stellar and gas disks have smoothly varying scale heights and velocity dispersions, contrary to assumptions of previous studies. Using the velocity dispersions, we find that the gravitational instability parameter Q follows a fairly uniform profile with radius and is greater than or equal to 1 across the star forming disk. The star formation law has a slope that is significantly different from those found in more face-on galaxy studies, both in deprojected and pixel-by-pixel plots. Midplane gas pressure based on the varying scale heights and velocity dispersions appears to roughly hold a power-law correlation with the midplane volume density ratio.

Massive stars can alter physical conditions and properties of their ambient interstellar dust grains via radiative heating and shocks. The HII regions in the Large Magellanic Cloud (LMC) offer ideal sites to study the stellar energy feedback effects on dust because stars can be resolved, and the galaxy's nearly face-on orientation allows us to unambiguously associate HII regions with their ionizing massive stars. The Spitzer Space Telescope survey of the LMC provides multi-wavelength (3.6 to 160 $\mu$m) photometric data of all HII regions. To investigate the evolution of dust properties around massive stars, we have analyzed spatially-resolved IR dust emission from two classical HII regions (N63 and N180) and two simple superbubbles (N70 and N144) in the LMC. We produce photometric spectral energy distributions (SEDs) of numerous small subregions for each region based on its stellar distributions and nebular morphologies. We use DustEM dust emission model fits to characterize the dust properties. Color-color diagrams and model fits are compared with the radiation field (estimated from photometric and spectroscopic surveys). Strong radial variations of SEDs can be seen throughout the regions, reflecting the available radiative heating. Emission from very small grains drastically increases at locations where the radiation field is the highest, while polycyclic aromatic hydrocarbons (PAHs) appear to be destroyed. PAH emission is the strongest in the presence of molecular clouds, provided that the radiation field is low.

We investigate the correlation between CO and HI emission in 18 nearby galaxies from the CARMA Survey Toward IR-Bright Nearby Galaxies (STING) at sub-kpc and kpc scales. Our sample, spanning a wide range in stellar mass and metallicity, reveals evidence for a metallicity dependence of the HI column density measured in regions exhibiting CO emission. Such a dependence is predicted by the equilibrium model of McKee & Krumholz, which balances H_2 formation and dissociation. The observed HI column density is often smaller than predicted by the model, an effect we attribute to unresolved clumping, although values close to the model prediction are also seen. We do not observe HI column densities much larger than predicted, as might be expected were there a diffuse HI component that did not contribute to H_2 shielding. We also find that the H_2 column density inferred from CO correlates strongly with the stellar surface density, suggesting that the local supply of molecular gas is tightly regulated by the stellar disk.

We have determined column densities of H I and/or H_2 for sight lines in the Magellanic Clouds from archival HST and FUSE spectra of H I Lyman-alpha and H_2 Lyman-band absorption. Together with some similar data from the literature, we now have absorption-based N(H I) and/or N(H_2) for 285 LMC and SMC sight lines (114 with a detection or limit for both species) -- enabling more extensive, direct, and accurate determinations of molecular fractions, gas-to-dust ratios, and elemental depletions in these two nearby, low-metallicity galaxies. For sight lines where the N(H I) estimated from 21 cm emission is significantly higher than the value derived from Lyman-alpha absorption (presumably due to emission from gas beyond the target stars), integration of the 21 cm profile only over the velocity range seen in Na I or H_2 absorption generally yields much better agreement. Conversely, N(21 cm) can be lower than N(Ly-alpha) by factors of 2--3 in some LMC sight lines -- suggestive of small-scale structure within the 21 cm beam(s) and/or some saturation in the emission. The mean gas-to-dust ratios obtained from N(H_tot)/E(B-V) are larger than in our Galaxy, by factors of 2.8--2.9 in the LMC and 4.1--5.2 in the SMC -- i.e., factors similar to the differences in metallicity. The N(H_2)/E(B-V) ratios are more similar in the three galaxies, but with considerable scatter within each galaxy. These data may be used to test models of the atomic-to-molecular transition at low metallicities and predictions of N(H_2) based on comparisons of 21 cm emission and the IR emission from dust.

We present an analysis of the relationship between molecular gas and current star formation rate surface density at sub-kpc and kpc scales in a sample of 14 nearby star-forming galaxies. Measuring the relationship in the bright, high molecular gas surface density ($\Shtwo\gtrsim$20 \msunpc) regions of the disks to minimize the contribution from diffuse extended emission, we find an approximately linear relation between molecular gas and star formation rate surface density, $\nmol\sim0.96\pm0.16$, with a molecular gas depletion time $\tdep\sim2.30\pm1.32$ Gyr. We show that, in the molecular regions of our galaxies there are no clear correlations between \tdep\ and the free-fall and effective Jeans dynamical times throughout the sample. We do not find strong trends in the power-law index of the spatially resolved molecular gas star formation law or the molecular gas depletion time across the range of galactic stellar masses sampled (\mstar $\sim$$10^{9.7}-10^{11.5}$ \msun). There is a trend, however, in global measurements that is particularly marked for low mass galaxies. We suggest this trend is probably due to the low surface brightness CO, and it is likely associated with changes in CO-to-H2 conversion factor.

This paper examines the detailed sub-milliarcsecond polarization properties of an individual SiO maser feature displaying a rotation in polarization electric vector position angle of approximately \pi/2 across the feature. Such rotations are a characteristic observational signature of circumstellar SiO masers detected toward a number of late-type, evolved stars. We employ a new calibration method for accurate circular VLBI polarimetry at millimeter wavelengths, to present the detailed Stokes {I,Q,U,V} properties for this feature. We analyze the fractional linear and circular polarization as a function of projected angular distance across the extent of the feature, and compare these measurements against several theoretical models proposed for sharp rotations of electric vector position angle in polarized SiO maser emission. We find that the rotation is most likely caused by the angle \theta between the line of sight and a projected magnetic field crossing the critical Van Vleck angle for maser propagation. The fractional linear polarization profile m_l(\theta) is well-fit by standard models for polarized maser transport, but we find less agreement for the fractional circular polarization profile m_c(\theta).

This study explores the effects of different assumptions and systematics on the determination of the local, spatially resolved star formation law. Using four star formation rate (SFR) tracers (H\alpha with azimuthally averaged extinction correction, mid-infrared 24 micron, combined H\alpha and mid-infrared 24 micron, and combined far-ultraviolet and mid-infrared 24 micron), several fitting procedures, and different sampling strategies we probe the relation between SFR and molecular gas at various spatial resolutions and surface densities within the central 6.5 kpc in the disk of NGC4254. We find that in the high surface brightness regions of NGC4254 the form of the molecular gas star formation law is robustly determined and approximately linear and independent of the assumed fraction of diffuse emission and the SFR tracer employed. When the low surface brightness regions are included, the slope of the star formation law depends primarily on the assumed fraction of diffuse emission. In such case, results range from linear when the fraction of diffuse emission in the SFR tracer is ~30% or less (or when diffuse emission is removed in both the star formation and the molecular gas tracer), to super-linear when the diffuse fraction is ~50% and above. We find that the tightness of the correlation between gas and star formation varies with the choice of star formation tracer. The 24 micron SFR tracer by itself shows the tightest correlation with the molecular gas surface density, whereas the H\alpha corrected for extinction using an azimuthally-averaged correction shows the highest dispersion. We find that for R<0.5R_25 the local star formation efficiency is constant and similar to that observed in other large spirals, with a molecular gas depletion time ~2 Gyr.

The synchrotron self-Compton (SSC) emission from Gamma-ray Burst (GRB) forward shock can extend to the very-high-energy (VHE; $E_\gamma > $100 GeV) range. Such high energy photons are rare and are attenuated by the cosmic infrared background before reaching us. In this work, we discuss the prospect to detect these VHE photons using the current ground-based Cherenkov detectors. Our calculated results are consistent with the upper limits obtained with several Cherenkov detectors for GRB 030329, GRB 050509B, and GRB 060505 during the afterglow phase. For 5 bursts in our nearby GRB sample (except for GRB 030329), current ground-based Cherenkov detectors would not be expected to detect the modeled VHE signal. Only for those very bright and nearby bursts like GRB 030329, detection of VHE photons is possible under favorable observing conditions and a delayed observation time of $\la$10 hours.

We intend to determine the type of circumburst medium and measure directly the initial Lorentz factor $\Gamma_0$ of GRB outflows. If the early X-ray afterglow lightcurve has a peak and the whole profile across the peak is consistent with the standard external shock model, the early rise profile of light curves can be used to differentiate whether the burst was born in interstellar medium (ISM) or in stellar wind. In the thin shell case, related to a sub-relativistic reverse shock, the peak time occurring after the end of the prompt emission, can be used to derive an accurate $\Gamma_0$, especially for the ISM case. The afterglow lightcurves for a flat electron spectrum $1<p<2$ have been derived analytically. In our GRB sample, we obtain $\Gamma_0 \sim 300$ for the bursts born in ISM. We did not find any good case for bursts born in stellar wind and behaving as a thin shell that can be used to constrain $\Gamma_0$ reliably.

Gamma-ray Bursts (GRBs) are among the potential extragalactic sources of very-high-energy (VHE) gamma-rays. We discuss the prospects of detecting VHE gamma-rays with current ground-based Cherenkov instruments during the afterglow phase. Using the fireball model, we calculate the synchrotron self-Compton (SSC) emission from forward-shock electrons. The modeled results are compared with the observational afterglow data taken with and/or the sensitivity level of ground-based VHE instruments (e.g. STACEE, H.E.S.S., MAGIC, VERITAS, and Whipple). We find that modeled SSC emission from bright and nearby bursts such as GRB 030329 are detectable by these instruments even with a delayed observation time of ~10 hours.

In the neutron-rich internal shocks model for gamma-ray bursts (GRBs), the Lorentz factors (LFs) of ion shells are variable, and so are the LFs of accompanying neutron shells. For slow neutron shells with a typical LF of approximate tens, the typical $\beta$-decay radius is $\sim 10^{14}-10^{15}$ cm. As GRBs last long enough [$T_{90}>14(1+z)$ s], one earlier but slower ejected neutron shell will be swept successively by later ejected ion shells in the range $\sim10^{13}-10^{15}$ cm, where slow neutrons have decayed significantly. Part of the thermal energy released in the interaction will be given to the electrons. These accelerated electrons will be mainly cooled by the prompt soft $\gamma-$rays and give rise to GeV emission. This kind of GeV emission is particularly important for some very long GRBs and is detectable for the upcoming satellite {\it Gamma-Ray Large Area Space Telescope} (GLAST).

We present the results from arcsecond resolution observations of various line transitions at 1.3 mm toward hypercompact HII region G28.20-0.04N. With the SMA data, we have detected and mapped the transitions in the CH$_{3}$CN, CO, $^{13}$CO, SO$_{2}$, OCS, and CH$_{3}$OH molecular lines as well as the radio recombination line H30$\alpha$. The observations and analysis indicate a hot core associated with G28.20-0.04N. The outflow and possible rotation are detected in this region.

This article presents a multi-wavelength study towards S87, based on a dataset of submillimeter/far-/mid-infrared (sub-mm/FIR/MIR) images and molecular line maps. The sub-mm continuum emission measured with JCMT/SCUBA reveals three individual clumps, namely, SMM1, SMM2, and SMM3. The MIR/FIR images obtained by the Spitzer Space Telescope indicate that both SMM1 and SMM3 harbor point sources. The J=1-0 transitions of CO, 13CO, C18O, and HCO+, measured with the 13.7m telescope of the Purple Mountain Observatory, exhibit asymmetric line profiles. Our analysis of spectral energy distributions (SEDs) shows that all of the three sub-mm clumps are massive (110--210 $M_{\odot}$), with average dust temperatures in the range ~20--40K. A multi-wavelength comparison convinces us that the asymmetric profiles of molecular lines should result from two clouds at slightly different velocities, and it further confirms that the star-forming activity in SMM1 is stimulated by a cloud-cloud collision. The stellar contents and SEDs suggest that SMM1 and SMM3 are high-mass and intermediate-mass star-forming sites respectively. However, SMM2 has no counterpart downwards 70 micron, which is likely to be a cold high-mass starless core. These results, as mentioned above, expose multiple phases of star formation in S87.

Using the IRAM 30 m telescope, a mapping survey in optically thick and thin lines was performed towards 46 high mass star-forming regions. The sample includes UC H{\sc ii} precursors and UC H{\sc ii} regions. Seventeen sources are found to show "blue profiles", the expected signature of collapsing cores. The excess of sources with blue over red profiles ([$N_{\rm blue}$ -- $N_{\rm red}$]/$N_{\rm total}$) is 29% in the HCO$^+$ $J$=1--0 line, with a probability of 0.6% that this is caused by random fluctuations. UC H{\sc ii} regions show a higher excess (58%) than UC H{\sc ii} precursors (17%), indicating that material is still accreted after the onset of the UC H{\sc ii} phase. Similar differences in the excess of blue profiles as a function of evolutionary state are not observed in low mass star-forming regions. Thus, if confirmed for high mass star-forming sites, this would point at a fundamental difference between low- and high-mass star formation. Possible explanations are inadequate thermalization, stronger influence of outflows in massive early cores, larger gas reserves around massive stellar objects or different trigger mechanisms between low- and high- mass star formation.