Gamma-ray bursts (GRBs) are known to be the most violent explosions in the universe, and a variety of correlations between observable GRB properties have been proposed in literature, but none of these correlations is valid for both long GRBs and short GRBs. In this paper we report the discovery of a universal correlation which is suitable for both long and short GRBs using three prompt emission properties of GRBs, i.e. the isotropic peak luminosity $L_{\rm iso}$, the peak energy of the time-integtated prompt emission spectrum $E_{\rm peak}$, and the "high signal" timescale $T_{\rm 0.45}$, $L_{\rm iso} \propto E_{\rm peak}^{1.94} T_{0.45}^{0.37}$. This universal correlation just involves properties of GRB prompt emission and does not require any information of afterglow phase, which can be used as a relatively unbiased redshift estimator. Here we use this correlation to estimate the pseudo-redshifts for short Gamma Ray Bursts and then use Lynden-Bell method to obtain a non-parametric estimate of their luminosity function and formation rate. The luminosity function is $\psi(L_0)\propto{L_0^{-0.63\pm{0.07}}}$ for dim SGRBs and $\psi(L_0)\propto{L_0^{-1.96\pm{0.28}}}$ for bright SGRBs, with the break point $6.95_{-0.76}^{+0.84}\times10^{50} erg/s$. The local formation rate of SGRBs is about 15 events $\rm Gpc^{-3}yr^{-1}$ . This universal correlation may have important implications for GRB physics, implying that the long and short GRBs should share similar radiation processes.

Luminosity of GRB 170817A is much lower than that of other sGRBs. The measurement of the superluminal movement of the radio afterglow emission confirms the presence of the relativistic jet, and the emission features can be well explained by the structured jet model. In this paper, we calculate the luminosity distribution of sGRBs and its evolution with redshift based on the structured (Gaussian) jet model, and find that the typical luminosity increase with redshift, for nearby sGRBs (such as for luminosity distance less than 200 Mpc) the typical gamma-ray luminosity is just around 10^47-10^48 erg s-1, which naturally explains the very low radiation luminosity of GRB 170817A. We derived the detection probability of sGRBs by Fermi-GBM and found that the expected detection rate of sGRBs is only about 1 yr-1 within the distance of several hundred Mpc. We explored the effect of the power-law index {\alpha} of the merger time distribution on the observed characteristics and found that it had little effect on the observed luminosity and viewing-angle distributions. However, it is very interesting that, for different values of {\alpha}, the distributions of the number of observed sGRBs are quite different, so it is possible to determine the value of {\alpha} through observed distributions of the number of sGRBs. We used the Bayesian method to make a quantitative analysis and found that the value of {\alpha} may be identified when the number of observed sGRBs with known redshifts is more than 200. Finally, we compare our results of gamma-ray luminosity distribution with sGRBs with known redshifts, and found that our results are consistent with the observation, which implies that our simulation results can reproduce the observed luminosity distribution well.

We have updated the Munich galaxy formation model, L-Galaxies, to follow the radial distributions of stars and atomic and molecular gas in galaxy discs. We include an H2-based star-formation law, as well as a detailed chemical-enrichment model with explicit mass-dependent delay times for SN-II, SN-Ia and AGB stars. Information about the star formation, feedback and chemical-enrichment histories of discs is stored in 12 concentric rings. The new model retains the success of its predecessor in reproducing the observed evolution of the galaxy population, in particular, stellar mass functions and passive fractions over the redshift range 0<=z<=3 and mass range 8<=log(M_*/Msun)<=12, the black hole-bulge mass relation at z=0, galaxy morphology as a function of stellar mass and the mass-metallicity relations of both stellar and gas components. In addition, its detailed modelling of the radial structure of discs allows qualitatively new comparisons with observation, most notably with the relative sizes and masses of the stellar, atomic and molecular components in discs. Good agreement is found with recent data. Comparison of results obtained for simulations differing in mass resolution by more than two orders of magnitude shows that all important distributions are numerically well converged even for this more detailed model. An examination of metallicity and surface-density gradients in the stars and gas indicates that our new model, with star formation, chemical enrichment and feedback calculated self-consistently on local disc scales, reproduces some but not all of the trends seen in recent many-galaxy IFU surveys.

We use semi-analytic galaxy catalogs based on two high-resolution cosmological $N$-body simulations, Millennium-WMAP7 and Millennium-II, to investigate the formation of the Local Group (LG) analogs. Unlike previous studies, we use the observed stellar masses to select the LG member (Milky Way (MW) and M31) analogs, and then impose constrains using the observed separation, isolation, and kinematics of the two main member galaxies. By comparing radial and low-ellipticity orbits between the MW and M31, we find higher tangential velocity results in higher total mass, which are 4.4$^{+2.4}_{-1.5}\times$10$^{12}\rm M_{\odot}$ and 6.6$^{+2.7}_{-1.5}\times$10$^{12}\rm M_{\odot}$ for radial and low-ellipticity orbits. The orbits also influence the individual mass distribution of MW and M31 analogs. For radial orbits, the typical host halo masses of the MW and M31 are 1.5$^{+1.4}_{-0.7}\times$10$^{12}\rm M_{\odot}$ and 2.5$^{+1.3}_{-1.1}\times$10$^{12}\rm M_{\odot}$; for low-ellipticity orbits, the masses are 2.5$^{+2.2}_{-1.4}\times$10$^{12}\rm M_{\odot}$ and 3.8$^{+2.8}_{-1.8}\times$10$^{12} \rm M_{\odot}$. The LG is located primarily in filaments with tails extending toward higher densities up to $\delta\sim4.5$. The dark matter velocity anisotropy parameters $\beta$ of both the MW and M31 analogs are close to zero in the center, increasing to 0.2--0.3 at 50--80 kpc and decreasing slowly outward. The slope is much flatter than computed from the MW satellites, and the amplitude is smaller than traced by halo stars. Values of $\beta$ from different tracers agree at $\sim$120 kpc where $\beta \sim$ 0.2. We also find that model predictions agree broadly with observations in the radial distribution and luminosity function of satellites around the MW and M31.

In the standard cosmological model, dark matter drives the structure formation and constructs potential wells within which galaxies may form. The baryon fraction in dark halos can reach the universal value (15.7%) in massive clusters and decreases rapidly as the mass of the system decreases. The formation of dwarf galaxies is sensitive both to baryonic processes and the properties of dark matter owing to the shallow potential wells in which they form. In dwarf galaxies in the Local Group, dark matter dominates the mass content even within their optical-light half-radii (r_e ~ 1 kpc). However, recently it has been argued that not all dwarf galaxies are dominated by dark matter. Here we report 19 dwarf galaxies that could consist mainly of baryons up to radii well beyond r_e, at which point they are expected to be dominated by dark matter. Of these, 14 are isolated dwarf galaxies, free from the influence of nearby bright galaxies and high dense environments. This result provides observational evidence that could challenge the formation theory of low-mass galaxies within the framework of standard cosmology. Further observations, in particular deep imaging and spatially-resolved kinematics, are needed to constrain the baryon fraction better in such galaxies.

We investigate the formation of ultra-diffuse galaxies (UDGs) using the Auriga high-resolution cosmological magneto-hydrodynamical simulations of Milky Way-sized galaxies. We identify a sample of $92$ UDGs in the simulations that match a wide range of observables such as sizes, central surface brightness, S\'{e}rsic indices, colors, spatial distribution and abundance. Auriga UDGs have dynamical masses similar to normal dwarfs. In the field, the key to their origin is a strong correlation present in low-mass dark matter haloes between galaxy size and halo spin parameter. Field UDGs form in dark matter haloes with larger spins compared to normal dwarfs in the field, in agreement with previous semi-analytical models. Satellite UDGs, on the other hand, have two different origins: $\sim 55\%$ of them formed as field UDGs before they were accreted; the remaining $\sim 45\%$ were normal field dwarfs that subsequently turned into UDGs as a result of tidal interactions.

The massive galaxy population above the characteristic Schechter mass M*~10^10.6 Msun contributes to about half of the total stellar mass in the local Universe. These massive galaxies usually reside in hot dark matter haloes above the critical shock-heating mass ~10^12 Msun, where the external cold gas supply to these galaxies is expected to be suppressed. When galaxies run out of their cold gas reservoir, they become dead and quiescent. Therefore, massive quiescent galaxies living in hot haloes are commonly believed to be gas-poor. Based on the data from SDSS, ALFALFA, GASS and COLD GASS surveys, here we show that the vast majority of the massive, quiescent, central disk galaxies in the nearby Universe have a remarkably large amount of cold atomic hydrogen gas, surprisingly similar to star-forming galaxies. Both star-forming and quiescent disk galaxies show identical symmetric double-horn HI spectra, indicating similar regularly rotating HI disks. Relative to their star-forming counterparts, massive quiescent central disk galaxies are quenched because of their significantly reduced molecular gas content, lower dust content, and lower star formation efficiency. Our findings reveal a new picture, which clearly demonstrates the detailed star-formation quenching process in massive galaxies and provides a stringent constraint on the physical mechanism of quenching.

We introduce a new halo/subhalo finder, HIKER (a Halo fInder based on KERnel-shift algorithm), which takes advantage of a machine learning method -- the mean-shift algorithm combined with the Plummer kernel function, to effectively locate density peaks corresponding to halos/subhalos in density field. Based on these density peaks, dark matter halos are identified as spherical overdensity structures, and subhalos are bound substructures with boundaries at their tidal radius. By testing HIKER code with mock halos, we show that HIKER performs excellently in recovering input halo properties. Especially, HIKER has higher accuracy in locating halo/subhalo centres than most halo finders. With cosmological simulations, we further show that HIKER reproduces the abundance of dark matter halos and subhalos quite accurately, and the HIKER halo/subhalo mass functions and $V_{max}$ functions are in good agreement with two widely used halo finders, SUBFIND and AHF.

Low mass galaxies are expected to be dark matter dominated even within their centrals. Recently two observations reported two dwarf galaxies in group environment with very little dark matter in their centrals. We explore the population and origins of dark-matter-deficient galaxies (DMDGs) in two state-of-the-art hydrodynamical simulations, the EAGLE and Illustris projects. For all satellite galaxies with $10^9<M_*<10^{10}$ M$_{\odot}$ in groups with $M_{200}>10^{13}$ M$_{\odot}$, we find that about $2.6\%$ of them in the EAGLE, and $1.5\%$ in the Illustris are DMDGs with dark matter fractions below $50\%$ inside two times half-stellar-mass radii. We demonstrate that DMDGs are highly tidal disrupted galaxies; and because dark matter has higher binding energy than stars, mass loss of the dark matter is much more rapid than stars in DMDGs during tidal interactions. If DMDGs were confirmed in observations, they are expected in current galaxy formation models.

It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $\Lambda$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $\Lambda$CDM simulations indicating that such a set up is not inconsistent with $\Lambda$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.

Space very long baseline interferometry (VLBI) has unique applications in high-resolution imaging of fine structure of astronomical objects and high-precision astrometry due to the key long space-Earth or space-space baselines beyond the Earth's diameter. China has been actively involved in the development of space VLBI in recent years. This review briefly summarizes China's research progress in space VLBI and the future development plan.

We developed a new semi-analytic galaxy formation model: Galaxy Assembly with Binary Evolution (GABE). For the first time we introduce binary evolution into semi-analytic models of galaxy formation by using Yunnan-II stellar population synthesis model, which includes various binary interactions. When implementing our galaxy formation model onto the merger trees extracted from the Millennium simulation, it can reproduce a large body of observational results. We find that in the local universe the model including binary evolution reduces the luminosity at optical and infrared wavelengths slightly, while increases the luminosity at ultraviolet wavelength significantly, especially in $F_{\rm UV}$ band. The resulting luminosity function does not change very much over SDSS optical bands and infrared band, but the predicted colors are bluer, especially when $F_{\rm UV}$ band is under consideration. The new model allows us to explore the physics of various high energy events related to the remnants of binary stars, e.g. type Ia supernovae, short gamma-ray bursts and gravitational wave events, and their relation with host galaxies in a cosmological context.

w-Projection is a wide-field imaging technique that is widely used in radio synthesis arrays. Processing the wide-field big data generated by the future Square Kilometre Array (SKA) will require significant updates to current methods to significantly reduce the time consumed on data processing. Data loading and gridding are found to be two major time-consuming tasks in w-projection. In this paper, we investigate two parallel methods of accelerating w-projection processing on multiple nodes: the hybrid Message Passing Interface (MPI) and Open Multi-Processing (OpenMP) method based on multicore Central Processing Units (CPUs) and the hybrid MPI and Compute Unified Device Architecture (CUDA) method based on Graphics Processing Units (GPUs). Both methods are successfully employed and operated in various computational environments, confirming their robustness. The experimental results show that the total runtime of both MPI + OpenMP and MPI + CUDA methods is significantly shorter than that of single-thread processing. MPI + CUDA generally shows faster performance when running on multiple nodes than MPI + OpenMP, especially on large numbers of nodes. The single-precision GPU-based processing yields faster computation than the double-precision processing; while the single- and doubleprecision CPU-based processing shows consistent computational performance. The gridding time remarkably increases when the support size of the convolution kernel is larger than 8 and the image size is larger than 2,048 pixels. The present research offers useful guidance for developing SKA imaging pipelines.

It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called `shape alignment' has been observed in isolated Milky Way-like galaxies but not in `Local Group'-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the {\it facing} region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the {\it away} region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at $\sim5\sigma$ levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

We compare galaxy morphology predicted by the Illustris hydrodynamical simulation and a semi-analytic model (SAM) grafted in the halo merger trees from the Illustris-Dark matter simulation. Morphology is classified according to the luminous profile and the kinematic bulge-to-total ratio for Illustris galaxy, and the bulge-to-total stellar mass ratio for SAM galaxy. For late-type galaxies in the Illustris catalogue, most of their counterparts in the SAM model have the same type, and the consistency between two models is higher for lower mass galaxies. For early-type (ET) galaxies in Illustris, the consistency is quite low for the counterparts except for most massive galaxies. By comparing in detail the growth histories of some matched galaxy pairs of Milky Way mass in Illustris and the SAM model, we notice two aspects of differences in determining galaxy morphology between the two galaxy formation implementations. First, in the SAM, major merger and frequent minor mergers result in the growth of bulges and turn the galaxy into ET, while bulge formation is not connected to mergers as tightly as in SAM for the Illustris galaxies. In addition, the satellite stellar mass can decrease significantly due to tidal stripping before merging into the central galaxy in Illustris, while it does not decrease in the SAM model. This results in less mergers with large (stellar) mass ratios in the Illustris simulation, and less effect of mergers on shaping galaxy morphology.

The Square Kilometre Array (SKA) will be the first low-frequency instrument with the capability to directly image the structures of the Epoch of Reionization (EoR). Indeed, deep imaging of the EoR over 5 targeted fields of 20 square degrees each has been selected as the highest priority science objective for SKA1, probably the most exciting scientific outcome to be expected for the SKA1-Low. Aiming at preparing for this highly challenging observation, we perform an extensive pre-selection of the 'quietest' and 'cleanest' candidate fields in the southern sky to be suited for deep imaging of the EoR using existing catalogs and observations over a broad band. The candidate fields should meet a number of strict criteria to avoid contamination from foreground structures and sources such as the Galactic plane, LMC/SMC, clusters of galaxies and bright radio sources. The candidate fields should also exhibit both the lowest average surface brightness and smallest variance to ensure uniformity and high quality deep imaging over the fields. We use two sets of field of view in our selection: the exactly same size of 20 square degrees as defined in current SKA1 documentation and 40 square degrees - an area of twice big as the EoR imaging fields. The choice of a larger area allows us to maximally reduce the environmental effect from radio structures/sources in the vicinity of the candidate field. Our selection yields a sample of 29 (9) `ideal' fields of 20 (40) square degrees in the southern sky that could be targeted for deep imaging of the EoR. We have launched a campaign to observe these candidate fields with some of the SKA pathfinders to set further constraints on these candidates. These pre-selected candidate fields can provide a valuable guidance for eventual selection of the five targeted fields for imaging the EoR structures with the future SKA1-low.

The correlation between the spins of dark matter halos and the large-scale structure (LSS) has been studied in great detail over a large redshift range, while investigations of galaxies are still incomplete. Motivated by this point, we use the state-of-the-art hydrodynamic simulation, Illustris-1, to investigate mainly the spin--LSS correlation of galaxies at redshift of $z=0$. We mainly find that the spins of low-mass, blue, oblate galaxies are preferentially aligned with the slowest collapsing direction ($e_3$) of the large-scale tidal field, while massive, red, prolate galaxy spins tend to be perpendicular to $e_3$. The transition from a parallel to a perpendicular trend occurs at $\sim10^{9.4} M_{\odot}/h$ in the stellar mass, $\sim0.62$ in the g-r color, and $\sim0.4$ in triaxiality. The transition stellar mass decreases with increasing redshifts. The alignment was found to be primarily correlated with the galaxy stellar mass. Our results are consistent with previous studies both in N-body simulations and observations. Our study also fills the vacancy in the study of the galaxy spin--LSS correlation at $z=0$ using hydrodynamical simulations and also provides important insight to understand the formation and evolution of galaxy angular momentum.

The Chinese Space Station Optical Survey (CSS-OS) is a major science project of the Space Application System of the China Manned Space Program. This survey is planned to perform both photometric imaging and slitless spectroscopic observations, and it will focus on different cosmological and astronomical goals. Most of these goals are tightly dependent on the accuracy of photometric redshift (photo-z) measurement, especially for the weak gravitational lensing survey as a main science driver. In this work, we assess if the current filter definition can provide accurate photo-z measurement to meet the science requirement. We use the COSMOS galaxy catalog to create a mock catalog for the CSS-OS. We compare different photo-z codes and fitting methods that using the spectral energy distribution (SED) template-fitting technique, and choose to use a modified LePhare code in photo-z fitting process. Then we investigate the CSS-OS photo-z accuracy in certain ranges of filter parameters, such as band position, width, and slope. We find that the current CSS-OS filter definition can achieve reasonably good photo-z results with sigma_z~0.02 and outlier fraction ~3%.

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system.

We use a volume-limited galaxy sample from the SDSS Data Release 7 to explore the dependence of galactic conformity on the large-scale environment, measured on $\sim$ 4 Mpc scales. We find that the star formation activity of neighbour galaxies depends more strongly on the environment than on the activity of their primary galaxies. In under-dense regions most neighbour galaxies tend to be active, while in over-dense regions neighbour galaxies are mostly passive, regardless of the activity of their primary galaxies. At a given stellar mass, passive primary galaxies reside in higher density regions than active primary galaxies, leading to the apparently strong conformity signal. The dependence of the activity of neighbour galaxies on environment can be explained by the corresponding dependence of the fraction of satellite galaxies. Similar results are found for galaxies in a semi-analytical model, suggesting that no new physics is required to explain the observed large-scale conformity.

We present our photometric monitoring of a flat spectrum radio quasar (FSRQ) 3C 454.3 at Yunnan observatories from 2006 to 2011. We find that the optical color of 3C 454.3 shows obvious redder-when-brighter trend, which reaches a saturation stage when the source is brighter than 15.15 mag at V band. We perform a simulation with multiple values of disk luminosity and spectral index to reproduce the magnitude-color diagram. The results show that the contamination caused by the disk radiation alone is difficult to produce the observed color variability. The variability properties during the outburst in December 2009 are also compared with $\gamma$-ray data derived from Fermi $\gamma$-ray space telescope. The flux variation of these two bands follow a linear relation with $F_{\gamma} \propto F_R^{1.14\pm0.07}$, which provides an observational evidence for external Compton process in 3C 454.3. Meanwhile, this flux correlation indicates that electron injection is the main mechanism for variability origin. We also explore the variation of the flux ratio $F_{\gamma}/F_R$ and the detailed structures in the lightcurves, and discuss some possible origins for the detailed variability behaviors.

We investigate the impact of different redshift distributions of random samples on the baryon acoustic oscillations (BAO) measurements of $D_V(z)r_\mathrm{d}^\mathrm{fid}/r_\mathrm{d}$ from the two-point correlation functions (2PCF) of galaxies in the Data Release 12 (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). Big surveys, such as BOSS, usually assign redshifts to the random samples by randomly drawing values from the measured redshift distributions of the data, which would necessarily introduce fiducial signals of fluctuations into the random samples, weakening the signals of BAO, if the cosmic variance cannot be ignored. We propose a smooth function of redshift distribution that fits the data well to populate the random galaxy samples. The resulting cosmological parameters match the input parameters of the mock catalogue very well. The significance of BAO signals has been improved by $0.33\sigma$ for a low-redshift (LOWZ) sample and by $0.03\sigma$ for a constant-stellar-mass (CMASS) sample, though the absolute values do not change significantly. Given the precision of the measurements of current cosmological parameters, it would be appreciated for the future improvements on the measurements of galaxy clustering.

We take advantage of the statistical power of the large-volume dark-matter-only Millennium simulation, combined with a sophisticated semi-analytic galaxy formation model, to explore whether the recently reported $z=3.7$ quiescent galaxy ZF-COSMOS-20115 (ZF; Glazebrook et al. 2017) can be accommodated in current galaxy formation models. In our model, a population of quiescent galaxies (QGs) with stellar masses and star formation rates comparable to those of ZF naturally emerges at redshifts $z<4$. There are two and five ZF analogues at the redshift $3.86$ and $3.58$ in the Millennium simulation volume, respectively. We demonstrate that, while the $z>3.5$ massive QGs are rare (about 2\% of the galaxies with the similar stellar masses), the existing AGN feedback model implemented in the semi-analytic galaxy formation model can successfully explain the formation of the high-redshift QGs as it does on their lower redshift counterparts.

The standard galaxy formation theory assumes that baryons and dark matter are initially well-mixed before becoming segregated due to radiative cooling. We use non-radiative hydrodynamical simulations to explicitly examine this assumption and find that baryons and dark matter can also be segregated because of different physics obeyed by gas and dark matter during the build-up of the halo. As a result, baryons in many haloes do not originate from the same Lagrangian region as the dark matter. When using the fraction of corresponding dark matter and gas particles in the initial conditions (the "paired fraction") as a proxy of the dark matter and gas segregation strength of a halo, on average about $25$ percent of the baryonic and dark matter of the final halo are segregated in the initial conditions. This is at odds with the assumption of the standard galaxy formation model. A consequence of this effect is that the baryons and dark matter of the same halo initially experience different tidal torques and thus their angular momentum vectors are often misaligned. The degree of the misalignment is largely preserved during later halo assembly and can be understood with the tidal torque theory. The result challenges the precision of some semi-analytical approaches which utilize dark matter halo merger trees to infer properties of gas associated to dark matter haloes.

A particular population of galaxies have drawn much interest recently, which are as faint as typical dwarf galaxies but have the sizes as large as $L^*$ galaxies, the so called "ultra-diffuse galaxie" (UDGs). The lack of tidal features of UDGs in dense environments suggests that their host halos are perhaps as massive as that of the Milky Way. On the other hand, galaxy formation efficiency should be much higher in the halos of such masses. Here we use the model galaxy catalog generated by populating two large simulations: the Millennium-II cosmological simulation and Phoenix simulations of 9 big clusters with the semi-analytic galaxy formation model. This model reproduces remarkably well the observed properties of UDGs in the nearby clusters, including the abundance, profile, color, and morphology, etc. We search for UDG candidates using the public data and find 2 UDG candidates in our Local Group and 23 in our Local Volume, in excellent agreement with the model predictions. We demonstrate that UDGs are genuine dwarf galaxies, formed in the halos of $\sim 10^{10}M_{\odot}$. It is the combination of the late formation time and high-spins of the host halos that results in the spatially extended feature of this particular population. The lack of tidal disruption features of UDGs in clusters can also be explained by their late infall-time.

We study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model massive galaxies are quenched by AGN feedback depending on both black hole and hot gas mass, and hence indirectly on stellar mass. In addition, satellite galaxies of any mass can be quenched by ram-pressure or tidal stripping of gas and through the suppression of gaseous infall. This combination of processes produces quenching efficiencies which depend on stellar mass, host halo mass, environment density, distance to group centre and group central galaxy properties in ways which agree qualitatively with observation. Some discrepancies remain in dense regions and close to group centres, where quenching still seems too efficient. In addition, although the mean stellar age of massive galaxies agrees with observation, the assumed AGN feedback model allows too much ongoing star formation at late times. The fact that both AGN feedback and environmental effects are stronger in higher density environments leads to a correlation between the quenching of central and satellite galaxies which roughly reproduces observed conformity trends inside haloes.

We study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

We present a detailed study of the formation of an inverse S-shaped filament prior to its eruption in active region NOAA 11884 from October 31 to November 2, 2013. In the initial stage, clockwise rotation of a small positive sunspot around the main negative trailing sunspot formed a curved filament. Then the small sunspot cancelled with negative magnetic flux to create a longer active-region filament with an inverse S-shape. At the cancellation site a brightening was observed in UV and EUV images and bright material was transferred to the filament. Later the filament erupted after cancellation of two opposite polarities under the upper part of the filament. Nonlinear force-free field (NLFFF) extrapolation of vector photospheric fields suggests that the filament may have a twisted structure, but this cannot be confirmed from the current observations.

The distribution of smaller satellite galaxies around large central galaxies has attracted attention because peculiar spatial and kinematic configurations have been detected in some systems. A particularly striking example of such behavior is seen in the satellite system of the Andromeda galaxy, where around 80\% are on the nearside of that galaxy, facing the Milky Way. Motivated by this departure from anisotropy, we examined the spatial distribution of satellites around pairs of galaxies in the SDSS. By stacking tens of thousands of satellites around galaxy pairs we found that satellites tend to bulge towards the other central galaxy, preferably occupying the space between the pair, rather than being spherically or axis-symmetrically distributed around each host. The bulging is a function of the opening angle examined and is fairly strong -- there are up to $\sim$10\% more satellites in the space between the pair, than expected from uniform. Consequently, it is a statistically very strong signal, being inconsistent with a uniform distribution at the 5$\sigma$ level. The possibility that the observed signal is the result of the overlap of two haloes with extended satellite distributions, is ruled out by testing this hypothesis by performing the same tests on isolated galaxies (and their satellites) artificially placed at similar separations. These findings highlight the unrelaxed and interacting nature of galaxies in pairs.

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

We have updated the Munich galaxy formation model to the Planck first-year cosmology, while modifying the treatment of baryonic processes to reproduce recent data on the abundance and passive fractions of galaxies from z= 3 down to z=0. Matching these more extensive and more precise observational results requires us to delay the reincorporation of wind ejecta, to lower the surface density threshold for turning cold gas into stars, to eliminate ram-pressure stripping in haloes less massive than ~10^14 Msun, and to modify our model for radio mode feedback. These changes cure the most obvious failings of our previous models, namely the overly early formation of low-mass galaxies and the overly large fraction of them that are passive at late times. The new model is calibrated to reproduce the observed evolution both of the stellar mass function and of the distribution of star formation rate at each stellar mass. Massive galaxies (M>10^11 [Msun]) assemble most of their mass before z=1 and are predominantly old and passive at z=0, while lower mass galaxies assemble later and, for M<10^9.5 (Msun), are still predominantly blue and star forming at z=0. This phenomenological but physically based model allows the observations to be interpreted in terms of the efficiency of the various processes that control the formation and evolution of galaxies as a function of their stellar mass, gas content, environment and time.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($\Lambda$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $\Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $\Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $\Lambda$CDM predictions. In fact, ~10% of $\Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ~90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ~10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a 'cD envelope'. These clouds typically originate from tidal stripping at z < 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ~30 per cent.

The accretion of satellites onto central galaxies along vast cosmic filaments is an apparent outcome of the anisotropic collapse of structure in our Universe. Numerical work (based on gravitational dynamics of N-body simulations) indicates that satellites are beamed towards hosts along preferred directions imprinted by the velocity shear field. Here we use the Sloan Digital Sky Survey to observationally test this claim. We construct 3D filaments and sheets and examine the relative position of satellite galaxies. A statistically significant alignment between satellite galaxy position and filament axis in observations is confirmed. We find a qualitatively compatible alignments by examining satellites and filaments similarly identified in the Millennium simulation, semi-analytical galaxy catalogue. We also examine the dependence of the alignment strength on galaxy properties such as colour, magnitude and (relative) satellite magnitude, finding that the alignment is strongest for the reddest and brightest central and satellite galaxies. Our results confirm the theoretical picture and the role of the cosmic web in satellite accretion. Furthermore our results suggest that filaments identified on larger scales can be reflected in the positions of satellite galaxies that are quite close to their hosts.

We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the "too big to fail" problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.

We investigate whether the satellite luminosity function (LF) of primary galaxies identified in the Sloan Digital Sky Survey (SDSS) depends on whether the host galaxy is in a filament or not. Isolated primary galaxies are identified in the SDSS spectroscopic sample while potential satellites (that are up to 4 magnitudes fainter than their hosts) are searched for in the much deeper photometric sample. Filaments are constructed from the galaxy distribution by the "Bisous" process. Isolated primary galaxies are divided into two subsamples: those in filaments and those not in filaments. We examine the stacked mean satellite LF of both the filament and non-filament sample and find that, on average, the satellite LFs of galaxies in filaments is significantly higher than those of galaxies not in filaments. The filamentary environment can increases the abundance of the brightest satellites ($M_\mathrm{sat.} < M_\mathrm{prim.} + 2.0$), by a factor of $\sim 2$ compared with non-filament isolated galaxies. This result is independent of primary galaxy magnitude although the satellite LF of galaxies in the faintest magnitude bin, is too noisy to determine if such a dependence exists. Since our filaments are extracted from a spectroscopic flux-limited sample, we consider the possibility that the difference in satellite LF is due to a redshift, colour or environmental bias, finding these to be insufficient to explain our result. The dependence of the satellite LF on the cosmic web suggests that the filamentary environment may have a strong effect on the efficiency of galaxy formation.

Cadmium Zinc Telluride and Cadmium Telluride are the detector materials of choice for the detection of X-rays in the X-ray energy band E >= 5keV with excellent spatial and spectral resolution and without cryogenic cooling. Owing to recent breakthroughs in grazing incidence mirror technology, next-generation hard X-ray telescopes will achieve angular resolution between 5 and 10 arc seconds - about an order of magnitude better than that of the NuSTAR hard X-ray telescope. As a consequence, the next generation of X-ray telescopes will require pixelated X-ray detectors with pixels on a grid with a lattice constant of <= 250um. Additional detector requirements include a low energy threshold of less than 5keV and an energy resolution of less than one keV. The science drivers for a high angular-resolution X-ray mission include studies and measurements of black hole spins, the cosmic evolution of super-massive black holes, active galactic nuclei feedback, and the behaviour of matter at very high densities. In this contribution, we report on our R&D studies with the goal to optimise small-pixel Cadmium Zinc Telluride and Cadmium Telluride detectors.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, neutron stars, and gamma-ray bursts. We designed, built and tested a X-ray polarimeter, X-Calibur, to be used in the focal plane of the balloon-borne InFOCuS grazing incidence X-ray telescope. X-Calibur combines a low-Z scatterer with a CZT detector assembly to measure the polarization of 20-80keV X-rays making use of the fact that polarized photons scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of ~80%. The X-Calibur detector assembly is completed, tested, and fully calibrated. The response to a polarized X-ray beam was measured successfully at the Cornell High Energy Synchrotron Source. This paper describes the design, calibration and performance of the X-Calibur polarimeter. In principle, a similar space-borne scattering polarimeter could operate over the broader 2-100keV energy band.

X-ray polarimetry has great scientific potential and new experiments, such as X-Calibur, PoGOLite, XIPE, and GEMS, will not only be orders of magnitude more sensitive than previous missions, but also provide the capability to measure polarization over a wide energy range. However, the measured spectra depend on the collection area, detector responses, and, in case of balloon-borne experiments, the absorption of X-rays in the atmosphere, all of which are energy dependent. Combined with the typically steep source spectra, this leads to significant biases that need to be taken into account to correctly reconstruct energy-resolved polarization properties. In this paper, we present a method based on an iterative unfolding algorithm that makes it possible to simultaneously reconstruct the energy spectrum and the polarization properties as a function of true photon energy. We apply the method to a simulated X-Calibur data set and show that it is able to recover both the energy spectrum and the energy-dependent polarization fraction.

Recent work has suggested that the amplitude of the size mass relation of massive early type galaxies evolves with redshift. Here we use a semi-analytical galaxy formation model to study the size evolution of massive early type galaxies. We find this model is able to reproduce the amplitude of present day amplitude and slope of the relation between size and stellar mass for these galaxies, as well as its evolution. The amplitude of this relation reflects the typical compactness of dark halos at the time when most of the stars are formed. This link between size and star formation epoch is propagated in galaxy mergers. Mergers of high or moderate mass ratio (less than 1:3) become increasingly important with increasing present day stellar mass for galaxies more massive than $10^{11.4}M_{\odot}$. At lower masses, low mass ratio mergers play a more important role. In situ star formation contribute more to the size growth than it does to stellar mass growth. We also find that, for ETGs identified at $z=2$, minor mergers dominate subsequent growth both for stellar mass and in size, consistent with earlier theoretical results.

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

The emptiness of the Local Void has been put forward as a serious challenge to the current standard paradigm of structure formation in LCDM. We use a high resolution cosmological N-body simulation, the Millennium-II run, combined with a sophisticated semi-analytic galaxy formation model, to explore statistically whether the local void is allowed within our current knowledge of galaxy formation in LCDM. We find that about 15 percent of the Local Group analogue systems (11 of 77) in our simulation are associated with nearby low density regions having size and 'emptiness' similar to those of the observed Local Void. This suggests that, rather than a crisis of the LCDM, the emptiness of the Local Void is indeed a success of the standard LCDM theory. The paucity of faint galaxies in such voids results from a combination of two factors: a lower amplitude of the halo mass function in the voids than in the field, and a lower galaxy formation efficiency in void haloes due to halo assembly bias effects. While the former is the dominated factor, the later also plays a sizable role. The halo assembly bias effect results in a stellar mass fraction 25 percent lower for void galaxies when compared to field galaxies with the same halo mass.

Subhalo abundance matching (SHAM) inserts galaxies into dark matter only simulations of the growth of cosmic structure in a way that requires minimal assumptions about galaxy formation. A galaxy is placed at the potential minimum of each distinct self-bound subhalo with a luminosity which is a monotonically increasing function of the maximum mass (or circular velocity) attained over the subhalo's earlier history. Galaxy and subhalo properties are linked by matching model and observed luminosity functions. Simulated structures can then be compared in detail with observation, for example, through galaxy correlation statistics, group catalogues, or galaxy-galaxy lensing. Robust astrophysical conclusions can be drawn from such a comparison only on scales which are unaffected by the numerical limitations of the simulation. Here we compare results for the Millennium Simulation (MS) with subhalos defined using the SUBFIND to those obtained applying identical analysis to the much higher resolution Millennium-II. Correlation statistics on scales between 200 kpc and 2 Mpc converge to within 20% only for subhalos with masses at infall corresponding to at least 1000 simulation particles in the MS. Numerically converged results can be obtained to much lower infall particle number, if galaxies are followed even after their associated subhalos have been tidally disrupted, as in most recent semi-analytic galaxy formation simulations. This allows robust comparison between simulation and observation over a wider dynamic range in mass than for a SHAM analysis which ignores such "orphan" galaxies.

Previous studies indicate that assembly bias effects are stronger for lower mass dark matter haloes. Here we make use of high resolution re-simulations of rich clusters and their surroundings from the Phoenix Project and a large volume cosmological simulation, the Millennium-II run, to quantify assembly bias effects on dwarf-sized dark matter haloes. We find that, in the regions around massive clusters, dwarf-sized haloes ($[10^9,10^{11}]\ms$) form earlier ($\Delta z \sim 2$ in redshift) and possess larger $V_{\rm max}$ ($\sim20%$) than the field galaxies. We find that this environmental dependence is largely caused by tidal interactions between the ejected haloes and their former hosts, while other large scale effects are less important. Finally we assess the effects of assembly bias on dwarf galaxy formation with a sophisticated semi-analytical galaxy formation model. We find that the dwarf galaxies near massive clusters tend to be redder ($\Delta(u-r) = 0.5$) and have three times as much stellar mass compared to the field galaxies with the same halo mass. These features should be seen with observational data.

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

We have updated our radially-resolved SAMs of galaxy formation, which track both the atomic and molecular gas phases of the ISM. The models are adapted from those of Guo et al. using similar methodology as in Fu et al. and are run on halo merger trees from the MS and MS II with the following main changes: (1) We adopt a simple star formation law where \Sigma_SFR \propto \Sigma_H2 (2) We inject the heavy elements produced by supernovae directly into the halo hot gas, instead of first mixing them with the cold gas in the disk. (3) We include radial gas inflows in disks using a model of the form v_inflow = \alpha r. The models are used to study the radial profiles of star formation rate and gas-phase metallicity in present-day galaxies. The \Sigma_H2 profiles in L* galaxies place strong constraints on inflow velocities, favouring models where v_inflow~7km/s at a galactocentric radius of 10kpc. Radial gas inflow has little influence on gas-phase and stellar metallicity gradients, which are affected much more strongly by the fraction of metals that are directly injected into the halo gas, rather than mixed with the cold gas. Metals ejected out of the galaxy in early epochs result in late infall of pre-enriched gas and flatter present-day gas-phase metallicity gradients. A prescription in which 80% of the metals are injected into the halo gas results in good fits to the flat observed metallicity gradients in galaxies with stellar masses greater than 10^10 M_sun, as well as the relations between gas-phase metallicity and sSFR in the outer parts of galactic disks. We examine the correlation between gas-phase metallicity gradient and global galaxy properties, finding that it is most strongly correlated with the bulge-to-total ratio of the galaxy. This is because gas is consumed when the bulge forms during galaxy mergers, and the gas-phase metallicity gradient is then set by newly-accreted gas.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

Using the technique of Angulo & White (2010) we scale the Millennium and Millennium-II simulations of structure growth in a LCDM universe from the cosmological parameters with which they were carried out (based on first-year results from the Wilkinson Microwave Anisotropy Probe, WMAP1) to parameters consistent with the seven-year WMAP data (WMAP7). We implement semi-analytic galaxy formation modelling on both simulations in both cosmologies to investigate how the formation, evolution and clustering of galaxies are predicted to vary with cosmological parameters. The increased matter density Omega_m and decreased linear fluctuation amplitude sigma8 in WMAP7 have compensating effects, so that the abundance and clustering of dark halos are predicted to be very similar to those in WMAP1 for z <= 3. As a result, local galaxy properties can be reproduced equally well in the two cosmologies by slightly altering galaxy formation parameters. The evolution of the galaxy populations is then also similar. In WMAP7, structure forms slightly later. This shifts the peak in cosmic star formation rate to lower redshift, resulting in slightly bluer galaxies at z=0. Nevertheless, the model still predicts more passive low-mass galaxies than are observed. For rp< 1Mpc, the z=0 clustering of low-mass galaxies is weaker for WMAP7 than for WMAP1 and closer to that observed, but the two cosmologies give very similar results for more massive galaxies and on large scales. At z>1 galaxies are predicted to be more strongly clustered for WMAP7. Differences in galaxy properties, including, clustering, in these two cosmologies are rather small up to redshift 3. Given that there are still considerable residual uncertainties in galaxy formation models, it is very difficult to distinguish WMAP1 from WMAP7 through observations of galaxy properties or their evolution.

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed.

In a semi-numerical model of reionization, the evolution of ionization fraction is simulated approximately by the ionizing photon to baryon ratio criterion. In this paper we incorporate a semi-analytical model of galaxy formation based on the Millennium II N-body simulation into the semi-numerical modeling of reionization. The semi-analytical model is used to predict the production of ionizing photons, then we use the semi-numerical method to model the reionization process. Such an approach allows more detailed modeling of the reionization, and also connects observations of galaxies at low and high redshifts to the reionization history. The galaxy formation model we use was designed to match the low-$z$ observations, and it also fits the high redshift luminosity function reasonably well, but its prediction on the star formation falls below the observed value, and we find that it also underpredicts the stellar ionizing photon production rate, hence the reionization can not be completed at $z \sim 6$ without taking into account some other potential sources of ionization photons. We also considered simple modifications of the model with more top heavy initial mass functions (IMF), with which the reionization can occur at earlier epochs. The incorporation of the semi-analytical model may also affect the topology of the HI regions during the EoR, and the neutral regions produced by our simulations with the semi-analytical model appeared less poriferous than the simple halo-based models.

We report on the optimization of the hard X-ray polarimeter X-Calibur for a high-altitude balloon-flight in the focal plane of the InFOC{\mu}S X-ray telescope from Fort Sumner (NM) in Fall 2013. X-Calibur combines a low-Z scintillator slab to Compton-scatter photons with a high-Z Cadmium Zinc Telluride (CZT) detector assembly to photo-absorb the scattered photons. The detector makes use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity and reaches a sensitivity close to the best theoretically possible. In this paper, we discuss the optimization of the design of the instrument based on Monte Carlo simulations of polarized and unpolarized X-ray beams and of the most important background components. We calculate the sensitivity of the polarimeter for the upcoming balloon flight from Fort Sumner and for additional longer balloon flights with higher throughput mirrors. We conclude by emphasizing that Compton polarimeters on satellite borne missions can be used down to energies of a few keV.

[abridged] We construct lightcones for the semi-analytic galaxy formation simulation of Guo et al. (2011) and make mock catalogues for comparison with deep high-redshift surveys. Photometric properties are calculated with two different stellar population synthesis codes (Bruzual & Charlot 2003; Maraston 2005) in order to study sensitivity to this aspect of the modelling. The catalogues are publicly available and include photometry for a large number of observed bands from 4000{\deg}A to 6{\mu}m, as well as rest-frame photometry and intrinsic properties of the galaxies. Guo et al. (2011) tuned their model to fit the low-redshift galaxy population but noted that at z > 1 it overpredicts the abundance of galaxies below the "knee" of the stellar mass function. Here we extend the comparison to deep galaxy counts in the B, i, J, K and IRAC 3.6{\mu}m, 4.5{\mu}m and 5.8{\mu}m bands, to the redshift distributions of K and 5.8{\mu}m selected galaxies, and to the evolution of rest-frame luminosity functions in the B and K bands. The B, i and J counts are well reproduced, but at longer wavelengths the overabundant high-redshift galaxies produce excess faint counts. The predicted redshift distributions for K and 5.8{\mu}m selected samples highlight the effect of emission from thermally pulsing AGB stars. The full treatment of Maraston (2005) predicts three times as many z~2 galaxies in faint 5.8{\mu}m selected samples as the model of Bruzual & Charlot (2003), whereas the two models give similar predictions for K-band selected samples. Although luminosity functions are adequately reproduced out to z~3 in rest-frame B, the same is true at rest-frame K only if TP-AGB emission is included, and then only at high luminosity. Fainter than L* the two synthesis models agree but overpredict the number of galaxies, another reflection of the overabundance of ~10^10M\odot model galaxies at z > 1.

We apply Monte Carlo Markov Chain (MCMC) methods to large-scale simulations of galaxy formation in a LambdaCDM cosmology in order to explore how star formation and feedback are constrained by the observed luminosity and stellar mass functions of galaxies. We build models jointly on the Millennium and Millennium-II simulations, applying fast sampling techniques which allow observed galaxy abundances over the ranges 7<log(M*/Msun)<12 and z=0 to z=3 to be used simultaneously as constraints in the MCMC analysis. When z=0 constraints alone are imposed, we reproduce the results of previous modelling by Guo et al. (2012), but no single set of parameters can reproduce observed galaxy abundances at all redshifts simultaneously, reflecting the fact that low-mass galaxies form too early and thus are overabundant at high redshift in this model. The data require the efficiency with which galactic wind ejecta are reaccreted to vary with redshift and halo mass quite differently than previously assumed, but in a similar way as in some recent hydrodynamic simulations of galaxy formation. We propose a specific model in which reincorporation timescales vary inversely with halo mass and are independent of redshift. This produces an evolving galaxy population which fits observed abundances as a function of stellar mass, B- and K-band luminosity at all redshifts simultaneously. It also produces a significant improvement in two other areas where previous models were deficient. It leads to present day dwarf galaxy populations which are younger, bluer, more strongly star-forming and more weakly clustered on small scales than before, although the passive fraction of faint dwarfs remains too high.

We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ~50, 000 galaxy systems. We find no evidence for any anistropy in the satellite galaxy distribution relative to the major axes of the primaries. We derive accurate projected number density profiles of satellites down to 4 magnitudes fainter than their primaries. We find the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and on the colour and concentration of their primaries.

We study two nearby, early-type galaxies, NGC4342 and NGC4291, that host unusually massive black holes relative to their low stellar mass. The observed black hole-to-bulge mass ratios of NGC4342 and NGC4291 are ~6.9% and ~1.9%, respectively, which significantly exceed the typical observed ratio of ~0.2%. As a consequence of the exceedingly large black hole-to-bulge mass ratios, NGC4342 and NGC4291 are ~5.1 sigma and ~3.4 sigma outliers from the M_BH - M_bulge scaling relation, respectively. In this paper, we explore the origin of the unusually high black hole-to-bulge mass ratio. Based on Chandra X-ray observations of the hot gas content of NGC4342 and NGC4291, we compute gravitating mass profiles, and conclude that both galaxies reside in massive dark matter halos, which extend well beyond the stellar light. The presence of dark matter halos around NGC4342 and NGC4291 and a deep optical image of the environment of NGC4342 indicate that tidal stripping, in which >90% of the stellar mass was lost, cannot explain the observed high black hole-to-bulge mass ratios. Therefore, we conclude that these galaxies formed with low stellar masses, implying that the bulge and black hole did not grow in tandem. We also find that the black hole mass correlates well with the properties of the dark matter halo, suggesting that dark matter halos may play a major role in regulating the growth of the supermassive black holes.

We examine how the atomic and molecular gas components of galaxies evolve to higher redshifts using the semi-analytic galaxy formation models of Fu et al. (2010) in which we track the surface density profiles of gas in disks. We adopt two different prescriptions based either on gas surface density and metallicity, or on interstellar pressure, to compute the molecular fraction as a function of radius in each disk. We demonstrate that the adopted star formation law determines how the {\em balance} between gas, stars and metals changes with time in the star-forming galaxy population, but does not influence the total mass in stars formed into galaxies at redshifts below $z\sim 2.5$. The redshift evolution of the mass-metallicity relation places strong constraints on the timescale over which cold gas is converted into stars in high redshift galaxies, and favours models where this remains constant at values around 1-2 Gyr. Future observations of the evolution of the average molecular-to-atomic gas ratio in galaxies as a function of stellar mass and redshift will constrain models of the atomic-to-molecular transition.

The 6 cm and 20 cm radio continuum properties of all 85 galaxies with reported 22 GHz H2O maser emission and luminosity distance D > 0.5 Mpc are studied. For the total of 55 targets for which both 6 cm and 20 cm measurements exist and for the subsample of 42 sources with masers related to active galactic nuclei (AGN), a spectral index could be determined from an assumed power-law dependence. The mean value of the resulting spectral index is in both cases 0.66+-0.07. Comparing radio properties of the maser galaxies with a sample of Seyferts without detected H2O maser, we find that (1) the spectral indices agree within the error limits, and (2) maser host galaxies have higher nuclear radio continuum luminosities, exceeding those of the comparison sample by factors of order 5. Only considering the subsample of galaxies with masers associated with AGN, there seems to be a trend toward rising maser luminosity with nuclear radio luminosity (both at 6 cm and 20 cm). However, when accounting for the Malmquist effect, the correlation weakens to a level, which is barely significant. Overall, the study indicates that nuclear radio luminosity is a suitable indicator to guide future AGN maser searches and to enhance detection rates, which are otherwise quite low (<10%).

X-ray polarimetry promises to give new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

We present measurements of projected autocorrelation functions w_p(r_p) for the stellar mass of galaxies and for their light in the U, B and V bands, using data from the third data release of the DEEP2 Galaxy Redshift Survey and the final data release of the Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass and light by comparing these to projected autocorrelations of dark matter estimated from the Millennium Simulations (MS) at z=1 and 0.07, the median redshifts of our galaxy samples. All of the autocorrelation and bias functions show systematic trends with spatial scale and waveband which are impressively similar at the two redshifts. This shows that the well-established environmental dependence of stellar populations in the local Universe is already in place at z=1. The recent MS-based galaxy formation simulation of Guo et al. (2011) reproduces the scale-dependent clustering of luminosity to an accuracy better than 30% in all bands and at both redshifts, but substantially overpredicts mass autocorrelations at separations below about 2 Mpc. Further comparison of the shapes of our stellar mass bias functions with those predicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation amplitude of sigma_8 0.8 rather than the sigma_8=0.9 assumed by the MS.

We study relations between stellar mass, star formation and gas-phase metallicity in a sample of 177,071 unique emission line galaxies from the SDSS-DR7, as well as in a sample of 43,767 star forming galaxies at z=0 from the cosmological semi-analytic model L-GALAXIES. We demonstrate that metallicity is dependent on star formation rate at fixed mass, but that the trend is opposite for low and for high mass galaxies. Low-mass galaxies that are actively forming stars are more metal-poor than quiescent low-mass galaxies. High-mass galaxies, on the other hand, have lower gas-phase metallicities if their star formation rates are small. Remarkably, the same trends are found for our sample of model galaxies. We find that massive model galaxies with low gas-phase metallicities have undergone a gas-rich merger in the past, inducing a starburst which exhausted their cold gas reservoirs and shut down star formation. This led to a gradual dilution in the gas-phase metallicities of these systems via accretion of gas. These model galaxies have lower-than-average gas-to-stellar mass ratios and higher-than-average central black hole masses. We confirm that massive galaxies with low gas-phase metallicities in our observational sample also have very massive black holes. We propose that accretion may therefore play a significant role in regulating the gas-phase metallicities of present-day massive galaxies.

We study the luminosity function of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For primaries of similar luminosity to the Milky Way and M31, we are able to stack as many as $\sim\negthinspace 20,000$ galaxy systems to obtain robust statistical results. We derive the satellite luminosity function extending almost 8 magnitudes fainter than the primary galaxy. We also determine how the satellite luminosity function varies with the luminosity, colour and concentration of the primary. We find that, in the mean, isolated primaries of comparable luminosity to the Milky Way and M31 contain about a factor of two fewer satellites brighter than $M_V=-14$ than the average of the Milky Way and M31.

We compare the galaxy population in the Virgo, Fornax, Coma and Perseus cluster to a state-of-the-art semi-analytic model, focusing on the regime of dwarf galaxies with luminosities from approximately 10^8 L_sun to 10^9 L_sun. We find that the number density profiles of dwarfs in observed clusters are reproduced reasonably well, and that the red fractions of model clusters provide a good match to Coma and Perseus. On the other hand, the red fraction among dwarf galaxies in Virgo is clearly lower than in model clusters. We argue that this is mainly caused by the treatment of environmental effects in the model. This explanation is supported by our finding that the colours of central ("field") dwarf galaxies are reproduced well, in contrast to previous claims. Finally, we find that the dwarf-to-giant ratio in model clusters is too high. This may indicate that the current model prescription for tidal disruption of faint galaxies is still not efficient enough.

We study the relation between molecular gas and star formation in a volume-limited sample of 222 galaxies from the COLD GASS survey, with measurements of the CO(1-0) line from the IRAM 30m telescope. The galaxies are at redshifts 0.025<z<0.05 and have stellar masses in the range 10.0<log(M*/Msun)<11.5. The IRAM measurements are complemented by deep Arecibo HI observations and homogeneous SDSS and GALEX photometry. A reference sample that includes both UV and far-IR data is used to calibrate our estimates of star formation rates from the seven optical/UV bands. The mean molecular gas depletion timescale, tdep(H2), for all the galaxies in our sample is 1 Gyr, however tdep(H2) increases by a factor of 6 from a value of ~0.5 Gyr for galaxies with stellar masses of 10^10 Msun to ~3 Gyr for galaxies with masses of a few times 10^11 Msun. In contrast, the atomic gas depletion timescale remains contant at a value of around 3 Gyr. This implies that in high mass galaxies, molecular and atomic gas depletion timescales are comparable, but in low mass galaxies, molecular gas is being consumed much more quickly than atomic gas. The strongest dependences of tdep(H2) are on the stellar mass of the galaxy (parameterized as log tdep(H2)= (0.36+/-0.07)(log M* - 10.70)+(9.03+/-0.99)), and on the specific star formation rate. A single tdep(H2) versus sSFR relation is able to fit both "normal" star-forming galaxies in our COLD GASS sample, as well as more extreme starburst galaxies (LIRGs and ULIRGs), which have tdep(H2) < 10^8 yr. Normal galaxies at z=1-2 are displaced with respect to the local galaxy population in the tdep(H2) versus sSFR plane and have molecular gas depletion times that are a factor of 3-5 times longer at a given value of sSFR due to their significantly larger gas fractions.

We are conducting COLD GASS, a legacy survey for molecular gas in nearby galaxies. Using the IRAM 30m telescope, we measure the CO(1-0) line in a sample of ~350 nearby (D=100-200 Mpc), massive galaxies (log(M*/Msun)>10.0). The sample is selected purely according to stellar mass, and therefore provides an unbiased view of molecular gas in these systems. By combining the IRAM data with SDSS photometry and spectroscopy, GALEX imaging and high-quality Arecibo HI data, we investigate the partition of condensed baryons between stars, atomic gas and molecular gas in 0.1-10L* galaxies. In this paper, we present CO luminosities and molecular hydrogen masses for the first 222 galaxies. The overall CO detection rate is 54%, but our survey also uncovers the existence of sharp thresholds in galaxy structural parameters such as stellar mass surface density and concentration index, below which all galaxies have a measurable cold gas component but above which the detection rate of the CO line drops suddenly. The mean molecular gas fraction MH2/M* of the CO detections is 0.066+/-0.039, and this fraction does not depend on stellar mass, but is a strong function of NUV-r colour. Through stacking, we set a firm upper limit of MH2/M*=0.0016+/-0.0005 for red galaxies with NUV-r>5.0. The average molecular-to-atomic hydrogen ratio in present-day galaxies is 0.3, with significant scatter from one galaxy to the next. The existence of strong detection thresholds in both the HI and CO lines suggests that "quenching" processes have occurred in these systems. Intriguingly, atomic gas strongly dominates in the minority of galaxies with significant cold gas that lie above these thresholds. This suggests that some re-accretion of gas may still be possible following the quenching event.

CdZnTe (CZT) is a wide bandgap II-VI semiconductor developed for the spectroscopic detection of X-rays and {\gamma}-rays at room temperature. The Swift Burst Alert Telescope is using an 5240 cm2 array of 2 mm thick CZT detectors for the detection of 15-150 keV X-rays from Gamma-Ray Bursts. We report on the systematic tests of thicker (\geq 0.5 cm) CZT detectors with volumes between 2 cm3 and 4 cm3 which are potential detector choices for a number of future X-ray telescopes that operate in the 10 keV to a few MeV energy range. The detectors contacted in our laboratory achieve Full Width Half Maximum energy resolutions of 2.7 keV (4.5%) at 59 keV, 3 keV (2.5%) at 122 keV and 4 keV (0.6%) at 662 keV. The 59 keV and 122 keV energy resolutions are among the world-best results for \geq 0.5 cm thick CZT detectors. We use the data set to study trends of how the energy resolution depends on the detector thickness and on the pixel pitch. Unfortunately, we do not find clear trends, indicating that even for the extremely good energy resolutions reported here, the achievable energy resolutions are largely determined by the properties of individual crystals. Somewhat surprisingly, we achieve the reported results without applying a correction of the anode signals for the depth of the interaction. Measuring the interaction depths thus does not seem to be a pre-requisite for achieving sub-1% energy resolutions at 662 keV.

We apply updated semi-analytic galaxy formation models simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II simulations. These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatments of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges and of gaseous and stellar disks, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over 5 orders of magnitude in stellar mass and 9 magnitudes in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionisation effects may not be needed to solve the "missing satellite" problem except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ~30kpc for massive galaxies. For M* < 6 x 10^10Msun, however, the model overpredicts clustering at scales below 1 Mpc, suggesting that the sigma_8 adopted in the simulations (0.9) is too high. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhalos (so-called orphans) are included. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ~10^10Msun galaxies beyond z~0.6, reflecting deficiencies in the way star-formation rates are modelled.

We report on Monte Carlo studies of the hard X-ray polarimeter X-Calibur. The polarimeter will be used in the focal plane of a grazing incidence hard X-ray telescope. It combines a low-Z Compton scatterer with a high-Z Cadmium Zinc Telluride (CZT) detector assembly to measure the polarization of 10 keV - 80 keV X-rays. X-Calibur makes use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. In contrast of competing designs, which use only a small fraction of the incoming X-rays, X-Calibur achieves a high detection efficiency of order unity. In this contributions, we discuss a Monte Carlo study which compares X-Calibur's polarimeteric performance achieved using different scattering materials (Scintillator, Be, LiH, Li), and calculate the sensitivity of X-Calibur when used with different balloon-borne and space-borne mirror assemblies.

We present cosmological hydrodynamical simulations of the formation of dwarf galaxies in a representative sample of haloes extracted from the Millennium-II Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show different mass assembly histories which are reflected in different star formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8 solar masses, consistent with other published simulations of galaxy formation in similar mass haloes. Our final objects have structures and stellar populations consistent with dwarf elliptical and dwarf irregular galaxies. However, in a Lambda CDM universe, 10^10 solar mass haloes must typically contain galaxies with much lower stellar mass than our simulated objects if they are to match observed galaxy abundances. The dwarf galaxies formed in our own and all other current hydrodynamical simulations are more than an order of magnitude more luminous than expected for haloes of this mass. We discuss the significance and possible implications of this result.

X-ray polarimetry promises to give qualitatively new information about high-energy sources. Examples of interesting source classes are binary black hole systems, rotation and accretion powered neutron stars, Microquasars, Active Galactic Nuclei and Gamma-Ray Bursts. Furthermore, X-ray polarimetry affords the possibility for testing fundamental physics, e.g. to observe signatures of light bending in the strong gravitational field of a black hole, to detect third order Quantum Electrodynamic effects in the magnetosphere of Magnetars, and to perform sensitive tests of Lorentz Invariance. In this paper we discuss scientific drivers of hard (>10 keV) X-ray polarimetry emphasizing how observations in the hard band can complement observations at lower energies (0.1 - 10 keV). Subsequently, we describe four different technical realizations of hard X-ray polarimeters suitable for small to medium sized space borne missions, and study their performance in the signal-dominated case based on Monte Carlo simulations. We end with confronting the instrument requirements for accomplishing the science goals with the capabilities of the four polarimeters.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We extend existing semi-analytic models of galaxy formation to track atomic and molecular gas in disk galaxies. Simple recipes for processes such as cooling, star formation, supernova feedback, and chemical enrichment of the stars and gas are grafted on to dark matter halo merger trees derived from the Millennium Simulation. Each galactic disk is represented by a series of concentric rings. We assume that surface density profile of infalling gas in a dark matter halo is exponential, with scale radius r_d that is proportional to the virial radius of the halo times its spin parameter $\lambda$. As the dark matter haloes grow through mergers and accretion, disk galaxies assemble from the inside out. We include two simple prescriptions for molecular gas formation processes in our models: one is based on the analytic calculations by Krumholz, McKee & Tumlinson (2008), and the other is a prescription where the H_2 fraction is determined by the kinematic pressure of the ISM. Motivated by the observational results of Leroy et al. (2008), we adopt a star formation law in which $\Sigma_{SFR}\propto\Sigma_{H_2}$ in the regime where the molecular gas dominates the total gas surface density, and $\Sigma_{SFR}\propto \Sigma_{gas}^2$ where atomic hydrogen dominates. We then fit these models to the radial surface density profiles of stars, HI and H_2 drawn from recent high resolution surveys of stars and gas in nearby galaxies. We explore how the ratios of atomic gas, molecular gas and stellar mass vary as a function of global galaxy scale parameters, including stellar mass, stellar surface density, and gas surface density. We elucidate how the trends can be understood in terms of three variables that determine the partition of baryons in disks: the mass of the dark matter halo, the spin parameter of the halo, and the amount of gas recently accreted from the external environment.

We introduced a semi-analytic model to study the detail of radial related processes on galaxy disks. Based on this recipe, we can calculate neutral and molecular gas in interstellar medium, and include the molecular phase into star formation recipes. From our model, we reproduced the radial surface density profiles and mass functions shown in the observational results. We can also give some predictions and explanations to the scale relations of neutral and molecular gas in disk galaxies.

For any assumed stellar Initial Mass Function, the Sloan Digital Sky Survey (SDSS) gives a precise determination of the stellar mass function of galaxies for 10^8 M_sun < M_* < 10^12 M_sun. Within the concordance LCDM cosmology, the Millennium simulations give a precise halo mass function for all halos within which galaxies can form. Under the plausible hypothesis that the stellar mass of a galaxy is an increasing function of the maximum mass ever attained by its halo, these combine to give halo mass as a function of stellar mass. The result agrees quite well with observational estimates of mean halo mass as a function of stellar mass from stacking analyses of the gravitational lensing signal and the satellite dynamics of SDSS galaxies. For M_* ~ 5.5 x 10^10 M_sun, the stellar mass usually assumed for the Milky Way, the implied halo mass is ~ 2 x 10^12 M_sun, consistent with most recent direct estimates and inferences from the MW/M31 Timing Argument. The fraction of the baryons associated with each halo which are present as stars in its central galaxy reaches a maximum of 20% at masses somewhat below that of the Milky Way, and falls rapidly at both higher and lower masses. These conversion efficiencies are lower than in almost all recent high-resolution simulations of galaxy formation, showing that these are not yet viable models for the formation of typical members of the galaxy population. When inserted in the Millennium-II Simulation, our derived relation between stellar mass and halo mass predicts a stellar mass autocorrelation function in excellent agreement with that measured directly in the SDSS. The implied Tully-Fisher relation also appears consistent with observation, suggesting that galaxy luminosity functions and Tully-Fisher relations can be reproduced simultaneously in a LCDM cosmology.

To shed light onto the circumnuclear environment of 22 GHz H2O maser galaxies, we have analyzed some of their multi-wavelength properties, including the far infrared luminosity (FIR), the luminosity of the [O III]\lambda5007 emission line, the nuclear X-ray luminosity, and the equivalent width of the neutral iron Ka emission line (EW (Ka)). Our statistical analysis includes a total of 85 sources, most of them harboring an active galactic nucleus (AGN). There are strong anti-correlations between EW (Ka) and two "optical thickness parameters", i.e. the ratios of the X-ray luminosity versus the presumably more isotropically radiated [O III] and far infrared (FIR) luminosities. Based on these anti-correlations, a set of quantitative criteria, EW (Ka) > 300eV, L_{2-10keV} < 2L_[O III] and L_{FIR} > 600L_{2-10keV} can be established for Compton-thick nuclear regions. 18 H2O maser galaxies belong to this category. There are no obvious correlations between the EW (Ka), the [O III] luminosity and the isotropic H2O maser luminosity. When comparing samples of Seyfert 2s with and without detected H2O maser lines, there seem to exist differences in EW (Ka) and the fraction of Compton-thick nuclei. This should be studied further. For AGN masers alone, there is no obvious correlation between FIR and H2O maser luminosities. However, including masers associated with star forming regions, a linear correlation is revealed. Overall, the extragalactic FIR-H2O data agree with the corresponding relation for Galactic maser sources, extrapolated by several orders of magnitude to higher luminosities.

We investigate large-angle scale temperature anisotropy in the Cosmic Microwave Background (CMB) with the Wilkinson Microwave Anisotropy Probe (WMAP) data and model the large-angle anomalies as the effect of the CMB quadrupole anisotropies caused by the local density inhomogeneities. The quadrupole caused by the local density inhomogeneities is different from the special relativity kinematic quadrupole. If the observer inhabits a strong inhomogeneous region, the local quadrupole should not be neglected. We calculate such local quadrupole under the assumption that there is a huge density fluctuation field in direction $(284^{\circ},74^{\circ})$, where the density fluctuation is $10^{-3}$, and its center is $\sim 112h^{-1} \rm {Mpc}$ away from us. After removing such mock signals from WMAP data, the power in quadrupole, $C_2$, increases from the range $(200\sim260\mu \rm{K^2})$ to $\sim1000\mu \rm{K^2}$. The quantity S, which is used to estimate the alignment between the quadrupole and the octopole, decreases from $(0.7\sim0.74)$ to $(0.31\sim0.37)$, while the model predict that $C_2=1071.5\mu \rm{K^2}$, $S=0.412$. So our local density inhomogeneity model can, in part, explain the WMAP low-$\ell$ anomalies.

We use a sample of 800 galaxies with H I mass measurements from the HyperLeda catalogue and optical photometry from the fourth data release of the Sloan Digital Sky Survey to calibrate a new photometric estimator of the H I to-stellar mass ratio for nearby galaxies. Our estimator, which is motivated by the Kennicutt-Schmidt star formation law, is log(G_{HI}/S) = -1.73238(g-r) + 0.215182mu_i - 4.08451, where mu_i is the i-band surface brighteness and g-r is the optical colour estimated from the g- and r-band Petrosian apparent agnitudes. This estimator has a scatter of sigma = 0.31 dex in log(G_{HI}/S), compared to sigma ~ 0.4 dex for previous estimators that were based on colour alone. We investigate whether the residuals in our estimate of log(G_{HI}/S) depend in a systematic way on a variety of different galaxy properties. We find no effect as a function of stellar mass or 4000A break strength, but there is a systematic effect as a function of the concentration index of the light. We then apply our estimator to a sample of 10^5 emission-line galaxies in the SDSS DR4 and derive an estimate of the H I mass function, which is in excellent agreement with recent results from H I blind surveys. Finally, we re-examine the well-known relation between gas-phase metallicity and stellar mass and ask whether there is a dependence on H I-to-stellar mass ratio, as predicted by chemical evolution models. We do find that gas-poor galaxies are more metal rich at fixed stellar mass. We compare our results with the semi-analytic models of De Lucia & Blaizot, which include supernova feedback, as well as the cosmological infall of gas.

We present an observational search for the possible correlation between cosmic microwave background (CMB) polarization map and soft X-ray background (SXRB) based on the ROSAT All-sky Survey (RASS) archive and WMAP five-year observations. This is motivated by the fact that some of the CMB polarization may arise from the scattering of CMB photons due to the free electrons generated by X-ray heating sources in the epoch of reionization. Detection of such a correlation allows one to study the role of X-ray heating in the process of reionization. However, the cross angular power spectrum of the CMB polarization and SXRB maps constructed from ROSAT RASS and WMAP five-year maps is consistent with no correlation. We attribute this negative detection to both the extremely weak signals and large instrumental noises. While a future search is needed with high sensitivity instruments for both CMB polarization and soft X-ray photons, our current results can still be used as a useful constraint on the effect of X-ray heating in the epoch of reionization.

Motivated by recent observational studies of the environment of z~6 QSOs, we have used the Millennium Run (MR) simulations to construct a very large (~20 deg^2) mock redshift survey of star-forming galaxies at z~6. We use this simulated survey to study the relation between density enhancements in the distribution of i-dropouts and Lya-emitters, and their relation to the most massive halos and protocluster regions at z~6. Our simulation predicts significant variations in surface density across the sky with some voids and filaments extending over scales of 1 degree, much larger than probed by current surveys. Approximately one third of all z~6 halos hosting i-dropouts brighter than z=26.5 mag (~M* at z=6) become part of z=0 galaxy clusters. i-dropouts associated with protocluster regions are found in regions where the surface density is enhanced on scales ranging from a few to several tens of arcminutes on the sky. We analyze two structures of i-dropouts and Lya-emitters observed with the Subaru Telescope and show that these structures must be the seeds of massive clusters-in-formation. In striking contrast, six z~6 QSO fields observed with HST show no significant enhancements in their i-dropout number counts. With the present data, we cannot rule out the QSOs being hosted by the most massive halos. However, neither can we confirm this widely used assumption. We conclude by giving detailed recommendations for the interpretation and planning of observations by current and future ground- and space based instruments that will shed new light on questions related to the large-scale structure at z~6.

We study model predictions for three high-redshift galaxy populations: Lyman break galaxies at z~3 (LBGs), optically selected star-forming galaxies at z~2 (BXs), and distant red galaxies at z~2 (DRGs).Our galaxy formation model simultaneously reproduces the abundances, redshift distributions and clustering of all three observed populations. The star formation rates (SFRs) of model LBGs and BXs are lower than those quoted for real samples, reflecting different initial mass functions and scatter in model dust properties. About 85% of model galaxies selected as DRGs are star-forming, with SFRs ranging up to 100 M_sun/yr. Model LBGs, BXs and DRGs together account for less than half of all star formation over the range 1.5<z<3.2. Model BXs have metallicities which agree roughly with observation, but model LBGs are only slightly more metal-poor, in disagreement with recent observational results. Model galaxies are predominantly disk-dominated. About 30% of model galaxies with M>10^{11}M_sun are classified as LBGs or BXs at the relevant redshifts, while 65% are classified as DRGs. Almost all model LBGs and BXs are central galaxies, but about a quarter of DRGs are satellites. Half of all LBG descendants at z=2 would be identified as BX's, but very few as DRGs. Clustering increases with decreasing redshift for descendants of all three populations, becoming stronger than that of L^* galaxies by z=0, when many have become satellite galaxies. Their stellar mass growth is dominated by star formation until z~1 and thereafter by mergers. Most LBGs and DRGs end up as red ellipticals, while BXs have a more varied fate. 99% of local galaxies with M>10^{11.5}M_sun are predicted to have at least one LBG/BX/DRG progenitor, and over 70% above 10^{11}M_sun. (abbreviated)

We introduce a model which uses semi-analytic techniques to trace formation and evolution of galaxy disks in their cosmological context. For the first time we model the growth of gas and stellar disks separately. In contrast to previous work we follow in detail the angular momentum accumulation history through the gas cooling, merging and star formation processes. Our model successfully reproduces the stellar mass--radius distribution and gas-to-stellar disk size ratio distribution observed locally. We also investigate the dependence of clustering on galaxy size and find qualitative agreement with observation. There is still some discrepancy at small scale for less massive galaxies, indicating that our treatment of satellite galaxies needs to be improved.

We use galaxy and dark halo data from the public database for the Millennium Simulation to study the growth of galaxies in the De Lucia et al. (2006) model for galaxy formation. Previous work has shown this model to reproduce many aspects of the systematic properties and the clustering of real galaxies, both in the nearby universe and at high redshift. It assumes the stellar masses of galaxies to increase through three processes, major mergers, the accretion of smaller satellite systems, and star formation. We show the relative importance of these three modes to be a strong function of stellar mass and of redshift. Galaxy growth through major mergers depends strongly on stellar mass, but only weakly on redshift. Except for massive systems, minor mergers contribute more to galaxy growth than major mergers at all redshifts and at all stellar masses. For galaxies significantly less massive than the Milky Way, star formation dominates the growth at all epochs. For galaxies significantly more massive than the Milky Way, growth through mergers is the dominant process at all epochs. At a stellar mass of $6\times 10^{10}M_\odot$, star formation dominates at $z>1$ and mergers at later times. At every stellar mass, the growth rates through star formation increase rapidly with increasing redshift. Specific star formation rates are a decreasing function of stellar mass not only at $z=0$ but also at all higher redshifts. For comparison, we carry out a similar analysis of the growth of dark matter halos. In contrast to the galaxies, growth rates depend strongly on redshift, but only weakly on mass. They agree qualitatively with analytic predictions for halo growth.

A solution of codimension 2 brane is found for which 4 dimensional Friedmann cosmology is recovered on the brane with time dependent tension, in the Einstein frame. The effective parameter $p/\rho$ of equation of state on the brane can be quintessence like, de Sitter like or phantom like, depending on integration constants of the solution.