Model Line Profiles Research Articles

Deriving X-Ray Line Profiles for Massive-star Winds from Momentum-conserving Dynamical Working Surface Solutions

We present a general procedure for deriving a line-profile model for massive-star X-ray spectra that captures the dynamics of the wind more directly. The basis of the model is the analytic solution to the problem of variable jets in Herbig–Haro objects given by Cantó et al. In deriving our model, we generalize this jet solution to include flows with a prescribed nonzero acceleration for the context of radiatively driven winds. We provide example line profiles generated from our model for the case of sinusoidal velocity and mass-ejection variations. The example profiles show the expected shape of massive-star X-ray emission lines, as well as interesting but complicated trends with the model parameters. This establishes the possibility that observed X-rays could be a result of temporal variations seeded at the wind base, rather than purely generated intrinsically within the wind volume, and can be described via a quantitative language that connects with the physical attributes of those variations, consistently with the downstream momentum-conserving nature of radiatively cooled shocked radial flows.

Towards inferring the geometry of kilonovae

ABSTRACT Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints, suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody, the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4–7 per cent).

A new hydrochloric acid 2-0 band analysis: Two temperature study

A new study of the 2−0 band of the hydrochloric acid molecule was performed by high resolution Fourier-transform absorption spectroscopy in the near infrared region. The spectra were measured at two different temperatures, 293 K and 315 K, for different pressures at each temperature. The spectral linewidths were analysed in a two-step procedure, being first performed by directly measuring the linewidth and second by fitting each spectral line to a model line profile, using Gaussian, Loretzian and Voigt profiles. A study of the profiles that best describe the spectral line fits is carried out in this work. The behavior of the spectral lines self-broadening and their corresponding self-induced shifts were studied for different values of rotational quantum numbers. The analysis was performed for both isotopes of the molecule and the self-broadening and self-shift coefficients were determined.

Infrared Spectroscopy of Be Stars: Influence of the Envelope Parameters on Brackett-Series Behaviour

The IR spectra of Be stars display numerous hydrogen recombination lines, constituting a great resource for obtaining information on the physical and dynamic structures of different regions within the circumstellar envelope. Nevertheless, this spectral region has not been analysed in depth, and there is a lack of synthetic spectra with which to compare observations. Therefore, we computed synthetic spectra with the HDUST code for different disc parameters. Here, we present our results on the spectral region that includes lines of the Brackett series. We discuss the dependence of the line series strengths on several parameters that describe the structure of the disc. We also compared model line profiles, fluxes, and EWs with observational data for two Be stars (MX Pup and π Aqr). Even though the synthetic spectra adequately fit our observations of both stars and allow us to constrain the parameters of the disc, there is a discrepancy with the observed data in the EW and flux measurements, especially in the case of MX Pup. It is possible that by including Brackett lines of higher terms or adding the analysis of other series, we may be able to better constrain the parameters of the observed disc.

Ionized gas kinematics and chemical abundances of low-mass star-forming galaxies at z ∼ 3

Context. Feedback from massive stars plays a crucial role in regulating the growth of young star-forming galaxies (SFGs) and in shaping their interstellar medium (ISM). This feedback contributes to the removal and mixing of metals via galactic outflows and to the clearance of neutral gas, which facilitates the escape of ionizing photons. Aims. Our goal is to study the impact of stellar feedback on the chemical abundances of the ISM in a sample of SFGs with strong emission lines at z ∼ 3. Methods. We selected 35 low-mass SFGs (7.9 < log(M⋆/M⊙) < 10.3) from deep spectroscopic surveys based on their CIII]λ1908 emission. We used new follow-up near-infrared (NIR) observations to examine their rest-optical emission lines and to identify ionized outflow signatures through broad emission line wings detected after Gaussian modeling of [OIII]λλ4959,5007 profiles. We characterized the gas-phase metallicity and carbon-to-oxygen (C/O) abundance of the galaxies using a Te-based method via the OIII]λ1666/[OIII]λ5007 ratio and photoionization models. Results. We find line ratios and rest-frame equivalent widths (EWs) characteristic of high-ionization conditions powered by massive stars. Our sample displays a mean rest-frame EW([OIII]λ5007) of ∼560 Å, while about 15% of the SFGs show EW([OIII]λλ4959,5007) > 1000 Å and EW(CIII]) > 5 Å, closely resembling those now seen in epoch of reionization (EoR) galaxies with the James Webb Space Telescope. We find high Te values, which imply low gas-phase metallicities 12+log(O/H) ∼ 7.5–8.5 (mean of 17% solar) and C/O abundances from 23% to 128% solar, with no apparent increasing trend with metallicity. Our sample follows the mass-metallicity relation at z ∼ 3, with some galaxies showing lower gas-phase metallicities. This results in significant deviations from the fundamental metallicity relation. From our [OIII]λλ4959,5007 line profile modeling, we find that 65% of our sample shows an outflow component, which is found both blue- or redshifted relative to the ionized gas systemic velocity, and the mean maximum velocities are vmax ∼ 280 km s−1. We find a weak correlation between vmax and the star formation rate surface density (ΣSFR) of vmax = (2.41 ± 0.03) × ΣSFR(0.06 ± 0.03). Moreover, we find that the mass-loading factor μ of our galaxy sample is typically lower than in more massive galaxies from the literature, but it is higher than in typical local dwarf galaxies. In the stellar mass range covered by our sample, we find that μ increases with ΣSFR. This suggests that for a given stellar mass, denser starbursts in low-mass galaxies produce stronger outflows. Our results complement the picture drawn by similar studies at lower redshift, suggesting that the removal of ionized gas in low-mass SFGs driven by stellar feedback is regulated by their stellar mass and by the strength and concentration of their star formation, that is, ΣSFR.

Astraeus – VIII. A new framework for Lyman-α emitters applied to different reionization scenarios

ABSTRACT We use the astraeus framework to investigate how the visibility and spatial distribution of Lyman-α (Lyα) emitters (LAEs) during reionization is sensitive to a halo mass-dependent fraction of ionizing radiation escaping from the galactic environment (fesc) and the ionization topology. To this end, we consider the two physically plausible bracketing scenarios of fesc increasing and decreasing with rising halo mass. We derive the corresponding observed Lyα luminosities of galaxies for three different analytic Lyα line profiles and associated Lyα escape fraction ($f_\mathrm{esc}^\mathrm{Ly\alpha }$) models: importantly, we introduce two novel analytic Lyα line profile models that describe the surrounding interstellar medium (ISM) as dusty gas clumps. They are based on parameterizing results from radiative transfer simulations, with one of them relating $f_\mathrm{esc}^\mathrm{Ly\alpha }$ to fesc by assuming the ISM of being interspersed with low-density tunnels. Our key findings are: (i) for dusty gas clumps, the Lyα line profile develops from a central to double peak dominated profile as a galaxy’s halo mass increases; (ii) LAEs are galaxies with $M_h\gtrsim 10^{10}\, {\rm \rm M_\odot }$ located in overdense and highly ionized regions; (iii) for this reason, the spatial distribution of LAEs is primarily sensitive to the global ionization fraction and only weakly in second-order to the ionization topology or a halo mass-dependent fesc; (iv) furthermore, as the observed Lyα luminosity functions reflect the Lyα emission from more massive galaxies, there is a degeneracy between the fesc-dependent intrinsic Lyα luminosity and the Lyα attenuation by dust in the ISM if fesc does not exceed $\sim 50~{{\ \rm per \, cent}}$.

The IAG spectral atlas of the spatially resolved Sun: Centre-to-limb observations

Context. Solar surface magneto-convection appears as granulation pattern that impacts spectral lines in terms of both shape and wavelength. Such induced effects also tend to vary over the observed solar disc because of the changing observation angle and, thus, the changing observation height as well. Centre-to-limb observations of the resolved Sun offer an insight into the variable spectral behaviour across different heliocentric observing positions, providing crucial information about limb darkening, convective velocities, and line profile variability relevant to radial velocity (RV) calculations. Thus, RV measurements and exoplanet transit spectroscopy depend on precise reference templates. Aims. We want to provide a spectroscopic centre-to-limb solar atlas at high spectral resolution and high-frequency accuracy. The atlas shall help improve the understanding of the solar atmosphere and convection processes. Methods. We performed high-resolution observations of the resolved quiet Sun with a Fourier transform spectrograph at the Institut für Astrophysik und Geophysik in Göttingen. Our dataset contains a wavelength range from 4200 Å to 8000 Å. We obtained 165 spectra in total, with a spectral resolution of Δν = 0.024 cm−1, corresponding to a resolving power R of 700 000 at ∼6000 Å. Results. We present a centre-to-limb solar atlas containing 14 heliocentric positions. To check for consistency, we investigated the Fe I 6175 Å absorption line and compared our line profiles with previous centre-to-limb observations and also with simulations. The line profile and also the bisector profiles are generally consistent with previous observations, but we have identified differences to model line profiles, especially close to the solar limb.

Dust masses for a large sample of core-collapse supernovae from optical emission line asymmetries: dust formation on 30-year time-scales

ABSTRACT Modelling the red–blue asymmetries seen in the broad emission lines of core-collapse supernovae (CCSNe) is a powerful technique to quantify total dust mass formed in the ejecta at late times (>5 yr after outburst) when ejecta dust temperatures become too low to be detected by mid-infrared (IR) instruments. Following our success in using the Monte Carlo radiative transfer code damocles to measure the dust mass evolution in SN 1987A and other CCSNe, we present the most comprehensive sample of dust mass measurements yet made with damocles, for CCSNe aged between 4 and 60 yr after outburst. Our sample comprises multi-epoch late-time optical spectra taken with the Gemini/Gemini Multi-Object Spectrographs (GMOS) and Very Large Telescope (VLT) X-Shooter spectrographs, supplemented by archival spectra. For the 14 CCSNe that we have modelled, we confirm a dust mass growth with time that can be fit by a sigmoid curve that is found to saturate beyond an age of ∼30 yr, at a mass of 0.23$^{+0.17}_{-0.12}$ M⊙. For an expanded sample including dust masses found in the literature for a further 11 CCSNe and six CCSN remnants, the dust mass at saturation is found to be 0.42$^{+0.09}_{-0.05}$ M⊙. Uncertainty limits for our dust masses were determined from a Bayesian analysis using the affine invariant Markov chain Monte Carlo ensemble sampler emcee with damocles. The best-fitting line profile models for our sample all required grain radii between 0.1 and 0.5 $\mu$m. Our results are consistent with CCSNe forming enough dust in their ejecta to significantly contribute to the dust budget of the Universe.

Accretion geometry of the black hole binary MAXI J1820+070 probed by frequency-resolved spectroscopy

ABSTRACT The geometry of the inner accretion flow in the hard and hard-intermediate states of X-ray binaries remains controversial. Using Neutron star Interior Composition Explorer observations of the black hole X-ray binary MAXI J1820+070 during the rising phase of its 2018 outburst, we study the evolution of the timing properties, in particular the characteristic variability frequencies of the prominent iron K α line. Using frequency-resolved spectroscopy, which is robust against uncertainties in the line profile modelling, we find that reflection occurs at large distances from the Comptonizing region in the bright hard state. During the hard-to-soft transition, the variability properties suggest that the reflector moves closer to the X-ray source. In parallel, the peak of the iron line shifts from 6.5 to ∼7 keV, becoming consistent with that expected of from a highly inclined disc extending close to the black hole. We additionally find significant changes in the dependence of the root-mean-square (rms) variability on both energy and Fourier frequency as the source softens. The evolution of the rms-energy dependence, the line profile, and the timing properties of the iron line as traced by the frequency-resolved spectroscopy all support the picture of a truncated disc/inner flow geometry.

Re-evaluation of ortho-para-dependence of self pressure-broadening in the ν1 + ν3 band of acetylene.

Optical frequency comb-referenced measurements of self pressure-broadened line profiles of the R(8) to R(13) lines in the ν1 + ν3 combination band of acetylene near 1.52 µm are reported. The analysis of the data found no evidence for a previously reported [Iwakuni et al., Phys. Rev. Lett. 117(14), 143902 (2016)] systematic alternation in self pressure-broadened line widths with the nuclear spin state of the molecule. The present work brought out the need for the use of an accurate line profile model and careful accounting for weak background absorptions due to hot band and lower abundance isotopomer lines. The data were adequately fit using the quadratic speed-dependent Voigt profile model, neglecting the small speed-dependent shift. Parameters describing the most probable and speed-dependent pressure-broadening, most probable shift, and the line strength were determined for each line. Detailed modeling of the results of Iwakuni et al. showed that their neglect of collisional narrowing due to the speed-dependent broadening term combined with the strongly absorbing data recorded and analyzed in transmission mode were the reasons for their results.

Chandra spectral measurements of the O supergiant ζ Puppis indicate a surprising increase in the wind mass-loss rate over 18 yr

ABSTRACT New long Chandra grating observations of the O supergiant ζ Pup show not only a brightening of the X-ray emission line flux of 13 per cent in the 18 yr since Chandra’s first observing cycle, but also clear evidence – at more than 4σ significance – of increased wind absorption signatures in its Doppler-broadened line profiles. We demonstrate this with non-parametric analysis of the profiles as well as Gaussian fitting and then use line-profile model fitting to derive a mass-loss rate of 2.47 ± 0.09 × 10−6${\mathrm{M_{\odot }~{\mathrm{y}r^{-1}}}}$, which is a 40 per cent increase over the value obtained from the cycle 1 data. The increase in the individual emission line fluxes is greater for short-wavelength lines than long-wavelength lines, as would be expected if a uniform increase in line emission is accompanied by an increase in the wavelength-dependent absorption by the cold wind in which the shock-heated plasma is embedded.

Measurement of pressure shift and broadening coefficients for Ne 3s[3/2]2→3p[5/2]3 transition in Ne and He using diode-laser absorption spectroscopy

Measurements of pressure broadening and shift coefficients for 640.2 nm Ne line were performed in an RF discharge without stabilization of neutral gas temperature. A diode laser with a short external cavity was used for metastable neon absorption spectroscopy in a natural mixture of neon isotopes. An appropriate fit function was constructed permitting deduction of Doppler and Lorentz components for the Voigt profile of the spectral line. To determine pressure broadening ξ and shift β coefficients for pure Ne and Ne:He mixture we used the information about Lorentz linewidth, position of the line center and gas temperature derived from spectral line shape fitted to a model line profile. The values of the collisional coefficients reduced to 300 K in units of 10-10 s-1 cm3, are ξNe-Ne = (1.9±0.2), ξNe-He = (3.4±0.2); βNe-Ne = -(0.7±0.1), βNe-He = (0.13±0.04).

Update of the Solar Lyα Profile Line Model

We present a modification of a model of solar cycle evolution of the solar Lyα line profile, along with a sensitivity study of interstellar neutral hydrogen to uncertainties in radiation pressure level. The line profile model, originally developed by Kowalska-Leszczynska et al., is parameterized by the composite solar Lyα flux, which recently was revised. We present modified parameters of the previously developed model of solar radiation pressure for neutral hydrogen and deuterium atoms in the heliosphere. The mathematical function used in the model, as well as the fitting procedure, remain unchanged. We show selected effects of the model modification on interstellar neutral H properties in the heliosphere, and we discuss the sensitivity of these quantities to uncertainties in the calibration of the composite Lyα series.

Modelling the asymmetries of the Sun’s radial p-mode line profiles

Context. The advent of space-borne missions has substantially increased the number and quality of the measured power spectrum of solar-like oscillators. It now allows for the p-mode line profiles to be resolved and facilitates an estimation of their asymmetry. The fact that this asymmetry can be measured for a variety of stars other than the Sun calls for a revisiting of acoustic mode asymmetry modelling. This asymmetry has been shown to be related to a highly localised source of stochastic driving in layers just beneath the surface. However, existing models assume a very simplified, point-like source of excitation. Furthermore, mode asymmetry could also be impacted by a correlation between the acoustic noise and the oscillating mode. Prior studies have modelled this impact, but only in a parametrised fashion, which deprives them of their predictive power. Aims. In this paper, we aim to develop a predictive model for solar radial p-mode line profiles in the velocity spectrum. Unlike the approach favoured by prior studies, this model is not described by free parameters and we do not use fitting procedures to match the observations. Instead, we use an analytical turbulence model coupled with constraints extracted from a 3D hydrodynamic simulation of the solar atmosphere. We then compare the resulting asymmetries with their observationally derived counterpart. Methods. We model the velocity power spectral density by convolving a realistic stochastic source term with the Green’s function associated with the radial homogeneous wave equation. We compute the Green’s function by numerically integrating the wave equation and we use theoretical considerations to model the source term. We reconstruct the velocity power spectral density and extract the line profile of radial p-modes as well as their asymmetry. Results. We find that stochastic excitation localised beneath the mode upper turning point generates negative asymmetry for ν < νmax and positive asymmetry for ν > νmax. On the other hand, stochastic excitation localised above this limit generates negative asymmetry throughout the p-mode spectrum. As a result of the spatial extent of the source of excitation, both cases play a role in the total observed asymmetries. By taking this spatial extent into account and using a realistic description of the spectrum of turbulent kinetic energy, both a qualitative and quantitative agreement can be found with solar observations performed by the GONG network. We also find that the impact of the correlation between acoustic noise and oscillation is negligible for mode asymmetry in the velocity spectrum.

Mirror, mirror on the outflow cavity wall

Aims. The inner regions of high-mass protostars are often invisible in the near-infrared, obscured by thick envelopes and discs. We aim to investigate the inner gaseous disc of IRAS 11101-5829 through scattered light from the outflow cavity walls. Methods. We observed the immediate environment of the high-mass young stellar object IRAS 11101-5829 and the closest knots of its jet, HH135-136, with the integral field unit VLT/SINFONI. We also retrieved archival data from the high-resolution long-slit spectrograph VLT/X-shooter. We analysed imaging and spectroscopic observations to discern the nature of the near-infrared CO emission. Results. We detect the first three bandheads of the υ = 2−0 CO vibrational emission for the first time in this object. It is coincident with continuum and Brγ emission and extends up to ~10 000 au to the north-east and ~10 000 au to the south-west. The line profiles have been modelled as a Keplerian rotating disc assuming a single ring in local thermodynamic equilibrium. The model output gives a temperature of ~3000 K, a CO column density of ~1 × 1022 cm−2, and a projected Keplerian velocity vK sin idisc ~ 25 km s−1, which is consistent with previous modelling in other high-mass protostars. In particular, the low value of vK sin idisc suggests that the disc is observed almost face-on, whereas the well-constrained geometry of the jet imposes that the disc must be close to edge-on. This apparent discrepancy is interpreted as the CO seen reflected in the mirror of the outflow cavity wall. Conclusions. From both jet geometry and disc modelling, we conclude that all the CO emission is seen through reflection by the cavity walls and not directly. This result implies that in the case of highly embedded objects, as for many high-mass protostars, line profile modelling alone might be deceptive and the observed emission could affect the derived physical and geometrical properties; in particular the inclination of the system can be incorrectly interpreted.

Broadening and Shift Coefficients of H2O Absorption Lines in the 8650–9020 cm–1 Spectral Region

The broadening and shift coefficients of H2O absorption lines are compared for Ar, He, H2, and N2 buffer gases. The broadening and shift coefficients were derived from the analysis of the absorption spectra recorded with an FTIR spectrometer in the spectral region 8650–9020 cm−1 with a spectral resolution of 0.01 cm–1. Using two model line profiles (Voigt and speed-dependent Voigt profiles) the parameters of the H2O absorption lines were calculated. It is shown that the speed-dependent Voigt profile provides better agreement with experimental data.

Pressure-dependent sensitivity of a single-pass methane detection system using a continuous-wave distributed feedback laser at 3270 nm.

A single-pass, compact, and low-power-consumption methane sensing system is presented, and its sensitivity is examined for a set of reference cells with pressures of 7.4, 74, and 740Torr and the same methane concentration at laboratory temperature of 23°C. A GaSb-based continuous-wave distributed feedback tunable diode laser at 3270nm and wavelength modulation spectroscopy were employed to collect 2f spectra in the ν3, R3 band of CH412. The collected 2f spectra were fitted to a modulated Voigt line profile model. An Allan-Werle variance analysis shows that the best detectivity, 4ppbm, is obtained for the highest pressure cell. At pressures of 74 and 7.4Torr the detectivities are 16 and 28ppbm, respectively. For these low-pressure cells, further sensitivity improvements were achieved using a large modulation depth (m>2.2) at the expense of spectral resolution.

CIELO-RGS: a catalog of soft X-ray ionized emission lines

Context.High-resolution X-ray spectroscopy has advanced our understanding of the hot Universe by revealing physical properties like kinematics, temperature, and abundances of the astrophysical plasmas. Despite technical and scientific achievements, the lack of scientific products at a level higher than count spectra is hampering complete scientific exploitation of high-quality data. This paper introduces the Catalog of Ionized Emission Lines Observed by the Reflection Grating Spectrometer (CIELO-RGS) onboard theXMM-Newtonspace observatory.Aims.The CIELO-RGS catalog aims to facilitate the exploitation of emission features in the public RGS spectra archive. In particular, we aim to analyze the relationship between X-ray spectral diagnostics parameters and measurements at other wavelengths. This paper focuses on the methodology of catalog generation, describing the automated line-detection algorithm.Methods.A moderate sample (∼2400 observations) of high-quality RGS spectra available atXMM-NewtonScience Archive is used as our starting point. A list of potential emission lines is selected based on a multi-scale peak-detection algorithm in a uniform and automated way without prior assumption on the underlying astrophysical model. The candidate line list is validated via spectral fitting with simple continuum and line profile models. We also compare the catalog content with published literature results on a small number of exemplary sources.Results.We generate a catalog of emission lines (1.2 × 104) detected in ∼1600 observations toward stars, X-ray binaries, supernovae remnants, active galactic nuclei, and groups and clusters of galaxies. For each line, we report the observed wavelength, broadening, energy and photon flux, equivalent width, and so on.

The distance, rotation, and physical parameters of ζ Pup

We scrutinize the Hipparcos parallax for the bright O supergiant zeta Pup, and confirm that the implied distance of 332+/-11pc appears to be reliable. We then review the implications for the star's physical parameters, and the consequences for the interpretation of P(phot), the 1.78-d photometric period. The equatorial rotation period is <3.7d (with 95% confidence), ruling out a proposed ~5.1d value. If the photometric period is the rotation period then i, the inclination of the rotation axis to the line of sight, is 33.2+/-1.8 degrees. The inferred mass, radius, and luminosity are securely established to be less than canonical values for the spectral type, and are not in agreement with single-star evolution models. The runaway status, rapid rotation, and anomalous physical properties are all indicative of an evolutionary history involving binary (or multiple-star) interaction. We perform simple starspot modelling to show that the low axial inclination required if P(rot) = 1.78d has testable spectroscopic consequences, which have not been identified in existing time series. If P(phot) is directly related to drivers of systematic, high-velocity stellar-wind variability (`discrete absorption components') in zeta Pup, antisolar differential rotation is required. Model line profiles calculated on that basis are at variance with observations.

Retrieval of Total Scattering Cross Sections of Molecules from Inhomogeneously Broadened Absorption Lines

A new method of retrieving quantitative information on the hard, soft, and diffraction collision’s frequencies from inhomogeneously pressure-broadened line profiles was proposed and tested. The essence of the method lies in the processing of recorded profiles using various profiles’ models containing these frequencies treated as adjustable parameters in different manners. Three self-broadened H2O absorption lines free of an instrumental function were processed with the use of four model line profiles. Retrieved set of pressure line broadening and narrowing constants allows finding out the sought frequencies. The estimation of the total scattering cross section, ∼630 Å2, in pure water vapor was made.