Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBreaking Beta: A comparison of mass modelling methods for spherical systems
65
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We apply four different mass modelling methods to a suite of publicly available mock data for spherical stellar systems. We focus on the recovery of the density and velocity anisotropy as a function of radius, using either line-of-sight velocity data only, or adding proper motion data. All methods perform well on isotropic and tangentially anisotropic mock data, recovering the density and velocity anisotropy within their 95% confidence intervals over the radial range 0.25 < R/Rhalf < 4, where Rhalf is the half light radius. However, radially-anisotropic mocks are more challenging. For line-of-sight data alone, only methods that use information about the shape of the velocity distribution function are able to break the degeneracy between the density profile and the velocity anisotropy to obtain an unbiased estimate of both. This shape information can be obtained through directly fitting a global phase space distribution function, by using higher order Virial Shape Parameters, or by assuming a Gaussian velocity distribution function locally, but projecting it self-consistently along the line of sight. Including proper motion data yields further improvements, and in this case, all methods give a good recovery of both the radial density and velocity anisotropy profiles.
rate research
Read More
Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and thoroughly tested. MAMPOSSt assumes a shape for the gravitational potential, but instead of postulating a shape for the distribution function in terms of energy and angular momentum, or supposing Gaussian line-of-sight velocity distributions, MAMPOSSt assumes a VA profile and a shape for the 3D velocity distribution, here Gaussian. MAMPOSSt requires no binning, differentiation, nor extrapolation of the observables. Tests on cluster-mass haloes from LambdaCDM cosmological simulations show that, with 500 tracers, MAMPOSSt is able to jointly recover the virial radius, tracer scale radius, dark matter scale radius and outer or constant VA with small bias (<10% on scale radii and <2% on the two other quantities) and inefficiencies of 10%, 27%, 48% and 20%, respectively. MAMPOSSt does not perform better when some parameters are frozen, and even worse when the virial radius is set to its true value, which appears to be the consequence of halo triaxiality. The accuracy of MAMPOSSt depends weakly on the adopted interloper removal scheme, including an efficient iterative Bayesian scheme that we introduce here, which can directly obtain the virial radius with as good precision as MAMPOSSt. Our tests show that MAMPOSSt with Gaussian 3D velocities is very competitive with, and up to 1000x faster than other methods. Hence, MAMPOSSt is a very powerful and rapid tool for the mass and anisotropy modeling of systems such as clusters and groups of galaxies, elliptical and dwarf spheroidal galaxies.
Sub-millimetre dust emission provides information on the physics of interstellar clouds and dust. Noise can produce anticorrelation between the colour temperature T_C and the spectral index beta. This must be separated from the intrinsic beta(T) relation of dust. We compare methods for the analysis of the beta(T) relation. We examine sub-millimetre observations simulated as simple modified black body emission or using 3D radiative transfer modelling. In addition to chi^2 fitting, we examine the results of the SIMEX method, basic Bayesian model, hierarchical models, and one method that explicitly assumes a functional form for beta(T). All methods exhibit some bias. Bayesian method shows significantly lower bias than direct chi^2 fits. The same is true for hierarchical models that also result in a smaller scatter in the temperature and spectral index values. However, significant bias was observed in cases with high noise levels. Beta and T estimates of the hierarchical model are biased towards the relation determined by the data with the highest S/N ratio. This can alter the recovered beta(T) function. With the method where we explicitly assume a functional form for the beta(T) relation, the bias is similar to the Bayesian method. In the case of an actual Herschel field, all methods agree, showing some degree of anticorrelation between T and beta. The Bayesian method and the hierarchical models can both reduce the noise-induced parameter correlations. However, all methods can exhibit non-negligible bias. This is particularly true for hierarchical models and observations of varying signal-to-noise ratios and must be taken into account when interpreting the results.
10,000 simulations of 1000-particle realisations of the same cluster are computed by direct force summation. Over three crossing times the original Poisson noise is amplified more than tenfold by self-gravity. The clusters fundamental dipole mode is strongly excited by Poisson noise, and this mode makes a major contribution to driving diffusion of stars in energy. The diffusive flow through action space is computed for the simulations and compared with the predictions of both local-scattering theory and the Balescu-Lenard (BL) equation. The predictions of local-scattering theory are qualitatively wrong because the latter neglects self-gravity. These results imply that local-scattering theory is of little value. Future work on cluster evolution should employ either N-body simulation or the BL equation. However, significant code development will be required to make use of the BL equation practicable.
We present high resolution H{sc i} 21cm Giant Meterwave Radio Telescope (GMRT) observations of the superthin galaxy FGC1540 with a spatial resolution of 10$$ $times$ 8$$ and a spectral resolution of 1.73 kms$^{-1}$ and an rms noise of 0.9 mJy per beam. We obtain its rotation curve as well as deprojected radial H{sc i} surface density profile by fitting a 3-dimensional tilted ring model directly to the H{sc i} data cubes by using the publicly-available software, Fully Automated Tirrific (FAT). We also present the rotation curve of FGC1540 derived from its optical spectroscopy study using the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. We use the rotation curve, the H{sc i} surface density profile together with Spitzer 3.6 $mu$m and the SDSS $i$--band data to construct the mass models for FGC1540. We find that both the Pseudo-isothermal (PIS), as well as Navarro-Frenk-White (NFW) dark matter (DM) halos, fit the observed rotation curve equally well. The PIS model indicates a compact dark matter halo ($R_{rm C}/R_{rm D}$ < 2), with the best-fitting core radius ($R_{rm C}$) approximately half the exponential stellar disc scale length ($R_{rm D}$), which is in agreement with the mass models of superthin galaxies studied earlier in the literature. Since the vertical thickness of the galactic stellar disc is determined by a balance between the net gravitational field and the velocity dispersion in the vertical direction, the compact dark matter halo may be primarily responsible in regulating the superthin vertical structure of the stellar disc in FGC1540 as was found in case of the superthin galaxy UGC7321.
Angular momentum, or spin, is a fundamental property of black holes (BHs), yet it is much more difficult to estimate than mass or accretion rate (for actively accreting systems). In recent years, high-quality X-ray observations have allowed for detailed measurements of the Fe K$alpha$ emission line, where relativistic line broadening allows constraints on the spin parameter (the X-ray reflection method). Another technique uses accretion disk models to fit the AGN continuum emission (the continuum-fitting, or CF, method). Although each technique has model-dependent uncertainties, these are the best empirical tools currently available and should be vetted in systems where both techniques can be applied. A detailed comparison of the two methods is also useful because neither method can be applied to all AGN. The X-ray reflection technique targets mostly local (z $lesssim$ 0.1) systems, while the CF method can be applied at higher redshift, up to and beyond the peak of AGN activity and growth. Here, we apply the CF method to two AGN with X-ray reflection measurements. For both the high-mass AGN, H1821+643, and the Seyfert 1, NGC 3783, we find a range in spin parameter consistent with the X-ray reflection measurements. However, the near-maximal spin favored by the reflection method for NGC 3783 is more probable if we add a disk wind to the model. Refinement of these techniques, together with improved X-ray measurements and tighter BH mass constraints, will permit this comparison in a larger sample of AGN and increase our confidence in these spin estimation techniques.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
