Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The disagreement between direct late-time measurements of the Hubble constant from the SH0ES collaboration, and early-universe measurements based on the $\Lambda$CDM model from the Planck collaboration might, at least in principle, be explained by new physics in the early universe. Recently, the application of the Effective Field Theory of Large-Scale Structure to the full shape of the power spectrum of the SDSS/BOSS data has revealed a new, rather powerful, way to measure the Hubble constant and the other cosmological parameters from Large-Scale Structure surveys. In light of this, we analyze two models for early universe physics, Early Dark Energy and Rock 'n' Roll, that were designed to significantly ameliorate the Hubble tension. Upon including the information from the full shape to the Planck, BAO, and Supernovae measurements, we find that the degeneracies in the cosmological parameters that were introduced by these models are well broken by the data, so that these two models do not significantly ameliorate the tension.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

The scale of ongoing and future electromagnetic surveys pose formidable challenges to classify astronomical objects. Pioneering efforts on this front include citizen science campaigns adopted by the Sloan Digital Sky Survey (SDSS). SDSS datasets have been recently used to train neural network models to classify galaxies in the Dark Energy Survey (DES) that overlap the footprint of both surveys. Herein, we demonstrate that knowledge from deep learning algorithms, pre-trained with real-object images, can be transferred to classify galaxies that overlap both SDSS and DES surveys, achieving state-of-the-art accuracy $\gtrsim99.6\%$. We demonstrate that this process can be completed within just eight minutes using distributed training. While this represents a significant step towards the classification of DES galaxies that overlap previous surveys, we need to initiate the characterization of unlabelled DES galaxies in new regions of parameter space. To accelerate this program, we use our neural network classifier to label over ten thousand unlabelled DES galaxies, which do not overlap previous surveys. Furthermore, we use our neural network model as a feature extractor for unsupervised clustering and find that unlabeled DES images can be grouped together in two distinct galaxy classes based on their morphology, which provides a heuristic check that the learning is successfully transferred to the classification of unlabelled DES images. We conclude by showing that these newly labeled datasets can be combined with unsupervised recursive training to create large-scale DES galaxy catalogs in preparation for the Large Synoptic Survey Telescope era.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.

In this paper, we report new limits on 21cm emission from cosmic reionization based on a 135-day observing campaign with a 64-element deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. This work extends the work presented in Parsons et al. (2014) with more collecting area, a longer observing period, improved redundancy-based calibration, optimal fringe-rate filtering, and improved power-spectral analysis using optimal quadratic estimators. The result is a new $2\sigma$ upper limit on $\Delta^{2}(k)$ of (22.4 mK)$^2$ in the range $0.15 < k < 0.5h\ {\rm Mpc}^{-1}$ at $z = 8.4$. This represents a three-fold improvement over the previous best upper limit. As we discuss in more depth in a forthcoming paper (Pober et al. 2015, in prep), this upper limit supports and extends previous evidence against extremely cold reionization scenarios. We conclude with a discussion of implications for future 21cm reionization experiments, including the newly funded Hydrogen Epoch of Reionization Array (HERA). $\textbf{The limits presented in this paper have been retracted: The erratum can be found in Appendix A.}$

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.

The Hydrogen Epoch of Reionization Array (HERA) is a radio interferometer aiming to detect the power spectrum of 21 cm fluctuations from neutral hydrogen from the Epoch of Reionization (EOR). Drawing on lessons from the Murchison Widefield Array (MWA) and the Precision Array for Probing the Epoch of Reionization (PAPER), HERA is a hexagonal array of large (14 m diameter) dishes with suspended dipole feeds. Not only does the dish determine overall sensitivity, it affects the observed frequency structure of foregrounds in the interferometer. This is the first of a series of four papers characterizing the frequency and angular response of the dish with simulations and measurements. We focus in this paper on the angular response (i.e., power pattern), which sets the relative weighting between sky regions of high and low delay, and thus, apparent source frequency structure. We measure the angular response at 137 MHz using the ORBCOMM beam mapping system of Neben et al. We measure a collecting area of 93 m^2 in the optimal dish/feed configuration, implying HERA-320 should detect the EOR power spectrum at z~9 with a signal-to-noise ratio of 12.7 using a foreground avoidance approach with a single season of observations, and 74.3 using a foreground subtraction approach. Lastly we study the impact of these beam measurements on the distribution of foregrounds in Fourier space.

We present a new method for interferometric imaging that is ideal for the large fields of view and compact arrays common in 21 cm cosmology. We first demonstrate the method with simulations for two very different low frequency interferometers, the Murchison Widefield Array (MWA) and the MIT Epoch of Reionization (MITEoR) Experiment. We then apply the method to the MITEoR data set collected in July 2013 to obtain the first northern sky map from 128 MHz to 175 MHz at about 2 degree resolution, and find an overall spectral index of -2.73+/-0.11. The success of this imaging method bodes well for upcoming compact redundant low-frequency arrays such as HERA. Both the MITEoR interferometric data and the 150 MHz sky map are publicly available at http://space.mit.edu/home/tegmark/omniscope.html.

We examine the degree to which observations of large-scale cosmic microwave background (CMB) polarization can shed light on the puzzling large-scale power modulation in maps of CMB anisotropy. We consider a phenomenological model in which the observed anomaly is caused by modulation of large-scale primordial curvature perturbations, and calculate Fisher information and error forecasts for future polarization data, constrained by the existing CMB anisotropy data. Because a significant fraction of the available information is contained in correlations with the anomalous temperature data, it is essential to account for these constraints. We also present a systematic approach to finding a set of normal modes that maximize the available information, generalizing the well-known Karhunen-Loeve transformation to take account of the constraints from the temperature data. A polarization map covering at least $\sim 60\%$ of the sky should be able to provide a $3\sigma$ detection of modulation at the level favored by the temperature data. A significant fraction of the information in such a data set is contained in the single mode that optimally encapsulates the signal due to temperature-polarization correlation.

We present an improved Global Sky Model (GSM) of diffuse galactic radio emission from 10 MHz to 5 THz, whose uses include foreground modeling for CMB and 21 cm cosmology. Our model improves on past work both algorithmically and by adding new data sets such as the Planck maps and the enhanced Haslam map. Our method generalizes the Principal Component Analysis approach to handle non-overlapping regions, enabling the inclusion of 29 sky maps with no region of the sky common to all. We also perform a blind separation of our GSM into physical components with a method that makes no assumptions about physical emission mechanisms (synchrotron, free-free, dust, etc). Remarkably, this blind method automatically finds five components that have previously only been found "by hand", which we identify with synchrotron, free-free, cold dust, warm dust, and the CMB anisotropy, with maps and spectra agreeing with previous work but in many cases with smaller error bars. The improved GSM is available online at github.com/jeffzhen/gsm2016.

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology.

We demonstrate that the observation of old strange quark stars (SQSs) can set important limits on the scattering cross sections $\sigma_q$ between the light quarks and the non-interacting scalar dark matter (DM). By analyzing a set of 1403 of solitary pulsarlike compact stars in the Milky Way, we find the old solitary pulsar PSR J1801-0857D can set the most stringent upper limits on $\sigma_q$ or the DM-proton scattering cross sections $\sigma_p$. By converting $\sigma_q$ into $\sigma_p$ based on effective operator analyses, we show the resulting $\sigma_p$ limit by assuming PSR J1801-0857D to be a SQS could be comparable with that of the current direct detection experiments but much weaker (by several orders of magnitude) than that obtained by assuming PSR J1801-0857D to be a neutron star (NS), which requires an extremely small $\sigma_p$ far beyond the limits of direct detection experiments. Our findings imply that the old pulsars are favored to be SQSs rather than NSs if the scalar DM were observed by future terrestrial experiments.

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

In order to study the "Cosmic Dawn" and the Epoch of Reionization with 21 cm tomography, we need to statistically separate the cosmological signal from foregrounds known to be orders of magnitude brighter. Over the last few years, we have learned much about the role our telescopes play in creating a putatively foreground-free region called the "EoR window." In this work, we examine how an interferometer's effects can be taken into account in a way that allows for the rigorous estimation of 21 cm power spectra from interferometric maps while mitigating foreground contamination and thus increasing sensitivity. This requires a precise understanding of the statistical relationship between the maps we make and the underlying true sky. While some of these calculations would be computationally infeasible if performed exactly, we explore several well-controlled approximations that make mapmaking and the calculation of map statistics much faster, especially for compact and highly redundant interferometers designed specifically for 21 cm cosmology. We demonstrate the utility of these methods and the parametrized trade-offs between accuracy and speed using one such telescope, the upcoming Hydrogen Epoch of Reionization Array, as a case study.

Isospin-violating dark matter (IVDM), which couples differently with protons and neutrons, provides a promising mechanism to ameliorate the tension among recent direct detection experiments. Assuming DM is non-interacting bosonic asymmetric IVDM, we investigate how the existence of old neutron stars limits the DM-proton scattering cross-section $\sigma_{\rm p}$, especially the effects of the isospin violating DM-nucleon interactions and the symmetry energy in the equation of state of isospin asymmetric nuclear matter. Our calculations are completely based on general relativity and especially the structure of neutron stars is obtained by solving the Tolman-Oppenheimer-Volkoff equations with nuclear matter equation of state constrained by terrestrial experiments. We find that, by considering the more realistic neutron star model rather than a simple uniform neutron sphere as usual, the $\sigma_{\rm p}$ bounds from old neutron stars can be varied by more than an order of magnitude depending on the specific values of the DM neutron-to-proton coupling ratio $f_{\rm n}/f_{\rm p}$, and they can be further varied by more than a factor of two depending on the density dependence of the symmetry energy. In particular, we demonstrate that the observed nearby isolated old neutron star PSR B1257+12 can set a very strong limit on $\sigma_{\rm p}$ for low-mass DM particles ($ \le 20 \, {\rm GeV}$) that reaches a sensitivity beyond the current best limits from direct detection experiments and disfavors the DM interpretation of previously-reported positive experimental results, including the IVDM.

A general formula of the symmetry energy for many-body interaction is proposed and the commonly used two-body interaction symmetry energy is recovered. Within Landau's theory (Lt), we generalize two equations of state (EoS) CCS$\delta$3 and CCS$\delta$5 to asymmetric nuclear matter. We assume that the density and density difference between protons and neutrons divided by their sum are order parameters. We use different EoS to study neutron stars by solving the TOV equations. We demonstrate that different EoS give different mass and radius relation for neutron stars even when they have exactly the same ground state (gs) properties ($E/A$, $\rho_0$, $K$, $S$, $L$ and $K_{sym}$). Furthermore, for one EoS we change $K_{sym}$ and fix all the other gs parameters. We find that for some $K_{sym}$ the EoS becomes unstable at high density even for neutron matter. This suggests that a neutron star (NS) can exist below and above the instability region but in different states: a quark gluon plasma (QGP) at high density and baryonic matter at low density. If the star's central density is in the instability region, then we associate these conditions to the occurrence of Supernovae (SN).

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

Isospin-violating dark matter (IVDM) provides a possible mechanism to ameliorate the tension among recent direct detection experiments. For IVDM, we demonstrate that the results of direct detection experiments based on neutron-rich target nuclei may depend strongly on the density dependence of the symmetry energy which is presently largely unknown and controls the neutron skin thickness that reflects the relative difference of neutron and proton form factors in the neutron-rich nuclei. In particular, using the neutron and proton form factors obtained from Skyrme-Hartree-Fock calculations by varying the symmetry energy within the uncertainty region set by the latest model-independent measurement of the neutron skin thickness of $^{208}$Pb from PREX experiment at JLab, we find that, for IVDM with neutron-to-proton coupling ratio fixed to $f_n/f_p=-0.7$, the form factor effect may enhance the sensitivity of Xe-based detectors (e.g., XENON100 and LUX) to the DM-proton cross section by a factor of $3$ in the DM mass region constrained by CMDS-II(Si) and even by more than an order of magnitude for heavy DM with mass larger than $80$ GeV, compared with the results using the empirical Helm form factor. Our results further indicate that the form factor effect can significantly modify the recoil spectrum of Xe-based detectors for heavy IVDM with $f_n/f_p=-0.7$.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.

Effects of the non-Newtonian gravity on properties of finite nuclei are studied by consistently incorporating both the direct and exchange contribution of the Yukawa potential in the Hartree-Fock approach using a well-tested Skyrme force for the strong interaction. It is shown for the first time that the strength of the Yukawa term in the non-Newtonian gravity is limited to $\log(|\alpha|)<1.75/[\lambda(\rm fm)]^{0.54} + 33.6$ within the length scale of $\lambda=1-10$ fm in order for the calculated properties of finite nuclei not to be in conflict with accurate experimental data available.

Waves play a crucial role in diagnosing the plasma properties of various structures in the solar corona and coronal heating. Slow magneto-acoustic (MA) waves are one of the important magnetohydrodynamic waves. In past decades, numerous slow MA waves were detected above the active regions and coronal holes, but rarely found elsewhere. Here, we investigate a `tornado'-like structure consisting of quasi-periodic streaks within a dark cavity at about 40--110 Mm above the quiet-Sun region on 2011 September 25. Our analysis reveals that these streaks are actually slow MA wave trains. The properties of these wave trains, including the phase speed, compression ratio, kinetic energy density, etc., are similar to those of the reported slow MA waves, except that the period of these waves is about 50 s, much shorter than the typical reported values (3--5 minutes).

In this talk, we report our recent study of constraining the non-Newtonian gravity at femtometer scale. We incorporate the Yukawa-type non-Newtonian gravitational potential consistently to the Skyrme functional form using the exact treatment for the direct contribution and density-matrix expansion method for the exchange contribution. The effects from the non-Newtonian potential on finite nuclei properties are then studied together with a well-tested Skyrme force. Assuming that the framework without non-Newtonian gravity can explain the binding energies and charge radii of medium to heavy nuclei within 2% error, we set an upper limit for the strength of the non-Newtonian gravitational potential at femtometer scale.

We explore effects of the light vector $U$-boson, which is weakly coupled to nucleons, on the transition density $\rho_{t}$ and pressure $P_{t}$ at the inner edge separating the liquid core from the solid crust of neutron stars. Three methods, i.e., the thermodynamical approach, the curvature matrix approach and the Vlasov equation approach are used to determine the transition density $\rho_{t}$ with the Skyrme effective nucleon-nucleon interactions. We find that the $\rho_{t}$ and $P_{t}$ depend on not only the ratio of coupling strength to mass squared of the $U$-boson $g^{2}/\mu ^{2}$ but also its mass $\mu $ due to the finite range interaction from the $U$-boson exchange. In particular, our results indicate that the $\rho_{t}$ and $P_{t}$ are sensitive to both $g^{2}/\mu ^{2}$ and $\mu $ if the $U$-boson mass $\mu $ is larger than about 2 MeV. Furthermore, we show that both $g^{2}/\mu ^{2}$ and $\mu $ can have significant influence on the mass-radius relation and the crustal fraction of total moment of inertia of neutron stars. In addition, we study the exchange term contribution of the $U$-boson based on the density matrix expansion method, and demonstrate that the exchange term effects on the nuclear matter equation of state as well as the $\rho_{t}$ and $P_{t}$ are generally negligible.

Several anomalies appear to be present in the large-angle cosmic microwave background (CMB) anisotropy maps of WMAP, including the alignment of large-scale multipoles. Models in which isotropy is spontaneously broken (e.g., by a scalar field) have been proposed as explanations for these anomalies, as have models in which a preferred direction is imposed during inflation. We examine models inspired by these, in which isotropy is broken by a multiplicative factor with dipole and/or quadrupole terms. We evaluate the evidence provided by the multipole alignment using a Bayesian framework, finding that the evidence in favor of the model is generally weak. We also compute approximate changes in estimated cosmological parameters in the broken-isotropy models. Only the overall normalization of the power spectrum is modified significantly.

In this paper, we present an automated system, which has the capability to catch and track solar limb prominences based on observations from EUV 304 passband. The characteristic parameters and their evolution, including height, position angle, area, length and brightness, are obtained without manual interventions. By applying the system to the STEREO-B/SECCHI/EUVI 304 data during 2007 April -2009 October, we obtain a total of 9477 well-tracked prominences and a catalog of these events available online at http://space.ustc.edu.cn/dreams/slipcat/. A detailed analysis of these prominences suggests that the system has a rather good performance. We have obtained several interesting statistical results based on the catalog. Most prominences appear below the latitude of 60 degrees and at the height of about 26 Mm above the solar surface. Most of them are quite stable during the period they are tracked. Nevertheless, some prominences have an upward speed of more than 100 km/s, and some others show significant downward and/or azimuthal speeds. There are strong correlations among the brightness, area and height. The expansion of a prominence is probably one major cause of its fading during the rising or erupting process.

The numerical studies of the interplanetary coupling between multiple magnetic clouds (MCs) are continued by a 2.5-dimensional ideal magnetohydrodynamic (MHD) model in the heliospheric meridional plane. The interplanetary direct collision (DC) / oblique collision (OC) between both MCs results from their same/different initial propagation orientations. Here the OC is explored in contrast to the results of the DC (Xiong et al., 2007). Both the slow MC1 and fast MC2 are consequently injected from the different heliospheric latitudes to form a compound stream during the interplanetary propagation. The MC1 and MC2 undergo contrary deflections during the process of oblique collision. Their deflection angles of $|\delta \theta_1|$ and $|\delta \theta_2|$ continuously increase until both MC-driven shock fronts are merged into a stronger compound one. The $|\delta \theta_1|$, $|\delta \theta_2|$, and total deflection angle $\Delta \theta$ ($\Delta \theta = |\delta \theta_1| + |\delta \theta_2|$) reach their corresponding maxima when the initial eruptions of both MCs are at an appropriate angular difference. Moreover, with the increase of MC2's initial speed, the OC becomes more intense, and the enhancement of $\delta \theta_1$ is much more sensitive to $\delta \theta_2$. The $|\delta\theta_1|$ is generally far less than the $|\delta\theta_2|$, and the unusual case of $|\delta\theta_1|\simeq|\delta\theta_2|$ only occurs for an extremely violent OC. But because of the elasticity of the MC body to buffer the collision, this deflection would gradually approach an asymptotic degree. Therefore, the deflection due to the OC should be considered for the evolution and ensuing geoeffectiveness of interplanetary interaction among successive coronal mass ejections (CMEs).

Numerical studies of the interplanetary "multiple magnetic clouds (Multi-MC)" are performed by a 2.5-dimensional ideal magnetohydrodynamic (MHD) model in the heliospheric meridional plane. Both slow MC1 and fast MC2 are initially emerged along the heliospheric equator, one after another with different time interval. The coupling of two MCs could be considered as the comprehensive interaction between two systems, each comprising of an MC body and its driven shock. The MC2-driven shock and MC2 body are successively involved into interaction with MC1 body. The momentum is transferred from MC2 to MC1. After the passage of MC2-driven shock front, magnetic field lines in MC1 medium previously compressed by MC2-driven shock are prevented from being restored by the MC2 body pushing. MC1 body undergoes the most violent compression from the ambient solar wind ahead, continuous penetration of MC2-driven shock through MC1 body, and persistent pushing of MC2 body at MC1 tail boundary. As the evolution proceeds, the MC1 body suffers from larger and larger compression, and its original vulnerable magnetic elasticity becomes stiffer and stiffer. So there exists a maximum compressibility of Multi-MC when the accumulated elasticity can balance the external compression. With respect to Multi-MC geoeffectiveness, the evolution stage is a dominant factor, whereas the collision intensity is a subordinate one. The magnetic elasticity, magnetic helicity of each MC, and compression between each other are the key physical factors for the formation, propagation, evolution, and resulting geoeffectiveness of interplanetary Multi-MC.

Numerical studies have been performed to interpret the observed "shock overtaking magnetic cloud (MC)" event by a 2.5 dimensional magnetohydrodynamic (MHD) model in heliospheric meridional plane. Results of an individual MC simulation show that the MC travels with a constant bulk flow speed. The MC is injected with very strong inherent magnetic field over that in the ambient flow and expands rapidly in size initially. Consequently, the diameter of MC increases in an asymptotic speed while its angular width contracts gradually. Meanwhile, simulations of MC-shock interaction are also presented, in which both a typical MC and a strong fast shock emerge from the inner boundary and propagate along heliospheric equator, separated by an appropriate interval. The results show that the shock firstly catches up with the preceding MC, then penetrates through the MC, and finally merges with the MC-driven shock into a stronger compound shock. The morphologies of shock front in interplanetary space and MC body behave as a central concave and a smooth arc respectively. The compression and rotation of magnetic field serve as an efficient mechanism to cause a large geomagnetic storm. The MC is highly compressed by the the overtaking shock. Contrarily, the transport time of incidental shock influenced by the MC depends on the interval between their commencements. Maximum geoeffectiveness results from that when the shock enters the core of preceding MC, which is also substantiated to some extent by a corresponding simplified analytic model. Quantified by $Dst$ index, the specific result gives that the geoeffectiveness of an individual MC is largely enhanced with 80% increment in maximum by an incidental shock.

Numerical studies of the interplanetary "shock overtaking magnetic cloud (MC)" event are continued by a 2.5 dimensional magnetohydrodynamic (MHD) model in heliospheric meridional plane. Interplanetary direct collision (DC)/oblique collision (OC) between an MC and a shock results from their same/different initial propagation orientations. For radially erupted MC and shock in solar corona, the orientations are only determined respectively by their heliographic locations. OC is investigated in contrast with the results in DC \citep{Xiong2006}. The shock front behaves as a smooth arc. The cannibalized part of MC is highly compressed by the shock front along its normal. As the shock propagates gradually into the preceding MC body, the most violent interaction is transferred sideways with an accompanying significant narrowing of the MC's angular width. The opposite deflections of MC body and shock aphelion in OC occur simultaneously through the process of the shock penetrating the MC. After the shock's passage, the MC is restored to its oblate morphology. With the decrease of MC-shock commencement interval, the shock front at 1 AU traverses MC body and is responsible for the same change trend of the latitude of the greatest geoeffectiveness of MC-shock compound. Regardless of shock orientation, shock penetration location regarding the maximum geoeffectiveness is right at MC core on the condition of very strong shock intensity. An appropriate angular difference between the initial eruption of an MC and an overtaking shock leads to the maximum deflection of the MC body. The larger the shock intensity is, the greater is the deflection angle. The interaction of MCs with other disturbances could be a cause of deflected propagation of interplanetary coronal mass ejection (ICME).

With an exponentially increasing amount of astronomical data, the complexity and dimension of astronomical data are likewise growing rapidly. Extracting information from such data becomes a critical and challenging problem. For example, some algorithms can only be employed in the low-dimensional spaces, so feature selection and feature extraction become important topics. Here we describe the difference between feature selection and feature extraction methods, and introduce the taxonomy of feature selection methods as well as the characteristics of each method. We present a case study comparing the performance and computational cost of different feature selection methods. For the filter method, ReliefF and fisher filter are adopted; for the wrapper method, improved CHAID, linear discriminant analysis (LDA), Naive Bayes (NB) and C4.5 are taken as learners. Applied on the sample, the result indicates that from the viewpoints of computational cost the filter method is superior to the wrapper method. Moreover, different learning algorithms combined with appropriate feature selection methods may arrive at better performance.