Multiscale cluster lens mass mapping – I. Strong lensing modelling     

E. Jullo  J.-P. Kneib 《Monthly notices of the Royal Astronomical Society》2009,395(3):1319-1332

We propose a novel technique to refine the modelling of galaxy cluster mass distribution using gravitational lensing. The idea is to combine the strengths of both 'parametric' and 'non-parametric' methods to improve the quality of the fit. We develop a multiscale model that allows sharper contrast in regions of higher density where the number of constraints is generally higher. Our model consists of (i) a multiscale grid of radial basis functions with physically motivated profiles and (ii) a list of galaxy-scale potentials at the location of the cluster member galaxies. This arrangement of potentials of different sizes allows us to reach a high resolution for the model with a minimum number of parameters. We apply our model to the well-studied cluster Abell 1689. We estimate the quality of our mass reconstruction with a Bayesian Monte Carlo Markov Chain sampler. For a selected subset of multiple images, we manage to halve the errors between the positions of predicted and observed images compared to previous studies. This is due to the flexibility of multiscale models at intermediate scale between cluster and galaxy scale. The software developed for this paper is part of the public lenstool package which can be found at http://www.oamp.fr/cosmology/lenstool .  相似文献   

Interpreting the transmission windows of distant quasars     

A. Maselli  A. Ferrara  S. Gallerani 《Monthly notices of the Royal Astronomical Society》2009,395(4):1925-1933

We propose the apparent shrinking criterion (ASC) to interpret the spatial extent, R _w, of transmitted flux windows in the absorption spectra of high- z quasars. The ASC can discriminate between the two regimes in which R _w corresponds either to the physical size,   R _{H  ii}   , of the quasar H  ii region or to the distance,   R ^max_w  , at which the transmitted flux drops to  =0.1  and a Gunn–Peterson (GP) trough appears. In the first case [H  ii region (HR) regime], one can determine the intergalactic medium mean H  i fraction,   x _H I  ; in the second [proximity region (PR) regime], the value of R _w allows one to measure the local photoionization rate and the local enhancement of the photoionization rate,  Γ_G  , due to nearby/intervening galaxies. The ASC has been tested against radiative transfer+smoothed particle hydrodynamics numerical simulations, and applied to a sample of 15 high-   z ( z > 5.8  ) quasar spectra. All sample quasars are found to be in the PR regime; hence, their observed spectral properties (inner flux profile, extent of transmission window) cannot reliably constrain the value of   x _{H  i}   . Four sample quasars show evidence for a local enhancement (up to 50 per cent) in the local photoionization rate possibly produced by a galaxy overdensity. We discuss the possible interpretations and uncertainties of this result.  相似文献   

Numerical simulations of hot halo gas in galaxy mergers     

Manodeep Sinha  Kelly Holley-Bockelmann 《Monthly notices of the Royal Astronomical Society》2009,397(1):190-207

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.  相似文献   

The impact of dust on the scaling properties of galaxy clusters     

Antonio C. da Silva  rea Catalano  Ludovic Montier  Etienne Pointecouteau  Joseph Lanoux  Martin Giard 《Monthly notices of the Royal Astronomical Society》2009,396(2):849-859

We investigate the effect of dust on the scaling properties of galaxy clusters based on hydrodynamic N -body simulations of structure formation. We have simulated five dust models plus radiative cooling and adiabatic models using the same initial conditions for all runs. The numerical implementation of dust was based on the analytical computations of Montier & Giard. We set up dust simulations to cover different combinations of dust parameters that make evident the effects of size and abundance of dust grains. Comparing our radiative plus dust cooling runs with a purely radiative cooling simulation, we find that dust has an impact on cluster scaling relations. It mainly affects the normalization of the scalings (and their evolution), whereas it introduces no significant differences in their slopes. The strength of the effect critically depends on the dust abundance and grain size parameters as well as on the cluster scaling. Indeed, cooling due to dust is effective in the cluster regime and has a stronger effect on the 'baryon driven' statistical properties of clusters such as   L _X– M , Y – M , S – M   scaling relations. Major differences, relative to the radiative cooling model, are as high as 25 per cent for the   L _X– M   normalization, and about 10 per cent for the Y – M and S – M normalizations at redshift zero. On the other hand, we find that dust has almost no impact on the 'dark matter driven'   T _mw– M   scaling relation. The effects are found to be dependent in equal parts on both dust abundances and grain size distributions for the scalings investigated in this paper. Higher dust abundances and smaller grain sizes cause larger departures from the radiative cooling (i.e. with no dust) model.  相似文献   

An implementation of radiative transfer in the cosmological simulation code gadget     

Margarita Petkova  Volker Springel 《Monthly notices of the Royal Astronomical Society》2009,396(3):1383-1403

We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH) simulation code gadget . It is based on a fast, robust and photon-conserving integration scheme where the radiation transport problem is approximated in terms of moments of the transfer equation and by using a variable Eddington tensor as a closure relation, following the Optically Thin Variable Eddington Tensor suggestion of Gnedin & Abel. We derive a suitable anisotropic diffusion operator for use in the SPH discretization of the local photon transport, and we combine this with an implicit solver that guarantees robustness and photon conservation. This entails a matrix inversion problem of a huge, sparsely populated matrix that is distributed in memory in our parallel code. We solve this task iteratively with a conjugate gradient scheme. Finally, to model photon sink processes we consider ionization and recombination processes of hydrogen, which is represented with a chemical network that is evolved with an implicit time integration scheme. We present several tests of our implementation, including single and multiple sources in static uniform density fields with and without temperature evolution, shadowing by a dense clump and multiple sources in a static cosmological density field. All tests agree quite well with analytical computations or with predictions from other radiative transfer codes, except for shadowing. However, unlike most other radiative transfer codes presently in use for studying re-ionization, our new method can be used on-the-fly during dynamical cosmological simulation, allowing simultaneous treatments of galaxy formation and the re-ionization process of the Universe.  相似文献   

A counterpart to the radial-orbit instability in triaxial stellar systems     

F. Antonini  R. Capuzzo-Dolcetta  D. Merritt 《Monthly notices of the Royal Astronomical Society》2009,399(2):671-682

  相似文献   

The influence of magnetic fields on neutrino-dominated accretion disc     

Yi Xie  Zhao-Yan Huang  Xiang-Fu Jia  Shao-Juan Fan  Fang-Fang Liu 《Monthly notices of the Royal Astronomical Society》2009,398(2):583-590

We consider the influence of magnetic fields on the model of neutrino-dominated accretion flows (NDAFs) for gamma-ray bursts (GRBs) via the assumption that the accretion rate of the disc is totally caused by the torque of the Lorentz force, i.e. the magnetic braking of large-scale magnetic fields and magnetic viscosity of small-scale magnetic fields. We calculate the structure, composition, luminosity of neutrino emission and the Poynting flux, and the rate of mass loss driven by neutrino heating or launched centrifugally by large-scale magnetic fields, based on the physical condition of the magnetized NDAFs. It is shown that the magnetized disc is favourable to interpret the diverse prompt emissions as well as the X-ray flares observed in the early afterglow of GRBs.  相似文献   

An MHD gadget for cosmological simulations     

K. Dolag  F. Stasyszyn 《Monthly notices of the Royal Astronomical Society》2009,398(4):1678-1697

Various radio observations have shown that the hot atmospheres of galaxy clusters are magnetized. However, our understanding of the origin of these magnetic fields, their implications on structure formation and their interplay with the dynamics of the cluster atmosphere, especially in the centres of galaxy clusters, is still very limited. In preparation for the upcoming new generation of radio telescopes (like Expanded Very Large Array, Low Wavelength Array, Low Frequency Array and Square Kilometer Array), a huge effort is being made to learn more about cosmological magnetic fields from the observational perspective. Here we present the implementation of magnetohydrodynamics (MHD) in the cosmological smoothed particle hydrodynamics (SPH) code gadget . We discuss the details of the implementation and various schemes to suppress numerical instabilities as well as regularization schemes, in the context of cosmological simulations. The performance of the SPH–MHD code is demonstrated in various one- and two-dimensional test problems, which we performed with a fully, three-dimensional set-up to test the code under realistic circumstances. Comparing solutions obtained using athena , we find excellent agreement with our SPH–MHD implementation. Finally, we apply our SPH–MHD implementation to galaxy cluster formation within a large, cosmological box. Performing a resolution study we demonstrate the robustness of the predicted shape of the magnetic field profiles in galaxy clusters, which is in good agreement with previous studies.  相似文献   

Cosmological systematics beyond nuisance parameters: form-filling functions     

T. D. Kitching  A. Amara  F. B. Abdalla  B. Joachimi  A. Refregier 《Monthly notices of the Royal Astronomical Society》2009,399(4):2107-2128

In the absence of any compelling physical model, cosmological systematics are often misrepresented as statistical effects and the approach of marginalizing over extra nuisance systematic parameters is used to gauge the effect of the systematic. In this article, we argue that such an approach is risky at best since the key choice of function can have a large effect on the resultant cosmological errors.
As an alternative we present a functional form-filling technique in which an unknown, residual, systematic is treated as such. Since the underlying function is unknown, we evaluate the effect of every functional form allowed by the information available (either a hard boundary or some data). Using a simple toy model, we introduce the formalism of functional form filling. We show that parameter errors can be dramatically affected by the choice of function in the case of marginalizing over a systematic, but that in contrast the functional form-filling approach is independent of the choice of basis set.
We then apply the technique to cosmic shear shape measurement systematics and show that a shear calibration bias of  | m ( z )| ≲ 10⁻³ (1 + z )^0.7  is required for a future all-sky photometric survey to yield unbiased cosmological parameter constraints to per cent accuracy.
A module associated with the work in this paper is available through the open source icosmo code available at http://www.icosmo.org .  相似文献   

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination     

Xiao-Gang Pan  Hai-Yin Zhou 《Chinese Astronomy and Astrophysics》2009,33(3):293-304

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.  相似文献