Parker Solar Probe (PSP), NASA's latest and closest mission to the Sun, is on a journey to investigate fundamental enigmas of the inner heliosphere. This paper reports initial observations made by the Solar Probe Analyzer for Ions (SPAN-I), one of the instruments in the Solar Wind Electrons Alphas and Protons (SWEAP) instrument suite. We address the presence of secondary proton beams in concert with ion-scale waves observed by FIELDS, the electromagnetic fields instrument suite. We show two events from PSP's 2nd orbit that demonstrate signatures consistent with wave-particle interactions. We showcase 3D velocity distribution functions (VDFs) measured by SPAN-I during times of strong wave power at ion-scales. From an initial instability analysis, we infer that the VDFs departed far enough away from local thermodynamic equilibrium (LTE) to provide sufficient free energy to locally generate waves. These events exemplify the types of instabilities that may be present and, as such, may guide future data analysis characterizing and distinguishing between different wave-particle interactions.

Electrostatic analyzers of different designs have been used since the earliest days of the space age, beginning with the very earliest solar wind measurements made by Mariner 2 en route to Venus in 1962. The Parker Solar Probe (PSP) mission, NASA's first dedicated mission to study the innermost reaches of the heliosphere, makes its thermal plasma measurements using a suite of instruments called the Solar Wind Electrons, Alphas, and Protons (SWEAP) investigation. SWEAP's electron Parker Solar Probe Analyzer (SPAN-E) instruments are a pair of top-hat electrostatic analyzers on PSP that are capable of measuring the electron distribution function in the solar wind from 2 eV to 30 keV. For the first time, in-situ measurements of thermal electrons provided by SPAN-E will help reveal the heating and acceleration mechanisms driving the evolution of the solar wind at the points of acceleration and heating, closer than ever before to the Sun. This paper details the design of the SPAN-E sensors and their operation, data formats, and measurement caveats from Parker Solar Probe's first two close encounters with the Sun.

The metallicity of strong HI systems, spanning from damped Lyman-alpha absorbers (DLAs) to Lyman-limit systems (LLSs) is explored between z = 5->0 using the EAGLE high-resolution cosmological hydrodynamic simulation of galaxy formation. The metallicities of LLSs and DLAs steadily increase with time in agreement with observations. DLAs are more metal rich than LLSs, although the metallicities in the LLS column density range (NHI = 10^17 -10^20 cm^-2) are relatively flat, evolving from a median HI-weighted metallicity of Z<10^-2 Zsol at z = 3 to ~10^-0.5 Zsol by z = 0. The metal content of HI systems tracks the increasing stellar content of the Universe, holding ~5% of the integrated total metals released from stars at z = 0. We also consider partial LLS (pLLS, NHI = 10^16-10^17 cm^-2) metallicities, and find good agreement with Wotta et al. (2016) for the fraction of systems above (37%) and below (63%) 0.1 Zsol. We also find a large dispersion of pLLS metallicities, although we do not reproduce the observed metallicity bimodality and instead we make the prediction that a larger sample will yield more pLLSs around 0.1 Zsol. We under-predict the median metallicity of strong LLSs, and predict a population of Z < 10^-3 Zsol DLAs at z > 3 that are not observed, which may indicate more widespread early enrichment in the real Universe compared to EAGLE.

The Radiation Assessment Detector (RAD), on board Mars Science Laboratory's (MSL) Curiosity rover, has been measuring ground level particle fluxes along with the radiation dose rate at the surface of Mars since August 2012. Similar to neutron monitors at Earth, RAD sees many Forbush decreases (FDs) in the galactic cosmic ray (GCR) induced surface fluxes and dose rates. These FDs are associated with coronal mass ejections (CMEs) and/or stream/corotating interaction regions (SIRs/CIRs). Orbiting above the Martian atmosphere, the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft has also been monitoring space weather conditions at Mars since September 2014. The penetrating particle flux channels in the Solar Energetic Particle (SEP) instrument onboard MAVEN can also be employed to detect FDs. For the first time, we study the statistics and properties of a list of FDs observed in-situ at Mars, seen both on the surface by MSL/RAD and in orbit detected by the MAVEN/SEP instrument. Such a list of FDs can be used for studying interplanetary CME (ICME) propagation and SIR evolution through the inner heliosphere. The magnitudes of different FDs can be well-fitted by a power-law distribution. The systematic difference between the magnitudes of the FDs within and outside the Martian atmosphere may be mostly attributed to the energy-dependent modulation of the GCR particles by both the pass-by ICMEs/SIRs and the Martian atmosphere.

We use the Aurora radiation-hydrodynamical simulations to study the mean free path (MFP) for hydrogen ionizing photons during the epoch of reionization. We directly measure the MFP by averaging the distance 1 Ry photons travel before reaching an optical depth of unity along random lines-of-sight. During reionization the free paths tend to end in neutral gas with densities near the cosmic mean, while after reionizaton the end points tend to be overdense but highly ionized. Despite the increasing importance of discrete, over-dense systems, the cumulative contribution of systems with $N_{\rm{HI}} \lesssim 10^{16.5}~{\rm cm^{-2}}$ suffices to drive the MFP at $z \approx 6$, while at earlier times higher column densities are more important. After reionization the typical size of HI systems is close to the local Jeans length, but during reionization it is much larger. The mean free path for photons originating close to galaxies, $\rm{MFP_{gal}}$, is much smaller than the cosmic MFP. After reionization this enhancement can remain significant up to starting distances of $\sim 1$ comoving Mpc. During reionization, however, $\rm{MFP_{gal}}$ for distances $\sim 10^2 - 10^3$ comoving kpc typically exceeds the cosmic MFP. These findings have important consequences for models that interpret the intergalactic MFP as the distance escaped ionizing photons can travel from galaxies before being absorbed and may cause them to under-estimate the required escape fraction from galaxies, and/or the required emissivity of ionizing photons after reionization.

We introduce a new suite of radiation-hydrodynamical simulations of galaxy formation and reionization called Aurora. The Aurora simulations make use of a spatially adaptive radiative transfer technique that lets us accurately capture the small-scale structure in the gas at the resolution of the hydrodynamics, in cosmological volumes. In addition to ionizing radiation, Aurora includes galactic winds driven by star formation and the enrichment of the universe with metals synthesized in the stars. Our reference simulation uses 2x512^3 dark matter and gas particles in a box of size 25 comoving Mpc/h with a force softening scale of at most 0.28 kpc/h. It is accompanied by simulations in larger and smaller boxes and at higher and lower resolution, employing up to 2x1024^3 particles, to investigate numerical convergence. All simulations are calibrated to yield simulated star formation rate (SFR) functions in close agreement with observational constraints at redshift z = 7 and to achieve reionization at z = 8.3, which is consistent with the observed optical depth to reionization. We focus on the design and calibration of the simulations and present some first results. The median stellar metallicities of low-mass galaxies at z = 6 are consistent with the metallicities of dwarf galaxies in the Local Group, which are believed to have formed most of their stars at high redshifts. After reionization, the mean photoionization rate decreases systematically with increasing resolution. This coincides with a systematic increase in the abundance of neutral hydrogen absorbers in the IGM.

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

High-resolution simulations of star-forming massive galactic discs have shown that clumps form with a characteristic baryonic mass in the range $10^7-10^8~M_{\odot}$, with a small tail exceeding $10^9~M_{\odot}$ produced by clump-clump mergers. This is in contrast with the observed kpc-size clumps with masses up to $10^{10}~M_{\odot}$ in high-redshift star-forming galaxies. In this paper we show that the comparison between simulated and observed star-forming clumps is hindered by limited observational spatial resolution and sensitivity. We post-process high-resolution hydrodynamical simulations of clumpy discs using accurate radiative transfer to model the effect of ionizing radiation from young stars and to compute H$\alpha$ emission maps. By comparing the intrinsic clump size and mass distributions with those inferred from convolving the H$\alpha$ maps with different gaussian apertures, we mimick the typical resolution used in observations. We found that with 100 pc resolution, mock observations can recover the intrinsic clump radii and stellar masses, in agreement with those found by lensing observations. Instead, using a 1 kpc resolution smears out individual clumps, resulting in their apparent merging. This causes significant overestimations of the clump radii and, therefore, masses derived using methods that use their observed sizes. We show that limited sensitivity can also force observations to significantly overestimate the clump masses. We conclude that a significant fraction of giant clumps detected in the observations may result from artificially inflated radii and masses, and that $\approx 100$ pc spatial resolution is required to capture correctly the physical characteristics of star-forming clumps if they are coherent structures produced by disc fragmentation.

We present a sample of 45 Damped Lyman-$\alpha$ system (DLA, \nhi\ $ \geq 2 \times 10^{20} {\rm cm}^{-2}$) counterparts (33 detections, 12 upper limits) which host gamma-ray bursts (GRB-DLAs) in order to investigate star-formation and metallicity within galaxies hosting DLAs. Our sample spans $z \sim 2-6$ and is nearly three times larger than any previously detected DLA counterparts survey based on quasar line-of-sight searches (QSO-DLAs). We report star formation rates (SFRs) from rest-frame UV photometry and SED modeling. We find that DLA counterpart SFRs are not correlated with either redshift or HI column density. Thanks to the combination of \hst\ and ground-based observations, we also investigate DLA host star-formation efficiency. Our GRB-DLA counterpart sample spans both higher efficiency and low efficiency star formation regions compared to the local Kennicutt-Schmidt relation, local star formation laws, and $z \sim 3$ cosmological simulations. We also compare the depletion times of our DLA hosts sample to other objects in the local Universe, our sample appears to deviate from the star formation efficiencies measured in local spiral and dwarf galaxies. Furthermore, we find similar efficiencies as local inner disks, SMC, and LBG outskirts. Finally, our enrichment time measurements show a spread of systems with under- and over-abundance of metals which may suggest that these systems had episodic star formation and a metal enrichment/depletion as a result of strong stellar feedback and/or metal inflow/outflow.

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

We study the distribution and evolution of highly ionised intergalactic metals in the Evolution and Assembly of Galaxies and their Environment (EAGLE) cosmological, hydrodynamical simulations. EAGLE has been shown to reproduce a wide range of galaxy properties while its subgrid feedback was calibrated without considering gas properties. We compare the predictions for the column density distribution functions (CDDFs) and cosmic densities of SiIV, CIV, NV, OVI and NeVIII absorbers with observations at redshift z = 0 to ~ 6 and find reasonable agreement, although there are some differences. We show that the typical physical densities of the absorbing gas increase with column density and redshift, but decrease with the ionization energy of the absorbing ion. The typical metallicity increases with both column density and time. The fraction of collisionally ionized metal absorbers increases with time and ionization energy. While our results show little sensitivity to the presence or absence of AGN feedback, increasing/decreasing the efficiency of stellar feedback by a factor of two substantially decreases/increases the CDDFs and the cosmic densities of the metal ions. We show that the impact of the efficiency of stellar feedback on the CDDFs and cosmic densities is largely due to its effect on the metal production rate. However, the temperatures of the metal absorbers, particularly those of strong OVI, are directly sensitive to the strength of the feedback.

The use of background quasars provides a powerful tool to probe the cool gas in the circum-galactic medium of foreground galaxies. Here, we present new observations with SINFONI and X-Shooter of absorbing-galaxy candidates at z=0.7-1. We report the detection with both instruments of the H-alpha emission line of one sub-DLA at z_abs=0.94187 with log N(HI)=19.38^+0.10_-0.15 towards SDSS J002133.27+004300.9. We estimate the star formation rate: SFR=3.6+/-2.2 solar masses per year in that system. A detailed kinematic study indicates a dynamical mass M_dyn=10^9.9+/-0.4 solar masses and a halo mass M_halo=10^11.9+/-0.5 solar masses. In addition, we report the OII detection with X-Shooter of another DLA at z_abs=0.7402 with log N(HI)=20.4+/-0.1 toward Q0052+0041 and an estimated SFR of 5.3+/-0.7 solar masses per year. Three other objects are detected in the continuum with X-Shooter but the nature and redshift of two of these objects are unconstrained due to the absence of emission lines, while the third object might be at the redshift of the quasar. We use the objects detected in our whole N(HI)-selected SINFONI survey to compute the metallicity difference between the galaxy and the absorbing gas, delta_HI(X), where a positive (negative) value indicates infall (outflow). We compare this quantity with the quasar line of sight alignment with the galaxy's major (minor) axis, another tracer of infall (outflow). We find that these quantities do not correlate as expected from simple assumptions. Additional observations are necessary to relate these two independent probes of gas flows around galaxies.

The observed high covering fractions of neutral hydrogen (HI) with column densities above $\sim 10^{17} \rm{cm}^{-2}$ around Lyman-Break Galaxies (LBGs) and bright quasars at redshifts z ~ 2-3 has been identified as a challenge for simulations of galaxy formation. We use the EAGLE cosmological, hydrodynamical simulation, which has been shown to reproduce a wide range of galaxy properties and for which the subgrid feedback was calibrated without considering gas properties, to study the distribution of HI around high-redshift galaxies. We predict the covering fractions of strong HI absorbers ($N_{\rm{HI}} \gtrsim 10^{17} \rm{cm}^{-2}$) inside haloes to increase rapidly with redshift but to depend only weakly on halo mass. For massive ($M_{200} \gtrsim 10^{12} {\rm M_{\odot}}$) halos the covering fraction profiles are nearly scale-invariant and we provide fitting functions that reproduce the simulation results. While efficient feedback is required to increase the HI covering fractions to the high observed values, the distribution of strong absorbers in and around halos of a fixed mass is insensitive to factor of two variations in the strength of the stellar feedback. In contrast, at fixed stellar mass the predicted HI distribution is highly sensitive to the feedback efficiency. The fiducial EAGLE simulation reproduces both the observed global column density distribution function of HI and the observed radial covering fraction profiles of strong HI absorbers around LBGs and bright quasars.

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

The density of the warm ionized gas in high-redshift galaxies is known to be higher than what is typical in local galaxies on similar scales. At the same time, the mean global properties of the high- and low-redshift galaxies are quite different. Here, we present a detailed differential analysis of the ionization parameters of 14 star-forming galaxies at redshift 2.6-3.4, compiled from the literature. For each of those high-redshift galaxies, we construct a comparison sample of low-redshift galaxies closely matched in specific star formation rate (sSFR) and stellar mass, thus ensuring that their global physical conditions are similar to the high-redshift galaxy. We find that the median log [OIII] 5007/ [OII] 3727 line ratio of the high-redshift galaxies is 0.5 dex higher than their local counterparts. We construct a new calibration between the [OIII] 5007/ [OII] 3727 emission line ratio and ionization parameter to estimate the difference between the ionization parameters in the high and low-redshift samples. Using this, we show that the typical density of the warm ionized gas in star-forming regions decreases by a median factor of $7.1^{+10.2}_{-5.4}$ from z ~ 3.3 to z ~ 0 at fixed mass and sSFR. We show that metallicity differences cannot explain the observed density differences. Because the high- and low-redshift samples are comparable in size, we infer that the relationship between star formation rate density and gas density must have been significantly less efficient at z ~2-3 than what is observed in nearby galaxies with similar levels of star formation activity.

We combine cosmological, hydrodynamical simulations with accurate radiative transfer corrections to investigate the relation between strong HI absorbers (N_HI >~ 10^17 /cm^2) and galaxies at redshift z = 3. We find a strong anti-correlation between the column density and the impact parameter that connects the absorber to the nearest galaxy. The median impact parameters for Lyman Limit (LL) and Damped Lyman-{\alpha} (DLA) systems are ~10 and ~1 proper kpc, respectively. If normalized to the size of the halo of the nearest central galaxy, the median impact parameters for LL and DLA systems become ~1 and ~10^-1 virial radii, respectively. At a given HI column density, the impact parameter increases with the mass of the closest galaxy, in agreement with observations. We predict most strong HI absorbers to be most closely associated with extremely low-mass galaxies, M_star < 10^8 M_sun and star formation rate <10^-1 M_sun/yr. We also find a correlation between the column density of absorbers and the mass of the nearest galaxy. This correlation is most pronounced for DLAs with N_HI > 10^21 /cm^2 which are typically close to galaxies with M_star >~ 10^9 M_sun. Similar correlations exist between column density and other properties of the associated galaxies such as their star formation rates, halo masses and HI content. The galaxies nearest to HI absorbers are typically far too faint to be detectable with current instrumentation, which is consistent with the high rate of (often unpublished) non-detections in observational searches for the galaxy counterparts of strong HI absorbers. Moreover, we predict that the detected nearby galaxies are typically not the galaxies that are most closely associated with the absorbers, thus causing the impact parameters, star formation rates and stellar masses of the observed counterparts to be biased high.

A substantial fraction of all ionizing photons originate from radiative recombinations. However, in radiative transfer calculations this recombination radiation is often assumed to be absorbed 'on-the-spot' because for most methods the computational cost associated with the inclusion of gas elements as sources is prohibitive. We present a new, CPU and memory efficient implementation for the transport of ionizing recombination radiation in the TRAPHIC radiative transfer scheme. TRAPHIC solves the radiative transfer equation by tracing photon packets at the speed of light and in a photon-conserving manner in spatially adaptive smoothed particle hydrodynamics simulations. Our new implementation uses existing features of the TRAPHIC scheme to add recombination radiation at no additional cost in the limit in which the fraction of the simulation box filled with radiation approaches 1. We test the implementation by simulating an HII region in photoionization equilibrium and comparing to reference solutions presented in the literature, finding excellent agreement. We apply our implementation to discuss the evolution of the HII region to equilibrium. We show that the widely used case A and B approximations yield accurate ionization profiles only near the source and near the ionization front, respectively. We also discuss the impact of recombination radiation on the geometry of shadows behind optically thick absorbers. We demonstrate that the shadow region may be completely ionized by the diffuse recombination radiation field and discuss the important role of heating by recombination radiation in the shadow region.

It is often assumed that local sources of ionizing radiation have little impact on the distribution of HI in the post-reionization Universe. While this is a good assumption for the IGM, analytic arguments suggest that local sources may typically be more important than the background radiation for high column density absorbers (N_HI > 10^17 /cm^2). We post-process cosmological simulations with accurate radiation transport to investigate the impact of local stellar sources on the HI distribution. We demonstrate that the limited numerical resolution and the simplified treatment of the ISM that are typical of current cosmological simulations provide significant challenges, but that many of the problems can be overcome by taking two steps. First, using ISM particles rather than stellar particles as sources results in a much better sampling of the source distribution. Second, by rescaling the source luminosities so that the amount of radiation escaping into the IGM agrees with that required to produce the observed background, many of the results become insensitive to errors in the predicted fraction of the radiation that escapes the immediate vicinity of the sources. By adopting this strategy and by varying the assumptions about the structure of the unresolved ISM, we conclude that we can robustly estimate the effect of local sources for column densities N_HI << 10^21 /cm^2. However, neither the escape fraction of ionizing radiation nor the effect of local sources on the abundance of N_HI >~ 10^21 systems can be predicted with confidence. We find that local stellar radiation is unimportant for N_HI << 10^17, but that it can affect Lyman Limit and DLA systems. For 10^18 < N_HI < 10^21 the impact of local sources increases with redshift. At z = 5 the abundance of absorbers with N_HI >> 10^17 is substantially reduced, but at z = 0 the reduction only becomes significant for N_HI >~ 10^21 /cm^2.

We use a set of cosmological simulations combined with radiative transfer calculations to investigate the distribution of neutral hydrogen in the post-reionization Universe. We assess the contributions from the metagalactic ionizing background, collisional ionization and diffuse recombination radiation to the total ionization rate at redshifts z=0-5. We find that the densities above which hydrogen self-shielding becomes important are consistent with analytic calculations and previous work. However, because of diffuse recombination radiation, whose intensity peaks at the same density, the transition between highly ionized and self-shielded regions is smoother than what is usually assumed. We provide fitting functions to the simulated photoionization rate as a function of density and show that post-processing simulations with the fitted rates yields results that are in excellent agreement with the original radiative transfer calculations. The predicted neutral hydrogen column density distributions agree very well with the observations. In particular, the simulations reproduce the remarkable lack of evolution in the column density distribution of Lyman limit and weak damped Ly\alpha\ systems below z = 3. The evolution of the low column density end is affected by the increasing importance of collisional ionization with decreasing redshift. On the other hand, the simulations predict the abundance of strong damped Ly\alpha\ systems to broadly track the cosmic star formation rate density.

We model the evolution of IR galaxies using a phenomenological approach to match the observed source counts at different IR wavelengths. We introduce a new algorithm for reproducing source counts based on direct integration of probability distributions rather than Monte-Carlo sampling. We construct a simple model for the evolution of the luminosity function and the colour distribution of IR galaxies which utilizes a minimum number of free parameters. Moreover we analyze how each of these parameters is constrained by observational data. The model is based on pure luminosity evolution and adopts the Dale & Helou SED templates. We find that the 850um source counts and their redshift distribution depend strongly on the shape of the luminosity evolution function, but only weakly on the details of the SEDs. We derive the best-fit evolutionary model using the 850um counts and redshift distribution as constraints. Moreover our best-fit shows a flattening of the faint end of the luminosity function towards high redshifts and requires a colour evolution which implies the typical dust temperatures of objects with the same luminosities to decrease with redshift. We compare our best-fit model to observed source counts at shorter and longer wavelengths which indicates our model reproduces the 70um and 1100um source counts remarkably well, but under-produces the counts at intermediate wavelengths. Analysis reveals that the discrepancy arises at low redshifts, indicating that revision of the adopted SED library towards lower dust temperatures (at a fixed infrared luminosity) is required. This modification is equivalent to a population of cold galaxies existing at low redshifts, as also indicated by recent Herschel results, which are underrepresented in IRAS sample. We show that the modified model successfully reproduces the source counts in a wide range of IR and submm wavelengths.

We use the weighted integral form of spherical Bessel functions, and introduce a new analytical set of complete and biorthogonal potential--density basis functions. The potential and density functions of the new set have finite central values and they fall off, respectively, similar to $r^{-(1+l)}$ and $r^{-(4+l)}$ at large radii where $l$ is the latitudinal quantum number of spherical harmonics. The lowest order term associated with $l=0$ is the perfect sphere of de Zeeuw. Our basis functions are intrinsically suitable for the modeling of three dimensional, soft-centred stellar systems and they complement the basis sets of Clutton-Brock, Hernquist & Ostriker and Zhao. We test the performance of our functions by expanding the density and potential profiles of some spherical and oblate galaxy models.

A set of bi-orthogonal potential-density basis functions is introduced to model the density and its associated gravitational field of three dimensional stellar systems. Radial components of our basis functions are weighted integral forms of spherical Bessel functions. We discuss the properties of our basis functions and demonstrate their shapes for the latitudinal Fourier number $l=2$.