This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The interaction between galaxies is believed to be the main origin of the peculiarities of galaxies, which disturbs not only the morphology but also their kinematics. These disturbed and asymmetric features are the indicators of galaxy interaction. We study the velocity field of the ionized gas of the paired galaxies in the SDSS-IV MaNGA IFU survey. Using the kinemetry package, we fit the velocity field of the ionized gas to quantify the degree of kinematic asymmetry. We find that the star formation rate (SFR) of the paired galaxies with high kinematic asymmetry is significantly enhanced even when the projected separation between the pair members is quite large (dp~100 kpc). On the contrary, no significant SFR enhancement is found for the paired galaxies with low kinematic asymmetry even when their projected separation is small (dp < 30 kpc). Moreover, we also find that the fraction of galaxies with high kinematic asymmetry is much higher in close pairs (dp < 30 kpc) than those with larger dp, which explains well the early statistical findings of the significant SFR enhancement in close pairs. Our new findings illustrate that the kinematic asymmetry is an excellent indicator of galaxy-galaxy interaction strength, which helps us better understand the merging stage of the observed galaxy pairs.

Compact group (CG) is a kind of special galaxy system where the galaxy members are separated at the distances of the order of galaxy size. The strong interaction between the galaxy members makes CGs ideal labs for studying the environmental effects on galaxy evolution. Traditional photometric selection algorithm biases against the CG candidates at low redshifts, while spectroscopic identification technique is affected by the spectroscopic incompleteness of sample galaxies and typically biases against the high redshift candidates. In this study, we combine these two methods and select CGs in the main galaxy sample of the Sloan Digital Sky Survey, where we also have taken the advantages of the complementary redshift measurements from the LAMOST spectral and GAMA surveys. We have obtained the largest and most complete CG samples to date. Our samples include 6,144 CGs and 8,022 CG candidates, which are unique in the studies of the nature of the CGs and the evolution of the galaxies inside.

Nebular emission lines are critical to measure physical properties in the ionized gas (e.g., metallicity, the SFR, or dust attenuation). They also account for a significant fraction of broadband fluxes, in particular at the highest redshifts, and therefore can strongly affect the determination of other physical properties, such as the stellar mass, which are crucial in shaping our understanding of galaxy formation and evolution. We investigate a sample of 51 LBGs at 3.0<z<3.8 with detected [OIII] line emissions and estimated the physical properties of these galaxies to examine the impact of including nebular emission data in the SED fitting. We used CIGALE to fit the rest-frame UV-to-NIR SEDs of these galaxies and their emission line data simultaneously. We ran CIGALE with and without the nebular model or the emission line data, and compare the results to show the importance of including the nebular emission line data in the SED fitting. We find that without the nebular model, the fitting overestimates the stellar mass due to the strong [OIII] lines that are redshifted to the Ks-band, consistent with previous results. The emission line data are necessary to constrain the nebular model in the SED fitting. Examining the Ks-band excess, which is mostly used to estimate the emissions of the [OIII]+Hb lines when there is no spectral data, we find that the estimation and observation are statistically consistent. However, the difference can reach up to more than 1 dex in some catastrophic cases, which shows the importance of obtaining spectroscopic measurements for these lines. We also estimate the equivalent width of the Hb absorption and find it negligible compared to the Hb emission. We conclude that line emission is important to constrain the nebular models and to obtain reliable estimates of the physical properties of galaxies. These data should be taken into account in the SED fitting.

We present a novel method of magnetohydrodynamics (MHD) within the smoothed particle hydrodynamics scheme using the Geometric Density average force expression (GDSPH). GDSPH has recently been shown to reduce the leading order errors and greatly improve the accuracy near density discontinuities, eliminating surface tension effects. Here, we extend the study to investigate how SPMHD benefits from this method. We implement ideal MHD in the Gasoline2 and Changa codes with both GDSPH and traditional SPH (TSPH) schemes. An constrained hyperbolic divergence cleaning scheme is employed to control the divergence error, and a switch for artificial resistivity with minimized dissipation is used. We test the codes with a large suite of MHD tests, and show that in all problems the results are comparable or improved over previous SPMHD implementations. While both GDSPH and TSPH perform well with relatively smooth or highly supersonic flows, GDSPH shows significant improvements in the presence of strong discontinuities and large dynamic scales. In particular, when applied to an astrophysical problem of the collapse of a magnetized cloud, GDSPH realistically captures the development of a magnetic tower and jet launching in the weak-field regime, and exhibit fast convergence with resolution, while TSPH failed to do so. Our new method shows qualitatively similar results to the ones from the meshless finite mass/volume (MFM/MFV) schemes within the Gizmo code, while remaining computationally less expensive.

Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback.

We measure the bivariate luminosity function (BLF) of galaxy pairs and use it to probe and characterize the galaxy-galaxy interaction between pair members. The galaxy pair sample is selected from the main galaxy sample of Sloan Digital Sky Survey and supplied with a significant number of redshifts from the LAMOST spectral and GAMA surveys. We find the BLFs depend on the projected distance $d_{\text{p}}$ between pair members. At large separation $d_{\text{p}} > 150 h^{-1}\ \text{kpc}$, the BLF degenerates into a luminosity function (LF) of single galaxies, indicating few interactions between pair members. At $100 h^{-1}\ \text{kpc} \leq d_{\text{p}} \leq 150 h^{-1}\ \text{kpc}$, the BLF starts to show the correlation between pair members, in the sense that the shape of the conditional luminosity function (CLF) of one member galaxy starts to depend on the luminosity of the other member galaxy. Specifically, the CLF with a brighter companion has a steeper faint-end slope, which becomes even more significant at $50 h^{-1}\ \text{kpc} \leq d_{\text{p}} \leq 100 h^{-1}\ \text{kpc}$. This behavior is consistent with the scenario, \textit{and also is the observational evidence}, that dynamic friction drives massive major merger pairs to merge more quickly. At close distance $d_{\text{p}} \leq 50 h^{-1}\ \text{kpc}$, besides the merging time-scale effect, the BLF also shows an overall brightening of $\Delta M_r \geq 0.04$ mag, which reveals the enhanced star formation of the close-pair phase. By combining another statistical conclusion that the star formation rate of late-type galaxies in close pairs is enhanced at a level of about 40\%, we further conclude that the average starburst time-scale of close pairs is as long as 0.4 Gyr.

The dominant feedback mechanism in low mass haloes is usually assumed to take the form of massive stars exploding as supernovae (SNe). We perform very high resolution cosmological zoom-in simulations of five dwarf galaxies to z = 4 with our mechanical SN feedback model. This delivers the correct amount of momentum corresponding to the stage of the SN remnant evolution resolved, and has been shown to lead to realistic dwarf properties in isolated simulations. We find that in 4 out of our 5 simulated cosmological dwarfs, SN feedback has insufficient impact resulting in excessive stellar masses, extremely compact sizes and central super-solar stellar metallicities. The failure of the SN feedback in our dwarfs is physical in nature within our model and is the result of the build up of very dense gas in the early universe due to mergers and cosmic inflows prior to the first SN occurring. We demonstrate that our results are insensitive to resolution (provided that it is high enough), details of the (spatially uniform) UV background and reasonable alterations within our star formation prescription. We therefore conclude that the ability of SNe to regulate dwarf galaxy properties is dependent on other physical processes, such as turbulent pressure support, clustering and runaway of SN progenitors and other sources of stellar feedback.

We present optical luminosity functions (LFs) of galaxies in the $^{0.1}g$, $^{0.1}r$, and $^{0.1}i$ bands, calculated using data in $\sim 40$ $deg^{2}$ sky area of LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in Southern Galactic Cap. Redshifts for galaxies brighter $r = 18.1$ were obtained mainly with LAMOST. In each band, LFs derived using both parametric and non-parametric maximum likelihood methods agree well with each other. In the $^{0.1}r$ band, our fitting parameters of the Schechter function are $\phi_{*}=(1.65\pm0.36)\times10^{-2}h^{3}Mpc^{-3}$, $M_{*}=-20.69\pm0.06$ mag, and $\alpha=-1.12\pm0.08$, in agreements with previous studies. Separate LFs are also derived for emission line galaxies and absorption line galaxies, respectively. The LFs of absorption line galaxies show a dip at $^{0.1}r \sim 18.5$ and can be well fitted by a double-Gaussian function, suggesting a bi-modality in passive galaxies.

The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

The star formation history (SFH) of galaxies allow us to investigate when galaxies formed their stars and assembled their mass. We can constrain the SFH with high level of precision from galaxies with resolved stellar populations, since we are able to discriminate between stars of different ages from the spectrum they emit. However, the relative importance of secular evolution (nature) over nurture is not yet clear, and separating the effects of interaction-driven evolution in the observed galaxy properties is not trivial. The aim of this study is to use MaNGA (Mapping Nearby Galaxies at APO) Integral Field Unit (IFU) data, in combination with multi-wavelength data, to constrain the SFH of nearby isolated galaxies. We present here the new techniques we are developing to constrain the SFH with high level of precision from Spectral Energy Distribution (SED) fitting. This study is part of a China-Chile collaboration program where we are applying these new techniques to investigate how galaxies formed and evolve in different environments.

The derivation of accurate stellar populations of galaxies is a non-trivial task because of the well-known age-metallicity degeneracy. We aim to break this degeneracy by invoking a chemical evolution model(CEM) for isolated disk galaxy, where its metallicity enrichment history(MEH) is modelled to be tightly linked to its star formation history(SFH). Our CEM has been successfully tested on several local group dwarf galaxies whose SFHs and MEHs have been both independently measured from deep color-magnitude diagrams of individual stars. By introducing the CEM into the stellar population fitting algorithm as a prior, we expect that the SFH of galaxies could be better constrained.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

We utilize the data from the Apache Point Observatory Galactic Evolution Experiment-2 (APOGEE-2) in the fourteenth data release of the Sloan Digital Sky Survey (SDSS) to calculate the line-of-sight velocity dispersion $\sigma_{1D}$ of a sample of old open clusters (age larger than 100\,Myr) selected from the Milky Way open cluster catalog of Kharchenko et al. (2013). Together with their $K_s$ band luminosity $L_{K_s}$, and the half-light radius $r_{h}$ of the most probable members, we find that these three parameters show significant pairwise correlations among each other. Moreover, a fundamental plane-{\it like} relation among these parameters is found for the oldest open clusters (age older than 1\,Gyr), $L_{K_s}\propto\sigma_{1D}^{0.82\pm0.29}\cdot r_h^{2.19\pm0.52}$ with $rms \sim\, 0.31$\,mag in the $K_s$ band absolute magnitude. The existence of this relation, which deviates significantly from the virial theorem prediction, implies that the dynamical structures of the old open clusters are quite similar, when survived from complex dynamical evolution to age older than 1 Gyr.

We outline a full-scale search for galaxies exhibiting double-peaked profiles of promi- nent narrow emission lines, motivated by the prospect of finding objects related to merging galaxies, and even dual active galactic nuclei candidates as by-product, from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) Data Re- lease 4. We assemble a large sample of 325 candidates with double-peaked or strong asymmetric narrow emission lines, with 33 objects therein appearing optically resolved dual-cored structures, close companions or signs of recent interaction on the Sloan Dig- ital Sky Survey images. A candidate from LAMOST (J074810.95+281349.2) is also stressed here based on the kinematic and spatial decompositions of the double-peaked narrow emission line target, with analysis from the cross-referenced Mapping Nearby Galaxies at the Apache Point Observatory (MaNGA) survey datacube. MaNGA en- ables us to constrain the origin of double peaks for these sources, and with the IFU data we infer that the most promising origin of double-peaked profiles for LAMOST J074810.95+281349.2 is the `Rotation Dominated + Disturbance' structure.

We use high-resolution cosmological hydrodynamical simulations of Milky Way-sized galaxies with varying supernovae feedback strengths and merger histories to investigate the formation of their gaseous halos and especially their hot ($>10^6$~K) X-ray luminous coronae. Our simulations predict the presence of significant hot gas in the halos as early as $z=3-4$, well before the halos ought to be able to sustain hot mode accretion in the conventional picture. The nascent coronae grow inside-out and initially do so primarily as a result of outflows from the central galaxies powered by merger-induced shock heating and strong supernovae feedback, both of which are elemental features of today's successful galaxy formation models. Furthermore, the outflows and the forming coronae also accelerate the transition from cold to hot mode accretion by contributing to the conditions for sustaining stable accretion shocks. They also disrupt the filamentary streams funneling cold gas onto the central galaxies by causing their mouths to fray into a broad delta, detach from the galaxies, and be pushed away to larger radii. And even though at early times the filaments repeatedly re-form, the hot gas and the outflows act to weaken the filaments and accelerate their ultimate disruption. Although galactic outflows are generally thought of as ejective feedback, we find that their action on the filaments suggests a preventive role as well.

We analyze the metal accumulation in dwarf and spiral galaxies by following the history of metal enrichment and outflows in a suite of twenty high-resolution simulated galaxies. These simulations agree with the observed stellar and gas-phase mass-metallicity relation, an agreement that relies on large fractions of the produced metals escaping into the CGM. For instance, in galaxies with Mvir ~ 1e9.5 -- 1e10 solar masses, we find that about ~ 85% of the available metals are outside of the galactic disk at z = 0, although the fraction decreases to a little less than half in Milky Way-mass galaxies. In many cases, these metals are spread far beyond the virial radius. We analyze the metal deficit within the ISM and stars in the context of previous work tracking the inflow and outflow of baryons. Outflows are prevalent across the entire mass range, as is reaccretion. We find that between 40 and 80% of all metals removed from the galactic disk are later reaccreted. The outflows themselves are metal enriched relative to the ISM by a factor of 0.2 dex because of the correspondence between sites of metal enrichment and outflows. As a result, the metal mass loading factor scales as eta_metals \propto v_circ^-0.91, a somewhat shallower scaling than the total mass loading factor. We analyze the simulated galaxies within the context of analytic chemical evolution models by determining their net metal expulsion efficiencies, which encapsulate the rates of metal loss and reaccretion. We discuss these results in light of the inflow and outflow properties necessary for reproducing the mass-metallicity relation.

It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation.

We present the Feedback Acting on Baryons in Large-scale Environments (FABLE) suite of cosmological hydrodynamical simulations of galaxies, groups and clusters. The simulations use the AREPO moving-mesh code with a set of physical models for galaxy formation based on the successful Illustris simulation, but with updated AGN and supernovae feedback models. This allows us to simultaneously reproduce the observed redshift evolution of the galaxy stellar mass function together with the stellar and gas mass fractions of local groups and clusters across a wide range of halo masses. Focusing on the properties of groups and clusters, we find very good agreement with a range of observed scaling relations, including the X-ray luminosity--total mass and gas mass relations as well as the total mass--temperature and Sunyaev-Zel'dovich flux--mass relations. Careful comparison of our results with scaling relations based on X-ray hydrostatic masses as opposed to weak lensing-derived masses reveals some discrepancies, which hint towards a non-negligible X-ray mass bias in observed samples. We further show that radial profiles of density, pressure and temperature of the simulated intracluster medium are in very good agreement with observations, in particular for $r > 0.3\ r_{500}$. In the innermost regions however we find too large entropy cores, which indicates that a more sophisticated modelling of the physics of AGN feedback may be required to accurately reproduce the observed populations of cool-core and non-cool-core clusters.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

While feedback from massive stars exploding as supernovae (SNe) is thought to be one of the key ingredients regulating galaxy formation, theoretically it is still unclear how the available energy couples to the interstellar medium and how galactic scale outflows are launched. We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections of thermal and/or kinetic energy, two parametrizations of delayed cooling feedback and a `mechanical' feedback scheme that injects the correct amount of momentum depending on the relevant scale of the SN remnant resolved. All schemes make use of individually time-resolved SN events. Adopting isolated disk galaxy setups at different resolutions, with the highest resolution runs reasonably resolving the Sedov-Taylor phase of the SN, we aim to find a physically motivated scheme with as few tunable parameters as possible. As expected, simple injections of energy overcool at all but the highest resolution. Our delayed cooling schemes result in overstrong feedback, destroying the disk. The mechanical feedback scheme is efficient at suppressing star formation, agrees well with the Kennicutt-Schmidt relation and leads to converged star formation rates and galaxy morphologies with increasing resolution without fine tuning any parameters. However, we find it difficult to produce outflows with high enough mass loading factors at all but the highest resolution, indicating either that we have oversimplified the evolution of unresolved SN remnants, require other stellar feedback processes to be included, require a better star formation prescription or most likely some combination of these issues.

We present a three-dimensional modeling of the Milky Way dust distribution by fitting the value-added star catalog of LAMOST spectral survey. The global dust distribution can be described by an exponential disk with scale-length of 3,192 pc and scale height of 103 pc. In this modeling, the Sun is located above the dust disk with a vertical distance of 23 pc. Besides the global smooth structure, two substructures around the solar position are also identified. The one located at $150^{\circ}<l<200^{\circ}$ and $-5^{\circ}<b<-30^{\circ}$ is consistent with the Gould Belt model of \citet{Gontcharov2009}, and the other one located at $140^{\circ}<l<165^{\circ}$ and $0^{\circ}<b<15^{\circ}$ is associated with the Camelopardalis molecular clouds.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We investigate the star formation history and the dust attenuation in the galaxy merger Mrk848. Thanks to the multiwavelength photometry from the ultraviolet (UV) to the infrared (IR), and MaNGA's integral field spectroscopy, we are able to study this merger in a detailed way. We divide the whole merger into the core and tail regions, and fit both the optical spectrum and the multi-band spectral energy distribution (SED) to models to obtain the star formation properties for each region respectively. We find that the color excess of stars in the galaxy $E(B-V)_s^{SED}$ measured with the multi-band SED fitting is consistent with that estimated both from the infrared excess (the ratio of IR to UV flux) and from the slope of the UV continuum. Furthermore, the reliability of the $E(B-V)_s^{SED}$ is examined with a set of mock SEDs, showing that the dust attenuation of the stars can be well constrained by the UV-to-IR broadband SED fitting. The dust attenuation obtained from optical continuum $E(B-V)_s^{spec}$ is only about half of $E(B-V)_s^{SED}$. The ratio of the $E(B-V)_s^{spec}$ to the $E(B-V)_g$ obtained from the Balmer decrement is consistent with the local value (around 0.5). The difference between the results from the UV-to-IR data and the optical data is consistent with the picture that younger stellar populations are attenuated by an extra dust component from the birth clouds compared to older stellar populations which are only attenuated by the diffuse dust. Both with the UV-to-IR SED fitting and the spectral fitting, we find that there is a starburst younger than 100~Myr in one of the two core regions, consistent with the scenario that the interaction-induced gas inflow can enhance the star formation in the center of galaxies.

We study the classification and composite spectra of galaxy in the fourth data release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). We select 40,182 spectra of galaxies from LAMOST DR4, which have photometric in- formation but no spectroscopic observations in the Sloan Digital Sky Survey(SDSS). These newly observed spectra are re-calibrated and classified into six classes, i.e. pas- sive, H{\alpha}-weak, star-forming, composite, LINER and Seyfert using the line intensity (H\b{eta}, [OIII]{\lambda}5007, H{\alpha} and [NII]{\lambda}6585). We also study the correlation between spectral classes and morphological types through three parameters: concentration index, (u - r) color, and D4000n index. We calculate composite spectra of high signal-to-noise ra- tio(S/N) for six spectral classes, and using these composites we pick out some features that can differentiate the classes effectively, including H\b{eta}, Fe5015, H{\gamma}A, HK, and Mg2 band etc. In addition, we compare our composite spectra with the SDSS ones and analyse their difference. A galaxy catalogue of 40,182 newly observed spectra (36,601 targets) and the composite spectra of the six classes are available online.

We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spiral galaxies may be conducive to producing a rich population of massive binary black holes. We consider here both detailed controlled simulations of spirals and cosmological simulations, as well as the current range of observed values for metallicity and star formation in outer disks. We find that outer HI disks contribute at least as much as dwarf galaxies do to the observed LIGO/Virgo detection rates. Identifying the host galaxies of merging massive black holes should provide constraints on cosmological parameters and insights into the formation channels of binary mergers.

We present a spectroscopic redshift catalog from the LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in the Southern Galactic Cap (SGC), which is designed to observe all sources (Galactic and extra-galactic) by using repeating observations with a limiting magnitude of $r=18.1~mag$ in two $20~deg^2$ fields. The project is mainly focusing on the completeness of LAMOST ExtraGAlactic Surveys (LEGAS) in the SGC, the deficiencies of source selection methods and the basic performance parameters of LAMOST telescope. In both fields, more than 95% of galaxies have been observed. A post-processing has been applied to LAMOST 1D spectrum to remove the majority of remaining sky background residuals. More than 10,000 spectra have been visually inspected to measure the redshift by using combinations of different emission/absorption features with uncertainty of $\sigma_{z}/(1+z)<0.001$. In total, there are 1528 redshifts (623 absorption and 905 emission line galaxies) in Field A and 1570 redshifts (569 absorption and 1001 emission line galaxies) in Field B have been measured. The results show that it is possible to derive redshift from low SNR galaxies with our post-processing and visual inspection. Our analysis also indicates that up to 1/4 of the input targets for a typical extra-galactic spectroscopic survey might be unreliable. The multi-wavelength data analysis shows that the majority of mid-infrared-detected absorption (91.3%) and emission line galaxies (93.3%) can be well separated by an empirical criterion of $W2-W3=2.4$. Meanwhile, a fainter sequence paralleled to the main population of galaxies has been witnessed both in $M_r$/$W2-W3$ and $M_*$/$W2-W3$ diagrams, which could be the population of luminous dwarf galaxies but contaminated by the edge-on/highly inclined galaxies ($\sim30\%$).

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Using the Eris zoom-in cosmological simulation of assembly of a Milky Way analog, we study chemical enrichment of stars due to accretion of metal-enriched gas from the interstellar medium during the Galaxy's development. We consider metal-poor and old stars in both Galactic halo and bulge and make use of stellar orbits, gas density and metallicity distributions in Eris. Assuming spherically symmetric Bondi-Hoyle accretion, we find that halo and bulge stars accrete metals at the rate of about 10^-24 solar mass per year and 10^-22 solar mass per year, respectively, at redshifts z < 3, but this accretion rate increases hundred-fold to about 10^-20 solar mass per year at higher redshifts due to increased gas density. Bulge and halo stars accrete similar amounts of metals at high redshifts as kinematically distinct bulge and halo are not yet developed at these redshifts and both sets of stars encounter similar metal distribution in the ISM on average. Accretion alone can enrich main-sequence stars up to [Fe/H] -2 in extreme cases.Median enrichment level due to accretion in these stars is about [Fe/H]~-6 to -5.Because accretion mostly takes place at high redshifts, it is alpha-enriched to [alpha/Fe]~0.5. We find that accretive metal enrichment is significant enough to affect the predicted metallicity distribution function of halo stars at [Fe/H] < -5.This suggests that attempts to infer the natal chemical environment of the most metal-poor stars from their observed enrichment today can be hindered due to metal accretion. Peculiar enrichment patterns such as those predicted to arise from pair-instability supernovae could help in disentangling natal and accreted metal content of stars.

We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major and minor mergers hosting at least one AGN, the fraction harbouring dual AGN is ~20-30 and ~1-10 per cent, respectively. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a shallow imaging survey will require very high angular resolution, whereas a deep imaging survey will be less resolution-dependent.

We confirm that the object DDO216-A1 is a substantial globular cluster at the center of Local Group galaxy DDO216 (the Pegasus dwarf irregular), using Hubble Space Telescope ACS imaging. By fitting isochrones, we find the cluster metallicity to be [M/H] = -1.6 +/-0.2, for reddening E(B-V) = 0.16 +/-0.02; the best-fit age is 12.3 +/-0.8 Gyr. There are ~30 RR Lyrae variables in the cluster; the magnitude of the fundamental mode pulsators gives a distance modulus of 24.77 +/-0.08 - identical to the host galaxy. The ratio of overtone to fundamental mode variables and their mean periods make DDO216-A1 an Oosterhoff Type I cluster. We find an I-band central surface brightness 20.85 +/-0.17 F814W mag per square arcsecond, a half-light radius of 3.1 arcsec (13.4 pc), and an absolute magnitude M814 = -7.90 +/-0.16 (approximately 10^5 solar masses). King models fit to the cluster give the core radius and concentration index, r_c = 2.1" +/-0.9" and c = 1.24 +/-0.39. The cluster is an "extended" cluster somewhat typical of some dwarf galaxies and the outer halo of the Milky Way. The cluster is projected <30 pc south of the center of DDO216, unusually central compared to most dwarf galaxy globular clusters. Analytical models of dynamical friction and tidal destruction suggest that it probably formed at a larger distance, up to ~1 kpc, and migrated inward. DDO216 has an unexceptional cluster specific frequency, S_N = 10. DDO216 is the lowest-luminosity Local Group galaxy to host a 10^5 solar mass globular cluster, and the only transition-type (dSph/dIrr) in the Local Group with a globular.

We study the integral Galactic extinction and reddening based on the galaxy catalog of the South Galactic Cap U-band Sky Survey (SCUSS), where $u$ band galaxy number counts and $u-r$ color distribution are used to derive the Galactic extinction and reddening respectively. We compare these independent statistical measurements with the reddening map of \citet{Schlegel1998}(SFD) and find that both the extinction and reddening from the number counts and color distribution are in good agreement with the SFD results at low extinction regions ($E(B-V)^{SFD}<0.12$ mag). However, for high extinction regions ($E(B-V)^{SFD}>0.12$ mag), the SFD map overestimates the Galactic reddening systematically, which can be approximated by a linear relation $\Delta E(B-V)= 0.43[E(B-V)^{SFD}-0.12$]. By combing the results of galaxy number counts and color distribution together, we find that the shape of the Galactic extinction curve is in good agreement with the standard $R_V=3.1$ extinction law of \cite{ODonnell1994}.

We present and analyze the possibility of using optical ${\it u}$-band luminosities to estimate star-formation rates (SFRs) of galaxies based on the data from the South Galactic Cap ${\it u }$ band Sky Survey (SCUSS), which provides a deep ${\it u}$-band photometric survey covering about 5000 $deg^2$ of the South Galactic Cap. Based on two samples of normal star-forming galaxies selected by the BPT diagram, we explore the correlations between ${\it u}$-band, H$\alpha$, and IR luminosities by combing SCUSS data with the Sloan Digital Sky Survey (SDSS) and ${\it Wide}$-${\it field\ Infrared\ Survey\ Explorer}$ (${\it WISE}$). The attenuation-corrected ${\it u}$-band luminosities are tightly correlated with the Balmer decrement-corrected H$\alpha$ luminosities with an rms scatter of $\sim$ 0.17 dex. The IR-corrected ${\it u }$ luminosities are derived based on the correlations between the attenuation of ${\it u}$-band luminosities and ${\it WISE}$ 12 (or 22) $\mu$m luminosities, and then calibrated with the Balmer-corrected H$\alpha$ luminosities. The systematic residuals of these calibrations are tested against the physical properties over the ranges covered by our sample objects. We find that the best-fitting nonlinear relations are better than the linear ones and recommended to be applied in the measurement of SFRs. The systematic deviations mainly come from the pollution of old stellar population and the effect of dust extinction; therefore, a more detailed analysis is needed in the future work.

We investigate the angular momentum evolution of four disk galaxies residing in Milky Way-sized halos formed in cosmological zoom-in simulations with various sub-grid physics and merging histories. We decompose these galaxies kinematically and photometrically, into their disk and bulge components. The simulated galaxies and their components lie on the observed sequences in the $j_*$--$M_*$ diagram relating the specific angular momentum and mass of the stellar component. We find that galaxies in low-density environments follow the relation $j_* \propto M_*^{\alpha}$ past major mergers, with $\alpha \sim 0.6$ in the case of strong feedback, when bulge-to-disk ratios are relatively constant, and $\alpha \sim 1.4$ in the other cases, when secular processes operate on shorter timescales. We compute the retention factors (i.e. the ratio of the specific angular momenta of stars and dark matter) for both disks and bulges and show that they vary relatively slowly after averaging over numerous but brief fluctuations. For disks, the retention factors are usually close to unity, while for bulges, they are a few times smaller. Our simulations therefore indicate that galaxies and their halos grow in a quasi-homologous way.

Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends sensitively on the progenitor's metallicity. Rare low-metallicity star formation during galaxy assembly can produce more detected compact binaries than typical star formation. Using detailed simulations of galaxy and chemical evolution, we determine how sensitively the compact binary populations of galaxies with similar present-day appearance depend on the details of their assembly. We also demonstrate by concrete example the extent to which dwarf galaxies overabundantly produce compact binary mergers, particularly binary black holes, relative to more massive galaxies. We discuss the implications for transient multimessenger astronomy with compact binary sources.

Active galactic nuclei (AGN) is one of the main drivers for transition from star-forming disk to passive spheroidal galaxies. However, the role of large-scale environment versus one-on-one interactions in triggering different types of AGN is still uncertain. We present a statistical study of the prevalence of the nuclear activity in isolated galaxies and physically bound isolated pairs. For the purpose of this study we considered optically and radio selected nuclear activity types. We aim to assess the effect of one-on-one interaction on the fraction of AGN and the role of their large-scale environment. To study the effect of one-on-one interaction on the fraction of AGN in isolated galaxy pairs, we compare with a sample of isolated galaxies homogeneously selected under the same isolation criterion. We examine the effect of the large-scale environment by comparing with control samples of single galaxies and galaxy pairs. In general we found no difference in the prevalence of optical AGN for the considered samples. For massive galaxies, the fraction of optical AGN in isolated galaxies is slightly higher than that in control samples. Also the fraction of passives in high mass isolated galaxies is smaller than in any other sample. Generally, there is no dependence on optical nuclear activity with local environment. On the other hand, we found evidence that radio AGN are strongly affected by the local environment. Optical AGN phenomenon is related to cold gas accretion, while radio AGN is related to hot gas accretion. In this context, there is more cold gas, fueling the central optical AGN, in isolated systems. Our results are in agreement with a scenario where cold gas accretion by secular evolution is the main driver of optical AGN, while hot gas accretion and one-on-one interactions are the main drivers of radio AGN activity.

We investigate the UV continuum slope $\alpha$ of a large quasar sample from SDSS DR7. By using specific continuum windows, we build two samples at lower ($0.71<z<1.19$) and higher ($1.90<z<3.15$) redshifts, which correspond to the continuum slopes at longer (NUV) and shorter (FUV) rest wavelength ranges respectively. Overall, the average continuum slopes are $-0.36$ and $-0.51$ for $\alpha_{\rm NUV}$ and $\alpha_{\rm FUV}$ with similar dispersions $\sigma_{\alpha} \sim 0.5$. For both samples, we confirm the luminosity dependence of the continuum slope, i.e., fainter quasars have redder spectra. We further find that both $\alpha_{\rm NUV}$ and $\alpha_{\rm FUV}$ have a common upper limit ($\sim 1/3$) which is almost independent of the quasar luminosity $L_{\rm bol}$. This finding implies that the intrinsic quasar continuum (or the bluest quasar), at any luminosity, obey the standard thin disk model. We propose that the other quasars with redder $\alpha$ are caused by the reddening from the dust {\it locally}. With this assumption, we employ the dust extinction scenario to model the observed $L_{\rm bol}-\alpha$ relation. We find that, a typical value of $E(B-V)\sim0.1$ to $0.3$ mag (depending on the types of extinction curve) of the quasar {\it local} dust is enough to explain the discrepancy of $\alpha$ between the observation ($\sim-0.5$) and the standard accretion disk model prediction ($\sim 1/3$).

Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We derive the distributions of time between gas crossing the virial radius and being accreted to the star forming region (which allows mixing within the corona), as well as the time between gas being accreted to the star forming region and then forming stars (which allows mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star forming region. Gas entering high mass galaxies is pre-enriched in low mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z~5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star forming region, with gas that is feeding the star forming region surpassing 0.1Z by z=0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulation (Mvir~2x10^{10}Msun, with M*~10^7Msun), by contrast, accretes primordial gas through the virial radius and onto the disc at all times. Much like the classical analytical solutions to the `G-dwarf problem', overproduction of low-metallicity stars is ameliorated by the inefficiency of star formation. Finally, gas outflow/metal removal rates from star forming regions as a function of galactic mass are presented.

The r-process nuclei are robustly synthesized in the material ejected during a neutron star binary merger (NSBM), as tidal torques transport angular momentum and energy through the outer Lagrange point in the form of a vast tidal tail. If NSBM are indeed solely responsible for the solar system r- process abundances, a galaxy like our own would require to host a few NSBM per million years, with each event ejecting, on average, about 5x10^{-2} M_sun of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than in ordinary supernovae, NSBM deposit a comparable amount of energy into the interstellar medium (ISM). In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution in a cosmological context, we perform a suite of three-dimensional hydrodynamic simulations with optically-thin radiative cooling of isolated NSBM ejecta expanding in environments with conditions adopted from Milky Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution that results from thermal instability in atomic gas is remarkably similar to that of a standard supernova blast wave. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBM blast waves. Among other quantities, we examine the radius, time, mass and kinetic energy content of the NSBM remnant at shell formation as well as the momentum injected to the ISM. We find that the shell formation epoch is attained when the swept-up mass is about 10^3 M_sun, at this point the mass fraction of r-process material is drastically enhanced up to two orders of magnitude in relation to a solar metallicity ISM.

We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $\beta$-model of Miller & Bregman (2013) based on OVII absorption and OVIII emission measurements. Warm--hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram--pressure stripping of the gas reservoirs in dwarf galaxy satellites.

The SCUSS is a deep $u$-band imaging survey in the south Galactic cap using the 2.3m Bok telescope. The survey observations were completed in the end of 2013, covering an area of about 5000 square degrees. We release the data in the region with an area of about 4000 deg$^2$ that is mostly covered by the Sloan digital sky survey. The data products contain calibrated single-epoch images, stacked images, photometric catalogs, and a catalog of star proper motions derived by Peng et al, 2015. The median seeing and magnitude limit ($5\sigma$) are about 2".0 and 23.2 mag, respectively. There are about 8 million objects having measurements of absolute proper motions. All the data and related documentations can be accessed through the SCUSS data release website of \url{http://batc.bao.ac.cn/Uband/data.html}.

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

We report the discovery of an isolated compact galaxy triplet SDSS J084843.45+164417.3, which is first detected by the LAMOST spectral survey and then confirmed by the spectroscopic observation of the BFOSC of the 2.16 meter telescope. It is found that this triplet is an isolated and extremely compact system, which has an aligned configuration and very small radial velocity dispersion. The member galaxies have similar colors and show marginal star formation activities. These results enhance the opinion that the compact triplets are well-evolved systems rather than the hierarchically forming structures. This occasional discovery reveals the limitations of the fiber spectral redshift surveys in studying such compact system, and declares the necessity of additional observations to complete the current redshift sample.

A small fraction($<10\%$) of SDSS main sample galaxies(MGs) have not been targeted with spectroscopy due to the the fiber collision effect. These galaxies have been compiled into the input catalog of the LAMOST extra-galactic survey and named as the complementary galaxy sample. In this paper, we introduce the project and the status of the spectroscopies of the complementary galaxies in the first two years of the LAMOST spectral survey(till Sep. of 2014). Moreover, we present a sample of 1,102 galaxy pairs identified from the LAMOST complementary galaxies and SDSS MGs, which are defined as that the two members have a projected distance smaller than 100 kpc and the recessional velocity difference smaller than 500 $\rm kms^{-1}$. Compared with the SDSS only selected galaxy pairs, the LAMOST-SDSS pairs take the advantages of not being biased toward large separations and therefor play as a useful supplement to the statistical studies of galaxy interaction and galaxy merging.

We present preliminary results of the quasar survey in Large Sky Area Multi- Object Fiber Spectroscopic Telescope (LAMOST) first data release (DR1), which includes pilot survey and the first year regular survey. There are 3921 quasars identified with reliability, among which 1180 are new quasars discovered in the survey. These quasars are at low to median redshifts, with highest z of 4.83. We compile emission line measurements around the H{\alpha}, H{\beta}, Mg II, and C IV regions for the new quasars. The continuum luminosities are inferred from SDSS photo- metric data with model fitting as the spectra in DR1 are non-flux-calibrated. We also compile the virial black hole mass estimates, and flags indicating the selec- tion methods, broad absorption line quasars. The catalog and spectra for these quasars are available online. 28% of the 3921 quasars are selected with optical- infrared colours independently, indicating that the method is quite promising in completeness of quasar survey. LAMOST DR1 and the on-going quasar survey will provide valuable data in the studies of quasars.

We present preliminary results of the spectral analysis on the radial distributions of the star formation history in both, a galaxy merger and a spiral isolated galaxy observed with MaNGA. We find that the central part of the isolated galaxy is composed by older stellar population ($\sim$2 Gyr) than in the outskirts ($\sim$7 Gyr). Also, the time-scale is gradually larger from 1 Gyr in the inner part to 3 Gyr in the outer regions of the galaxy. In the case of the merger, the stellar population in the central region is older than in the tails, presenting a longer time-scale in comparison to central part in the isolated galaxy. Our results are in agreement with a scenario where spiral galaxies are built from inside-out. In the case of the merger, we find evidence that interactions enhance star formation in the central part of the galaxy.

The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic Cap, using the 90Prime wide-field imager on the 2.3m Bok telescope at Kitt Peak. The survey observations started in 2010 and ended in 2013. The final survey area is about 5000 deg2 with a median 5-sigma point source limiting magnitude of about 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS photometric catalog with those of the SDSS and CFHTLS.

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

We study the effect of sub-grid physics, galaxy mass, structural parameters and resolution on the fragmentation of gas-rich galaxy discs into massive star forming clumps. The initial conditions are set up with the aid of the ARGO cosmological hydrodynamical simulation. Blast-wave feedback does not suppress fragmentation, but reduces both the number of clumps and the duration of the unstable phase. Once formed, bound clumps cannot be destroyed by our feedback model. Widespread fragmentation is promoted by high gas fractions and low halo concentrations. Yet giant clumps $M > 10^8 M_{\odot}$ lasting several hundred Myr are rare and mainly produced by clump-clump mergers. They occur in massive discs with maximum rotational velocities $V_{max} > 250$ km/s at $z \sim 2$, at the high mass end of the observed galaxy population at those redshifts. The typical gaseous and stellar masses of clumps in all runs are in the range $\sim 10^7-10^8 M_{\odot}$ for galaxies with disc mass in the range $10^{10}-8\times 10^{10} M_{\odot}$. Clumps sizes are usually in the range $100-400$ pc, in agreement with recent clump observations in lensed high-z galaxies. \\ We argue that many of the giant clumps identified in observations are not due to in-situ fragmetation, or are the result of blending of smaller structures owing to insufficient resolution. Using an analytical model describing local collapse inside spiral arms, we can predict the characteristic gaseous masses of clumps at the onset of fragmentation ($\sim 3-5 \times 10^7 M_{\odot}$) quite accurately, while the conventional Toomre mass overestimates them. Due to their moderate masses, clumps which migrate to the centre have marginal effect on bulge growth.

We develop a template-fit method to automatically identify and classify late-type K and M dwarfs in spectra from the LAMOST. A search of the commissioning data, acquired in 2009-2010, yields the identification of 2612 late-K and M dwarfs. The template fit method also provides spectral classification to half a subtype, classifies the stars along the dwarf-subdwarf metallicity sequence, and provides improved metallicity/gravity information on a finer scale. The automated search and classification is performed using a set of cool star templates assembled from the Sloan Digital Sky Survey spectroscopic database. We show that the stars can be efficiently classified despite shortcomings in the LAMOST commissioning data which include bright sky lines in the red. In particular we find that the absolute and relative strengths of the critical TiO and CaH molecular bands around 7000A are cleanly measured, which provides accurate spectral typing from late-K to mid-M, and makes it possible to estimate metallicities in a way that is more efficient and reliable than with the use of spectral indices or spectral-index based parameters such as zeta. Most of the cool dwarfs observed by LAMOST are found to be metal-rich dwarfs. We use a calibration of spectral type to absolute magnitude and estimate spectroscopic distances for all the stars; we also recover proper motions from the SUPERBLINK and PPMXL catalogs. Our analysis of the estimated transverse motions suggests a mean velocity and standard deviation for the UVW components of velocity to be: U=-9.8 km/s; V=-22.8 km/s; W=-7.9 km/s. The resulting values are general agreement with previous reported results, which yields confidence in our spectral classification and spectroscopic distance estimates, and illustrates the potential for using LAMOST spectra of K and M dwarfs for investigating the chemo-kinematics of the local Galactic disk and halo.

Based on SDSS g, r and SCUSS (South Galactic Cap of u-band Sky Survey) $u$ photometry, we develop a photometric calibration for estimating the stellar metallicity from $u-g$ and $g-r$ colors by using the SDSS spectra of 32,542 F- and G-type main sequence stars, which cover almost $3700$ deg$^{2}$ in the south Galactic cap. The rms scatter of the photometric metallicity residuals relative to spectrum-based metallicity is $0.14$ dex when $g-r<0.4$, and $0.16$ dex when $g-r>0.4$. Due to the deeper and more accurate magnitude of SCUSS $u$ band, the estimate can be used up to the faint magnitude of $g=21$. This application range of photometric metallicity calibration is wide enough so that it can be used to study metallicity distribution of distant stars. In this study, we select the Sagittarius (Sgr) stream and its neighboring field halo stars in south Galactic cap to study their metallicity distribution. We find that the Sgr stream at the cylindrical Galactocentric coordinate of $R\sim 19$ kpc, $\left| z\right| \sim 14$ kpc exhibits a relative rich metallicity distribution, and the neighboring field halo stars in our studied fields can be modeled by two-Gaussian model, with peaks respectively at [Fe/H]$=-1.9$ and [Fe/H]$=-1.5$.

We use the Herschel SPIRE color-color diagram to study the spectral energy distribution (SED) and the redshift estimation of high-z galaxies. We compiled a sample of 57 galaxies with spectroscopically confirmed redshifts and SPIRE detections in all three bands at $z=2.5-6.4$, and compared their average SPIRE colors with SED templates from local and high-z libraries. We find that local SEDs are inconsistent with high-z observations. The local calibrations of the parameters need to be adjusted to describe the average colors of high-z galaxies. For high-z libraries, the templates with an evolution from z=0 to 3 can well describe the average colors of the observations at high redshift. Using these templates, we defined color cuts to divide the SPIRE color-color diagram into different regions with different mean redshifts. We tested this method and two other color cut methods using a large sample of 783 Herschel-selected galaxies, and find that although these methods can separate the sample into populations with different mean redshifts, the dispersion of redshifts in each population is considerably large. Additional information is needed for better sampling.

A sample of 70 E+A galaxies are selected from 37, 206 galaxies in the second data release (DR2) of Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) according to the criteria for E+A galaxies defined by Goto, and each of these objects is further visually identified. In this sample, most objects are low redshift E+A galaxies with z < 0.25, and locate in the high latitude sky area with magnitude among 14 to 18 mag in g, r and i bands. A stellar population analysis for the whole sample indicates that the E+A galaxies are characterized by both young and old stellar populations (SPs), and the metal-rich SPs have relatively higher contributions than the metal-poor ones. Additionally, a morphological classification for these objects is performed based on the images taken from the Sloan Digital Sky Survey (SDSS).

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

We investigate the luminosity and redshift dependence of the quasar continuum by means of composite spectrum using a large non-BAL radio-quiet quasar sample drawn from the Sloan Digital Sky Survey. Quasar continuum slopes in the UV-Opt band are measured at two different wavelength ranges, i.e., $\alpha_{\nu12}$ ($1000\sim 2000 \rm\AA$) and $\alpha_{\nu24}$ ($2000 \sim 4000 \rm\AA$) derived from power law fitting. Generally, the UV spectra slope becomes harder (higher $\alpha_{\nu}$) towards higher bolometric luminosity. On the other hand, when quasars are further grouped into luminosity bins, we find both $\alpha_{\nu12}$ and $\alpha_{\nu24}$ show significant anti-correlation with redshift (i.e., quasar continuum becomes redder towards higher redshift). We suggest that the cosmic dust extinction is very likely the cause of this observed $\alpha_\nu-z$ relation. We build a simple cosmic dust extinction model to quantify the observed reddening tendency and find an effective dust density $n\sigma_v \sim 10^{-5}h~\rm Mpc^{-1}$ at $z<1.5$. The other possibilities that could produce such a reddening effect have also been discussed.

We investigate the production sites and the enrichment history of $r$-process elements in the Galaxy, as traced by the [Eu/Fe] ratio, using the high resolution, cosmological zoom-in simulation `Eris'. At $z=0$, Eris represents a close analog to the Milky Way, making it the ideal laboratory to understand the chemical evolution of our Galaxy. Eris formally traces the production of oxygen and iron due to Type-Ia and Type-II supernovae. We include in post-processing the production of $r$-process elements from compact binary mergers. Unlike previous studies, we find that the nucleosynthetic products from compact binary mergers can be incorporated into stars of very low metallicity and at early times, even with a minimum delay time of 100 Myr. This conclusion is relatively insensitive to modest variations in the merger rate, minimum delay time, and the delay time distribution. By implementing a first-order prescription for metal-mixing, we can further improve the agreement between our model and the data for the chemical evolution of both [$\alpha$/Fe] and [Eu/Fe].We argue that compact binary mergers could be the dominant source of $r$-process nucleosynthesis in the Galaxy.

Reliable halo mass estimation for a given galaxy system plays an important role both in cosmology and galaxy formation studies. Here we set out to find the way that can improve the halo mass estimation for those galaxy systems with limited brightest member galaxies been observed. Using four mock galaxy samples constructed from semi-analytical formation models, the subhalo abundance matching method and the conditional luminosity functions, respectively, we find that the luminosity gap between the brightest and the subsequent brightest member galaxies in a halo (group) can be used to significantly reduce the scatter in the halo mass estimation based on the luminosity of the brightest galaxy alone. Tests show that these corrections can significantly reduce the scatter in the halo mass estimations by $\sim 50\%$ to $\sim 70\%$ in massive halos depending on which member galaxies are considered. Comparing to the traditional ranking method, we find that this method works better for groups with less than five members, or in observations with very bright magnitude cut.

Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent.

The South Galactic Cap $u$-band Sky Survey (SCUSS) provides a deep $u$-band imaging of about 5000 deg$^2$ in south Galactic cap. It is about 1.5 mag deeper than the SDSS $u$-band. In this paper we evaluate the capability of quasar selection using both SCUSS and SDSS data, based on considerations of the deep SCUSS $u$-band imaging and two-epoch $u$-band variability. We find that the combination of the SCUSS $u$-band and the SDSS $griz$ band allows us to select more faint quasars and more quasars at redshift around 2.2 than the selection only with the SDSS $ugriz$ data. Quasars have significant $u$-band variabilities. The fraction of quasars with large two-epoch variability is much higher than that of stars. The selection by variability can select both low-redshift quasars with ultraviolet excess and mid-redshift ($2 < z <3.5$) quasars where quasar selection by optical colors is inefficient. The above two selections are complementary and make full use of the SCUSS u-band advantages.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

We present novel 3D multi-scale SPH simulations of gas-rich galaxy mergers between the most massive galaxies at $z \sim 8 - 10$, designed to scrutinize the direct collapse formation scenario for massive black hole seeds proposed in \citet{mayer+10}. The simulations achieve a resolution of 0.1 pc, and include both metallicity-dependent optically-thin cooling and a model for thermal balance at high optical depth. We consider different formulations of the SPH hydrodynamical equations, including thermal and metal diffusion. When the two merging galaxy cores collide, gas infall produces a compact, optically thick nuclear disk with densities exceeding $10^{-10}$ g cm$^3$. The disk rapidly accretes higher angular momentum gas from its surroundings reaching $\sim 5$ pc and a mass of $\gtrsim 10^9$ $M_{\odot}$ in only a few $10^4$ yr. Outside $\gtrsim 2$ pc it fragments into massive clumps. Instead, supersonic turbulence prevents fragmentation in the inner parsec region, which remains warm ($\sim 3000-6000$ K) and develops strong non-axisymmetric modes that cause prominent radial gas inflows ($> 10^4$ $M_{\odot}$ yr$^{-1}$), forming an ultra-dense massive disky core. Angular momentum transport by non-axisymmetric modes should continue below our spatial resolution limit, quickly turning the disky core into a supermassive protostar which can collapse directly into a massive black hole of mass $10^8-10^9$ $M_{\odot}$ via the relativistic radial instability. Such a "cold direct collapse"' explains naturally the early emergence of high-z QSOs. Its telltale signature would be a burst of gravitational waves in the frequency range $10^{-4} - 10^{-1}$ Hz, possibly detectable by the planned eLISA interferometer.

We examine gas accretion and subsequent star formation in representative galaxies from the McMaster Unbiased Galaxy Simulations (Stinson et al. 2010). Accreted gas is bimodal with a natural temperature division at $10^5$ K, near the peak of the cooling curve. Cold-mode accretion dominates inflows at early times, creating a peak in total accretion at redshift z=2-4 and declining exponentially below z$\sim$2. Hot-mode accretion peaks near z=1-2 and declines gradually. Hot-mode exceeds cold-mode accretion at z$\sim$1.8 for all four galaxies rather than when the galaxy reaches a characteristic mass. Cold-mode accretion can fuel immediate star formation, while hot-mode accretion preferentially builds a large, hot gas reservoir in the halo. Late-time star formation relies on reservoir gas accreted 2-8 Gyr prior. Thus, the reservoir allows the star formation rate to surpass the current overall gas accretion rate. Stellar feedback cycles gas from the interstellar medium back into the hot reservoir. Stronger feedback results in more gas cycling, gas removal in a galactic outflow and less star formation overall, enabling simulations to match the observed star formation history. For lower mass galaxies in particular, strong feedback can delay the star formation peak to z=1-2 from the accretion peak at z=2-4.

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

We present more results from a fully cosmological LCDM simulation of a group of isolated dwarf galaxies that has been shown to reproduce the observed stellar mass and cold gas content, resolved star formation histories, and metallicities of dwarfs in the Local Volume. Here we investigate the energetics and timetable of the cusp-core transformation. As suggested by previous work, supernova-driven gas outflows remove dark matter (DM) cusps and create kpc-size cores in all systems having a stellar mass above 1e6 Msun. The "DM core mass removal efficiency" -- dark mass ejected per unit stellar mass -- ranges today from a few to a dozen, and increases with decreasing host mass. Because dwarfs form the bulk of their stars prior to redshift 1 and the amount of work required for DM heating and core formation scales approximately as Mvir^{5/3}, the unbinding of the DM cusp starts early and the formation of cored profiles is not as energetically onerous as previously claimed. DM particles in the cusp typically migrate to 2-3 core radii after absorbing a few percent of the energy released by supernovae. The present-day slopes of the inner dark matter mass profiles, Gamma=d log M/dlog R=2.5-3, of the simulated "Bashful" and "Doc" dwarfs are similar to those measured in the luminous Fornax and Sculptor dwarf spheroidals. None of the simulated galaxies has a circular velocity profile exceeding 20 km/s in the inner 1 kpc, implying that supernova feedback is key to solve the "too-big-to-fail" problem for Milky Way subhalos.

We report the discovery of 28 candidate high-velocity stars (HVSs) at heliocentric distances of less than 3 kpc, based on the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data Release 1. Our sample of HVS candidates covers a much broader color range than the equivalent ranges discussed in previous studies and comprises the first and largest sample of HVSs in the solar neighborhood. The sufficiently accurate observed and derived parameters for all candidates allow us to ascertain their nature as genuine HVSs, while a subset of 12 objects represents the most promising candidates. Our results also highlight the great potential of discovering statistically large numbers of HVSs of different spectral types in LAMOST survey data. This will ultimately enable us to achieve a better understanding of the nature of Galactic HVSs and their ejection mechanisms, and to constrain the structure of the Galaxy.

LAMOST has released more than two million spectra, which provide the opportunity to search for double-peaked narrow emission line (NEL) galaxies and AGNs. The double-peaked narrow-line profiles can be well modeled by two velocity components, respectively blueshifted and redshifted with respect to the systemic recession velocity. This paper presents 20 double-peaked NEL galaxies and AGNs found from LAMOST DR1 using a search method based on multi-gaussian fit of the narrow emission lines. Among them, 10 have already been published by other authors, either listed as genuine double-peaked NEL objects or as asymmetric NEL objects, the remaining 10 being first discoveries. We discuss some possible origins for double-peaked narrow-line features, as interaction between jet and narrow line regions, interaction with companion galaxies and black hole binaries. Spatially resolved optical imaging and/or follow-up observations in other spectral bands are needed to further discuss the physical mechanisms at work.

Based on SDSS and South Galactic Cap of u-band Sky Survey (SCUSS) early data, we use star counts method to estimate the Galactic structure parameters in an intermediate latitude with 10,180 main-sequence (MS) stars in absolute magnitude interval of $4 \leq M_r \leq 13$. We divide the absolute magnitude into five intervals:$4 \leq M_r < 5$, $5 \leq M_r < 6$, $6 \leq M_r < 8$, $8 \leq M_r < 10$, $10 \leq M_r \leq 13$, and estimate the Galactic structure parameters in each absolute magnitude interval to explore their possible variation with the absolute magnitude. Our study shows the parameters depend on absolute magnitude. For the thin disk, the intrinsic faint MS stars have large local space density and they tend to stay close to the Galactic plane. A plausible explanation is that faint MS stars with long lifetime experience long gravitational interaction time result in a short scaleheight. However, for the thick disk, the parameters show a complex trend with absolute magnitude, which may imply the complicated original of the thick disk. For the halo, the intrinsic faint MS stars have large local density and small axial ratio, which indicate a flattened inner halo and a more spherical outer halo.

We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in a fully cosmological Lambda cold dark matter context as part of high-resolution smoothed particle hydrodynamic simulations. We performed bulge-disc decompositions of the galaxies using artificial observations, in order to conduct a fair comparison to observations. We also dynamically decomposed the galaxies and compared the star formation histories of the bulges to those of the entire galaxies. These star formation histories showed that the bulges were primarily formed before z = 1 and during periods of rapid star formation. Both galaxies have large amounts of early star formation, which is likely related to the relatively high bulge-to-disc ratios also measured for them. Unlike almost all previous cosmological simulations, the realistically concentrated bulges of these galaxies do not lead to unphysically high rotational velocities, causing them to naturally lie along the observed Tully-Fisher relation.

We compare the central mass distribution of galaxies simulated with three different models of the interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to 10^4K, additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H_2) with shielding of atomic and molecular hydrogen, in addition to metal line cooling. In order to analyze the effect of these models, we follow the evolution of four field galaxies with V_peak < 120 km/s to a redshift of zero using high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological LCDM context. The spiral galaxies produced in simulations with either primordial cooling or H_2 physics have realistic, rising rotation curves. In contrast, the simulations with metal line cooling and otherwise similar feedback and star formation produced galaxies with the peaked rotation curves typical of most previous LCDM simulations of spiral galaxies. The less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling or H_2 are linked to changes in the angular momentum distribution of the baryons. These galaxies had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback-heated gas cools more slowly than when metal line cooling is included and so requires less energy to be expelled. When H_2 is included, the accompanying shielding produces large amounts of clumpy, cold gas where H_2 forms. Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material. (abridged)

We examine the statistical properties of the brightest group galaxies (BGGs) using a complete spectroscopic sample of groups/clusters of galaxies selected from the Data Release 7 of the Sloan Digital Sky Survey. We test whether BGGs and other bright members of groups are consistent with an ordered population among the total population of group galaxies. We find that the luminosity distributions of BGGs do not follow the predictions from the order statistics (OS). The average luminosities of BGGs are systematically brighter than OS predictions. On the other hand, by properly taking into account the brightening effect of the BGGs, the luminosity distributions of the second brightest galaxies are in excellent agreement with the expectations of OS. The brightening of BGGs relative to the OS expectation is consistent with a scenario that the BGGs on average have over-grown about 20 percent masses relative to the other member galaxies. The growth ($\Delta M$) is not stochastic but correlated with the magnitude gap ($G_{1,2}$) between the brightest and the second brightest galaxy. The growth ($\Delta M$) is larger for the groups having more prominent BGGs (larger $G_{1,2}$) and averagely contributes about 30 percent of the final $G_{1,2}$ of the groups of galaxies.

We use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk.

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift 3, using the "Eris" suite of cosmological "zoom-in" simulations. The reference run adopts a blastwave scheme for supernova feedback that produces galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. Synthetic spectra through the multiphase CGM produce interstellar absorption line strengths of Lya, CII, CIV, SiII, and SiIV as a function of galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high-redshift by Steidel et al. (2010). Only about one third of all the gas within R_vir is outflowing. The fraction of sightlines within one virial radius that intercept optically thick material is 27%, in agreement with recent observations by Rudie et al. (2012). Such optically thick absorption is shown to trace inflowing "cold" streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar, give origin to strong (log N > 13) CII absorption with a covering factor of 22% within R_vir and 10% within 2 R_vir. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large OVI halo, with a typical column density log N>14 and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. (2011). Our zoom-in simulations of this single system appear to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings. (Abridged)

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Ly\alpha absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

We present a detailed study of the metal-enriched circum-galactic medium of a massive galaxy at z=3 using "ErisMC", a new cosmological hydrodynamic "zoom-in" simulations of a disk galaxy with mass comparable to the Milky Way. The run adopts a blastwave scheme for supernova feedback that generates galactic outflows without explicit wind particles, a star formation recipe based on a high gas density threshold, and high temperature metal cooling. ErisMC's main progenitor at z=3 resembles a "Lyman break" galaxy of mass M_vir=2.4e11 M_sun, virial radius R_vir=48 kpc, and star formation rate 18 M_sun/yr, and its metal-enriched CGM extends to 200 (physical) kpc from its center. Approximately 41, 9, and 50 percent of all gas-phase metals at z=3 are locked in a hot (T> 3e5 K), warm (3e5 >T> 3e4 K), and cold (T< 3e4 K) medium, respectively. We identify three sources of heavy elements: 1) the main host, responsible for 60% of all the metals found within 3R_vir; 2) its satellite progenitors, responsible for 28% of all the metals within 3R_vir, and for only 5% of those beyond 3R_vir; and nearby dwarfs, which give origin to 12% of all the metals within 3R_vir and 95% of those beyond 3R_vir. Late (z<5) galactic "superwinds" account for only 9% of all the metals observed beyond 2R_vir, the bulk having been released at redshifts 5< z < 8 by early star formation and outflows. In the CGM, lower overdensities are typically enriched by `older', colder metals. Heavy elements are accreted onto Eris along filaments via low-metallicity cold inflows, and are ejected hot via galactic outflows at a few hundred km/s. The outflow mass-loading factor is of order unity for the main halo, but can exceed 10 for nearby dwarfs. We stress that our "zoom-in" simulation focuses on the CGM of a single massive system and cannot describe the enrichment history of the IGM as a whole. (abridged)

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<\alpha<67^\circ$ and $42^\circ<\delta<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

The disk galaxy UGC8802 has high neutral gas content and a flat profile of star formation rate compared to other disk galaxies with similar stellar mass. It also shows a steep metallicity gradient. We construct a chemical evolution model to explore its growth history by assuming its disk grows gradually from continuous gas infall, which is shaped by a free parameter -- the infall-peak time. By adopting the recently observed molecular surface density related star formation law, we show that a late infall-peak time can naturally explain the observed high neutral gas content, while an inside-out disk formation scenario can fairly reproduce the steep oxygen abundance gradient. Our results show that most of the observed features of UGC8802 can be well reproduced by simply `turning the knob' on gas inflow with one single parameter, which implies that the observed properties of gas-rich galaxies could also be modelled in a similar way.

Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust grain formation of H2, shielding of HI and H2, and photodissociation of H2 by Lyman-Werner radiation from nearby stellar populations. Because this model does not assume equilibrium abundances, it is particularly well suited for simulations that model low-metallicity environments, such as dwarf galaxies and the early Universe. We further introduce an explicit link between star formation and local H2 abundance. This link limits star formation to "star-forming regions," represented by areas with abundant H2. With this implementation, we determine the effect of H2 on star formation in a cosmological simulation of a dwarf galaxy. This simulation is the first cosmological simulation with non-equilibrium H2 abundances to be integrated to a redshift of zero or to include efficient SN feedback. We find that our simulations are consistent with the observed Tully-Fisher, global Kennicutt-Schmidt, and resolved Kennicutt-Schmidt relations. We find that the inclusion of shielding of both the atomic and molecular hydrogen and, to a lesser extent, the additional cooling from H2 at temperatures between 200 and 5000 K increases the amount of cold gas in the galaxies. The changes to the ISM result in an increased amount of cold, dense gas in the disk of the galaxy and the formation of a clumpier interstellar media (ISM). The explicit link between star formation and H2 and the clumpier ISM results in a bluer galaxy with a greater spatial distribution of star formation at a redshift of zero. (abridged)

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.

A clear prediction of the Cold Dark Matter model is the existence of cuspy dark matter halo density profiles on all mass scales. This is not in agreement with the observed rotation curves of spiral galaxies, challenging on small scales the otherwise successful CDM paradigm. In this work we employ high resolution cosmological hydro-dynamical simulations to study the effects of dissipative processes on the inner distribution of dark matter in Milky-Way like objects (M~1e12 Msun). Our simulations include supernova feedback, and the effects of the radiation pressure of massive stars before they explode as supernovae. The increased stellar feedback results in the expansion of the dark matter halo instead of contraction with respect to N-body simulations. Baryons are able to erase the dark matter cuspy distribution creating a flat, cored, dark matter density profile in the central several kpc of a massive Milky-Way like halo. The profile is well fit by a Burkert profile, with fitting parameters consistent with the observations. In addition, we obtain flat rotation curves as well as extended, exponential stellar disk profiles. While the stellar disk we obtain is still partially too thick to resemble the MW thin disk, this pilot study shows that there is enough energy available in the baryonic component to alter the dark matter distribution even in massive disc galaxies, providing a possible solution to the long standing problem of cusps vs. cores.

We investigate the formation of the first massive black holes in high redshift galaxies, with the goal of providing insights to which galaxies do or do not host massive black holes. We adopt a novel approach to forming seed black holes in galaxy halos in cosmological SPH+N-body simulations. The formation of massive black hole seeds is dictated directly by the local gas density, temperature, and metallicity, and motivated by physical models of massive black hole formation. We explore seed black hole populations as a function of halo mass and redshift, and examine how varying the efficiency of massive black hole seed formation affects the relationship between black holes and their hosts. Seed black holes tend to form in halos with mass between 10^7 and 10^9 Msun, and the formation rate is suppressed around z = 5 due to the diffusion of metals throughout the intergalactic medium. We find that the time of massive black hole formation and the occupation fraction of black holes are a function of the host halo mass. By z = 5, halos with mass M_halo > 3 x 10^9 Msun host massive black holes regardless of the efficiency of seed formation, while the occupation fraction for smaller halos increases with black hole formation efficiency. Our simulations explain why massive black holes are found in some bulgeless and dwarf galaxies, but we also predict that their occurrence becomes rarer and rarer in low-mass systems.

We perform SPH+N-body cosmological simulations of massive disk galaxies, including a formalism for black hole seed formation and growth, and find that satellite galaxies containing supermassive black hole seeds are often stripped as they merge with the primary galaxy. These events naturally create a population of "wandering" black holes that are the remnants of stripped satellite cores; galaxies like the Milky Way may host 5 -- 15 of these objects within their halos. The satellites that harbor black hole seeds are comparable to Local Group dwarf galaxies such as the Small and Large Magellanic Clouds; these galaxies are promising candidates to host nearby intermediate mass black holes. Provided that these wandering black holes retain a gaseous accretion disk from their host dwarf galaxy, they give a physical explanation for the origin and observed properties of some recently discovered off-nuclear ultraluminous X-ray sources such as HLX-1.

We introduce a simple model to explore the star formation histories of disk galaxies. We assume that the disk origins and grows by continuous gas infall. The gas infall rate is parametrized by the Gaussian formula with one free parameter: infall-peak time $t_p$. The Kennicutt star formation law is adopted to describe how much cold gas turns into stars. The gas outflow process is also considered in our model. We find that, at given galactic stellar mass $M_*$, model adopting late infall-peak time $t_p$ results in blue colors, low metallicity, high specific star formation rate and high gas fraction, while gas outflow rate mainly influences the gas-phase metallicity and star formation efficiency mainly influences the gas fraction. Motivated by the local observed scaling relations, we construct a mass-dependent model by assuming low mass galaxy has later infall-peak time $t_p$ and larger gas outflow rate than massive systems. It is shown that this model can be in agreement with not only the local observations, but also the observed correlations between specific star formation rate and galactic stellar mass $SFR/M_* \sim M_*$ at intermediate redshift $z<1$. Comparison between the Gaussian-infall model and exponential-infall model is also presented. It shows that the exponential-infall model predicts higher star formation rate at early stage and lower star formation rate later than that of Gaussian-infall. Our results suggest that the Gaussian infall rate may be more reasonable to describe the gas cooling process than the exponential infall rate, especially for low-mass systems.

We examine the distribution of axis ratios of a large sample of disk galaxies hosting type 2 AGNs selected from the Sloan Digital Sky Survey and compare it with a well-defined control sample of non-active galaxies. We find them significantly different, where the type 2 AGNs show both an excess of edge-on objects and deficit of round objects. This systematical bias can not be explained by a nuclear obscurer oriented randomly with respect to the stellar disk. However, a nuclear obscurer coplanar with the stellar disk also does not fit the data very well. By assuming that the nuclear obscurer having an opening angle of ~60 degree, we find the observed axis ratio distribution can be nicely reproduced by a mean tilt angle of ~30 degree between the nuclear obscurer and the stellar disk.

Fully cosmological, high resolution N-Body + SPH simulations are used to investigate the chemical abundance trends of stars in simulated stellar halos as a function of their origin. These simulations employ a physically motivated supernova feedback recipe, as well as metal enrichment, metal cooling and metal diffusion. As presented in an earlier paper, the simulated galaxies in this study are surrounded by stellar halos whose inner regions contain both stars accreted from satellite galaxies and stars formed in situ in the central regions of the main galaxies and later displaced by mergers into their inner halos. The abundance patterns ([Fe/H] and [O/Fe]) of halo stars located within 10 kpc of a solar-like observer are analyzed. We find that for galaxies which have not experienced a recent major merger, in situ stars at the high [Fe/H] end of the metallicity distribution function are more [alpha/Fe]-rich than accreted stars at similar [Fe/H]. This dichotomy in the [O/Fe] of halo stars at a given [Fe/H] results from the different potential wells within which in situ and accreted halo stars form. These results qualitatively match recent observations of local Milky Way halo stars. It may thus be possible for observers to uncover the relative contribution of different physical processes to the formation of stellar halos by observing such trends in the halo populations of the Milky Way, and other local L* galaxies.

We report the discovery of eight new quasars in one extragalactic field (five degree centered at RA=$08^h58^m08.2^s$, Dec=$01^o32'29.7''$) with the LAMOST commissioning observations on December 18, 2009. These quasars, with $i$ magnitudes from 16.44 to 19.34 and redshifts from 0.898 to 2.773, were not identified in the SDSS spectroscopic survey, though six of them with redshifts less than 2.5 were selected as quasar targets in SDSS. Except one source without near-IR $Y$-band data, seven of these eight new quasars meet a newly proposed quasar selection criterion involving both near-IR and optical colors. Two of them were found in the 'redshift desert' for quasars ($z$ from 2.2 to 3), indicating that the new criterion is efficient for recovering the missing quasars with similar optical colors as stars. Although LAMOST met some problems during the commissioning observations, we were still able to identify other 38 known SDSS quasars in this field, with $i$ magnitudes from 16.24 to 19.10 and redshifts from 0.297 to 4.512. Our identifications imply that a substantial fraction of quasars may be missing in the previous quasar surveys. The implication of our results to the future LAMOST quasar survey is discussed.

A study of the IGM metal enrichment using a series of SPH simulations is presented, employing metal cooling and turbulent diffusion of metals and thermal energy. An adiabatic feedback mechanism was adopted where gas cooling was prevented to generate galactic winds without explicit wind particles. The simulations produced a cosmic star formation history (SFH) that is broadly consistent with observations until z $\sim$ 0.5, and a steady evolution of the universal neutral hydrogen fraction ($\Omega_{\rm H I}$). At z=0, about 40% of the baryons are in the warm-hot intergalactic medium (WHIM), but most metals (80%-90%) are locked in stars. At higher z the proportion of metals in the IGM is higher due to more efficient loss from galaxies. The IGM metals primarily reside in the WHIM throughout cosmic history. The metallicity evolution of the gas inside galaxies is broadly consistent with observations, but the diffuse IGM is under enriched at z $\sim$ 2.5. Galactic winds most efficiently enrich the IGM for halos in the intermediate mass range $10^{10}$M$_{\sun}$ - $10^{11}$ M$_{\sun}$. At the low mass end gas is prevented from accreting onto halos and has very low metallicities. At the high mass end, the fraction of halo baryons escaped as winds declines along with the decline of stellar mass fraction of the galaxies. This is likely because of the decrease in star formation activity and in wind escape efficiency. Metals enhance cooling which allows WHIM gas to cool onto galaxies and increases star formation. Metal diffusion allows winds to mix prior to escape, decreasing the IGM metal content in favour of gas within galactic halos and star forming gas. Diffusion significantly increases the amount of gas with low metallicities and changes the density-metallicity relation.

The redshift range from 2.2 to 3, is known as the 'redshift desert' of quasars because quasars with redshift in this range have similar optical colors as normal stars and are thus difficult to be found in optical sky surveys. A quasar candidate, SDSS J085543.40-001517.7, which was selected by a recently proposed criterion involving near-IR $Y-K$ and optical $g-z$ colors, was identified spectroscopically as a new quasar with redshift of 2.427 by the LAMOST commissioning observation in December 2009 and confirmed by the observation made with the NAOC/Xinglong 2.16m telescope in March 2010. This quasar was not targeted in the SDSS spectroscopic survey because it locates in the stellar locus of the optical color-color diagrams, while it is clearly separated from stars in the $Y-K$ vs. $g-z$ diagram. Comparing with other SDSS quasars we found this new quasar with $i$ magnitude of 16.44 is apparently the brightest one in the redshift range from 2.3 to 2.7. From the spectral properties we derived its central black hole mass as $(1.4\sim3.9) \times 10^{10} M_\odot$ and the bolometric luminosity as $3.7\times 10^{48}$ \ergs, which indicates that this new quasar is intrinsically very bright and belongs to the most luminous quasars in the universe. Our identification supports that quasars in the redshift desert can be found by the quasar selection criterion involving the near-IR colors. More missing quasars are expected to be recovered by the future LAMOST spectroscopic surveys, which is important to the study of the cosmological evolution of quasars at redshift higher than 2.2.

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy initial mass function at high metallicity and high redshift. We examine whether this hypothesis has any testable observational consequences. First we determine, using a set of hydrodynamic galaxy formation simulations, that the CMB temperature floor would have influenced the majority of stars formed at redshifts between z=3 and 6, and probably even to higher redshift. Five signatures of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current simulations; and (5) an enhancement in the abundances of alpha elements such as O and Mg at metallicities -2 < [Fe/H] < -0.5. Observations are not presently able to either confirm or rule out the presence of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular clusters could provide an explanation for the observed bimodality of their metallicity distribution.

The Tully-Fisher relation of spiral galaxies shows notable dependence on morphological types, with earlier type spirals having systematically lower luminosity at fixed maximum rotation velocity $V_{max}$. This decrement of luminosity is more significant in shorter wavelengths. By modeling the rotation curve and stellar population of different morphological type spiral galaxies in combination, we find the $V_{max}$ of spiral galaxies is weakly dependent on the morphological type, whereas the difference of the stellar population originating from the bulge disk composition effect mainly account for the morphological type dependence of the Tully-Fisher relation.

The formation of brown dwarfs (BDs) due to the fragmentation of proto-stellar disks undergoing pairwise encounters was investigated. High resolution allowed the use of realistic initial disk models where both the vertical structure and the local Jeans mass were resolved. The results show that objects with masses ranging from giant planets to low mass stars can form during such encounters from initially stable disks. The parameter space of initial spin-orbit orientations and the azimuthal angles for each disk was explored. An upper limit on the initial Toomre Q value of ~2 was found for fragmentation to occur. Depending on the initial configuration, shocks, tidal-tail structures and mass inflows were responsible for the condensation of disk gas. Retrograde disks were generally more likely to fragment. When the interaction timescale was significantly shorter than the disks' dynamical timescales, the proto-stellar disks tended to be truncated without forming objects. The newly-formed objects had masses ranging from 0.9 to 127 Jupiter masses, with the majority in the BD regime. They often resided in star-BD multiples and in some cases also formed hierarchical orbiting systems. Most of them had large angular momenta and highly flattened, disk-like shapes. The objects had radii ranging from 0.1 to 10 AU. The disk gas was assumed to be locally isothermal, appropriate for the short cooling times in extended proto-stellar disks, but not for condensed objects. An additional case with explicit cooling that reduced to zero for optically thick gas was simulated to test the extremes of cooling effectiveness and it was still possible to form objects in this case. Detailed radiative transfer is expected to lengthen the internal evolution timescale for these objects, but not to alter our basic results.

We study the chemical evolution of the disks of the Milky Way (MW) and of Andromeda (M31), in order to reveal common points and differences between the two major galaxies of the Local group. We use a large set of observational data for M31, including recent observations of the Star Formation Rate (SFR) and gas profiles, as well as stellar metallicity distributions along its disk. We show that, when expressed in terms of the corresponding disk scale lengths, the observed radial profiles of MW and M31 exhibit interesting similarities, suggesting the possibility of a description within a common framework. We find that the profiles of stars, gas fraction and metallicity of the two galaxies, as well as most of their global properties, are well described by our model, provided the star formation efficiency in M31 disk is twice as large as in the MW. We show that the star formation rate profile of M31 cannot be fitted with any form of the Kennicutt-Schmidt law (KS Law) for star formation. We attribute those discrepancies to the fact that M31 has undergone a more active star formation history, even in the recent past, as suggested by observations of a "head-on" collision with the neighboring M32 galaxy about 200 Myr ago. The MW has most probably undergone a quiescent secular evolution, making possible a fairly successful description with a simple model. If M31 is more typical of spiral galaxies, as recently suggested by Hammer et al. (2007), more complex models, involving galaxy interactions, will be required for the description of spirals.

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 investigate the radial color gradients of galactic disks using a sample of about 20,000 face-on spiral galaxies selected from the fourth data release of the Sloan Digital Sky Survey (SDSS-DR4). We combine galaxies with similar concentration, size and luminosity to construct composite galaxies, and then measure their color profiles by stacking the azimuthally averaged radial color profiles of all the member galaxies. Except for the smallest galaxies (R_{50}<3 kpc), almost all galaxies show negative disk color gradients with mean g-r gradient G_{gr}=-0.006 mag kpc^{-1} and r-z gradient G_{rz}=-0.018 mag kpc^{-1}. The disk color gradients are independent of the morphological types of galaxies and strongly dependent on the disk surface brightness \mu_{d}, with lower surface brightness galactic disks having steeper color gradients. We quantify the intrinsic correlation between color gradients and surface brightness as G_{gr}=-0.011\mu_{d}+0.233 and G_{rz}=-0.015\mu_{d}+0.324. These quantified correlations provide tight observational constraints on the formation and evolution models of spiral galaxies.

We use the ROSAT All-Sky Survey to study the X-ray properties of a sample of 625 groups and clusters of galaxies selected from the Sloan Digital Sky Survey. We stack clusters with similar velocity dispersions and investigate whether their average X-ray luminosities and surface brightness profiles vary with the radio activity level of their central galaxies. We find that at a given value of $\sigma$, clusters with a central radio AGN have more concentrated X-ray surface brightness profiles, larger central galaxy masses, and higher X-ray luminosities than clusters with radio-quiet central galaxies. The enhancement in X-ray luminosity is more than a factor of two, is detected with better than 6$\sigma$ significance, and cannot be explained by X-ray emission from the radio AGN itself. This difference is largely due to a subpopulation of radio-quiet, high velocity dispersion clusters with low mass central galaxies. These clusters are underluminous at X-ray wavelengths when compared to otherwise similar clusters where the central galaxy is radio-loud, more massive, or both.

Using a samples of 61506 spiral galaxies selected from the SDSS DR2, we examine the luminosity function (LF) of spiral galaxies with different inclination angles. We find that the characteristic luminosity of the LF, $L^*$, decreases with increasing inclination, while the faint-end slope, $\alpha$, depends only weakly on it. The inclination-dependence of the LF is consistent with that expected from a simple model where the optical depth is proportional to the cosine of the inclination angle, and we use a likelihood method to recover both the coefficient in front of the cosine, $\gamma$, and the LF for galaxies viewed face-on. The value of $\gamma$ is quite independent of galaxy luminosity in a given band, and the values of $\gamma$ obtained in this way for the 5 SDSS bands give an extinction curve which is a power law of wavelength ($\tau\propto\lambda^{-n}$), with a power index $n=0.96\pm0.04$. Using the dust extinction for galaxies obtained by Kauffmann et al. (2003), we derive an `extinction-corrected' luminosity function for spiral galaxies. Dust extinction makes $M^*$ dimmer by about 0.5 magnitudes in the $z$-band, and about 1.2 magnitudes in the $u$- band. Since our analysis is based on a sample where selection effects are well under control, the dimming of edge-on galaxies relative to face-on galaxies is best explained by assuming that galaxy disks are optically thick in dust absorptions.

The formation of brown dwarfs via encounters between proto-stars has been confirmed with high-resolution numerical simulations with a restricted treatment of the thermal conditions. The new results indicate that young brown dwarfs (BDs) formed this way are disk-like and often reside in multiple systems. The newly-formed proto-BDs disks are up to 18 AU in size and spin rapidly making small-scale bipolar outflows, fragmentation and the possible formation of planetary companions likely as have recently been observed for BDs. The object masses range from 2 to 73 Jupiter masses, distributed in a manner consistent with the observed sub-stellar initial mass function. The simulations usually form multiple BDs on eccentric orbits about a star. One such system was hierarchical, a BD binary in orbit around a star, which may explain recently observed hierarchical systems. One third of the BDs were unbound after a few thousand years and interactions among orbiting BDs may eject more or add to the number of binaries. Improvements over prior work include resolution down to a Jupiter mass, self-consistent models of the vertical structure of the initial disks and careful attention to avoid artificial fragmentation.

We use a large sample of galaxies drawn from the Sloan Digital Sky Survey (SDSS) and Two Micro All Sky Survey (2MASS) to present Color-Magnitude Relations (CMRs) for late-type galaxies in both optical and optical-infrared bands. A sample from SDSS Data Release 4 (DR4) is selected to investigate the optical properties. Optical-infrared colors are estimated from a position matched sample of DR4 and 2MASS, in which the photometric aperture mismatch between these two surveys is carefully corrected. It is shown that, after correcting the dust attenuation, the optical colors for faint galaxies (i.e. $M_r > -21$) have very weak correlation with the luminosity, while the optical colors for bright galaxies (i.e. $M_r < -21$) are redder for more luminous galaxies. All (optical, optical-infrared and infrared) colors show similar but stronger correlations with stellar mass than with absolute magnitude. The optical colors correlate more strongly with stellar mass surface density than with stellar mass, while optical-infrared and infrared colors show stronger correlations with stellar mass. By comparing the observed colors of our sample galaxies with the colors predicted by stellar population synthesis model, we find that massive late-type galaxies have older and higher metallicity stellar population than less massive galaxies. This suggests that CMRs for late-type galaxies are trends defined by the combination of stellar mean age and metallicity. Moreover, our results suggest that the stellar mean metallicity of late-type galaxy is mainly determined by its stellar mass, while the star formation history is mainly regulated by the stellar mass surface density.

We use the ROSAT All Sky Survey (RASS) to study the soft X-ray properties of a homogeneous sample of 46,420 quasars selected from the third data release of the Sloan Digital Sky Survey (SDSS). Optical luminosities, both at rest-frame 2500\AA ($L_{2500}$) and in [OIII] ($L_{[\rm{OIII}]}$) span more than three orders of magnitude, while redshifts range over $0.1<z<5.4$. We detect 3366 quasars directly in the observed 0.1--2.4 keV band. Sub-samples of radio-loud and radio-quiet objects (RLQs and RQQs) are obtained by cross-matching with the FIRST catalogue. We study the distribution of X-ray luminosity as a function of optical luminosity, redshift and radio power using both individual detections and stacks of complete sets of similar quasars. At every optical luminosity and redshift $\log L_{2\kev}$ is, to a good approximation, normally distributed with dispersion $\sim 0.40$, at least brightwards of the median X-ray luminosity. This median X-ray luminosity of quasars is a power law of optical luminosity with index $\sim 0.53$ for $L_{2500}$ and $\sim 0.30$ for $L_{[\rm{OIII}]}$. RLQs are systematically brighter than RQQs by about a factor of 2 at given optical luminosity. The zero-points of these relations increase systematically with redshift, possibly in different ways for RLQs and RQQs. Evolution is particularly strong at low redshift and if the optical luminosity is characterised by $L_{[\rm{OIII}]}$. At low redshift and at given $L_{[\rm{OIII}]}$ the soft X-ray emission from type II AGN is more than 100 times weaker than that from type I AGN.

We have examined some basic properties of damped Ly$\alpha$ systems(DLAs) by semi-analytic model. We assume that DLA hosts are disk galaxies whose mass function is generated by Press-Schechter formulism at redshift 3. Star formation and chemical evolution undergo in the disc. We select modelled DLAs according to their observational criterion by Monte Carlo simulation using random line of sights and disk inclinations. The DLA ages are set to be 1 to 3 Gyr. By best-fitting the predicted metallicity distribution to the observed ones, we get the effective yield for DLAs about $0.25Z_{\odot}$. On the basis of this constrain, we further compared our model predictions with observations at redshift 3 in the following items: number density; gas content; HI frequency distribution; star formation rate density; relationship between metallicity and HI column density. We found that the predicted number density at redshift 3 agree well with the observed value, but the gas content $\Omega_{DLA}$ is about 3 times larger than observed since our model predicts more DLA systems with higher column density. The frequency distribution at higher HI column density is quite consistent with observation while some difference exists at lower HI end. The predicted star formation rate density contributed by DLAs is consistent with the most recent observations. Meanwhile, the connection between DLAs and Lyman Break galaxies(LBGs) is discussed by comparing their UV luminosity functions which shows that the DLAs host galaxies are much fainter than LBGs. However, there is a discrepancy between model prediction and observation in the correlation between metallicity and HI column density for DLAs. Further investigations are needed for the star formation mode at high redshift environments.

Based on the disk galaxy formation theory within the framework of standard LCDM hierarchical picture, we selected modelled DLAs, according to their observational criterion, by Monte Carlo simulation with the random inclinations being considered, to examine their observed properties. By best-fitting the predicted metallicity distribution to the observed ones, we get the effective yield for DLAs about 0.25Z_sun, which is comparable to those for SMC and LMC. And the predicted distribution is the same as that of observation at the significant level higher than 60%. The predicted column density distribution of modelled DLAs is compared with the observed ones with the corresponding number density, gas content being discussed. We found that the predicted number density n(z) at redshift 3 agree well with the observed value, but the gas content \Omega_DLA is about 3 times larger than observed since our model predicts more DLA systems with higher column density. It should be noted that the predicted star formation rate density contributed by DLAs is consistent with the most recent observations if the star formation timescale in DLAs is assumed to be 1 to 3 Gyr. Meanwhile, the connection between DLAs and LBGs is discussed by comparing their UV luminosity functions which shows that the DLAs host galaxies are much fainter than LBGs. We also predict that only few percent of DLAs can host LBGs which is also consistent with current observations. However, there is a discrepancy between model prediction and observation in the correlation between metallicity and HI column density for DLAs. We suggest that this could result from either the inadequacy of Schmidt-type star formation law at high redshift, the diversities of DLA populations, or the model limitations.

abridged: We use a complete sample of about 140,000 galaxies from the Sloan Digital Sky Survey (SDSS) to study the size distribution of galaxies and its dependence on their luminosity, stellar mass, and morphological type. The large SDSS database provides statistics of unprecedented accuracy. For each type of galaxy, the size distribution at given luminosity (or stellar mass) is well described by a log-normal function, characterized by its median $\bar{R}$ and dispersion $\sigma_{\ln R}$. For late-type galaxies, there is a characteristic luminosity at $M_{r,0}\sim -20.5$ (assuming $h=0.7$) corresponding to a stellar mass $M_0\sim 10^{10.6}\Msun$. Galaxies more massive than $M_0$ have $\bar{R}\propto M^{0.4}$ and $\sigma_{\ln R}\sim 0.3$, while less massive galaxies have $\bar{R}\propto M^{0.15}$ and $\sigma_{\ln R}\sim 0.5$. For early-type galaxies, the $\bar{R}$ - $M$ relation is significantly steeper, $\bar{R}\propto M^{0.55}$, but the $\sigma_{\ln R}$ - $M$ relation is similar to that of late-type galaxies. Faint red galaxies have sizes quite independent of their luminosities.

Current theory of disk galaxy formation is used to study fundamental-plane (FP) type of relations for disk galaxies. We examine how the changes in model parameters affect these relations and explore the possibility of using such relations to constrain theoretical models. The distribution of galaxy disks in the space of their fundamental properties are predicted to be concentrated in a plane, with the Tully-Fisher (TF) relation (a relation between luminosity $L$ and maximum rotation velocity $V_m$) being an almost edge-on view. Using rotation velocities at larger radii generally leads to larger TF scatter. In searching for a third parameter, we find that both the disk scale-length $R_d$ (or surface brightness) and the rotation-curve shape are correlated with the TF scatter. The FP relation in the $(\Log L, \Log V_m, \Log R_d)$-space obtained from the theory is $L\propto R_d^{\alpha'} V_m^{\beta'}$, with ${\alpha'}\sim 0.50$ and ${\beta'}\sim 2.60$, consistent with the preliminary result we obtain from observational data. Among the model parameters we probe, variation in any of them can generate significant scatter in the TF relation, but the effects of the spin parameter and halo concentration can be reduced significantly by introducing $R_d$ while the scatter caused by varying $m_d$ (the ratio between disk mass and halo mass) is most effectively reduced by introducing the parameters which describes the rotation-curve shape. The TF and FP relations combined should therefore provide useful constraints on models of galaxy formation.