Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGaussPy+: A fully automated Gaussian decomposition package for emission line spectra
137
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values.
rate research
Read More
We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.
We present a machine learning package for the classification of periodic variable stars. Our package is intended to be general: it can classify any single band optical light curve comprising at least a few tens of observations covering durations from weeks to years, with arbitrary time sampling. We use light curves of periodic variable stars taken from OGLE and EROS-2 to train the model. To make our classifier relatively survey-independent, it is trained on 16 features extracted from the light curves (e.g. period, skewness, Fourier amplitude ratio). The model classifies light curves into one of seven superclasses - Delta Scuti, RR Lyrae, Cepheid, Type II Cepheid, eclipsing binary, long-period variable, non-variable - as well as subclasses of these, such as ab, c, d, and e types for RR Lyraes. When trained to give only superclasses, our model achieves 0.98 for both recall and precision as measured on an independent validation dataset (on a scale of 0 to 1). When trained to give subclasses, it achieves 0.81 for both recall and precision. In order to assess classification performance of the subclass model, we applied it to the MACHO, LINEAR, and ASAS periodic variables, which gave recall/precision of 0.92/0.98, 0.89/0.96, and 0.84/0.88, respectively. We also applied the subclass model to Hipparcos periodic variable stars of many other variability types that do not exist in our training set, in order to examine how much those types degrade the classification performance of our target classes. In addition, we investigate how the performance varies with the number of data points and duration of observations. We find that recall and precision do not vary significantly if the number of data points is larger than 80 and the duration is more than a few weeks. The classifier software of the subclass model is available from the GitHub repository (https://goo.gl/xmFO6Q).
We describe the development of automated emission line detection software for the Fiber Multi-Object Spectrograph (FMOS), which is a near-infrared spectrograph fed by $400$ fibers from the $0.2$ deg$^2$ prime focus field of view of the Subaru Telescope. The software, FIELD (FMOS software for Image-based Emission Line Detection), is developed and tested mainly for the FastSound survey, which is targeting H$alpha$ emitting galaxies at $z sim 1.3$ to measure the redshift space distortion as a test of general relativity beyond $z sim 1$. The basic algorithm is to calculate the line signal-to-noise ratio ($S/N$) along the wavelength direction, given by a 2-D convolution of the spectral image and a detection kernel representing a typical emission line profile. A unique feature of FMOS is its use of OH airglow suppression masks, requiring the use of flat-field images to suppress noise around the mask regions. Bad pixels on the detectors and pixels affected by cosmic-rays are efficiently removed by using the information obtained from the FMOS analysis pipeline. We limit the range of acceptable line-shape parameters for the detected candidates to further improve the reliability of line detection. The final performance of line detection is tested using a subset of the FastSound data; the false detection rate of spurious objects is examined by using inverted frames obtained by exchanging object and sky frames. The false detection rate is $< 1$% at $S/N > 5$, allowing an efficient and objective emission line search for FMOS data at the line flux level of $gtrsim 1.0 times 10^{-16}$[erg/cm$^2$/s].
The hyperfine transitions of the ground-rotational state of the hydroxyl radical (OH) have emerged as a versatile tracer of the diffuse molecular interstellar medium. We present a novel automated Gaussian decomposition algorithm designed specifically for the analysis of the paired on-source and off-source optical depth and emission spectra of these transitions. In contrast to existing automated Gaussian decomposition algorithms, AMOEBA (Automated MOlecular Excitation Bayesian line-fitting Algorithm) employs a Bayesian approach to model selection, fitting all 4 optical depth and 4 emission spectra simultaneously. AMOEBA assumes that a given spectral feature can be described by a single centroid velocity and full width at half-maximum, with peak values in the individual optical depth and emission spectra then described uniquely by the column density in each of the four levels of the ground-rotational state, thus naturally including the real physical constraints on these parameters. Additionally, the Bayesian approach includes informed priors on individual parameters which the user can modify to suit different data sets. Here we describe AMOEBA and evaluate its validity and reliability in identifying and fitting synthetic spectra with known parameters.
High-resolution optical integral field units (IFUs) are rapidly expanding our knowledge of extragalactic emission nebulae in galaxies and galaxy clusters. By studying the spectra of these objects -- which include classic HII regions, supernova remnants, planetary nebulae, and cluster filaments -- we are able to constrain their kinematics (velocity and velocity dispersion). In conjunction with additional tools, such as the BPT diagram, we can further classify emission regions based on strong emission-line flux ratios. LUCI is a simple-to-use python module intended to facilitate the rapid analysis of IFU spectra. LUCI does this by integrating well-developed pre-existing python tools such as astropy and scipy with new machine learning tools for spectral analysis (Rhea et al. 2020). Furthermore, LUCI provides several easy-to-use tools to access and fit SITELLE data cubes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
