اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدImproving the full spectrum fitting method: accurate convolution with Gauss-Hermite functions
284
0
0.0
(
0
)
تأليف
Michele Cappellari
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
I start by providing an updated summary of the penalized pixel-fitting (pPXF) method, which is used to extract the stellar and gas kinematics, as well as the stellar population of galaxies, via full spectrum fitting. I then focus on the problem of extracting the kinematic when the velocity dispersion $sigma$ is smaller than the velocity sampling $Delta V$, which is generally, by design, close to the instrumental dispersion $sigma_{rm inst}$. The standard approach consists of convolving templates with a discretized kernel, while fitting for its parameters. This is obviously very inaccurate when $sigma<Delta V/2$, due to undersampling. Oversampling can prevent this, but it has drawbacks. Here I present a more accurate and efficient alternative. It avoids the evaluation of the under-sampled kernel, and instead directly computes its well-sampled analytic Fourier transform, for use with the convolution theorem. A simple analytic transform exists when the kernel is described by the popular Gauss-Hermite parametrization (which includes the Gaussian as special case) for the line-of-sight velocity distribution. I describe how this idea was implemented in a significant upgrade to the publicly available pPXF software. The key advantage of the new approach is that it provides accurate velocities regardless of $sigma$. This is important e.g. for spectroscopic surveys targeting galaxies with $sigmallsigma_{rm inst}$, for galaxy redshift determinations, or for measuring line-of-sight velocities of individual stars. The proposed method could also be used to fix Gaussian convolution algorithms used in todays popular software packages.
قيم البحث
اقرأ أيضاً
We investigate the precision of the ages and metallicities of 21,000 mock simple stellar populations (SSPs) determined through full-spectrum fitting. The mock SSPs cover an age range of 6.8 $<$ log (age/yr) $<$ 10.2, for three wavelength ranges in th
e optical regime, using both Padova and MIST isochrone models. Random noise is added to the model spectra to achieve S/N ratios between 10 to 100 per wavelength pixel. We find that for S/N $geq$ 50, this technique can yield ages of SSPs to an overall precision of $Delta,mbox{log(age/yr)} sim 0.1$ for ages in the ranges 7.0 $leq$ log (age/yr) $leq$ 8.3 and 8.9 $leq$ log (age/yr) $leq$ 9.4. For the age ranges of 8.3 $leq$ log (age/yr) $leq$ 8.9 and log (age/yr) $geq$ 9.5, which have significant flux contributions from asymptotic giant branch (AGB) and red giant branch (RGB) stars, respectively, the age uncertainty rises to about $pm 0.3$ dex. The precision of age and metallicity estimation using this method depends significantly on the S/N and the wavelength range used in the fitting. We quantify the systematic differences in age predicted by the MIST and Padova isochrone models, due to their different assumptions about stellar physics in various important (i.e., luminous) phases of stellar evolution, which needs to be taken in consideration when comparing ages of star clusters obtained using these popular models. Knowing the strengths and limitations of this technique is crucial in interpreting the results obtained for real star clusters and for deciding the optimal instrument setup before performing the observations.
We propose an alternative approach to the construction of fitting functions to the nonlinear matter power spectrum extracted from $N$-body simulations based on the relative matter power spectrum $delta(k,a)$, defined as the fractional deviation in th
e absolute matter power spectrum produced by a target cosmology away from a reference $Lambda$CDM prediction. From the computational perspective, $delta(k,a)$ is fairly insensitive to the specifics of the simulation settings, and numerical convergence at the 1%-level can be readily achieved without the need for huge computing capacity. Furthermore, $delta(k,a)$ exhibits several interesting properties that enable a piece-wise construction of the full fitting function, whereby component fitting functions are sought for single-parameter variations and then multiplied together to form the final product. Then, to obtain 1%-accurate absolute power spectrum predictions for any target cosmology only requires that the community as a whole invests in producing one single ultra-precise reference $Lambda$CDM absolute power spectrum, to be combined with the fitting function to produce the desired result. To illustrate the power of this approach, we have constructed the fitting function RelFit using only five relatively inexpensive $w$CDM simulations (box length $L=256 h^{-1}$Mpc, $N=1024^3$ particles, initialised at $z_i=49$). In a 6-parameter space spanning ${omega_m,A_s,n_s,w,omega_b,h}$, the output relative power spectra of RelFit are consistent with the predictions of the CosmicEmu emulator to 1% or better for a wide range of cosmologies up to $ksimeq 10$/Mpc. Thus, our approach could provide an inexpensive and democratically accessible route to fulfilling the 1%-level accuracy demands of the upcoming generation of large-scale structure probes, especially in the exploration of non-standard or exotic cosmologies on nonlinear scales.
A recent paper by Ge et al. performs a series of experiments with two full spectral fitting codes, pPXF and starlight, finding that the two yield consistent results when the input spectrum is not heavily reddened. For E(B-V) > 0.2, however, they clai
m starlight leads to severe biases in the derived properties. Counterintuitively, and at odds with previous simulations, they find that this behaviour worsens significantly as the signal-to-noise ratio of the input spectrum increases. This communication shows that this is entirely due to an A_V < 1 mag condition imposed while initializing the Markov chains in the code. This choice is normally irrelevant in real-life galaxy work but can become critical in artificial experiments. Alleviating this usually harmless initialization constraint changes the Ge et al. results completely, as was explained to the authors before their publication. We replicate their spectral fitting experiments, finding much smaller biases. Furthermore both bias and scatter in the derived properties all converge as S/N increases, as one would expect. We also show how the very output of the code provides ways of diagnosing anomalies in the fits. The code behaviour has been documented in careful and extensive experiments in the literature, but the biased analysis of Ge et al. is just not representative of starlight at all.
We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are det
ected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 {mu}m allows us to fit a single temperature greybody with temperature 80.3pm0.6K and dust emissivity index 1.73pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clumps bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.
In the Gaia era, understanding the effects of perturbations of the Galactic disc is of major importance in the context of dynamical modelling. In this theoretical paper, we extend previous work in which, making use of the epicyclic approximation, the
linearized Boltzmann equation had been used to explicitly compute, away from resonances, the perturbed distribution function of a Galactic thin disc population in the presence of a non-axisymmetric perturbation of constant amplitude. Here we improve this theoretical framework in two distinct ways in the new code that we present. First, we use better estimates for the action-angle variables away from quasi-circular orbits, computed from the AGAMA software, and we present an efficient routine to numerically re-express any perturbing potential in these coordinates with an accuracy well below the percent level. The use of more accurate action estimates allows us to identify resonances such as the outer 1:1 bar resonance at larger azimuthal velocities than the outer Lindblad resonance, and to extend our previous theoretical results well above the Galactic plane, where we explicitly show how they differ from the epicyclic approximation. In particular, the displacement of resonances in velocity space as a function of height can in principle constrain the 3D structure of the Galactic potential. Second, we allow the perturbation to be time-dependent, thereby allowing us to model the effect of transient spiral arms or of a growing bar. The theoretical framework and tools presented here will be useful for a thorough analytical dynamical modelling of the complex velocity distribution of disc stars as measured by past and upcoming Gaia data releases.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
