We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$.

We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galaxy (LRG) and a "bright galaxy" (BG) sample ($r<21$). We have found 286 new high probability candidates, 133 from the LRG sample and 153 from the BG sample. We have then ranked these candidates based on a value that combines the CNN likelihood to be a lens and the human score resulting from visual inspection (P-value) and we present here the highest 82 ranked candidates with P-values $\ge 0.5$. All these high-quality candidates have obvious arc or point-like features around the central red defector. Moreover, we define the best 26 objects, all with scores P-values $\ge 0.7$ as a "golden sample" of candidates. This sample is expected to contain very few false positives and thus it is suitable for follow-up observations. The new lens candidates come partially from the the more extended footprint adopted here with respect to the previous analyses, partially from a larger predictive sample (also including the BG sample). These results show that machine learning tools are very promising to find strong lenses in large surveys and more candidates that can be found by enlarging the predictive samples beyond the standard assumption of LRGs. In the future, we plan to apply our CNN to the data from next-generation surveys such as the Large Synoptic Survey Telescope, Euclid, and the Chinese Space Station Optical Survey.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

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

Galaxy properties are known to correlate most tightly with the galaxy effective stellar velocity dispersion $\sigma_{\rm e}$. Here we look for {\em additional} trends at fixed $\sigma_{\rm e}$ using 1339 galaxies ($M_\ast \gtrsim 6\times10^9$ M$_\odot$) with different morphologies in the MaNGA (DR14) sample with integral-field spectroscopy data. We focus on the gradients ($\gamma_{\rm rms} \equiv \sigma(R_{\rm e}/4)/\sigma_{\rm e}$) of the stellar root-mean-square velocity ($V_{\rm rms} \equiv \sqrt{V^2 + \sigma^2}$), which we show traces the total mass density gradient $\gamma_{\rm tot}$ derived from dynamical models and, more weakly, the bulge fraction. We confirm that $\gamma_{\rm rms}$ increases with $\sigma_{\rm e}$, age and metallicity. We additionally find that these correlations still exist at fixed $\sigma_{\rm e}$, where galaxies with larger $\gamma_{\rm rms}$ are found to be older and more metal-rich. It means that mass density gradients contain information of the stellar population which is not fully accounted for by $\sigma_{\rm e}$. This result puts an extra constraint on our understanding of galaxy quenching. We compare our results with galaxies in the IllustrisTNG hydrodynamical simulations and find that, at fixed $\sigma_{\rm e}$, similar trends exist with age, the bulge fraction, and the total mass density slope but, unlike observations, no correlation with metallicity can be detected in the simulations.

The streaming instability (SI) is a mechanism to aerodynamically concentrate solids in protoplanetary disks and facilitate the formation of planetesimals. Recent numerical modeling efforts have demonstrated the increasing complexity of the initial mass distribution of planetesimals. To better constrain this distribution, we conduct SI simulations including the self-gravity with hitherto the highest resolution. To subsequently identify all of the self-bound clumps, we develop a new clump-finding tool, PLanetesimal ANalyzer (\texttt{PLAN}). We then apply a maximum likelihood estimator to fit a suite of parameterized models with different levels of complexity to the simulated mass distribution. To determine which models are best-fitting and statistically robust, we apply three model selection criteria with different complexity penalties. We find that the initial mass distribution of clumps is not universal regarding both the functional forms and parameter values. Our model selection criteria prefer models different from those previously considered in the literature. Fits to multi-segment power law models break to a steeper distribution above masses close to 100 km collapsed planetesimals, similar to observed Kuiper Belt size distributions. We find evidence for a turnover in the low mass end of the planetesimal mass distribution in our high resolution run. Such a turnover is expected for gravitational collapse, but had not previously been reported.

We develop an analytic mass model for lensing galaxies, based on a broken power-law (BPL) density profile, which is a power-law profile with a mass deficit or surplus in the central region. Under the assumption of an elliptically symmetric surface mass distribution, the deflection angle and magnification can be evaluated analytically for this new model. We compute the theoretical prediction for various quantities, including the volume and surface mass density profiles of the galaxies, and the aperture and luminosity-weighted line-of-sight velocity dispersions, and compare them to those measured from the Illustris simulation. We find an excellent agreement between our model prediction and the simulation, which validates our modeling. The high efficiency and accuracy of our model manifests itself as a promising tool for studying properties of galaxies with strong lensing.

Kilonovae are optical flashes produced in the aftermath of neutron star-neutron star mergers (NNMs) or neutron star-black hole mergers (NBMs). In this work, we use the Millennium Simulation, combined with a semi-analytic galaxy formation model--GABE (Galaxy Assembly with Binary Evolution) which adopts binary stellar population synthesis models, to explore the cosmic event rate of kilonovae, and the properties of their host galaxies in a cosmological context. We find that model with supernova kick velocity of 0 km/s fits the observation best, in agreement with the exception of some formation channels of binary neutron star. This indicates that NNMs prefer to originate from binary systems with low kick velocities. With V$_{\rm kick}$=0 km/s, the cosmic event rate of NNMs and NBMs at z=0 are 283 Gpc$^{-3}$yr$^{-1}$ and 91 Gpc$^{-3}$yr$^{-1}$, respectively, marginally consistent with the constraint from LVC GWTC-1. For Milky Way-mass galaxies, we predict the NNM rate is $25.7^{+59.6}_{-7.1}$ Myr$^{-1}$, which is also in good agreement with the observed properties of binary neutron stars in the Milky Way. Taking all the NNMs into account in the history of Milky Way-mass galaxies, we find that the averaged r-process elements yield with A>79 in a NNM and NBM event should be 0.01 M$_{\odot}$ to be consistent with observation. We conclude that NGC 4993, the host galaxy of GW170817, is a typical host galaxy for NNMs. However, generally NNMs and NBMs tend to reside in young, blue, star-forming, late-type galaxies, with stellar mass and gaseous metallicity distribution peaking at $10^{10.65}$ M$_{\odot}$ and 8.72-8.85, respectively. By studying kilonovae host galaxies in the cosmological background, it is promising to constrain model details better when we have more events in the forthcoming future. (abridged)

We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse. Turbulence impedes the formation of dense particle clumps by acting as an effective turbulent diffusivity, but it can also promote planetesimal formation by concentrating solids, for example in zonal flows. We quantify the effect of turbulent diffusivity using numerical simulations of the streaming instability in small local domains, forced with velocity perturbations that establish approximately Kolmogorov-like fluid turbulence. We find that planetesimal formation is strongly suppressed by turbulence once velocity fluctuations exceed a threshold value of $\delta v^2 \simeq 10^{-3.5} - 10^{-3} c_s^2$. Turbulence whose strength is just below the threshold reduces the rate of solid material being converted into bound clumps. The main effect of turbulence is to thicken the mid-plane solid layer. Our results are thus consistent with planetesimal formation requiring a mid-plane solid-to-gas ratio $\epsilon \gtrsim 0.5$. We describe a method for tracking bound clumps in our simulations, and use this method to construct a mass function of planetesimals that is measured for each clump shortly after its collapse. Adopting this definition of the initial mass, instead of measuring masses at a fixed time, reduces planetesimal masses by a factor of three. For models in which planetesimals form, we show that the initial mass function is well-described by a broken power law, whose parameters are robust to the inclusion and strength of imposed turbulence. Turbulence in protoplanetary disks is likely to substantially exceed the threshold for planetesimal formation at radii where temperatures $T \gtrsim 10^3 \ {\rm K}$ lead to thermal ionization. Fully in situ planetesimal and planet formation may therefore not be viable for the closest-in exoplanets.

A new approach is presented to compute the seismic normal modes of a fully heterogeneous, rotating planet. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method (FMM). To solve the resulting quadratic eigenvalue problem (QEP), the approach utilizes extended Lanczos vectors forming a subspace computed from a non-rotating planet -- with the shape of boundaries of a rotating planet and accounting for the centrifugal potential -- to reduce the dimension of the original problem significantly. The subspace is guaranteed to be contained in the space of functions to which the seismic normal modes belong. The reduced system can further be solved with a standard eigensolver. The computational accuracy is illustrated using all the modes with relative small meshes and also tested against standard perturbation calculations relative to a standard Earth model. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Mars models studying the effects of heterogeneity on a large range of scales.

Recent research shows that the spectral lag is closely related to the spectral evolution in GRBs. In this paper, we study the spectral lag for a radiating jet shell with a high-energy cut-off radiation spectrum. For the jet shell with a cut-off power-law spectrum, the spectral lag monotonically increases with the photon energy and levels off at a certain photon energy. It is the same for the jet shell with a Band cut-off spectrum (Bandcut). However, a turn-over from the positive lags to negative lags appears in the high-energy range for the jet shell with Bandcut, which is very similar to that observed in GRB~160625B. The dependence of the spectral lags on the spectral shape/evolution are studied in details. In addition, the spectral lag behavior observed in GRB~160625B is naturally reproduced based on our theoretical outcome.

The streaming instability concentrates solid particles in protoplanetary disks, leading to gravitational collapse into planetesimals. Despite its key role in producing particle clumping and determining critical length scales in the instability's linear regime, the influence of the disk's radial pressure gradient on planetesimal properties has not been examined in detail. Here, we use streaming instability simulations that include particle self-gravity to study how the planetesimal initial mass function depends on the radial pressure gradient. Fitting our results to a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, we find a single value $p \approx 1.6$ describes simulations in which the pressure gradient varies by $\gtrsim 2$. An exponentially truncated power-law provides a significantly better fit, with a low mass slope of $p^\prime \approx 1.3$ that weakly depends on the pressure gradient. The characteristic truncation mass is found to be $\sim M_G = 4 \pi^5 G^2 \Sigma_p^3 / \Omega^4$. We exclude the cubic dependence of the characteristic mass with pressure gradient suggested by linear considerations, finding instead a linear scaling. These results strengthen the case for a streaming-derived initial mass function that depends at most weakly on the aerodynamic properties of the disk and participating solids. A simulation initialized with zero pressure gradient---which is {\em not} subject to the streaming instability---also yields a top-heavy mass function but with modest evidence for a different shape. We discuss the consistency of the theoretically predicted mass function with observations of Kuiper Belt planetesimals, and describe implications for models of early stage planet formation..

The Chinese Space Station Optical Survey (CSS-OS) is a planned full sky survey operated by the Chinese Space Station Telescope (CSST). It can simultaneously perform the photometric imaging and spectroscopic slitless surveys, and will probe weak and strong gravitational lensing, galaxy clustering, individual galaxies and galaxy clusters, active galactic nucleus (AGNs), and so on. It aims to explore the properties of dark matter and dark energy and other important cosmological problems. In this work, we focus on two main CSS-OS scientific goals, i.e. the weak gravitational lensing (WL) and galaxy clustering surveys. We generate the mock CSS-OS data based on the observational COSMOS and zCOSMOS catalogs. We investigate the constraints on the cosmological parameters from the CSS-OS using the Markov Chain Monte Carlo (MCMC) method. The intrinsic alignments, galaxy bias, velocity dispersion, and systematics from instrumental effects in the CSST WL and galaxy clustering surveys are also included, and their impacts on the constraint results are discussed. We find that the CSS-OS can improve the constraints on the cosmological parameters by a factor of a few (even one order of magnitude in the optimistic case), compared to the current WL and galaxy clustering surveys. The constraints can be further enhanced when performing joint analysis with the WL, galaxy clustering, and galaxy-galaxy lensing data. Therefore, the CSS-OS is expected to be a powerful survey for exploring the Universe. Since some assumptions may be still optimistic and simple, it is possible that the results from the real survey could be worse. We will study these issues in details with the help of simulations in the future.

We study the sudden optical and ultraviolet (UV) brightening of 1ES 1927+654, which until now was known as a narrow-line active galactic nucleus (AGN). 1ES 1927+654 was part of the small and peculiar class of "true Type-2" AGN, which lack broad emission lines and line-of-sight obscuration. Our high-cadence spectroscopic monitoring captures the appearance of a blue, featureless continuum, followed several weeks later by the appearance of broad Balmer emission lines. This timescale is generally consistent with the expected light travel time between the central engine and the broad-line emission region in (persistent) broad-line AGN. Hubble Space Telescope spectroscopy reveals no evidence for broad UV emission lines (e.g., CIV1549, CIII]1909, MgII2798), probably owing to dust in the broad-line emission region. To the best of our knowledge, this is the first case where the lag between the change in continuum and in broad-line emission of a "changing-look" AGN has been temporally resolved. The nature and timescales of the photometric and spectral evolution disfavor both a change in line-of-sight obscuration and a change of the overall rate of gas inflow as driving the drastic spectral transformations seen in this AGN. Although the peak luminosity and timescales are consistent with those of tidal disruption events seen in inactive galaxies, the spectral properties are not. The X-ray emission displays a markedly different behavior, with frequent flares on timescales of hours to days, and will be presented in a companion publication.

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

By comparing the dynamical and lensing masses of early-type lens galaxies, one can constrain both the cosmological parameters and the density profiles of galaxies. We explore the constraining power on cosmological parameters and the effect of the lens mass model in this method with 161 galaxy-scale strong lensing systems, which is currently the largest sample with both high resolution imaging and stellar dynamical data. We assume a power-law mass model for the lenses, and consider three different parameterizations for $\gamma$ (i.e., the slope of the total mass density profile) to include the effect of the dependence of $\gamma$ on redshift and surface mass density. When treating $\delta$ (i.e., the slope of the luminosity density profile) as a universal parameter for all lens galaxies, we find the limits on the cosmological parameter $\Omega_m$ are quite weak and biased, and also heavily dependent on the lens mass model in the scenarios of parameterizing $\gamma$ with three different forms. When treating $\delta$ as an observable for each lens, the unbiased estimate of $\Omega_m$ can be obtained only in the scenario of including the dependence of $\gamma$ on both the redshift and the surface mass density, that is $\Omega_m = 0.381^{+0.185}_{-0.154}$ at 68\% confidence level in the framework of a flat $\Lambda$CDM model. We conclude that the significant dependencies of $\gamma$ on both the redshift and the surface mass density, as well as the intrinsic scatter of $\delta$ among the lenses, need to be properly taken into account in this method.

A critical step toward the emergence of planets in a protoplanetary disk consists in accretion of planetesimals, bodies 1-1000 km in size, from smaller disk constituents. This process is poorly understood partly because we lack good observational constraints on the complex physical processes that contribute to planetesimal formation. In the outer solar system, the best place to look for clues is the Kuiper belt, where icy planetesimals survived to this day. Here we report evidence that Kuiper belt planetesimals formed by the streaming instability, a process in which aerodynamically concentrated clumps of pebbles gravitationally collapse into 100-km-class bodies. Gravitational collapse was previously suggested to explain the ubiquity of equal-size binaries in the Kuiper belt. We analyze new hydrodynamical simulations of the streaming instability to determine the model expectations for the spatial orientation of binary orbits. The predicted broad inclination distribution with 80% of prograde binary orbits matches the observations of trans-Neptunian binaries. The formation models which imply predominantly retrograde binary orbits can be ruled out. Given its applicability over a broad range of protoplanetary disk conditions, it is expected that the streaming instability seeded planetesimal formation also elsewhere in the solar system, and beyond.

A Rayleigh-Ritz with Continuous Galerkin method based approach is presented to compute the normal modes of a planet in the presence of an essential spectrum. The essential spectrum is associated with a liquid outer core. The presence of a liquid outer core requires the introduction of a mixed Continuous Galerkin finite-element approach. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid domains. The relevant generalized eigenvalue problem is solved by a combination of several highly parallel, computationally efficient methods. Self-gravitation is treated as an N-body problem and the relevant gravitational potential is evaluated directly and efficiently utilizing the fast multipole method. The computational experiments are performed on constant elastic balls and the isotropic version of the preliminary reference earth model (PREM) for validation. Our proposed algorithm is illustrated in fully heterogeneous models including one combined with crust 1.0.

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single epoch images, coadded images, coadded source catalogs, and associated products and services assembled over the first three years of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a photometric precision of < 1% in all bands, and an astrometric precision of 151 mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1 includes nearly 400M distinct astronomical objects detected in ~10,000 coadd tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures. Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive coadd image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

Recently, Rastall gravity is undergoing a significant surge in popularity. We obtain a power-law total mass-density profile for the inner region (within several effective radius) of early-type galaxies (ETGs) from the space-time structures which are described by the static spherically-symmetric solutions of Rastall gravity under the assumption of perfect fluid matter. We find that in the inner region of ETGs, the Rastall dimensionless parameter $\beta=\kappa\lambda$ determines the mass distribution. We then use 118 galaxy-galaxy strong gravitational lensing systems to constrain the Rastall dimensionless parameter $\beta$. We find that the mean value of $\beta$ for total 118 ETGs is $\beta=0.163\pm0.001$(68\% CL) with a minor intrinsic scatter of $\delta=0.020\pm 0.001$. Our work observationally illustrates the physical meaning of the Rastall dimensionless parameter in galaxy scale. From the Newtonian approximation of Rastall gravity, we also find that an absolute isothermal mass distribution for ETGs is not allowed in the framework of Rastall gravity.

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1.

The interstellar medium is a key ingredient that governs star formation in galaxies. We present a detailed study of the infrared (~ 1-500 micron) spectral energy distributions of a large sample of 193 nearby (z ~ 0.088) luminous infrared galaxies (LIRGs) covering a wide range of evolutionary stages along the merger sequence. The entire sample has been observed uniformly by 2MASS, WISE, Spitzer, and Herschel. We perform multi-component decomposition of the spectra to derive physical parameters of the interstellar medium, including the intensity of the interstellar radiation field and the mass and luminosity of the dust. We also constrain the presence and strength of nuclear dust heated by active galactic nuclei. The radiation field of LIRGs tends to have much higher intensity than in quiescent galaxies, and it increases toward advanced merger stages as a result of central concentration of the interstellar medium and star formation. The total gas mass is derived from the dust mass and the galaxy stellar mass. We find that the gas fraction of LIRGs is on average ~ 0.3 dex higher than that of main-sequence star-forming galaxies, rising moderately toward advanced merger stages. All LIRGs have star formation rates that place them above the galaxy star formation main sequence. Consistent with recent observations and numerical simulations, the global star formation efficiency of the sample spans a wide range, filling the gap between normal star-forming galaxies and extreme starburst systems.

We naturally obtain a power-law total mass-density profile for the inner regions (within several effective radius) of early-type galaxies(ETGs) from space-time structures, which are described by the static spherically-symmetric solutions of Einstein equations with the perfect fluid matters. We point out that, in the inner region of an early type galaxy, the average value of the equation of state $\omega$ determines the power-law index $\gamma$ of the total mass-density profile given by \cite{2006ApJ...649..599K}. We then use 119 galaxy-galaxy strong gravitational lensing systems from SLACS, BELLS, BELLS GALLERY and SL2S, to constrain the equation of state $\omega$ for the inner regions of ETGs. We find that the average value of $\omega$ for totally 119 ETGs is $\omega=-0.315\pm0.073$ (68\% CL), which is very close to $-1/3$. If $\omega=-1/3$ is the best value to describe the equation of state of the inner regions of the dark matter in galaxies, the deviation from $-1/3$ could represent the mixture of baryonic matter and dark matter in galaxies.

More than one hundred galaxy-scale strong gravitational lens systems have been found by searching for the emission lines coming from galaxies with redshifts higher than the lens galaxies. Based on this spectroscopic-selection method, we introduce the deep Residual Networks (ResNet, a kind of deep Convolutional Neural Networks) to search for the galaxy-Ly$\alpha$ emitter (LAE) lens candidates by recognizing the Ly$\alpha$ emission lines coming from high redshift galaxies ($2 < z < 3$) in the spectra of early-type galaxies (ETGs) at middle redshift ($z\sim 0.5$). The spectra of the ETGs come from the Data Release 12 (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey \uppercase\expandafter{\romannumeral3} (SDSS-\uppercase\expandafter{\romannumeral3}). In this paper, we first build a 28 layers ResNet model, and then artificially synthesize 150,000 training spectra, including 140,000 spectra without Ly$\alpha$ lines and 10,000 ones with Ly$\alpha$ lines, to train the networks. After 20 training epochs, we obtain a near-perfect test accuracy at 0.9954. The corresponding loss is 0.0028 and the completeness is 93.6\%. We finally apply our ResNet model to our predictive data with 174 known lens candidates. We obtain 1232 hits including 161 of the 174 known candidates (92.5\% discovery rate). Apart from the hits found in other works, our ResNet model also find 536 new hits. We then perform several subsequent selections on these 536 hits and present 5 most believable lens candidates.

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

We measure the total mass-density profiles out to three effective radii for a sample of 63 $z \sim 0.5$, massive early-type galaxies (ETGs) acting as strong gravitational lenses through a joint analysis of lensing and stellar dynamics. The compilation is selected from three galaxy-scale strong-lens samples including the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS), BELLS for GALaxy-Ly$\alpha$ EmitteR sYstems Survey, and Strong Lensing Legacy Survey (SL2S). Utilizing the wide source-redshift coverage (0.8--3.5) provided by these three samples, we build a statistically significant ensemble of massive ETGs for which robust mass measurements can be achieved within a broad range of Einstein radii up to three effective radii. Characterizing the three-dimensional total mass-density distribution by a power-law profile as $\rho \propto r^{-\gamma}$, we find that the average logarithmic density slope for the entire sample is $\langle\gamma\rangle=2.000_{-0.032}^{+0.033}$ ($68\%$CL) with an intrinsic scatter $\delta=0.180_{-0.028}^{+0.032}$. Further parameterizing $\langle\gamma\rangle$ as a function of redshift $z$ and ratio of Einstein radius to effective radius $R_{ein}/R_{eff}$, we find the average density distributions of these massive ETGs become steeper at larger radii and later cosmic times with magnitudes $\mathrm{d} \langle\gamma\rangle / \mathrm{d}z = -0.309_{-0.160}^{+0.166}$ and $\mathrm{d} \langle\gamma\rangle / \mathrm{d} \log_{10} \frac{R_{ein}}{R_{eff}} = 0.194_{-0.083}^{+0.092}$.

The Streaming Instability (SI) is a mechanism to concentrate solids in protoplanetary disks. Nonlinear particle clumping from the SI can trigger gravitational collapse into planetesimals. To better understand the numerical robustness of the SI, we perform a suite of vertically-stratified 3D simulations with fixed physical parameters known to produce strong clumping. We vary the numerical implementation, namely the computational domain size and the vertical boundary conditions (vBCs), comparing newly-implemented outflow vBCs to the previously-used periodic and reflecting vBCs. We find strong particle clumping by the SI is mostly independent of the vBCs. However, peak particle densities are higher in larger simulation domains due to a larger particle mass reservoir. We report SI-triggered zonal flows, i.e., azimuthally-banded radial variations of gas pressure. These structures have low amplitudes, insufficient to halt particle radial drift, confirming that particle trapping in gas pressure maxima is not the mechanism of the SI. We find that outflow vBCs produce artificially large gas outflow rates at vertical boundaries. However, the outflow vBCs reduce artificial reflections at vertical boundaries, allowing more particle sedimentation, and showing less temporal variation and better convergence with box size. The radial spacing of dense particle filaments is $\sim0.15$ gas scale heights ($H$) for all vBCs, which sets the feeding zone for planetesimal growth in self-gravitating simulations. Our results validate the use of the outflow vBCs in SI simulations, even with vertical boundaries close ($\leq 0.4H$) to the disk midplane. Overall, our study demonstrates the numerical robustness of nonlinear particle clumping by the SI.

Schwinn et al. (2017) have recently compared the abundance and distribution of massive substructures identified in a gravitational lensing analysis of Abell 2744 by Jauzac et al. (2016) and N-body simulation and found no cluster in {\Lambda}CDM simulation that is similar to Abell 2744. Schwinn et al.(2017) identified the measured projected aperture masses with the actual masses associated with subhaloes in the MXXL N-body simulation. We have used the high resolution Phoenix cluster simulations to show that such an identification is incorrect: the aperture mass is dominated by mass in the body of the cluster that happens to be projected along the line-of-sight to the subhalo. This enhancement varies from factors of a few to factors of more than 100, particularly for subhaloes projected near the centre of the cluster. We calculate aperture masses for subhaloes in our simulation and compare them to the measurements for Abell 2744. We find that the data for Abell 2744 are in excellent agreement with the matched predictions from {\Lambda}CDM. We provide further predictions for aperture mass functions of subhaloes in idealized surveys with varying mass detection thresholds.

We investigate the discrepancy between the two-dimensional projected lensing mass and the dynam- ical mass for an ensemble of 97 strong gravitational lensing systems discovered by the Sloan Lens ACS (SLACS) Survey, the BOSS Emission-Line Lens Survey (BELLS), and the BELLS for GALaxy-Ly{\alpha} EmitteR sYs- tems (BELLS GALLERY) Survey. We fit the lensing data to obtain the Einstein mass, and use the velocity dispersion of the lensing galaxies provided by the Sloan Digital Sky Survey (SDSS) to get the projected dynamical mass within the Einstein radius by assuming the power-law mass approximation. The discrepancy is found to be obvious and quantified by Bayesian analysis. For the singular isothermal sphere (SIS) mass model, we obtain that the Einstein mass is 20.7% more than the dynamical mass, and the discrepancy increases with the redshift of the lensing galaxies. For more general power-law mass model, the discrepancy still exists within 1{\sigma} credible region. The main reason of the the discrepancy could be mass contamination, including all invisible masses along the line of sight. In addition, the measurement errors and the approximation of the mass models could lead to part of the discrepancy.

We study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

The spatial distribution of satellite galaxies encodes rich information of the structure and assembly history of galaxy clusters. In this paper, we select a redMaPPer cluster sample in SDSS Stripe 82 region with $0.1 \le z \le 0.33$, $20<\lambda<100$ and $P_{\rm cen}>0.7$. Using the high-quality weak lensing data from CS82 Survey, we constrain the mass profile of this sample. Then we compare directly the mass density profile with the satellite number density profile. We find that the total mass and number density profiles have the same shape, both well fitted by an NFW profile. The scale radii agree with each other within 1$\sigma$ error ($r_{\rm s,gal}=0.34_{-0.03}^{+0.04}$Mpc vs $r_{\rm s}=0.37_{-0.10}^{+0.15}$Mpc ).

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

Prior to the arrival of Cassini-Huygens, aerosol production in Titan's atmosphere was believed to begin in the stratosphere where chemical processes are predominantly initiated by FUV radiation. However, measurements taken by Cassini UVIS and CAPS indicate that haze formation initiates in the thermosphere where there is a greater flux of EUV photons and energetic particles available to initiate chemical reactions, including the destruction of N2. The discovery of previously unpredicted nitrogen species in Titan's atmosphere by Cassini INMS indicates that nitrogen participates in the chemistry to a much greater extent than was appreciated before Cassini. The degree of nitrogen in the haze is important for understanding the diversity of molecules present in Titan's atmosphere and on its surface. We have conducted a series of simulation experiments using either spark discharge or FUV photons to initiate chemistry in CH4/N2 gas mixtures (0.01% CH4/99.99% N2 to 10% CH4/90% N2). We obtained in situ real-time measurements using an HR-ToF-AMS to measure the particle composition as a function of particle size and a PIT-MS to measure the composition of gas phase products. These two techniques allow us to investigate the effect of energy source and initial CH4 concentration on the degree of nitrogen incorporation in both the gas and solid phase products. The results presented here confirm that FUV photons produce not only solid phase nitrogen bearing products but also gas phase nitrogen species. We find that in both the gas and solid phase, nitrogen is found in nitriles rather than amines and both the gas phase and solid phase products are composed primarily of molecules with a low degree of aromaticity. The UV experiments reproduce the absolute abundances measured in Titan's stratosphere for a number of gas phase species including C4H2, C6H6, HCN, CH3CN, HC3N, and C2H5CN. [Edited for length]

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2<z<0.4$ and $A=6.61\pm0.75, B=-0.15\pm0.05$ for $0.4<z<0.6$. These amplitudes and slopes are completely consistent with predictions from recent simulations. We also measure the stellar-to-halo mass ratio for our samples, and find results consistent with previous measurements from lensing and other techniques.

Planetesimals may form from the gravitational collapse of dense particle clumps initiated by the streaming instability. We use simulations of aerodynamically coupled gas-particle mixtures to investigate whether the properties of planetesimals formed in this way depend upon the sizes of the particles that participate in the instability. Based on three high resolution simulations that span a range of dimensionless stopping time $6 \times 10^{-3} \leq \tau \leq 2$ no statistically significant differences in the initial planetesimal mass function are found. The mass functions are fit by a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, with $p=1.5-1.7$ and errors of $\Delta p \approx 0.1$. Comparing the particle density fields prior to collapse, we find that the high wavenumber power spectra are similarly indistinguishable, though the large-scale geometry of structures induced via the streaming instability is significantly different between all three cases. We interpret the results as evidence for a near-universal slope to the mass function, arising from the small-scale structure of streaming-induced turbulence.

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

Simulations are expected to be the powerful tool to investigate the baryon effects on dark matter (DM) halos. Recent high resolution, cosmological hydrodynamic simulations (\citealt{Cintio14}, DC14) predict that the inner density profiles of DM halos depend systematically on the ratio of stellar to DM mass ($M_{\ast}/M_{\rm halo}$) which is thought to be able to provide good fits to the observed rotation curves of galaxies. The DC14 profile is fitted from the simulations which are confined to $M_{\rm halo}\le 10^{12}M_{\sun}$, in order to investigate the physical processes that may affect all halos, we extrapolate it to much larger halo mass, including that of galaxy clusters. The inner slope of DC14 profile is flat for low halo mass, it approaches 1 when the halo mass increases towards $10^{12}M_{\sun}$ and decreases rapidly after that mass. We use DC14 profile for lenses and find that it predicts too few lenses compared with the most recent strong lensing observations SQLS (\citealt{Inada12}). We also calculate the strong lensing probabilities for a simulated density profile which continues the halo mass from the mass end of DC14 ($\sim 10^{12}M_{\sun}$) to the mass that covers the galaxy clusters (\citealt{Schaller15}, Schaller15), and find that this Schaller15 model predict too many lenses compared with other models and SQLS observations. Interestingly, Schaller15 profile has no core, however, like DC14, the rotation curves of the simulated halos are in excellent agreement with observational data. Furthermore, we show that the standard two-population model SIS+NFW cannot match the most recent SQLS observations for large image separations.

We present a new tool for colour-magnitude diagram (CMD) studies, $Powerful~CMD$. This tool is built on the basis of the advanced stellar population synthesis (ASPS) model, in which single stars, binary stars, rotating stars, and star formation history have been taken into account. Via $Powerful~CMD$, the distance modulus, colour excess, metallicity, age, binary fraction, rotating star fraction, and star formation history of star clusters can be determined simultaneously from observed CMDs. The new tool is tested via both simulated and real star clusters. Five parameters of clusters NGC6362, NGC6652, NGC6838 and M67 are determined and compared to other works. It is shown that this tool is useful for CMD studies, in particular for those with the data of the Hubble Space Telescope (HST). Moreover, we find that the inclusion of binaries in theoretical stellar population models may lead to smaller colour excess compared to the case of single star population models.

We present Canada-France-Hawaii Telescope (CFHT) MegaCam observations of a galaxy-quasar strong gravitational lens system, SDSS J1640+1932. This system, located at z=0.195 (foreground elliptical galaxy) and z=0.778 (background quasar), was first visually identified by us in the Sloan Digital Sky Survey (SDSS) database. Our CFHT imaging with an angular resolution of 0.7$^{\prime \prime}$ clearly resolves 4 lensed images and a nearly complete Einstein ring. Modeling the system with a singular isothermal ellipsoid (SIE) total mass distribution, we find an Einstein radius of ${2.49^{\prime \prime}}_{-0.049}^{+0.063}$ enclosing a inferred mass of $7.25_{-0.29}^{+0.37}\times10^{11} M_{\odot}$. The quasar and its host galaxy have been magnified by a factor of 23, and the time delay relative to the leading image is determined to be 23.4-25.2 days. These parameters vary minimally when our model is fitted to the {\it{g}}-, {\it{r}}- or {\it{i}}-band images.

We perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

Planet searches using the radial velocity technique show a paucity of companions to solar-type stars within ~5 AU in the mass range of ~10 - 80 M$_{\text{Jup}}$. This deficit, known as the brown dwarf desert, currently has no conclusive explanation. New substellar companions in this region help asses the reality of the desert and provide insight to the formation and evolution of these objects. Here we present 10 new brown dwarf and two low-mass stellar companion candidates around solar-type stars from the Multi-object APO Radial-Velocity Exoplanet Large-Area Survey (MARVELS) of the Sloan Digital Sky Survey III (SDSS-III). These companions were selected from processed MARVELS data using the latest University of Florida Two Dimensional (UF2D) pipeline, which shows significant improvement and reduction of systematic errors over previous pipelines. The 10 brown dwarf companions range in mass from ~13 to 76 M$_{\text{Jup}}$ and have orbital radii of less than 1 AU. The two stellar companions have minimum masses of ~98 and 100 M$_{\text{Jup}}$. The host stars of the MARVELS brown dwarf sample have a mean metallicity of [Fe/H] = 0.03 $\pm$ 0.08 dex. Given our stellar sample we estimate the brown dwarf occurrence rate around solar-type stars with periods less than ~300 days to be ~0.56%.

The defining characteristic of the cold dark matter (CDM) hypothesis is the presence of a very large number of low-mass haloes, too small to have made a visible galaxy. Other hypotheses for the nature of the dark matter, such as warm dark matter (WDM), predict a much smaller number of such low-mass haloes. Strong lensing systems offer the possibility of detecting small-mass haloes through the distortions they induce in the lensed image. Here we show that the main contribution to the image distortions comes from haloes along the line of sight rather than subhaloes in the lens as has normally been assumed so far. These interlopers enhance the differences between the predictions of CDM and WDM models. We derive the total perturber mass function, including both subhaloes and interlopers, and show that measurements of approximately 20 strong lens systems with a detection limit of $M_{\rm low}=10^7 h^{-1} M_{\odot}$ would distinguish (at 3 sigma) between CDM and a WDM model consisting of 7 keV sterile neutrinos such as those required to explain the recently detected 3.5 keV X-ray emission line from the centres of galaxies and clusters.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

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

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

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

As the first paper in a series on the study of the galaxy-galaxy lensing from Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our image processing pipeline that corrects the systematics primarily introduced by the Point Spread Function (PSF). Using this pipeline, we processed SDSS DR7 imaging data in $r$ band and generated a background galaxy catalog containing the shape information of each galaxy. Based on our own shape measurements of the galaxy images from SDSS DR7, we extract the galaxy-galaxy (GG) lensing signals around foreground spectroscopic galaxies binned in different luminosity and stellar mass. The overall signals are in good agreement with those obtained by \citet{Mandelbaum2005, Mandelbaum2006} from the SDSS DR4. The results in this paper with higher signal to noise ratio is due to the larger survey area than SDSS DR4, confirm that more luminous/massive galaxies bear stronger GG lensing signal. We also divide the foreground galaxies into red/blue and star forming/quenched subsamples and measured their GG lensing signals, respectively. We find that, at a specific stellar mass/luminosity, the red/quenched galaxies have relatively stronger GG lensing signals than their counterparts especially at large radii. These GG lensing signals can be used to probe the galaxy-halo mass relations and their environmental dependences in the halo occupation or conditional luminosity function framework.

The cold dark matter (CDM) cosmological model unambigously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of warm particles (WDM), then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs.Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of $10^7 h^{-1} M_{\odot}$ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters.

We study the formation of planetesimals in protoplanetary disks from the gravitational collapse of solid over-densities generated via the streaming instability. To carry out these studies, we implement and test a particle-mesh self-gravity module for the Athena code that enables the simulation of aerodynamically coupled systems of gas and collisionless self-gravitating solid particles. Upon employment of our algorithm to planetesimal formation simulations, we find that (when a direct comparison is possible) the Athena simulations yield predicted planetesimal properties that agree well with those found in prior work using different numerical techniques. In particular, the gravitational collapse of streaming-initiated clumps leads to an initial planetesimal mass function that is well-represented by a power-law, dN/dM ~ M^{-p}, with p = 1.6 +/- 0.1, which equates to a differential size distribution dN/dR ~ R^{-q}, with q = 2.8 +/- 0.1. We find no significant trends with resolution from a convergence study of up to 512^3 grid zones and Npar = 1.5x10^8 particles. Likewise, the power-law slope appears indifferent to changes in the relative strength of self-gravity and tidal shear, and to the time when (for reasons of numerical economy) self-gravity is turned on, though the strength of these claims is limited by small number statistics. For a typically assumed radial distribution of minimum mass solar nebula solids (assumed here to have dimensionless stopping time \tau = 0.3), our results support the hypothesis that bodies on the scale of large asteroids or Kuiper Belt Objects could have formed as the high-mass tail of a primordial planetesimal population.

We use the shear catalog from the CFHT Stripe-82 Survey to measure the subhalo masses of satellite galaxies in redMaPPer clusters. Assuming a Chabrier Initial Mass Function (IMF) and a truncated NFW model for the subhalo mass distribution, we find that the sub-halo mass to galaxy stellar mass ratio increases as a function of projected halo-centric radius $r_p$, from $M_{\rm sub}/M_{\rm star}=4.43^{+ 6.63}_{- 2.23}$ at $r_p \in [0.1,0.3]$ $h^{-1}Mpc$ to $M_{\rm sub}/M_{\rm star}=75.40^{+ 19.73}_{- 19.09}$ at $r_p \in [0.6,0.9]$ $h^{-1}Mpc$. We also investigate the dependence of subhalo masses on stellar mass by splitting satellite galaxies into two stellar mass bins: $10<\log(M_{\rm star}/M_{\rm sun})<10.5$ and $11<\log(M_{\rm star}/M_{\rm sun})<12$. The best-fit subhalo mass of the more massive satellite galaxy bin is larger than that of the less massive satellites: $\log(M_{\rm sub}/M_{\rm sun})=11.14 ^{+ 0.66 }_{- 0.73}$ ($M_{\rm sub}/M_{\rm star}=19.5^{+19.8}_{-17.9}$) versus $\log(M_{\rm sub}/M_{\rm sun})=12.38 ^{+ 0.16 }_{- 0.16}$ ($M_{\rm sub}/M_{\rm star}=21.1^{+7.4}_{-7.7}$).

We assess the effectiveness of the Jeans-Anisotropic-MGE (JAM) technique with a state-of-the-art cosmological hydrodynamic simulation, the Illustris project. We perform JAM modelling on 1413 simulated galaxies with stellar mass M^* > 10^{10}M_{sun}, and construct an axisymmetric dynamical model for each galaxy. Combined with a Markov Chain Monte Carlo (MCMC) simulation, we recover the projected root-mean-square velocity (V_rms) field of the stellar component, and investigate constraints on the stellar mass-to-light ratio, M^*/L, and the fraction of dark matter f_{DM} within 2.5 effective radii (R_e). We find that the enclosed total mass within 2.5 R_e is well constrained to within 10%. However, there is a degeneracy between the dark matter and stellar components with correspondingly larger individual errors. The 1 sigma scatter in the recovered M^*/L is 30-40% of the true value. The accuracy of the recovery of M^*/L depends on the triaxial shape of a galaxy. There is no significant bias for oblate galaxies, while for prolate galaxies the JAM-recovered stellar mass is on average 18% higher than the input values. We also find that higher image resolutions alleviate the dark matter and stellar mass degeneracy and yield systematically better parameter recovery.

We explore the Minkowski functionals of weak lensing convergence map to distinguish between $f(R)$ gravity and the general relativity (GR). The mock weak lensing convergence maps are constructed with a set of high-resolution simulations assuming different gravity models. It is shown that the lensing MFs of $f(R)$ gravity can be considerably different from that of GR because of the environmentally dependent enhancement of structure formation. We also investigate the effect of lensing noise on our results, and find that it is likely to distinguish F5, F6 and GR gravity models with a galaxy survey of $\sim3000$ degree$^2$ and with a background source number density of $n_g=30~{\rm arcmin}^{-2}$, comparable to an upcoming survey dark energy survey (DES). We also find that the $f(R)$ signal can be partially degenerate with the effect of changing cosmology, but combined use of other observations, such as the cosmic microwave background (CMB) data, can help break this degeneracy.

We derived constraints on cosmological parameters using weak lensing peak statistics measured on the $\sim130~{\rm deg}^2$ of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). This analysis demonstrates the feasibility of using peak statistics in cosmological studies. For our measurements, we considered peaks with signal-to-noise ratio in the range of $\nu=[3,6]$. For a flat $\Lambda$CDM model with only $(\Omega_{\rm m}, \sigma_8)$ as free parameters, we constrained the parameters of the following relation $\Sigma_8=\sigma_8(\Omega_{\rm m}/0.27)^{\alpha}$ to be: $\Sigma_8=0.82 \pm 0.03 $ and $\alpha=0.43\pm 0.02$. The $\alpha$ value found is considerably smaller than the one measured in two-point and three-point cosmic shear correlation analyses, showing a significant complement of peak statistics to standard weak lensing cosmological studies. The derived constraints on $(\Omega_{\rm m}, \sigma_8)$ are fully consistent with the ones from either WMAP9 or Planck. From the weak lensing peak abundances alone, we obtained marginalised mean values of $\Omega_{\rm m}=0.38^{+0.27}_{-0.24}$ and $\sigma_8=0.81\pm 0.26$. Finally, we also explored the potential of using weak lensing peak statistics to constrain the mass-concentration relation of dark matter halos simultaneously with cosmological parameters.

It has been realized that the shock wave effects play an important role in neutrino oscillations during the supernova explosion. In recent years, with the development of simulations about supernova explosion, we have a better understanding about the density profiles and the shock waves in supernovae than before. It has been shown that the appearance of shock waves not only varies with time, but is also affected by the mass of the supernova. When the mass of the supernova happens to be in a certain range (e.g. it equals 10.8 times the mass of the sun), there might be a reverse shock wave, another sudden change of density except the forward shock wave, emerging in the supernova. In addition, there are some other time-dependent changes of density profiles in different supernova models. Because of these complex density profiles, the expression of the crossing probability at the high resonance, $P_H$, which we used previously would be no longer applicable. In order to get more accurate and reasonable results, we use the data of density profiles in three different supernova models obtained from simulations to study the variations of $P_s$ (the survival probability of $\nu_e\rightarrow\nu_e$), as well as $P_c$ (the conversion probability of $\nu_x\rightarrow\nu_e$). It is found that the mass of the supernova does make a difference on the behavior of $P_s$, and affects $P_c$ at the same time. With the results of $P_s$ and $P_c$, we can estimate the number of $\nu_e$ remained after they go through the matter in the supernova.

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

The New Vacuum Solar Telescope (NVST) is a 1 meter vacuum solar telescope that aims to observe the fine structures on the Sun. The main tasks of NVST are high resolution imaging and spectral observations, including the measurements of solar magnetic field. NVST is the primary ground-based facility of Chinese solar community in this solar cycle. It is located by the Fuxian Lake of southwest China, where the seeing is good enough to perform high resolution observations. In this paper, we first introduce the general conditions of Fuxian Solar Observatory and the primary science cases of NVST. Then, the basic structures of this telescope and instruments are described in detail. Finally, some typical high resolution data of solar photosphere and chromosphere are also shown.

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 select satellite galaxies from the galaxy group catalog constructed with the SDSS spectroscopic galaxies and measure the tangential shear around these galaxies with source catalog extracted from CFHT/MegaCam Stripe-82 Survey to constrain the mass of subhalos associated with them. The lensing signal is measured around satellites in groups with masses in the range [10^{13}, 5x10^{14}]h^{-1}M_{sun}, and is found to agree well with theoretical expectation. Fitting the data with a truncated NFW profile, we obtain an average subhalo mass of log M_{sub}= 11.68 \pm 0.67 for satellites whose projected distances to central galaxies are in the range [0.1, 0.3] h^{-1}Mpc, and log M_{sub}= 11.68 \pm 0.76 for satellites with projected halo-centric distance in [0.3, 0.5] h^{-1}Mpc. The best-fit subhalo masses are comparable to the truncated subhalo masses assigned to satellite galaxies using abundance matching and about 5 to 10 times higher than the average stellar mass of the lensing satellite galaxies.

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

We report the discovery of a candidate brown dwarf or a very low mass stellar companion (MARVELS-5b) to the star HIP 67526 from the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The radial velocity curve for this object contains 31 epochs spread over 2.5 years. Our Keplerian fit using a Markov Chain Monte Carlo approach, reveals that the companion has an orbital period of $90.2695^{+0.0188}_{-0.0187}$ days, an eccentricity of $0.4375 \pm 0.0040$ and a semi-amplitude of $2948.14^{+16.65}_{-16.55}$ m s$^{-1}$. Using additional high-resolution spectroscopy, we find the host star has an effective temperature $T_{\rm{eff}}=6004 \pm 34$ K, a surface gravity $\log g$ [cgs] $=4.55 \pm 0.17$ and a metallicity [Fe/H] $=+0.04 \pm 0.06$. The stellar mass and radius determined through the empirical relationship of Torres et al. (2010), yields 1.10$\pm$0.09 $M_{\sun}$ and 0.92$\pm$0.19 $R_{\sun}$. The minimum mass of MARVELS-5b is $65.0 \pm 2.9 M_{Jup}$, indicating that it is likely to be either a brown dwarf or a very low mass star, thus occupying a relatively sparsely-populated region of the mass function of companions to solar-type stars. The distance to this system is 101$\pm$10 pc from the astrometric measurements of Hipparcos. No stellar tertiary is detected in the high-contrast images taken by either FastCam lucky imaging or Keck adaptive optics imaging, ruling out any star with mass greater than 0.2$M_{\sun}$ at a separation larger than 40 AU.

With galaxy groups constructed from the Sloan Digital Sky Survey (SDSS), we analyze the expected galaxy-galaxy lensing signals around satellite galaxies residing in different host haloes and located at different halo-centric distances. We use Markov Chain Monte Carlo (MCMC) method to explore the potential constraints on the mass and density profile of subhaloes associated with satellite galaxies from SDSS-like surveys and surveys similar to the Large Synoptic Survey Telescope (LSST). Our results show that for SDSS-like surveys, we can only set a loose constraint on the mean mass of subhaloes. With LSST-like surveys, however, both the mean mass and the density profile of subhaloes can be well constrained.

Binary stars are common in star clusters and galaxies, but the detailed effects of binary evolution are not taken into account in some colour-magnitude diagram (CMD) studies. This paper studies the CMDs of twelve globular clusters via binary-star stellar populations. The observational CMDs of star clusters are compared to those of binary-star populations, and then the stellar metallicities, ages, distances and reddening values of these star clusters are obtained. The paper also tests the different effects of binary and single stars on CMD studies. It is shown that binaries can fit the observational CMDs of the sample globular clusters better compared to single stars. This suggests that the effects of binary evolution should be considered when modeling the CMDs and stellar populations of star clusters and galaxies.

Binary stars are common in star clusters and galaxies, but the detailed effects of binary evolution are not taken into account in some colour-magnitude diagram (CMD) studies. This paper studies the CMDs of twelve globular clusters via binary-star stellar populations. The observational CMDs of star clusters are compared to those of binary-star populations, and then the stellar metallicities, ages, distances and reddening values of these star clusters are obtained. The paper also tests the different effects of binary and single stars on CMD studies. It is shown that binaries can fit the observational CMDs of the sample globular clusters better compared to single stars. This suggests that the effects of binary evolution should be considered when modeling the CMDs and stellar populations of star clusters and galaxies.

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

We investigate the potential of exploiting the Sunyaev-Zeldovich effect (SZE) to study the properties of hot gas in galaxy groups. It is shown that, with upcoming SZE surveys, one can stack SZE maps around galaxy groups of similar halo masses selected from large galaxy redshift surveys to study the hot gas in halos represented by galaxy groups. We use various models for the hot halo gas to study how the expected SZE signals are affected by gas fraction, equation of state, halo concentration, and cosmology. Comparing the model predictions with the sensitivities expected from the SPT, ACT and Planck surveys shows that a SPT-like survey can provide stringent constraints on the hot gas properties for halos with masses M ~> 10^{13} h^{-1}Msun. We also explore the idea of using the cross correlation between hot gas and galaxies of different luminosity to probe the hot gas in dark matter halos without identifying galaxy groups to represent dark halos. Our results show that, with a galaxy survey as large as the Sloan Digital Sky Survey and with the help of the conditional luminosity function (CLF) model, one can obtain stringent constraints on the hot gas properties in halos with masses down to 10^{13} h^{-1}Msun. Thus, the upcoming SZE surveys should provide a very promising avenue to probe the hot gas in relatively low-mass halos where the majority of L*-galaxies reside.

Context. Absorption line indices are widely used to determine the stellar population parameters such as age and metallicity of galaxies, but it is not easy to obtain the line indices of some distant galaxies that have colours available. Aims. This paper investigates the correlations between absorption line indices and colours. Methods. A few statistical fitting methods are mainly used, via both the observational data of Sloan Digital Sky Survey and a widely used theoretical stellar population model. Results. Some correlations between widely used absorption line indices and ugriz colours are found from both observational data of early-type galaxies and a theoretical simple stellar population model. In particular, good correlations between colours and widely used absorption line indices such as Dn(4000), HgammaA, HgammaF, HdeltaA, Mg1, Mg2, and Mgb, are shown in this paper. Conclusions. Some important absorption line indices of early-type galaxies can be estimated from their colours using correlations between absorption line indices and colours. For example, age-sensitive absorption line indices can be estimated from (u-r) or (g-r) colours and metallicity-sensitive ones from (u - z) or (g - z). This is useful for studying the stellar populations of distant galaxies, especially for statistical investigations.

Context. Absorption line indices are widely used to determine the stellar population parameters such as age and metallicity of galaxies, but it is not easy to obtain the line indices of some distant galaxies that have colours available. Aims. This paper investigates the correlations between absorption line indices and colours. Methods. A few statistical fitting methods are mainly used, via both the observational data of Sloan Digital Sky Survey and a widely used theoretical stellar population model. Results. Some correlations between widely used absorption line indices and ugriz colours are found from both observational data of early-type galaxies and a theoretical simple stellar population model. In particular, good correlations between colours and widely used absorption line indices such as Dn(4000), HgammaA, HgammaF, HdeltaA, Mg1, Mg2, and Mgb, are shown in this paper. Conclusions. Some important absorption line indices of early-type galaxies can be estimated from their colours using correlations between absorption line indices and colours. For example, age-sensitive absorption line indices can be estimated from (u-r) or (g-r) colours and metallicity-sensitive ones from (u - z) or (g - z). This is useful for studying the stellar populations of distant galaxies, especially for statistical investigations.

We use galaxy groups selected from the Sloan Digital Sky Survey (SDSS) together with mass models for individual groups to study the galaxy-galaxy lensing signals expected from galaxies of different luminosities and morphological types. We compare our model predictions with the observational results obtained from the SDSS by Mandelbaum et al. (2006) for the same samples of galaxies. The observational results are well reproduced in a $\Lambda$CDM model based on the WMAP 3-year data, but a $\Lambda$CDM model with higher $\sigma_8$, such as the one based on the WMAP 1-year data,significantly over-predicts the galaxy-galaxy lensing signal. We model, separately, the contributions to the galaxy-galaxy lensing signals from different galaxies: central versus satellite, early-type versus late-type, and galaxies in halos of different masses. We also examine how the predicted galaxy-galaxy lensing signal depends on the shape, density profile, and the location of the central galaxy with respect to its host halo.

The kinematics of satellite galaxies reflect the masses of the extended dark matter haloes in which they orbit, and thus shed light on the mass-luminosity relation (MLR) of their corresponding central galaxies. In this paper we select a large sample of centrals and satellites from the Sloan Digital Sky Survey (SDSS) and measure the kinematics (velocity dispersions) of the satellite galaxies as a function of the $r$-band luminosity of the central galaxies. Using the analytical framework presented in Paper I, we use these data to infer {\it both} the mean and the scatter of the MLR of central galaxies, carefully taking account of selection effects and biases introduced by the stacking procedure. As expected, brighter centrals on average reside in more massive haloes. In addition, we find that the scatter in halo masses for centrals of a given luminosity, $\sigma_{\log M}$, also increases with increasing luminosity. As we demonstrate, this is consistent with $\sigma_{\log L}$, which reflects the scatter in the conditional probability function $P(L_c|M)$, being independent of halo mass. Our analysis of the satellite kinematics yields $\sigma_{\log L}=0.16\pm0.04$, in excellent agreement with constraints from clustering and group catalogues, and with predictions from a semi-analytical model of galaxy formation. We thus conclude that the amount of stochasticity in galaxy formation, which is characterized by $\sigma_{\log L}$, is well constrained, is independent of halo mass, and is in good agreement with current models of galaxy formation.

Galaxy clustering and galaxy-galaxy lensing probe the connection between galaxies and their dark matter haloes in complementary ways. On one hand, the halo occupation statistics inferred from the observed clustering properties of galaxies are degenerate with the adopted cosmology. Consequently, different cosmologies imply different mass-to-light ratios for dark matter haloes. On the other hand, galaxy-galaxy lensing yields direct constraints on the actual mass-to-light ratios and it can be used to break this degeneracy, and thus to constrain cosmological parameters. In this paper we establish the link between galaxy luminosity and dark matter halo mass using the conditional luminosity function (CLF). We constrain the CLF parameters using the galaxy luminosity function and the luminosity dependence of the correlation lengths of galaxies. The resulting CLF models are used to predict the galaxy-galaxy lensing signal. For a cosmology with $(\Omega_{\rm m},\sigma_8)=(0.238,0.734)$, the model accurately fits the galaxy-galaxy lensing data obtained from the SDSS. For a comparison cosmology with $(\Omega_{\rm m},\sigma_8)=(0.3,0.9)$, however, we can accurately fit the luminosity function and clustering properties of the galaxy population, but the model predicts mass-to-light ratios that are too high, resulting in a strong overprediction of the galaxy-galaxy lensing signal. We conclude that the combination of galaxy clustering and galaxy-galaxy lensing is a powerful probe of the galaxy-dark matter connection, with the potential to yield tight constraints on cosmological parameters. Since this method mainly probes the mass distribution on non-linear scales, it is complementary to constraints obtained from the galaxy power-spectrum, which mainly probes the large-scale (linear) matter distribution.