Synthetic Spectra Research Articles

From theory to practice: Unlocking the distribution of capacitive times in electrochemical impedance spectroscopy

The distribution of relaxation times (DRT) method is a non-parametric approach for analyzing electrochemical impedance spectroscopy (EIS) data. However, we must be careful when using the DRT method on electrochemical systems with blocking electrodes, such as those encountered in batteries and supercapacitors. This is because, at low frequencies, the asymptotic behavior of the DRT model cannot capture unbounded impedances. To address this issue, we explore the distribution of capacitive times (DCT), a method that, despite being developed decades ago, is still not widely used. In this work, we detail the theoretical underpinnings of the DCT, deriving DCT-specific analytical formulae based on several standard impedance models. We also draw parallels between DCT and DRT and show how these two methods differ in capturing timescales and peaks, elucidating the scenarios where DCT can serve as a viable alternative should the DRT not be applicable. Additionally, we develop a novel method featuring two deep neural networks for DCT deconvolution. We systematically tested this method using a diverse set of synthetic and actual EIS spectra to ensure its efficacy and reliability. Overall, this article seeks to expand the scope of non-parametric approaches for EIS data analysis, particularly to systems characterized by blocking electrodes.

The EDIBLES survey

Context. Small linear carbon chain radicals such as C2 and C3 act as both the building blocks and dissociation fragments of larger carbonaceous species. Their rotational excitation traces the temperature and density of local environments. However, these homo-nuclear di- and triatomic species are only accessible through their electronic and vibrational features because they lack a permanent dipole moment, and high signal-to-noise ratio data are necessary as the result of their generally low abundances in the interstellar medium (ISM). Aims. In order to improve our understanding of small carbonaceous species in the ISM, we carried out a sensitive survey of C2 and C3 using the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES) dataset. We also expanded our searches to C4, C5, and the 13C12C isotopologue in the most molecule-rich sightlines. Methods. We fitted synthetic spectra generated following a physical excitation model to the C2 (2-0) Phillips band to obtain the C2 column density (N) as well as the kinetic temperature (Tkin) and number density (n) of the host cloud. The C3 molecule was measured through its à − $ \tilde X$ (000-000) electronic origin band system. We simulated the excitation of this band with a double-temperature Boltzmann distribution. Results. We present the largest combined survey of C2 and C3 to date in which the individual transitions can be resolved. In total, we detected C2 in 51 velocity components along 40 sightlines, and C3 in 31 velocity components along 27 sightlines. Further analysis confirms the two molecules are detected in the same velocity components. We find a very good correlation between N(C2) and N(C3) with a Pearson correlation coefficient r = 0.93 and an average N(C2)/N(C3) ratio of 15.5± 1.4. A comparison with the behaviour of the C2 diffuse interstellar bands (DIBs) shows that there are no clear differences among sightlines with and without detections of C2 and C3. This is in direct contrast to the better-studied non-C2 DIBs, which have reduced strengths in molecule-rich environments, consistent with the idea that the C2 DIBs are indeed a distinguishable DIB family. We also identify, for the first time, the Q(2), Q(3), and Q(4) transitions of the 13C12C (2-0) Phillips band in the stacked average spectrum of molecule-rich sightlines, and estimate the isotopic ratio of carbon 12C/13C to be 79±8, consistent with literature results. At this stage it is not yet possible to identify these transitions in individual sightlines. Our search for the C4 and C5 optical bands was unsuccessful; even in stacked spectra no unambiguous identification could be made.

Experimental characterization of the active and passive fast-ion H-alpha emission in W7-X using FIDASIM

This paper presents the first results from the analysis of Balmer-alpha spectra at Wendelstein 7-X which contain the broad charge exchange emission from fast-ions. The measured spectra are compared to synthetic spectra predicted by the FIDASIM code, which has been supplied with the 3D magnetic fields from VMEC, 5D fast-ion distribution functions from ASCOT, and a realistic Neutral Beam Injection geometry including beam particle blocking elements. Detailed modeling of the beam emission shows excellent agreement between measured beam emission spectra and predictions. In contrast, modeling of beam halo radiation and Fast-Ion H-Alpha signals (FIDA) is more challenging due to strong passive contributions. While about 50% of the halo radiation can be attributed to passive signals from edge neutrals, the FIDA emission—in particular for an edge-localized line of sights—is dominated by passive emission. This is in part explained by high neutral densities in the plasma edge and in part by edge-born fast-ion populations as demonstrated by detailed modeling of the edge fast-ion distribution.

He abundance in NGC 1850 A and B: Are we observing the early stage of the formation of multiple populations in a stellar cluster?

ABSTRACT We present the results for a sample of B stars in the Large Magellanic Cloud young double stellar cluster NGC 1850 A and NGC 1850 B, as observed with the integral-field spectrograph at the Very Large Telescope, the Multi Unit Spectroscopic Explorer (MUSE). We compare the observed equivalent widths (EWs) of four He lines (4922, 5015, 6678, and 7065 Å) with those determined from synthetic spectra computed with different He mass fractions (Y = 0.25, 0.27, 0.30, and 0.35) with the code synspec, which takes into account the non-local thermodynamic equilibrium effect. From this comparison, we determine the He mass fraction of the B stars, finding a distribution that is not homogeneous. The stars can be divided in three groups: He-weak (Y < 0.24) and He-normal (0.24 ≤ Y ≤ 0.26) stars, belonging to the main sequence of NGC 1850 A, and He-rich stars (0.33 ≤ Y ≤ 0.38), situated in the main sequence associated with NGC 1850 B. We analyse the stellar rotation as possibly being responsible for the anomalous features of the He lines in the He-rich stars. We provide a simple analysis of the differences between the observed EWs and those obtained from theoretical models with different rotation velocities (Vsini = 0 and 250 km s–1). The resolution of the MUSE spectra does not allow us to obtain a conclusive result; however, our analysis support the He-enhanced hypothesis.

On the Identification of N-rich Metal-poor Field Stars with Future Chinese Space Station Telescope

During the long term evolution of globular clusters (GCs), some member stars are lost to the field. The recently found nitrogen-rich (N-rich) metal-poor field stars are promising candidates of these GC escapees, since N enhancement is the fingerprint of chemically enhanced populations in GCs. In this work, we discuss the possibility of identifying N-rich metal-poor field stars with the upcoming Chinese Space Station Telescope (CSST). We focus on the main survey camera with NUV, u, g, r, i, z, y filters and slitless spectrograph with a resolution about 200. The combination of UV sensitive equipment and prominent N-related molecular lines in the UV band bodes well for the identification: the color–color diagram of (u − g) versus (g − r) is capable of separating N-rich field stars from normal halo stars, if metallicity can be estimated without using the information on u-band photometry. Besides, the synthetic spectra show that a signal-to-noise ratio of 10 is sufficient to identify N-rich field stars. In the near future, a large sample of N-rich field stars found by CSST, combined with state-of-the-art N-body simulations will be crucial to deciphering GC-Galaxy co-evolution.

Minimum Noise Fraction Analysis of TGO/NOMAD LNO Channel High-Resolution Nadir Spectra of Mars

NOMAD is a suite of spectrometers on the board of the ESA-Roscosmos Trace Gas Orbiter (TGO) spacecraft and is capable of investigating the Martian environment at very high spectral resolution in the ultraviolet–visible and infrared spectral ranges by means of three separate channels: UVIS (0.2–0.65 μm), LNO (2.2–3.8 μm), and SO (2.3–4.3 μm). Among all channels, LNO is the only one operating at infrared wavelengths in nadir-viewing geometry, providing information on the whole atmospheric column and on the surface. Unfortunately, the LNO data are characterized by an overall low level of signal-to-noise ratio (SNR), limiting their contribution to the scientific objectives of the TGO mission. In this study, we assess the possibility of enhancing LNO nadir data SNR by applying the Minimum Noise Fraction (MNF), a well-known algorithm based on the Principal Components technique that has the advantage of providing transform eigenvalues ordered with increasing noise. We set up a benchmark process on an ensemble of synthetic spectra in order to optimize the algorithm specifically for LNO datasets. We verify that this optimization is limited by the presence of spectral artifacts introduced by the MNF itself, and the maximum achievable SNR is dependent on the scientific purpose of the analysis. MNF application study cases are provided to LNO data subsets in the ranges 2.627–2.648 μm and 2.335–2.353 μm (spectral orders 168 and 189, respectively) covering absorption features of gaseous H2O and CO and CO2 ice, achieving a substantial enhancement in the quality of the observations, whose SNR increases up to a factor of 10. While such an enhancement is still not enough to enable the investigation of spectral features of faint trace gases (in any case featured in orders whose spectral calibration is not fully reliable, hence preventing the application of the MNF), interesting perspectives for improving retrieval of both atmospheric and surface features from LNO nadir data are implied.

Search for spectroscopic binaries using rotational velocities in five open clusters observed by ESO

ABSTRACT In this paper we detect double-lined spectroscopic binaries (SB2) in five open clusters: NGC 2243, NGC 2420, NGC 3532, NGC 6253, and NGC 6705 (M 11) using a method based on high values of the projected rotational velocity when they are fitted with single star spectral model. Observed spectra were obtained from the European Southern Observatory archive. The method was validated on the sets of synthetic spectra for the single and binary stars. It is able to reliably select spectroscopic binaries without confusing them with single stars, if the components in a binary rotate slowly and the radial velocity separation is sufficiently high. We found 60 SB2 candidates: two in NGC 2243, eight each in NGC 2420 and NGC 3532, 17 in NGC 6253, and 25 in NGC 6705. Comparison with literature confirms 18 of them, thus we found 42 new SB2 candidates.

Examining quasar absorption-line analysis methods: the tension between simulations and observational assumptions key to modelling clouds

ABSTRACT A key assumption in quasar absorption-line studies of the circumgalactic medium (CGM) is that each absorption component maps to a spatially isolated ‘cloud’ structure that has single valued properties (e.g. density, temperature, metallicity). We aim to assess and quantify the degree of accuracy underlying this assumption. We used adaptive mesh refinement hydrodynamic cosmological simulations of two z = 1 dwarf galaxies and generated synthetic quasar absorption-line spectra of their CGM. For the Si ii λ1260 transition, and the C iv λλ1548, 1550 and O vi λλ1031, 1037 fine-structure doublets, we objectively determined which gas cells along a line of sight (LOS) contribute to detected absorption. We implemented a fast, efficient, and objective method to define individual absorption components in each absorption profile. For each absorption component, we quantified the spatial distribution of the absorbing gas. We studied a total of 1302 absorption systems containing a total of 7755 absorption components. 48 per cent of Si ii, 68 per cent of C iv, and 72 per cent of O vi absorption components arise from two or more spatially isolated ‘cloud’ structures along the LOS. Spatially isolated ‘cloud’ structures were most likely to have cloud–cloud LOS separations of 0.03Rvir (1.3 kpc), 0.11Rvir (4.8 kpc), and 0.13Rvir (5.6 kpc) for Si ii, C iv, and O vi, respectively. There can be very little overlap between multiphase gas structures giving rise to absorption components. If our results reflect the underlying reality of how absorption lines record CGM gas, they place tension on current observational analysis methods as they suggest that component-by-component absorption-line formation is more complex than is assumed and applied for chemical-ionization modelling.

The First Y Dwarf Data from JWST Show that Dynamic and Diabatic Processes Regulate Cold Brown Dwarf Atmospheres

The James Webb Space Telescope (JWST) is now observing Y dwarfs, the coldest known brown dwarfs, with effective temperatures T eff ≲ 475 K. The first published observations provide important information: not only is the atmospheric chemistry out of equilibrium, as previously known, but the pressure–temperature profile is not in the standard adiabatic form. The rapid rotation of these Jupiter-size, isolated, brown dwarfs dominates the atmospheric dynamics, and thermal and compositional changes disrupt convection. These processes produce a colder lower atmosphere, and a warmer upper atmosphere, compared to a standard adiabatic profile. Leggett et al. presented empirical models where the pressure–temperature profile was adjusted so that synthetic spectra reproduced the 1 ≲ λ (μm) ≲ 20 spectral energy distributions of brown dwarfs with 260 ≤ T eff (K) ≤ 540. We show that spectra generated by these models fit the first JWST Y dwarf spectrum better than standard-adiabat models. Unexpectedly, there is no 4.3 μm PH3 feature in the JWST spectrum and atmospheres without phosphorus better reproduce the 4 μm flux peak. Our analysis of new JWST photometry indicates that the recently discovered faint secondary of the WISE J033605.05-014350AB system has T eff ≈ 295 K, making it the first dwarf in the significant luminosity gap between the 260 K WISE J085510.83-071442.5, and all other known Y dwarfs. The adiabat-adjusted disequilibrium-chemistry models are recommended for analyses of all brown dwarfs cooler than 600 K, and a grid is publicly available. Photometric color transformations are provided in an appendix.

Day–Night Transport-induced Chemistry and Clouds on WASP-39b: Gas-phase Composition

JWST has recently detected the first robust photochemical product on an exoplanet: sulfur dioxide (SO2) on WASP-39b. The data from the NIRISS instrument also reveal signs of partial coverage of clouds. Most of the previous studies have focused on interpreting spectral data with 1D models. To explore how the chemical species and cloud particles are altered by global circulation, we applied a 2D photochemical model and a 2D microphysical cloud model separately to postprocess the thermal and dynamical structures simulated by a 3D general circulation model of WASP-39b. We found that SO2 produced by photochemistry on the dayside can be transported to the nightside owing to the efficient replenishment of horizontal transport. The morning–evening limb differences in methane (CH4) abundances predicted by the 1D models disappeared after horizontal transport is included. Similarly, the inclusion of horizontal transport also reduced the limb differences in SO2. Our modeling results suggest that the fast zonal wind results in minimal or negligible limb asymmetry in composition. Based on the synthetic spectra generated by our 2D atmosphere simulations, we propose that observing SO2 absorption in the emission spectra of WASP-39b at different phases may offer opportunities to probe the horizontal quenching process of photochemical products. We will focus on the gas-phase chemistry in this paper and leave the results regarding clouds in the subsequent paper as part of the series.

Surface chemical composition of single WNh stars

Context. Wolf–Rayet (WR) stars of the WNh category contain a significant fraction of hydrogen at their surface. They can be hydrogen-burning, very massive stars or stars in a post-main sequence phase of evolution. Also, WNh stars are sometimes not included in population synthesis models. Aims. We aim to better characterise the properties of single WNh stars in the Galaxy and the Magellanic Clouds. In particular, we want to constrain their surface chemistry beyond the hydrogen content by determining the helium, carbon, and nitrogen surface abundances. Methods. We perform a spectroscopic analysis of 22 single WNh stars. We fit their ultraviolet and/or optical spectra using synthetic spectra computed with the code CMFGEN. We determine the main stellar parameters (temperature, luminosity, mass-loss rates) and the surface H, He, C, and N mass fractions. We investigate the ability of current evolutionary models to reproduce all parameters at the same time. Results. We find that all WNh stars show the signatures of CNO-cycle material at their surface: they are carbon-depleted and nitrogen-rich. A clear trend of higher nitrogen content at higher metallicity is observed, as expected. The amount of hydrogen (X) varies significantly from one star to another, independently of luminosity. Values of X larger than 0.4 are not exceptional. The majority of Galactic WNh stars can be explained by evolutionary models, provided sufficient fine-tuning of the input parameters of evolutionary calculations. At lower metallicity, most stars escape predictions from evolutionary models. This has been noted in the literature but constraints on the surface nitrogen content exacerbate this severe issue. Conclusions. Our study highlights the need to refine the treatment of WR stars in both stellar evolution and population synthesis models.

Ground motion simulations of the 1887 Pinal de Amoles (MW = 5.53) and the 1989 Landa de Matamoros (MW = 4.94) earthquakes, Querétaro, Mexico

Ground motions characteristics of the 1887 Pinal de Amoles (MW = 5.53) and the 1989 Landa de Matamoros (MW = 4.94) Querétaro earthquakes were studied. For this purpose, the specific barrier model (SBM) was implemented in the context of the stochastic ground motion representation of earthquakes. The SBM was calibrated by modeling the peak ground acceleration (PGA) and the Fourier acceleration spectra (FAS) of the 11 May 2016 (MW = 4.80) Guadalajara earthquake. The calibration results showed that the SBM and selected source and path parameters were able to reproduce synthetic PGA and spectra in good agreement with observed waveforms; in particular, a stress drop of 160 bars fits the results better. Several rupture scenarios for the 1887 and 1989 events were performed by obtaining synthetic seismograms considering different values of stress drop (10, 30, 50, 100, 150, 200, 300, and 500 bars). For a typical stress-drop value of 100 bars, the PGA at distances of 10 and 200 km are in the following intervals: 2.57–295 cm/s2 and 1.72–186 cm/s2 for the Pinal de Amoles and Landa de Matamoros earthquakes, respectively. The results obtained in this study contribute toward a better understanding of strong ground motions in the region.

Estimating stellar parameters and identifying very metal-poor stars for low-resolution spectra (R ∼ 200)

Abstract Very metal-poor (VMP, [Fe/H]<-2.0) stars serve as invaluable repositories of insights into the nature and evolution of the first-generation stars formed in the early galaxy. The upcoming China Space Station Telescope (CSST) will provide us with a large amount of spectral data that may contain plenty of VMP stars, and thus it is crucial to determine the stellar atmospheric parameters ( $T_{\textrm{eff}}$ , $\log$ g, and [Fe/H]) for low-resolution spectra similar to the CSST spectra ( $R\sim 200$ ). This study introduces a novel two-dimensional Convolutional Neural Network (CNN) model, comprised of three convolutional layers and two fully connected layers. The model’s proficiency is assessed in estimating stellar parameters, particularly metallicity, from low-resolution spectra ( $R \sim 200$ ), with a specific focus on enhancing the search for VMP stars within the CSST spectral data. We mainly use 10 008 spectra of VMP stars from LAMOST DR3, and 16 638 spectra of non-VMP stars ([Fe/H]>-2.0) from LAMOST DR8 for the experiments and apply random forest and support vector machine methods to make comparisons. The resolution of all spectra is reduced to $R\sim200$ to match the resolution of the CSST, followed by pre-processing and transformation into two-dimensional spectra for input into the CNN model. The validation and practicality of this model are also tested on the MARCS synthetic spectra. The results show that using the CNN model constructed in this paper, we obtain Mean Absolute Error (MAE) values of 99.40 K for $T_{\textrm{eff}}$ , 0.22 dex for $\log$ g, 0.14 dex for [Fe/H], and 0.26 dex for [C/Fe] on the test set. Besides, the CNN model can efficiently identify VMP stars with a precision rate of 94.77%, a recall rate of 93.73%, and an accuracy of 95.70%. This paper powerfully demonstrates the effectiveness of the proposed CNN model in estimating stellar parameters for low-resolution spectra ( $R\sim200$ ) and recognizing VMP stars that are of interest for stellar population and galactic evolution work.

Searching for Novel Chemistry in Exoplanetary Atmospheres Using Machine Learning for Anomaly Detection

The next generation of telescopes will yield a substantial increase in the availability of high-quality spectroscopic data for thousands of exoplanets. The sheer volume of data and number of planets to be analyzed greatly motivate the development of new, fast, and efficient methods for flagging interesting planets for reobservation and detailed analysis. We advocate the application of machine learning (ML) techniques for anomaly (novelty) detection to exoplanet transit spectra, with the goal of identifying planets with unusual chemical composition and even searching for unknown biosignatures. We successfully demonstrate the feasibility of two popular anomaly detection methods (local outlier factor and one-class support vector machine) on a large public database of synthetic spectra. We consider several test cases, each with different levels of instrumental noise. In each case, we use receiver operating characteristic curves to quantify and compare the performance of the two ML techniques.

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Measurements of the turbulent density wavenumber spectrum, δnˆe(k⊥) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, νe∗<1 , Te/Ti>1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core (ρ ≈ 0.7) over a broad wavenumber range, 0.5⩽k⊥⩽16 cm−1 ( 0.1⩽k⊥ρs⩽5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, Ps(k⊥) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δnˆe(k⊥) . The synthetic DBS Ps(k⊥) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δnˆe(k⊥) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δnˆe(k⊥) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay (k −0.6) for k < 3.5 cm−1, with k −2.6 at intermediate-k ( 3.5⩽k⩽8.5 cm−1), and rapid decay (k −9.4) for k > 8.5 cm−1. Scans of physics parameters using TGLF suggest that the normalized ∇Te scale-length, R/LTe , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k) while R/LTe -driven TEM/ETG-type modes (intermediate/high k) are marginally sub-dominant.

Simulated Radio and Neutrino Imaging of a Microquasar

Microquasar stellar systems emit electromagnetic radiation and high-energy particles. Thanks to their location within our own galaxy, they can be observed in high detail. Still, many of their inner workings remain elusive; hence, simulations, as the link between observations and theory, are highly useful. In this paper, both high-energy particle and synchrotron radio emissions from simulated microquasar jets are calculated using special relativistic imaging. A finite ray speed imaging algorithm is employed on hydrodynamic simulation data, producing synthetic images seen from a stationary observer. A hydrodynamical model is integrated in the above emission models. Synthetic spectra and maps are then produced that can be compared to observations from detector arrays. As an application, the model synthetically observes microquasars during an episodic ejection at two different spatio-temporal scales: one on the neutrino emission region scale and the other on the synchrotron radio emission scale. The results are compared to the sensitivity of existing detectors.

Study of Electron Impact Excitation of Na-like Kr Ion for Impurity Seeding Experiment in Large Helical Device

In this work, a krypton gas impurity seeding experiment was conducted in a Large Helical Device. Emission lines from the Na-like Kr ion in the extreme ultraviolet wavelength region, such as 22.00 nm, 17.89 nm, 16.51 nm, 15.99 nm, and 14.08 nm, respective to 2p63p(2P1/2o)−2p63s(2S1/2), 2p63p(2P3/2o)−2p63s(2S1/2), 2p63d(2D3/2)−2p63p(2P3/2o), 2p63d(2D5/2)−2p63p(2P3/2o), and 2p63d(2D3/2)−2p63p(2P1/2o) transitions, are observed. In order to generate a theoretical synthetic spectrum, an extensive calculation concerning the excitation of the Kr25+ ion through electron impact was performed for the development of a suitable plasma model. For this, the relativistic multiconfiguration Dirac–Hartree–Fock method was employed along with its extension to the relativistic configuration interaction method to compute the relativistic bound-state wave functions and excitation energies of the fine structure levels using the General Relativistic Atomic Structure Package-2018. In addition, another set of calculations was carried out utilizing the relativistic many-body perturbation theory and relativistic configuration interaction methods integrated within the Flexible Atomic Code. To investigate the reliability of our findings, the results of excitation energies, transition probabilities, and weighted oscillator strengths of different dipole-allowed transitions obtained from these different methods are presented and compared with the available data. Further, the detailed electron impact excitation cross-sections and their respective rate coefficients are obtained for various fine structure resolved transitions using the fully relativistic distorted wave method. Rate coefficients, calculated using the Flexible Atomic Code for population and de-population kinetic processes, are integrated into the collisional-radiative plasma model to generate a theoretical spectrum. Further, the emission lines observed from the Kr25+ ion in the impurity seeding experiment were compared with the present plasma model spectrum, demonstrating a noteworthy overall agreement between the measurement and the theoretical synthetic spectrum.

Recovering non-Maxwellian particle velocity distribution functions from collective Thomson-scattered spectra

Collective optical Thomson scattering (TS) is a diagnostic commonly used to characterize plasma parameters. These parameters are typically extracted by a fitting algorithm that minimizes the difference between a measured scattered spectrum and an analytic spectrum calculated from the velocity distribution function (VDF) of the plasma. However, most existing TS analysis algorithms assume that the VDFs are Maxwellian, and applying an algorithm that makes this assumption does not accurately extract the plasma parameters of a non-Maxwellian plasma due to the effect of non-Maxwellian deviations on the TS spectra. We present new open-source numerical tools for forward modeling analytic spectra from arbitrary VDFs and show that these tools are able to more accurately extract plasma parameters from synthetic TS spectra generated by non-Maxwellian VDFs compared to standard TS algorithms. Estimated posterior probability distributions of fits to synthetic spectra for a variety of example non-Maxwellian VDFs are used to determine uncertainties in the extracted plasma parameters and show that correlations between parameters can significantly affect the accuracy of fits in plasmas with non-Maxwellian VDFs.

Influence of the optical Fe II quasi-continuum on measuring the spectral parameters of active galactic nuclei

Aims. We explore the influence of optical Fe II quasi-continuum on the measured spectral parameters in the λλ4150−5500 Å range for the spectra of Type 1 active galactic nuclei (AGNs). Methods. We assume that the broad line region is composed of two sub-regions: the very broad line region (VBLR) and the intermediate line region (ILR). We constructed a large set of synthetic AGN spectra by taking different portions of the VBLR and ILR contributions, where initially the VBLR and ILR model spectra were constructed on the basis of prototypes of two observed spectra with dominant VBLR (i.e. ILR) emission. To investigate the influence of the optical Fe II quasi-continuum on the AGN measured spectral parameters, we fit the power-law continuum and emission lines in a set of model spectra, as commonly done for observed AGN spectra. We then compared the spectral parameters obtained after the fitting procedure with those of the model. Results. We find that the optical Fe II quasi-continuum can be very strong in the case of spectra with strong and very broad Fe II lines and it is difficult to fully separate it from the power-law continuum. This gives the effect of a slightly underestimated Hβ width and underestimated fluxes of the Hβ and Fe II lines, while the continuum flux is then slightly overestimated. The most affected spectral parameters are the line equivalent widths (EWs), especially EW Fe II, which may be strongly underestimated. We discuss the possible underlying physics in the quasar main sequence, as implied by the results of our spectral modelling. We find that the set of AGN model spectra assuming different ILR and VBLR contributions can aptly reproduce the quasar main sequence, that is, the full width at half maximum Hβ versus Fe II/Hβ anti-correlation, where both parameters in this anti-correlation are strongly dependent on the ILR and VBLR contribution rate.

Unique distant classical Cepheid OGLE GD-CEP-1353 with anomalously high abundances of s- and r-process elements

Aims. While looking for recently discovered distant Cepheids with an interesting chemical composition, we noticed one star (OGLE GD-CEP-1353) with extremely large equivalent widths of spectral lines of heavy elements. The aim of this work is to perform an abundance analysis, and to find a possible explanation for the found chemical anomaly. Methods. Quantitative analysis of the equivalent widths and synthetic spectrum synthesis were used to derive abundances in this star. Both local and nonlocal thermodynamic equilibrium (LTE and NLTE) approximations were used in our analysis. Results. Abundances of 28 chemical elements from carbon to thorium were derived. While light and iron peak elements show abundances typical for distant Cepheids (located in the outer disk), the s-process elements are overabundant about one dex. r-process elements are slightly less overabundant. This makes the star a unique Cepheid of our Galaxy.