We study the properties of magnetic fields in the circumgalactic medium (CGM) of $z<1$ galaxies by correlating Faraday rotation measures of $\sim 1,000$ high-redshift radio sources with the foreground galaxy number density estimated from the DESI Legacy Imaging Surveys. This method enables us to extract signals of rotation measures contributed by intervening gas around multiple galaxies. Our results show that there is no detectable correlation between the distribution of rotation measures and the number of foreground galaxies, contrary to several previous results. Utilizing the non-detection signals, we estimate $3\sigma$ upper limits to the rotation measures from the CGM of $\sim 20 \rm \ rad/m^{2}$ within~50 kpc and $\sim 10 \rm \ rad / m^{2}$ at separations of $100$ kpc. By adopting a column density distribution of ionized gas obtained from absorption line measurements, we further estimate the strengths of coherent magnetic fields parallel to the line of sight of $<\rm 2 \ \mu G$ in the CGM. We show that the estimated upper limits of rotation measures and magnetic field strengths are sufficient to constrain outputs of recent galaxy magneto-hydrodynamic simulations. Finally, we discuss possible causes for the inconsistency between our results and previous works.

We study the evolution of galaxies and the circumgalactic medium (CGM) through cosmic time by correlating $\sim50,000$ MgII absorbers, tracers of cool gas ($\sim10^{4}$~K), detected in the SDSS quasar spectra with galaxies detected in the DESI Legacy Imaging Surveys. By doing so, we extract the properties of galaxies associated with absorbers from redshift 0.4 to 1.3 with effectively $\sim 15,000$ pairs and explore the covering fraction of MgII absorbers as a function of galaxy type, stellar mass, impact parameter, and redshift. We find that the gas covering fraction increases with stellar mass of galaxies by $\sim M_{*}^{0.4}$. However, after we normalize the impact parameter by the virial radius of dark matter halos, the gas profiles around galaxies with masses ranging from $10^{9}$ to $10^{11} \, M_{\odot}$ become weakly dependent on stellar mass. In addition, the gas distribution depends on galaxy type: the covering fraction within $0.3 \, r_{vir}$ around star-forming galaxies is 2-4 times higher than that around passive galaxies at all redshifts. We find that the covering fraction of strong absorbers ($W_{\lambda2796}>1$ Angstrom) around both types of galaxies evolves significantly with redshift, similarly to the evolution of star-formation rate of galaxies, while such an evolution is not detected for weak absorbers ($W_{\lambda2796}<1$ Angstrom). We quantify the HI mass traced by strong absorbers and find that the gas mass around galaxies evolves consistently with the star-formation rate of galaxies. This result suggests that the properties of galaxies and their CGM co-evolve through cosmic time. Finally, we discuss the origins of strong absorbers around passive galaxies and argue that its redshift evolution may trace the star-formation activity of satellite galaxies.

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

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

We develop a semi-analytic model to explore the physical properties of cool pressure-confined circumgalactic clouds with mass ranging from $10$ to $10^{8} \, \rm M_{\odot}$ in a hot diffuse halo. We consider physical effects that control the motion and mass loss of the clouds, and estimate the lifetime and the observed properties of individual cool gas clouds inferred from the CLOUDY simulation. Our results show that the cool pressure-confined gas clouds have physical properties consistent with absorption line systems with neutral hydrogen column densities $N_{\rm HI}\geq10^{18.5} \rm cm^{-2}$ such as strong metal absorbers, sub-DLAs, and DLAs. The cool circumgalactic clouds are transient due to evaporation and recycling and therefore a constant replenishment is needed to maintain the cool CGM. We further model the ensemble properties of the cool CGM with clouds originated from outflows, inflows, or/and in-situ formation with a range of initial cloud mass function and velocity distribution. We find that only with a certain combination of parameters, an outflow model can broadly reproduce three cool gas properties around star-forming galaxies simultaneously: the spatial distribution, down-the-barrel outflow absorption, and gas velocity dispersion. Both a constant insitu model and gas inflow model can reproduce the observed covering fractions of high $\rm N_{HI}$ gas around passive galaxies but they fail to reproduce sufficient number of low $\rm N_{HI}$ systems. The limitations and the failures of the current models are discussed. Our results illustrate that semi-analytic modeling is a promising tool to understand the physics of the cool CGM which is usually unresolved by state-of-the-art cosmological hydrodynamic simulations.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

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

We study the distribution and kinematics of the cool circumgalactic medium (CGM) of emission line galaxies (ELGs) traced by metal absorption lines. Using about 200,000 ELGs from SDSS-IV eBOSS and half a million background quasars from SDSS, we measure the median absorption strength of MgII and FeII lines in quasar spectra for impact parameters ranging from 10 kpc to 1 Mpc. For comparison we measure the same quantity around luminous red galaxies (LRGs). On scales greater than 100 kpc both ELGs and LRGs exhibit similar absorption profiles. However, metal absorption is 5-10 times stronger around ELGs on smaller scales. The metal absorption strength is anisotropic, with an excess along the minor-axis of the galaxies, indicating an outflow origin of the absorbing gas. The ratio between the velocity dispersion of the cool CGM and that of its host dark matter halo is about one for ELGs but about half for LRGs. These results show that the dichotomy of galaxy types is reflected in both the density distribution and kinematics of the CGM traced by metal absorption lines. Our results provide strong evidence that the CGM of ELGs is enriched by gas outflows generated by star formation.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.

We investigate the limitations of statistical absorption measurements with the SDSS optical spectroscopic surveys. We show that changes in the data reduction strategy throughout different data releases have led to a better accuracy at long wavelengths, in particular for sky line subtraction, but a degradation at short wavelengths with the emergence of systematic spectral features with an amplitude of about one percent. We show that these features originate from inaccuracy in the fitting of modeled F-star spectra used for flux calibration. The best-fit models for those stars are found to systematically over-estimate the strength of metal lines and under-estimate that of Lithium. We also identify the existence of artifacts due to masking and interpolation procedures at the wavelengths of the hydrogen Balmer series leading to the existence of artificial Balmer $\alpha$ absorption in all SDSS optical spectra. All these effects occur in the rest-frame of the standard stars and therefore present Galactic longitude variations due to the rotation of the Galaxy. We demonstrate that the detection of certain weak absorption lines reported in the literature are solely due to calibration effects. Finally, we discuss new strategies to mitigate these issues.

Metal abundance and its evolution are studied for Mg II quasar absorption line systems from their weak, unsaturated spectral lines using stacked spectra from the archived data of Sloan Digital Sky Survey. They show an abundance pattern that resembles that of Galactic halo or Small Magellanic Cloud, with metallicity [Z/H] showing an evolution from redshift z=2 to 0.5: metallicity becomes approximately solar or even larger at z~0. We show that the evolution of the metal abundance traces the cumulative amount of the hydrogen fuel consumed in star formation in galaxies. With the aid of a spectroscopic simulation code, we infer the median gas density of the cloud to be roughly 0.3 $\rm cm^{-3}$, with which the elemental abundance in various ionization stages, in particular C I, is consistently explained. This gas density implies that the size of the Mg II clouds is of the order of 0.03 kpc, which suggests that individual Mg II clouds around a galaxy are of a baryonic mass typically $\rm 10^{3} \, M_{\odot}$. This means that Mg II clouds are numerous and `foamy', rather than the large entity that covers a sizable fraction of galaxies with a single cloud.

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

We characterize the luminosity functions of galaxies residing in $z\sim0$ groups and clusters over the broadest ranges of luminosity and mass reachable by the Sloan Digital Sky Survey. Our measurements cover four orders of magnitude in luminosity, down to about $M_r=-12$ mag or $L=10^7\,L_\odot$, and three orders of magnitude in halo mass, from $10^{12}$ to $10^{15} \, {\rm M}_\odot$. We find a characteristic scale, $M_r\sim-18$ mag or $L\sim10^9\, L_\odot$, below which the slope of the luminosity function becomes systematically steeper. This trend is present for all halo masses and originates mostly from red satellites. This ubiquitous faint-end upturn suggests that it is formation, rather than halo-specific environmental effect, that plays a major role in regulating the stellar masses of faint satellites. We show that the satellite luminosity functions can be described in a simple manner by a double Schechter function with amplitudes scaling with halo mass over the entire range of observables. Combining these conditional luminosity functions with the dark matter halo mass function, we accurately recover the entire field luminosity function over 10 visual magnitudes and reveal that satellite galaxies dominate the field luminosity function at magnitudes fainter than $-17$. We find that the luminosity functions of blue and red satellite galaxies show distinct shapes and we present estimates of the stellar mass fraction as a function of halo mass and galaxy type. Finally, using a simple model, we demonstrate that the abundances and the faint-end slopes of blue and red satellite galaxies can be interpreted in terms of their formation history, with two distinct modes separated by some characteristic time.

We use star, galaxy and quasar spectra taken by the Sloan Digital Sky Survey to map out the distribution of diffuse interstellar bands (DIBs) induced by the Milky Way. We show that, after carefully removing the intrinsic spectral energy distribution of each source, it is possible to measure statistical flux fluctuations at the 1e-3 level, detect more than 20 DIBs and measure their strength as a function of position on the sky. We create a map of DIB absorption covering about 5000 square degrees and measure correlations with various tracers of the interstellar medium: atomic and molecular hydrogen, dust and polycyclic aromatic hydrocarbons (PAHs). After recovering known correlations, we show that each DIB has a different dependence on atomic and molecular hydrogen: while they are all positively correlated with N(HI), they exhibit a range of behaviours with N(H2) showing positive, negative or no correlation. We show that a simple parametrization involving only N(HI) and N(H2) applied to all the DIBs is sufficient to reproduce a large collection of observational results reported in the literature: it allows us to naturally describe the relations between DIB strength and dust reddening (including the so-called skin effect), the related scatter, DIB pair-wise correlations & families, the affinity for $\sigma/\zeta$-type environments and other correlations related to molecules. Our approach allows us to characterize DIB dependencies in a simple manner and provides us with a metric to characterize the similarity between different DIBs.

We map out calcium II & sodium I absorption (Fraunhofer H, K & D lines) induced by both the interstellar medium and the circumgalactic medium of the Milky Way. Our measurements cover more than $9000$ deg$^2$ and make use of about $300,000$ extragalactic spectra from the Sloan Digital Sky Survey. We present absorption maps for these two species and then compare their distributions to those of neutral hydrogen and dust. We show that the abundance of Na I with respect to neutral hydrogen stays roughly constant in different environments, while that of Ca II decreases with hydrogen column density. Studying how these tracers vary as a function of velocity, we show that, on average, the N(Na I)/N(Ca II) ratio decreases at higher velocity with respect to the local standard of rest, similar to the local Routly-Spitzer effect but seen on Galactic scale. We show that it is likely caused by higher gas/dust density at lower velocity. Finally, we show that Galactic Ca II and Na I absorption needs to be taken into account for precision photometry and, more importantly, for photometric redshift estimation with star forming galaxies. Our maps of Ca II and Na I absorption are publicly available.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We explore the distribution of cool (~$10^4$K) gas around galaxies and its dependence on galaxy properties. By cross-correlating about 50,000 MgII absorbers with millions of sources from the SDSS (optical), WISE (IR), and GALEX (UV) surveys we effectively extract about 2,000 galaxy-absorber pairs at z~0.5 and probe relations between absorption strength and galaxy type, impact parameter and azimuthal angle. We find that cool gas traced by MgII absorbers exists around both star-forming and passive galaxies with a similar incidence rate on scales greater than 100 kpc but each galaxy type exhibits a different behavior on smaller scales: MgII equivalent width does not correlate with the presence of passive galaxies whereas stronger MgII absorbers tend to be found in the vicinity of star-forming galaxies. This effect is preferentially seen along the minor axis of these galaxies, suggesting that some of the gas is associated with outflowing material. In contrast, the distribution of cool gas around passive galaxies is consistent with being isotropic on the same scales. We quantify the average excess MgII equivalent width $<\delta W_{0}^{\rm MgII}>$ as a function of galaxy properties and find $<\delta W_0^{\rm MgII}>\propto SFR^{1.2}, sSFR^{0.5}$ and $M_\ast^{0.4}$ for star-forming galaxies. This work demonstrates that the dichotomy between star-forming and passive galaxies is reflected in the CGM traced by low-ionized gas. We also measure the covering fraction of MgII absorption and find it to be about 2-10 times higher for star-forming galaxies than passive ones within 50 kpc. We estimate the amount of neutral gas in the halo of $<\log M_\ast/{\rm M_\odot}>$~10.8 galaxies to be a few x$10^9 \rm M_\odot$ for both types of galaxies. Finally, we find that correlations between absorbers and sources detected in the UV and IR lead to physical trends consistent with those measured in the optical.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

This paper is a review on the observational Hubble parameter data that have gained increasing attention in recent years for their illuminating power on the dark side of the universe --- the dark matter, dark energy, and the dark age. Currently, there are two major methods of independent observational H(z) measurement, which we summarize as the "differential age method" and the "radial BAO size method". Starting with fundamental cosmological notions such as the spacetime coordinates in an expanding universe, we present the basic principles behind the two methods. We further review the two methods in greater detail, including the source of errors. We show how the observational H(z) data presents itself as a useful tool in the study of cosmological models and parameter constraint, and we also discuss several issues associated with their applications. Finally, we point the reader to a future prospect of upcoming observation programs that will lead to some major improvements in the quality of observational H(z) data.

The universe, with large-scale homogeneity, is locally inhomogeneous, clustering into stars, galaxies and larger structures. Such property is described by the smoothness parameter $\alpha$ which is defined as the proportion of matter in the form of intergalactic medium. If we take consideration of the inhomogeneities in small scale, there should be modifications of the cosmological distances compared to a homogenous model. Dyer and Roeder developed a second-order ordinary differential equation (D-R equation) that describes the angular diameter distance-redshift relation for inhomogeneous cosmological models. Furthermore, we may obtain the D-R equation for observational $H(z)$ data (OHD). The density-parameter $\Omega_{\rm M}$, the state of dark energy $\omega$, and the smoothness-parameter $\alpha$ are constrained by a set of OHD in a spatially flat $\Lambda$CDM universe as well as a spatially flat XCDM universe. By using of $\chi^2$ minimization method we get $\alpha=0.81^{+0.19}_{-0.20}$ and $\Omega_{\rm M}=0.32^{+0.12}_{-0.06}$ at $1\sigma$ confidence level. If we assume a Gaussian prior of $\Omega_{\rm M}=0.26\pm0.1$, we get $\alpha=0.93^{+0.07}_{-0.19}$ and $\Omega_{\rm M}=0.31^{+0.06}_{-0.05}$. For XCDM model, $\alpha$ is constrained to $\alpha\geq0.80$ but $\omega$ is weakly constrained around -1, where $\omega$ describes the equation of the state of the dark energy ($p_{\rm X}=\omega\rho_{\rm X}$). We conclude that OHD constrains the smoothness parameter more effectively than the data of SNe Ia and compact radio sources.

We use the redshift Hubble parameter $H(z)$ data derived from relative galaxy ages, distant type Ia supernovae (SNe Ia), the Baryonic Acoustic Oscillation (BAO) peak, and the Cosmic Microwave Background (CMB) shift parameter data, to constrain cosmological parameters in the Undulant Universe. We marginalize the likelihood functions over $h$ by integrating the probability density $P\propto e^{-\chi^2/2}$. By using the Markov Chain Monte Carlo (MCMC) technique, we obtain the best fitting results and give the confidence regions on the $b-\Omega_{\rm m0}$ plane. Then we compare their constraints. Our results show that the $H(z)$ data play a similar role with the SNe Ia data in cosmological study. By presenting the independent and joint constraints, we find that the BAO and CMB data play very important roles in breaking the degeneracy compared with the $H(z)$ and SNe Ia data alone. Combined with the BAO or CMB data, one can improve the constraints remarkably. The SNe Ia data sets constrain $\Omega_{\rm m0}$ much tighter than the $H(z)$ data sets, but the $H(z)$ data sets constrain $b$ much tighter than the SNe Ia data sets. All these results show that the Undulant Universe approaches the $\Lambda \rm$CDM model. We expect more $H(z)$ data to constrain cosmological parameters in future.

Recent observations show that our universe is accelerating by dark energy, so it is important to investigate the thermodynamical properties of it. The Undulant Universe is a model with equation of state $\omega(a)=-\cos(b\ln a)$ for dark energy, where we show that there neither the event horizon nor the particle horizon exists. However, as a boundary of keeping thermodynamical properties, the apparent horizon is a good holographic screen. The Universe has a thermal equilibrium inside the apparent horizon, so the Unified First Law and the Generalized Second Law of thermodynamics are satisfied. As a thermodynamical whole, the evolution of the Undulant Universe behaves very well in the current phase. However, when considering the unification theory, the failure of conversation law at the epoch of the matter dominated or near singularity need some more consideration for the form of the Undulant Universe.

We employ an analytical approach to investigate the signatures of Baryon Acoustic Oscillations (BAOs) on the convergence power spectrum of weak lensing by large scale structure. It is shown that the BAOs wiggles can be found in both of the linear and nonlinear convergence power spectra of weak lensing at about $40\le l\le600$, but they are weaker than that of matter power spectrum. Although the statistical error for LSST are greatly smaller than that of CFHT and SNAP survey especially at about $30<l<300$, they are still larger than the their maximum variations of BAOs wiggles. Thus, the detection of BAOs with the ongoing and upcoming surveys such as LSST, CFHT and SNAP survey confront a technical challenge.

We use the current weak lensing data to constrain the equation of state of dark energy $w$ and the total mass of massive neutrinos $\sum m_{\nu}$. The constraint on $w$ would be weak if only the current weak lensing data are used. With the addition of other observational data such as the type Ia supernovae, baryon acoustic oscillation, and the high redshift Hubble parameter data $H(z)$ derived from relative galaxy ages to break the degeneracy, the result is significantly improved. For the pure $w$CDM model without massive neutrinos, we find $w=-1.00^{+0.10}_{-0.12}$. For the $w$CDM model with the massive neutrino component, we show that the constraint on $w$ is almost unchanged, there is very little degeneracy between $w$ and $\sum m_{\nu}$. After marginalizing over other parameters, we obtain the probability distribution function of $\sum m_{\nu}$, and find that the upper limit is $\sum m_{\nu} \leq 0.8$ eV at 95.5% confidence level for the combined data sets. Our constraints of $w$ and $\sum m_{\nu}$ are both compatible and comparable with the constraints obtained from the WMAP 5-year data.