Line Wings Research Articles

Optimising the Hα index for the identification of activity signals in FGK stars

Context. The Ca II H&K and Hα lines are two of the most used activity diagnostics for detecting stellar activity signals in the optical regime, and for inferring possible false positives in exoplanet detection with the radial velocity method. The flux in the two lines is known to follow the solar activity cycle, and to correlate well with sunspot number and other activity diagnostics. However, for other stars, the flux in these lines is known to have a wide range of correlations, increasing the difficulty in the interpretation of the signals observed with the Hα line. Aims. In this work we investigate the effect of the Hα bandpass width on the correlation between the Ca II and Hα indices with the aim of improving the Hα index to better identify and model the signals coming from activity variability. Methods. We used a sample of 152 FGK dwarfs observed with HARPS for more than 13 yr with enough cadence to be able to detect rotational modulations and cycles in activity proxies. We calculated the Ca II and Hα activity indices using a range of bandwidths for Hα between 0.1 and 2.0 Å. We studied the correlation between the indices’ time series at long and short timescales, and analysed the impact of stellar parameters, activity level, and variability on the correlations. Results. The correlation between Ca II and Hα, both at short and long timespans, is maximised when using narrow Hα bandwidths, with a maximum at 0.6 Å. For some inactive stars, as the activity level increases, the flux in the Hα line core increases, while the flux in the line wings decreases as the line becomes shallower and broader. The balance between these fluxes can cause stars to show the negative correlations observed in the literature when using a wide bandwidth on Hα. These anti-correlations may become positive correlations if using the 0.6 Å bandwidth. We demonstrate that rotationally modulated signals observed in SCa II, which appear flat or noisy when using 1.6 Å on SHα, can become more evident if a 0.6 Å bandpass is used instead. Low activity variability appears to be a contributing factor for the cases of weak or no correlations. Conclusions. Calculating the Hα index using a bandpass of 0.6 Å maximises the correlation between Ca II and Hα, both at short and long timescales. On the other hand, the use of the broader 1.6 Å, generally used in exoplanet detection to identify stellar activity signals, degrades the signal by including the flux in the line wings. In view of these results, we strongly recommend the use of a 0.6 Å bandwidth when computing the Hα index for the identification of activity rotational modulation and magnetic cycle signals in solar-type stars.

Optimal Spectral Lines for Measuring Chromospheric Magnetic Fields

This paper identifies spectral lines from X-ray to IR wavelengths which are optimally suited to measuring vector magnetic fields as high as possible in the solar atmosphere. Instrumental and Earth's atmospheric properties, as well as solar abundances, atmospheric properties, and elementary atomic physics are considered without bias toward particular wavelengths or diagnostic techniques. While narrowly focused investigations of individual lines have been reported in detail, no assessment of the comparative merits of all lines has ever been published. Although in the UV, on balance the Mg+ h and k lines near 2800 Å are optimally suited to polarimetry of plasma near the base of the solar corona. This result was unanticipated, given that longer-wavelength lines offer greater sensitivity to the Zeeman effect. While these lines sample optical depths photosphere to the coronal base, we argue that cores of multiple spectral lines provide a far more discriminating probe of magnetic structure as a function of optical depth than the core and inner wings of a strong line. Thus, together with many chromospheric lines of Fe+ between 2585 Å and the h line at 2803 Å, this UV region promises new discoveries concerning how the magnetic fields emerge, heat, and accelerate plasma as they battle to dominate the force and energy balance within the poorly understood chromosphere.

An efficient Monte Carlo model for the slowing down of photoelectrons. Application to H-[formula omitted] in exoplanet atmospheres

Photoelectrons, the fast electrons produced in the photoionization of planetary atmospheres, drive transformations in the atmospheric gas that are often inhibited by energy considerations for thermal electrons. The transformations include excitation and ionization of atoms and molecules, which affect the detectability of these gases and constrain the fraction of incident stellar radiation that transforms into heat. To gain insight into these important questions, we build a Monte Carlo model that solves the slowing down of photoelectrons in a gas with arbitrary amounts of H and He atoms and thermal electrons. Our novel multi-score scheme differs from similar tools in that it efficiently handles rare collisional channels, as in the case of low-abundance excited atoms that undergo superelastic and inelastic collisions. The model is validated and its performance demonstrated. Further, we investigate whether photoelectrons might affect the population of the excited hydrogen H(2) detected at some exoplanet atmospheres by transmission spectroscopy in the H-α line. For the ultra-hot Jupiter HAT-P-32b, we find that photoelectron-driven excitation of H(2) is inefficient at the pressures probed by the line core but becomes significant (yet sub-dominant) deeper in the atmosphere where the line wings form. The contribution of photoelectrons to the destruction of H(2) either by collisional deexcitation or ionization is entirely negligible, a conclusion likely to hold for exoplanet atmospheres at large. Importantly, photoelectrons dominate the gas ionization at the altitudes probed by the H-α line, a fact that will likely affect, even if indirectly, the population of H(2) and other tracers such as metastable helium. Future modeling of these excited levels should incorporate photoelectron-driven ionization.

The water vapor foreign continuum in the 8100-8500 cm−1 spectral range

In the Earth's atmosphere, the foreign absorption continuum of water vapor is due to the interaction of water molecules with other atmospheric gases (mostly N2 and O2). Following our study of the self-continuum in the high energy edge of the 1.25 µm transparency window (Korovela et al, J. Quant. Spectrosc. Radiat. Transfer 286 (2022) 108206), we report here the first room temperature measurements of water vapor foreign continuum in the same spectral region. Foreign continuum cross-sections, CF, are derived at several selected spectral points between 8120 and 8500 cm−1 for humidified nitrogen, humidified oxygen and humidified air (10000 ppm of H2O). The absorption signal is measured by cavity ring-down spectroscopy (CRDS) using pressure ramps up to 750 Torr. While the measured total continuum absorption follows nicely the expected quadratic dependence versus the total pressure, the uncertainty on the retrieved weak foreign continuum is strongly affected by other contributions which have to be subtracted from the measured absorption (far wings of the resonance lines, O2 collision induced absorption, self-continuum, Rayleigh scattering). The obtained H2O-air and H2O-N2CF cross-section values are found to be comparable while the H2O-O2CF values appear to be significantly smaller. Overall, the retrieved CF values for H2O-air mixture validate the MT_CKD model in the considered region.

Detection of Na and K in the Atmosphere of the Hot Jupiter HAT-P-1b with P200/DBSP

We present a new optical transmission spectrum of the hot Jupiter HAT-P-1b based on two transits observed with the Double Spectrograph (DBSP) on the Palomar 200-inch (P200) telescope. The DBSP transmission spectrum, covering a wavelength range from 3250 to 10007 Å, is consistent with that observed with the Hubble Space Telescope (HST), but the former has a finer spectral resolution. The DBSP spectrum alone reveals the presence of a pressure-broadened line wing for Na, the line cores for both Na and K, and tentative evidence for H2O. We obtain consistent results from the spectral retrieval analyses performed on the DBSP-only data set and the DBSP, HST, and Spitzer combined data set. Our retrievals suggest a mostly clear atmosphere for HAT-P-1b, with a cloud coverage of % that is dominated by enhanced haze. We derive subsolar abundances for Na, K, and C, and subsolar to solar for O. Future observations with James Webb Space Telescope and ground-based high-resolution spectrographs should be able to not only confirm the presence of these species but also stringently constrain the formation and migration pathways for HAT-P-1b.

Dark Off-limb Gap: Manifestation of a Temperature Minimum and the Dynamic Nature of the Chromosphere

We study off-limb emission of the lower solar atmosphere using high-resolution imaging spectroscopy in the Hβ and Ca ii 8542 Å lines obtained with the CHROMospheric Imaging Spectrometer (CHROMIS) and the CRisp Imaging SpectroPolarimeter (CRISP) on the Swedish 1 m Solar Telescope. The Hβ line-wing images show the dark intensity gap between the photospheric limb and chromosphere, which is absent in the Ca ii images. We calculate synthetic spectra of the off-limb emissions with the RH code in one-dimensional spherical geometry and find good agreement with the observations. The analysis of synthetic line profiles shows that the gap in the Hβ line-wing images maps the temperature minimum region between the photosphere and chromosphere due to the well-known opacity and emissivity gap of Balmer lines in this layer. However, the observed gap is detected farther from the line core in the outer line-wing positions than in the synthetic profiles. We found that an increased microturbulence in the model chromosphere is needed to reproduce the dark gap in the outer line wing, suggesting that the observed Hβ gap is the manifestation of the temperature minimum and the dynamic nature of the solar chromosphere. The temperature minimum produces a small enhancement in synthetic Ca ii line-wing intensities. Observed off-limb Ca ii line-wing emissions show similar enhancement below the temperature minimum layer near the edge of the photospheric limb.

Discovery of Line Pressure Broadening and Direct Constraint on Gas Surface Density in a Protoplanetary Disk

The gas surface density profile of protoplanetary disks is one of the most fundamental physical properties to understanding planet formation. However, it is challenging to determine the surface density profile observationally, because the H2 emission cannot be observed in low-temperature regions. We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) archival data of the 12CO J = 3 − 2 line toward the protoplanetary disk around TW Hya and discovered extremely broad line wings due to the pressure broadening. In conjunction with a previously reported optically thin CO isotopologue line, the pressure broadened line wings enabled us to directly determine the midplane gas density for the first time. The gas surface density at ∼5 au from the central star reaches ∼103 g cm−2, which suggests that the inner region of the disk has enough mass to form a Jupiter-mass planet. Additionally, the gas surface density drops at the inner cavity by ∼2 orders of magnitude compared to outside the cavity. We also found a low CO abundance of ∼10−6 with respect to H2, even inside the CO snow line, which suggests conversion of CO to less volatile species. Combining our results with previous studies, the gas surface density jumps at r ∼ 20 au, suggesting that the inner region (3 < r < 20 au) might be the magnetorotational instability dead zone. This study sheds light on the direct gas surface density constraint without assuming the CO/H2 ratio using ALMA.

Characterising the orbit and circumstellar environment of the high-mass binary MWC 166 A

Context. Stellar evolution models are highly dependent on accurate mass estimates, especially for highly massive stars in the early stages of stellar evolution. The most direct method for obtaining model-independent stellar masses is derivation from the orbit of close binaries. Aims. Our aim was to derive the first astrometric plus radial velocity orbit solution for the single-lined spectroscopic binary star MWC 166 A, based on near-infrared interferometry over multiple epochs and ∼100 archival radial velocity measurements, and to derive fundamental stellar parameters from this orbit. A supplementary aim was to model the circumstellar activity in the system from K band spectral lines. Methods. The data used include interferometric observations from the VLTI instruments GRAVITY and PIONIER, as well as the MIRC-X instrument at the CHARA Array. We geometrically modelled the dust continuum to derive relative astrometry at 13 epochs, determine the orbital elements, and constrain individual stellar parameters at five different age estimates. We used the continuum models as a base to examine differential phases, visibilities, and closure phases over the Br γ and He I emission lines in order to characterise the nature of the circumstellar emission. Results. Our orbit solution suggests a period of P = 367.7 ± 0.1 d, approximately twice as long as found with previous radial velocity orbit fits. We derive a semi-major axis of 2.61 ± 0.04 au at d = 990 ± 50 pc, an eccentricity of 0.498 ± 0.001, and an orbital inclination of 53.6 ± 0.3°. This allowed the component masses to be constrained to M1 = 12.2 ± 2.2 M⊙ and M2 = 4.9 ± 0.5 M⊙. The line-emitting gas was found to be localised around the primary and is spatially resolved on scales of ∼11 stellar radii, where the spatial displacement between the line wings is consistent with a rotating disc. Conclusions. The large spatial extent and stable rotation axis orientation measured for the Br γ and He I line emission are inconsistent with an origin in magnetospheric accretion or boundary-layer accretion, but indicate an ionised inner gas disc around this Herbig Be star. We observe line variability that could be explained either with generic line variability in a Herbig star disc or V/R variations in a decretion disc scenario. We have also constrained the age of the system, with relative flux ratios suggesting an age of ∼(7 ± 2)×105 yr, consistent with the system being composed of a main-sequence primary and a secondary still contracting towards the main-sequence stage.

Testing Lyα Emission-line Reconstruction Routines at Multiple Velocities in One System

The 1215.67 Å H i Lyα emission line dominates the ultraviolet flux of low-mass stars, including the majority of known exoplanet hosts. Unfortunately, strong attenuation by the interstellar medium (ISM) obscures the line core in most stars, requiring the intrinsic Lyα flux to be reconstructed based on fits to the line wings. We present a test of the widely used Lyα emission-line reconstruction code lyapy using phase-resolved, medium-resolution STIS G140M observations of the close white dwarf–M dwarf binary EG UMa. The Doppler shifts induced by the binary orbital motion move the Lyα emission line in and out of the region of strong ISM attenuation. Reconstructions of each spectrum should produce the same Lyα profile regardless of phase, under the well-justified assumption that there is no intrinsic line variability between observations. Instead, we find that the reconstructions underestimate the Lyα flux by almost a factor of 2 for the lowest velocity, most attenuated spectrum, due to a degeneracy between the intrinsic Lyα and ISM profiles. Our results imply that many stellar Lyα fluxes derived from G140M spectra reported in the literature may be underestimated, with potential consequences for, for example, estimates of extreme-ultraviolet stellar spectra and ultraviolet inputs into simulations of exoplanet atmospheres.

Large Interferometer For Exoplanets (LIFE)

Context.An important future goal in exoplanetology is to detect and characterize potentially habitable planets. Concepts for future space missions have already been proposed: from a large UV-optical-infrared space mission for studies in reflected light, to the Large Interferometer for Exoplanets (LIFE) for analyzing the thermal portion of the planetary spectrum. Using nulling interferometry, LIFE will allow us to constrain the radius and effective temperature of (terrestrial) exoplanets, as well as provide unique information about their atmospheric structure and composition.Aims.We explore the potential of LIFE for characterizing emission spectra of Earth at various stages of its evolution. This allows us (1) to test the robustness of Bayesian atmospheric retrieval frameworks when branching out from a modern Earth scenario while still remaining in the realm of habitable (and inhabited) exoplanets, and (2) to refine the science requirements for LIFE for the detection and characterization of habitable, terrestrial exoplanets.Methods.We performed Bayesian retrievals on simulated spectra of eight different scenarios, which correspond to cloud-free and cloudy spectra of four different epochs of the evolution of the Earth. Assuming a distance of 10 pc and a Sun-like host star, we simulated observations obtained with LIFE using its simulator LIFEsim,considering all major astrophysical noise sources.Results.With the nominal spectral resolution(R= 50) and signal-to-noise ratio (assumed to be S/N = 10 at 11.2 μm), we can identify the main spectral features of all the analyzed scenarios (most notably CO2, H2O, O3, and CH4). This allows us to distinguish between inhabited and lifeless scenarios. Results suggest that O3 and CH4in particular yield an improved abundance estimate by doubling the S/N from 10 to 20. Neglecting clouds in the retrieval still allows for a correct characterization of the atmospheric composition. However, correct cloud modeling is necessary to avoid biases in the retrieval of the correct thermal structure.Conclusions.From this analysis, we conclude that the baseline requirements for R and S/N are sufficient for LIFE to detect O3and CH4in the atmosphere of an Earth-like planet with an O2abundance of around 2% in volume mixing ratio. Doubling the S/N would allow a clearer detection of these species at lower abundances. This information is relevant in terms of the LIFE mission planning. We also conclude that cloud-free retrievals of cloudy planets can be used to characterize the atmospheric composition of terrestrial habitable planets, but not the thermal structure of the atmosphere. From the inter-model comparison performed, we deduce that differences in the opacity tables (caused by, e.g., a different line wing treatment) may be an important source of systematic errors.

The transfer of polarized radiation in resonance lines with partial frequency redistribution, J-state interference, and arbitrary magnetic fields

Aims. We present the theoretical framework and numerical methods we have implemented to solve the problem of the generation and transfer of polarized radiation in spectral lines without assuming local thermodynamical equilibrium, while accounting for scattering polarization, partial frequency redistribution (due to both the Doppler effect and elastic collisions), J-state interference, and hyperfine structure. The resulting radiative transfer code allows one to model the impact of magnetic fields of an arbitrary strength and orientation through the Hanle, incomplete Paschen-Back, and magneto-optical effects. We also evaluate the suitability of a series of approximations for modeling the scattering polarization in the wings of strong resonance lines at a much lower computational cost, which is particularly valuable for the numerically intensive case of three-dimensional radiative transfer. Methods. We examine the suitability of the considered approximations by using our radiative transfer code to model the Stokes profiles of the Mg II h &amp; k lines and of the H I Lyman-α line in magnetized one-dimensional models of the solar atmosphere. Results. Neglecting Doppler redistribution in the scattering processes that are unperturbed by elastic collisions (i.e., treating them as coherent in the observer’s frame) produces a negligible error in the scattering polarization wings of the Mg II resonance lines and a minor one in the Lyman-α wings, although it is unsuitable to model the cores of these lines. For both lines, the scattering processes that are perturbed by elastic collisions only give a significant contribution to the intensity component of the emissivity. Neglecting collisional as well as Doppler redistribution (so that all scattering processes are coherent) represents a rough but suitable approximation for the wings of the Mg II resonance lines, but a very poor one for the Lyman-α wings. The magnetic sensitivity in the scattering polarization wings of the considered lines can be modeled by accounting for the magnetic field in only the ηI and ρV coefficients of the Stokes-vector transfer equation (i.e., using the zero-field expression for the emissivity).

Discovery of optical and infrared accretion disc wind signatures in the black hole candidate MAXI J1348–630

MAXI J1348–630 is a low mass X-ray binary discovered in 2019 during a bright outburst. During this event, the system experienced both hard and soft states following the standard evolution. We present multi-epoch optical and near-infrared spectroscopy obtained with X-shooter at the Very Large Telescope. Our dataset includes spectra taken during the brightest phases of the outburst as well as the decay towards quiescence. We study the evolution of the main emission lines, paying special attention to the presence of features commonly associated with accretion disc winds, such as blueshifted absorptions, broad emission line wings and flat-top profiles. We find broad emission line wings in Hα during the hard-to-soft transition and blueshifted absorption troughs at ∼ − 500 km s−1 in Hβ, He I–5876, Hα and Paβ during the bright soft-intermediate state. In addition, flat-top profiles are seen throughout the outburst. We interpret these observables as signatures of a cold (i.e., optical-to-infrared) accretion disc wind present in the system. We discuss the properties of the wind and compare them with those seen in other X-ray transients. In particular, the wind velocity that we observe is low when compared to those of other systems, which might be a direct consequence of the relatively low binary inclination, as suggested by several observables. This study strengthens the hypothesis that cold winds are a common feature in low mass X-ray binaries and that they can also be detected in low inclination objects via high-quality optical and infrared spectroscopy.

COCOPLOT: COlor COllapsed PLOTting software Using colour to view 3D data as a 2D image

ABSTRACT Most modern solar observatories deliver data products formatted as 3D spatiotemporal data cubes, that contain additional, higher dimensions with spectral and/or polarimetric information. This multidimensional complexity presents a major challenge when browsing for features of interest in several dimensions simultaneously. We developed the COlor COllapsed PLOTting (COCOPLOT) software as a quick-look and context image software, to convey spectral profile or time evolution from all the spatial pixels (x, y) in a 3D [nx, ny, nλ] or [nx, ny, nt] data cube as a single image, using colour. This can avoid the need to scan through many wavelengths, creating difference and composite images when searching for signals satisfying multiple criteria. Filters are generated for the red, green, and blue channels by selecting values of interest to highlight in each channel, and their weightings. These filters are combined with the data cube over the third dimension axis to produce an nx × ny × 3 cube displayed as one true colour image. Some use cases are presented for data from the Swedish 1-m Solar Telescope and Interface Region Imaging Spectrograph, including H α solar flare data, a comparison with k-means clustering for identifying asymmetries in the Ca ii K line and off-limb coronal rain in IRIS C ii slit-jaw images. These illustrate identification by colour alone using COCOPLOT of locations including line wing or central enhancement, broadening, wing absorption, and sites with intermittent flows or time-persistent features. COCOPLOT is publicly available in both idl and python.

The MOSDEF-LRIS survey: connection between galactic-scale outflows and the properties of z ∼ 2 star-forming galaxies

ABSTRACT We investigate the conditions that facilitate galactic-scale outflows using a sample of 155 typical star-forming galaxies at z ∼ 2 drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. The sample includes deep rest-frame UV spectroscopy from the Keck Low-Resolution Imaging Spectrometer (LRIS), which provides spectral coverage of several low-ionization interstellar (LIS) metal absorption lines and Lyα emission. Outflow velocities are calculated from the centroids of the LIS absorption and/or Lyα emission, as well as the highest velocity component of the outflow from the blue wings of the LIS absorption lines. Outflow velocities are found to be marginally correlated or independent of galaxy properties, such as star-formation rate (SFR) and star-formation rate surface density (ΣSFR). Outflow velocity scales with SFR as a power-law with index 0.24, which suggests that the outflows may be primarily driven by mechanical energy generated by supernovae explosions, as opposed to radiation pressure acting on dusty material. On the other hand, outflow velocity and ΣSFR are not significantly correlated, which may be due to the limited dynamic range of ΣSFR probed by our sample. The relationship between outflow velocity and ΣSFR normalized by stellar mass (ΣsSFR), as a proxy for gravitational potential, suggests that strong outflows (e.g. &amp;gt; 200 km s−1) become common above a threshold of log(ΣsSFR/$\rm {yr}^{-1}\ \rm {kpc}^{-2}$) ∼ −11.3, and that above this threshold, outflow velocity uncouples from ΣsSFR. These results highlight the need for higher resolution spectroscopic data and spatially resolved imaging to test the driving mechanisms of outflows predicted by theory.

Spectrum synthesis for radially pulsating stars with shocked atmospheres

ABSTRACT spec_puls describes a suite of computer programs to simulate the emergent spectrum from a radially pulsating star. It combines a Christy-type non-linear pulsation code with classical stellar atmosphere codes. The principal aim is to interpret the dynamical spectrum of the radially pulsating extreme helium star V652 Herculis (V652 Her), which shows a strong shock at minimum radius. The components are general enough to treat other classes of radial pulsation. The theoretical spectrum from a shocked pulsation model shows line doubling, with the blue component emerging at standstill velocity and accelerating blueward. The doubling phase depends on line depth and parent ion. The behaviour of line cores post-shock points to a drop in the ionization temperature, although the gas temperature in the model remains high. Shock compression leads to phase-dependent strengthening of Stark-broadened line wings, with the far wings responding first. With velocity-, temperature-, and pressure-sensitive diagnostics, detailed tomography of the pulsation-driven shock in V652 Her seems possible. Even when no shock is present, the dynamical spectrum is significantly different from a model in hydrostatic equilibrium. Using the quasi-static approximation (e.g. at maximum radius) may lead to a considerable underestimate of the star’s mean effective temperature and surface gravity.

A New Method for Direct Measurement of Isotopologue Ratios in Protoplanetary Disks: A Case Study of the 12CO/13CO Ratio in the TW Hya Disk

Planetary systems are thought to be born in protoplanetary disks. Isotope ratios are a powerful tool for investigating the material origin and evolution from molecular clouds to planetary systems via protoplanetary disks. However, it is challenging to measure the isotope (isotopologue) ratios, especially in protoplanetary disks, because the emission lines of major species are saturated. We developed a new method to overcome these challenges by using optically thin line wings induced by thermal broadening. As a first application of the method, we analyzed two carbon monoxide isotopologue lines, 12CO 3–2 and 13CO 3–2, from archival observations of a protoplanetary disk around TW Hya with the Atacama Large Millimeter/submillimeter Array. The 12CO/13CO ratio was estimated to be 21 ± 5 at disk radii of 70–110 au, which is significantly smaller than the value observed in the local interstellar medium, ∼69. It implies that an isotope exchange reaction occurs in a low-temperature environment with C/O > 1. In contrast, it is suggested that 12CO/13CO is higher than ∼84 in the outer disk (r > 130 au), which can be explained by the difference in the binding energy of the isotopologues on dust grains and the CO gas depletion processes. Our results imply that the gas-phase 12CO/13CO can vary by a factor of >4 even inside a protoplanetary disk and therefore can be used to trace material evolution in disks.

Absorption of methane broadened by carbon dioxide in the 3.3 μm spectral region: From line centers to the far wings

This work studied the infrared absorption of methane broadened by carbon dioxide, which can contribute to the radiative budget of CO2-dominated atmospheres. Fourier transform absorption spectra of CH4 perturbed by CO2 were recorded in the 3.3 μm spectral region, at room temperature and total pressures ranging from 3 to 25 bars. These experimental data were modeled using a theoretical approach taking collision-induced line mixing into account. Comparisons between measured and calculated spectra demonstrate that the proposed model is capable of accurately representing the absorption of methane broadened by collisions with CO2, from line centers to the far wings. For practical applications, this rigorous spectral modeling was used to derive a simple χ-factor model to represent the spectral shape of CO2-broadened CH4 line wings. Comparisons with experimental values show that, at room temperature, the proposed χ-factor reproduces the measured spectra to within 10% in the band wings where the absorption is mainly due to the far line wings, while the usual Lorentzian profile leads to relative differences several orders of magnitude larger. This line-shape correction was successfully validated through comparisons with heavily saturated spectra of the ν3 band of methane recorded at sub-atmospheric pressures.

The atmospheric continuum in the “terahertz gap” region (15–700 cm−1): Review of experiments at SOLEIL synchrotron and modeling

The results of our recent water vapor related continuum studies in the range of the “terahertz gap” at the AILES beam line of the SOLEIL synchrotron using a Fourier transform spectrometer are reviewed and summarized. The continuum of pure water vapor and its mixture with N2, O2 and dry air were investigated in the 15–700 cm−1 range in conditions close to the atmospheric ones. Analysis of the frequency and temperature dependence of the self-continuum revealed a significant contribution of a stable water dimer to the observed absorption, which is not taken into account by the MT_CKD semi-empirical model of the continuum widely used for atmospheric applications. We propose a physically based approach to continuum modeling, which explicitly accounts for the contribution of dimers and far wings of monomer resonance lines. We applied the suggested approach for the water vapor self-continuum modeling in the H2O rotational band and supported it by comparing with the SOLEIL data. The frequency dependence of the relative contributions of far wings, stable and metastable dimers to the self-continuum is evaluated. Prospects for further developments of the atmospheric continuum modeling are considered.

The Atmospheric Response to High Nonthermal Electron-beam Fluxes in Solar Flares. II. Hydrogen-broadening Predictions for Solar Flare Observations with the Daniel K. Inouye Solar Telescope

Abstract Redshifted components of chromospheric emission lines in the hard X-ray impulsive phase of solar flares have recently been studied through their 30 s evolution with the high resolution of the Interface Region Imaging Spectrograph. Radiative-hydrodynamic flare models show that these redshifts are generally reproduced by electron-beam-generated chromospheric condensations. The models produce large ambient electron densities, and the pressure broadening of the hydrogen Balmer series should be readily detected in observations. To accurately interpret the upcoming spectral data of flares with the DKIST, we incorporate nonideal, nonadiabatic line-broadening profiles of hydrogen into the RADYN code. These improvements allow time-dependent predictions for the extreme Balmer line wing enhancements in solar flares. We study two chromospheric condensation models, which cover a range of electron-beam fluxes (1 − 5 × 1011 erg s−1 cm−2) and ambient electron densities (1 − 60 × 1013 cm−3) in the flare chromosphere. Both models produce broadening and redshift variations within 10 s of the onset of beam heating. In the chromospheric condensations, there is enhanced spectral broadening due to large optical depths at Hα, Hβ, and Hγ, while the much lower optical depth of the Balmer series H12−H16 provides a translucent window into the smaller electron densities in the beam-heated layers below the condensation. The wavelength ranges of typical DKIST/ViSP spectra of solar flares will be sufficient to test the predictions of extreme hydrogen wing broadening and accurately constrain large densities in chromospheric condensations.

APEX and NOEMA observations of H2S in nearby luminous galaxies and the ULIRG Mrk 231

Context. In order to understand the evolution and feedback of active galactic nuclei (AGN) and star formation, it is important to use molecular lines as probes of physical conditions and chemistry. Aims. We use H2S to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium in luminous infrared galaxies. Specifically, our aim is to search for evidence of shock enhancement of H2S related to galactic-scale mechanical feedback processes such as outflows. Methods. Using the APEX single-dish telescope, we have observed the 110–101 transition of ortho-H2S at 168 GHz towards the centres of 12 nearby luminous infrared galaxies. We have also observed the same line towards the ultra-luminous infrared galaxy Mrk 231 with the NOEMA interferometer. Results. We detected H2S towards NGC 253, NGC 1068, NGC 3256, NGC 4418, NGC 4826, NGC 4945, Circinus, M 83, and Mrk 231. Upper limits were obtained for NGC 1097, NGC 1377, and IC 860. We also detected line emission from HCN 2–1 in all galaxies in the APEX survey as well as HCO+, HNC, CH3CN, CH3OH, H2CS, HOC+, and SO in several of the sample galaxies. Mrk 231 has a rich 2 mm molecular spectrum and, in addition to H2S, we detect emission from HC3N, CH3OH, HC18O+, C2S, and CH3CCH. Four galaxies show elevated H2S emission relative to HCN: Circinus, NGC 3256, NGC 4826, and NGC 4418. We suggest that the high line ratios are caused by elevated H2S abundances in the dense gas. However, we do not find any clear connection between the H2S/HCN line intensity ratio and the presence (or speed) of molecular outflows in the sample galaxies. Therefore, H2S abundances do not seem to be globally affected by the large-scale outflows. In addition, the H2S/HCN line ratio is not enhanced in the line wings compared to the line core in Mrk 231. This suggests that H2S abundances do not increase in the dense gas in the outflow. However, we do find that the H2S and HCN luminosities (LH2S and LHCN) correlate well with the total molecular gas mass in the outflow, Moutflow(H2), in contrast to LCO and LHCO+. We also find that the line luminosity of H2S correlates with the total infrared luminosity in a similar way as that of H2O. Conclusions. We do not find any evidence of H2S abundance enhancements in the dense gas due to galactic-scale outflows in our sample galaxies, nor in the high-resolution study of Mrk 231. We discuss possible mechanisms behind the suggested H2S abundance enhancements in NGC 4418, Circinus, NGC 3256, and NGC 4826. These include radiative processes (for example X-rays or cosmic rays) or smaller-scale shocks. Further high-resolution and multi-transition studies are required to determine the cause behind the elevated H2S emission in these galaxies. We suggest that LH2S serves as a tracer of the dense gas content, similar to LHCN, and that the correlation between LH2S and Moutflow(H2) implies a relation between the dense gas reservoir and the properties and evolution of the molecular feedback. This potential link requires further study since it holds important keys to our understanding of how the properties of molecular outflows relate to those of their host galaxies. Finally, the similar infrared-correlation coefficients between H2S and H2O may indicate that they originate in the same regions in the galaxy: warm gas in shocks or irradiated by star formation or an AGN.