We study prospects of constraining the primordial magnetic field (PMF) and its evolution during the dark ages and cosmic dawn in light of EDGES 21-cm signal. Our analysis has been carried out on a `colder IGM' background which is one of the promising avenues to interpret the EDGES signal. We consider the dark matter-baryon interactions for the excess cooling. We find that the colder IGM suppresses both the residual free electron fraction and the coupling coefficient between the ionised and neutral components. The Compton heating also gets affected in colder IGM background. Consequently, the IGM heating rate due to the PMF enhances compared to the standard scenario. Thus, a significant fraction of the magnetic energy, for $B_0 \lesssim 0.5 \, {\rm nG}$, gets transferred to the IGM and the magnetic field decays at a much faster rate compared to the simple $(1+z)^2$ scaling during the dark ages and cosmic dawn. This low PMF is an unlikely candidate for explaining the rise of the EDGES absorption signal at lower redshift. We also see that the PMF and DM-baryon interaction together introduces a plateau-like feature in the redshift evolution of the IGM temperature. We find that the upper limit on the PMF depends on the underlying DM-baryon interaction. Higher PMF can be allowed when the interaction cross-section is higher and/or the DM particle mass is lower. Our study shows that the PMF with $B_0$ up to $\sim 0.4 \, {\rm nG}$, which is ruled out in the standard model, can be allowed if DM-baryon interaction with suitable cross-section and DM mass is considered.

We present semi-analytical models of high redshift Lyman-{\alpha} emitters (LAEs) in order to constrain the star formation efficiency in those galaxies. Our supernova feedback induced star formation model along with Sheth-Tormman halo mass function correctly reproduces the shape, amplitude and the redshift evolution of UV and Lyman-{\alpha} luminosity functions of LAEs in the redshift range z = 2 to 7.3. We show that the fraction of Lyman-{\alpha} emitting galaxies increases with increasing redshifts reaching to unity just after the reionisation. However, we show that star formation efficiency in those LAEs does not show any redshift evolution within the uncertainty in available observations. This would have significant repercussion on the reionisation of the intergalactic medium.

We fit the recently updated UV luminosity functions (LF) of high-$z$ ($1.5\le z\le8.0$) galaxies using our semi-analytical models of galaxy formation that take into account various feedback processes. In order to reproduce the overall redshift evolution we require the efficiency of converting gas into stars to decrease with decreasing redshift. Even for $z\ge 6$, our models require supernovae (SNe) feedback to reproduce the observed LF suggesting the prevalence of galactic winds that could have polluted the inter-galactic medium even at very high redshifts. The observed LF in the low luminosity end for $z<2.5$ shows an upward turn. In our models we reproduce this trend using passively evolving population of galaxies. Measuring stellar mass, age and metallicity of these galaxies using multi-band spectral energy distribution (SED) fitting will place strong constraints on the existence of such galaxies. While Active galactic Nuclei (AGN) feedback is essential to reproduce the LF at high luminosity end for $z<4$, it may not needed in the case of $z\ge 6$. We show that the expected turn around in the LF due to cooling criteria occurs at luminosity much lower than what has been probed with present day observations. With future deep observations, that can measure the LF more accurately, we will be able to distinguish between different modes of SNe feedback and get insights into the physical processes that drive the galaxy evolution.

We present semi-analytical models of galactic outflows in high redshift galaxies driven by both hot thermal gas and non-thermal cosmic rays. Thermal pressure alone may not sustain a large scale outflow in low mass galaxies (i.e $M\sim 10^8$~M$_\odot$), in the presence of supernovae (SNe) feedback with large mass loading. We show that inclusion of cosmic ray pressure allows outflow solutions even in these galaxies. In massive galaxies for the same energy efficiency, cosmic ray driven winds can propagate to larger distances compared to pure thermally driven winds. On an average gas in the cosmic ray driven winds has a lower temperature which could aid detecting it through absorption lines in the spectra of background sources. Using our constrained semi-analytical models of galaxy formation (that explains the observed UV luminosity functions of galaxies) we study the influence of cosmic ray driven winds on the properties of the intergalactic medium (IGM) at different redshifts. In particular, we study the volume filling factor, average metallicity, cosmic ray and magnetic field energy densities for models invoking atomic cooled and molecular cooled halos. We show that the cosmic rays in the IGM could have enough energy that can be transferred to the thermal gas in presence of magnetic fields to influence the thermal history of the intergalactic medium. The significant volume filling and resulting strength of IGM magnetic fields can also account for recent $\gamma$-ray observations of blazars.

Studying the cosmic dawn and the epoch of reionization through the redshifted 21 cm line are among the major science goals of the SKA1. Their significance lies in the fact that they are closely related to the very first stars in the universe. Interpreting the upcoming data would require detailed modelling of the relevant physical processes. In this article, we focus on the theoretical models of reionization that have been worked out by various groups working in India with the upcoming SKA in mind. These models include purely analytical and semi-numerical calculations as well as fully numerical radiative transfer simulations. The predictions of the 21 cm signal from these models would be useful in constraining the properties of the early galaxies using the SKA data.

We present an improved photometric redshift estimator code, CuBAN$z$, that is publicly available at https://goo.gl/fpk90V}{https://goo.gl/fpk90V. It uses the back propagation neural network along with clustering of the training set, which makes it more efficient than existing neural network codes. In CuBAN$z$, the training set is divided into several self learning clusters with galaxies having similar photometric properties and spectroscopic redshifts within a given span. The clustering algorithm uses the color information (i.e. $u-g$, $g-r$ etc.) rather than the apparent magnitudes at various photometric bands as the photometric redshift is more sensitive to the flux differences between different bands rather than the actual values. Separate neural networks are trained for each cluster using all possible colors, magnitudes and uncertainties in the measurements. For a galaxy with unknown redshift, we identify the closest possible clusters having similar photometric properties and use those clusters to get the photometric redshifts using the particular networks that were trained using those cluster members. For galaxies that do not match with any training cluster, the photometric redshifts are obtained from a separate network that uses entire training set. This clustering method enables us to determine the redshifts more accurately. SDSS Stripe 82 catalog has been used here for the demonstration of the code. For the clustered sources with redshift range $z_{\rm spec}<0.7$, the residual error ($\langle (z_{{\rm spec}}-z_{{\rm phot}})^2\rangle^{1/2} $) in the training/testing phase is as low as 0.03 compared to the existing ANNz code that provides residual error on the same test data set of 0.05. Further, we provide a much better estimate of the uncertainty of the derived photometric redshift.

We present results of our time variability studies of Mg II and Al III absorption lines in a sample of 22 Low Ionization Broad Absorption Line QSOs (LoBAL QSOs) at 0.2 <= zem <= 2.1 using the 2m telescope at IUCAA Girawali Observatory over a time-scale of 10 days to 7.69 years in the QSO's rest frame. Spectra are analysed in conjunction with photometric light curves from Catalina Real-Time Transient Survey. Long time-scale (i.e >= 1 year) absorption line variability is seen in 8 cases (36% systems) while only 4 of them (i.e 18% systems) show variability over short time-scales (i.e < 1 year). We notice a tendency of highly variable LoBAL QSOs to have high ejection velocity, low equivalent width and low redshift. The detection rate of variability in LoBAL QSOs showing Fe fine-structure lines (FeLoBAL QSOs) is less than that seen in non-Fe LoBAL QSOs. Absorption line variability is more frequently detected in QSOs having continuum dominated by Fe emission lines compared to rest of the QSOs. Confirming these trends with a bigger sample will give vital clues for understanding the physical distinction between different BAL QSO sub-classes. We correlate the absorption line variability with various parameters derived from continuum light curves and find no clear correlation between continuum flux and absorption line variabilities. However, sources with large absorption line variability also show large variability in their light curves. We also see appearance/disappearance of absorption components in 2 cases and clear indications for profile variations in 4 cases. The observed variability can be best explained by a combination of process driven by continuum variations and clouds transiting across the line of sight.

We present a semi-analytical model of high redshift galaxy formation. In our model the star formation inside a galaxy is regulated by the feedback from supernova (SNe) driven outflows. We derive a closed analytical form for star formation rate in a single galaxy taking account of the SNe feedback in a self-consistent manner. We show that our model can explain the observed correlation between the stellar mass and the circular velocity of galaxies from dwarf galaxies to massive galaxies of $10^{12} M_\odot$. For small mass dwarf galaxies additional feedback other than supernova feedback is needed to explain the spread in the observational data. Our models reproduce the observed 3-D fundamental correlation between the stellar mass, gas phase metallicity and star formation rate in galaxies establishing that the SNe feedback plays a major role in building this relation. Further, the observed UV luminosity functions of Lyman-Break galaxies (LBGs) are well explained by our feedback induced star formation model for a vast redshift range of $1.5 \le z \le 8$. In particular, the flattening of the luminosity functions at the low luminosity end naturally arises due to our explicit SNe feedback treatment.

We present a physically motivated semi-analytic model to understand the clustering of high redshift LBGs. We show that the model parameters constrained by the observed luminosity function, can be used to predict large scale (\theta > 80 arcsec) bias and angular correlation function of galaxies. These predictions are shown to reproduce the observations remarkably well. We then adopt these model parameters to calculate the halo occupation distribution (HOD) using the conditional mass function. The halo model using this HOD is shown to provide a reasonably good fit to the observed clustering of LBGs at both large (\theta>80") and small (\theta < 10") angular scales for z=3-5 and several limiting magnitudes. However, our models underpredict the clustering amplitude at intermediate angular scales, where quasi-linear effects are important. The average mass of halos contributing to the observed clustering is found to be 6.2 x 10^{11} M_\odot and the characteristic mass of a parent halo hosting satellite galaxies is 1.2 \times 10^{12} M_\odot for a limiting absolute magnitude of -20.5 at z=4. For a given threshold luminosity these masses decrease with increasing z and at any given z these are found to increase with increasing value of threshold luminosity. We find that approximately 40 % of the halos above a minimum mass M_{min}, can host detectable central galaxies and about 5-10 % of these halos are likely to also host a detectable satellite. The satellites form typically a dynamical timescale prior to the formation of the parent halo. The small angular scale clustering is due to central-satellite pairs and is quite sensitive to changes in the duration of star formation in a halo. The present data favor star formation in a halo lasting typically for a few dynamical time-scales. Our models also reproduce different known trends between parameters related to star formation.

Laboratory experiments measuring neutrino oscillations, indicate small mass differences between different mass eigenstates of neutrinos. The absolute mass scale is however not determined, with at present the strongest upper limits coming from astronomical observations rather than terrestrial experiments. The presence of massive neutrinos suppresses the growth of perturbations below a characteristic mass scale, thereby leading to a decreased abundance of collapsed dark matter halos. Here we show that this effect can significantly alter the predicted luminosity function (LF) of high redshift galaxies. In particular we demonstrate that a stringent constraint on the neutrino mass can be obtained using the well measured galaxy LF and our semi-analytic structure formation models. Combining the constraints from the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) data with the LF data at z = 4, we get a limit on the sum of the masses of 3 degenerate neutrinos \Sigma m_\nu < 0.52 eV at the 95 % CL. The additional constraints using the prior on Hubble constant strengthens this limit to \Sigma m_\nu < 0.29 eV at the 95 % CL. This neutrino mass limit is a factor of order 4 improvement compared to the constraint based on the WMAP7 data alone, and as stringent as known limits based on other astronomical observations. As different astronomical measurements may suffer from different set of biases, the method presented here provides a complementary probe of \Sigma m_\nu . We suggest that repeating this exercise on well measured luminosity functions over different redshift ranges can provide independent and tighter constraints on \Sigma m_\nu .

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

We study winds in high redshift galaxies driven by a relativistic cosmic ray (proton) component in addition to the hot thermal gas component. Cosmic rays (CRs) are likely to be efficiently generated in supernova shocks inside galaxies. We obtain solutions of such CR driven free winds in a gravitational potential of the NFW form, relevant to galaxies. Cosmic rays naturally provide the extra energy and/or momentum input to the system, needed for a transonic wind solution in a gas with adiabatic index $\gamma=5/3$. We show that CRs can effectively drive winds even when the thermal energy of the gas is lost due to radiative cooling. These wind solutions predict an asymptotic wind speed closely related to the circular velocity of the galaxy. Furthermore, the mass outflow rate per unit star formation rate (eta_w) is predicted to be ~ 0.2-0.5 for massive galaxies, with masses $M \sim 10^{11}-10^{12} M_\odot$. We show eta_w to be inversely proportional to the square of the circular velocity. Magnetic fields at the $\mu$G levels are also required in these galaxies to have a significant mass loss. A large eta_w for small mass galaxies implies that CR driven outflows could provide a strong negative feedback to the star formation in dwarf galaxies. Further, our results will also have important implications to the metal enrichment of the IGM. These conclusions are applicable to the class of free wind models where the source region is confined to be within the sonic point.

We present a semi-analytical model of star formation which explains simultaneously the observed UV luminosity function of high redshift Lyman break galaxies (LBGs) and luminosity functions of Lyman-alpha emitters. We consider both models that use the Press-Schechter (PS) and Sheth-Tormen (ST) halo mass functions to calculate the abundances of dark matter halos. The Lyman-alpha luminosity functions at z < 4 are well reproduced with only <10% of the LBGs emitting Lyman-alpha lines with rest equivalent width greater than the limiting equivalent width of the narrow band surveys. However, the observed luminosity function at z > 5 can be reproduced only when we assume that nearly all LBGs are Lyman-alpha emitters. Thus it appears that 4 < z < 5 marks the epoch when a clear change occurs in the physical properties of the high redshift galaxies. As Lyman-alpha escape depends on dust and gas kinematics of the inter stellar medium (ISM), this could mean that on an average the ISM at z > 5 could be less dusty, more clumpy and having more complex velocity field. All of these will enable easier escape of the Lyman-alpha photons. At z > 5 the observed Lyman-alpha luminosity function are well reproduced with the evolution in the halo mass function along with very minor evolution in the physical properties of high redshift galaxies. In particular, upto z = 6.5, we do not see the effect of evolving inter galactic medium (IGM) opacity on the Lyman-alpha escape from these galaxies.

The form of the halo mass function is a basic ingredient in any semi-analytical galaxy formation model. We study the existing forms of the mass functions in the literature and compare their predictions for semi-analytical galaxy formation models. Two methods are used in the literature to compute the net formation rate of halos, one by simply taking the derivative of the halo mass function and the other using the prescription due to Sasaki (1994). For the Press-Schechter (PS) mass function, we compare various model predictions, using these two methods. However, as the Sasaki formalism cannot be easily generalized for other mass functions, we use the derivative while comparing model predictions of different mass functions. We show that the reionization history and UV luminosity function of Lyman break galaxies (LBGs) predicted by the PS mass function differs from those using any other existing mass function, like Sheth-Tormen (ST) mass function.In particular the reionization efficiency of molecular cooled halos has to be substantially reduced when one uses the ST and other mass functions obtained from the simulation instead of the PS mass function. Using $\chi^2$-minimization, we find that the observed UV luminosity functions of LBGs at $3.0\le z\le 7.4$ are better reproduced by models using the ST mass function compared to models that use the PS mass function. On the other hand, the volume filling factor of the metals expelled from the galaxies through supernovae driven outflows differs very little between models with different mass functions. It depends on the way we treat merging outflows. We also show that the porosity weighted average quantities related to the outflow are not very sensitive to the differences in the halo mass function.

We present semi-analytic models of galactic outflows, constrained by available observations on high redshift star formation and reionization. Galactic outflows are modeled in a manner akin to models of stellar wind blown bubbles. Large scale outflows can generically escape from low mass halos (M<10^9 M_sun) for a wide range of model parameters but not from high mass halos (M> 10^{11} M_sun). The gas phase metallicity of the outflow and within the galaxy are computed. Ionization states of different metal species are calculated and used to examine the detectability of metal lines from the outflows. The global influence of galactic outflows is also investigated. Models with only atomic cooled halos significantly fill the IGM at z~3 with metals (with -2.5>[Z/Z_sun]>-3.7), the actual extent depending on the efficiency of winds, the IMF, the fractional mass that goes through star formation and the reionization history of the universe. In these models, a large fraction of outflows at z~3 are supersonic, hot (T> 10^5 K) and have low density, making metal lines difficult to detect. They may also result in significant perturbations in the IGM gas on scales probed by the Lyman-alpha forest. On the contrary, models including molecular cooled halos with a normal mode of star formation can potentially volume fill the universe at z> 8 without drastic dynamic effects on the IGM, thereby setting up a possible metallicity floor (-4.0<[Z/Z_sun]<-3.6). Interestingly, molecular cooled halos with a ``top-heavy'' mode of star formation are not very successful in establishing the metallicity floor because of the additional radiative feedback, that they induce. (Abridged)

We present a self-consistent, semi-analytical LCDM model of star formation and reionization. For the cosmological parameters favored by the WMAP data, our models consistently reproduce the optical depth to reionization, redshift of reionization and the observed luminosity functions (LF) and hence the star formation rate (SFR) density at 3<z<6 for a reasonable range of model parameters. While simple photoionization feedback produces the correct shape of LF at z = 6, for $z = 3$ we need additional feedback that suppresses star formation in halos with 10<log(M/M_\odot)<11. Models with prolonged continuous star formation activities are preferred over those with short bursts as they are consistent with the existence of a Balmer break in considerable fraction of observed galaxies even at z~6. The halo number density evolution from the standard LCDM structure formation model that fits LF up to z=6 is consistent with the upper limits on z~7 LF and source counts at 8<z<12 obtained from the HUDF observations without requiring any dramatic change in the nature of star formation. However, to reproduce the observed LF at 6<z<10, obtained from the near-IR observations around strong lensing clusters, we need a strong evolution in the IMF, reddening correction and the mode of star formation at z>8. Molecular cooled halos, which may be important for reionizing the universe, are not detectable in the present deep field observations. However, their presence and contribution to reionization can be inferred indirectly from the redshift evolution of the LF at 6<z< 12.Accurately measuring the LF at high z can be used to understand the nature of star formation in the dark ages and probe the history of reionization. (Abridged)

Supernova explosions in the early star forming galaxies will accelerate cosmic rays (CRs). CRs are typically confined in the collapsed objects for a short period before escaping into the intergalactic medium (IGM). Galactic outflows can facilitate this escape by advecting CRs into the IGM. An outflow that results in a termination shock can also generate more CRs. We show that the CR protons from the above processes can significantly affect the thermal history of the IGM. Within plausible range of parameters, cosmic ray heating can compensate for adiabatic cooling and explain the measured IGM temperature at redshifts z between 2 to 4, even with early reionization.