We present a complementarity study of gravitational waves and double Higgs production in the $4b$ channel, exploring the gauge singlet scalar extension of the SM. This new physics extension serves as a simplified benchmark model that realizes a strongly first-order electroweak phase transition necessary to generate the observed baryon asymmetry in the universe. In calculating the signal-to-noise ratio of the gravitational waves, we incorporate the effect of the recently discovered significant suppression of the gravitational wave signals from sound waves for strong phase transitions, make sure that supercooled phase transitions do complete and adopt a bubble wall velocity that is consistent with a successful electroweak baryogenesis by solving the velocity profiles of the plasma. The high-luminosity LHC sensitivity to the singlet scalar extension of the SM is estimated using a shape-based analysis of the invariant $4b$ mass distribution. We find that while the region of parameter space giving detectable gravitational waves is shrunk due to the new gravitational wave simulations, the qualitative complementary role of gravitational waves and collider searches remain unchanged.

We perform a systematic search for high-redshift ($z >$ 1.5) extreme variability quasars (EVQs) using repeat spectra from the Sixteenth Data Release of Sloan Digital Sky Survey, which provides a baseline spanning up to $\sim$18 yrs in the observed frame. We compile a sample of 348 EVQs with a maximum continuum variability at rest frame 1450 Angstrom of more than 100% (i.e., $\delta$V $\equiv$ (Max$-$Min)/Mean $>$1). Among those EVQs, we discover 23 new Changing-Look Active Galactic Nuclei (CLAGNs) with CL behaviors in different UV lines, e.g., CIV, CIII], SiIV, and MgII. Importantly, spurious CLAGNs caused by SDSS problematic spectral flux calibration, e.g., fiber drop issue, are rejected. Compared with normal EVQs, we find that CLAGNs have systematically lower continuum luminosity and broader emission line widths, hence lower Eddington ratios. Statistical tests of EVQ/CLAGN properties suggest that CLAGNs are likely to be a subset of EVQs with less efficient accretion rather than a distinctive AGN population. Moreover, we find that the difference composite spectra (i.e., bright$-$faint) of EVQs and CLAGNs are both similar to those of normal quasars with a power-law slope of $\alpha_{\rm \lambda,diff}\sim-2$, further supporting that they belong to the same population with analogous intrinsic variability mechanisms, e.g., rapid accretion rate change with a possible state transition. Finally, we also confirm the anti-breathing of CIV (i.e., line width increases as luminosity increases) in EVQs, and find that in addition to $\sim$ 0.4 dex systematic uncertainty in single-epoch CIV virial black hole mass estimates, an extra scatter of $\sim$ 0.3 dex will be introduced by extreme variability.

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$\beta$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{{\rm 2-12\,keV}}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$\beta$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$\sigma$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

We present the design and performance of broadband and tunable infrared-blocking filters for millimeter and sub-millimeter astronomy composed of small scattering particles embedded in an aerogel substrate. The ultra-low-density (typically < 150 mg/cm^3) aerogel substrate provides an index of refraction as low as 1.05, removing the need for anti-reflection coatings and allowing for broadband operation from DC to above 1 THz. The size distribution of the scattering particles can be tuned to provide a variable cutoff frequency. Aerogel filters with embedded high-resistivity silicon powder are being produced at 40-cm diameter to enable large-aperture cryogenic receivers for cosmic microwave background polarimeters, which require large arrays of sub-Kelvin detectors in their search for the signature of an inflationary gravitational-wave background.

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 report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

The 21~cm intensity mapping experiments promise to measure the baryon acoustic oscillations (BAO) through the mapping of H{\small I} gas at the post-reionization epoch. Its success depends on the assumption that H{\small I} gas traces the matter density distribution {\it linearly} on large scales. However, both nonlinear halo clustering and nonlinear effects modulating H{\small I} gas in halos may spoil this. We employ three approaches to generate the mock H{\small I} density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of H{\small I} linearity is valid at the scale corresponding to the first BAO peak, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and H{\small I} content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the H{\small I} bias is scale-independent down to small scales. We also find that the linear H{\small I} bias scales approximately linearly with redshift for $z\le 3$.

We present accurate measurements of the total HI mass in dark matter halos of different masses at z ~ 0, by stacking the HI spectra of entire groups from the Arecibo Fast Legacy ALFA Survey. The halos are selected from the optical galaxy group catalog constructed from the Sloan Digital Sky Survey DR7 Main Galaxy sample, with reliable measurements of halo mass and halo membership. We find that the HI-halo mass relation is not a simple monotonic function, as assumed in several theoretical models. In addition to the dependence of halo mass, the total HI gas mass shows strong dependence on the halo richness, with larger HI masses in groups with more members at fixed halo masses. Moreover, halos with at least three member galaxies in the group catalog have a sharp decrease of the HI mass, potentially caused by the virial halo shock-heating and the AGN feedback. The dominant contribution of the HI gas comes from the central galaxies for halos of $M<10^{12.5}h^{-1}M_{\odot}$, while the satellite galaxies dominate over more massive halos. Our measurements are consistent with a three-phase formation scenario of the HI-rich galaxies. The smooth cold gas accretion is driving the HI mass growth in halos of $M<10^{11.8}h^{-1}M_{\odot}$, with late-forming halos having more HI accreted. The virial halo shock-heating and AGN feedback will take effect to reduce the HI supply in halos of $10^{11.8}h^{-1}M_{\odot}<M<10^{13}h^{-1}M_{\odot}$. The HI mass in halos more massive than $10^{13}h^{-1}M_{\odot}$ generally grows by mergers, with the dependence on halo richness becoming much weaker.

Axion-like-particles (ALPs) emitted from the core of magnetars can convert to photons in the strong magnetic field of the magnetosphere. We study such emissions in the soft gamma-ray range from 300 keV to 1 MeV, where the ALP spectrum peaks and astrophysical backgrounds from resonant Compton upscattering are expected to be suppressed. Using published quiescent soft-gamma flux upper limits in 6 Magnetars obtained with $\textit{CGRO}$ COMPTEL, $\textit{INTEGRAL}$ SPI/IBIS/ISGRI and the $\textit{Fermi}$ Gamma Ray Burst Monitor (GBM), we put limits on the ALP-photon coupling obtained from conversions, assuming that ALP emission from the core is just sub-dominant to bounds from neutrino cooling. For core temperatures of $10^9$ K, the constraints on the ALP-photon coupling coming from magnetars 1E 2259+586 and J170849.0-400910 are better than the current limits obtained from the CAST experiment. We provide a detailed study of the dependence of our results on the magnetar core temperature. Our results motivate a program of studying quiescent soft gamma ray emission from magnetars in the 300 keV - 1 MeV band with the $\textit{Fermi}$-GBM.

We report the discovery of SDSS J0849+1114 as the first known triple Type 2 Seyfert nucleus. It represents three active black holes that are identified from new spatially resolved optical slit spectroscopy using the Dual Imaging Spectrograph on the 3.5 m telescope at the Apache Point Observatory. We also present new complementary observations including the Hubble Space Telescope Wide Field Camera 3 U- and Y-band imaging, Chandra Advanced CCD Imaging Spectrometer S-array X-ray 0.5--8 keV imaging spectroscopy, and NSF Karl G. Jansky Very Large Array radio 9.0 GHz imaging in its most extended A configuration. These comprehensive multiwavelength observations, when combined together, strongly suggest that all three nuclei are active galactic nuclei. While they are now still at kiloparsec-scale separations, where the host-galaxy gravitational potential dominates, the black holes may evolve into a bound triple system in $\lesssim$2 Gyr. These triple merger systems may explain the overly massive stellar cores that have been observed in some elliptical galaxies such as M87, which are expected to be unique gravitational wave sources. Similar systems may be more common in the early universe, when galaxy mergers are thought to have been more frequent.

Periodic quasars are candidates for binary supermassive black holes (BSBHs) efficiently emitting low frequency gravitational waves. Recently, $\sim$150 candidates were identified from optical synoptic surveys. However, they may be false positives caused by stochastic quasar variability given the few cycles covered (typically 1.5). To independently test the binary hypothesis, we search for evidence of truncated or gapped circumbinary accretion discs (CBDs) in their spectral energy distributions (SEDs). Our work is motivated by CBD simulations that predict flux deficits as cutoffs from central cavities opened by secondaries or notches from minidiscs around both BHs. We find that candidate periodic quasars show SEDs similar to those of control quasars matched in redshift and luminosity. While seven of 138 candidates show a blue cutoff in the IR-optical-UV SED, six of which may represent CBDs with central cavities, the red SED fraction is similar to that in control quasars, suggesting no correlation between periodicity and SED anomaly. Alternatively, dust reddening may cause red SEDs. The fraction of extremely radio-loud quasars, e.g., blazars (with $R>100$), is tentatively higher than that in control quasars (at 2.5$\sigma$). Our results suggest that, assuming most periodic candidates are robust, IR-optical-UV SEDs of CBDs are similar to those of accretion discs of single BHs, if the periodicity is driven by BSBHs; the higher blazar fraction may signal precessing radio jets. Alternatively, most current candidate periodic quasars identified from few-cycle light curves may be false positives. Their tentatively higher blazar fraction and lower Eddington ratios may both be caused by selection biases.

The broad Mg II line in quasars has distinct variability properties compared with broad Balmer lines: it is less variable, and usually does not display a "breathing" mode, the increase in the average cloud distance when luminosity increases. We demonstrate that these variability properties of Mg II can be reasonably well explained by simple Locally Optimally Emitting Cloud (LOC) photoionization models, confirming earlier photoionization results. In the fiducial LOC model, the Mg II-emitting gas is on average more distant from the ionizing source than the H$\alpha$/H$\beta$ gas, and responds with a lower amplitude to continuum variations. If the broad-line region (BLR) is truncated at a physical radius of $\sim 0.3$ pc (for a $10^{8.5}M_{\odot}$ BH accreting at Eddington ratio of 0.1), most of the Mg II flux will always be emitted near this outer boundary and hence will not display breathing. These results indicate that reverberation mapping results on broad Mg II, while generally more difficult to obtain due to the lower line responsivity, can still be used to infer the Mg II BLR size and hence black hole mass. But it is possible that Mg II does not have a well defined intrinsic BLR size-luminosity relation for individual quasars, even though a global one for the general population may still exist. The dramatic changes in broad H$\alpha$/H$\beta$ emission in the observationally-rare changing-look quasars are fully consistent with photoionization responses to extreme continuum variability, and the LOC model provides natural explanations for the persistence of broad Mg II in changing-look quasars defined on H$\alpha$/H$\beta$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies.

Dark matter exists in our Universe, but its nature remains mysterious. The remarkable sensitivity of the Laser Interferometer Gravitational-Wave Observatory (LIGO) may be able to solve this mystery. A good dark matter candidate is the ultralight dark photon. Because of its interaction with ordinary matter, it induces displacements on LIGO mirrors that can lead to an observable signal. In a study that bridges gravitational wave science and particle physics, we perform a direct dark matter search using data from LIGO's first (O1) data run, as opposed to an indirect search for dark matter via its production of gravitational waves. We demonstrate an achieved sensitivity on squared coupling as $\sim 4 \times 10^{-45}$, in a $U(1)_{\rm B}$ dark photon dark matter mass band around $m_{\rm A} \sim 4 \times 10^{-13}$eV. Substantially improved search sensitivity is expected during the coming years of continued data taking by LIGO and other gravitational wave detectors in a growing global network.

Changing-Look (CL) is a rare phenomenon of Active Galactic Nuclei (AGNs) that exhibit emerging or disappearing broad lines accompanied by continuum variations on astrophysically short timescales ($\lesssim$ 1 yr to a few decades). While previous studies have found Balmer-line (broad H$\alpha$ and/or H$\beta$) CL AGNs, the broad Mg II line is persistent even in dim states. No unambiguous Mg II CL AGN has been reported to date. We perform a systematic search of Mg II CL AGNs using multi-epoch spectra of a special population of Mg II-emitters (characterized by strong broad Mg II emission with little evidence for AGN from other normal indicators such as broad H$\alpha$ and H$\beta$ or blue power-law continua) from the Fourteenth Data Release of the Sloan Digital Sky Survey. We present the discovery of the first unambiguous case of an Mg II CL AGN, SDSS J152533.60+292012.1 (at redshift $z$ = 0.449), which is turning off within rest-frame 286 days. The dramatic diminishing of Mg II equivalent width (from 110 $\pm$ 26 Angstrom to being consistent with zero), together with little optical continuum variation ($\Delta V_{\rm max-min}$ $=$ 0.17 $\pm$ 0.05 mag) coevally over $\sim$ 10 years, rules out dust extinction or a tidal disruption event. Combined with previously known H$\beta$ CL AGNs, we construct a sequence that represents different temporal stages of CL AGNs. This CL sequence is best explained by the photoionization model of Guo et al. (2019). In addition, we present two candidate turn-on Mg II CL AGNs and a sample of 361 Mg II-emitters for future Mg II CL AGN searches.

We present the relation between the sphaleron energy and the gravitational wave signals from a first order electroweak phase transition. The crucial ingredient is the scaling law between the sphaleron energy at the temperature of the phase transition and that at zero temperature. We estimate the baryon number preservation criterion, and observe that for a sufficiently strong phase transition, it is possible to probe the electroweak sphaleron using measurements of future space-based gravitational wave detectors.

We present Chandra ACIS-S X-ray imaging spectroscopy for five dual active galactic nucleus (AGN) candidates. Our targets were drawn from a sample of 1286 [O III]-selected AGN pairs systematically selected from the Sloan Digital Sky Survey Seventh Data Release. Each of the targets contains two nuclei separated by $\sim$3-9 kpc in projection, both of which are optically classified as Type 2 (obscured) AGNs based on diagnostic ratios of the narrow emission lines. Combined with independent, empirical star formation rate estimates based on the host-galaxy stellar continua, the new Chandra X-ray observations allow us to evaluate the dual-AGN hypothesis for each merging system. We confirm two (SDSS J0907+5203 and SDSS J1544+0446) of the five targets as bona-fide dual AGNs. For the other three targets, the existing data are consistent with the dual-AGN scenario, but we cannot rule out the possibility of stellar/shock heating and/or one AGN ionizing both gaseous components in the merger. The average X-ray-to-[O III] luminosity ratio in our targets seems to be systematically smaller than that observed in single AGNs but is higher than that seen in dual AGNs selected from AGNs with double-peaked narrow emission lines. We suggest that the systematically smaller X-ray-to-[O III] luminosity ratio observed in dual AGNs than in single AGNs is due to a high nuclear gas column likely from strong merger-induced inflows. Unlike double-peaked-[O III]-selected dual AGNs, the new sample selected from resolved galaxy pairs are not subject to the orientation bias caused by the double-peak line-of-sight velocity splitting selection, which also contributes to lowering the X-ray-to-[O III] luminosity ratio.

Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required.

We present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution.

We present a study of a sample of 223 radio-loud quasars (up to redshift $<$0.3) in order to investigate their spectral properties. Twenty-six of these radio-loud quasars are identified as Flat Spectrum Radio Quasars (FSRQs) and fifty-four are identified as Steep Spectrum Radio Quasars (SSRQs) based on their radio spectral index. We study the [O III] line properties of these quasars to investigate the origin and properties of blue wings (shift of the profile towards lower wavelengths) and blue outliers (shift of the whole spectroscopic feature). Most of the quasars show blue wings with velocities up to 420 km $s^{-1}$. We find that around 17\% of the quasars show outliers with velocities spanning 419 to -315 km $s^{-1}$. Finally, we revisit the $\it M_{\rm BH} - \sigma$ relation of our sample using [S II]$\lambda$6716, 6731 and [O III] linewidths as surrogates for stellar velocity dispersions, $\sigma$, to investigate their location on the $\it M_{\rm BH} - \sigma$ relation for quiescent galaxies. Due to the strong blending of [S II] with $\rm H_{\alpha}$, we could estimate $\sigma_{[\rm SII]}$ of only 123 quasars. We find that the radio-loud quasars do not show a relationship between $\it M_{\rm BH}$ and $\sigma_{\rm [SII]/[OIII]}$ up to a redshift of 0.3, although they cluster around the local relation. We find an overall offset of 0.12$\pm$0.05 dex of our sample of radio-loud quasars from the $\it M_{\rm BH} - \sigma$ relation of quiescent galaxies. Quasars in our highest redshift bin (z=0.25-0.3) show a deviation of $\sim$0.33 $\pm$ 0.06 dex with respect to the local relation.

The large-scale redshift-space distortion (RSD) in galaxy clustering can probe $f\sigma_8$, a combination of the cosmic structure growth rate and the matter fluctuation amplitude, which can constrain dark energy models and test theories of gravity. While the RSD on small scales (e.g. a few to tens of $h^{-1}{\rm Mpc}$) can further tighten the $f\sigma_8$ constraints, galaxy assembly bias, if not correctly modelled, may introduce systematic uncertainties. Using a mock galaxy catalogue with built-in assembly bias, we perform a preliminary study on how assembly bias may affect the $f\sigma_8$ inference. We find good agreement on scales down to 8--9$h^{-1}{\rm Mpc}$ between a $f\sigma_8$ metric from the redshift-space two-point correlation function with the central-only mock catalogue and that with the shuffled catalogue free of assembly bias, implying that $f\sigma_8$ information can be extracted on such scales even with assembly bias. We then apply the halo occupation distribution (HOD) and three subhalo clustering and abundance matching (SCAM) models to model the redshift-space clustering with the mock. Only the SCAM model based on $V_{\rm peak}$ (used to create the mock) can reproduce the $f\sigma_8$ metric, and the other three could not. However, the $f\sigma_8$ metrics determined from central galaxies from all the models are able to match the expected one down to 8$h^{-1}{\rm Mpc}$. Our results suggest that halo models with no or incorrect assembly bias prescription could still be used to model the RSD down to scales of $\sim 8 h^{-1}{\rm Mpc}$ to tighten the $f\sigma_8$ constraint, with a sample of central galaxies or with a flexible satellite occupation prescription.

We investigate the abundance, small-scale clustering and galaxy-galaxy lensing signal of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). To this end, we present new measurements of the redshift and stellar mass dependence of the lensing properties of the galaxy sample. We analyse to what extent models assuming the Planck18 cosmology fit to the number density and clustering can accurately predict the small-scale lensing signal. In qualitative agreement with previous BOSS studies at redshift $z \sim 0.5$ and with results from the Sloan Digital Sky Survey, we find that the expected signal at small scales ($0.1 < r_{\rm p} < 3 \, h^{-1} \mathrm{Mpc}$) is higher by $\sim 25\%$ than what is measured. Here, we show that this result is persistent over the redshift range $0.1 < z < 0.7$ and for galaxies of different stellar masses. If interpreted as evidence for cosmological parameters different from the Planck CMB findings, our results imply $S_8 = \sigma_8 \sqrt{\Omega_{\rm m} / 0.3} = 0.744 \pm 0.015$, whereas $S_8 = 0.832 \pm 0.013$ for Planck18. However, in addition to being in tension with CMB results, such a change in cosmology alone does not accurately predict the lensing amplitude at larger scales. Instead, other often neglected systematics like baryonic feedback or assembly bias are likely contributing to the small-scale lensing discrepancy. We show that either effect alone, though, is unlikely to completely resolve the tension. Ultimately, a combination of the two effects in combination with a moderate change in cosmological parameters might be needed.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

Changing look active-galatic-nuclei (CL AGNs) can yield considerable insight into accretion physics as well as the co-evolution of black holes and their host galaxies. A large sample of these CL AGNs is essential to achieve the latter goal. We propose an effective method to select CL candidates from spectroscopic quasar catalogs using the mid-infrared (MIR) variability information derived from ALL-WISE/NEOWISE data releases. Our primary selection criteria include both a large amplitude flux variation and a transition of MIR color from an AGN to a normal galaxy. A pilot spectroscopic follow-up of 7 candidates among about 300 candidates selected from Sloan Digital Sky Survey low-redshift (z<0.5) AGN sample results in 6 new turn-off CL AGNs. We expect to obtain hundreds of CL AGNs once full spectroscopic follow-up of the sample is carried out.

A class of dark photon dark matter models with ultralight masses would lead to oscillation of a test body through a coupling with baryons or $B-L$ charge. This periodical oscillation of an observer results in swing of a star's apparent position due to the effect of aberration of light, which could be probed with high-precision astrometry observations of stars in the Milky Way. We propose to use the observations of stellar positions of a number of stars by Gaia to search for this kind of dark photon dark matter. We show that this astrometry method is able to give promising sensitivities to search for the dark photon dark matter in the mass range of $10^{-23}\sim10^{-21}$ eV.

We propose to use gravitational waves from extreme mass ratio inspirals (EMRI), composed of a boson star and a supermassive black hole in the center of galaxies, as a new method to search for boson stars. Gravitational waves from EMRI have the advantage of being long-lasting within the frequency band of future space-based interferometer gravitational wave detectors and can accumulate large signal-to-noise ratio (SNR) for very sub-solar mass boson stars. Compared to gravitational waves from boson star binaries, which fall within the LIGO band, we find that much larger ranges of the mass and compactness of boson stars, as well as the underlying particle physics parameter space, can be probed by EMRI. We take tidal disruption of the boson stars into account and distinguish those which dissolve before the inner-most-stable-circular-orbit (ISCO) and those which dissolve after it. Due to the large number of cycles recorded, EMRIs can lead to a very precise mass determination of the boson star and distinguish it from standard astrophysical compact objects in the event of a signal. Tidal effects in inspiralling binary systems, as well as possible correlated electromagnetic signals, can also serve as potential discriminants.

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1<z<4.5 and is flux-limited to i_PSF<21.7 without any other cuts on quasar properties. The main sample characterization includes: 1) spectral measurements of the continuum and broad emission lines for individual objects from the coadded first-season spectroscopy in 2014; 2) identification of broad and narrow absorption lines in the spectra; 3) optical variability properties for continuum and broad lines from multi-epoch spectroscopy. We provide improved systemic redshift estimates for all quasars, and demonstrate the effects of signal-to-noise ratio on the spectral measurements. We compile measured properties for all 849 quasars along with supplemental multi-wavelength data for subsets of our sample from other surveys. The SDSS-RM sample probes a diverse range in quasar properties, and shows well detected continuum and broad-line variability for many objects from first-season monitoring data. The compiled properties serve as the benchmark for follow-up work based on SDSS-RM data. The spectral fitting tools are made public along with this work.

We present a dedicated complementarity study of gravitational wave and collider measurements of the simplest extension of the Higgs sector: the singlet scalar augmented Standard Model. We study the following issues: (i) the electroweak phase transition patterns admitted by the model, and the proportion of parameter space for each pattern; (ii) the regions of parameter space that give detectable gravitational waves at future space-based detectors; and (iii) the current and future collider measurements of di-Higgs production, as well as searches for a heavy weak diboson resonance, and how these searches interplay with regions of parameter space that exhibit strong gravitational wave signals. We carefully investigate the behavior of the normalized energy released during the phase transition as a function of the model parameters, address subtle issues pertaining to the bubble wall velocity, and provide a description of different fluid velocity profiles. On the collider side, we identify the subset of points that are most promising in terms of di-Higgs and weak diboson production studies while also giving detectable signals at LISA, setting the stage for future benchmark points that can be used by both communities.

We study the evolution of star-forming galaxies with $10^{10} M_\odot<M_*<10^{11.6} M_\odot$ over the redshift range of 0.7<z<1.2 using the emission line galaxies (ELGs) in the extended Baryon Oscillation Spectroscopic Survey (eBOSS). By applying the incomplete conditional stellar mass function (ICSMF) model proposed in Guo et al., we simultaneously constrain the sample completeness, the stellar--halo mass relation (SHMR) and the quenched galaxy fraction. We obtain the intrinsic stellar mass functions for star-forming galaxies in the redshift bins of 0.7<z<0.8, 0.8<z<0.9, 0.9<z<1.0 and 1.0<z<1.2, as well as the stellar mass function for all galaxies in the redshift bin of 0.7<z<0.8. We find that the eBOSS ELG sample only selects about 1%-10% of the star-forming galaxy population at the different redshifts, with the lower redshift samples to be more complete. There is only weak evolution in the SHMR of the ELGs from z=1.2 to z=0.7, as well as the intrinsic galaxy stellar mass functions for lower-mass galaxies of $M_*<10^{11} M_\odot$. There is significant decrease of the stellar mass function for star-forming galaxies with redshift at the massive end. Our best-fitting models show that the central ELGs at these redshifts live in halos of mass $M\sim10^{12} M_\odot$ while the satellite ELGs occupy slightly more massive halos of $M\sim10^{12.6} M_\odot$. The average satellite fraction of the observed ELGs varies from 13% to 17%, with the galaxy bias increasing from 1.1 to 1.4 from z=0.7 to 1.2.

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.

Quasars whose broad emission lines show temporal, bulk radial velocity (RV) shifts have been proposed as candidate sub-parsec (sub-pc), binary supermassive black holes (BSBHs). We identified a sample of 16 BSBH candidates based on two-epoch spectroscopy among 52 quasars with significant RV shifts over a few rest-frame years. The candidates showed consistent velocity shifts independently measured from two broad lines (H$\beta$ and H$\alpha$ or Mg${\rm \,II}$) without significant changes in the broad-line profiles. Here in the third paper of the series, we present further third- and fourth-epoch spectroscopy for 12 of the 16 candidates for continued RV tests, spanning $\sim$5--15 yr in the quasars' rest frames. Cross-correlation analysis of the broad H$\beta$ calibrated against [O${\rm\,III}]\,\lambda 5007$ suggests that 5 of the 12 quasars remain valid as BSBH candidates. They show broad H$\beta$ RV curves that are consistent with binary orbital motion without significant changes in the broad line profiles. Their broad H$\alpha$ (or Mg${\rm \,II}$) lines display RV shifts that are either consistent with or smaller than those seen in broad H$\beta$. The RV shifts can be explained by a $\sim$0.05--0.1 pc BSBH with an orbital period of $\sim$40--130 yr, assuming a mass ratio of 0.5--2 and a circular orbit. However, the parameters are not well constrained given the few epochs that sample only a small portion of the hypothesized binary orbital cycle. The apparent occurrence rate of sub-pc BSBHs is $\lesssim$13$\pm$5% among all SDSS quasars, with no significant difference in the subsets with and without single-epoch broad line velocity offsets. Dedicated long-term spectroscopic monitoring is still needed to further confirm or reject these BSBH candidates.

Observations have shown that UV/optical variation amplitude of quasars depend on several physi- cal parameters including luminosity, Eddington ratio, and likely also black hole mass. Identifying new factors which correlate with the variation is essential to probe the underlying physical processes. Combining ~ten years long quasar light curves from SDSS stripe 82 and X-ray data from Stripe 82X, we build a sample of X-ray detected quasars to investigate the relation between UV/optical variation amplitude ($\sigma_{rms}$) and X-ray loudness. We find that quasars with more intense X-ray radiation (com- pared to bolometric luminosity) are more variable in UV/optical. Such correlation remains highly significant after excluding the effect of other parameters including luminosity, black hole mass, Ed- dington ratio, redshift, rest-frame wavelength (i.e., through partial correlation analyses). We further find the intrinsic link between X-ray loudness and UV/optical variation is gradually more prominent on longer timescales (up to 10 years in the observed frame), but tends to disappear at timescales < 100 days. This suggests a slow and long-term underlying physical process. The X-ray reprocessing paradigm, in which UV/optical variation is produced by a variable central X-ray emission illuminating the accretion disk, is thus disfavored. The discovery points to an interesting scheme that both the X-ray corona heating and UV/optical variation is quasars are closely associated with magnetic disc turbulence, and the innermost disc turbulence (where corona heating occurs) correlates with the slow turbulence at larger radii (where UV/optical emission is produced).

We formulate a model of the conditional colour-magnitude distribution (CCMD) to describe the distribution of galaxy luminosity and colour as a function of halo mass. It consists of two populations of different colour distributions, dubbed pseudo-blue and pseudo-red, respectively, with each further separated into central and satellite galaxies. We define a global parameterization of these four colour-magnitude distributions and their dependence on halo mass, and we infer parameter values by simultaneously fitting the space densities and auto-correlation functions of 79 galaxy samples from the Sloan Digital Sky Survey defined by fine bins in the colour-magnitude diagram (CMD). The model deprojects the overall galaxy CMD, revealing its tomograph along the halo mass direction. The bimodality of the colour distribution is driven by central galaxies at most luminosities, though at low luminosities it is driven by the difference between blue centrals and red satellites. For central galaxies, the two pseudo-colour components are distinct and orthogonal to each other in the CCMD: at fixed halo mass, pseudo-blue galaxies have a narrow luminosity range and broad colour range, while pseudo-red galaxies have a narrow colour range and broad luminosity range. For pseudo-blue centrals, luminosity correlates tightly with halo mass, while for pseudo-red galaxies colour correlates more tightly (redder galaxies in more massive haloes). The satellite fraction is higher for redder and for fainter galaxies, with colour a stronger indicator than luminosity. We discuss the implications of the results and further applications of the CCMD model.

We explore the evolution of the time variability (in the optical $g$-band and on timescales of weeks to years) of SDSS Stripe 82 quasars along the quasar main sequence. A parent sample of $1004$ quasars within $0.5\leq z \leq 0.89$ are used for our statistical studies, we then make subsamples from our parent sample: a subsample of $246$ quasars with similar luminosities, and a subsample of $399$ quasars with similar Rfe (i.e., the ratio of the equivalent width of FeII within $4435$--$4685\ \mathrm{\AA}$ to that of \Hbeta). We find the variability amplitude decreases with luminosity ($L_{\mathrm{bol}}$). The anti-correlation between the variability amplitude and Rfe is weak but statistically significant. The characteristic timescale, $\tau$, correlates mostly with quasar luminosity, its dependence on Rfe is statistically insignificant. After controlling luminosity and Rfe, the high- and low-FWHM samples have similar structure functions. These results support the framework that Rfe is governed by Eddington ratio and FWHM of Hbeta is mostly determined by orientation. We then provide new empirical relations between variability parameters and quasar properties (i.e., luminosity and Rfe). Our new relations are consistent with the scenario that quasar variability is driven by the thermal fluctuations in the accretion disk, $\tau$ seems to correspond to the thermal timescale. From our new relations, we find the short-term variability is mostly sensitive to $L_{\mathrm{bol}}$. Basing on this, we propose that quasar short-term (a few months) variability might be a new type of "Standard Candle" and can be adopted to probe cosmology.

Kiloparsec-scale dual active galactic nuclei (AGNs) are active supermassive black hole pairs co-rotating in galaxies with separations of less than a few kpc. Expected to be a generic outcome of hierarchical galaxy formation, their frequency and demographics remain uncertain. We have carried out an imaging survey with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) of AGNs with double-peaked narrow [O III] emission lines. HST/WFC3 offers high image quality in the near-infrared (NIR) to resolve the two stellar nuclei, and in the optical to resolve [O III] from ionized gas in the narrow-line regions. This combination has proven to be key in sorting out alternative scenarios. With HST/WFC3 we are able to explore a new population of close dual AGNs at more advanced merger stages than can be probed from the ground. Here we show that the AGN SDSS J0924+0510, which had previously shown two stellar bulges, contains two spatially distinct [O III] regions consistent with a dual AGN. While we cannot completely exclude cross-ionization from a single central engine, the nearly equal ratios of [O III] strongly suggest a dual AGN with a projected angular separation of 0."4, corresponding to a projected physical separation of $r_p$ = 1 kpc at redshift z = 0.1495. This serves as a proof of principle for combining high-resolution NIR and optical imaging to identify close dual AGNs. Our result suggests that studies based on low-resolution and/or low-sensitivity observations may miss close dual AGNs and thereby may underestimate their occurrence rate on $\lesssim$ kpc scales.

We propose a novel method to constrain the missing fraction of galaxies using galaxy clustering measurements in the galaxy conditional stellar mass function (CSMF) framework, which is applicable to surveys that suffer significantly from sample selection effects. The clustering measurements, which are not sensitive to the random sampling (missing fraction) of galaxies, are widely used to constrain the stellar-halo mass relation (SHMR). By incorporating a missing fraction (incompleteness) component into the CSMF model (ICSMF), we use the incomplete stellar mass function and galaxy clustering to simultaneously constrain the missing fractions and the SHMRs. Tests based on mock galaxy catalogs with a few typical missing fraction models show that this method can accurately recover the missing fraction and the galaxy SHMR, and hence provides us reliable measurements of the galaxy stellar mass functions. We then apply it to the Baryon Oscillation Spectroscopic Survey (BOSS) over the redshift range of 0.1<z<0.8 for galaxies of $M_*>10^{11}M_\odot$. We find the sample completeness for BOSS is over 80% at z<0.6, but decreases at higher redshifts to about 30%. After taking these completeness factors into account, we provide accurate measurements of the stellar mass functions for galaxies with $10^{11}M_\odot<M_*<10^{12}M_\odot$, as well as the SHMRs, over the redshift range 0.1<z<0.8 in this largest galaxy redshift survey.

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

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

NGC 5548 has been intensively monitored by the AGN Space Telescope and Optical Reverberation Mapping collaboration. Approximately after half of the light curves, the correlation between the broad emission lines and the lag-corrected ultraviolet continua becomes weak. This anomalous behavior is accompanied by an increase of soft X-ray emission. We propose a simple model to understand this anomalous behavior, i.e., the corona might fall down, thereby increasing the covering fraction of the inner disk. Therefore, X-ray and extreme ultraviolet emission suffer from spectral variations. The ultraviolet continua variations are driven by both X-ray and extreme ultraviolet variations. Consequently, the spectral variability induced by the falling corona would dilute the correlation between the broad emission lines and the ultraviolet continua. Our model can explain many additional observational facts, including the dependence of the anomalous behavior on velocity and ionization energy. We also show that the time lag and correlation between the X-ray and the ultraviolet variations change as NGC 5548 displays the anomalous behavior. The time lag is dramatically longer than the expectation from disk reprocessing if the anomalous behavior is properly excluded. During the anomalous state, the time lag approaches the light-travel timescale of disk reprocessing albeit with a much weaker correlation. We speculate that the time lag in the normal state is caused by reprocessing of the broad line region gas. As NGC 5548 enters the abnormal state, the contribution of the broad line region gas is smaller; the time lag reflects disk reprocessing. We also discuss alternative scenarios.

Because of vast volume of data being produced by today's scientific simulations and experiments, lossy data compressor allowing user-controlled loss of accuracy during the compression is a relevant solution for significantly reducing the data size. However, lossy compressor developers and users are missing a tool to explore the features of scientific datasets and understand the data alteration after compression in a systematic and reliable way. To address this gap, we have designed and implemented a generic framework called Z-checker. On the one hand, Z-checker combines a battery of data analysis components for data compression. On the other hand, Z-checker is implemented as an open-source community tool to which users and developers can contribute and add new analysis components based on their additional analysis demands. In this paper, we present a survey of existing lossy compressors. Then we describe the design framework of Z-checker, in which we integrated evaluation metrics proposed in prior work as well as other analysis tools. Specifically, for lossy compressor developers, Z-checker can be used to characterize critical properties of any dataset to improve compression strategies. For lossy compression users, Z-checker can detect the compression quality, provide various global distortion analysis comparing the original data with the decompressed data and statistical analysis of the compression error. Z-checker can perform the analysis with either coarse granularity or fine granularity, such that the users and developers can select the best-fit, adaptive compressors for different parts of the dataset. Z-checker features a visualization interface displaying all analysis results in addition to some basic views of the datasets such as time series. To the best of our knowledge, Z-checker is the first tool designed to assess lossy compression comprehensively for scientific datasets.

Multi-hydrogenated species with proper symmetry properties can present different spin configurations, and thus exist under different spin symmetry forms, labeled as para and ortho for two-hydrogen molecules. We investigated here the ortho-to-para ratio (OPR) of H$_2$Cl$^+$ in the light of new observations performed in the z=0.89 absorber toward the lensed quasar PKS 1830-211 with the Atacama Large Millimeter/submillimeter Array (ALMA). Two independent lines of sight were observed, to the southwest (SW) and northeast (NE) images of the quasar, with OPR values found to be $3.15 \pm 0.13$ and $3.1 \pm 0.5$ in each region, respectively, in agreement with a spin statistical weight of 3:1. An OPR of 3:1 for a molecule containing two identical hydrogen nuclei can refer to either a statistical result or a high-temperature limit depending on the reaction mechanism leading to its formation. It is thus crucial to identify rigorously how OPRs are produced in order to constrain the information that these probes can provide. To understand the production of the H$_2$Cl$^+$ OPR, we undertook a careful theoretical study of the reaction mechanisms involved with the aid of quasi-classical trajectory calculations on a new global potential energy surface fit to a large number of high-level ab initio data. Our study shows that the major formation reaction for H$_2$Cl$^+$ produces this ion via a hydrogen abstraction rather than a scrambling mechanism. Such a mechanism leads to a 3:1 OPR, which is not changed by destruction and possible thermalization reactions for H$_2$Cl$^+$ and is thus likely to be the cause of observed 3:1 OPR ratios, contrary to the normal assumption of scrambling.

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.

Studies have shown that UV/optical light curves of quasars can be described with the prevalent damped random walk (DRW, also known as Ornstein-Uhlenbeck process) model. A white noise power spectral density (PSD) is expected at low frequency in this model, however, direct observational constraint to the low frequency PSD slope is hard due to limited lengths of the light curves available. Meanwhile, quasars show too large scatter in their DRW parameters to be attributed to the uncertainties in the measurements and the dependence of variation to known physical factors. In this work we present simulations showing that, if the low frequency PSD deviates from DRW, the red noise leakage can naturally produce large scatter in variation parameters measured from simulated light curves. The steeper the low frequency PSD slope is, the larger scatter we expect. Based on the observations of SDSS Stripe 82 quasars, we find the low frequency PSD slope should be no steeper than -1.3. The actual slope could be flatter, which consequently requires that quasar variabilities should be influenced by other unknown factors. We speculate that magnetic field and/or metallicity could be such additional factors.

Primordial black hole (PBH), which can be naturally produced in the early universe, remains a promising dark matter candidate . It can merge with a supermassive black hole (SMBH) in the center of a galaxy and generate gravitational wave (GW) signals in the favored frequency region of LISA-like experiments. In this work, we initiate the study on the event rate calculation for such extreme mass ratio inspirals (EMRI). Including the sensitivities of various proposed GW detectors, we find that such experiments offer a novel and outstanding tool to test the scenario where PBH constitutes (fraction of) dark matter. The PBH energy density fraction of DM ($f_\text{PBH}$) could potentially be explored as small as $10^{-3} \sim 10^{-4}$. Further, LISA has the capability to search for PBH mass upto $10^{-2} \sim 10^{-1} M_\odot$. Other proposed GW experiments can probe lower PBH mass regime.

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

We study the galaxy clustering dependence on the $\left[\mathrm{O\scriptsize{II}}\right]$ emission line luminosity in the SDSS DR7 Main galaxy sample at mean redshift $z\sim0.1$. We select volume-limited samples of galaxies with different $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity thresholds and measure their projected, monopole and quadrupole two-point correlation functions. We model these observations using the 1$h^{-1}\rm{Gpc}$ MultiDark Planck cosmological simulation and generate light-cones with the SUrvey GenerAtoR algorithm. To interpret our results, we adopt a modified (Sub)Halo Abundance Matching scheme, accounting for the stellar mass incompleteness of the emission line galaxies. The satellite fraction constitutes an extra parameter in this model and allows to optimize the clustering fit on both small and intermediate scales (i.e. $r_p\lesssim 30h^{-1}\rm{Mpc})$, with no need of any velocity bias correction. We find that, in the local Universe, the $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity correlates with all the clustering statistics explored and with the galaxy bias. This latter quantity correlates more strongly with the SDSS $r$-band magnitude than $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity. In conclusion, we propose a straightforward method to produce reliable clustering models, entirely built on the simulation products, which provides robust predictions of the typical ELG host halo masses and satellite fraction values. The SDSS galaxy data, MultiDark mock catalogues and clustering results are made publicly available.

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

We use the hydrodynamical galaxy formation simulations from the Illustris suite to study the origin and properties of galaxy velocity bias, i.e., the difference between the velocity distributions of galaxies and dark matter inside halos. We find that galaxy velocity bias is a decreasing function of the ratio of galaxy stellar mass to host halo mass. In general, central galaxies are not at rest with respect to dark matter halos or the core of halos, with a velocity dispersion above 0.04 times that of the dark matter. The central galaxy velocity bias is found to be mostly caused by the close interactions between the central and satellite galaxies. For satellite galaxies, the velocity bias is related to their dynamical and tidal evolution history after being accreted onto the host halos. It depends on the time after the accretion and their distances from the halo centers, with massive satellites generally moving more slowly than the dark matter. The results are in broad agreements with those inferred from modeling small-scale redshift-space galaxy clustering data, and the study can help improve models of redshift-space galaxy clustering.

We study in this work a scenario that the universe undergoes a two step phase transition with the first step happened to the dark matter sector and the second step being the transition between the dark matter and the electroweak vacuums, where the barrier between the two vacuums, that is necessary for a strongly first order electroweak phase transition (EWPT) as required by the electroweak baryogenesis mechanism, arises at the tree-level. We illustrate this idea by working with the standard model (SM) augmented by a scalar singlet dark matter and an extra scalar singlet which mixes with the SM Higgs boson. We study the conditions for such pattern of phase transition to occur and especially for the strongly first order EWPT to take place, as well as its compatibility with the basic requirements of a successful dark matter, such as observed relic density and constraints of direct detections. We further explore the discovery possibility of this pattern EWPT by searching for the gravitational waves generated during this process in spaced based interferometer, by showing a representative benchmark point of the parameter space that the generated gravitational waves fall within the sensitivity of eLISA, DECIGO and BBO.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

Based on recent $Herschel$ results, the ortho-to-para ratio (OPR) of NH$_2$ has been measured towards the following high-mass star-forming regions: W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4), and G34.3+0.1. The OPR at thermal equilibrium ranges from the statistical limit of three at high temperatures to infinity as the temperature tends toward zero, unlike the case of H$_{2}$. Depending on the position observed along the lines-of-sight, the OPR was found to lie either slightly below the high temperature limit of three (in the range $2.2-2.9$) or above this limit ($\sim3.5$, $\gtrsim 4.2$, and $\gtrsim 5.0$). In low temperature interstellar gas, where the H$_{2}$ is para-enriched, our nearly pure gas-phase astrochemical models with nuclear-spin chemistry can account for anomalously low observed NH$_2$-OPR values. We have tentatively explained OPR values larger than three by assuming that spin thermalization of NH$_2$ can proceed at least partially by H-atom exchange collisions with atomic hydrogen, thus increasing the OPR with decreasing temperature. In this paper, we present quasi-classical trajectory calculations of the H-exchange reaction NH$_2$ + H, which show the reaction to proceed without a barrier, confirming that the H-exchange will be efficient in the temperature range of interest. With the inclusion of this process, our models suggest both that OPR values below three arise in regions with temperatures $\gtrsim20-25$~K, depending on time, and values above three but lower than the thermal limit arise at still lower temperatures.

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

We present measurements of the clustering of galaxies as a function of their stellar mass in the Baryon Oscillation Spectroscopic Survey. We compare the clustering of samples using 12 different methods for estimating stellar mass, isolating the method that has the smallest scatter at fixed halo mass. In this test, the stellar mass estimate with the smallest errors yields the highest amplitude of clustering at fixed number density. We find that the PCA stellar masses of Chen etal (2012) clearly have the tightest correlation with halo mass. The PCA masses use the full galaxy spectrum, differentiating them from other estimates that only use optical photometric information. Using the PCA masses, we measure the large-scale bias as a function of Mgal for galaxies with logMgal>=11.4, correcting for incompleteness at the low-mass end of our measurements. Using the abundance-matching ansatz to connect dark matter halo mass to stellar mass, we construct theoretical models of b(Mgal) that match the same stellar mass function but have different amounts of scatter in stellar mass at fixed halo mass, sigma_logM. Using this approach, we find sigma_logM=0.18^{+0.01}_{-0.02}. This value includes both intrinsic scatter as well as random errors in the stellar masses. To partially remove the latter, we use repeated spectra to estimate statistical errors on the stellar masses, yielding an upper limit to the intrinsic scatter of 0.16 dex.

Faint red galaxies in the Sloan Digital Sky Survey show a puzzling clustering pattern in previous measurements. In the two-point correlation function (2PCF), they appear to be strongly clustered on small-scales, indicating a tendency to reside in massive haloes as satellite galaxies. However, their weak clustering on large scales suggests that they are more likely to be found in low mass haloes. The interpretation of the clustering pattern suffers from the large sample variance in the 2PCF measurements, given the small volume of the volume-limited sample of such faint galaxies. We introduce a method to improve the clustering measurements of faint galaxies by making a full use of a flux-limited sample to obtain volume-limited measurements with an increased effective volume. In the improved 2PCF measurements, the fractional uncertainties on large-scales drop by more than 40 per cent, and the strong contrast between small-scale and large-scale clustering amplitudes seen in previous work is no longer prominent. From halo occupation distribution modelling of the measurements, we find that a considerable fraction of faint red galaxies to be satellites in massive haloes, a senario supported by the strong covariance of small-scale 2PCF measurements and the relative spatial distribution of faint red galaxies and luminous galaxies. However, the satellite fraction is found to be degenerate with the slope of the distribution profile of satellites in inner haloes. We compare the modelling results with semi-analytic model predictions and discuss the implications.

We have developed an analytical method based on forward-modeling techniques to characterize the high-mass end of the red sequence (RS) galaxy population at redshift $z\sim0.55$, from the DR10 BOSS CMASS spectroscopic sample, which comprises $\sim600,000$ galaxies. The method, which follows an unbinned maximum likelihood approach, allows the deconvolution of the intrinsic CMASS colour-colour-magnitude distributions from photometric errors and selection effects. This procedure requires modeling the covariance matrix for the i-band magnitude, g-r colour and r-i colour using Stripe 82 multi-epoch data. Our results indicate that the error-deconvolved intrinsic RS distribution is consistent, within the photometric uncertainties, with a single point ($<0.05~{\rm{mag}}$) in the colour-colour plane at fixed magnitude, for a narrow redshift slice. We have computed the high-mass end ($^{0.55}M_i \lesssim -22$) of the $^{0.55}i$-band RS Luminosity Function (RS LF) in several redshift slices within the redshift range $0.52<z<0.63$. In this narrow redshift range, the evolution of the RS LF is consistent, within the uncertainties in the modeling, with a passively-evolving model with $\Phi_* = (7.248 \pm 0.204) \times10^{-4}$ Mpc$^{-3}$ mag$^{-1}$, fading at a rate of $1.5\pm0.4$ mag per unit redshift. We report RS completeness as a function of magnitude and redshift in the CMASS sample, which will facilitate a variety of galaxy-evolution and clustering studies using BOSS. Our forward-modeling method lays the foundations for future studies using other dark-energy surveys like eBOSS or DESI, which are affected by the same type of photometric blurring/selection effects.

We have collected near-infrared to X-ray data of 20 multi-epoch heavily reddened SDSS quasars to investigate the physical mechanism of reddening. Of these, J2317+0005 is found to be a UV cutoff quasar. Its continuum, which usually appears normal, decreases by a factor 3.5 at 3000{\AA}, compared to its more typical bright state during an interval of 23 days. During this sudden continuum cut-off, the broad emission line fluxes do not change, perhaps due to the large size of the Broad Line Region (BLR), r > 23 / (1+z) days. The UV continuum may have suffered a dramatic drop out. However, there are some difficulties with this explanation. Another possibility is that the intrinsic continuum did not change, but was temporarily blocked out, at least towards our line of sight. As indicated by X-ray observations, the continuum rapidly recovers after 42 days. A comparison of the bright state and dim states would imply an eclipse by a dusty cloud with a reddening curve having a remarkably sharp rise shortward of 3500{\AA}. Under the assumption of being eclipsed by a Keplerian dusty cloud, we characterized the cloud size with our observations, however, which is a little smaller than the 3000\AA\ continuum-emitting size inferred from accretion disk models. Therefore, we speculate this is due to a rapid outflow or inflow with a dusty cloud passing through our line-of-sight to the center.

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $\chi^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

We investigated the optical/ultraviolet (UV) color variations for a sample of 2169 quasars based on multi-epoch spectroscopy in the Sloan Digital Sky Survey (SDSS) data release seven (DR7) and data release nine (DR9). To correct the systematic difference between DR7 and DR9 due to the different instrumental setup, we produced a correction spectrum by using a sample of F-stars observed both in DR7 and DR9. The correction spectrum was then applied to quasars when comparing the spectra of DR7 with DR9. In each object, the color variation was explored by comparing the spectral index of the continuum power-law fit on the brightest spectrum with the faintest one, and also by the shape of their difference spectrum. In 1876 quasars with consistent color variations from two methods, we found that most sources (1755, $\sim 94\%$) show bluer-when-brighter (BWB) trend, and the redder-when-brighter (RWB) trend is only detected in 121 objects ($\sim 6\%$). The common BWB trend is supported by the bluer composite spectrum constructed from bright spectra than that from faint spectra, and also by the blue composite difference spectrum. The correction spectrum is proved to be highly reliable by comparing the composite spectrum from corrected DR9 and original DR7 spectra. Assuming that the optical/UV variability is triggered by fluctuations, RWB trend can likely be explained if the fluctuations occur firstly at outer disk region, and the inner disk region has not fully responded yet when the fluctuation being propagated inward. In contrast, the common BWB trend implies that the fluctuations are likely more often happening firstly in inner disk region.

Small- and intermediate-scale galaxy clustering can be used to establish the galaxy-halo connection to study galaxy formation and evolution and to tighten constraints on cosmological parameters. With the increasing precision of galaxy clustering measurements from ongoing and forthcoming large galaxy surveys, accurate models are required to interpret the data and extract relevant information. We introduce a method based on high-resolution N-body simulations to accurately and efficiently model the galaxy two-point correlation functions (2PCFs) in projected and redshift spaces. The basic idea is to tabulate all information of haloes in the simulations necessary for computing the galaxy 2PCFs within the framework of halo occupation distribution or conditional luminosity function. It is equivalent to populating galaxies to dark matter haloes and using the mock 2PCF measurements as the model predictions. Besides the accurate 2PCF calculations, the method is also fast and therefore enables an efficient exploration of the parameter space. As an example of the method, we decompose the redshift-space galaxy 2PCF into different components based on the type of galaxy pairs and show the redshift-space distortion effect in each component. The generalizations and limitations of the method are discussed.

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times.

We utilize $\Lambda$CDM halo occupation models of galaxy clustering to investigate the evolving stellar mass dependent clustering of galaxies in the PRIsm MUlti-object Survey (PRIMUS) and DEEP2 Redshift Survey over the past eight billion years of cosmic time, between $0.2<z<1.2$. These clustering measurements provide new constraints on the connections between dark matter halo properties and galaxy properties in the context of the evolving large-scale structure of the universe. Using both an analytic model and a set of mock galaxy catalogs, we find a strong correlation between central galaxy stellar mass and dark matter halo mass over the range $M_\mathrm{halo}\sim10^{11}$-$10^{13}~h^{-1}M_\odot$, approximately consistent with previous observations and theoretical predictions. However, the stellar-to-halo mass relation (SHMR) and the mass scale where star formation efficiency reaches a maximum appear to evolve more strongly than predicted by other models, including models based primarily on abundance-matching constraints. We find that the fraction of satellite galaxies in haloes of a given mass decreases significantly from $z\sim0.5$ to $z\sim0.9$, partly due to the fact that haloes at fixed mass are rarer at higher redshift and have lower abundances. We also find that the $M_1/M_\mathrm{min}$ ratio, a model parameter that quantifies the critical mass above which haloes host at least one satellite, decreases from $\approx20$ at $z\sim0$ to $\approx13$ at $z\sim0.9$. Considering the evolution of the subhalo mass function vis-\`{a}-vis satellite abundances, this trend has implications for relations between satellite galaxies and halo substructures and for intracluster mass, which we argue has grown due to stripped and disrupted satellites between $z\sim0.9$ and $z\sim0.5$.

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 present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

We present the measurements and modelling of the projected and redshift-space clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. For a volume-limited luminous red galaxy sample in the redshift range of $0.48<z<0.55$, we perform halo occupation distribution modelling of the small- and intermediate-scale ($0.1$--$60h^{-1}{\rm {Mpc}}$) projected and redshift-space two-point correlation functions, with an accurate model built on high resolution $N$-body simulations. To interpret the measured redshift-space distortions, the distribution of galaxy velocities must differ from that of the dark matter inside haloes of $\sim 10^{13}$--$10^{14}h^{-1}{\rm M_{\odot}}$, i.e. the data require the existence of galaxy velocity bias. Most notably, central galaxies on average are not at rest with respect to the core of their host haloes (defined by the inner 25% of particles around the halo potential minimum), but rather move around it with a 1D velocity dispersion of $0.22^{+0.03}_{-0.04}$ times that of the dark matter, implying a spatial offset from the centre at the level of $\lesssim$1% of the halo virial radius. The luminous satellite galaxies move more slowly than the dark matter, with velocities $0.86^{+0.08}_{-0.03}$ times those of the dark matter, which suggests that the velocity and spatial distributions of these satellites cannot both be unbiased. The constraints mainly arise from the Fingers-of-God effect at nonlinear scales and the smoothing to the Kaiser effect in the translinear regime; the robustness of the results is demonstrated by a variety of tests. In addition, no clear evidence is found for a strong luminosity dependence of the velocity bias. We discuss the implications of the existence of galaxy velocity bias for investigations of galaxy formation and cosmology.

In a sample of 60 quasars selected from Sloan Digital Sky Survey (SDSS) with at least six-epoch spectroscopy, we investigate the variability of emission lines and continuum luminosity at various aspects. A strong anti-correlation between the variability and continuum luminosity at 2500 {\AA} is found for the sample, which is consistent with previous works. In individual sources, we find that half of the sample objects have a trend of being bluer-when-brighter trend (BWB), while the remaining half exhibit redder-when-brighter trend (RWB). Although the mechanism for RWB is unclear, the effects of host galaxy contribution due to seeing variations can not be completely ruled out. As expected from photoionization model, the positive correlations between the broad emission line and continuum luminosity are found in most individual sources, as well as for the whole sample. We confirm the Baldwin effect in most individual objects and the whole sample, while a negative Baldwin effect is also found in several quasars, which can be at least partly (if not all) due to the host galaxy contamination. We find positive correlations between the broad emission line luminosity and line width in most individual quasars, as well as the whole sample, implying a more variable line base than the line core.

Rate constants for the N + C2 reaction have been measured in a continuous supersonic flow reactor over the range 57 K to 296 K by the relative rate technique employing the N + OH - H + NO reaction as a reference. Excess concentrations of atomic nitrogen were produced by the microwave discharge method and C2 and OH radicals were created by the in-situ pulsed laser photolysis of precursor molecules C2Br4 and H2O2 respectively. In parallel, quantum dynamics calculations were performed based on an accurate global potential energy surfaces for the three lowest lying quartet states of the C2N molecule. The 14A" potential energy surface is barrierless, having two deep potential wells corresponding to the NCC and CNC intermediates. Both the experimental and theoretical work show that the rate constant decreases to low temperature, although the experimentally measured values fall more rapidly than the theoretical ones except at the lowest temperatures. Astrochemical simulations indicate that this reaction could be the dominant source of CN in dense interstellar clouds.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

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.

The three-point correlation function (3PCF) for galaxies provides an opportunity to measure the non-Gaussianity generated from nonlinear structure formation and also probes information about galaxy formation and evolution that is generally not available from the two-point correlation function (2PCF). We measure the 3PCF of the Sloan Digital Sky Survey DR7 main sample galaxies in both redshift and projected spaces on scales up to 40Mpc/h. We explore the dependence of the 3PCF on galaxy stellar mass and color in order to constrain the formation and evolution for galaxies of different properties. The study of the dependence on these properties also helps better constrain the relation between galaxy stellar mass and color and the properties of their hosting dark-matter halos. We focus on the study of the reduced 3PCF, Q, defined as the ratio between the 3PCF and the sum of the products of the 2PCFs. We find a very weak stellar mass dependence of Q in both redshift and projected spaces. On small scales, more massive galaxies tend to have slightly higher amplitudes of Q. The shape dependence of Q is also weak on these small scales, regardless of stellar mass and color. The reduced 3PCF has a strong color dependence for the low-mass galaxies, while no significant dependence on color is found for the high-mass galaxies. Low-mass red galaxies have higher amplitudes and stronger shape dependence of the reduced 3PCF than the blue galaxies, implying that these low-mass red galaxies tend to populate filamentary structures. The linear galaxy bias model fails to interpret the color dependence of Q, emphasizing the importance of a nonvanishing quadratic bias parameter in the correct modeling of the galaxy color distribution.

We study the clustering of galaxies, as a function of their colour, from Data Release Ten (DR10) of the SDSS-III Baryon Oscillation Spectroscopic Survey. We select 122,967 galaxies with 0.43 < z < 0.7 into a "Blue" sample and 131,969 into a "Red" sample based on k+e corrected (to z=0.55) r-i colours and i band magnitudes. The samples are chosen to each contain more than 100,000 galaxies, have similar redshift distributions, and maximize the difference in clustering amplitude. The Red sample has a 40% larger bias than the Blue (b_Red/b_Blue = 1.39+-0.04), implying the Red galaxies occupy dark matter halos with an average mass that is 0.5 log Mo greater. Spherically averaged measurements of the correlation function, \xi 0, and the power spectrum are used to locate the position of the baryon acoustic oscillation (BAO) feature of both samples. Using \xi 0, we obtain distance scales, relative to our reference LCDM cosmology, of 1.010+-0.027 for the Red sample and 1.005+-0.031 for the Blue. After applying reconstruction, these measurements improve to 1.013+/-0.020 for the Red sample and 1.008+-0.026 for the Blue. For each sample, measurements of \xi 0 and the second multipole moment, \xi 2, of the anisotropic correlation function are used to determine the rate of structure growth, parameterized by f\sigma 8. We find f\sigma 8,Red = 0.511+-0.083, f\sigma 8,Blue = 0.509+/-0.085, and f\sigma 8,Cross = 0.423+-0.061 (from the cross-correlation between the Red and Blue samples). We use the covariance between the bias and growth measurements obtained from each sample and their cross-correlation to produce an optimally-combined measurement of f\sigma 8,comb = 0.443+-0.055. In no instance do we detect significant differences in distance scale or structure growth measurements obtained from the Blue and Red samples.

3C 66A was monitored by the BATC(Beijing-Arizona-Taipei-Connecticut) telescope from 2005 to 2008,1994 observations were obtained on 89 nights. Detailed research and analysis was performed on these observations in this paper. A long term burst occurred in the whole light curve. No intra-day variability was claimed in our campaign by intra-night light curve analysis. Time lag of shorter wavelenth preceding longer wavelength was shown by correlation analysis. The results showed that the optical spectral shape turned flatter when the source brightened, and the spectral variability indicator was bigger on shorter time-scale as determined by the color indices variation analysis.

We compiled a sample of 60 quasars with spectroscopy on at least six epochs from the Sloan Digital Sky Survey (SDSS) to study the variabilities of the spectral shape, the continuum and the emission lines luminosity. In this paper, we present the results of SDSS J030639.57+000343.1. We found a strong anti-correlation between the continuum luminosity at $5100 \AA$ and the spectral index, implying a bluer-when-brighter trend. The luminosity of the broad $\rm H_\alpha$ line is proportion to the continuum luminosity at $5100 \AA$. Correspondingly, we did not find strong correlation between the equivalent width of broad $\rm H\alpha$ and the continuum luminosity, i.e. no baldwin effect of broad $\rm H\alpha$ in this source.

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We present the measurement of the two-point cross-correlation function (CCF) of 8,198 Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) quasars and 349,608 DR10 CMASS galaxies from the Baryonic Oscillation Spectroscopic Survey (BOSS) at redshift <z>~0.5 (0.3<z<0.9). The cross-correlation function can be reasonably well fit by a power-law model xi_QG(r)=(r/r_0)^(-gamma) on projected scales of r_p=2-25 Mpc/h with r_0=6.61+-0.25 Mpc/h and gamma=1.69+-0.07. We estimate a quasar linear bias of b_Q=1.38+-0.10 at <z>=0.53 from the CCF measurements. This linear bias corresponds to a characteristic host halo mass of ~4x10^12 M_sun/h, compared to ~10^13 M_sun/h characteristic host halo mass for CMASS galaxies. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be db_Q/dlogL=0.20+-0.34 or 0.11+-0.32 (depending on different luminosity divisions) for quasar luminosities -23.5>M_i(z=2)>-25.5, implying a weak luminosity dependence of quasar clustering for the bright end of the quasar population at <z>~0.5. We compare our measurements with theoretical predictions, Halo Occupation Distribution (HOD) models and mock catalogs. These comparisons suggest quasars reside in a broad range of host halos, and the host halo mass distributions significantly overlap with each other for quasars at different luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF equally well, but with different outcomes such as the satellite fraction and host halo mass distribution. These results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD approach and the need for additional information in populating quasars in dark matter halos with HOD. [Abridged]

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We analyze the optical spectra of massive (log M*/Msun > 11.4) radio-loud galaxies at z~0.2 and z~0.6. By comparing stellar population parameters of these radio-loud samples with radio-quiet control samples, we investigate how the presence of a radio-emitting jet relates to the recent star formation history of the host galaxy. We also investigate how the emission-line properties of the radio galaxies evolve with redshift by stacking their spectra. Our main results are the following. (1) Both at low and at high redshift, half as many radio-loud as radio-quiet galaxies have experienced significant star formation in the past Gyr. (2) The Balmer absorption line properties of massive galaxies that have experienced recent star formation show that star formation occurred as a burst in many of these systems. (3) Both the radio and the emission-line luminosity of radio AGN evolve significantly with redshift. However, radio galaxies with similar stellar population parameters, have similar emission-line properties both at high- and at low-redshift. These results suggest that massive galaxies experience cyclical episodes of gas accretion, star formation and black hole growth, followed by the production of a radio jet that shuts down further activity. The behaviour of galaxies with log M*/Msun > 11.4 is the same at z = 0.6 as it is at z = 0.2, except that higher redshift galaxies experience more star formation and black hole growth and produce more luminous radio jets during each accretion cycle.

In fiber-fed galaxy redshift surveys, the finite size of the fiber plugs prevents two fibers from being placed too close to one another, limiting the ability of studying galaxy clustering on all scales. We present a new method for correcting such fiber collision effects in galaxy clustering statistics based on spectroscopic observations. Our method makes use of observations in tile overlap regions to measure the contributions from the collided population, and to therefore recover the full clustering statistics. The method is rooted in solid theoretical ground and is tested extensively on mock galaxy catalogs. We demonstrate that our method can well recover the projected and the full three-dimensional redshift-space two-point correlation functions on scales both below and above the fiber collision scale, superior to the commonly used nearest neighbor and angular correction methods. We discuss potential systematic effects in our method. The statistical correction accuracy of our method is only limited by sample variance, which scales down with (the square root of) the volume probed. For a sample similar to the final SDSS-III BOSS galaxy sample, the statistical correction error is expected to be at the level of 1% on scales 0.1--30Mpc/h for the two-point correlation functions. The systematic error only occurs on small scales, caused by non-perfect correction of collision multiplets, and its magnitude is expected to be smaller than 5%. Our correction method, which can be generalized to other clustering statistics as well, enables more accurate measurements of full three-dimensional galaxy clustering on all scales with galaxy redshift surveys. (abridged)

We analyze the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of 3,275 deg^2. We identify, characterize, and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the `CMASS' sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120h^-1Mpc or k < 0.01hMpc^-1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. We use mock galaxy catalogs that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on correlation function measurements and that this systematic error is negligible for the spherically averaged correlation function. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature (Anderson et al. 2012), constrain cosmological models using the full shape of the correlation function (Sanchez et al. 2012), and measure the rate of structure growth (Reid et al. 2012). (abridged)

We present a comprehensive study of 250,000 galaxies targeted by the Baryon Oscillation Spectroscopic Survey (BOSS) up to z ~ 0.7 with the specific goal of identifying and characterising a population of galaxies that has evolved without significant merging. We compute a likelihood that each BOSS galaxy is a progenitor of the Luminous Red Galaxies (LRGs) sample, targeted by SDSS-I/II up z ~ 0.5, by using the fossil record of LRGs and their inferred star-formation histories, metallicity histories and dust content. We determine merger rates, luminosity growth rates and the evolution of the large-scale clustering between the two surveys, and we investigate the effect of using different stellar population synthesis models in our conclusions. We demonstrate that our sample is slowly evolving (of the order of 2 +/- 1.5% per Gyr by merging). Our conclusions refer to the bright and massive end of the galaxy population, with Mi0.55 < -22, and M* > 1E11.2 Msolar, corresponding roughly to 95% and 40% of the LRGs and BOSS galaxy populations, respectively. Our analysis further shows that any possible excess of flux in BOSS galaxies, when compared to LRGs, from potentially unresolved targets at z ~ 0.55 must be less than 1% in the r0.55-band (approximately equivalent to the g-band in the rest-frame of galaxies at z=0.55). When weighting the BOSS galaxies based on the predicted properties of the LRGs, and restricting the analysis to the reddest BOSS galaxies, we find an evolution of the large-scale clustering that is consistent with dynamical passive evolution, assuming a standard cosmology. We conclude that our likelihoods give a weighted sample that is as clean and as close to passive evolution (in dynamical terms, i.e. no or negligible merging) as possible, and that is optimal for cosmological studies.

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

We investigate the effects of potential sources of systematic error on the angular and photometric redshift, z_phot, distributions of a sample of redshift 0.4 < z < 0.7 massive galaxies whose selection matches that of the Baryon Oscillation Spectroscopic Survey (BOSS) constant mass sample. Utilizing over 112,778 BOSS spectra as a training sample, we produce a photometric redshift catalog for the galaxies in the SDSS DR8 imaging area that, after masking, covers nearly one quarter of the sky (9,913 square degrees). We investigate fluctuations in the number density of objects in this sample as a function of Galactic extinction, seeing, stellar density, sky background, airmass, photometric offset, and North/South Galactic hemisphere. We find that the presence of stars of comparable magnitudes to our galaxies (which are not traditionally masked) effectively remove area. Failing to correct for such stars can produce systematic errors on the measured angular auto-correlation function, w, that are larger than its statistical uncertainty. We describe how one can effectively mask for the presence of the stars, without removing any galaxies from the sample, and minimize the systematic error. Additionally, we apply two separate methods that can be used to correct the systematic errors imparted by any parameter that can be turned into a map on the sky. We find that failing to properly account for varying sky background introduces a systematic error on w. We measure w, in four z_phot slices of width 0.05 between 0.45 < z_phot < 0.65 and find that the measurements, after correcting for the systematic effects of stars and sky background, are generally consistent with a generic LambdaCDM model, at scales up to 60 degrees. At scales greater than 3 degrees and z_phot > 0.5, the magnitude of the corrections we apply are greater than the statistical uncertainty in w.

The variations of gravity were measured with a high precision LaCoste-Romberg D gravimeter during a total solar eclipse to investigate the effect of solar eclipse on the gravitational field. The observed anomaly $(7.0 \pm 2.7) \times 10^{-8}$ m/s$^2$ during the eclipse implies that there may be a shielding property of gravitation.

By making use of Duan-Ge's decomposition theory of gauge potential and the Duan's topological current theory proposed by Prof. Duan Yi-Shi, we study a two component superfluid Bose condensed system, which is supposed being realized in the interior of neutron stars in the form of a coexistent neutron superfluid and protonic superconductor. We propose that this system possesses vortex lines. The topological charge of the vortex lines are characterized by the Hopf indices and the Brower degrees of $\phi$-mapping.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

We study whether the bias factors of galaxies can be unbiasedly recovered from their power spectra and bispectra. We use a set of numerical N-body simulations and construct large mock galaxy catalogs based upon the semi-analytical model of Croton et al. (2006). We measure the reduced bispectra for galaxies of different luminosity, and determine the linear and first nonlinear bias factors from their bispectra. We find that on large scales down to that of the wavenumber k=0.1h/Mpc, the bias factors b1 and b2 are nearly constant, and b1 obtained with the bispectrum method agrees very well with the expected value. The nonlinear bias factor b2 is negative, except for the most luminous galaxies with M<-23 which have a positive b2. The behavior of b2 of galaxies is consistent with the b2 mass dependence of their host halos. We show that it is essential to have an accurate estimation of the dark matter bispectrum in order to have an unbiased measurement of b1 and b2. We also test the analytical approach of incorporating halo occupation distribution to model the galaxy power spectrum and bispectrum. The halo model predictions do not fit the simulation results well on the precision requirement of current cosmological studies.

We use a set of numerical N-body simulations to study the large-scale behavior of the reduced bispectrum of dark matter and compare the results with the second-order perturbation theory and the halo models for different halo mass functions. We find that the second-order perturbation theory (PT2) agrees with the simulations fairly well on large scales of k<0.05 h/Mpc, but it shows a signature of deviation as the scale goes down. Even on the largest scale where the bispectrum can be measured reasonably well in our simulations, the inconsistency between PT2 and the simulations appears for the colinear triangle shapes. For the halo model, we find that it can only serve as a qualitative method to help study the behavior of Q on large scales and also on relatively small scales. The failure of second-order perturbation theory will also affect the precise determination of the halo models, since they are connected through the 3-halo term in the halo model. The 2-halo term has too much contribution on the large scales, which is the main reason for the halo model to overpredict the bispectrum on the large scales. Since neither of the models can provide a satisfying description for the bispectrum on scales of about 0.1h/Mpc for the requirement of precision cosmology, we release the reduced bispectrum of dark matter on a large range of scales for future analytical modeling of the bispectrum.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

We study the layer of electrons on bare strange star surfaces, taking the Dirac exchange-energy into account. Because electrons are fermions, the electron wave function must be of exchange-antisymmetry. The Dirac exchange-energy originates, consequently, from the exchange-antisymmetry of electron wave functions. This consideration may result in changing the electron distribution and the electric field on the surface of bare strange star. The strong magnetic field effect on the structures of the electrospheres is also discussed.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

Many quasi-simultaneous optical observations of 17 blazars are obtained from previous papers published over the last 19 years in order to investigate the spectral slope variability and understand the radiation mechanism of blazars. The long-period dereddened optical spectral slopes are calculated. We analyse the average spectral slope distribution, which suggests that the spectra of flat spectrum radio quasars (FSRQs) and high energy peaked BL Lac objects (HBLs) are probably deformed by other emission components. The average spectral slopes of low energy peaked BL Lac objects(LBLs), which scatter around 1.5, show a good accordance with the synchrotron self-Compton (SSC) loss-dominated model. We present and discuss the variability between the spectral slope and optical luminosity. The spectra of all HBLs and LBLs get flatter when they turn brighter, while for FSRQs this trend does not exist or may even be reversed. This phenomenon may imply that there is a thermal contribution to the optical spectrum for FSRQs. For the FSRQ 1156+295, there is a hint that the slope gets flatter at both the brightest and faintest states. Our result shows that three subclasses locate in different regions in the pattern of slope variability indicator versus average spectral slope. The relativistic jet mechanism is supported by the significant correlation between the optical Doppler factor and the average spectral slope.

A new single-dynamical-scalar-field model of dark energy is proposed, in which either higher derivative terms nor structures of extra dimension are needed. With the help of a fixed background vector field, the parameter for the effective equation of state of dark energy may cross $w=-1$ in the evolution of the universe. After suitable choice of the potential, the crossing $w=-1$ and transition from decelerating to accelerating occur at $z\approx 0.2$ and $z\approx 1.7$, respectively.

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.