Dark matter haloes play a fundamental role in cosmological structure formation. The most common approach to model their assembly mechanisms is through N-body simulations. In this work we present an innovative pathway to predict dark matter halo formation from the initial density field using a Deep Learning algorithm. We implement and train a Deep Convolutional Neural Network (DCNN) to solve the task of retrieving Lagrangian patches from which dark matter halos will condense. The volumetric multi-label classification task is turned into a regression problem by means of the euclidean distance transformation. The network is complemented by an adaptive version of the watershed algorithm to form the entire protohalo identification pipeline. We show that splitting the segmentation problem into two distinct sub-tasks allows for training smaller and faster networks, while the predictive power of the pipeline remains the same. The model is trained on synthetic data derived from a single full N-body simulation and achieves deviations of ~10% when reconstructing the dark matter halo mass function at z=0. This approach represents a promising framework for learning highly non-linear relations in the primordial density field. As a practical application, our method can be used to produce mock dark matter halo catalogues directly from the initial conditions of N-body simulations.

The ever-increasing sensitivity of the network of gravitational-wave detectors has resulted in the accelerated rate of detections from compact binary coalescence systems in the third observing run of Advanced LIGO and Advanced Virgo. Not only has the event rate increased, but also the distances to which phenomena can be detected, leading to a rise in the required sky volume coverage to search for counterparts. Additionally, the improvement of the detectors has resulted in the discovery of more compact binary mergers involving neutron stars, revitalizing dedicated follow-up campaigns. While significant effort has been made by the community to optimize single telescope observations, using both synoptic and galaxy-targeting methods, less effort has been paid to coordinated observations in a network. This is becoming crucial, as the advent of gravitational-wave astronomy has garnered interest around the globe, resulting in abundant networks of telescopes available to search for counterparts. In this paper, we extend some of the techniques developed for single telescopes to a telescope network. We describe simple modifications to these algorithms and demonstrate them on existing network examples. These algorithms are implemented in the open-source software \texttt{gwemopt}, used by some follow-up teams, for ease of use by the broader community.

One of the most computationally demanding aspects of the hydrodynamical modelling of Astrophysical phenomena is the transport of energy by radiation or relativistic particles. Physical processes involving energy transport are ubiquitous and of capital importance in many scenarios ranging from planet formation to cosmic structure evolution, including explosive events like core collapse supernova or gamma-ray bursts. Moreover, the ability to model and hence understand these processes has often been limited by the approximations and incompleteness in the treatment of radiation and relativistic particles. The DIAPHANE project has focused in developing a portable and scalable library that handles the transport of radiation and particles (in particular neutrinos) independently of the underlying hydrodynamic code. In this work, we present the computational framework and the functionalities of the first version of the DIAPHANE library, which has been successfully ported to three different smoothed-particle hydrodynamic codes, GADGET2, GASOLINE and SPHYNX. We also present validation of different modules solving the equations of radiation and neutrino transport using different numerical schemes.

The dearth of dwarf galaxies in the local universe is hard to reconcile with the large number of low mass haloes expected within the concordance $\Lambda$CDM paradigm. In this paper we perform a systematic evaluation of the uncertainties affecting the measurement of DM halo abundance using galaxy kinematics. Using a large sample of dwarf galaxies with spatially resolved kinematic data we derive a correction to obtain the observed abundance of galaxies as a function of their halo maximum circular velocity from the line-of-sight velocity function in the Local Volume. This estimate provides a direct means of comparing the predictions of theoretical models and simulations (including nonstandard cosmologies and novel galaxy formation physics) to the observational constraints. The new "galactic $V_{max}$" function is steeper than the line-of-sight velocity function but still shallower than the theoretical CDM expectation, showing that some unaccounted physical process is necessary to reduce the abundance of galaxies and/or drastically modify their density profiles compared to CDM haloes. Using this new galactic $V_{max}$ function, we investigate the viability of baryonic solutions such as feedback-powered outflows and photoevaporation of gas from an ionising radiation background. At the 3-$\sigma$ confidence level neither energetic feedback nor photoevaporation are effective enough to reconcile the disagreement. In the case of maximum baryonic effects, the theoretical estimate still deviates significantly from the observations for $V_{max} < 60$ km/s. CDM predicts at least 1.8 times more galaxies with $V_{max} = 50$ km/s and 2.5 times more than observed at $30$ km/s. Recent hydrodynamic simulations seem to resolve the discrepancy but disagree with the properties of observed galaxies with resolved kinematics. (abridged)

The observed number of dwarf galaxies as a function of rotation velocity is significantly smaller than predicted by the standard model of cosmology. This discrepancy cannot be simply solved by assuming strong baryonic feedback processes, since they would violate the observed relation between maximum circular velocity ($v_{\rm max}$) and baryon mass of galaxies. A speculative but tantalising possibility is that the mismatch between observation and theory points towards the existence of non-cold or non-collisionless dark matter (DM). In this paper, we investigate the effects of warm, mixed (i.e warm plus cold), and self-interacting DM scenarios on the abundance of dwarf galaxies and the relation between observed HI line-width and maximum circular velocity. Both effects have the potential to alleviate the apparent mismatch between the observed and theoretical abundance of galaxies as a function of $v_{\rm max}$. For the case of warm and mixed DM, we show that the discrepancy disappears, even for luke-warm models that evade stringent bounds from the Lyman-$\alpha$ forest. Self-interacting DM scenarios can also provide a solution as long as they lead to extended ($\gtrsim 1.5$ kpc) dark matter cores in the density profiles of dwarf galaxies. Only models with velocity-dependent cross sections can yield such cores without violating other observational constraints at larger scales.

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

We explore fundamental properties of the distribution of low mass dark matter halos within the cosmic web using warm dark matter (WDM) and cold dark matter (CDM) cosmological simulations. Using self abundance-matched mock galaxy catalogs, we show that the distribution of dwarf galaxies in a WDM universe, wherein low mass halo formation is heavily suppressed, is nearly indistinguishable to that of a CDM universe whose low mass halos are not seen because galaxy formation is suppressed below some threshold halo mass. However, if the scatter between dwarf galaxy luminosity and halo properties is large enough, low mass CDM halos would sometimes host relatively bright galaxies thereby populating CDM voids with the occasional isolated galaxy and reducing the numbers of completely empty voids. Otherwise, without high mass to light scatter, all mock galaxy clustering statistics that we consider--the auto-correlation function, the numbers and radial profiles of satellites, the numbers of isolated galaxies, and the PDF of small voids--are nearly identical in CDM and WDM. WDM voids are neither larger nor emptier than CDM voids, when constructed from abundance-matched halo catalogs. It is thus a challenge to determine whether the CDM problem of the over-abundance of small halos with respect to the number density of observed dwarf galaxies has a cosmological solution or an astrophysical solution. However, some clues about the dark matter particle and the scatter between the properties of dwarf galaxies and their dark matter halo hosts might be found in the cosmic web of galaxies in future surveys of the local volume.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

Future large-scale structure surveys of the Universe will aim to constrain the cosmological model and the true nature of dark energy with unprecedented accuracy. In order for these surveys to achieve their designed goals, they will require predictions for the nonlinear matter power spectrum to sub-percent accuracy. Through the use of a large ensemble of cosmological N-body simulations, we demonstrate that if we do not understand the uncertainties associated with simulating structure formation, i.e. knowledge of the `true' simulation parameters, and simply seek to marginalize over them, then the constraining power of such future surveys can be significantly reduced. However, for the parameters {n_s, h, Om_b, Om_m}, this effect can be largely mitigated by adding the information from a CMB experiment, like Planck. In contrast, for the amplitude of fluctuations sigma8 and the time-evolving equation of state of dark energy {w_0, w_a}, the mitigation is mild. On marginalizing over the simulation parameters, we find that the dark-energy figure of merit can be degraded by ~2. This is likely an optimistic assessment, since we do not take into account other important simulation parameters. A caveat is our assumption that the Hessian of the likelihood function does not vary significantly when moving from our adopted to the 'true' simulation parameter set. This paper therefore provides strong motivation for rigorous convergence testing of N-body codes to meet the future challenges of precision cosmology.

The nature of structure formation around the particle free streaming scale is still far from understood. Many attempts to simulate hot, warm, and cold dark matter cosmologies with a free streaming cutoff have been performed with cosmological particle-based simulations, but they all suffer from spurious structure formation at scales below their respective free streaming scales -- i.e. where the physics of halo formation is most affected by free streaming. We perform a series of high resolution numerical simulations of different WDM models, and develop an approximate method to subtract artificial structures in the measured halo mass function. The corrected measurements are then used to construct and calibrate an extended Press-Schechter (EPS) model with sharp-$k$ window function and adequate mass assignment. The EPS model gives accurate predictions for the low redshift halo mass function of CDM and WDM models, but it significantly under-predicts the halo abundance at high redshifts. By taking into account the ellipticity of the initial patches and connecting the characteristic filter scale to the smallest ellipsoidal axis, we are able to eliminate this inconsistency and obtain an accurate mass function over all redshifts and all dark matter particle masses covered by the simulations. As an additional application we use our model to predict the microhalo abundance of the standard neutralino-CDM scenario and we give the first quantitative prediction of the mass function over the full range of scales of CDM structure formation.

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of LCDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sufficient statistical precision, but lacked robust verification of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than ~1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may affect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.

We explore the cosmological halo-to-halo scatter of the distribution of mass within dark matter halos utilizing a well-resolved statistical sample of clusters from the cosmological Millennium simulation. We find that at any radius, the spherically-averaged dark matter density of a halo (corresponding to the "smooth-component") and its logarithmic slope are well-described by a Gaussian probability distribution. At small radii (within the scale radius), the density distribution is fully determined by the measured Gaussian distribution in halo concentrations. The variance in the radial distribution of mass in dark matter halos is important for the interpretation of direct and indirect dark matter detection efforts. The scatter in mass profiles imparts approximately a 25 percent cosmological uncertainty in the dark matter density at the Solar neighborhood and a factor of ~3 uncertainty in the expected Galactic dark matter annihilation flux. The aggregate effect of halo-to-halo profile scatter leads to a small (few percent) enhancement in dark matter annihilation background if the Gaussian concentration distribution holds for all halo masses versus a 10 percent enhancement under the assumption of a log-normal concentration distribution. The Gaussian nature of the cluster profile scatter implies that the technique of "stacking" halos to improve signal to noise should not suffer from bias.

We explore the clustering properties of high redshift dark matter halos, focusing on halos massive enough to host early generations of stars or galaxies at redshift 10 and greater. Halos are extracted from an array of dark matter simulations able to resolve down to the "mini-halo" mass scale at redshifts as high as 30, thus encompassing the expected full mass range of halos capable of hosting luminous objects and sources of reionization. Halo clustering on large-scales agrees with the Sheth, Mo & Tormen halo bias relation within all our simulations, greatly extending the regime where large-scale clustering is confirmed to be "universal" at the 10-20% level (which means, for example, that 3sigma halos of cluster mass at z=0 have the same large-scale bias with respect to the mass distribution as 3sigma halos of galaxy mass at z=10). However, on small-scales, the clustering of our massive halos (> ~10^9 Msun/h) at these high redshifts is stronger than expected from comparisons with small-scale halo clustering extrapolated from lower redshifts. This implies "non-universality" in the scale-dependence of halo clustering, at least for the commonly used parameterizations of the scale-dependence of bias that we consider. We provide a fit for the scale-dependence of bias in our results. This study provides a basis for using extraordinarily high redshift galaxies (redshift ~10) as a probe of cosmology and galaxy formation at its earliest stages. We show also that mass and halo kinematics are strongly affected by finite simulation volumes. This suggests the potential for adverse affects on gas dynamics in hydrodynamic simulations of limited volumes, such as is typical in simulations of the formation of the "first stars", though further study is warranted.

The dark matter halo mass function is a key repository of cosmological information over a wide range of mass scales, from individual galaxies to galaxy clusters. N-body simulations have established that the friends-of-friends (FOF) mass function has a universal form to a surprising level of accuracy (< 10%). The high-mass tail of the mass function is exponentially sensitive to the amplitude of the initial density perturbations, the mean matter density parameter, Omega_m, and to the dark energy controlled late-time evolution of the density field. Observed group and cluster masses, however, are usually stated in terms of a spherical overdensity (SO) mass which does not map simply to the FOF mass. Additionally, the widely used halo models of structure formation -- and halo occupancy distribution descriptions of galaxies within halos -- are often constructed exploiting the universal form of the FOF mass function. This again raises the question of whether FOF halos can be simply related to the notion of a spherical overdensity mass. By employing results from Monte Carlo realizations of ideal Navarro-Frenk-White (NFW) halos and N-body simulations, we study the relationship between the two definitions of halo mass. We find that the vast majority of halos (80-85%) in the mass-range 10^{12.5}-10^{15.5} M_sun/h indeed allow for an accurate mapping between the two definitions (~ 5%), but only if the halo concentrations are known. Nonisolated halos fall into two broad classes: those with complex substructure that are poor fits to NFW profiles and those ``bridged'' by the (isodensity-based)FOF algorithm. A closer investigation of the bridged halos reveals that the fraction of these halos and their satellite mass distribution is cosmology dependent. (abridged)

The majority of X-ray burst sources do not display a burst rate that increases with luminosity as expected, but this is seen in the two clocked bursters XB1323-619 and GS1826-24. We present a detailed investigation of these two sources which in the case of the first source, spans 18 years. Based on measurements of the burst rate, X-ray luminosity, the alpha-parameter and the two time constants generally present in the burst decays, we demonstrate the importance of the rp nuclear burning process. A detailed comparison with theory shows that although the burst rate in each source agrees well with the theoretical value, there is a difference of more than a factor of 5 in the burst rate at a given luminosity between the sources. We show that the main reason for this is that the two sources have substantially different emitting areas on the neutron star in non-burst emission, a factor often neglected. Variation of this area may explain the inverse relation of burst rate with luminosity in the majority of burst sources.

We construct mock galaxy catalogues to analyse clustering properties of a Lambda cold dark matter (LCDM) universe within a cosmological dark matter simulation of sufficient resolution to resolve structure down to the scale of dwarfs. We show that there is a strong age-clustering correlation for objects likely to host luminous galaxies, which includes the satellite halo (subhalo) population. Older mock galaxies are significantly more clustered in our catalog, which consists of satellite haloes as well as the central peaks of discrete haloes, selected solely by peak circular velocity. This age dependence is caused mainly by the age-clustering relation for discrete haloes, recently found by Gao et al., acting mostly on field members, combined with the tendency for older mock galaxies to lie within groups and clusters, where galaxy clustering is enhanced. Our results suggest that the clustering age dependence is manifested in real galaxies. At small scales (less than ~5 Mpc/h), the very simple assumption that galaxy colour depends solely on halo age is inconsistent with the strength of the observed clustering colour trends, where red galaxies become increasingly more clustered than blue galaxies toward smaller scales, suggesting that luminosity weighted galaxy ages do not closely trace the assembly epoch of their dark matter hosts. The age dependence is present but is much weaker for satellite haloes lying within groups and clusters than for the global population.

We use an array of high-resolution N-body simulations to determine the mass function of dark matter haloes at redshifts 10-30. We develop a new method for compensating for the effects of finite simulation volume that allows us to find an approximation to the true ``global'' mass function. By simulating a wide range of volumes at different mass resolution, we calculate the abundance of haloes of mass 10^{5-12} Msun/h. This enables us to predict accurately the abundance of the haloes that host the sources that reionize the universe. In particular, we focus on the small mass haloes (>~10^{5.5-6} Msun/h) likely to harbour population III stars where gas cools by molecular hydrogen emission, early galaxies in which baryons cool by atomic hydrogen emission at a virial temperature of ~10^4K (10^{7.5-8} Msun/h), and massive galaxies that may be observable at redshift ~10. When we combine our data with simulations that include high mass halos at low redshift, we find that the best fit to the halo mass function depends not only on linear overdensity, as is commonly assumed in analytic models, but also upon the slope of the linear power spectrum at the scale of the halo mass. The Press-Schechter model gives a poor fit to the halo mass function in the simulations at all epochs; the Sheth-Tormen model gives a better match, but still overpredicts the abundance of rare objects at all times by up to 50%. Finally, we consider the consequences of the recently released WMAP 3-year cosmological parameters. These lead to much less structure at high redshift, reducing the number of z=10 ``mini-haloes'' by more than a factor of two and the number of z=30 galaxy hosts by more than four orders of magnitude. Code to generate our best-fit halo mass function may be downloaded from http://icc.dur.ac.uk/Research/PublicDownloads/genmf_readme.html

We model gas cooling in high-resolution N-body simulations in order to investigate the formation of the first generation of stars. We follow a region of a LCDM universe especially selected to contain a rich cluster by the present day. The properties of the dark haloes that form in these sub-solar mass resolution simulations are presented in a companion paper by Gao et al. The first gas clouds able to cool by molecular hydrogen line emission collapse at extremely high redshift, z ~ 47, when the mass of the dark halo is 2.4 x 10^5 Msun/h. By z ~ 30, a substantial population of haloes are capable of undergoing molecular hydrogen cooling although their ability to form stars is dependent on the efficiency of feedback processes such as dissociating Lyman-Werner radiation. The mass of the main halo grows extremely rapidly and, by z ~ 36, its virial temperature has reached 10^4K, at which point gas cooling becomes dominated by more effective atomic line processes. By z ~ 30, a small ``group'' of such potential galaxies will have formed unless prevented from doing so by feedback processes. By this redshift, massive (>~ 100 Msun) population III stars are able to ionise gas well beyond their own host halo and neighbouring HII regions can percolate to form an ionized superbubble. Such patches would be too widely separated to contribute significantly to reionisation at this time. The large number density of early cooling haloes in the pre-reionised universe raises the exciting prospect that this ultra-early generation of stars may be observable as gamma-ray bursts or supernovae.

We use cosmological LCDM numerical simulations to model the evolution of the substructure population in sixteen dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10^11 to 10^14 Msun at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation: dn/dv=(1/8)(v_cmax/v_cmax,host)^-4. Halo to halo variance in the VDF is a factor of roughly two to four. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for r >~ 0.3r_vir, where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease toward the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias toward more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power law cosmologies of P ~ k^n, and demonstrate that a steeper spectral index, n, results in significantly less substructure.

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

We present results from analysis of the observation of the dipping, quasi-periodic bursting LMXB XB 1323-619 made with XMM-Newton in Jan., 2003. A number of absorption lines were discovered in the spectrum, notably at 6.70 and 6.98 keV which we identify with scattering by ions Fe XXV and Fe XXVI. Such features have been seen in other dipping sources, but their origin has not been understood. Curve of growth analysis provided a consistent solution assuming collisional ionization with kT=45 keV, close to the electron temperature we previously determined for the ADC in this source. We thus propose that the absorption lines in the dipping LMXB are produced in the ADC. Spectral evolution in dipping was well-described by the progressive covering model which we have previously shown to give very good explanations of many dipping sources. We discuss the proposal of Boirin et al. (2004b), that spectral evolution in the dipping LMXB may be explained by subjecting the continuum to a highly ionized absorber. This would require a decrease in X-ray intensity by a factor of ~3 in dipping at energies where photoelectric absorption is not effective (40-100 keV) which previous analysis of BeppoSAX data on several sources rules out, e.g. in XB 1323-619, any decrease was less than 10+/-10%. We find a remarkable linear increase in the rate of X-ray bursts with time over the 14-year period since 1989, and a systematic non-linear increase in source luminosity (L). The linear variation of burst rate with L shows that the burst rate is proportional to mass accretion rate, and if continued implies that the gap between bursts will become zero on January 11, 2008. In reality we expect the source to undergo a transition from X-ray bursting to X-ray flaring confirming that flaring is unstable nuclear burning.

We review evidence that the census of Milky Way satellites similar to those known may be incomplete at low latitude due to obscuration and in the outer halo due to a decreasing sensitivity to dwarf satellites with distance. We evaluate the possible impact that incompleteness has on comparisons with substructure models by estimating corrections to the known number of dwarfs using empirical and theoretical models. If we assume that the true distribution of Milky Way satellites is uniform with latitude, then we estimate a 33% incompleteness in the total number of dwarfs due to obscuration at low latitude. Similarly, if we suppose that the radial distribution of Milky Way satellites matches that of M31, or that of the oldest sub-halos or the most massive sub-halos in a simulation, we estimate a total number of Milky Way dwarfs ranging from 1 -- 3 times the known population. Although the true level of incompleteness is quite uncertain, the fact that our extrapolations yield average total numbers of MW dwarfs that are realistically 1.5 -- 2 times the known population, shows that incompleteness needs to be taken seriously when comparing to models of dwarf galaxy formation. Interestingly, the radial distribution of the oldest sub-halos in a Lambda+CDM simulation of a Milky Way-like galaxy possess a close match to the observed distribution of M31's satellites, which suggests that reionization may be an important factor controlling the observability of sub-halos. We also assess the prospects for a new SDSS search for Milky Way satellites to constrain the possible incompleteness in the outer halo.

We present galaxy circular velocity functions (GCVFs) for 34 low redshift (z < 0.15) clusters identified in the Sloan Digital Sky Survey (SDSS), for fifteen clusters drawn from dark matter simulations of hierarchical structure growth in a LambdaCDM cosmology, and for ~29,000 SDSS field galaxies. We find that the observed and simulated cluster GCVFs take the form of a power law. The cumulative GCVFs of the simulated clusters are very similar across a wide range of cluster masses, provided individual subhalo circular velocities are scaled by the circular velocities of the parent cluster. Once all sources of error are accounted for, the intrinsic scatter in the cumulative, scaled observed cluster GCVF is consistent with the simulations. The slope of the observed cluster GCVF is -2.5, independent of cluster velocity dispersion. The average slope of the simulated GCVFs is somewhat steeper, although formally consistent given the errors. Using our highest resolution cluster, we find that the effects of baryons on galaxy rotation curves is to flatten the simulated cluster GCVF into better agreement with observations. Finally, we find that the field GCVF deviates significantly from a power law, being flatter than the cluster GCVF at circular velocities less than 200 km/s, and steeper at circular velocities greater than 200 km/s.

We use a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 $h^{-1}$Mpc volume containing 80 million particles. Our low redshift results allow us to probe low $\sigma$ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the $10^{10}-10^{11} \msun$ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z $\simeq$ 6) QSOs, which are thought to be hosted by haloes at 5$\sigma$ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.

The optical expansion parallax of NGC 6543 has been detected and measured using two epochs of HST images separated by a time baseline of only three years. We have utilized three separate methods of deriving the angular expansion of bright fiducials, the results of which are in excellent agreement. We combine our angular expansion estimates with spectroscopically obtained expansion velocities to derive a distance to NGC 6543 of 1001$\pm$269 pc. The deduced kinematic age of the inner bright core of the nebula is 1039$\pm$259 years; however, the kinematic age of the polar caps that surround the core is larger - perhaps the result of deceleration or earlier mass ejection. The morphology and expansion patterns of NGC 6543 provide insight into a complex history of axisymmetric, interacting stellar mass ejections.