The impact of the streaming between baryons and dark matter on the first structures has been actively explored by recent studies. We investigate how the key results are affected by two popular approximations. One is to implement the streaming by accounting for only the relative motion while assuming "baryons trace dark matter" spatially at the initialization of simulation. This neglects the smoothing on the gas density taking place before the initialization. In our simulation initialized at $z_i=200$, it overestimates the gas density power spectrum up to 40\% at $k\approx10^2~h~\mbox{Mpc}^{-1}$ at $z=20$. Halo mass ($M_h$) and baryonic fraction in halos ($f_{b,h}$) are also overestimated, but the relation between the two remains unchanged. The other approximation tested is to artificially amplify the density/velocity fluctuations in the cosmic mean density to simulate the first minihalos that form in overdense regions. This gives a head start to the halo growth while the subsequent growth is similar to that in the mean density. The growth in a true overdense region, on the other hand, is accelerated gradually in time. For example, raising $\sigma_8$ by 50\% effectively transforms $z\rightarrow\sqrt{1.5}z$ in the halo mass growth history while in 2-$\sigma$ overdensity, the growth is accelerated by a constant in redshift: $z\rightarrow{z+4.8}$. As a result, halos are more grown in the former than in the latter before $z\approx27$ and vice versa after. The $f_{b,h}$-$M_h$ relation is unchanged in those cases as well, suggesting that the Pop III formation rate for a given $M_h$ is insensitive to the tested approximations.

Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

Spicules are rapidly evolving fine-scale jets of magnetized plasma in the solar chromosphere. It remains unclear how these prevalent jets originate from the solar surface and what role they play in heating the solar atmosphere. Using the Goode Solar Telescope at the Big Bear Solar Observatory, we observed spicules emerging within minutes of the appearance of opposite-polarity magnetic flux around dominant-polarity magnetic field concentrations. Data from the Solar Dynamics Observatory showed subsequent heating of the adjacent corona. The dynamic interaction of magnetic fields (likely due to magnetic reconnection) in the partially ionized lower solar atmosphere appears to generate these spicules and heat the upper solar atmosphere.

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.

Chromospheric activities prior to solar flares provide important clues to solar flare initiation, but are as yet poorly understood. We report a significant and rapid H$\alpha$ line broadening before the solar flare SOL2011-09-29T18:08, that was detected using the unprecedented high-resolution H$\alpha$ imaging spectroscopy with the Fast Imaging Solar Spectrograph (FISS) installed on the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory.The strong H$\alpha$ broadening extends as a blue excursion up to -4.5 A and as a red excursion up to 2.0 A, which implies a mixture of velocities in the range of -130 km/s to 38 km/s by applying the cloud model, comparable to the highest chromospheric motions reported before. The H$\alpha$ blueshifted broadening lasts for about 4 minutes, and is temporally and spatially correlated with the start of a rising filament, which is later associated with the main phase of the flare as detected by the Atmosphere Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The potential importance of this H$\alpha$ blueshifted broadening as a preflare chromospheric activity is briefly discussed within the context of the two-step eruption model.

The decline in abundance of Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$\alpha$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$\alpha$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$\alpha$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 5\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$\alpha$ peak. We discuss our results in the context of the interpretation of observations.

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.

Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $\Lambda$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other.

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

A very large dynamic range with simultaneous capture of both large- and small-scales in the simulations of cosmic structures is required for correct modelling of many cosmological phenomena, particularly at high redshift. This is not always available, or when it is, it makes such simulations very expensive. We present a novel sub-grid method for modelling low-mass ($10^5\,M_\odot\leq M_{\rm halo}\leq 10^9\,M_\odot$) haloes, which are otherwise unresolved in large-volume cosmological simulations limited in numerical resolution. In addition to the deterministic halo bias that captures the average property, we model its stochasticity that is correlated in time. We find that the instantaneous binned distribution of the number of haloes is well approximated by a log-normal distribution, with overall amplitude modulated by this "temporal correlation bias". The robustness of our new scheme is tested against various statistical measures, and we find that temporally correlated stochasticity generates mock halo data that is significantly more reliable than that from temporally uncorrelated stochasticity. Our method can be applied for simulating processes that depend on both the small- and large-scale structures, especially for those that are sensitive to the evolution history of structure formation such as the process of cosmic reionization. As a sample application, we generate a mock distribution of medium-mass ($ 10^{8} \leq M/M_{\odot} \leq 10^{9}$) haloes inside a 500 Mpc$\,h^{-1}$, $300^3$ grid simulation box. This mock halo catalogue bears a reasonable statistical agreement with a halo catalogue from numerically-resolved haloes in a smaller box, and therefore will allow a very self-consistent sets of cosmic reionization simulations in a box large enough to generate statistically reliable data.

The Near-IR Imaging Spectropolarimeter (NIRIS) is a polarimeter that is installed at the New Solar Telescope at Big Bear Solar Observatory. This instrument takes advantages of the highest spatial resolution and flux. The primary mirror is an on-axis type, so it was our interest to evaluate its contribution to the crosstalk among the Stokes parameters since we could not put our calibration optics before the mirror. We would like to present our efforts to compensate for the crosstalk among Stokes profiles caused by the relay optics from the telescope to the detector. The overall data processing pipeline is also introduced.

In this study we report detailed observations of magnetic environment at four footpoints of two warm coronal loops observed on 5 May 2016 in NOAA AR 12542 (Loop I) and 17 Dec 2015 in NOAA AR 12470 (Loop II). These loops were connecting a plage region with sunspot periphery (Loop I) and a sunspot umbra (Loop II). We used Solar Dynamics Observatory (SDO) and Goode Solar Telescope (GST) data to describe the phenomenon and understand its causes. The study indicates loop brightening episodes were associated with magnetic flux emergence and cancellation processes observed in SDO's Helioseismic and Magnetic Imager (HMI) and GST's Near InfraRed Imaging Spectrapolarimeter (NIRIS) data. The observed activity was driven by magnetic reconnection between small-scale emerging dipoles and large-scale pre-existing fields, suggesting that the reconnection occurred in the lower chromosphere at the edge of an extended plage region, where the loops were rooted. We suggest that plasma, evaporated during these reconnection events, gradually filled the loops and as it cooled the visible density front propagated from one footpoint of the loop to another at a rate of 90-110 km s$^{-1}$. This study also indicates that at least some of the bright loops seen in SDO Atmospheric Imaging Assembly images rooted in sunspot umbra may be heated due to magnetic activity taking place at the remote (non-sunspot) footpoint.

We present a systematic study of the cosmic variance that existed in the formation of first stars and galaxies. We focus on the cosmic variance induced by the large-scale density and velocity environment engraved at the epoch of recombination. The density environment is predominantly determined by the dark-matter overdensity, and the velocity environment by the dark matter-baryon streaming velocity. Toward this end, we introduce a new cosmological initial condition generator BCCOMICS, which solves the quasi-linear evolution of small-scale perturbations under the large-scale density and streaming-velocity environment and generates the initial condition for dark matter and baryons, as either particles or grid data at a specific redshift. We also describe a scheme to simulate the formation of first galaxies inside density peaks and voids, where a local environment is treated as a separate universe. The resulting cosmic variance in the minihalo number density and the amount of cooling mass are presented as an application. Density peaks become a site for enhanced formation of first galaxies, which compete with the negative effect from the dark matter-baryon streaming velocity on structure formation.

The photospheric response to solar flares, also known as coronal back reaction, is often observed as sudden flare-induced changes in vector magnetic field and sunspot motions. However, it remains obscure whether evolving flare ribbons, the flare signature closest to the photosphere, are accompanied by changes in vector magnetic field therein. Here we explore the relationship between the dynamics of flare ribbons in the chromosphere and variations of magnetic fields in the underlying photosphere, using high-resolution off-band H-alpha images and near-infrared vector magnetograms of the M6.5 flare on 2015 June 22 observed with the 1.6 m Goode Solar Telescope. We find that changes of photospheric fields occur at the arrival of the flare ribbon front, thus propagating analogously to flare ribbons. In general, the horizontal field increases and the field lines become more inclined to the surface. When ribbons sweep through regions that undergo a rotational motion, the fields transiently turn more vertical with decreased horizontal field and inclination angle, and then restore and/or become more horizontal than before the ribbon arrival. The ribbon propagation decelerates near the sunspot rotation center, where the vertical field becomes permanently enhanced. Similar magnetic field changes are discernible in magnetograms from the Helioseismic and Magnetic Imager (HMI), and an inward collapse of coronal magnetic fields is inferred from the time sequence of non-linear force-free field models extrapolated from HMI magnetograms. We conclude that photospheric fields respond nearly instantaneously to magnetic reconnection in the corona.

We present unprecedented high-resolution TiO images and Fe I 1565 nm spectropolarimetric data of two light bridges taken by the 1.6-m Goode Solar Telescope at Big Bear Solar Observatory. In the first light bridge (LB1), we find striking knot-like dark structures within the central dark lane. Many dark knots show migration away from the penumbra along the light bridge. The sizes, intensity depressions and apparent speeds of their proper motion along the light bridges of 33 dark knots identified from the TiO images are mainly in the ranges of 80$\sim$200~km, 30\%$\sim$50\%, and 0.3$\sim$1.2~km~s$^{-1}$, respectively. In the second light bridge (LB2), a faint central dark lane and striking transverse intergranular lanes were observed. These intergranular lanes have sizes and intensity depressions comparable to those of the dark knots in LB1, and also migrate away from the penumbra at similar speeds. Our observations reveal that LB2 is made up of a chain of evolving convection cells, as indicated by patches of blue shift surrounded by narrow lanes of red shift. The central dark lane generally corresponds to blueshifts, supporting the previous suggestion of central dark lanes being the top parts of convection upflows. In contrast, the intergranular lanes are associated with redshifts and located at two sides of each convection cell. The magnetic fields are stronger in intergranular lanes than in the central dark lane. These results suggest that these intergranular lanes are manifestations of convergent convective downflows in the light bridge. We also provide evidence that the dark knots observed in LB1 may have a similar origin.

We study five sequential solar flares (SOL2015-08-07) occurring in Active Region 12396 observed with the Goode Solar Telescope (GST) at the BBSO, complemented by IRIS and SDO observations. The main flaring region is an arch filament system (AFS) consisting of multiple bundles of dark filament threads enclosed by semi-circular flare ribbons. We study the magnetic configuration and evolution of the active region by constructing coronal magnetic field models based on SDO/HMI magnetograms using two independent methods, i.e., the nonlinear force-free field (NLFFF) extrapolation and the flux rope insertion method. The models consist of multiple flux ropes with mixed signs of helicity, i.e., positive (negative) in the northern (southern) region, which is consistent with the GST observations of multiple filament bundles. The footprints of quasi-separatrix layers (QSLs) derived from the extrapolated NLFFF compare favorably with the observed flare ribbons. An interesting double-ribbon fine structure located at the east border of the AFS is consistent with the fine structure of the QSL's footprint. Moreover, magnetic field lines traced along the semi-circular footprint of a dome-like QSL surrounding the AFS are connected to the regions of significant helicity and Poynting flux injection. The maps of magnetic twist show that positive twist became dominant as time progressed, which is consistent with the injection of positive helicity before the flares. We hence conclude that these circular shaped flares are caused by 3D magnetic reconnection at the QSLs associated with the AFS possessing mixed signs of helicity.

As one of the most violent eruptions on the Sun, flares are believed to be powered by magnetic reconnection. The fundamental physics involving the release, transfer and deposition of energy have been studied extensively. Taking advantage of the unprecedented resolution provided by the 1.6-m Goode Solar Telescope, here we show a sudden rotation of vector magnetic fields, about 12$^{\circ}$-20$^{\circ}$ counterclockwise, associated with a flare. Unlike the permanent changes reported previously, the azimuth-angle change is transient and co-spatial/temporal with H$\alpha$ emission. The measured azimuth angle becomes closer to that in potential fields suggesting untwist of flare loops. The magnetograms were obtained in the near infrared at 1.56~$\mu$m, which is minimally affected by flare emission and no intensity profile change was detected. We believe that these transient changes are real and discuss the possible explanations in which the high energy electron beams or $Alfv\acute{e}n$ waves play a crucial role.

Solar Active Region (AR) 12673 is the most flare productive AR in the solar cycle 24. It produced four X-class flares including the X9.3 flare on 06 September 2017 and the X8.2 limb event on 10 September. Sun and Norton (2017) reported that this region had an unusual high rate of flux emergence, while Huang et al. (2018) reported that the X9.3 flare had extremely strong white-light flare emissions. Yang at al. (2017) described the detailed morphological evolution of this AR. In this report, we focus on usual behaviors of the light bridge (LB) dividing the delta configuration of this AR, namely the strong magnetic fields (above 5500 G) in the LB and apparent photospheric twist as shown in observations with a 0.1 arcsec spatial resolution obtained by the 1.6m telescope at Big Bear Solar Observatory.

We study the impact of thermal inflation on the formation of cosmological structures and present astrophysical observables which can be used to constrain and possibly probe the thermal inflation scenario. These are dark matter halo abundance at high redshifts, satellite galaxy abundance in the Milky Way, and fluctuation in the 21-cm radiation background before the epoch of reionization. The thermal inflation scenario leaves a characteristic signature on the matter power spectrum by boosting the amplitude at a specific wavenumber determined by the number of e-foldings during thermal inflation ($N_{\rm bc}$), and strongly suppressing the amplitude for modes at smaller scales. For a reasonable range of parameter space, one of the consequences is the suppression of minihalo formation at high redshifts and that of satellite galaxies in the Milky Way. While this effect is substantial, it is degenerate with other cosmological or astrophysical effects. The power spectrum of the 21-cm background probes this impact more directly, and its observation may be the best way to constrain the thermal inflation scenario due to the characteristic signature in the power spectrum. The Square Kilometre Array (SKA) in phase 1 (SKA1) has sensitivity large enough to achieve this goal for models with $N_{\rm bc}\gtrsim 26$ if a 10000-hr observation is performed. The final phase SKA, with anticipated sensitivity about an order of magnitude higher, seems more promising and will cover a wider parameter space.

Solar flares are generally believed to be powered by free magnetic energy stored in the corona, but the build up of coronal energy alone may be insufficient for the imminent flare occurrence. The flare onset mechanism is a critical but less understood problem, insights into which could be gained from small-scale energy releases known as precursors, which are observed as small pre-flare brightenings in various wavelengths, and also from certain small-scale magnetic configurations such as the opposite polarity fluxes, where magnetic orientation of small bipoles is opposite to that of the ambient main polarities. However, high-resolution observations of flare precursors together with the associated photospheric magnetic field dynamics are lacking. Here we study precursors of a flare using unprecedented spatiotemporal resolution of the 1.6 m New Solar Telescope, complemented by novel microwave data. Two episodes of precursor brightenings are initiated at a small-scale magnetic channel (a form of opposite polarity fluxes) with multiple polarity inversions and enhanced magnetic fluxes and currents, lying near the footpoints of sheared magnetic loops. The low-atmospheric origin of these precursor emissions is corroborated by microwave spectra. We propose that the emerging magnetic channel field interacts with the sheared arcades to cause precursor brightenings at the main flare core region. These high-resolution results provide evidence of low-atmospheric small-scale energy release and possible relationship to the onset of the main flare.

Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\bar\delta$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\bar\delta]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization.

Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.

Cosmic reionization is thought to be primarily fueled by the first generations of galaxies. We examine their stellar and gaseous properties, focusing on the star formation rates and the escape of ionizing photons, as a function of halo mass, redshift, and environment using the full suite of the {\it Renaissance Simulations} with an eye to provide better inputs to global reionization simulations. This suite, carried out with the adaptive mesh refinement code Enzo, is unprecedented in terms of their size and physical ingredients. The simulations probe overdense, average, and underdense regions of the universe of several hundred comoving Mpc$^3$, each yielding a sample of over 3,000 halos in the mass range $10^7 - 10^{9.5}~\Ms$ at their final redshifts of 15, 12.5, and 8, respectively. In the process, we simulate the effects of radiative and supernova feedback from 5,000 to 10,000 metal-free (Population III) stars in each simulation. We find that halos as small as $10^7~\Ms$ are able to form stars due to metal-line cooling from earlier enrichment by massive Population III stars. However, we find such halos do not form stars continuously. Using our large sample, we find that the galaxy-halo occupation fraction drops from unity at virial masses above $10^{8.5}~\Ms$ to $\sim$50\% at $10^8 ~\Ms$ and $\sim$10\% at $10^7~\Ms$, quite independent of redshift and region. Their average ionizing escape fraction is $\sim$5\% in the mass range $10^8 - 10^9~\Ms$ and increases with decreasing halo mass below this range, reaching 40--60\% at $10^7~\Ms$. Interestingly, we find that the escape fraction varies between 10--20\% in halos with virial masses $\sim 3 \times 10^9~\Ms$. Taken together, our results confirm the importance of the smallest galaxies as sources of ionizing radiation contributing to the reionization of the universe.

We study the dynamical effect of relative velocities between dark matter and baryonic fluids, which remained supersonic after the epoch of recombination. The impact of this supersonic motion on the formation of cosmological structures was first formulated by Tseliakhovich & Hirata (2010), in terms of the linear theory of small-scale fluctuations coupled to large-scale, relative velocities in mean-density regions. In their formalism, they limited the large-scale density environment to be those of the global mean density. We improve on their formulation by allowing variation in the density environment as well as the relative velocities. This leads to a new type of coupling between large-scale and small-scale modes. We find that the small-scale fluctuation grows in a biased way: faster in the overdense environment and slower in the underdense environment. We also find that the net effect on the global power spectrum of the density fluctuation is to boost its overall amplitude from the prediction by Tseliakhovich & Hirata (2010). Correspondingly, the conditional mass function of cosmological halos and the halo bias parameter are both affected in a similar way. The discrepancy between our prediction and that by Tseliakhovich & Hirata (2010) is significant, and therefore the related cosmology and high-redshift astrophysics should be revisited. The mathematical formalism of this study can be used for generating cosmological initial conditions of small-scale perturbations in generic, overdense (underdense) background patches.

Observations of the epoch of reionization give us clues about the nature and evolution of the sources of ionizing photons, or early stars and galaxies. We present a new suite of structure formation and radiative transfer simulations from the PRACE4LOFAR project designed to investigate whether the mechanism of radiative feedback, or the suppression of star formation in ionized regions from UV radiation, can be inferred from these observations. Our source halo mass extends down to $10^8 M_\odot$, with sources in the mass range $10^8$ to $10^9 M_\odot$ expected to be particularly susceptible to feedback from ionizing radiation, and we vary the aggressiveness and nature of this suppression. Not only do we have four distinct source models, we also include two box sizes (67 Mpc and 349 Mpc), each with two grid resolutions. This suite of simulations allows us to investigate the robustness of our results. All of our simulations are broadly consistent with the observed electron-scattering optical depth of the cosmic microwave background and the neutral fraction and photoionization rate of hydrogen at $z\sim6$. In particular, we investigate the redshifted 21-cm emission in anticipation of upcoming radio interferometer observations. We find that the overall shape of the 21-cm signal and various statistics are robust to the exact nature of source suppression, the box size, and the resolution. There are some promising model discriminators in the non-Gaussianity and small-scale power spectrum of the 21-cm signal.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

We present estimates of the nonlinear bias of cosmological halo formation, spanning a wide range in the halo mass from $\sim 10^{5} M_\odot$ to $\sim 10^{12} M_\odot$, based upon both a suite of high-resolution cosmological N-body simulations and theoretical predictions. The halo bias is expressed in terms of the mean bias and stochasticity as a function of local overdensity ($\delta$), under different filtering scales, which is realized as the density of individual cells in uniform grids. The sampled overdensities span a range wide enough to provide the fully nonlinear bias effect on the formation of haloes. A strong correlation between $\delta$ and halo population overdensity $\delta_h$ is found, along with sizable stochasticity. We find that the empirical mean halo bias matches, with good accuracy, the prediction by the peak-background split method based on the excursion set formalism, as long as the empirical, globally-averaged halo mass function is used. Consequently, this bias formalism is insensitive to uncertainties caused by varying halo identification schemes, and can be applied generically. We also find that the probability distribution function of biased halo numbers has wider distribution than the pure Poisson shot noise, which is attributed to the sub-cell scale halo correlation. We explicitly calculate this correlation function and show that both overdense and underdense regions have positive correlation, leading to stochasticity larger than the Poisson shot noise in the range of haloes and halo-collapse epochs we study.

We present a novel method to implement time-delayed propagation of radiation fields in cosmological radiative transfer simulations. Time-delayed propagation of radiation fields requires construction of retarded-time fields by tracking the location and lifetime of radiation sources along the corresponding light-cones. Cosmological radiative transfer simulations have, until now, ignored this "light-cone effect" or implemented ray-tracing methods that are computationally demanding. We show that radiative transfer calculation of the time-delayed fields can be easily achieved in numerical simulations when periodic boundary conditions are used, by calculating the time-discretized retarded-time Green's function using the Fast Fourier Transform (FFT) method and convolving it with the source distribution. We also present a direct application of this method to the long-range radiation field of Lyman-Werner band photons, which is important in the high-redshift astrophysics with first stars.

We present our prediction for the local 21 cm differential brightness temperature ($\delta T_{b}$) from a set of strongly clustered sources of Population III (Pop III) and II (Pop II) objects in the early Universe, by a numerical simulation of their formation and radiative feedback. These objects are located inside a highly biased environment, which is a rare, high-density peak ("Rarepeak") extending to $\sim7$ comoving Mpc. We study the impact of ultraviolet (UV) and X-ray photons on the intergalactic medium (IGM) and the resulting $\delta T_{b}$, when Pop III stars are assumed to emit X-ray photons by forming X-ray binaries very efficiently. We parameterize the rest-frame spectral energy distribution (SED) of X-ray photons, which regulates X-ray photon-trapping, IGM-heating, secondary Lyman-alpha pumping and the resulting morphology of $\delta T_{b}$. A combination of emission ($\delta T_{b}>0$) and absorption ($\delta T_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($\delta\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($\theta\sim10'$) object for 21 cm observation at $13\lesssim z\lesssim17$, which is found to be detectable as a single object by SKA with integration time of $\sim1000$ hours. Power spectrum analysis by some of the SKA precursors (LOFAR, MWA, PAPER) of such rare peaks is found difficult due to the rarity of these peaks, and the contribution only by these rare peaks to the total power spectrum remains subdominant compared to that by all astrophysical sources.

We present an overview of the theory of high-redshift star and X-ray source formation, and how they affect the 21-cm background. Primary focus is given to Lyman alpha pumping and X-ray heating mechanisms at cosmic dawn, opening a new observational window for high-redshift astrophysics by generating sizable fluctuations in the 21-cm background. We describe observational prospects for power spectrum analysis and 3D tomography (imaging) of the signature of these early astrophysical sources by SKA1-LOW and SKA2.

Due to their long mean free path, X-rays are expected to have an important impact on cosmic reionization by heating and ionizing the IGM on large scales, especially after simulations have suggested that Population III stars may form in pairs at redshifts as high as 20 - 30. We use the Pop III distribution and evolution from a self-consistent cosmological radiation hydrodynamics simulation of the formation of the first galaxies and a simple Pop III X-ray binaries model to estimate their X-ray output in a high density region larger than 100 comoving (Mpc)$^3$. We then combine three different methods --- ray tracing, a one-zone model, and X-ray background modeling --- to investigate the X-ray propagation, intensity distribution, and long term effects on the IGM thermal and ionization state. The efficiency and morphology of photo-heating and photo-ionization are dependent on the photon energies. The sub-keV X-rays only impact the IGM near the sources while the keV photons contribute significantly to the X-ray background and heat and ionize the IGM smoothly. The X-rays just below 1 keV are most efficient in heating and ionizing the IGM. We find that the IGM might be heated to over 100 K by $z=10$ and the high density source region might reach 10$^4$ K, limited by atomic hydrogen cooling. This may be important for predicting the 21-cm neutral hydrogen signals. But, on the other hand, the free electrons from X-ray ionizations are not enough to contribute significantly to the optical depth of CMB to the Thomson scattering.

A novel method to characterize the topology of the early-universe intergalactic medium during the epoch of cosmic reionization is presented. The 21-cm radiation background from high redshift is analyzed through the calculation of the 2-dimensional (2D) genus. The radiative transfer of hydrogen-ionizing photons and ionization-rate equations are calculated in a suite of numerical simulations under various input parameters. The 2D genus is calculated from the mock 21-cm images of the high-redshift Universe. We construct the 2D genus curve by varying the threshold differential brightness temperature, and compare this to the 2D genus curve of the underlying density field. We find that (1) the 2D genus curve reflects the evolutionary track of cosmic reionization and (2) the 2D genus curve can discriminate between certain reionization scenarios and thus indirectly probe the properties of radiation-sources. Choosing the right beam shape of a radio antenna is crucial for this analysis. To this end, the Square Kilometre Array (SKA) is found to be a suitable apparatus for this analysis in terms of sensitivity, even though some deterioration of the data for this purpose is unavoidable under the planned size of the antenna core.

Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has the largest impact at two different stages of reionization: one when reionization is $\sim 20\%$ and other when it is $\sim 80\%$ completed. We find a factor of $\sim 4$ amplification of the power spectrum at the largest scale available in our simulations. We do not find any significant anisotropy in the 21-cm power spectrum due to the light cone effect. We argue that for the power spectrum to become anisotropic, the light cone effect would have to make the ionized bubbles significantly elongated or compressed along the line-of-sight, which would require extreme reionization scenarios. We also calculate the two-point correlation functions parallel and perpendicular to the line-of-sight and find them to differ. Finally, we calculate an optimum frequency bandwidth below which the light cone effect can be neglected when extracting power spectra from observations. We find that if one is willing to accept a $10 \%$ error due to the light cone effect, the optimum frequency bandwidth for $k= 0.056 \, \rm{Mpc}^{-1}$ is $\sim 7.5$ MHz. For $k = 0.15$ and $0.41 \, \rm{Mpc}^{-1}$ the optimum bandwidth is $\sim 11$ and $\sim 16$ MHz respectively.

We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $\Delta z$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged)

An essential parameter for models of coronal heating and fast solar wind acceleration that rely on the dissipation of MHD turbulence is the characteristic energy-containing length $\lambda_{\bot}$ of the squared velocity and magnetic field fluctuations ($u^2$ and $b^2$) transverse to the mean magnetic field inside a coronal hole (CH) at the base of the corona. The characteristic length scale defines directly the heating rate. We use a time series analysis of solar granulation and magnetic field measurements inside two CHs obtained with the New Solar Telescope (NST) at Big Bear Solar Observatory. A data set for transverse magnetic fields obtained with the Solar Optical Telescope/Spectro-Polarimeter (SOT/SP) aboard {\it Hinode} spacecraft was utilized to analyze the squared transverse magnetic field fluctuations $b_t^2$. Local correlation tracking (LCT) was applied to derive the squared transverse velocity fluctuations $u^2$. We find that for $u^2$-structures, Batchelor integral scale $\lambda$ varies in a range of 1800 - 2100 km, whereas the correlation length $\varsigma$ and the $e$-folding length $L$ vary between 660 and 1460 km. Structures for $b_t^2$ yield $\lambda \approx 1600$ km, $\varsigma \approx 640$ km, and $L \approx 620$ km. An averaged (over $\lambda, \varsigma$, and $L$) value of the characteristic length of $u^2$-fluctuations is 1260$\pm$500 km, and that of $b_t^2$ is 950$\pm$560 km. The characteristic length scale in the photosphere is approximately 1.5-50 times smaller than that adopted in previous models (3-30$\times10^3$ km). Our results provide a critical input parameter for current models of coronal heating and should yield an improved understanding of fast solar wind acceleration.

The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized).

We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization.

The Swift Gamma-ray Burst (GRB) observatory responds to GRB triggers with optical observations in ~ 100 s, but cannot respond faster than ~ 60 s. While some ground-based telescopes respond quickly, the number of sub-60 s detections remains small. In mid- to late-2013, the Ultra-Fast Flash Observatory-Pathfinder is to be launched on the Lomonosov spacecraft to investigate early optical GRB emission. This pathfinder mission is necessarily limited in sensitivity and event rate; here we discuss a next generation rapid-response space observatory. We list science topics motivating our instruments, those that require rapid optical-IR GRB response, including: A survey of GRB rise shapes/times, measurements of optical bulk Lorentz factors, investigation of magnetic dominated (vs. non-magnetic) jet models, internal vs. external shock origin of prompt optical emission, the use of GRBs for cosmology, and dust evaporation in the GRB environment. We also address the impacts of the characteristics of GRB observing on our instrument and observatory design. We describe our instrument designs and choices for a next generation observatory as a second instrument on a low-earth orbit spacecraft, with a 120 kg instrument mass budget. Restricted to relatively modest mass and power, we find that a coded mask X-ray camera with 1024 cm2 of detector area could rapidly locate about 64 GRB triggers/year. Responding to the locations from the X-ray camera, a 30 cm aperture telescope with a beam-steering system for rapid (~ 1 s) response and a near-IR camera should detect ~ 29 GRB, given Swift GRB properties. Am additional optical camera would give a broadband optical-IR slope, allowing dynamic measurement of dust extinction at the source, for the first time.

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes.

We used photospheric intensity images and magnetic field measurements from the New Solar Telescope in Big Bear and Helioseismic Magnetic Imager on board Solar Dynamics Observatory (SDO) to study the the effect that the new small-scale emerging flux induces on solar granulation. We report that emerging flux appears to leave different types of footprint on solar granulation: i) diffuse irregular patches of increased brightness, ii) well defined filament-like structures and accompanied bright points, and iii) bright point-like features that appear inside granules. We suggest that the type of the footprint depends on the intensity of emerging fields. Stronger fields, emerging as a part of large magnetic structure, create on the solar surface a well defined filamentary pattern with bright points at the ends of the filaments, while weak turbulent fields are associated with bright patches inside the host granule.

Observations of redshifted 21-cm radiation from neutral hydrogen during the epoch of reionization (EoR) are considered to constitute the most promising tool to probe that epoch. One of the major goals of the first generation of low frequency radio telescopes is to measure the 3D 21-cm power spectrum. However, the 21-cm signal could evolve substantially along the line of sight (LOS) direction of an observed 3D volume, since the received signal from different planes transverse to the LOS originated from different look-back times and could therefore be statistically different. Using numerical simulations we investigate this so-called light cone effect on the spherically averaged 3D 21-cm power spectrum. For this version of the power spectrum, we find that the effect mostly `averages out' and observe a smaller change in the power spectrum compared to the amount of evolution in the mean 21-cm signal and its rms variations along the LOS direction. Nevertheless, changes up to 50% at large scales are possible. In general the power is enhanced/suppressed at large/small scales when the effect is included. The cross-over mode below/above which the power is enhanced/suppressed moves toward larger scales as reionization proceeds. When considering the 3D power spectrum we find it to be anisotropic at the late stages of reionization and on large scales. The effect is dominated by the evolution of the ionized fraction of hydrogen during reionization and including peculiar velocities hardly changes these conclusions. We present simple analytical models which explain qualitatively all the features we see in the simulations.

The peculiar velocity of the intergalactic gas responsible for the cosmic 21cm background from the epoch of reionization and beyond introduces an anisotropy in the three-dimensional power spectrum of brightness temperature fluctuations. Measurement of this anisotropy by future 21cm surveys is a promising tool for separating cosmology from 21cm astrophysics. However, previous attempts to model the signal have often neglected peculiar velocity or only approximated it crudely. This paper re-examines the effects of peculiar velocity on the 21cm signal in detail, improving upon past treatment and addressing several issues for the first time. (1) We show that properly accounting for finite optical depth eliminates the unphysical divergence of 21cm brightness temperature in overdense regions of the IGM found by previous work that employed the usual optically-thin approximation. (2) The approximation made previously to circumvent the diverging brightness temperature problem by capping velocity gradient can misestimate the power spectrum on all scales. (3) The observed power spectrum in redshift-space remains finite even in the optically-thin approximation if one properly accounts for the redshift-space distortion. However, results that take full account of finite optical depth show that this approximation is only accurate in the limit of high spin temperature. (4) The linear theory for redshift-space distortion results in ~30% error in the observationally relevant wavenumber range, at the 50% ionized epoch. (5) We describe and test two numerical schemes to calculate the 21cm signal from reionization simulations to incorporate peculiar velocity effects in the optically-thin approximation accurately. One is particle-based, the other grid-based, and while the former is most accurate, we demonstrate that the latter is computationally more efficient and can achieve sufficient accuracy. [Abridged]

We studied flux emergence events of sub-granular scale in a solar active region. New Solar Telescope (NST) of Big Bear Solar Observatory made it possible to clearly observe the photospheric signature of flux emergence with very high spatial (0".11 at 7057{\AA}) and temporal (15 s) resolution. From TiO observations with the pixel scale of 0".0375, we found several elongated granule-like features (GLFs) stretching from the penumbral filaments of a sunspot at a relatively high speed of over 4 km s-1. After a slender arched darkening appeared at a tip of a penumbral filament, a bright point (BP) developed and quickly moved away from the filament forming and stretching a GLF. The size of a GLF was approximately 0.5" wide and 3" long. The moving BP encountered nearby structures after several minutes of stretching, and a well-defined elongated shape of a GLF faded away. Magnetograms from SDO/HMI and NST/IRIM revealed that those GLFs are photospheric indicators of small-scale flux emergence, and their disappearance is related to magnetic cancellation. From two well-observed events, we describe detailed development of the sub-structures of GLFs, and different cancellation processes that each of the two GLFs underwent.

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

The Ultra-Fast Flash Observatory (UFFO) Pathfinder is a payload on the Russian Lomonosov satellite, scheduled to be launched in November 2011. The Observatory is designed to detect early UV/Optical photons from Gamma-Ray Bursts (GRBs). There are two telescopes and one main data acquisition system: the UFFO Burst Alert & Trigger Telescope (UBAT), the Slewing Mirror Telescope (SMT), and the UFFO Data Acquisition (UDAQ) system. The UDAQ controls and manages the operation and communication of each telescope, and is also in charge of the interface with the satellite. It will write the data taken by each telescope to the NOR flash memory and sends them to the satellite via the Bus-Interface system (BI). It also receives data from the satellite including the coordinates and time of an external trigger from another payload, and distributes them to two telescopes. These functions are implemented in field programmable gates arrays (FPGA) for low power consumption and fast processing without a microprocessor. The UDAQ architecture, control of the system, and data flow will be presented.

The Ultra-Fast Flash Observatory (UFFO) pathfinder is a space mission devoted to the measurement of Gamma-Ray Bursts (GRBs), especially their early light curves which will give crucial information on the progenitor stars and central engines of the GRBs. It consists of two instruments: the UFFO Burst Alert & Trigger telescope (UBAT) for the detection of GRB locations and the Slewing Mirror Telescope (SMT) for the UV/optical afterglow observations, upon triggering by UBAT. The UBAT employs a coded-mask {\gamma}/X-ray camera with a wide field of view (FOV), and is comprised of three parts: a coded mask, a hopper, and a detector module (DM). The UBAT DM consists of a LYSO scintillator crystal array, multi-anode photo multipliers, and analog and digital readout electronics. We present here the design and fabrication of the UBAT DM, as well as its preliminary test results.

The Ultra-Fast Flash Observatory-Pathfinder (UFFO-P) is to be launched onboard Lomonosov spacecraft in November 2011. It is to measure early UV/Optical photons from Gamma Ray Bursts (GRBs). Slewing Mirror Telescope (SMT) is one of two instruments designed for detection of UV/Optical images of the GRBs. SMT is a Ritchey-Chr\'etien telescope of 100 mm in diameter with a motorized slewing mirror at the entrance providing 17\times17 arcmin2 in Field of View (FOV) and 4 arcsec in pixel resolution. Its sky coverage can be further expanded up to 35 degrees in FOV by tilting a motorized slewing mirror. All mirrors were fabricated to about RMS 0.02 waves in wave front error (WFE) and 84.7% (in average reflectivity) over 200nm~650nm range. SMT was aligned to RMS 0.05 waves in WFE (test wavelength 632.8nm). From the static gravity test result, SMT optics system is expected to survive during launch. The technical details of SMT assembly and laboratory performance test results are reported.

The Ultra-Fast Flash Observatory (UFFO) is a new space-based experiment to observe Gamma-Ray Bursts (GRBs). GRBs are the most luminous electromagnetic events in the universe and occur randomly in any direction. Therefore the UFFO consists of two telescopes; UFFO Burst Alert & Trigger Telescope (UBAT) to detect GRBs using a wide field-of-view (FOV), and a Slewing Mirror Telescope (SMT) to observe UV/optical events rapidly within the narrow, targeted FOV. The SMT is a Ritchey-Chretien telescope that uses a motorized mirror system and an Intensified Charge-Coupled Device (ICCD). When the GRB is triggered by the UBAT, the SMT receives the position information and rapidly tilts the mirror to the target. The ICCD start to take the data within a second after GRB is triggered. Here we give the details about the SMT readout electronics that deliver the data.

Hundreds of gamma-ray burst (GRB) optical light curves have been measured since the discovery of optical afterglows. However, even after nearly 7 years of operation of the Swift Observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at short time scales associated with prompt emissions and progenitors. Because of this lack of sub-minute data, the characteristics of the rise phase of optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching sub-minute and sub-second timescales in a small spacecraft observatory. Rather than slewing the entire spacecraft to aim the optical instrument at the GRB position, we use rapidly moving mirror to redirect our optical beam. As a first step, we employ motorized slewing mirror telescope (SMT), which can point to the event within 1s, in the UFFO Pathfinder GRB Telescope onboard the Lomonosov satellite to be launched in Nov. 2011. UFFO's sub-minute measurements of the optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, provide the first early and detailed measurements of fast-rise GRB optical light curves, and help verify the prospect of GRB as a new standard candle. We will describe the science and the mission of the current UFFO Pathfinder project, and our plan of a full-scale UFFO-100 as the next step.

We detected 2.8 bright points (BPs) per Mm$^2$ in the Quiet Sun (QS) with the New Solar Telescope (NST) at Big Bear Solar Observatory; using the TiO 705.68 nm spectral line, at an angular resolution ~ 0.1'' to obtain 30 min data sequence. Some BPs formed knots that were stable in time and influenced the properties of the granulation pattern around them. The observed granulation pattern within ~ 3'' of knots presents smaller granules than those observed in a normal granulation pattern; i.e., around the knots a suppressed convection is detected. Observed BPs covered ~ 5% of the solar surface and were not homogeneously distributed. BPs had an average size of 0.22'', they were detectable for 4.28 min in average, and had an averaged contrast of 0.1% in the deep red TiO spectral line.

With the ever increasing influx of high resolution images of the solar surface obtained at a multitude of wavelengths, various processes occurring at small spatial scales have become a greater focus of our attention. Complex small-scale magnetic fields have been reported that appear to have enough stored to heat the chromosphere. While significant progress has been made in understanding small-scale phenomena, many specifics remain elusive. We present here a detailed study of a single event of disappearance of a magnetic dipole and associated chromospheric activity. Based on New Solar Telescope H$\alpha$ data and {\it Hinode} photospheric line-of-sight magnetograms and Ca II H images we report the following. 1) Our analysis indicates that even very small dipoles (elements separated by about 0\arcsec.5 or less) may reach the chromosphere and trigger non-negligible chromospheric activity. 2) Careful consideration of the magnetic environment where the new flux is deposited may shed light on the details of magnetic flux removal from the solar surface. We argue that the apparent collision and disappearance of two opposite polarity elements may not necessarily indicate their cancellation (i.e., reconnection, emergence of a "U" tube or submergence of $ \Omega $ loops). In our case, the magnetic dipole disappeared by reconnecting with overlying large-scale inclined plage fields. 3) Bright points seen in off-band H$\alpha$ images are very well-correlated with the Ca II H bright points, which in turn are co-spatial with G-band bright points. We further speculate that, in general, H$\alpha$ bright points are expected be co-spatial with photospheric BPs, however, a direct comparison is needed to refine their relationship.

Surface photometry of the quiet Sun has achieved an angular resolution of $0".1$ with the New Solar Telescope at Big Bear Solar Observatory revealing that a disproportionate fraction of the oscillatory events appear above observed bright point-like structures. During the tracking of these structures, we noted that the more powerful oscillatory events are cospatial with them, indicating that observed flux tubes may be the source of many observed oscillatory events.

We present some preliminary results from a series of extremely large, high-resolution N-body simulations of the formation of early nonlinear structures. We find that the high-z halo mass function is inconsistent with the Sheth-Tormen mass function, which tends to over-estimate the abundance of rare halos. This discrepancy is in rough agreement with previous results based on smaller simulations. We also show that the number density of minihaloes is correlated with local matter density, albeit with a significant scatter that increases with redshift, as minihaloes become increasingly rare. The average correlation is in rough agreement with a simple analytical extended Press-Schechter model, but can differ by up to factor of 2 in some regimes.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

The first, self-consistent calculations are presented of the cosmological, H2-dissociating UV background produced during the epoch of reionization by the sources of reionization. Large-scale radiative transfer simulations of reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down to dwarf galaxies of mass ~10^8 Msun, identified by very high-resolution N-body simulations, including the self-regulating effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the inhomogeneous radiation background in the Lyman-Werner bands of H2 between 11 and 13.6 eV. On average, the intensity of this Lyman-Werner background is found to rise to the threshold level at which dissociation suppresses H2 cooling and star formation inside minihalos, long before reionization is complete. Spatial variations in the Lyman-Werner background are found which result from the clustering of sources associated with large-scale structure formation, such that intensity fluctuations correlate with matter density fluctuations. As a result, the Lyman-Werner background rises to the threshold level for H2 suppression earlier in the vicinity of the reionization sources and their H II regions.

The first, self-consistent calculations of the cosmological H_2 dissociating UV background produced during the epoch of reionization (EOR) by the sources of reionization are presented. Large-scale radiative transfer simulations of reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down to dwarf galaxies of mass ~ 10^8 solar mass, identified by very high-resolution N-body simulations, including the self-regulating effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the inhomogeneous background in the Lyman-Werner band of H_2 between 11 and 13.6 eV.

We study the effect of starlight from the first stars on the ability of other minihaloes in their neighborhood to form additional stars. The question of what the dynamical consequences were for these target minihaloes, of their exposure to the ionizing and dissociating starlight from the Pop III star requires a full radiation-hydrodynamics calculation. Towards this end, we have performed a series of detailed, 1D, radiation-hydrodynamical simulations. We have varied the distance to the source (and, hence, the flux) and the mass and evolutionary stage of the target haloes to quantify this effect. We find: (1) trapping of the I-front and its transformation from R-type to D-type, preceded by a shock front; (2) photoevaporation of the ionized gas (i.e. all gas originally located outside the trapping radius); (3) formation of an H_2 precursor shell which leads the I-front, stimulated by partial photoionization; and (4) the shock-induced formation of H_2 in the minihalo neutral core when the shock speeds up and partially ionizes the gas. The fate of the neutral core is mostly determined by the response of the core to this shock front, which leads to molecular cooling and collapse that, when compared to the same halo without external radiation, is either: (a) expedited, (b) delayed, (c) unaltered, or (d) reversed or prevented, depending upon the flux (i.e. distance to the source) and the halo mass and evolutionary stage. Roughly speaking, most haloes that were destined to cool, collapse, and form stars in the absence of external radiation are found to do so even when exposed to the first Pop III star in their neighborhood, while those that would not have done so are still not able to. (abridged)

New results on the anisotropy of the cosmic microwave background (CMB) and its polarization based on the first 3 years of data from the Wilkinson Microwave Anisotropy Probe (WMAP) have revised the electron scattering optical depth downward from tau_es=0.17+0.08-0.07 to tau_es=0.09+/-0.03. This implies a shift of the effective reionization redshift from z_r~17 to z_r~11. Previous attempts to explain the high redshift of reionization inferred from the WMAP 1 year data have led to widespread speculation that the sources of reionization must have been much more efficient than those associated with the star formation observed at low redshift. This is consistent, for example, with the suggestion that early star formation involved massive, Population III stars that early on produced most of the ionizing radiation escaping from halos. It is therefore tempting to interpret the new WMAP results as implying that we can now relax those previous high demands on the efficiency of the sources of reionization and perhaps even turn the argument around as evidence against such high efficiency. We show that this is not the case, however. The new WMAP results also find that the primordial density fluctuation power spectrum has a lower amplitude, sigma_8, and departs substantially from the scale-invariant spectrum. We show that these effects combine to cancel the impact of the later reionization implied by the new value of tau_es on the required ionizing efficiency per collapsed baryon. The delay of reionization is surprisingly well matched by a comparable delay (by a factor of ~1.4 in scale factor) in the formation of the halos responsible for reionization.

The H atoms inside minihalos (i.e. halos with virial temperatures T_vir < 10^4 K, in the mass range roughly from 10^4 M_sun to 10^8 M_sun) during the cosmic dark ages in a LambdaCDM universe produce a redshifted background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm signal from individual halos of different mass and redshift and the evolving mass function of minihalos to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation at high redshift (z>8) to compute the mean brightness temperature of this background from both minihalos and the intergalactic medium (IGM) prior to the onset of Ly-alpha radiative pumping. We find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse shocked gas in the IGM.

N-body simulations find a universal structure for the halos which result from the nonlinear growth of Gaussian-random-noise density fluctuations in the CDM universe. This talk summarized our attempts to derive and explain this universal structure by analytical approximation and simplified models. As an example, we show here that a 1D spherical infall model involving a fluid approximation derived from the Boltzmann equation can explain not only the halo density profile but its phase-space density profile, as well.

We have derived the first, fully-cosmological, similarity solutions for CDM halo formation in the presence of nongravitational collisionality, which provides an analytical theory of the effect of the self-interacting dark matter (SIDM) hypothesis on halo density profiles. Collisions transport heat inward, which produces a constant-density core, while continuous infall pumps energy into the halo to stabilize the core against gravothermal catastrophe. These solutions improve upon earlier attempts to model the formation and evolution of SIDM halos, offer deeper insight than existing N-body experiments, and yield a more precise determination of the dependence of halo density profile on the value of the CDM self-interaction cross section. Different solutions arise for different values of the dimensionless collisionality parameter Q = s rho_b r_v \~ r_v/l_mfp, where s is the scattering cross section per unit mass, rho_b is the cosmic mean matter density, r_v is halo virial radius and l_mfp is the collision mean free path. The maximum flattening of central density occurs for an intermediate value of Q, Q_th, at which the halo is maximally relaxed to isothermality. The density profiles with constant-density cores preferred by dwarf and LSB rotation curves are best fit by the maximally-flattened (Q=Q_th) solution. If we assume that dwarfs and LSB galaxies formed at their typical collapse epoch in LCDM, then the value of s which makes Q=Q_th is s ~ 200 cm^{2}/g, much higher than previous estimates, s ~ 0.5-5 cm^{2}/g, based on N-body experiments. If s is independent of collision velocity, then the same value s ~ 200 cm^{2}/g would make Q>Q_th for clusters, which typically formed only recently, resulting in relatively less flattening of their central density profile and a smaller core.

The H atoms inside minihalos (i.e. halos with virial temperatures T_vir < 10^4 K, in the mass range roughly from 10^4 M_sun to 10^8 M_sun) during the cosmic dark ages in a LCDM universe produce a redshifted background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm signal from individual halos of different mass and redshift and the evolving mass function of minihalos to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation at high redshift (z > 8) to compute the mean brightness temperature of this background from both minihalos and the intergalactic medium (IGM) prior to the onset of Ly-alpha radiative pumping. We find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse shocked gas in the IGM.

The origin of the cosmic gamma-ray background at 1-20 MeV remains a mystery. We show that gamma-ray emission accompanying annihilation of 20 MeV dark matter particles explains most of the observed signal. Our model satisfies all of the current observational constraints, and naturally provides the origin of "missing" gamma-ray background at 1-20 MeV and 511 keV line emission from the Galactic center. We conclude that gamma-ray observations support the existence of 20 MeV dark matter particles. Improved measurements of the gamma-ray background in this energy band undoubtedly test our proposal.

The observed cosmic gamma-ray background at ~ MeV has often been attributed to Type Ia supernovae (SNIa). Since SNIa is close to a standard candle, one can calculate the gamma-ray intensity of SNIa integrated over redshifts fairly accurately, once the evolution of the SNIa rate is known. The latest SNIa rate measured at z < ~1.6 (Dahlen et al. 2004) indicates that the previous calculations of the gamma-ray background consistently overestimated the SNIa rate. With the new rate, we find that the SNIa contribution is an order of magnitude smaller than observed, and thus new population(s) of sources should be invoked.

Motivated by a recent detection of 511 keV photons from the center of our Galaxy, we calculate the spectrum of the soft gamma-ray background of the redshifted 511 keV photons from cosmological halos. Annihilation of dark matter particles into electron-positron pairs makes a substantial contribution to the gamma-ray background. Mass of such dark matter particles must be <~ 100 MeV so that resulting electron-positron pairs are on-relativistic. On the other hand, we show that in order for the annihilation not to exceed the observed background, the dark matter mass needs to be >~ 20 MeV. We include the contribution from the active galactic nuclei and supernovae. The halo substructures may increase the lower bound to >~ 60 MeV.

Dark-matter halos are the scaffolding around which galaxies and clusters are built. They form when the gravitational instability of primordial density fluctuations causes regions which are denser than average to slow their cosmic expansion, recollapse, and virialize. Understanding the equilibrium structure of these halos is thus a prerequisite for understanding galaxy and cluster formation. Numerical N-body simulations of structure formation from Gaussian-random-noise initial conditions in the CDM universe find a universal internal structure for halos. Objects as different in size and mass as dwarf spheroidal galaxies and galaxy clusters are predicted to have halos with the same basic structure when properly rescaled, independent of halo mass, of the shape of the power spectrum of primordial density fluctuations, and of the cosmological background parameters. This remarkable universality is a fundamental prediction of the CDM model, but our knowledge is limited to the ``empirical'' N-body simulation results, with little analytical understanding. We summarize here our attempts to fill this gap, in an effort to derive and give physical insight to the numerical results and extend them beyond the range of numerical simulation: (1) Simulated halos which form from highly simplified initial conditions involving gravitational instability in a cosmological pancake show that many of the universal properties of CDM halos are generic to cosmological gravitational collapse and do not require Gaussian-random-noise density fluctuations or hierarchical clustering. (2) A fluid approximation derived from the Boltzmann equation yields an analytical theory of halo dynamics which can explain many of the N-body results ... (abridged)

We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their dynamics, which take proper account of the collisional interaction of SIDM particles, based on a fluid approximation derived from the Boltzmann equation. These similarity solutions are relevant to galactic and cluster halo formation in the cold dark matter (CDM) model. Different solutions arise for different values of the dimensionless collisionality parameter, Q=sigma*rho_b*r_s, where sigma is the SIDM particle scattering cross section per unit mass, rho_b is the cosmic mean density, and r_s is the shock radius. For all solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of r_s, pumped by cosmological infall. Accordingly, core collapse must in general be delayed until infall becomes negligible, contrary to previous analyses based on isolated halos, which predict core collapse in a Hubble time. Our solutions agree with N-body simulations, which match observed galactic rotation curves if sigma=[0.56 - 5.6]cm^2/g, implying Q=[6.2e-7 - 3.6e-5]. Similar profiles also arise for Q=[1.37e-2 - 1.7e-1], or sigma=[1.2e4 - 2.7e4]cm^2/g, a regime not previously simulated.

We use the universal mass accretion history recently reported for simulations of halo formation in the cold dark matter model (CDM) to analyze the formation and growth of a single halo. We derive the time-dependent density profile three different ways, based upon three approximations of successively greater realism: equilibrium, radial orbits, and a fluid approximation. For the equilibrium model, the density profile is well-fit by either an NFW or Moore profile over a limited range of radii and scale factors. For the radial orbit model, we find profiles which are generally steeper than the NFW profile, with an inner logarithmic slope approaching -2, consistent with a purely radial collisionless system. In the fluid approximation, we find good agreement with the NFW and Moore profiles for radii resolved by N-body simulations (r/r_200 > 0.01), and an evolution of concentration parameter nearly identical to that found in N-body simulations. The evolving structure of cosmological halos is therefore best understood as the effect of a time-varying rate of mass infall on a smoothly distributed, isotropic, collisionless fluid.

We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions taking account of the collisional interaction of SIDM particles. This interaction results in a thermal conductivity that heats the halo core and flattens its density profile. These similarity solutions are relevant to galactic and cluster halo formation in the CDM model. We assume an initial mass profile dM/M M^{-eps}, as in the familiar secondary infall model. If eps=1/6, SIDM halos will evolve self-similarly, with a cold, supersonic infall terminated by a strong accretion shock. Different solutions arise for different values of the collisionality parameter, Q= sigma rho_b r_s, where sigma is the scattering cross section, rho_b is the cosmic mean density, and r_s is the shock radius. For all these solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of r_s. We find two different regimes for these solutions: 1) for Q \leq Q_{th}, the core density decreases and core size increases as Q increases; 2) for Q \geq Q_{th}, the core density increases and core size decreases as Q increases. Our similarity solutions are in agreement with previous N-body simulations of SIDM halos, which correspond to the low-Q regime, if Q=[8.4e-4 - 4.9e-2]Q_{th} (low-Q), or sigma=[0.56-5.6]cm^2/g. As Q=\infty, our similarity solution aquires a central density cusp, in agreement with some simulation results which used an ordinary collisional fluid to approximate the effects of SIDM collisionality. When Q=[18.6-231]Q_{th} or sigma=[1.2e4 - 2.71e4]cm^2/g, for which we find flat-density cores comparable to those of the observationally acceptable low-Q solutions, has not previously been identified. Further study of this regime is warranted.