Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

The formation of the first stars out of metal-free gas appears to result in stars at least an order of magnitude more massive than in the present-day case. We here consider what controls the transition from a primordial to a modern initial mass function. It has been proposed that this occurs when effective metal line cooling occurs at a metallicity threshold of Z/Z_sun > 10^{-3.5}. We study the influence of low levels of metal enrichment on the cooling and collapse of initially ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations with particle splitting. Our initial conditions represent protogalaxies forming within a previously ionized H ii region that has not yet had time to cool and recombine. These differ considerably from those used in simulations predicting a metallicity threshold, where the gas was initially cold and only partially ionized. In the centrally condensed potential that we study here, a wide variety of initial conditions for the gas yield a monolithic central collapse. Our models show no fragmentation during collapse to number densities as high as 10^5 cm^{-3}, for metallicities reaching as high as 10^{-1} Z_sun in one rotating case, far above the threshold suggested by previous work. Rotation allows for the formation of gravitationally stable gas disks over large fractions of the local Hubble time. Turbulence slows the growth of the central density slightly, but both spherically symmetric and turbulent initial conditions collapse and form a single sink particle. We therefore argue that fragmentation at moderate density depends on the initial conditions for star formation more than on the metal abundances present.

We present results of hydrodynamic simulations of star formation triggered by cloud-cloud collisions. During the early stages of star formation, low-mass objects form by gravitational instabilities in protostellar discs. A number of these low-mass objects are in the sub-stellar mass range, including a few objects of planetary mass. The disc instabilities that lead to the formation of low-mass objects in our simulations are the product of disc-disc interactions and/or interactions between the discs and their surrounding gas.

We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using an SPH code which satisfies the Jeans condition by invoking On-the-Fly Particle Splitting. Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor-Ebert, spheres. Collisions are characterised by M_0 (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius), and M (Mach Number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling, and (iii) adiabatic heating in optically thick protostellar fragments. The efficiency of star formation is found to vary between 10% and 30% in the different collisions studied and it appears to increase with decreasing M_0, and/or decreasing b, and/or increasing M. For b<0.5 collisions produce shock compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high, 1 to 5 x 10^{-5} M_sun yr^{-1}, for the first 1 to 3 x 10^4 yrs. The densities in the filaments, n >~ 5 x 10^5 cm^{-3}, are sufficient that they could be mapped in NH_3 or CS line radiation, in nearby star formation regions.

We combine bright XMM data with the Chandra Deep Field South observations in order to explore the behavior of the intrinsic AGN absorption, as a function of redshift and luminosity.Our sample consists of 359 sources selected in the hard 2-8 keV band, spanning the flux range 6\times10^{-16}-$3\times10^{-13} erg s^-1 cm^-2 with a high rate of spectroscopic or photometric redshift completeness (100 and 85 per cent respectively for the Chandra and XMM data. We derive the column density values using X-ray spectral fits. We find that the fraction of obscured AGN falls with increasing luminosity in agreement with previous findings. The fraction of obscured AGN shows an apparent increase at high redshifts (z>2). Simulations show that this effect can be most probably attributed to the fact that at high redshifts the column densities are overestimated.

In this paper we serendipitously identify X-ray cluster candidates using XMM-Newton archival observations complemented by 5-band optical photometric follow-up observations (r~23 mag) as part of the X-ray Identification (XID) programme. Our sample covers an area of ~2.1 sq. deg (15 XMM-Newton fields) and comprises a total of 21 (19 serendipitous + 2 target) extended X-ray sources to the limit f(0.5-2keV) ~ 6x10^{-15} cgs with a high probability (> 99.9%) of being extended on the XMM-Newton images. Of the 21 cluster candidates 7 are spectroscopically confirmed in the literature. Exploiting the optical data available for these fields we discover that 68% of the X-ray cluster candidates are associated with optical galaxy overdensities. We also attempt to constrain the redshifts of our cluster candidates using photometric methods. We thus construct the photometric redshift distribution of galaxies in the vicinity of each X-ray selected cluster candidate and search for statistically significant redshift peaks against that of the background distribution of field galaxies. Comparison of the photometric with spectroscopic redshift estimates for the confirmed clusters suggest that our simple method is robust out to z~0.5. For clusters at higher z, deeper optical data are required to estimate reliable photometric redshifts. Finally, using the sample of the 19 serendipitous X-ray selected cluster candidates we estimate their surface density down to f(0.5-2keV) ~ 6x10^{-15} cgs and find it to be in fair agreement with previous and recent studies.

(Abridged) In this paper we present photometric redshift estimates for a sample of X-ray selected sources detected in the wide field (~2 deg^2), bright [f_{X} (0.5-8 keV)~10^{-14} cgs] XMM-Newton/2dF survey. Unlike deeper X-ray samples comprising a large fraction of sources with colours dominated by the host galaxy, our bright survey primarily probes the QSO X-ray population. Therefore photometric redshift methods employing both galaxy and QSO templates need to be used. We employ the photometric redshift technique of Hatziminaoglou, Mathez & Pello (2000) using 5-band photometry from the SDSS. We separate our X-ray sources according to their optical profile to point-like and extended. We apply QSO and galaxy templates to the point-like and extended sources respectively. X-ray sources associated with Galactic stars are identified and discarded from our point-like sample on the basis of their low X-ray--to--optical flux ratio and their broad band colours that are best fit by stellar templates. Comparison of our results with spectroscopic redshifts available, allows calibration of our method and estimation of the photometric redshift accuracy. For ~70 per cent of the point-like sources photometric redshifts are correct within dz <= 0.3 (or ~75 per cent have dz/(1+z) <= 0.2), and the rms scatter is estimated to be sigma_z = 0.30. For the optically extended objects the photometric redshifts work only in the case of red (g - r > 0.5 mag) sources yielding dz <= 0.15 and dz/(1+z) <= 0.2 for 73 and 93 per cent respectively. However, we find that the above photometric redshift technique does not work in the case of extended sources with blue colours (g - r < 0.5): such sources cannot be fit successfully by QSO or galaxy templates, or any linear combination of the two.

We report the discovery of a new, non-eclipsing SU UMa-type dwarf nova, HS 2219+1824. Photometry obtained in quiescence (V ~ 17.5) reveals a double-humped light curve from which we derive an orbital period of ~ 86.2 min. Additional photometry obtained during a superoutburst reaching V ~ 12.0 clearly shows superhumps with a period of ~ 89.05 min. The optical spectrum contains double-peaked Balmer and HeI emission lines from the accretion disc as well as broad absorption troughs of Hbeta, Hgamma, and Hdelta from the white dwarf primary star. Modelling of the optical spectrum implies a white dwarf temperature of 13000 K <~ Twd <~ 17000 K, a distance of 180 pc <~ d <~ 230 pc, and suggests that the spectral type of the donor star is later than M5. Phase-resolved spectroscopy obtained during quiescence reveals a narrow Halpha emission line component which has a radial velocity amplitude and phase consistent with an origin on the secondary star, possibly on the irradiated hemisphere facing the white dwarf. This constitutes the first detection of line emission from the secondary star in a quiescent SU UMa star.

We present an analysis of the X-ray spectral properties of 61 hard X-ray selected (2-8 keV) sources from the bright ($f_{2-8 keV}>10^{-14}$ \funits) {\it XMM-Newton}/2dF survey. This comprises of 9 {\it XMM-Newton} pointings in the North Galactic Pole region ($\sim 1.6 \rm deg^2$) and overlaps with the SDSS, 2QZ and 2dFGRS surveys. Our sources contribute about 50 per cent of the 2-10 keV X-ray background down to the flux limit of $10^{-14}$ cgs. The hardness ratio distribution of the sample suggests a deficit of heavily absorbed sources. A spectral fit to the co-added total source spectrum yields a steep photon index, $\Gamma=1.83^{+0.04}_{-0.05}$. All but 8 sources have optical counterparts down to the SDSS photometric limit of $r\approx22.5$. Spectroscopic identifications exist for 34 sources. The vast majority are associated with Broad-Line (BL) AGN (24 sources) while only 7 present narrow or no emission lines. Five sources are probably associated with Galactic stars. Finally, for another 17 probable AGN we present photometric redshifts. The combined spectrum of the 24 spectroscopically identified BL AGN is steep ($\Gamma=2.02^{+0.04}_{-0.05}$), while that of the 7 AGN, which do not present broad lines is flatter with $\Gamma=1.64^{+0.11}_{-0.11}$. The spectrum of the 8 optically unidentified sources is flat with $\Gamma\approx 1.1$. Spectral fits to the individual BL AGN reveal large absorption (rest-frame column density $\rm>10^{22} cm^{-2}$) in only two cases. The individual spectra of the NL AGN present significant evidence for even a moderate absorption ($\rm 3\times 10^{21} cm^{-2}$) in only one case.

We report follow-up observations of a new white dwarf/red dwarf binary HS2237+8154, identified as a blue variable star from the Hamburg Quasar Survey. Ellipsoidal modulation observed in the R band as well as the Halpha radial velocity variations measured from time-resolved spectroscopy determine the orbital period to be Porb=178.10+-0.08min. The optical spectrum of HS2237+8154 is well described by a combination of a Teff=11500+-1500K white dwarf (assuming log g=8) and a dM3.5+-0.5 secondary star. The distance implied from the flux scaling factors of both stellar components is d=105+-25pc. Combining the constraints obtained from the radial velocity of the secondary and from the ellipsoidal modulation, we derive a binary inclination of i~50deg-70deg and stellar masses of Mwd=0.47-0.67Msun and Msec=0.2-0.4Msun. All observations imply that the secondary star must be nearly Roche-lobe filling. Consequently, HS2237+8154 may be either a pre-cataclysmic variable close to the start of mass transfer, or (considering its orbital period) a cataclysmic variable that terminated mass transfer and entered the period gap, or a hibernating nova.

The technique of estimating redshifts using photometric rather than spectroscopic observations has recently received great attention due to its simplicity and the accuracy of the results obtained. In this work, we estimate photometric redshifts for an X-ray selected QSO sample. This is the first time this technique is applied on such a sample. We first calculate the accuracy of the results obtained by comparing photometric to spectroscopic redshifts for a sub-sample of our QSO sample: for the majority (~67%) of the objects in this sub-sample, photometric redshift estimates are correct within Dz<0.3. We then derive the photometric redshift distribution for the whole QSO sample. In the future, we expect to use the photometric redshift distribution in order to derive the distributions of properties such as the Hardness Ratio and hence the hydrogen column density, the luminosity function etc. As an example, we estimate here the dependence of the Hardness Ratio of the QSO sample on photometric redshift.

We have invented a new algorithm to use with self-gravitating SPH Star Formation codes. The new method is designed to enable SPH simulations to self-regulate their numerical resolution, i.e. the number of SPH particles; the latter is calculated using the Jeans condition (Bate & Burkert 1997) and the local hydrodynamic conditions of the gas. We apply our SPH with Particle Splitting code to cloud-cloud collision simulations. Chapter 2 lists the properties of our standard SPH code. Chapter 3 discusses the efficiency of the standard code as this is applied to simulations of rotating, uniform clouds with m=2 density perturbations. Chapter 4 [astro-ph/0203057] describes the new method and the tests that it has successfully been applied to. It also contains the results of the application of Particle Splitting to the case of rotating clouds as those of Chapter 3, where, with great computational efficiency, we have reproduced the results of FD codes and SPH simulations with large numbers of particles. Chapter 5 gives a detailed account of the cloud-cloud collisions studied, starting from a variety of initial conditions produced by altering the cloud mass, cloud velocity and the collision impact parameter. In the majority of the cases studied, the collisions produced filaments (similar to those observed in ammonia in nearby Star Forming Regions) or networks of filaments; groups of protostellar cores have been produced by fragmentation of the filaments. The accretion rates at these cores are comparable to those of Class 0 objects. Due to time-step constraints the simulations stop early in their evolution. The star formation efficiency of this mechanism is extrapolated in time and is found to be 10-20%.

We describe and demonstrate a method for increasing the resolution locally in a Smoothed Particle Hydrodynamic (SPH) simulation, by splitting particles. We show that in simulations of self-gravitating collapse (of the sort which are presumed to occur in star formation) the method is stable, and affords great savings in computer time and memory. When applied to the standard Boss & Bodenheimer test -- which has been shown to depend critically on fulfilment of the Jeans Condition -- the results are comparable both with those obtained using Adaptive Mesh Refinement, and with those obtained using a standard high-resolution SPH simulation, but they are achieved with considerably less computational resource. Further development and testing is required before the method can safely be applied to more general flows.

We analyse and compare the clustering of young stars in Chamaeleon I and Taurus. We compute the mean surface-density of companion stars \bar{N} as a function of angular displacement \theta from each star. We then fit \bar{N}(\theta) with two simultaneous power laws. For Chamaeleon I, the exponents of the power laws are 1.97 and 0.28, with the elbow at ~0.011 degrees. For Taurus, we obtain 2.02 and 0.87, with the elbow at ~0.013 deg. For both star clusters the observational data make quite large systematic excursions from the best fitting curve in the binary regime. These excursions may be attributable to evolutionary effects of the types discussed recently by Nakajima et al. and Bate et al. In the clustering regime the data conform to the best fitting curve very well. We also calculate the box-dimensions for the two star clusters. However, the limited dynamic range makes these estimates simply descriptors of the large-scale clustering, and not admissible evidence for fractality. We also propose two algorithms for objectively generating maps of constant stellar surface-density in young star clusters. Such maps are useful for comparison with molecular-line and dust-continuum maps of star-forming clouds, and with the results of numerical simulations. These maps retain information which is suppressed in the evaluation of \bar{N}(\theta). Analyzing such maps can be used to discriminate between fractal structure and single-level clustering, and to determine the degree of central condensation in small-N clusters. Algorithm I uses a universal smoothing length, and therefore has a restricted dynamic range, but it is implicitly normalized. Algorithm II uses a local smoothing length, which gives it much greater dynamic range, but it has to be normalized explicitly.