We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas with a single temperature component cannot provide a consistent model to match the observational data. A simple two-phase model assuming a hot component and an ionized warm component can accommodate all the three observations. The total amount of the gas in individual halos is inferred to be comparable to the universal baryon fraction corresponding to the halo mass. The inferred temperature of the hot component is comparable to the halo virial temperature. The fraction of the hot component increases from $(15-40)\%$ for $10^{12.5}{M}_\odot$ halos to $(40-60)\%$ for $10^{14.5}{M}_\odot$ halos, where the ranges reflect uncertainties in the assumed gas density profile. Our results suggest that a significant fraction of the halo gas is in a non-thermalized component with temperature much lower than the virial temperature.

Stellar shells around galaxies could provide precious insights into their assembly history. However, their formation mechanism remains poorly empirically constrained. We present MUSE@VLT data of the most prominent outer shell of NGC474, to constrain its formation history. The stellar shell spectrum is clearly detected, with a signal-to-noise ratio of around 65 pix$^{-1}$. We use a full spectral fitting method to determine the line-of-sight velocity and the age and metallicity of the shell and associated point-like sources within the MUSE field of view. We detect six globular cluster (GC) candidates and eight planetary nebula (PN) candidates which are all kinematically associated to the stellar shell. We show that the shell has an intermediate metallicity, [M/H] = $-0.83^{+0.12}_{-0.12}$ and a possible alpha-enrichement, [alpha/Fe] ~ 0.3. This metallicity and the number of PNe are consistent with a progenitor galaxy more massive than 10^9 Msun. The photometry of the shell suggest that it is possibly composed of stars from the low metallicity outskirts of an intermediate mass (~1:10) companion. We show that at least two globular cluster candidates are quite young, with ages below 1.5 Gyr, which gives a lower-limit for the merger age. We note that spectroscopic data on the center of NGC~474 shows the presence of a young (~1~Gyr) stellar population. We conclude that the shell formation event is likely to be an intermediate-mass merger which happened around 1.5~Gyr ago and which has collaterally triggered the formation of massive, up to ~ 6 x 10^4 Msun, star cluster, and a nuclear starburst in NGC474.

We use the IllustrisTNG cosmological hydrodynamical simulation to study the formation of ultra-diffuse galaxies (UDGs) in galaxy clusters. We supplement the simulations with a realistic mass-size relation for galaxies at the time of infall into the cluster, as well as an analytical model to describe the tidally-induced evolution of their stellar mass, velocity dispersion and size. The model assumes "cuspy" NFW halos and, contrary to recent claims, has no difficulty reproducing the observed number of UDGs in clusters. Our results further suggest that the UDG population consists of a mixture of "normal" low surface brightness galaxies such as those found in the field ("born" UDGs, or B-UDGs), as well as a distinct population that owe their large size and low surface brightness to the effects of cluster tides ("tidal", or T-UDGs). The simulations indicate that T-UDGs entered the cluster earlier and should be more prevalent than B-UDGs near the cluster centres. T-UDGs should also have, at given stellar mass, lower velocity dispersion, higher metallicities, and lower dark matter content than B-UDGs. Our results suggest that systems like DF-44 are consistent with having been born as UDGs, while others such as DF2, DF4 and VLSB-D are possibly extreme T-UDG examples.

The MATLAS deep imaging survey has uncovered a plethora of dwarf galaxies in the low density environment it has mapped. A fraction of them are unusually extended and have a low-surface brightness. Among these so-called ultra-diffuse galaxies, a few seem to host an excess of globular clusters. With the integral-field unit spectrograph MUSE we have observed one of these galaxies - MATLAS J15052031+0148447 (MATLAS-2019) - located towards the nearby group NGC 5846 and measured its systemic velocity,age, and metallicity, and that of its globular clusters candidates. For the stellar body of MATLAS-2019 we derive a metallicity of -1.33+0.19-0.01 dex and an age of 11.2+1.8-0.8 Gyr. For some of the individual GCs and the stacked GC population, we derive consistent ages and metallicities. From the 11 confirmed globular clusters and using a Markov Chain Monte Carlo approach we derived a dynamical mass-to-light ratio of 4.2+8.6-3.4M/L. This is at the lower end of the luminosity-mass scaling relation defined by the Local Group dwarf galaxies. Furthermore, we couldn't confirm nor reject the possibility of a rotational component of the GC system. If present, this would further modify the inferred mass. Follow-up observations of the globular cluster population and of the stellar body of the galaxy are needed to assess whether this galaxy is lacking dark matter like it was suggested for the pair of dwarf galaxies in the field of NGC 1052, or if this is a miss-interpretation arising from systematic uncertainties of the method commonly used for these systems and the large uncertainties of the individual globular cluster velocities.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9$^{-1.8}_{+1.4}$~\kms, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major-axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of $\sigma_{\mathrm{int}}=10.6^{+3.9}_{-2.3}$~\kms. The velocity dispersion $\sigma_{\rm{DF2}\star}$(\re) for the stellar body within one effective radius is $10.8^{-4.0}_{+3.2}$~\kms. Considering various sources of systemic uncertainties this central value varies between 5 and 13~\kms, and we conservatively report a 95\% confidence upper limit to the dispersion within one \re\ of 21~\kms. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC~1052-DF2 that have been reported in the recent literature.

We survey globular clusters (GCs) in M85 using $ugi$-band images of a $1^{\circ} \times 1^{\circ}$ field obtained with the MegaCam at the 3.6 m Canada-France-Hawaii Telescope. We identify 1318 GC candidates with 20.0 mag $< g_0 <$ 23.5 mag in the entire survey region. Their radial number density profile is well fit by a S{\'e}rsic profile with $n$ = 2.58$^{+0.43}_{-0.33}$ and effective radius $R_{\rm e,GCS}$ = 4$\rlap{.}{'}$14 (= 22 kpc), showing that the candidates at $R < 20'$ are mostly genuine GCs in M85. We estimate the total number of GCs, $N$(total) = $1216^{+82}_{-50}$, and the specific frequency, $S_N = 1.41^{+0.10}_{-0.06}$. The overall color distribution of the GCs in M85 is bimodal, but the GCs in the central region at $R < 2'$ do not show a bimodal distribution clearly. The radial number density profile and surface number density map of the blue GCs (BGCs) show more extended structures than those of the red GCs (RGCs). The spatial distributions of both BGCs and RGCs are elongated, similar to that of the galaxy stellar light. The number fraction of the RGCs in the central region is much smaller compared to those in other early-type galaxies of similar luminosity. The mean $(g-i)_0$ color of the RGCs in M85 is about 0.1 mag bluer than typical values for other Virgo early-type galaxies of similar luminosity, indicating that a significant fraction of the RGCs in M85 may be younger than typical GCs. These results indicate that M85 might have undergone a major wet merger recently.

NGC 1052-DF2, an ultra diffuse galaxy (UDG), has been the subject of intense debate. Its alleged absence of dark matter, and the brightness and number excess of its globular clusters (GCs) at an initially assumed distance of 20Mpc, suggested a new formation channel for UDGs. We present the first systematic spectroscopic analysis of both the stellar body and the GCs (six of which were previously known, and one newly confirmed member) of this galaxy using MUSE@VLT. Even though NGC 1052-DF2 does not show any spatially extended emission lines we report the discovery of three planetary nebulae (PNe). We conduct full spectral fitting on the UDG and the stacked spectra of all GCs. The UDG's stellar population is old, 8.9$\pm$1.5 Gyr, metal-poor, with [M/H] = $-$1.07$\pm$0.12 with little or no $\alpha$-enrichment. The stacked spectrum of all GCs indicates a similar age of 8.9$\pm$1.8 Gyr, but lower metallicity, with [M/H] = $-$1.63$\pm$0.09, and similarly low $\alpha$-enrichment. There is no evidence for a variation of age and metallicity in the GC population with the available spectra. The significantly more metal-rich stellar body with respect to its associated GCs, the age of the population, its metallicity and alpha enrichment, are all in line with other dwarf galaxies. NGC 1052-DF2 thus falls on the same empirical mass-metallicity relation as other dwarfs, for the full distance range assumed in the literature. We find that both debated distance estimates (13 and 20 Mpc) are similarly likely, given the three discovered PNe.

Using deep, high resolution optical imaging from the Next Generation Virgo Cluster Survey we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses $10^{5}\lesssim M_{*}/M_{\odot} \lesssim10^{12}$. The nucleation fraction reaches a peak value $f_{n}\approx90\%$ for $M_{*} \approx 10^{9} M_{\odot}$ galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as $M_{*} \approx5\times10^{5} M_{\odot}$. Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at $M_{NSC}\approx7\times10^{5} M_{\odot}$, a factor 3-4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to $M_{NSC}/M_{*}\approx3.6\times10^{-3}$ for $M_{*} \approx 5\times10^{9} M_{\odot}$ galaxies. Nuclei in both more and less massive galaxies are much more prominent: $M_{NSC}\propto M_{*}^{0.46}$ at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including a very similar occupation distribution and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$.

As we demonstrated in Paper I, the quenched fractions of central and satellite galaxies as function of halo mass are extremely similar, as long as one controls for stellar mass. The same holds for the quenched fractions as a function of central velocity dispersion, which is tightly correlated with black hole mass, as long as one controls for both stellar and halo mass. Here we use mock galaxy catalogs constructed from the latest semi-analytic model, L-GALAXIES, and the state-of-the-art hydrodynamical simulation, EAGLE, to investigate whether these models can reproduce the trends seen in the data. We also check how the group finder used to identify centrals and satellites impacts our results. We find that L-GALAXIES fails to reproduce the trends. The predicted quenched fraction of central galaxies increases sharply with halo mass around $10^{12.5}h^{-1}M_{\odot}$ and with black hole mass around $\sim10^{6.5}M_{\odot}$, while the predicted quenched fraction of satellites increases with both halo and black hole masses gradually. In contrast, centrals and satellites in EAGLE follow almost the same trend as seen in the data. We discuss the implications of our results for how feedback processes regulate galaxy quenching.

Ultra-diffuse galaxies (UDGs) are unusual galaxies with low luminosities, similar to classical dwarf galaxies, but sizes up to $\sim\!5$ larger than expected for their mass. Some UDGs have large populations of globular clusters (GCs), something unexpected in galaxies with such low stellar density and mass. We have carried out a comprehensive study of GCs in both UDGs and classical dwarf galaxies at comparable stellar masses using HST observations of the Coma cluster. We present new imaging for 33 Dragonfly UDGs with the largest effective radii ($>2$ kpc), and additionally include 15 UDGs and 54 classical dwarf galaxies from the HST/ACS Coma Treasury Survey and the literature. Out of a total of 48 UDGs, 27 have statistically significant GC systems, and 11 have candidate nuclear star clusters. The GC specific frequency ($S_N$) varies dramatically, with the mean $S_N$ being higher for UDGs than for classical dwarfs. At constant stellar mass, galaxies with larger sizes (or lower surface brightnesses) have higher $S_N$, with the trend being stronger at higher stellar mass. At lower stellar masses, UDGs tend to have higher $S_N$ when closer to the center of the cluster, i.e., in denser environments. The fraction of UDGs with a nuclear star cluster also depends on environment, varying from $\sim\!40$\% in the cluster core, where it is slightly lower than the nucleation fraction of classical dwarfs, to $\lesssim20\%$ in the outskirts. Collectively, we observe an unmistakable diversity in the abundance of GCs, and this may point to multiple formation routes.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

We investigate the quenching properties of central and satellite galaxies, utilizing the halo masses and central-satellite identifications from the SDSS galaxy group catalog of Yang et al. We find that the quenched fractions of centrals and satellites of similar stellar masses have similar dependence on host halo mass. The similarity of the two populations is also found in terms of specific star formation rate and 4000 \AA\ break. The quenched fractions of centrals and satellites of similar masses show similar dependencies on bulge-to-total light ratio, central velocity dispersion and halo-centric distance in halos of given halo masses. The prevalence of optical/radio-loud AGNs is found to be similar for centrals and satellites at given stellar masses. All these findings strongly suggest that centrals and satellites of similar masses experience similar quenching processes in their host halos. We discuss implications of our results for the understanding of galaxy quenching.

M85 is a peculiar S0 galaxy in Virgo and is a well-known merger remnant. In this paper, we present the first spectroscopic study of globular clusters (GCs) in M85. We obtain spectra for 21 GC candidates and the nucleus of M85 using the Gemini Multi-Object Spectrograph on the Gemini North 8.1 m telescope. From their radial velocities, 20 of the GCs are found to be members of M85. We find a strong rotation signal of the M85 GC system with a rotation amplitude of 235 km s$^{-1}$. The rotation axis of the GC system has a position angle of about 161$^{\circ}$, which is 51$\rlap{.}{^\circ}$5 larger than that of the stellar light. The rotation-corrected radial velocity dispersion of the GC system is estimated to be $\sigma_{\rm r,cor} = $ 160 km s$^{-1}$. The rotation parameter $\Omega R_{\rm icor}/\sigma_{\rm r,cor}$ of the GC system is derived to be 1.47$^{+1.05}_{-0.48}$, which is one of the largest among known early-type galaxies. The ages and metallicities of the GCs, which show the same trend as the results based on Lick indices, are derived from full spectrum fitting (ULySS). About a half of the GCs are an intermediate-age population of which the mean age is $\sim$ 3.7 $\pm$ 1.9 Gyr, having a mean [Fe/H] value of --0.26. The other half are old and metal-poor. These results suggest that M85 experienced a wet merging event about 4 Gyr ago, forming a significant population of star clusters. The strong rotational feature of the GC system can be explained by an off-center major merging.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

We examine the thermal energy contents of the intergalactic medium (IGM) over three orders of magnitude in both mass density and gas temperature using thermal Sunyaev-Zel'dovich effect (tSZE). The analysis is based on {\it Planck} tSZE map and the cosmic density field, reconstructed for the SDSS DR7 volume and sampled on a grid of cubic cells of $(1h^{-1}{\rm Mpc})^3$, together with a matched filter technique employed to maximize the signal-to-noise. Our results show that the pressure - density relation of the IGM is roughly a power law given by an adiabatic equation of state, with an indication of steepening at densities higher than about $10$ times the mean density of the universe. The implied average gas temperature is $\sim 10^4\,{\rm K}$ in regions of mean density, $\rho_{\rm m} \sim {\overline\rho}_{\rm m}$, increasing to about $10^5\,{\rm K}$ for $\rho_{\rm m} \sim 10\,{\overline\rho}_{\rm m}$, and to $>10^{6}\,{\rm K}$ for $\rho_{\rm m} \sim 100\,{\overline\rho}_{\rm m}$. At a given density, the thermal energy content of the IGM is also found to be higher in regions of stronger tidal fields, likely due to shock heating by the formation of large scale structure and/or feedback from galaxies and AGNs. A comparison of the results with hydrodynamic simulations suggests that the current data can already provide interesting constraints on galaxy formation.

We present Keck/DEIMOS spectroscopy of globular clusters (GCs) around the ultra-diffuse galaxies (UDGs) VLSB-B, VLSB-D, and VCC615 located in the central regions of the Virgo cluster. We spectroscopically identify 4, 12, and 7 GC satellites of these UDGs, respectively. We find that the three UDGs have systemic velocities ($V_{sys}$) consistent with being in the Virgo cluster, and that they span a wide range of velocity dispersions, from $\sim 16$ to $\sim 47$ km/s, and high dynamical mass-to-light ratios within the radius that contains half the number of GCs ($ 407^{+916}_{-407}$, $21^{+15}_{-11}$, $60^{+65}_{-38}$, respectively). VLSB-D shows possible evidence for rotation along the stellar major axis and its $V_{sys}$ is consistent with that of the massive galaxy M84 and the center of the Virgo cluster itself. These findings, in addition to having a dynamically and spatially ($\sim 1$ kpc) off-centered nucleus and being extremely elongated, suggest that VLSB-D could be tidally perturbed. On the contrary, VLSB-B and VCC615 show no signals of tidal deformation. Whereas the dynamics of VLSB-D suggest that it has a less massive dark matter halo than expected for its stellar mass, VLSB-B and VCC615 are consistent with a $\sim 10^{12}$ M$_{\odot}$ dark matter halo. Although our samples of galaxies and GCs are small, these results suggest that UDGs may be a diverse population, with their low surface brightnesses being the result of very early formation, tidal disruption, or a combination of the two.

We derive stellar population parameters for a representative sample of ultracompact dwarfs (UCDs) and a large sample of massive globular clusters (GCs) with stellar masses $\gtrsim$ 10$^{6}$ $M_{\odot}$ in the central galaxy M87 of the Virgo galaxy cluster, based on model fitting to the Lick-index measurements from both the literature and new observations. After necessary spectral stacking of the relatively faint objects in our initial sample of 40 UCDs and 118 GCs, we obtain 30 sets of Lick-index measurements for UCDs and 80 for GCs. The M87 UCDs have ages $\gtrsim$ 8 Gyr and [$\alpha$/Fe] $\simeq$ 0.4 dex, in agreement with previous studies based on smaller samples. The literature UCDs, located in lower-density environments than M87, extend to younger ages and smaller [$\alpha$/Fe] (at given metallicities) than M87 UCDs, resembling the environmental dependence of the Virgo dE nuclei. The UCDs exhibit a positive mass-metallicity relation (MZR), which flattens and connects compact ellipticals at stellar masses $\gtrsim$ 10$^{8}$ $M_{\odot}$. The Virgo dE nuclei largely follow the average MZR of UCDs, whereas most of the M87 GCs are offset towards higher metallicities for given stellar masses. The difference between the mass-metallicity distributions of UCDs and GCs may be qualitatively understood as a result of their different physical sizes at birth in a self-enrichment scenario or of galactic nuclear cluster star formation efficiency being relatively low in a tidal stripping scenario for UCD formation. The existing observations provide the necessary but not sufficient evidence for tidally stripped dE nuclei being the dominant contributors to the M87 UCDs.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

We apply a halo-based group finder to four large redshift surveys, the 2MRS, 6dFGS, SDSS and 2dFGRS, to construct group catalogs in the low-redshift Universe. The group finder is based on that of Yang et al. but with an improved halo mass assignment so that it can be applied uniformly to various redshift surveys of galaxies. Halo masses are assigned to groups according to proxies based on the stellar mass/luminosity of member galaxies. The performances of the group finder in grouping galaxies according to common halos and in halo mass assignments are tested using realistic mock samples constructed from hydrodynamical simulations and empirical models of galaxy occupation in dark matter halos. Our group finder finds $\sim 94\%$ of the correct true member galaxies for $90-95\%$ of the groups in the mock samples; the halo masses assigned by the group finder are un-biased with respect to the true halo masses, and have a typical uncertainty of $\sim0.2\,{\rm dex}$. The properties of group catalogs constructed from the observational samples are described and compared with other similar catalogs in the literature.

A matched filter technique is applied to the Planck all-sky Compton y-parameter map to measure the thermal Sunyaev-Zel'dovich (tSZ) effect produced by galaxy groups of different halo masses selected from large redshift surveys in the low-z Universe. Reliable halo mass estimates are available for all the groups, which allows us to bin groups of similar halo masses to investigate how the tSZ effect depends on halo mass over a large mass range. Filters are simultaneously matched for all groups to minimize projection effects. We find that the integrated y-parameter and the hot gas content it implies are consistent with the predictions of the universal pressure profile model only for massive groups above $10^{14}\,{\rm M}_\odot$, but much lower than the model prediction for low-mass groups. The halo mass dependence found is in good agreement with the predictions of a set of simulations that include strong AGN feedback, but simulations including only supernova feedback significantly over predict the hot gas contents in galaxy groups. Our results suggest that hot gas in galaxy groups is either effectively ejected or in phases much below the virial temperatures of the host halos.

We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.

Small distortions in the images of Einstein rings or giant arcs offer the exciting prospect of detecting dark matter haloes or subhaloes of mass below $10^9$M$_{\odot}$, most of which are too small to have made a visible galaxy. A very large number of such haloes are predicted to exist in the cold dark matter model of cosmogony; in contrast other models, such as warm dark matter, predict no haloes below a mass of this order which depends on the properties of the warm dark matter particle. Attempting to detect these small perturbers could therefore discriminate between different kinds of dark matter particles, and even rule out the cold dark matter model altogether. Globular clusters in the lens galaxy also induce distortions in the image which could, in principle, contaminate the test. Here, we investigate the population of globular clusters in six early type galaxies in the Virgo cluster. We find that the number density of globular clusters of mass $\sim10^6$M$_{\odot}$ is comparable to that of the dark matter perturbers (including subhaloes in the lens and haloes along the line-of-sight). We show that the very different degrees of mass concentration in globular clusters and dark matter haloes result in different lensing distortions. These are detectable with milli-arcsecond resolution imaging which can distinguish between globular cluster and dark matter halo signals.

We show that the ratio between the stellar mass of central galaxy and the mass of its host halo, $f_c \equiv M_{*,c}/M_{\rm h}$, can be used as an observable proxy of halo assembly time, in that galaxy groups with higher $f_c$ assembled their masses earlier. Using SDSS groups of Yang et al., we study how $f_c$ correlates with galaxy properties such as color, star formation rate, metallicity, bulge to disk ratio, and size. Central galaxies of a given stellar mass in groups with $f_c>0.02$ tend to be redder in color, more quenched in star formation, smaller in size, and more bulge dominated, as $f_c$ increases. The trends in color and star formation appear to reverse at $f_c<0.02$, reflecting a down-sizing effect that galaxies in massive halos formed their stars earlier although the host halos themselves assembled later (lower $f_c$). No such reversal is seen in the size of elliptical galaxies, suggesting that their assembly follows halo growth more closely than their star formation. Satellite galaxies of a given stellar mass in groups of a given halo mass tend to be redder in color, more quenched in star formation and smaller in size as $f_c$ increases. For a given stellar mass, satellites also tend to be smaller than centrals. The trends are stronger for lower mass groups. For groups more massive than $\sim 10^{13}{\rm M}_\odot$, a weak reversed trend is seen in color and star formation. The observed trends in star formation are qualitatively reproduced by an empirical model based on halo age abundance matching, but not by a semi-analytical model tested here.

We present the results from a wide-field spectroscopic survey of globular clusters (GCs) in the Virgo Cluster. We obtain spectra for 201 GCs and 55 ultracompact dwarfs (UCDs) using the Hectospec on the Multiple Mirror Telescope, and derive their radial velocities. We identify 46 genuine intracluster GCs (IGCs), not associated with any Virgo galaxies, using the 3D GMM test on the spatial and radial velocity distribution.They are located at the projected distance 200 kpc $\lesssim$ R $\lesssim$ 500 kpc from the center of M87. The radial velocity distribution of these IGCs shows two peaks, one at $v_{\rm r}$ = 1023 km s$^{-1}$ associated with the Virgo main body, and another at $v_{\rm r}$ = 36 km s$^{-1}$ associated with the infalling structure. The velocity dispersion of the IGCs in the Virgo main body is $\sigma_{\rm{GC}} \sim$ 314 km s$^{-1}$, which is smoothly connected to the velocity dispersion profile of M87 GCs, but much lower than that of dwarf galaxies in the same survey field, $\sigma_{\rm{dwarf}} \sim$ 608 km s$^{-1}$. The UCDs are more centrally concentrated on massive galaxies, M87, M86, and M84. The radial velocity dispersion of the UCD system is much smaller than that of dwarf galaxies. Our results confirm the large-scale distribution of Virgo IGCs indicated by previous photometric surveys. The color distribution of the confirmed IGCs shows a bimodality similar to that of M87 GCs. This indicates that most IGCs are stripped off from dwarf galaxies and some from massive galaxies in the Virgo.

Globular clusters (GCs) are some of the most visible tracers of the merging and accretion history of galaxy halos. Metal-poor GCs, in particular, are thought to arrive in massive galaxies largely through dry, minor merging events, but it is rare to see a direct connection between GCs and visible stellar streams. NGC 474 is a post-merger early-type galaxy with dramatic fine structures made of concentric shells and radial streams that have been more clearly revealed by deep imaging. We present a study of GCs in NGC 474 to better establish the relationship between merger-induced fine structure and the GC system. We find that many GCs are superimposed on visible streams and shells, and about 35% of GCs outside $3R_{\rm e,galaxy}$ are located in regions of fine structure. The spatial correlation between the GCs and fine structure is significant at the 99.9% level, showing that this correlation is not coincidental. The colors of the GCs on the fine structures are mostly blue, and we also find an intermediate-color population that is dominant in the central region, and which will likely passively evolve to have colors consistent with a traditional metal-rich GC population. The association of the blue GCs with fine structures is direct confirmation that many metal-poor GCs are accreted onto massive galaxy halos through merging events, and that progenitors of these mergers are sub-L* galaxies.

We compare predictions of a number of empirical models and numerical simulations of galaxy formation to the conditional stellar mass functions (CSMF)of galaxies in groups of different masses obtained recently by Lan et al. to test how well different models accommodate the data. The observational data clearly prefer a model in which star formation in low-mass halos changes behavior at a characteristic redshift $z_c\sim 2$. There is also tentative evidence that this characteristic redshift depends on environment, becoming $z_c\sim 4$ in regions that eventually evolve into rich clusters of galaxies. The constrained model is used to understand how galaxies form and evolve in dark matter halos, and to make predictions for other statistical properties of the galaxy population, such as the stellar mass functions of galaxies at high $z$, the star formation and stellar mass assembly histories in dark matter halos. A comparison of our model predictions with those of other empirical models shows that different models can make vastly different predictions, even though all of them are tuned to match the observed stellar mass functions of galaxies.

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

We present an analysis of high-quality photometry for globular clusters (GCs) in the Virgo cluster core region, based on data from the Next Generation Virgo Cluster Survey (NGVS) pilot field, and in the Milky Way (MW) based on VLT/X-Shooter spectrophotometry. We find significant discrepancies in color-color diagrams between sub-samples from different environments, confirming that the environment has a strong influence on the integrated colors of GCs. GC color distributions along a single color are not sufficient to capture the differences we observe in color-color space. While the average photometric colors become bluer with increasing radial distance to the cD galaxy M87, we also find a relation between the environment and the slope and intercept of the color-color relations. A denser environment seems to produce a larger dynamic range in certain color indices. We argue that these results are not due solely to differential extinction, IMF variations, calibration uncertainties, or overall age/metallicity variations. We therefore suggest that the relation between the environment and GC colors is, at least in part, due to chemical abundance variations, which affect stellar spectra and stellar evolution tracks. Our results demonstrate that stellar population diagnostics derived from model predictions which are calibrated on one particular sample of GCs may not be appropriate for all extragalactic GCs. These results advocate a more complex model of the assembly history of GC systems in massive galaxies that goes beyond the simple bimodality found in previous decades.

Diffuse star clusters (DSCs) are old and dynamically hot stellar systems that have lower surface brightness and more extended morphology than globular clusters (GCs). Using the images from HST/ACS Fornax Cluster Survey, we find that 12 out of 43 early-type galaxies (ETGs) in the Fornax cluster host significant numbers of DSCs. Together with literature data from the HST/ACS Virgo Cluster Survey, where 18 out of 100 ETGs were found to host DSCs, we systematically study the relationship of DSCs with GCs, and their host galaxy environment. Two DSC hosts are post-merger galaxies, with most of the other hosts either having low mass or showing clear disk components. We find that while the number ratio of DSCs to GCs is nearly constant in massive galaxies, the DSC-to-GC ratio becomes systematically higher in lower mass hosts. This suggests that DSCs may be more efficient at forming (or surviving) in low density environments. DSC hosts are not special either in their position in the cluster, or in the galactic color-magnitude diagram. Why some disk and low-mass galaxies host DSCs while others do not is still a puzzle, however. The mean ages of DSC hosts and non-hosts are similar at similar masses, implying that formation efficiency, rather than survival, is the reason behind different DSC number fractions in early-type galaxies.

Observations of nearby galaxy clusters at low surface brightness have identified galaxies with low luminosities, but sizes as large as L* galaxies, leading them to be dubbed "ultra-diffuse galaxies" (UDGs). The survival of UDGs in dense environments like the Coma cluster suggests that UDGs could reside in much more massive dark halos. We report the detection of a substantial population of globular clusters (GCs) around a Coma UDG, Dragonfly 17 (DF17). We find that DF17 has a high GC specific frequency of S_N=26+/-13. The GC system is extended, with an effective radius of 12"+/-2", or 5.6+/-0.9 kpc at Coma distance, 70% larger than the galaxy itself. We also estimate the mean of the GC luminosity function to infer a distance of 97 (+17/-14) Mpc, providing redshift-independent confirmation that one of these UDGs is in the Coma cluster. The presence of a rich GC system in DF17 indicates that, despite its low stellar density, star formation was intense enough to form many massive star clusters. If DF17's ratio of total GC mass to total halo mass is similar to those in other galaxies, then DF17 has an inferred total mass of ~10^11 solar masses, only ~10% the mass of the Milky Way, but extremely dominated by dark matter, with M/L_V~1000. We suggest that UDGs like DF17 may be "pure stellar halos", i.e., galaxies that formed their stellar halo components, but then suffered an early cessation in star formation that prevented the formation of any substantial central disk or bulge.

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$\alpha$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

We present a survey of star clusters in the halo of IC 10, a starburst galaxy in the Local Group based on Subaru R band images and NOAO Local Group Survey UBVRI images. We find five new star clusters. All these star clusters are located far from the center of IC 10, while previously known star clusters are mostly in the main body. Interestingly the distribution of these star clusters shows an asymmetrical structure elongated along the east and south-west direction. We derive UBVRI photometry of 66 star clusters including these new star clusters as well as previously known star clusters. Ages of the star clusters are estimated from the comparison of their UBVRI spectral energy distribution with the simple stellar population models. We find that the star clusters in the halo are all older than 1 Gyr, while those in the main body have various ages from very young (several Myr) to old (>1 Gyr). The young clusters (<10 Myr) are mostly located in the H{\alpha} emission regions and are concentrated on a small region at 2' in the south-east direction from the galaxy center, while the old clusters are distributed in a wider area than the disk. Intermediate-age clusters (~100 Myr) are found in two groups. One is close to the location of the young clusters and the other is at ~4' from the location of the young clusters. The latter may be related with past merger or tidal interaction.

We present spectroscopy of 28 SNR candidates as well as one H II region in M81, and two SNR candidates in M82. Twenty six out of the M81 candidates turn out to be genuine SNRs, and two in M82 may be shocked condensations in the galactic outflow or SNRs. The distribution of [N II]/H{\alpha} ratios of M81 SNRs is bimodal. M81 SNRs are divided into two groups in the spectral line ratio diagrams: an [O III]-strong group and an [O III]-weak group. The latter have larger sizes, and may have faster shock velocity. [N II]/H{\alpha} ratios of the SNRs show a strong correlation with [S II]/H{\alpha} ratios. They show a clear radial gradient in [N II]/H{\alpha} and [S II]/H{\alpha} ratios: dLog ([N II]/H{\alpha})/dLog R = -0.018 {\pm} 0.008 dex/kpc and dLog ([S II]/H{\alpha})/dLog R = -0.016 {\pm} 0.008 dex/kpc where R is a deprojected galactocentric distance. We estimate the nitrogen and oxygen abundance of the SNRs from the comparison with shock-ionization models. We obtain a value for the nitrogen radial gradient, dLog(N/H)/dLogR = -0.023 {\pm} 0.009 dex/kpc, and little evidence for the gradient in oxygen. This nitrogen abundance shows a few times flatter gradient than those of the planetary nebulae and H II regions. We find that five SNRs are matched with X-ray sources. Their X-ray hardness colors are consistent with thermal SNRs.

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

Previous measurements of the composition of Ultra-High Energy Cosmic Rays(UHECRs) made by the High Resolution Fly's Eye(HiRes) and Pierre Auger Observatory(PAO) are seemingly contradictory, but utilize different detection methods, as HiRes was a stereo detector and PAO is a hybrid detector. The five year Telescope Array(TA) Middle Drum hybrid composition measurement is similar in some, but not all, respects in methodology to PAO, and good agreement is evident between data and a light, largely protonic, composition when comparing the measurements to predictions obtained with the QGSJetII-03 and QGSJet-01c models. These models are also in agreement with previous HiRes stereo measurements, confirming the equivalence of the stereo and hybrid methods. The data is incompatible with a pure iron composition, for all models examined, over the available range of energies. The elongation rate and mean values of Xmax are in good agreement with Pierre Auger Observatory data. This analysis is presented using two methods: data cuts using simple geometrical variables and a new pattern recognition technique.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We have searched for intermediate-scale anisotropy in the arrival directions of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky using data collected over a 5 year period by the surface detector of the Telescope Array experiment. We report on a cluster of events that we call the hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has a Li-Ma statistical significance of 5.1$\sigma$, and is centered at R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about 19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of events of 5.1$\sigma$ significance, appearing by chance in an isotropic cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$\sigma$).

We present new formulae for the mass functions of the clusters and the isolated clusters with non Gaussian initial conditions. For this study, we adopt the Extended Zel'dovich (EZL) model as a basic framework, focusing on the case of primordial non-Gaussianity of the local type whose degree is quantified by a single parameter $f_{nl}$. By making a quantitative comparison with the N-body results, we first demonstrate that the EZL formula with the constant values of three fitting parameters still works remarkably well for the local $f_{nl}$ case. We also modify the EZL formula to find an analytic expression for the mass function of isolated clusters which turns out to have only one fitting parameter other than the overall normalization factor and showed that the modified EZL formula with a constant value of the fitting parameter matches excellently the N-body results with various values of $f_{nl}$ at various redshifts. Given the simplicity of the generalized EZL formulae and their good agreements with the numerical results, we finally conclude that the EZL mass functions of the massive clusters and isolated clusters should be useful as an analytic guideline to constrain the scale dependence of the primordial non-Gaussianity of the local type.

We present an analytic formula for the supercluster mass function which is constructed by modifying the extended Zel'dovich model for the halo mass function. The formula has two characteristic parameters whose best-fit values are determined by fitting to the numerical results from N-body simulations for the standard LambdaCDM cosmology. The parameters are found to be independent of redshifts and robust against variation of the key cosmological parameters. Under the assumption that the same formula for the supercluster mass function is valid for non-standard cosmological models, we show that the relative abundance of the rich superclusters should be a powerful indicator of any deviation of the real universe from the prediction of the standard LambdaCDM model.

We present a spectroscopic study of the four extended star clusters (ESCs) in NGC 6822 based on the data obtained with Gemini Multi-Object Spectrograph (GMOS) on the Gemini-South 8.1m telescope. Their radial velocities derived from the spectra range from $-61.2 \pm 20.4$ km s$^{-1}$ (for C1) to $-115.34 \pm 57.9$ km s$^{-1}$ (for C4) and, unlike the intermediate age carbon stars, they do not display any sign of systematic rotation around NGC 6822. The ages and metallicities derived using the Lick indices show that the ESCs are old ($\geq 8$ Gyr) and metal poor ([Fe/H]$\lesssim -1.5$). NGC 6822 is found to have both metal poor ([Fe/H]$\approx -2.0$) and metal rich ([Fe/H]$\approx -0.9$) star clusters within 15 arcmin (2 kpc) from the center, while only metal poor clusters are observed in the outer halo with $r \geq 20$ arcmin (2.6 kpc). The kinematics, old ages, and low metallicities of ESCs suggest that ESCs may have accreted into the halo of NGC 6822. Based on the velocity distribution of ESCs, we have determined the total mass and the mass-to-light ratio of NGC 6822: $M_{N6822} = 7.5^{+4.5}_{-0.1} \times 10^{9} M_{\odot}$ and $(M/L)_{N6822} = 75^{+45}_{-1} (M/L)_{\odot}$. It shows that NGC 6822 is one of the most dark matter dominated dwarf galaxies in the Local Group.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

We present a study of Ks band luminosity evolution of the asymptotic giant branch (AGB) population in simple stellar systems using star clusters in the Large Magellanic Cloud (LMC). We determine physical parameters of LMC star clusters including center coordinates, radii, and foreground reddenings. Ages of 83 star clusters are derived from isochrone fitting with the Padova models, and those of 19 star clusters are taken from the literature. The AGB stars in 102 star clusters with log(age) = 7.3 - 9.5 are selected using near-infrared color magnitude diagrams based on 2MASS photometry. Then we obtain the Ks band luminosity fraction of AGB stars in these star clusters as a function of ages. The Ks band luminosity fraction of AGB stars increases, on average, as age increases from log(age) ~ 8.0, reaching a maximum at log(age) ~ 8.5, and it decreases thereafter. There is a large scatter in the AGB luminosity fraction for given ages, which is mainly due to stochastic effects. We discuss this result in comparison with five simple stellar population models. The maximum Ks band AGB luminosity fraction for bright clusters is reproduced by the models that expect the value of 0.7 - 0.8 at log(age) = 8.5 - 8.7. We discuss the implication of our results with regard to the study of size and mass evolution of galaxies.

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.

We present a measurement of the energy spectrum of ultra-high-energy cosmic rays performed by the Telescope Array experiment using monocular observations from its two new FADC-based fluorescence detectors. After a short description of the experiment, we describe the data analysis and event reconstruction procedures. Since the aperture of the experiment must be calculated by Monte Carlo simulation, we describe this calculation and the comparisons of simulated and real data used to verify the validity of the aperture calculation. Finally, we present the energy spectrum calculated from the merged monocular data sets of the two FADC-based detectors, and also the combination of this merged spectrum with an independent, previously published monocular spectrum measurement performed by Telescope Array's third fluorescence detector (Abu-Zayyad {\it et al.}, {Astropart. Phys.} 39 (2012), 109). This combined spectrum corroborates the recently published Telescope Array surface detector spectrum (Abu-Zayyad {\it et al.}, {Astrophys. Journ.} 768 (2013), L1) with independent systematic uncertainties.

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays with primary energies above 1.6 x 10^(18) eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 x 10^(18) eV and a steepening at 5.4 x 10^(19) eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of ultra-high energy cosmic ray surface detector data, that involves generating a complete simulation of ultra-high energy cosmic rays striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the "thinning" approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

We present a photometric study of star clusters in the nearby starburst galaxy M82 based on the UBVI, YJ and H band Hubble Space Telescope images. We find 1105 star clusters with V<23 mag. 1070 of them are located in the disk region, while 35 star clusters are in the halo region. The star clusters in the disk are composed of a dominant blue population with a color peak at (B-V)0 ~ 0.45, and a weaker red population. The luminosity function of the disk clusters shows a power-law distribution with a power law index alpha=-2.04+-0.03, and the scale height of their distribution is h_z=9".64+-0".40 (164 +- 7 pc), similar to that of the stellar thin disk of M82. We have derived the ages of ~630 star clusters using the spectral energy distribution fit method by comparing UBVI(YJ)H band photometric data with the simple stellar population models. The age distribution of the disk clusters shows that the most dominant cluster population has ages ranging from 100 Myr to 1 Gyr, with a peak at about 500 Myr. This suggests that M82 has undergone a disk-wide star formation about 500 Myr ago, probably through the interaction with M81. The brightest star clusters in the nuclear region are much brighter than those in other regions, indicating that more massive star clusters are formed in the denser environments. On the other hand, the colors of the halo clusters are similar to those of globular clusters in the Milky Way Galaxy, and their ages are estimated to be older than 1 Gyr. These are probably genuine old globular clusters in M82.

A new way to construct a fitting formula for the halo mass function is presented. Our formula is expressed as a solution to the modified Jedamzik matrix equation that automatically satisfies the normalization constraint. The characteristic parameters expressed in terms of the linear shear eigenvalues are empirically determined by fitting the analytic formula to the numerical results from the high-resolution N-body simulation and found to be independent of scale, redshift and background cosmology. Our fitting formula with the best-fit parameters is shown to work excellently in the wide mass-range at various redshifts: The ratio of the analytic formula to the N-body results departs from unity by up to 10% and 5% over 10^{11}<= M/(M_sun/h)<= 5x10^{15} at z=0,\ 0.5 and 1 for the FoF-halo and SO-halo cases, respectively.

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the distributions of the events in right ascension and declination are compatible with an isotropic distribution in all three sets. We then compare with previously reported clustering of the UHECR events at small angular scales. No significant clustering is found in the TA data. We then check the events with E>57 EeV for correlations with nearby active galactic nuclei. No significant correlation is found. Finally, we examine all three sets for correlations with the large-scale structure of the Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS hypothesis), while the event set with E>10 EeV is compatible with isotropy and is not compatible with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this set compatible with the LSS hypothesis.

The Telescope Array (TA) experiment, located in the western desert of Utah,USA, is designed for observation of extensive air showers from extremely high energy cosmic rays. The experiment has a surface detector array surrounded by three fluorescence detectors to enable simultaneous detection of shower particles at ground level and fluorescence photons along the shower track. The TA surface detectors and fluorescence detectors started full hybrid observation in March, 2008. In this article we describe the design and technical features of the TA surface detector.

We report the discovery of two new globular clusters in the remote halos of M81 and M82 in the M81 Group based on Hubble Space Telescope archive images. They are brighter than typical globular clusters (MV = -9.34 mag for GC-1 and M_V = -10.51 mag for GC-2), and much larger than known globular clusters with similar luminosity in the MilkyWay Galaxy and M81. Radial surface brightness profiles for GC-1 and GC-2 do not show any features of tidal truncation in the outer part. They are located much farther from both M81 and M82 in the sky, compared with previously known star clusters in these galaxies. Color-magnitude diagrams of resolved stars in each cluster show a well-defined red giant branch (RGB), indicating that they are metal-poor and old. We derive a low metallicity with [Fe/H] $\simeq -2.3$ and an old age ~14 Gyr for GC-2 from the analysis of the absorption lines in its spectrum in the Sloan Digital Sky Survey in comparison with the simple stellar population models. The I-band magnitude of the tip of the RGB for GC-2 is 0.26 mag fainter than that for the halo stars in the same field, showing that GC-2 is ~400 kpc behind the M81 halo along our line of sight. The deprojected distance to GC-2 from M81 is much larger than any other known globular clusters in the local universe. This shows that GC-2 is the most isolated globular cluster in the local universe.

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

Since 2007, the Telescope Array (TA) experiment, based in Utah, USA, has been observing ultra high energy cosmic rays to understand their origins. The experiment involves a surface detector (SD) array and three fluorescence detector (FD) stations. FD stations, installed surrounding the SD array, measure the air fluorescence light emitted from extensive air showers (EASs) for precise determination of their energies and species. The detectors employed at one of the three FD stations were relocated from the High Resolution Fly's Eye experiment. At the other two stations, newly designed detectors were constructed for the TA experiment. An FD consists of a primary mirror and a camera equipped with photomultiplier tubes. To obtain the EAS parameters with high accuracies, understanding the FD optical characteristics is important. In this paper, we report the characteristics and installation of new FDs and the performances of the FD components. The results of the monitored mirror reflectance during the observation time are also described in this report.

We present the mid-infrared (MIR) observation of a nearby galaxy cluster, Abell 2255 by the AKARI space telescope. Using the AKARI's continuous wavelength coverage between 3-24 micron and the wide field of view, we investigate the properties of cluster member galaxies to see how the infall of the galaxies, the cluster substructures, and the cluster-cluster merger influence their evolution. We show that the excess of MIR (11 micron) flux is a good indicator to discriminate galaxies at different evolutionary stages, and divide galaxies into three classes accordingly : strong MIR-excess (N3-S11>0.2) galaxies that include both unobscured and obscured star-forming galaxies, weak MIR-excess (-2.0<N3-S11<-1.2) galaxies that are quiescent, old (>5 Gyr) galaxies where the MIR emission arises mainly from the circumstellar dust around AGB stars, and intermediate MIR-excess (-1.2<N3-S11<0.2) galaxies in between the two classes that are less than a few Gyrs old past the prime star formation activity. With the MIR-excess diagnostics, we investigate how local and cluster-scale environments affect the individual galaxies. We derive the total star formation rate of ~130 Msun/yr for A2255 using the strong MIR-excess galaxies, which is consistent with other clusters of galaxies at similar redshifts and with similar masses. We find no strong evidence that supports enhanced star formation neither inside the cluster nor in the substructure region. The intermediate MIR-excess galaxies, representing galaxies in transition from star-forming galaxies to quiescent galaxies, are located preferentially at the medium density region or cluster substructures. Our findings suggest that galaxies are being transformed from star-forming galaxies into red, quiescent galaxies from the infall region through near the core, which can be well-explained by the ram-pressure stripping as previous simulation suggests.

In the last decade we witness an advent of new types of dwarf stellar systems in cluding ultra-compact dwarfs, ultra-faint dwarf spheroidals, and exotic globular clusters, breaking the old simple paradigm for dwarf galaxies and globular clusters. These objects become more intriguing, and understanding of these new findings be comes more challenging. Recently we discovered a new type of large scale structure in the Virgo cluster of galaxies: it is composed of globular clusters. Globular clusters in Virgo are found wandering between galaxies (intracluster globular clusters) as well as in galaxies. These intracluster globular clusters fill a significant fraction in the area of the Virgo cluster and they are dominated by blue globular clusters. These intracluster globular clusters may be closely related with the first dwarf galaxies in the universe.

Many problems of cold dark matter models such as the cusp problem and the missing satellite problem can be alleviated, if galactic halo dark matter particles are ultra-light scalar particles and in Bose-Einstein condensate (BEC), thanks to a characteristic length scale of the particles. We show that this finite length scale of the dark matter can also explain the recently observed common central mass of the Milky Way satellites ($\sim 10^7 M_\odot$) independent of their luminosity, if the mass of the dark matter particle is about $10^{-22} eV$.

We present the {\it AKARI} InfraRed Camera (IRC) imaging observation of early-type galaxies in A2218 at z $\simeq$ 0.175. Mid-infrared (MIR) emission from early-type galaxies traces circumstellar dust emission from AGB stars or/and residual star formation. Including the unique imaging capability at 11 and 15 $\mu$m, our {\it AKARI} data provide an effective way to investigate MIR properties of early-type galaxies in the cluster environment. Among our flux-limited sample of 22 red sequence early-type galaxies with precise dynamical and line strength measurements ($<$ 18 mag at 3 $\mu m$), we find that at least 41% have MIR-excess emission. The $N3-S11$ versus $N3$ (3 and 11 $\mu$m) color-magnitude relation shows the expected blue sequence, but the MIR-excess galaxies add a red wing to the relation especially at the fainter end. A SED analysis reveals that the dust emission from AGB stars is the most likely cause for the MIR-excess, with low level of star formation being the next possible explanation. The MIR-excess galaxies show a wide spread of $N3-S11$ colors, implying a significant spread (2--11 Gyr) in the estimated mean ages of stellar populations. We study the environmental dependence of MIR-excess early-type galaxies over an area out to a half virial radius ($\sim$1 Mpc). We find that the MIR-excess early-type galaxies are preferentially located in the outer region. From these evidences, we suggest that the fainter, MIR-excess early-type galaxies have just joined the red sequence, possibly due to the infall and subsequent morphological/spectral transformation induced by the cluster environment.

We suggest that the dark matter model based on Bose Einstein condensate or scalar field can resolve the apparently contradictory behaviors of dark matter in the Abell 520 and the Bullet cluster. During a collision of two galaxies in the cluster, if initial kinetic energy of the galaxies is large enough, two dark matter halos pass each other in a soliton-like way as observed in the Bullet cluster. If not, the halos merge due to the tiny repulsive interaction among dark matter particles as observed in the Abell 520. This idea can also explain the origin of the dark galaxy and the galaxy without dark matter.