We have performed microwave diagnostics of the magnetic field strengths in solar flare loops based on the theory of gyrosynchrotron emission. From Nobeyama Radioheliograph observations of three flare events at 17 and 34 GHz, we obtained the degree of circular polarization and the spectral index of microwave flux density, which were then used to map the magnetic field strengths in post-flare loops. Our results show that the magnetic field strength typically decreases from ~800 G near the loop footpoints to ~100 G at a height of 10--25 Mm. Comparison of our results with magnetic field modeling using a flux rope insertion method is also discussed. Our study demonstrates the potential of microwave imaging observations, even at only two frequencies, in diagnosing the coronal magnetic field of flaring regions.

Filaments are very common physical phenomena on the Sun and are often taken as important proxies of solar magnetic activities. The study of filaments has become a hot topic in the space weather research. For a more comprehensive understanding of filaments, especially for an understanding of solar activities of multiple solar cycles, it is necessary to perform a combined multifeature analysis by constructing a data set of multiple solar cycle data. To achieve this goal, we constructed a centennial data set that covers the H$\alpha$ data from five observatories around the world. During the data set construction, we encountered varieties of problems, such as data fusion, accurate determination of the solar edge, classifying data by quality, dynamic threshold, and so on, which arose mainly due to multiple sources and a large time span of data. But fortunately, these problems were well solved. The data set includes seven types of data products and eight types of feature parameters with which we can implement the functions of data searching and statistical analyses. It has the characteristics of better continuity and highly complementary to space observation data, especially in the wavelengths not covered by space observations, and covers many solar cycles (including more than 60 yr of high-cadence data). We expect that this new comprehensive data set as well as the tools will help researchers to significantly speed up their search for features or events of interest, for either statistical or case study purposes, and possibly help them get a better and more comprehensive understanding of solar filament mechanisms.

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

We explore the cosmological multi-tracer synergies between a DESI-like Emission Line Galaxy (ELG) distribution and a Tianlai-like 21cm intensity map. We use simulated maps generated from a particle simulation in the light-cone volume (Horizon Run 4), sky-trimmed and including the effects of foreground contamination, its removal, and instrument noise. We first validate how the foreground residual affects the recovered 21cm signal by putting different levels of foreground contamination into the 21cm maps. We find that the contamination cannot be ignored in the angular auto-power spectra of HI even when it is small, but has no influence on the accuracy of the angular cross-power spectra between HI and galaxies. In the foreground-cleaned map case, as information is lost in the cleaning procedure, there is also a bias in the cross-power spectrum. However, we found that the bias from the cross-power spectrum is scale-independent, which is easily parameterised as part of the model, while the offset in the HI auto-power spectrum is non-linear. In particular, we tested that the cross power also benefits from the cancellation of the bias in the power spectrum measurement that is induced by the instrument noise, which changes the shape of the auto-power spectra but leaves the cross-power unaffected. We then modelled the angular cross-power spectra to fit the BAO feature in broadband shape of the angular cross power spectrum, including contamination from the residual foreground and the effect of instrument noise. We forecast a constraint on the angular diameter distance $D_\mathrm{A}$ for the Tianlai Pathfinder redshift $0.775<z<1.03$, giving a distance measurement with a precision of 3.3% at that redshift.

A photometric redshift sample of Luminous Red Galaxies (hereafter LRGs) obtained from The DECam Legacy Survey (DECaLS) is analysed to probe cosmic distances by exploiting the wedge approach of the two-point correlation function. Although the cosmological information is highly contaminated by the uncertainties existing in the photometric redshifts from the galaxy map, an angular diameter distance can be probed at the perpendicular configuration in which the measured correlation function is minimally contaminated. An ensemble of wedged correlation functions selected up to a given threshold based on having the least contamination was studied in the previous work (arXiv:1903.09651v2 [astro-ph.CO]) using simulations, and the extracted cosmological information was unbiased within this threshold. We apply the same methodology for analysing the LRG sample from DECaLS which will provide the optical imaging for targeting two-thirds of the DESI footprint and measure the angular diameter distances at $z=0.69$ and $z=0.87$ to be $D_{A}(0.697)=(1499 \pm 77\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ and $D_{A}(0.874)=(1680 \pm 109\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ with a fractional error of 5.14% and 6.48% respectively. We obtain a value of $H_{0}=67.59\pm5.52$ km/s/Mpc which supports the $H_0$ measured by all other BAO results and is consistent with $\Lambda$CDM model.

White-light flares (WLFs), first observed in 1859, refer to a type of solar flares showing an obvious enhancement of the visible continuum emission. This type of enhancement often occurs in most energetic flares, and is usually interpreted as a consequence of efficient heating in the lower solar atmosphere through non-thermal electrons propagating downward from the energy release site in the corona. However, this coronal-reconnection model has difficulty in explaining the recently discovered small WLFs. Here we report a C2.3 white-light flare, which are associated with several observational phenomena: fast decrease in opposite-polarity photospheric magnetic fluxes, disappearance of two adjacent pores, significant heating of the lower chromosphere, negligible increase of hard X-ray flux, and an associated U-shaped magnetic field configuration. All these suggest that this white-light flare is powered by magnetic reconnection in the lower part of the solar atmosphere rather than by reconnection higher up in the corona.

The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate HI intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the HI sky for the frequency range 700-800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21cm distribution on large-scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information.

As an introduction of a kinematic survey of Magellanic Cloud (MC) star clusters, we report on the dynamical masses and mass-to-light ($M/L$) ratios of NGC 419 (SMC) and NGC 1846 (LMC). We have obtained more than one hundred high-resolution stellar spectra in and around each cluster using the multi-object spectrograph M2FS on the $Magellan$/Clay Telescope. Line-of-sight velocities and positions of the stars observed in each cluster were used as input to an expectation-maximization algorithm used to estimate cluster membership probabilities, resulting in samples of 46 and 52 likely members ($P_{M}\geq 50$%) in NGC 419 and NGC 1846, respectively. This process employed single-mass King models constrained by the structural parameters of the clusters and provided self-consistent dynamical mass estimates for both clusters. Our best-fit results show that NGC 419 has a projected central velocity dispersion of $2.44^{+0.37}_{-0.21}\ {\rm km\,s^{-1}}$, corresponding to a total mass of $7.6^{+2.5}_{-1.3}\times10^4\ {\rm M}_{\odot}$ and $V$-band $M/L$ ratio of $0.22^{+0.08}_{-0.05}$ in solar units. For NGC 1846, the corresponding results are $2.04^{+0.28}_{-0.24}\ {\rm km\,s^{-1}}$, $5.4^{+1.5}_{-1.4}\times10^4\ {\rm M}_{\odot}$ and $0.32^{+0.11}_{-0.11}$. The mean metallicities of NGC 419 and NGC 1846 are found to be $\rm [Fe/H]=-0.84\pm0.19$ and $-0.70\pm0.08$, respectively, based on the spectra of likely cluster members. We find marginal statistical evidence of rotation in both clusters, though in neither cluster does rotation alter our mass estimates significantly. We critically compare our findings with those of previous kinematic studies of these two clusters in order to evaluate the consistency of our observational results and analytic tools.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

The galaxy distributions along the line-of-sight are significantly contaminated by the uncertainty on redshift measurements obtained through multiband photometry, which makes it difficult to get cosmic distance information measured from baryon acoustic oscillations, or growth functions probed by redshift distortions. We investigate the propagation of the uncertainties into large scale clustering by exploiting all known estimators, and propose the wedge approach as a promising analysis tool to extract cosmic distance information still remaining in the photometric galaxy samples. We test our method using simulated galaxy maps with photometric uncertainties of $\sigma_{0} =\left(0.01, 0.02, 0.03\right)$. The measured anisotropy correlation function $\xi$ is binned into the radial direction of $s$ and the angular direction of $\mu$, and the variations of $\xi(s,\mu)$ with perpendicular and radial cosmic distance measures of $D_A$ and $H^{-1}$ are theoretically estimated by an improved RSD model. Although the radial cosmic distance $H^{-1}$ is unable to be probed from any of the three photometric galaxy samples, the perpendicular component of $D_A$ is verified to be accurately measured even after the full marginalisation of $H^{-1}$. We measure $D_A$ with approximately 6% precision which is nearly equivalent to what we can expect from spectroscopic DR12 CMASS galaxy samples.

We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in $^{12}$CO, and 29 (52%) are also detected in $^{13}$CO above a 3$\sigma$ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661km s$^{-1}$, and $\sim$54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H$_2$ conversion factor of $\langle\alpha\rangle=1.8^{+1.3}_{-0.8}M_\odot$(K km s$^{-1}$pc$^{2}$)$^{-1}$; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of $400-600$Myr for the (U)LIRGs, compared to the 1.3Gyr for local spiral galaxies. Finally, we re-examine the relationship between the $^{12}$CO/$^{13}$CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems.

We present the first results of a high-resolution Karl G. Jansky Very Large Array (VLA) imaging survey of luminous and ultra-luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). From the full sample of 68 galaxies, we have selected 25 LIRGs that show resolved extended emission at sufficient sensitivity to image individual regions of star-formation activity beyond the nucleus.~With wideband radio continuum observations, which sample the frequency range from $3-33$ GHz, we have made extinction-free measurements of the luminosities and spectral indicies for a total of 48 individual star-forming regions identified as having de-projected galactocentric radii ($r_{G}$) that lie outside the 13.2$\mu$m core of the galaxy.~The median $3-33$ GHz spectral index and 33 GHz thermal fraction measured for these "extranuclear" regions is $-0.51 \pm 0.13$ and $65 \pm 11\%$ respectively.~These values are consistent with measurements made on matched spatial scales in normal star-forming galaxies, and suggests that these regions are more heavily-dominated by thermal free-free emission relative to the centers of local ULIRGs.~Further, we find that the median star-formation rate derived for these regions is $\sim 1 M_{\odot}$ yr$^{-1}$, and when we place them on the sub-galactic star-forming main sequence of galaxies (SFMS), we find they are offset from their host galaxies' globally-averaged specific star-formation rates (sSFRs).~We conclude that while nuclear starburst activity drives LIRGs above the SFMS, extranuclear star-formation still proceeds in a more extreme fashion relative to what is seen in local spiral galaxies.

The mapping of galaxy clustering from real space to redshift space introduces the anisotropic property to the measured galaxy density power spectrum in redshift space, known as the redshift space distortion (RSD) effect. The mapping formula is intrinsically non-linear, which is complicated by the higher order polynomials due to indefinite orders of cross correlations between density and velocity fields, and the Finger--of--God (FoG) effect due to the randomness of the galaxy peculiar velocity field. In previous works, we have verified the robustness of advanced TNS mapping formula in our hybrid RSD model in dark matter case, where the halo bias models are not taken into account for the halo mapping formula in redshift space. Using 100 realizations of halo catalogs in N-body simulations, we find that our halo RSD model with the known halo bias model and the effective FoG function accurately predicts the halo power spectrum measurements, within 1$\sim$2% accuracy up to $k\sim 0.2h$/Mpc, depending on different halo masses and redshifts.

A moderately strong vector repulsion between quarks in dense quark matter is needed to explain how a quark core can support neutron stars heavier than two solar masses. We study this repulsion, parametrized by a four-fermion interaction with coupling g_V, in terms of non-perturbative gluon exchange in QCD in the Landau gauge. Matching the energy of quark matter, g_V n_q^2 (where n_q is the number density of quarks) with the quark exchange energy calculated in QCD with a gluon propagator parametrized by a finite gluon mass m_g and a frozen coupling alpha_s, at moderate quark densities, we find that gluon masses m_g in the range 200 - 600 MeV and alpha_s = 2 - 4 lead to a g_V consistent with neutron star phenomenology. Estimating the effects of quark masses and a color-flavor-locked (CFL) pairing gap, we find that g_V can be well approximated by a flavor-symmetric, decreasing function of density. We briefly discuss similar matchings for the isovector repulsion and for the pairing attraction.

The ``force-free" magnetosphere of a Kerr black-hole in the axisymmetric steady-state possesses two integral functions of the stream-function $\Psi$; the field-line angular-velocity (FLAV) $\Omega_{\rm F} (\Psi)$, and the current-function $I(\Psi)$. For the zero-angular-momentum-observer's FLAV, i.e., $\Omega_{F {\omega}}\equiv\Omega_{\rm F}-\omega$, the iso-rotation law breaks down; $\Omega_{\rm F {\omega}}$ vanishes at the null-surface $\cal{S}_{\rm N}$ with $\omega=\Omega_{\rm F}$, where $\omega$ is the frame-dragging angular-velocity, and this gives rise to a subsequent breakdown of the force-free and freezing-in conditions at $\cal{S}_{\rm N}$, and the particle-velocity and the electric-current must vanish at $\cal{S}_{\rm N}$. These impose strong constraints on a possible gap to be constructed at $\cal{S}_{\rm N}$ between the force-free domains, $\cal{D}_{\rm out}$ and $\cal{D}_{\rm in}$. Current circuits, $\cal{C}_{\rm out}$ and $\cal{C}_{\rm in}$, must be closed in each domain, with EMFs $\cal{E}_{\rm out}$ and $\cal{E}_{\rm in}$, respectively, oppositely directed at the inductive membrane $\cal{S}_{\rm N}$ and with resistances $I_{\rm (out)}$ and $I_{\rm (in)}$ at the resistive membranes $\cal{S}_{{\rm ff}\infty}$ and $\cal{S}_{{\rm ffH}}$. The particles pair-produced in $\cal{G}_{\rm N}$ must be ``zero-angular-momentum-particles", which will be dense enough to pin down the poloidal field-lines ${\bf B}_p$ threading $\cal{G}_{\rm N}$, to ensure magnetisation of $\cal{G}_{\rm N}$. The zero-angular-momentum-state of $\cal{G}_{\rm N}$ makes it possible in $\cal{D}_{\rm in}$ that the outgoing-flux $I_{\rm (in)}$ of {\em positive}-angular-momentum from the hole is equivalent to the ingoing-flux $I^{\rm (in)}$ of {\em negative}-angular-momentum from $\cal{G}_{\rm N}$, i.e., $I_{\rm (in)}(\Psi)=-I^{\rm (in)}(\Psi)$. A twin-pulsar model is proposed.

In the preliminary trajectory design of the multi-target rendezvous problem, a model that can quickly estimate the cost of the orbital transfer is essential. The estimation of the transfer time using solar sail between two arbitrary orbits is difficult and usually requires to solve an optimal control problem. Inspired by the successful applications of the deep neural network in nonlinear regression, this work explores the possibility and effectiveness of mapping the transfer time for solar sail from the orbital characteristics using the deep neural network. Furthermore, the Monte Carlo Tree Search method is investigated and used to search the optimal sequence considering a multi-asteroid exploration problem. The obtained sequences from preliminary design will be solved and verified by sequentially solving the optimal control problem. Two examples of different application backgrounds validate the effectiveness of the proposed approach.

We present an empirical stellar spectra library created using spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5. This library represents a uniform data set ranging from 3750 through 8500K in effective temperature (T_eff), from -2.5 through +1.0 dex in metallicity ([Fe/H]), and from 0 to 5.0 dex in gravity (log g). The spectra in the library have resolution R~1800, with well-calibrated fluxes and rest-framed wavelengths. Using a large number of red stars observed by LAMOST, we generated denser K type templates to fill in data missing from current empirical spectral libraries, particularly the late K type. For K giants, we calibrated the spectroscopic surface gravities against the asteroseismic surface gravities. To verify the reliability of the parameters labeled for this library, we performed an internal cross-validation by using a chi^2 minimization method to interpolate parameters of each individual spectrum using the remaining spectra in the library. We obtained precisions of 41 K, 0.11 dex, and 0.05 dex for T_eff, log g, and [Fe/H], respectively, which means the templates are labeled with correct stellar parameters. Through external comparisons, we confirmed that measurements of the stellar parameters through this library can achieve accuracies of approximately 125K in T_eff, 0.1 dex in [Fe/H] and 0.20 dex in log g without systematic offset. This empirical library is useful for stellar parameter measurements because it has large parameter coverage and full wavelength coverage from 3800 to 8900A.

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

Using observations by the Solar Dynamics Observatory from June 2010 to December 2017, we have performed the first statistical investigation of circular-ribbon flares (CFs) and examined the white-light emission in these CFs. We find 90 CFs occurring in 36 active regions (ARs), including 8 X-class, 34 M-class, 48 C- and B-class flares. The occurrence rate of white-light flares (WLFs) is 100\% (8/8) for X-class CFs, $\sim$62\% (21/34) for M-class CFs, and $\sim$8\% (4/48) for C- and B-class CFs. Sometimes we observe several CFs in a single AR, and nearly all of them are WLFs. Compared to normal CFs, CFs with white-light enhancement tend to have a shorter duration, smaller size, stronger electric current and more complicated magnetic field. We find that for X-class WLFs, the white-light enhancement is positively correlated with the flare class, implying that the white-light enhancement is largely determined by the amount of released energy. However, there is no such correlation for M- and C-class WLFs, suggesting that other factors such as the time scale, spatial scale and magnetic field complexity may play important roles in the generation of white-light emission if the released energy is not high enough.

We investigate the electron-positron pair cascade taking place in the magnetosphere of a rapidly rotating black hole. Because of the spacetime frame dragging, the Goldreich-Julian charge density changes sign in the vicinity of the event horizon, which leads to an occurrence of a magnetic-field aligned electric field, in the same way as the pulsar outer-magnetospheric accelerator. In this lepton accelerator, electrons and positrons are accelerated in the opposite directions, to emit copious gamma-rays via the curvature and inverse-Compton processes. We examine a stationary pair cascade, and show that a stellar-mass black hole moving in a gaseous cloud can emit a detectable very-high-energy flux, provided that the black hole is extremely rotating and that the distance is less than about 1 kpc. We argue that the gamma-ray image will have a point-like morphology, and demonstrate that their gamma-ray spectra have a broad peak around 0.01-1 GeV and a sharp peak around 0.1 TeV, that the accelerators become most luminous when the mass accretion rate becomes about 0.01% of the Eddington rate, and that the predicted gamma-ray flux little changes in a wide range of magnetospheric currents. An implication of the stability of such a stationary gap is discussed.

We report on the status of an ongoing effort to develop arrays of horn-coupled, polarization-sensitive microwave kinetic inductance detectors (MKIDs) that are each sensitive to two spectral bands between 125 and 280 GHz. These multi-chroic MKID arrays are tailored for next-generation, large-detector-count experiments that are being designed to simultaneously characterize the polarization properties of both the cosmic microwave background (CMB) and Galactic dust emission. We present our device design and describe laboratory-based measurement results from two 23-element prototype arrays. From dark measurements of our first engineering array we demonstrated a multiplexing factor of 92, showed the resonators respond to bath temperature changes as expected, and found that the fabrication yield was 100%. From our first optically loaded array we found the MKIDs respond to millimeter-wave pulses, additional optical characterization measurements are ongoing. We end by discussing our plans for scaling up this technology to kilo-pixel arrays over the next two years.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

Solid evidence of magnetic reconnection is rarely reported within sunspots, the darkest regions with the strongest magnetic fields and lowest temperatures in the solar atmosphere. Using the world's largest solar telescope, the 1.6-meter Goode Solar Telescope, we detect prevalent reconnection through frequently occurring fine-scale jets in the H${\alpha}$ line wings at light bridges, the bright lanes that may divide the dark sunspot core into multiple parts. Many jets have an inverted Y-shape, shown by models to be typical of reconnection in a unipolar field environment. Simultaneous spectral imaging data from the Interface Region Imaging Spectrograph show that the reconnection drives bidirectional flows up to 200~km~s$^{-1}$, and that the weakly ionized plasma is heated by at least an order of magnitude up to $\sim$80,000 K. Such highly dynamic reconnection jets and efficient heating should be properly accounted for in future modeling efforts of sunspots. Our observations also reveal that the surge-like activity previously reported above light bridges in some chromospheric passbands such as the H${\alpha}$ core has two components: the ever-present short surges likely to be related to the upward leakage of magnetoacoustic waves from the photosphere, and the occasionally occurring long and fast surges that are obviously caused by the intermittent reconnection jets.

The observed power spectrum in redshift space appears distorted due to the peculiar motion of galaxies, known as redshift-space distortions (RSD). While all the effects in RSD are accounted for by the simple mapping formula from real to redshift spaces, accurately modeling redshift-space power spectrum is rather difficult due to the non-perturbative properties of the mapping. Still, however, a perturbative treatment may be applied to the power spectrum at large-scales, and on top of a careful modeling of the Finger-of-God effect caused by the small-scale random motion, the redshift-space power spectrum can be expressed as a series of expansion which contains the higher-order correlations of density and velocity fields. In our previous work [JCAP 8 (Aug., 2016) 050], we provide a perturbation-theory inspired model for power spectrum in which the higher-order correlations are evaluated directly from the cosmological $N$-body simulations. Adopting a simple Gaussian ansatz for Finger-of-God effect, the model is shown to quantitatively describe the simulation results. Here, we further push this approach, and present an accurate power spectrum template which can be used to estimate the growth of structure as a key to probe gravity on cosmological scales. Based on the simulations, we first calibrate the uncertainties and systematics in the pertrubation theory calculation in a fiducial cosmological model. Then, using the scaling relations, the calibrated power spectrum template is applied to a different cosmological model. We demonstrate that with our new template, the best-fitted growth functions are shown to reproduce the fiducial values in a good accuracy of 1 \% at $k<0.18 \hompc$ for cosmologies with different Hubble parameters.

We present observations of an M5.7 white-light flare (WLF) associated with a small filament eruption in NOAA active region 11476 on 2012 May 10. During this flare, a circular flare ribbon appeared in the east and a remote brightening occurred in the northwest of the active region. Multi-wavelength data are employed to analyze the WLF, including white light (WL), ultraviolet, extreme ultraviolet, hard X-ray (HXR) and microwave. A close spatial and temporal relationship between the WL, HXR and microwave emissions is found in this WLF. However, the peak time of the WL emission lagged that of the HXR and microwave emissions by about 1-2 minutes. Such a result tends to support the back-warming mechanism for the WL emission. Interestingly, the enhanced WL emission occurred at the two footpoints of the filament. Through forced and potential field extrapolations, we find that the three-dimensional magnetic field in the flare region has a fan-spine feature and that a flux rope lies under the dome-like field structure. We describe the entire process of flare evolution into several steps, each producing the sequent brightening below the filament, the circular flare ribbons and the WL enhancement respectively. We suggest that reconnection between the magnetic field of the filament and the overlying magnetic field or reconnection within the flux rope leads to the WL enhancement.

There are two goals in this study. One is to investigate how frequently white-light flares (WLFs) occur in a flare-productive active region (NOAA active region 11515). The other is to investigate the relationship between WLFs and magnetic transients (MTs). We use the full-disk continuum filtergrams and line-of-sight magnetograms taken by SDO/HMI to identify WLFs and MTs, respectively. Images taken by SDO/AIA are also used to show the morphology of the flares in the upper atmosphere. We found at least 20 WLFs out of a total of 70 flares above C class (28.6%) in NOAA active region 11515 during its passage across the solar disk. Each of these WLFs occurred in a small region, with a short duration of about 5 minutes. The enhancement of white-light continuum intensity is usually small, with an average enhancement of 8.1%. The 20 WLFs observed were found along an unusual configuration of the magnetic field characterized by a narrow ribbon of negative field. Furthermore, the WLFs were found to be accompanied by MTs, with radical changes in magnetic field strength (or even a sign reversal) observed during the flare. In contrast, there is no obvious signature of MTs in those 50 flares without white-light enhancements. Our results suggest that WLFs occur much more frequently than what was previously thought, with most WLFs being fairly weak enhancements. This may explain why WLFs are not frequently reported. Our observations also suggest that MTs and WLFs are closely related and appear co-spatial and co-temporal, when considering HMI data. A larger enhancement of WL emission is often accompanied by a larger change of the line-of-sight component of the unsigned magnetic field. Considering the close relationship between MTs and WLFs, many previously reported flares with MTs may be WLFs.

We review the equation of state of matter in neutron stars from the solid crust through the liquid nuclear matter interior to the quark regime at higher densities. We focus in detail on the question of how quark matter appears in neutron stars, and how it affects the equation of state. After discussing the crust and liquid nuclear matter in the core we briefly review aspects of microscopic quark physics relevant to neutron stars, and quark models of dense matter based on the Nambu--Jona-Lasinio framework, in which gluonic processes are replaced by effective quark interactions. We turn then to describing equations of state useful for interpretation of both electromagnetic and gravitational observations, reviewing the emerging picture of hadron-quark continuity in which hadronic matter turns relatively smoothly, with at most only a weak first order transition, into quark matter with increasing density. We review construction of unified equations of state that interpolate between the reasonably well understood nuclear matter regime at low densities and the quark matter regime at higher densities. The utility of such interpolations is driven by the present inability to calculate the dense matter equation of state in QCD from first principles. As we review, the parameters of effective quark models -- which have direct relevance to the more general structure of the QCD phase diagram of dense and hot matter -- are constrained by neutron star mass and radii measurements, in particular favoring large repulsive density-density and attractive diquark pairing interactions. We describe the structure of neutron stars constructed from the unified equations of states with crossover. Lastly we present the current equations of state -- called "QHC18" for quark-hadron crossover -- in a parametrized form practical for neutron star modeling.

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

We investigate the acceleration of electrons and positrons by magnetic-field-aligned electric fields in the polar funnel of an accreting black hole (BH). Applying the pulsar outer-gap theory to BH magnetospheres, we find that such a lepton accelerator arises in the immediate vicinity of the event horizon due to frame-dragging, and that their gamma-ray luminosity increases with decreasing accretion rate. Furthermore, we demonstrate that the gamma-ray flux is enhanced along the rotation axis by more than an order of magnitude if the BH spin increases from $a=0.90M$ to $a=0.9999M$. As a result, if a ten-solar-mass, almost-maximally rotating BH is located within 3 kpc, when its accretion rate is between 0.005% and 0.01% of the Eddington rate, its high-energy flare becomes detectable with the Fermi/Large Area Telescope, provided that the flare lasts longer than 1.2 months and that we view the source nearly along the rotation axis. In addition, its very-high-energy flux is marginally detectable with the Cherenkov Telescope Array, provided that the flare lasts longer than a night and that our viewing angle is about 45 degrees with respect to the rotation axis.

We report new spectroscopic observations obtained with the Michigan/Magellan Fiber System of 308 red giants (RGs) located in two fields near the photometric center of the bar of the Large Magellanic Cloud. This sample consists of 131 stars observed in previous studies (in one field) and 177 newly-observed stars (in the second field) selected specifically to more reliably establish the metallicity and age distributions of the bar. For each star, we measure its heliocentric line-of-sight velocity, surface gravity and metallicity from its high-resolution spectrum (effective temperatures come from photometric colors). The spectroscopic Hertzsprung-Russell diagrams---modulo small offsets in surface gravities---reveal good agreement with model isochrones. The mean metallicity of the 177-RG sample is $\rm [Fe/H]=-0.76\pm0.02$ with a metallicity dispersion $\sigma=0.28\pm0.03$. The corresponding metallicity distribution---corrected for selection effects---is well fitted by two Gaussian components: one metal-rich with a mean $-0.66\pm0.02$ and a standard deviation $0.17\pm0.01$, and the other metal-poor with $-1.20\pm0.24$ and $0.41\pm0.06$. The metal-rich and metal-poor populations contain approximately 85% and 15% of stars, respectively. We also confirm the velocity dispersion in the bar center decreases significantly from $31.2\pm4.3$ to $18.7\pm1.9$ km s$^{-1}$ with increasing metallicity over the range $-2.09$ to $-0.38$. Individual stellar masses are estimated using the spectroscopic surface gravities and the known luminosities. We find that lower mass hence older RGs have larger metallicity dispersion and lower mean metallicity than the higher-mass, younger RGs. The estimated masses, however, extend to implausibly low values ($\rm \sim 0.1~M_{\odot}$) making it impossible to obtain an absolute age-metallicity or age distribution of the bar.

Aluminum lumped-element kinetic inductance detectors (LEKIDs) sensitive to millimeter-wave photons have been shown to exhibit high quality factors, making them highly sensitive and multiplexable. The superconducting gap of aluminum limits aluminum LEKIDs to photon frequencies above 100 GHz. Manganese-doped aluminum (Al-Mn) has a tunable critical temperature and could therefore be an attractive material for LEKIDs sensitive to frequencies below 100 GHz if the internal quality factor remains sufficiently high when manganese is added to the film. To investigate, we measured some of the key properties of Al-Mn LEKIDs. A prototype eight-element LEKID array was fabricated using a 40 nm thick film of Al-Mn deposited on a 500 {\mu}m thick high-resistivity, float-zone silicon substrate. The manganese content was 900 ppm, the measured $T_c = 694\pm1$ mK, and the resonance frequencies were near 150 MHz. Using measurements of the forward scattering parameter $S_{21}$ at various bath temperatures between 65 and 250 mK, we determined that the Al-Mn LEKIDs we fabricated have internal quality factors greater than $2 \times 10^5$, which is high enough for millimeter-wave astrophysical observations. In the dark conditions under which these devices were measured, the fractional frequency noise spectrum shows a shallow slope that depends on bath temperature and probe tone amplitude, which could be two-level system noise. The anticipated white photon noise should dominate this level of low-frequency noise when the detectors are illuminated with millimeter-waves in future measurements. The LEKIDs responded to light pulses from a 1550 nm light-emitting diode, and we used these light pulses to determine that the quasiparticle lifetime is 60 {\mu}s.

We develop a method to identify the spectroscopic signature of unresolved L-dwarf ultracool companions, which compares the spectra of candidates and their associated control stars using spectral ratio differences and residual spectra. We present SpeX prism-mode spectra (0.7-2.5 micron) for a pilot sample of 111 mid M dwarfs, including 28 that were previously identified as candidates for unresolved ultracool companionship (a sub-sample from Cook et al. 2016; paper 1) and 83 single M dwarfs that were optically colour-similar to these candidates (which we use as `control stars'). We identify four candidates with evidence for near-infrared excess. One of these (WISE J100202.50+074136.3) shows strong evidence for an unresolved L dwarf companion in both its spectral ratio difference and its residual spectra, two most likely have a different source for the near-infrared excess, and the other may be due to spectral noise. We also establish expectations for a null result (i.e. by searching for companionship signatures around the M dwarf control stars), as well as determining the expected outcome for ubiquitous companionship (as a means of comparison with our actual results), using artificially generated unresolved M+L dwarf spectra. The results of these analyses are compared to those for the candidate sample, and reasonable consistency is found. With a full follow-up programme of our candidates sample from Cook et al., we might expect to confirm up to 40 such companions in the future, adding extensively to the known desert population of M3-M5 dwarfs.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

The task of flux calibration for LAMOST (Large sky Area Multi-Object Spectroscopic Telescope) spectra is difficult due to many factors. For example, the lack of standard stars, flat fielding for large field of view, and variation of reddening between different stars especially at low galactic latitudes etc. Poor selection, bad spectral quality, or extinction uncertainty of standard stars not only might induce errors to the calculated spectral response curve (SRC), but also might lead to failures in producing final 1D spectra. In this paper, we inspected spectra with Galactic latitude |b|>=60 degree and reliable stellar parameters, determined through the LAMOST Stellar Parameter Pipeline (LASP), to study the stability of the spectrograph. To guarantee the selected stars had been observed by each fiber, we selected 37,931 high quality exposures of 29,000 stars from LAMOST DR2, and more than 7 exposures for each fiber. We calculated the SRCs for each fiber for each exposure, and calculated the statistics of SRCs for spectrographs with both the fiber variations and time variations. The result shows that the average response curve of each spectrograph (henceforth ASPSRC) is relatively stable with statistical errors <= 10%. From the comparison between each ASPSRC and the SRCs for the same spectrograph obtained by 2D pipeline, we find that the ASPSRCs are good enough to use for the calibration. The ASPSRCs have been applied to spectra which were abandoned by LAMOST 2D pipeline due to the lack of standard stars, increasing the number of LAMOST spectra by 52,181 in DR2. Comparing those same targets with SDSS, the relative flux differences between SDSS spectra and that of LAMOST spectra with the ASPSRC method are less than 10%, which underlines that the ASPSRC method is feasible for LAMOST flux calibration.

A sample of the LAMOST spectra of the early type M0-M3 dwarfs is compared with the Kepler observations. It is found that M dwarfs with strong chromospheric emission in $H_{\alpha}$ have large flare activity in general. The rotational periods derived from the Kepler measurements have close correlations with the sizes of the flares, the power-law distribution index and the equivalent widths of the $H_{\alpha}$ emission. A clear trend exists for higher magnetic activities being detected in faster rotating M dwarfs (rotation periods $<$ 20 day).

The mapping of dark matter clustering from real space to redshift space introduces the anisotropic property to the measured density power spectrum in redshift space, known as the redshift space distortion effect. The mapping formula is intrinsically non-linear, which is complicated by the higher order polynomials due to indefinite cross correlations between the density and velocity fields, and the Finger-of-God effect due to the randomness of the peculiar velocity field. Whilst the full higher order polynomials remain unknown, the other systematics can be controlled consistently within the same order truncation in the expansion of the mapping formula, as shown in this paper. The systematic due to the unknown non-linear density and velocity fields is removed by separately measuring all terms in the expansion directly using simulations. The uncertainty caused by the velocity randomness is controlled by splitting the FoG term into two pieces, 1) the "one-point" FoG term being independent of the separation vector between two different points, and 2) the "correlated" FoG term appearing as an indefinite polynomials which is expanded in the same order as all other perturbative polynomials. Using 100 realizations of simulations, we find that the Gaussian FoG function with only one scale-independent free parameter works quite well, and that our new mapping formulation accurately reproduces the observed 2-dimensional density power spectrum in redshift space at the smallest scales by far, up to $k\sim 0.2h$Mpc, considering the resolution of future experiments.

We report on the development of scalable prototype microwave kinetic inductance detector (MKID) arrays tailored for future multi-kilo-pixel experiments that are designed to simultaneously characterize the polarization properties of both the cosmic microwave background (CMB) and Galactic dust emission. These modular arrays are composed of horn-coupled, polarization-sensitive MKIDs, and each pixel has four detectors: two polarizations in two spectral bands between 125 and 280 GHz. A horn is used to feed each array element, and a planar orthomode transducer, composed of two waveguide probe pairs, separates the incoming light into two linear polarizations. Diplexers composed of resonant-stub band-pass filters separate the radiation into 125 to 170 GHz and 190 to 280 GHz pass bands. The millimeter-wave power is ultimately coupled to a hybrid co-planar waveguide microwave kinetic inductance detector using a novel, broadband circuit developed by our collaboration. Electromagnetic simulations show the expected absorption efficiency of the detector is approximately 90%. Array fabrication will begin in the summer of 2016.

The anisotropic galaxy clustering of large scale structure observed by the Baryon Oscillation Spectroscopic Survey Data Release 11 is analyzed to probe the sum of neutrino mass in the small $m_\nu < 1eV$ limit in which the early broadband shape determined before the last scattering surface is immune from the variation of $m_\nu$. The signature of $m_\nu$ is imprinted on the altered shape of the power spectrum at later epoch, which provides an opportunity to access the non--trivial $m_\nu$ through the measured anisotropic correlation function in redshift space (hereafter RSD instead of Redshift Space Distortion). The non--linear RSD corrections with massive neutrinos in the quasi linear regime are approximately estimated using one-loop order terms computed by tomographic linear solutions. We suggest a new approach to probe $m_\nu$ simultaneously with all other distance measures and coherent growth functions, exploiting this deformation of the early broadband shape of the spectrum at later epoch. If the origin of cosmic acceleration is unknown, $m_\nu$ is poorly determined after marginalising over all other observables. However, we find that the measured distances and coherent growth functions are minimally affected by the presence of mild neutrino mass. Although the standard model of cosmic acceleration is assumed to be the cosmological constant, the constraint on $m_\nu$ is little improved. Interestingly, the measured CMB distance to the last scattering surface sharply slices the degeneracy between the matter content and $m_\nu$, and the hidden $m_\nu$ is excavated to be $m_\nu=0.19^{+0.28}_{-0.17} eV$ which is different from massless neutrino more than 68% confidence.

Many previous studies have shown that magnetic fields as well as sunspot structures present rapid and irreversible changes associated with solar flares. In this paper we first use five X-class flares observed by SDO/HMI to show that not only the magnetic fields and sunspot structures do show rapid, irreversible changes but also these changes are closely related, both spatially and temporally. The magnitudes of the correlation coefficients between the temporal variations of horizontal magnetic field and sunspot intensity are all larger than 0.90, with a maximum value of 0.99 and an average value of 0.96. Then using four active regions in quiescent times, three observed and one simulated, we show that in sunspot penumbra regions there also exists a close correlation between sunspot intensity and horizontal magnetic field strength, in addition to the well-known one between sunspot intensity and normal magnetic field strength. Connecting these two observational phenomena, we show that the sunspot structure change and the magnetic field change are the two facets of the same phenomena of solar flares, one change might be induced by the change of the other due to a linear correlation between sunspot intensity and magnetic field strength out of a local force balance.

We study the singularity of the congruences for both timelike and null geodesic curves using the expansion of the early anisotropic Bianchi type I Universe. In this paper, we concentrate on the influence of the shear of the timelike and null geodesic congruences in the early Universe. Under some natural conditions, we derive the Raychaudhuri type equation for the expansion and the shear-related equations. Recently, scientists working on the LIGO (Laser Interferometer Gravitational-Wave Observatory) have shown many possibilities to observing the anisotropy of the primordial gravitational wave background radiation. We deduce the evolution equation for the shear that may be responsible for those observational results.

The cosmic distance can be precisely determined using a `standard ruler' imprinted by primordial baryon acoustic oscillation (hereafter BAO) in the early Universe. The BAO at the targeted epoch is observed by analysing galaxy clustering in redshift space (hereafter RSD) for which a theoretical formulation is not yet fully understood, and thus makes this methodology unsatisfactory. The BAO analysis following a full RSD modelling is contaminated by systematic uncertainties due to a non-linear smearing effects such as non-linear corrections and by random viral velocity of galaxies. However, the BAO can be probed independently of RSD contamination using the BAO peak positions located in the 2D anisotropic correlation function. A new methodology is presented to measure peak positions, to test whether it is also contaminated by the same systematics in RSD, and to provide the radial and transverse cosmic distances determined by the 2D BAO peak positions. We find that in our model independent anisotropic clustering analysis we can obtain about $2\%$ and $5\%$ constraints on $D_A$ and $H^{-1}$ respectively with current BOSS data, which is competitive with other analysis.

The ultra-faint dwarf galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high spectral resolution using the Michigan/Magellan Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high levels of r-process material in Ret 2 (mean [Eu/Fe]=+1.69+/-0.05) found in three of these stars (mean [Fe/H]=-2.88+/-0.10). The abundances closely match the r-process pattern found in the well-studied metal-poor halo star CS22892-052. Such r-process-enhanced stars have not been found in any other ultra-faint dwarf galaxy, though their existence has been predicted by at least one model. The fourth star in Ret 2 ([Fe/H]=-3.42+/-0.20) contains only trace amounts of Sr ([Sr/Fe]=-1.73+/-0.43) and no detectable heavier elements. One r-process enhanced star is also enhanced in C (natal [C/Fe]=+1.1). This is only the third such star known, which suggests that the nucleosynthesis sites leading to C and r-process enhancements are decoupled. The r-process-deficient star is enhanced in Mg ([Mg/Fe]=+0.81+/-0.14), and the other three stars show normal levels of alpha-enhancement (mean [Mg/Fe]=+0.34+/-0.03). The abundances of other alpha and Fe-group elements closely resemble those in ultra-faint dwarf galaxies and metal-poor halo stars, suggesting that the nucleosynthesis that led to the large r-process enhancements either produced no light elements or produced light-element abundance signatures indistinguishable from normal supernovae.

We delineate the properties of QCD matter at baryon density n_B=1-10n_0 (n_0: nuclear saturation density), through the construction of neutron star equations of state that satisfy the neutron star mass-radius constraints as well as physical conditions on the speed of sound. The QCD matter is described in the 3-window modeling: at n_B < 2n_0 purely nuclear matter; at n_B > 5n_0 percolated quark matter; and at 2n_0 < n_B < 5n_0 matter intermediate between these two which are constructed by interpolation. Using a schematic quark model with effective interactions inspired from hadron and nuclear physics, we analyze the strength of interactions necessary to describe observed neutron star properties. Our finding is that the interactions should remain as strong as in the QCD vacuum, indicating that gluons at n_B =1-10 n_0 remain non-perturbative even after quark matter formation.

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

Measurement of the local dark matter density plays an important role in both Galactic dynamics and dark matter direct detection experiments. However, the estimated values from previous works are far from agreeing with each other. In this work, we provide a well-defined observed sample with 1427 G \& K type main-sequence stars from the LAMOST spectroscopic survey, taking into account selection effects, volume completeness, and the stellar populations. We apply a vertical Jeans equation method containing a single exponential stellar disk, a razor thin gas disk, and a constant dark matter density distribution to the sample, and obtain a total surface mass density of $\rm {78.7 ^{+3.9}_{-4.7}\ M_{\odot}\ pc^{-2}}$ up to 1 kpc and a local dark matter density of $0.0159^{+0.0047}_{-0.0057}\,\rm M_{\odot}\,\rm pc^{-3}$. We find that the sampling density (i.e. number of stars per unit volume) of the spectroscopic data contributes to about two-thirds of the uncertainty in the estimated values. We discuss the effect of the tilt term in the Jeans equation and find it has little impact on our measurement. Other issues, such as a non-equilibrium component due to perturbations and contamination by the thick disk population, are also discussed.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

The apparent anisotropies of the galaxy clustering in observable redshift space provide a unique opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales via the Alcock--Paczynski effect and redshift-space distortions. While the improved theoretical models have been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear scales, the applicability of these models is still limited in the presence of the non--perturbative smearing effect caused by the randomness of the relative velocities. Although the cosmological constraint from the anisotropic clustering will be certainly improved with a more elaborate theoretical model, we here consider an alternative way by using the statistical power of both the power spectrum and bispectrum at large scales. Based on the Fisher matrix analysis, we estimate the benefit of combining the power spectra and bispectra, finding that the constraints on the cosmic expansion and growth of structure will be improved by a factor of two. This compensates for the loss of constraining power using the power spectrum alone due to the randomness of the relative velocities.

We have identified outflows and bubbles in the Taurus molecular cloud based on the $\sim 100$ deg$^2$ Five College Radio Astronomy Observatory $^{12}$CO(1-0) and $^{13}$CO(1-0) maps and the Spitzer young stellar object catalogs. In the main 44 deg$^2$ area of Taurus we found 55 outflows, of which 31 were previously unknown. We also found 37 bubbles in the entire 100 deg$^2$ area of Taurus, all of which had not been found before. The total kinetic energy of the identified outflows is estimated to be $\bf \sim 3.9 \times 10^{45}$ erg, which is \textbf{1\%} of the cloud turbulent energy. The total kinetic energy of the detected bubbles is estimated to be $\sim 9.2 \times 10^{46}$ erg, which is 29\% of the turbulent energy of Taurus. The energy injection rate from outflows is $\bf \sim 1.3 \times 10^{33}~\rm erg\ s^{-1}$, \textbf{0.4 - 2 times} the dissipation rate of the cloud turbulence. The energy injection rate from bubbles is $\sim 6.4 \times 10^{33}$ erg s$^{-1}$, \textbf{2 - 10 times} the turbulent dissipation rate of the cloud. The gravitational binding energy of the cloud is $\bf \sim 1.5 \times 10^{48}$ {\bf erg}, \textbf{385} and 16 times the energy of outflows and bubbles, respectively. We conclude that neither outflows nor bubbles can \textbf{provide enough energy to balance the overall gravitational binding energy and the turbulent energy of Taurus. However,} in the current epoch, stellar feedback is sufficient to maintain the observed turbulence in Taurus.

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

It is generally believed that the evolution of magnetic helicity has a close relationship with solar activity. Before the launch of SDO, earlier studies have mostly used MDI/SOHO line of sight magnetograms and assumed that magnetic fields are radial when calculating magnetic helicity injection rate from photospheric magnetograms. However, this assumption is not necessarily true. Here we use the vector magnetograms and line of sight magnetograms, both taken by HMI/SDO, to estimate the effects of non-radial magnetic field on measuring magnetic helicity injection rate. We find that: 1) The effect of non-radial magnetic field on estimating tangential velocity is relatively small; 2) On estimating magnetic helicity injection rate, the effect of non-radial magnetic field is strong when active regions are observed near the limb and is relatively small when active regions are close to disk center; 3) The effect of non-radial magnetic field becomes minor if the amount of accumulated magnetic helicity is the only concern.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

We present a sample of metallic-line star (Am) candidates from the Large sky Area Multi-Object fiber Spectroscopic Telescope Data Release one (LAMOST DR1). According to the characteristic of under-abundance of calcium and overabundance of iron element of Am stars, we propose an empirical separation curve derived from line indices of Ca II K-line and iron lines we choose for low resolution spectra. 3537 Am candidates are ultimately selected from more than 30,000 stars which are classified as A-type or early-F stars by both LAMOST pipeline and visual inspection. Then we make some analysis on this sample and finally provide a list of these Am candidates with 10 relevant parameters. Comparing with other catalogues, Am candidates selected from LAMOST DR1 are much fainter on the whole. Obviously, our list is an important complementary to already known bright Am catalogues, and it offers valuable material for the research on this type of chemically peculiar stars.

We construct an equation of state for massive neutron stars based on quantum chromodynamics phenomenology. Our primary purpose is to delineate the relevant ingredients of equations of state that simultaneously have the required stiffness and satisfy constraints from thermodynamics and causality. These ingredients are: (i) a repulsive density-density interaction, universal for all flavors; (ii) the color-magnetic interaction active from low to high densities; (iii) confining effects, which become increasingly important as the baryon density decreases; (iv) non-perturbative gluons, which are not very sensitive to changes of the quark density. We use the following "3-window" description: At baryon densities below about twice normal nuclear density, 2n_0, we use the Akmal-Pandharipande-Ravenhall (APR) equation of state, and at high densities, > (4-7)n_0, we use the three-flavor Nambu-Jona-Lasinio (NJL) model supplemented by vector and diquark interactions. In the transition density region, we smoothly interpolate the hadronic and quark equations of state in the chemical potential-pressure plane. Requiring that the equation of state approach APR at low densities, we find that the quark pressure in non-confining models can be larger than the hadronic pressure, unlike in conventional equations of state. We show that consistent equations of state of stiffness sufficient to allow massive neutron stars are reasonably tightly constrained, suggesting that gluon dynamics remains non-perturbative even at baryon densities ~10n_0.

We analyze the clustering of large scale structure in the Universe in a model independent method, accounting for anisotropic effects along and transverse to the line of sight. The Baryon Oscillation Spectroscopy Survey Data Release 11 provides a large sample of 690,000 galaxies, allowing determination of the Hubble expansion H, angular distance D_A, and growth rate G_T at an effective redshift of z=0.57. After careful bias and convergence studies of the effects from small scale clustering, we find that cutting transverse separations below 40 Mpc/h delivers robust results while smaller scale data leads to a bias due to unmodelled nonlinear and velocity effects. The converged results are in agreement with concordance LCDM cosmology, general relativity, and minimal neutrino mass, all within the 68% confidence level. We also present results separately for the northern and southern hemisphere sky, finding a slight tension in the growth rate -- potentially a signature of anisotropic stress, or just covariance with small scale velocities -- but within 68% CL.

Baryon acoustic oscillations (BAO), known as one of the largest cosmological objects, is now recognized as standard cosmological tool to measure geometric distances via the Alcock-Paczynski effect, by which the observed BAO exhibits characteristic anisotropies in addition to the redshift distortions. This implies that once we know the correct distances to the observed BAO, the tip points of baryon acoustic peaks in the anisotropic correlation function of galaxies, $\xi(\sigma,\pi)$, can form a great circle (hereafter 2D BAO circle) in the $\sigma$ and $\pi$ plane, where $\sigma$ and $\pi$ are the separation of galaxy pair parallel and perpendicular to the line-of-sight, respectively. This 2D BAO circle remains unchanged under the variations of the unknown galaxy bias and/or coherent motion, while it varies transversely and radially with respect to the variations of $D_A$ and $H^{-1}$, respectively. Hereby the ratio between transverse distance $D_A$ and the radial distance $H^{-1}$ reproduces the intrinsic shape of 2D BAO circle, which is {\it a priori} given by the known broadband shape of spectra. All BAO peaks of $\xi(\sigma,\pi)$ are precisely calculated with the improved theoretical model of redshift distortion. We test this broadband Alcock--Paczynski method using BOSS--like mock catalogues. The transverse and radial distances are probed in precision of several percentage fractional errors, and the coherent motion is observed to match with the fiducial values accurately.

Our observations of the Universe are fundamentally anisotropic, with data from galaxies separated transverse to the line of sight coming from the same epoch while that from galaxies separated parallel to the line of sight coming from different times. Moreover, galaxy velocities along the line of sight change their redshift, giving redshift space distortions. We perform a full two-dimensional anisotropy analysis of galaxy clustering data, fitting in a substantially model independent manner the angular diameter distance D_A, Hubble parameter H, and growth rate ddelta/dln a without assuming a dark energy model. The results demonstrate consistency with LCDM expansion and growth, hence also testing general relativity. We also point out the interpretation dependence of the effective redshift z_eff, and its cosmological impact for next generation surveys.

We study microlensing light curves by a triple lens, in particular, by a primary star plus two planets. A four-fold degeneracy is confirmed in the light curves, similar to the close and wide degeneracy found in a double lens. Furthermore, we derive a set of equations for triple-lens in the external shear approximation. By using these external shear equations, we identify two kinds of continuous degeneracies which may confused double and triple lenses, i.e. the continuous external shear degeneracy among triple-lens systems and the double-triple lens degeneracy. These degeneracies are particularly important in high magnification events, and thus some caution needs to be applied when one infers the fraction of stars hosting multiple planets from microlensing. We study the dependence of the degeneracies on the lensing parameters (e.g., source trajectory) and give recipes about how the degeneracies should be explored with real data.

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

We present a set of white dwarf-main sequence (WDMS) binaries identified spectroscopically from the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST, also called the Guo Shou Jing Telescope) pilot survey. We develop a color selection criteria based on what is so far the largest and most complete Sloan Digital Sky Survey (SDSS) DR7 WDMS binary catalog and identify 28 WDMS binaries within the LAMOST pilot survey. The primaries in our binary sample are mostly DA white dwarfs except for one DB white dwarf. We derive the stellar atmospheric parameters, masses, and radii for the two components of 10 of our binaries. We also provide cooling ages for the white dwarf primaries as well as the spectral types for the companion stars of these 10 WDMS binaries. These binaries tend to contain hot white dwarfs and early-type companions. Through cross-identification, we note that nine binaries in our sample have been published in the SDSS DR7 WDMS binary catalog. Nineteen spectroscopic WDMS binaries identified by the LAMOST pilot survey are new. Using the 3$\sigma$ radial velocity variation as a criterion, we find two post-common-envelope binary candidates from our WDMS binary sample.

We present a spectroscopic catalog of 67082 M dwarfs from the LAMOST pilot survey. For each spectrum of the catalog, spectral subtype, radial velocity, equivalent width of H${\alpha}$, a number of prominent molecular band indices and the metal sensitive parameter $\zeta$ are provided . Spectral subtype have been estimated by a remedied Hammer program (Original Hammer: Covey et al. 2007), in which indices are reselected to obtain more accurate auto-classified spectral subtypes. All spectra in this catalog have been visually inspected to confirm the spectral subtypes. Radial velocities have been well measured by our developed program which uses cross-correlation method and estimates uncertainty of radial velocity as well. We also examine the magnetic activity properties of M dwarfs traced by H${\alpha}$ emission line. The molecular band indices included in this catalog are temperature or metallicity sensitive and can be used for future analysis of the physical properties of M dwarfs. The catalog is available on the website \url{http://sciwiki.lamost.org/MCatalogPilot/}.

The cosmological power spectrum of the coherent matter flow is measured exploiting an improved prescription for the apparent anisotropic clustering pattern in redshift space. New statistical analysis is presented to provide an optimal observational platform to link the improved redshift distortion theoretical model to future real datasets. The statistical power as well as robustness of our method are tested against 60 realizations of 8 Gpc/h^3 dark matter simulation maps mocking the precision level of upcoming wide--deep surveys. We showed that we can accurately extract the velocity power spectrum up to quasi linear scales of k~0.1 h/Mpc at z = 0.35 and up to k~0.15 h/Mpc at higher redshifts within a couple of percentage precision level. Our understanding of redshift space distortion is proved to be appropriate for precision cosmology, and our statistical method will guide us to righteous path to meet the real world.

We investigate the direct determination of expansion history using redshift distortions without plugging into detailed cosmological parameters. The observed spectra in redshift space include a mixture of information: fluctuations of density-density and velocity-velocity spectra, and distance measures of perpendicular and parallel components to the line of sight. Unfortunately it is hard to measure all the components simultaneously without any specific prior assumption. Common prior assumptions include a linear/quasi-linear model of redshift distortions or a model for the shape of the power spectra, which eventually breaks down on small scales at later epochs where nonlinear structure formation disturbs coherent growth. The degeneracy breaking, between the effect of cosmic distances and redshift distortions for example, depends on the prior we assume. An alternative approach is to utilize the cosmological principle inscribed in the heart of the Friedmann-Lematre-Robertson-Walker (hereafter FLRW) universe, that is, the specific relation between the angular diameter distance and the Hubble parameter, in this degeneracy breaking. We show that utilizing this FLRW prior early in the step of distinguishing the distance effect from redshift distortions helps us improve the detectability of power spectra and distance measures with no leaning on a combination of other experiments.

We show that in the single component situation all perturbation variables in the comoving gauge are conformally invariant to all perturbation orders. Generally we identify a special time slicing, the uniform-conformal transformation slicing, where all perturbations are again conformally invariant to all perturbation orders. We apply this result to the delta N formalism, and show its conformal invariance.

New statistical method is proposed to coherently combine Baryon Acoustic Oscillation statistics (BAO) and peculiar velocity measurements exploiting decomposed density-density and velocity-velocity spectra in real space from the observed redshift distortions in redshift space, 1) to achieve stronger dark energy constraints, \sigma(w)=0.06 and \sigma(w_a)=0.20, which are enhanced from BAO or velocity measurements alone, and 2) to cross-check consistency of dark energy constraints from two different approaches; BAO as geometrical measurements and peculiar velocity as large scale structure formation observables. In addition to those advantages, as power spectra decomposition procedure is free from uncertainty of galaxy bias, this simultaneous fitting is an optimal method to extract cosmological parameters without any pre-assumption about galaxy bias.

We present new measurements of the coherent motion of galaxies based on observations of the large-scale redshift-space distortions seen in the two-dimensional two-point correlation function of Luminous Red Galaxies in Data Release Seven of the Sloan Digital Sky Survey. We have developed a new methodology for estimating these coherent motions, which is less dependent on the details of galaxy bias and of the cosmological model to explain the late-time acceleration of the expansion of the Universe. We measure a one-dimensional velocity dispersion of galaxies on large-scales of \sigma_v=3.01^{+0.45}_{-0.46} Mpc/h and \sigma_v=3.69^{+0.47}_{-0.47} \mpcoh$ at a mean redshift of z=0.25 and 0.38 respectively. These values are fully consistent with predictions for a WMAP7-normalised LCDM Universe and inconsistent at confidence of 3.8\sigma with a Dvali-Gabadadze-Porrati (DGP) model for the Universe. We can convert the units of these $\sigma_v$ measurements to 270^{+40}_{-41} km/s and 320^{+41}_{-41} km/s respectively (assuming a $\Lambda$CDM universe), which are lower that expected based on recent low redshift (z<0.2) measurements of the peculiar velocity field (or "bulk flows"). It is difficult to directly compare these measurements as they cover different redshift ranges and different areas of the sky. However, one possible cosmological explanation for this discrepancy is that our Galaxy is located in unusually over, or under, dense region of the Universe.

Using the latest physical modeling and constrained by the most recent data, we develop a phenomenological parameterized model of the contributions to intensity and polarization maps at millimeter wavelengths from external galaxies and Sunyaev-Zeldovich effects. We find such modeling to be necessary for estimation of cosmological parameters from Planck data. For example, ignoring the clustering of the infrared background would result in a bias in n_s of 7 sigma. We show that the simultaneous marginalization over a full foreground model can eliminate such biases, while increasing the statistical uncertainty in cosmological parameters by less than 20%. The small increases in uncertainty can be significantly reduced with the inclusion of higher-resolution ground-based data. The multi-frequency analysis we employ involves modeling 46 total power spectra and marginalization over 17 foreground parameters. We show that we can also reduce the data to a best estimate of the CMB power spectra, and just two principal components (with constrained amplitudes) describing residual foreground contamination.

We explore the complementarity of weak lensing and galaxy peculiar velocity measurements to better constrain modifications to General Relativity. We find no evidence for deviations from GR on cosmological scales from a combination of peculiar velocity measurements (for Luminous Red Galaxies in the Sloan Digital Sky Survey) with weak lensing measurements (from the CFHT Legacy Survey). We provide a Fisher error forecast for a Euclid-like space-based survey including both lensing and peculiar velocity measurements, and show that the expected constraints on modified gravity will be at least an order of magnitude better than with present data, i.e. we will obtain 5% errors on the modified gravity parametrization described here. We also present a model--independent method for constraining modified gravity parameters using tomographic peculiar velocity information, and apply this methodology to the present dataset.

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

We test General Relativity (GR) using current cosmological data: the cosmic microwave background (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs-Wolfe (ISW) effect from the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation of supernovae Type Ia (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the weak lensing (WL) data from CFHTLS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear and mildly non-linear scales. We first test a model where the effective Newton's constant, mu, and the ratio of the two gravitational potentials, eta, transit from the GR value to another constant at late times; in this case, we find that standard GR is fully consistent with the combined data. The strongest constraint comes from the ISW effect which would arise from this gravitational transition; the observed ISW signal imposes a tight constraint on a combination of mu and eta that characterizes the lensing potential. Next, we consider four pixels in time and space for each function mu and eta, and perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes are consistent with GR within the errors. Only one eigenmode shows a 2-sigma deviation from the GR prediction, which is likely to be due to a systematic effect. However, the detection of such a deviation demonstrates the power of our time- and scale-dependent PCA methodology when combining observations of structure formation and expansion history to test GR.

We develop an optimized technique to extract density--density and velocity--velocity spectra out of observed spectra in redshift space. The measured spectra of the distribution of halos from redshift distorted mock map are binned into 2--dimensional coordinates in Fourier space so as to be decomposed into both spectra using angular projection dependence. With the threshold limit introduced to minimize nonlinear suppression, the decomposed velocity--velocity spectra are reasonably well measured up to scale k=0.07 h/Mpc, and the measured variances using our method are consistent with errors predicted from a Fisher matrix analysis. The detectability is extendable to k\sim 0.1 h/Mpc with more conservative bounds at the cost of weakened constraint.

We present here a new parameterization for the bulk motions of galaxies and clusters (in the linear regime) that can be measured statistically from the shape and amplitude of the two-dimensional two-point correlation function. We further propose the one-dimensional velocity dispersion (v_p) of the bulk flow as a complementary measure of redshift-space distortions, which is model-independent and not dependent on the normalisation method. As a demonstration, we have applied our new methodology to the C4 cluster catalogue constructed from Data Release Three (DR3) of the Sloan Digital Sky Survey. We find v_p=270^{+433}km/s (also consistent with v_p=0) for this cluster sample (at z=0.1), which is in agreement with that predicted for a WMAP5-normalised LCDM model (i.e., v_p(LCDM=203km/s). This measurement does not lend support to recent claims of excessive bulk motions (\simeq1000 km/s) which appear in conflict with LCDM, although our large statistical error cannot rule them out. From the measured coherent evolution of v_p, we develop a technique to re-construct the perturbed potential, as well as estimating the unbiased matter density fluctuations and scale--independent bias.

Next generation surveys will observe the large-scale structure of the Universe with unprecedented accuracy. This will enable us to test the relationships between matter over-densities, the curvature perturbation and the Newtonian potential. Any large-distance modification of gravity or exotic nature of dark energy modifies these relationships as compared to those predicted in the standard smooth dark energy model based on General Relativity. In linear theory of structure growth such modifications are often parameterised by virtue of two functions of space and time that enter the relation of the curvature perturbation to, first, the matter over-density, and second, the Newtonian potential. We investigate the predictions for these functions in Brans-Dicke theory, clustering dark energy models and interacting dark energy models. We find that each theory has a distinct path in the parameter space of modified growth. Understanding these theoretical priors on the parameterisations of modified growth is essential to reveal the nature of cosmic acceleration with the help of upcoming observations of structure formations.

Dark energy in General Relativity is typically non-interacting with other matter. However, it is possible that the dark energy interacts with the dark matter, and in this case, the dark matter can violate the universality of free fall (the weak equivalence principle). We show that some forms of the dark sector interaction do not violate weak equivalence. For those interactions that do violate weak equivalence, there are no available laboratory experiments to probe this violation for dark matter. But cosmology provides a test for violations of the equivalence principle between dark matter and baryons -- via a test for consistency of the observed galaxy velocities with the Euler equation.

Measuring the statistics of galaxy peculiar velocities using redshift-space distortions is an excellent way of probing the history of structure formation. Because galaxies are expected to act as test particles within the flow of matter, this method avoids uncertainties due to an unknown galaxy density bias. We show that the parameter combination measured by redshift-space distortions, $f\sigma_8^{\rm mass}$ provides a good test of dark energy models, even without the knowledge of bias or $\sigma_8^{\rm mass}$ required to extract $f$ from this measurement (here $f$ is the logarithmic derivative of the linear growth rate, and $\sigma_8^{\rm mass}$ is the root-mean-square mass fluctuation in spheres with radius $8h^{-1}$Mpc). We argue that redshift-space distortion measurements will help to determine the physics behind the cosmic acceleration, testing whether it is related to dark energy or modified gravity, and will provide an opportunity to test possible dark energy clumping or coupling between dark energy and dark matter. If we can measure galaxy bias in addition, simultaneous measurement of both the overdensity and velocity fields can be used to test the validity of equivalence principle, through the continuity equation.

Using the linear theory of perturbations in General Relativity, we express a set of consistency relations that can be observationally tested with current and future large scale structure surveys. We then outline a stringent model-independent program to test gravity on cosmological scales. We illustrate the feasibility of such a program by jointly using several observables like peculiar velocities, galaxy clustering and weak gravitational lensing. After addressing possible observational or astrophysical caveats like galaxy bias and redshift uncertainties, we forecast in particular how well one can predict the lensing signal from a cosmic shear survey using an over-lapping galaxy survey. We finally discuss the specific physics probed this way and illustrate how $f(R)$ gravity models would fail such a test.

Much of the science case for the next generation of deep, wide-field optical/infrared surveys has been driven by the further study of dark energy. This is a laudable goal (and the subject of a companion white paper by Zhan et al.). However, one of the most important lessons of the current generation of surveys is that the interesting science questions at the end of the survey are quite different than they were when the surveys were being planned. The current surveys succeeded in this evolving terrain by being very general tools that could be applied to a number of very fundamental measurements. Likewise, the accessibility of the data enabled the broader cosmological and astronomical community to generate more science than the survey collaborations could alone. With that in mind, we should consider some of the basic physical and cosmological questions that surveys like LSST and JDEM-Wide will be able to address.

We construct a consistency test of General Relativity (GR) on cosmological scales. This test enables us to distinguish between the two alternatives to explain the late-time accelerated expansion of the universe, that is, dark energy models based on GR and modified gravity models without dark energy. We derive the consistency relation in GR which is written only in terms of observables - the Hubble parameter, the density perturbations, the peculiar velocities and the lensing potential. The breakdown of this consistency relation implies that the Newton constant which governs large-scale structure is different from that in the background cosmology, which is a typical feature in modified gravity models. We propose a method to perform this test by reconstructing the weak lensing spectrum from measured density perturbations and peculiar velocities. This reconstruction relies on Poisson's equation in GR to convert the density perturbations to the lensing potential. Hence any inconsistency between the reconstructed lensing spectrum and the measured lensing spectrum indicates the failure of GR on cosmological scales. The difficulties in performing this test using actual observations are discussed.

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. In this paper we present a formalism for forecasting the constraints on the growth of structure which would arise in an idealized survey. This Fisher matrix based formalism can be used to study the power and aid in the design of future surveys.

We study a braneworld model in which a phantom-like behaviour occurs with only cold dark matter and a cosmological constant, due to a large distance modification of gravity. With the addition of curvature, the geometrical tests are not strict enough to rule out models in which gravity is modified significantly on large scales. We show that this degeneracy in the parameter space is broken by the structure formation tests, such as the integrated Sachs-Wolfe effect, which can probe general relativity on large scales.

The solar corona is a complex system, with nonisothermal plasma and being in the self-gravitating field of the Sun. So the corona plasma is not only a nonequilibrium system but also a nonextensive one. We estimate the parameter of describing the degree of nonextensivity of the corona plasma and study the generalization of the solar wind speed theory in the framework of nonextensive statistical mechanics. It is found that, when use Chapman's corona model (1957) as the radial distribution of the temperature in the corona, the nonextensivity reduces the gas pressure outward and thus leads a significant deceleration effect on the radial speed of the solar wind.

In this paper, we study the large scale structure formation of the normal branch in DGP model (Dvail, Gabadadze and Porrati brane world model) by applying the scaling method developed by Sawicki, Song and Hu for solving the coupled perturbed equations of motion of on-brane and off-brane. There is detectable departure of perturbed gravitational potential from LCDM even at the minimal deviation of the effective equation of state w_eff below -1. The modified perturbed gravitational potential weakens the integrated Sachs-Wolfe effect which is strengthened in the self-accelerating branch DGP model. Additionally, we discuss the validity of the scaling solution in the de Sitter limit at late times.

Models which accelerate the expansion of the universe through the addition of a function of the Ricci scalar f(R) leave a characteristic signature in the large-scale structure of the universe at the Compton wavelength scale of the extra scalar degree of freedom. We search for such a signature in current cosmological data sets: the WMAP cosmic microwave background (CMB) power spectrum, SNLS supernovae distance measures, the SDSS luminous red galaxy power spectrum, and galaxy-CMB angular correlations. Due to theoretical uncertainties in the nonlinear evolution of f(R) models, the galaxy power spectrum conservatively yields only weak constraints on the models despite the strong predicted signature in the linear matter power spectrum. Currently the tightest constraints involve the modification to the integrated Sachs-Wolfe effect from growth of gravitational potentials during the acceleration epoch. This effect is manifest for large Compton wavelengths in enhanced low multipole power in the CMB and anti-correlation between the CMB and tracers of the potential. They place a bound on the Compton wavelength of the field be less than of order the Hubble scale.

The self-accelerating braneworld model (DGP) can be tested from measurements of the expansion history of the universe and the formation of structure. Current constraints on the expansion history from supernova luminosity distances, the CMB, and the Hubble constant exclude the simplest flat DGP model at about 3sigma. The best-fit open DGP model is, however, only a marginally poorer fit to the data than flat LCDM. Its substantially different expansion history raises structure formation challenges for the model. A dark-energy model with the same expansion history would predict a highly significant discrepancy with the baryon oscillation measurement due the high Hubble constant required and a large enhancement of CMB anisotropies at the lowest multipoles due to the ISW effect. For the DGP model to satisfy these constraints new gravitational phenomena would have to appear at the non-linear and cross-over scales respectively. A prediction of the DGP expansion history in a region where the phenomenology is well understood is that high-redshift galaxies should be substantially correlated with the CMB through the ISW effect. This correlation should provide a sharp test of the DGP model in the future.

We develop a scaling ansatz for the master equation in Dvali, Gabadadze, Porrati cosmologies, which allows us to solve the equations of motion for perturbations off the brane during periods when the on-brane evolution is scale-free. This allows us to understand the behavior of the gravitational potentials outside the horizon at high redshifts and close to the horizon today. We confirm that the results of Koyama and Maartens are valid at scales relevant for observations such as galaxy-ISW correlation. At larger scales, there is an additional suppression of the potential which reduces the growth rate even further and would strengthen the ISW effect.

We study the evolution of linear cosmological perturbations in f(R) models of accelerated expansion in the physical frame where the gravitational dynamics are fourth order and the matter is minimally coupled. These models predict a rich and testable set of linear phenomena. For each expansion history, fixed empirically by cosmological distance measures, there exists two branches of f(R) solutions that are parameterized by B propto d^2 f/dR^2. For B<0, which include most of the models previously considered, there is a short-timescale instability at high curvature that spoils agreement with high redshift cosmological observables. For the stable B>0 branch, f(R) models can reduce the large-angle CMB anisotropy, alter the shape of the linear matter power spectrum, and qualitatively change the correlations between the CMB and galaxy surveys. All of these phenomena are accessible with current and future data and provide stringent tests of general relativity on cosmological scales.

With a wariness of Occam's razor awakened by the discovery of cosmic acceleration, we abandon the usual assumption of zero mean curvature and ask how well it can be determined by planned surveys. We also explore the impact of uncertain mean curvature on forecasts for the performance of planned dark energy probes. We find that weak lensing and photometric baryon acoustic oscillation data, in combination with CMB data, can determine the mean curvature well enough that the residual uncertainty does not degrade constraints on dark energy. We also find that determinations of curvature are highly tolerant of photometric redshift errors.

The large scale structure formation in the standard Einstein gravity at later time is uniquely determined by the expansion history H(z) measured by the geometrical factor. A possible departure from the expected standard structure formation is detectable by growth structure test. We design new deviation parameters to test standard gravity at large scales and estimate the expected errors by weak lensing experiment. The effective Newtonian constant G_N at the present time can be constrained around 10% accuracy by measuring the difference curvature perturbations (\Phi-\Psi)/2 within 2% accuracy.

The cosmic acceleration was discovered in one of the brane-based models. We are interested in discriminating this model from the dark energy by tomographic cosmic shear. Growth factors are different in the two models when one adjusts parameters to get nearly identical H(z). The two models could be distinguished with independent determinations of both geometrical factors and the growth factor. We introduce new parameterizations to separate the influence of geometry and the influence of growth on cosmic shear maps. We find that future observations will be able to distinguish between both models.

An alternative mechanism that dims high redshift supernovae without cosmic acceleration utilizes an oscillation of photons into a pseudo-scalar particle during transit. Since angular diameter distance measures are immune to the loss of photons, this ambiguity in interpretation can be resolved by combining CMB acoustic peak measurements with the recent baryon oscillation detection in galaxy power spectra. This combination excludes a non-accelerating dark energy species at the 4sigma level regardless of the level of the pseudo-scalar coupling. While solutions still exist with substantial non-cosmological dimming of supernovae, they may be tested with future improvement in baryon oscillation experiments.

Small distortions in the observed shapes of distant galaxies, a cosmic shear due to gravitational lensing, can be used to simultaneously determine the distance-redshift relation, r(z), and the density contrast growth factor, g(z). Both of these functions are sensitive probes of the acceleration. Their simultaneous determination allows for a consistency test and provides sensitivity to physics beyond the standard dark energy paradigm.

Weak lensing observations and supernova observations, combined with CMB observations, can both provide powerful constraints on dark energy properties. Considering statistical errors only, we find luminosity distances inferred from 2000 supernovae and large-scale (l < 1000) angular power spectra inferred from redshift-binned cosmic shear maps place complementary constraints on w_0 and w_a where w(z) = w_0 + w_a(a-1). Further, each set of observations can constrain higher-dimensional parameterizations of w(z); we consider eigenmodes of the w(z) error covariance matrix and find such datasets can each constrain the amplitude of about 5 w(z) eigenmodes. We also consider another parameterization of the dark energy.

We study how parameter error forecasts for tomographic cosmic shear observations are affected by sky coverage, density of source galaxies, inclusion of CMB experiments, simultaneou fitting of non--dark energy parameters, and the parametrization of the history of the dark energy equation-of-state parameter w(z). We find tomographic shear-shear power spectra on large angular scales (l<1000) inferred from all-sky observations, in combination with Planck, can achieve sigma(w0)=0.06 and sigma(wa)=0.09 assuming the equation-of-state parameter is given by w(z)=w0+wa(1-a(z)) and that nine other matter content and primordial power spectrum parameters are simultaneously fit. Taking parameters other than w0, wa and Omegam to be completely fixed by the CMB we find errors on w0 and wa that are only 10% and 30% better respectively, justifying this common simplifying assumption. We also study `dark energy tomography': reconstuction of w(z) assumed to be constant within each of five independent redshift bins. With smaller-scale information included by use of the Jain & Taylor ratio statistic we find sigma(wi)<0.1 for all five redshift bins and sigma(wi)<0.02 for both bins at z<0.8. Finally, addition of cosmic shear can also reduce errors on quantities already determined well by the CMB. We find the sum of neutrino masses can be determined to +-0.013eV and that the primordial power specrum power-law index, ns, as well as dns/dlnk, can be determined more than a factor of two better than by Planck alone. These improvements may be highly valuable since the lower bound on the sum of neutrino masses is 0.06eV as inferred from atmospheric neutrino oscillations, and slow-roll models of inflation predict non-zero dns/dlnk at the forecasted error levels when |n_S-1|>0.04.

An effective field theory technique that combines the standard nuclear physics approach and chiral perturbation theory is applied to the $hen$ process, ${}^{3}{He}+n\to {}^{4}{He}+ \gamma$. For the initial and final nuclear states, high-precision wave functions are generated via the variational Monte Carlo method using the Argonne $v_{14}$ potential and Urbana VIII trinucleon interactions, while the relevant transition operators are calculated up to ${\cal O}(Q^4)$ in HB$\chi$PT. The imposition of the renormalization condition that the magnetic moments of ${}^{3}{He}$ and ${}^{3}{H}$ be reproduced allows us to carry out a parameter-free calculation of the $hen$ cross section. The result, $\sigma=(60\pm 3\pm 1) \mu b$, is in reasonable agreement with the experimental values, $(54\pm 6) \mu b$ and $(55\pm 3) \mu b$. This agreement demonstrates the validity of the calculational method previously used for estimating the reaction rate of the solar $hep$ process.

Distortions of Cosmic Microwave Background (CMB) temperature and polarization maps caused by gravitational lensing, observable with high angular resolution and high sensitivity, can be used to measure the neutrino mass. Assuming two massless species and one with mass m_nu we forecast sigma(m_nu) = 0.15 eV from the Planck satellite and sigma(m_nu)=0.04 eV from observations with twice the angular resolution and about 20 times the sensitivity. A detection is likely at this higher sensitivity since the observation of atmospheric neutrino oscillations require mass-squared differences greater than about (0.04 eV)^2.

We study the detectability, given CMB polarization maps, of departures from the inflationary consistency equation, r \equiv T/S \simeq -5 n_T, where T and S are the tensor and scalar contributions to the quadrupole variance, respectively. The consistency equation holds if inflation is driven by a slowly-rolling scalar field. Departures can be caused by: 1) higher-order terms in the expansion in slow-roll parameters, 2) quantum loop corrections or 3) multiple fields. Higher-order corrections in the first two slow-roll parameters are undetectably small. Loop corrections are detectable if they are nearly maximal and r \ga 0.1. Large departures (|\Delta n_T| \ga 0.1) can be seen if r \ga 0.001. High angular resolution can be important for detecting non-zero r+5n_T, even when not important for detecting non-zero r.

Holographic considerations may provide a glimpse of quantum gravity beyond what is currently accessible by other means. Here we apply holography to inflationary cosmology. We argue that the appropriate holographic bound on the total entropy in the inflationary perturbations is given by the {\it difference} of the area of the apparent horizon in the slow-roll approximation and its exact value near the end of inflation. This implies that the effective field theory of the inflaton weakly coupled to gravity requires modifications below the Planck scale. For a conventional model of inflation this scale is of order of the GUT scale, $\Lambda_{UV} \la 10^{16} GeV$, but could be considerably lower. Signatures of such new physics could show up in the CMB.

The intervening large--scale structure distorts cosmic microwave background (CMB) anisotropies via gravitational lensing. The same large--scale structure, traced by dusty star--forming galaxies, also induces anisotropies in the far--infrared background (FIRB). We investigate the resulting inter--dependence of the FIRB and CMB with a halo model for the FIRB. In particular, we calculate the cross--correlation between the lensing potential and the FIRB. The lensing potential can be quadratically estimated from CMB temperature and/or polarization maps. We show that the cross--correlation can be measured with high signal--to--noise with data from the Planck Surveyor. We discuss how such a measurement can be used to understand the nature of FIRB sources and their relation to the distribution of dark matter.

We show that the polarization of the cosmic microwave background can be used to detect gravity waves from inflation if the energy scale of inflation is above 3.2 times 10^15 GeV. These gravity waves generate polarization patterns with a curl, whereas (to first order in perturbation theory) density perturbations do not. The limiting ``noise'' arises from the second--order generation of curl from density perturbations, or rather residuals from its subtraction. We calculate optimal sky coverage and detectability limits as a function of detector sensitivity and observing time.