首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 484 毫秒
1.
We present a new Monte Carlo algorithm to generate merger trees describing the formation history of dark matter haloes. The algorithm is a modification of the algorithm of Cole et al. used in the galform semi-analytic galaxy formation model. As such, it is based on the Extended Press–Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the Λ cold dark matter Millennium N -body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases, we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees .  相似文献   

2.
Within the context of constraining an expansion of the dark energy equation of state   w ( z ),  we show that the eigendecomposition of Fisher matrices is sensitive to both the maximum order of the expansion and the basis set choice. We investigate the Fisher matrix formalism in the case that a particular function is expanded in some basis set. As an example we show results for an all-sky weak lensing tomographic experiment. We show that the set of eigenfunctions is not unique and that the best constrained functions are only reproduced accurately at very higher order   N ≳ 100  , a top-hat basis set requires an even higher order. We show that the common approach used for finding the marginalized eigenfunction errors is sensitive to the choice of  non- w ( z )  parameters and priors. The eigendecomposition of Fisher matrices is a potentially useful tool that can be used to determine the predicted accuracy with which an experiment could constrain   w ( z )  . It also allows for the reconstruction of the redshift sensitivity of the experiment to changes in   w ( z )  . However, the technique is sensitive to both the order and the basis set choice. Publicly available code is available as part of icosmo at http://www.icosmo.org .  相似文献   

3.
We present a set of gas-phase Planck mean and Rosseland mean opacity tables applicable for simulations of star and planet formation, stellar evolution and disc modelling at various metallicities in hydrogen-rich environments. The tables are calculated for gas temperatures between 1000 and 10 000 K and total hydrogen number densities between 102 and  1017 cm−3  . The carbon-to-oxygen ratio is varied from 0.43 to well above 2.0, the nitrogen-to-oxygen ratio between 0.14 and 100.0. The tables are calculated for a range of metallicities down to  [M/H]'= log  N M/ N H=−7.0  . We demonstrate how the mean opacities and the abundances of the opacity species vary with C/O, N/O and  [M/H]'  . We use the element abundances from Grevesse et al., and we provide additional tables for the oxygen abundance value from Caffau et al. All tables will be available online at http://star-www.st-and.ac.uk/~ch80/datasources.html .  相似文献   

4.
We present an online catalog for dwarfs with spectral types M, L and T. The catalog currently contains 893 objects listed in order of right ascension. For most of them the spectra can be downloaded from the catalog in ASCII or FITS formats. The web address is http://www.iac.es/galeria/ege/catalogo spectral/.  相似文献   

5.
We present a library of 1654 high-resolution stellar spectra, with a sampling of 0.3 Å and covering the wavelength range from 3000 to 7000 Å. The library was computed with the latest improvements in stellar atmospheres, incorporating non-local thermodynamic equilibrium (non-LTE) line-blanketed models for hot, massive  ( T eff≥ 27 500 K)  and line-blanketed models for cool  (3000 ≤ T eff≤ 4500 K)  stars. The total coverage of the grid is  3000 K ≤ T eff≤ 55 000 K  and  −0.5 ≤ log g ≤ 5.5  , for four chemical abundance values: twice solar, solar, half solar and 1/10 solar. Evolutionary synthesis models using this library are presented in a companion paper. We tested the general behaviour of the library by calculating and comparing equivalent widths of numerous H and He  i lines, and some of the commonly used metallic indices. We also compared the library with the empirical libraries STELIB and Indo-US. The full set of the synthetic stellar spectra is available from our websites ( http://www.iaa.csic.es/~rosa and http://www.astro.iag.usp.br/~lucimara/library.htm ).  相似文献   

6.
The He, C, N, and O abundances in more than 120 planetary nebulae (PNe) of our Galaxy and the Magellanic Clouds have been redetermined by analyzing new PNe observations. The characteristics of PNe obtained by modeling their spectra have been used to compile a new catalog of parameters for Galactic and extragalactic PNe, which is accessible at http://www.astro.spbu.ru/staff/afk/GalChemEvol.html. The errors in the parameters of PNe and their elemental abundances related to inaccuracies in the observational data have been analyzed. The He abundance is determined with an accuracy of 0.06 dex, while the errors in the C, N, and O abundances are 0.1–0.2 dex. Taking into account the inaccuracies in the corrections for the ionization stages of the elements whose lines are absent in the PNe spectra increases the errors in the He abundance to 0.1 dex and in the C, N, and O abundances to 0.2–0.3 dex. The elemental abundances in PNe of various Galactic subsystems and the Magellanic Clouds have been analyzed. This analysis suggests that the Galactic bulge objects are similar to type II PNe in Peimbert’s classification, whose progenitor stars belong to the thin-disk population with ages of at least 4–6 Gyr. A similarity between the elemental abundances in PNe of the Magellanic Clouds and the Galactic halo has been established.  相似文献   

7.
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 .  相似文献   

8.
As galaxy surveys become larger and more complex, keeping track of the completeness, magnitude limit and other survey parameters as a function of direction on the sky becomes an increasingly challenging computational task. For example, typical angular masks of the Sloan Digital Sky Survey contain about   N = 300 000  distinct spherical polygons. Managing masks with such large numbers of polygons becomes intractably slow, particularly for tasks that run in time     with a naive algorithm, such as finding which polygons overlap each other. Here we present a 'divide-and-conquer' solution to this challenge: we first split the angular mask into pre-defined regions called 'pixels', such that each polygon is in only one pixel, and then perform further computations, such as checking for overlap, on the polygons within each pixel separately. This reduces     tasks to     , and also reduces the important task of determining in which polygon(s) a point on the sky lies from     to     , resulting in significant computational speedup. Additionally, we present a method to efficiently convert any angular mask to and from the popular healpix format. This method can be generically applied to convert to and from any desired spherical pixelization. We have implemented these techniques in a new version of the mangle software package, which is freely available at http://space.mit.edu/home/tegmark/mangle/ , along with complete documentation and example applications. These new methods should prove quite useful to the astronomical community, and since mangle is a generic tool for managing angular masks on a sphere, it has the potential to benefit terrestrial mapmaking applications as well.  相似文献   

9.
  相似文献   

10.
This paper is primarily an investigation of whether the 'optimal extraction' techniques used in CCD spectroscopy can be applied to imaging photometry. It is found that using such techniques provides a gain of around 10 per cent in signal-to-noise ratio over normal aperture photometry. Formally, it is shown to be equivalent to profile fitting, but offers advantages of robust error estimation, freedom from bias introduced by mis-estimating the point spread function, and convenience. In addition some other techniques are presented, which can be applied to profile fitting, aperture photometry and the 'optimal' photometry. Code implementing these algorithms is available at http://www.astro.keele.ac.uk/~timn/.  相似文献   

11.
Many papers are devoted to the prediction of radiation conditions on board of a spacecraft (Pichkhadze et al., 2004; Khamidullina et al., 2008; 2012), and a number of software systems for corresponding calculations have been developed: the US information system CREME96 (https://creme.isde.vander-bilt.edu/); European SPENVIS (http://www.spenvis.oma.be/intro.php); Russian SEREIS (Kuznetsov et al., 2001; Model’ kosmosa, 2007) and COSRAD (http://cosrad.sinp.msu.ru/manual.html; Kuznetsov et al., 2011) based on the models of the radiation environment in near-Earth space (Bashkirov et al., 1998; Nymmik, 2004; Model’ kosmosa, 2007; Kuznetsov et al., 2011). In this paper we propose a simple calculation algorithm of short-term (for a few days) forecasting of dynamics of the radiation dose on the International Space Station (ISS) in radiation environment undisturbed by solar proton events. This algorithm does not use radiation environment models and detailed ballistic calculations, while it uses data of the onboard radiation monitoring system (RMS) and empirical relations, obtained for ISS orbital motion.  相似文献   

12.
We present the single stellar population (SSP) synthesis results of our new synthetic stellar atmosphere models library with a spectral sampling of 0.3 Å, covering the wavelength range from 3000 to 7000 Å for a wide range of metallicities (twice solar, solar, half solar and 1/10 solar). The stellar library is composed of 1650 spectra computed with the latest improvements in stellar atmospheres. In particular, it incorporates non-local thermodynamic equilibrium (LTE) line-blanketed models for hot  ( T eff≥ 27 500 K)  , and LTE line-blanketed models (Phoenix) for cool  (3000 ≤ T eff≤ 4500 K)  stars. Because of the high spectral resolution of this library, evolutionary synthesis models can be used to predict the strength of numerous weak absorption lines and the evolution of the profiles of the strongest lines over a wide range of ages. The SSP results have been calculated for ages from 1 Myr to 17 Gyr using the stellar evolutionary tracks provided by the Geneva and Padova groups. For young stellar populations, our results have a very detailed coverage of high-temperature stars with similar results for the Padova and Geneva isochrones. For intermediate and old stellar populations, our results, once degraded to a lower resolution, are similar to the ones obtained by other groups (limitations imposed by the stellar evolutionary physics notwidthstanding). The limitations and advantages of our models for the analysis of integrated populations are described. The full set of the stellar library and the evolutionary models are available for retrieval at the websites http://www.iaa.csic.es/~rosa and http://www.iaa.csic.es/~mcs/sed@ , or on request from the first two authors.  相似文献   

13.
We analyse the redshift space topology and geometry of the nearby Universe by computing the Minkowski functionals of the Updated Zwicky Catalogue (UZC). The UZC contains the redshifts of almost 20 000 galaxies, is 96 per cent complete to the limiting magnitude m Zw=15.5, and includes the Center for Astrophysics (CfA) Redshift Survey (CfA2). From the UZC we can extract volume-limited samples reaching a depth of 70  h −1 Mpc before sparse sampling dominates. We quantify the shape of the large-scale galaxy distribution by deriving measures of planarity and filamentarity from the Minkowski functionals. The nearby Universe shows a large degree of planarity and a small degree of filamentarity. This quantifies the sheet-like structure of the Great Wall, which dominates the northern region (CfA2N) of the UZC. We compare these results with redshift space mock catalogues constructed from high-resolution N -body simulations of two cold dark matter (CDM) models with either a decaying massive neutrino ( τ CDM) or a non-zero cosmological constant (ΛCDM). We use semi-analytic modelling to form and evolve galaxies in these dark matter‐only simulations. We are thus able, for the first time, to compile redshift space mock catalogues which contain galaxies, along with their observable properties, rather than dark matter particles alone. In both models the large-scale galaxy distribution is less coherent than the observed distribution, especially with regard to the large degree of planarity of the real survey. However, given the small volume of the region studied, this disagreement can still be a result of cosmic variance, as shown by the agreement between the ΛCDM model and the southern region of CfA2.  相似文献   

14.
  相似文献   

15.
N -body simulations of cold dark matter (CDM) have shown that, in this hierarchical structure formation model, dark matter halo properties, such as the density profile, the phase-space density profile, the distribution of axial ratio, the distribution of spin parameter and the distribution of internal specific angular momentum, follow 'universal' laws or distributions. Here, we study the properties of the first generation of haloes in a hot dark matter (HDM) dominated universe, as an example of halo formation through monolithic collapse. We find all these universalities to be present in this case also. Halo density profiles are very well fit by the Navarro, Frenk & White profile over two orders of magnitude in mass. The concentration parameter depends on mass as   c ∝ M 0.2  , reversing the dependence found in a hierarchical CDM universe. However, the concentration–formation time relation is similar in the two cases: earlier forming haloes tend to be more concentrated than their later forming counterparts. Halo formation histories are also characterized by two phases in the HDM case: an early phase of rapid accretion followed by slower growth. Furthermore, there is no significant difference between the HDM and CDM cases concerning the statistics of other halo properties: the phase-space density profile; the velocity anisotropy profile; the distribution of shape parameters; the distribution of spin parameter and the distribution of internal specific angular momentum are all similar in the two cases. Only substructure content differs dramatically. These results indicate that mergers do not play a pivotal role in establishing the universalities, thus contradicting models which explain them as consequences of mergers.  相似文献   

16.
We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f 0. Our method uses importance sampling to find the sampling DF f s that minimizes the mean-square formal errors in a given set of projections of the DF f 0. The result is a multimass N -body realization of the galaxy model in which 'interesting' regions of phase space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher mass particles.
As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N -body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor of ∼100 reduction in the variance at the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N -body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N -body simulations closer to the collisionless ideal.  相似文献   

17.
We use an extremely large volume  (2.4  h −3 Gpc3)  , high-resolution N -body simulation to measure the higher order clustering of dark matter haloes as a function of mass and internal structure. As a result of the large simulation volume and the use of a novel 'cross-moment' counts-in-cells technique which suppresses discreteness noise, we are able to measure the clustering of haloes corresponding to rarer peaks than was possible in previous studies; the rarest haloes for which we measure the variance are 100 times more clustered than the dark matter. We are able to extract, for the first time, halo bias parameters from linear up to fourth order. For all orders measured, we find that the bias parameters are a strong function of mass for haloes more massive than the characteristic mass   M *  . Currently, no theoretical model is able to reproduce this mass dependence closely. We find that the bias parameters also depend on the internal structure of the halo up to fourth order. For haloes more massive than   M *  , we find that the more concentrated haloes are more weakly clustered than the less concentrated ones. We see no dependence of clustering on concentration for haloes with masses   M < M *  ; this is contrary to the trend reported in the literature when segregating haloes by their formation time. Our results are insensitive to whether haloes are labelled by the total mass returned by the friends-of-friends group finder or by the mass of the most massive substructure. This implies that our conclusions are not an artefact of the particular choice of group finding algorithm. Our results will provide important input to theoretical models of galaxy clustering.  相似文献   

18.
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−3 (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 .  相似文献   

19.
We describe a made-to-measure (M2M) algorithm for constructing N -particle models of stellar systems from observational data (χ2M2M), extending earlier ideas by Syer & Tremaine. The algorithm properly accounts for observational errors, is flexible, and can be applied to various systems and geometries. We implement this algorithm in a parallel code nmagic and carry out a sequence of tests to illustrate its power and performance. (i) We reconstruct an isotropic Hernquist model from density moments and projected kinematics and recover the correct differential energy distribution and intrinsic kinematics. (ii) We build a self-consistent oblate three-integral maximum rotator model and compare how the distribution function is recovered from integral field and slit kinematic data. (iii) We create a non-rotating and a figure rotating triaxial stellar particle model, reproduce the projected kinematics of the figure rotating system by a non-rotating system of the same intrinsic shape, and illustrate the signature of pattern rotation in this model. From these tests, we comment on the dependence of the results from χ2M2M on the initial model, the geometry, and the amount of available data.  相似文献   

20.
The excitation mechanism of solar five-minute oscillations is studied in the present paper. We calculated the non-adiabatic oscillations of low- and intermediate-degree (l = 1 − 25) g4-p39 modes for the Sun. Both the thermodynamic and dynamic couplings are taken into account by using our non-local and time-dependent theory of convection. The results show that all the lowfrequencyf- and p-modes with periods P > 5.4 min are pulsationally unstable, while the coupling between convection and oscillations is neglected. However, when the convection coupling is taken into account, all the g- and low-frequency f- and p-modes with periods longer than 16 minutes (except the low-degree p1-modes) and the high frequency p-modes with periods shorter than 3 minutes become stable, and the intermediate-frequency p-modes with period from 3 to 16 minutes are pulsationally unstable. The pulsation amplitude growth rates depend only on the frequency and almost do not depend on l. They achieve the maximum at ν 3700 μHz (or P 270 sec). The coupling between convection and oscillations plays a key role for stabilization of low-frequency f- and p-modes and excitation of intermediate-frequency p-modes. We propose that the solar 5-minute oscillations are not caused by any single excitation mechanism, but they are resulted from the combined effect of “regular” coupling between convection and oscillations and turbulent stochastic excitation. For low- and intermediatefrequency p-modes, the coupling between convection and oscillations dominates; while for high-frequency modes, stochastic excitation dominates.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号