Showing 1–35 of 35 results for author: Mellah, I E

Whispering in the dark: faint X-ray emission from black holes with OB star companions

Authors: Koushik Sen, Ileyk El Mellah, Norbert Langer, Xiao-Tian Xu, Martin Quast, Daniel Pauli

Abstract: Context. Recent astrometric and spectroscopic surveys of OB stars have revealed a few stellar-mass black holes (BHs) with orbital periods as low as 10 days. No X-ray counterpart has been detected, due to the absence of a radiatively efficient accretion disk around the BH. Yet, dissipative processes in the hot, dilute and strongly magnetized plasma around the BH (so-called BH corona) can still lead… ▽ More Context. Recent astrometric and spectroscopic surveys of OB stars have revealed a few stellar-mass black holes (BHs) with orbital periods as low as 10 days. No X-ray counterpart has been detected, due to the absence of a radiatively efficient accretion disk around the BH. Yet, dissipative processes in the hot, dilute and strongly magnetized plasma around the BH (so-called BH corona) can still lead to non-thermal X-ray emission (e.g. synchrotron). Aims. We determine the X-ray luminosity distribution from BH+OB star binaries up to orbital periods of a few thousand days. Methods. We use detailed binary evolution models computed with MESA for initial primary masses of 10-90 $M_{\odot}$ and orbital periods from 1-3000 d. The X-ray luminosity is computed for a broad range of radiative efficiencies. Results. We show that particle acceleration through magnetic reconnection can heat the BH corona. A substantial fraction of the gravitational potential energy from the accreted plasma is converted into non-thermal X-ray emission. Our population synthesis analysis predicts at least 28 (up to 72) BH+OB star binaries in the Large Magellanic Cloud (LMC) to produce X-ray luminosity above 10$^{31}$ erg$\,$s$^{-1}$, observable through focused Chandra observations. We identify a population of SB1 systems in the LMC and HD96670 in the Milky Way comprising O stars with unseen companions of masses above 2.3 $M_{\odot}$ that aligns well with our predictions. The predicted luminosities of the OB companions to these X-ray-emitting BHs are 10$^{4.5-5.5}$ $L_{\odot}$. Conclusions. These results make the case for long-time exposure in X-rays of the stellar-mass BH candidates identified around OB stars. It will constrain the underlying population of X-ray-faint BHs, the evolution from single to double degenerate binaries, and the progenitors of gravitational wave mergers. (Abridged) △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: Submitted to A&A on May 31, 2024. Comments and feedback are very welcome

ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA

Authors: S. H. J. Wallstrom, T. Danilovich, H. S. P. Muller, C. A. Gottlieb, S. Maes, M. Van de Sande, L. Decin, A. M. S. Richards, A. Baudry, J. Bolte, T. Ceulemans, F. De Ceuster, A. de Koter, I. El Mellah, M. Esseldeurs, S. Etoka, D. Gobrecht, E. Gottlieb, M. Gray, F. Herpin, M. Jeste, D. Kee, P. Kervella, T. Khouri, E. Lagadec , et al. (13 additional authors not shown)

Abstract: The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but… ▽ More The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but a comparative study of the intricate interplay between chemistry and physics is still difficult because observational details such as frequencies and angular resolutions are rarely comparable. Aiming to overcome these deficiencies, ATOMIUM is an ALMA Large Programme to study the physics and chemistry of the circumstellar envelopes of a diverse set of oxygen-rich evolved stars under homogeneous observing conditions at three angular resolutions between ~0.02"-1.4". Here we summarize the molecular inventory of these sources, and the correlations between stellar parameters and molecular content. Seventeen oxygen-rich or S-type asymptotic giant branch (AGB) and red supergiant (RSG) stars have been observed in several tunings with ALMA Band 6, targeting a range of molecules to probe the circumstellar envelope and especially the chemistry of dust formation close to the star. We systematically assigned the molecular carriers of the spectral lines and measured their spectroscopic parameters and the angular extent of the emission of each line from integrated intensity maps. Across the ATOMIUM sample, we detect 291 transitions of 24 different molecules and their isotopologues. This includes several first detections in oxygen-rich AGB/RSG stars: PO v=1, SO2 v1=1 and v2=2, and several high energy H2O transitions. We also find several first detections in S-type AGB stars: vibrationally excited HCN v2=2,3 and SiS v=4,5,6, as well as first detections of the molecules SiC, AlCl, and AlF in W Aql... △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages plus appendices, forthcoming publication in Astronomy & Astrophysics

The first mm detection of a neutron star high-mass X-ray binary

Authors: J. van den Eijnden, L. Sidoli, M. Diaz Trigo, N. Degenaar, I. El Mellah, F. Fürst, V. Grinberg, P. Kretschmar, S. Martínez-Núñez, J. C. A. Miller-Jones, K. Postnov, T. D. Russell

Abstract: Neutron stars accreting from OB supergiants are often divided between persistently and transiently accreting systems, called Supergiant X-ray Binaries (SgXBs) and Supergiant Fast X-ray Transients (SFXTs). This dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, binary orbit, or magnetic field configuration, but direct observational eviden… ▽ More Neutron stars accreting from OB supergiants are often divided between persistently and transiently accreting systems, called Supergiant X-ray Binaries (SgXBs) and Supergiant Fast X-ray Transients (SFXTs). This dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, binary orbit, or magnetic field configuration, but direct observational evidence for these hypotheses remains sparse. To investigate their stellar winds, we present the results of pilot 100-GHz observations of one SFXT and one SgXB with the Northern Extended Millimetre Array. The SFXT, IGR J18410-0535, is detected as a point source at $63.4 \pm 9.6$ $μ$Jy, while the SgXB, IGR J18410-0535 remains undetected. Radio observations of IGR J18410-0535 imply a flat or inverted low-frequency spectrum, arguing for wind emission and against non-thermal flaring. Due to the uncertain SFXT distance, however, the observations do not necessarily imply a difference between the wind properties of the SFXT and SgXB. We compare the mm constraints with other HMXBs and isolated OB supergiants, before considering how future mm campaigns can constrain HMXB wind properties by including X-ray measurements. Specifically, we discuss caveats and future steps to successfully measure wind mass loss rates and velocities in HMXBs with coordinated mm, radio, and X-ray campaigns. △ Less

Submitted 7 September, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS Letters. 7 pages, 1 figure

Reconnection-driven flares in 3D black hole magnetospheres -- A scenario for hot spots around Sagittarius A*

Authors: I. El Mellah, B. Cerutti, B. Crinquand

Abstract: Low-luminosity supermassive and stellar-mass black holes (BHs) may be embedded in a collisionless and highly magnetized plasma. They show non-thermal flares indicative of efficient dissipative processes in the vicinity of the BH. During NIR flares from the supermassive BH Sagittarius A* (Sgr A*), GRAVITY detected circular motion and polarization evolution which suggest the presence of transient sy… ▽ More Low-luminosity supermassive and stellar-mass black holes (BHs) may be embedded in a collisionless and highly magnetized plasma. They show non-thermal flares indicative of efficient dissipative processes in the vicinity of the BH. During NIR flares from the supermassive BH Sagittarius A* (Sgr A*), GRAVITY detected circular motion and polarization evolution which suggest the presence of transient synchrotron-emitting hot spots moving around the BH. We study 3D reconnecting current layers in the magnetosphere of spinning BHs to determine whether plasma-loaded flux ropes formed near the event horizon could reproduce the hot spot observations and help constraining the BH spin. We perform global 3D particle-in-cell simulations in Kerr spacetime of a pair plasma embedded in a strong large-scale magnetic field originating in a disk in prograde Keplerian rotation. A cone-shaped current layer develops which surrounds the twisted open magnetic field lines threading the event horizon. Magnetic field lines coupling the disk to the BH inflate and reconnect a few gravitational radii above the disk. Particles accelerate and accumulate in a few rotating macroscopic flux ropes. Once flux ropes detach, they propagate in the current layer following what appears as a rapidly opening spiral when seen face-on. A single flux rope carries enough relativistic particles to emit synchrotron radiation at levels suitable to reproduce the flares' peak-luminosity of Sgr A* but it quickly fades away as it flows away. Our kinematic analysis of flux ropes' motion favors a BH spin of 0.65 to 0.8 for Sgr A*. The flares' duration require that the underlying magnetic loop seeded in the disk mid-plane has a finite lifetime and azimuthal extent. In this scenario, the hot spot corresponds to a spinning arc along which multiple reconnection sites power the net emission as flux ropes episodically detach. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Submitted to A&A, comments welcome

Journal ref: A&A 677, A67 (2023)

Stellar wind variability in Cygnus X-1 from high-resolution excess variance spectroscopy with Chandra

Authors: Lucia K. Härer, Michael L. Parker, Ileyk El Mellah, Victoria Grinberg, Ralf Ballhausen, Zsofi Igo, Amy Joyce, Jörn Wilms

Abstract: Stellar winds of massive stars are known to be driven by line absorption of UV photons, a mechanism which is prone to instabilities, causing the wind to be clumpy. The clumpy structure hampers wind mass-loss estimates, limiting our understanding of massive star evolution. The wind structure also impacts accretion in high-mass X-ray binary (HMXB) systems. We analyse the wavelength-dependent variabi… ▽ More Stellar winds of massive stars are known to be driven by line absorption of UV photons, a mechanism which is prone to instabilities, causing the wind to be clumpy. The clumpy structure hampers wind mass-loss estimates, limiting our understanding of massive star evolution. The wind structure also impacts accretion in high-mass X-ray binary (HMXB) systems. We analyse the wavelength-dependent variability of X-ray absorption in the wind to study its structure. Such an approach is possible in HMXBs, where the compact object serves as an X-ray backlight. We probe different parts of the wind by analysing data taken at superior and inferior conjunction. We apply excess variance spectroscopy to study the wavelength-dependent soft X-ray variability of the HMXB Cygnus X-1 in the low/hard spectral state. Excess variance spectroscopy quantifies the variability of an object above the statistical noise as a function of wavelength, which allows us to study the variability of individual spectral lines. As one of the first studies, we apply this technique to high-resolution gratings spectra provided by Chandra, accounting for various systematic effects. The frequency dependence is investigated by changing the time binning. The strong orbital phase dependence we observe in the excess variance is consistent with column density variations predicted by a simple model for a clumpy wind. We identify spikes of increased variability with spectral features found by previous spectroscopic analyses of the same data set, most notably from silicon in over-dense clumps in the wind. In the silicon line region, the variability power is redistributed towards lower frequencies, hinting at increased line variability in large clumps. In prospect of the microcalorimetry missions that are scheduled to launch within the next decade, excess variance spectra present a promising approach to constrain the wind structure. △ Less

Submitted 21 September, 2023; v1 submitted 27 April, 2023; originally announced April 2023.

Comments: Accepted by A&A

Journal ref: A&A 680, A72 (2023)

Observing the onset of the accretion wake in Vela X-1

Authors: C. M. Diez, V. Grinberg, F. Fürst, I. El Mellah, M. Zhou, A. Santangelo, S. Martínez-Núñez, R. Amato, N. Hell, P. Kretschmar

Abstract: High-Mass X-ray Binaries (HMXBs) offer a unique opportunity for the investigation of accretion onto compact objects and of wind structure in massive stars. A key source for such studies is the bright neutron star HMXB Vela X-1 whose convenient physical and orbital parameters facilitate the analysis and in particular enable studies of the wind structure in HMXBs. Here, we analyse simultaneous XMM-N… ▽ More High-Mass X-ray Binaries (HMXBs) offer a unique opportunity for the investigation of accretion onto compact objects and of wind structure in massive stars. A key source for such studies is the bright neutron star HMXB Vela X-1 whose convenient physical and orbital parameters facilitate the analysis and in particular enable studies of the wind structure in HMXBs. Here, we analyse simultaneous XMM-Newton and NuSTAR observations at $φ_{\mathrm{orb}} \approx$ 0.36-0.52 and perform time-resolved spectral analysis down to the pulse period of the neutron star, based on our previous NuSTAR-only results. For the first time, we are able to trace the onset of the wakes in a broad 0.5-78 keV range with a high-time resolution of $\sim$283 s and compare to theoretical predictions. We observe a clear rise of the absorption column density of the stellar wind $N_{\mathrm{H,1}}$ starting at orbital phase $\sim$0.44, corresponding to the wake structure entering our line of sight towards the neutron star, together with local extrema throughout the observation possibly associated with clumps or other structures in the wind. Periods of high absorption reveal the presence of multiple fluorescent emission lines of highly ionised species, mainly in the soft X-ray band between 0.5 and 4 keV, indicating photoionisation of the wind. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 18 pages, 18 figures

Journal ref: A&A 674, A147 (2023)

Clumpy wind studies and the non-detection of cyclotron line in OAO 1657-415

Authors: Pragati Pradhan, Carlo Ferrigno, Biswajit Paul, Enrico Bozzo, Ileyk El Mellah, David P. Huenemoerder, James F. Steiner, Victoria Grinberg, Felix Furst, Chandreyee Maitra, Patrizia Romano, Peter Kretschma, Jamie Kennea, Deepto Chakrabarty

Abstract: Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass loss rates. In High Mass X-ray Binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform spectro-timing analysis of the HMXB OAO 1657-415 with non-simultaneous NuSTAR and NICER observations. We compute the… ▽ More Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass loss rates. In High Mass X-ray Binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform spectro-timing analysis of the HMXB OAO 1657-415 with non-simultaneous NuSTAR and NICER observations. We compute the hardness ratio from the energy-resolved light curves, and using an adaptive rebinning technique, we thus select appropriate time segments to search for rapid spectral variations on timescales of a few hundreds to thousands of seconds. Column density and intensity of Iron K$α$ line were strongly correlated, and the recorded spectral variations were consistent with accretion from a clumpy wind. We also illustrate a novel framework to measure clump sizes, masses in HMXBs more accurately based on absorption measurements and orbital parameters of the source. We then discuss the limitations posed by current X-ray spacecrafts in such measurements and present prospects with future X-ray missions. We find that the source pulse profiles show a moderate dependence on energy. We identify a previously undetected dip in the pulse profile visible throughout the NuSTAR observation near spin phase 0.15 possibly caused by intrinsic changes in accretion geometry close to the neutron star. We do not find any evidence for the debated cyclotron line at $\sim$ 36\,keV in the time-averaged or the phase-resolved spectra with NuSTAR. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 32 pages, 11 figures in main text, 7 figures in Appendix, Accepted for publication in ApJ

Magnetic reconnection plasmoid model for Sagittarius A* flares

Authors: N. Aimar, A. Dmytriiev, F. H. Vincent, I. El Mellah, T. Paumard, G. Perrin, A. Zech

Abstract: Sagittarius A*, the supermassive black hole at the center of our galaxy, exhibits episodic near-infrared flares. The recent monitoring of three such events by the GRAVITY instrument has shown that some flares are associated with orbital motions in the close environment of the black hole with super Keplerian velocity. We develop a semi-analytic model of Sagittarius~A* flares based on an ejected lar… ▽ More Sagittarius A*, the supermassive black hole at the center of our galaxy, exhibits episodic near-infrared flares. The recent monitoring of three such events by the GRAVITY instrument has shown that some flares are associated with orbital motions in the close environment of the black hole with super Keplerian velocity. We develop a semi-analytic model of Sagittarius~A* flares based on an ejected large plasmoid, inspired by recent particle-in-cell global simulations of black hole magnetospheres. We model the infrared astrometric and photometric signatures associated to this model. We consider a spherical large plasmoid ejected along a conical orbit around the black hole. This plasmoid is assumed to be formed by successive mergers of smaller plasmoids produced through magnetic reconnection. Non-thermal electrons are injected in the plasmoid. We compute the evolution of the electron-distribution under the influence of synchrotron cooling. We solve the radiative transfer problem and transport the radiation along null geodesics of the Schwarzschild spacetime. We also take into account the quiescent radiation of the accretion flow, on top of which the flare evolves. For the first time, we successfully account for the astrometric and flux variations of the GRAVITY data with a flare model that incorporates an explicit modeling of the emission mechanism. We find good agreement between the prediction of our model and the recent data. In particular, the azimuthal velocity is set by the magnetic field line it belongs to, which is anchored in the inner parts of the accretion flow, hence the super-Keplerian motion. The astrometric track is also shifted with respect to the center of mass due to the quiescent radiation, in agreement with the difference measured with the GRAVITY data. These results support the picture of magnetic reconnection as a viable model for Sagittarius~A* infrared flares. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 20 pages, 14 figures, accepted in A&A

Journal ref: A&A 672, A62 (2023)

The VLT/SPHERE view of the ATOMIUM cool evolved star sample. I. Overview: Sample characterization through polarization analysis

Authors: M. Montargès, E. Cannon, A. de Koter, T. Khouri, E. Lagadec, P. Kervella, L. Decin, I. McDonald, W. Homan, L. B. F. M. Waters, R. Sahai, C. A. Gottlieb, J. Malfait, S. Maes, B. Pimpanuwat, M. Jeste, T. Danilovich, F. De Ceuster, M. Van de Sande, D. Gobrecht, S. H. J. Wallström, K. T. Wong, I. El Mellah, J. Bolte, F. Herpin , et al. (10 additional authors not shown)

Abstract: Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution dat… ▽ More Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust. Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $π^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 22 pages, 15 figures, 5 tables

Journal ref: A&A 671, A96 (2023)

First 3D Radiation-Hydrodynamic Simulations of Wolf-Rayet Winds

Authors: N. Moens, L. G. Poniatowski, L. Hennicker, J. O. Sundqvist, I. El Mellah, N. D. Kee

Abstract: Classical Wolf Rayet (WR) stars are direct supernova progenitors undergoing vigorous mass-loss. Understanding the dense and fast outflows of such WR stars is thus crucial for understanding advanced stages of stellar evolution, the dynamical feedback of massive stars on their environments, and characterizing the distribution of black hole masses. In this paper, we develop first time-dependent, mult… ▽ More Classical Wolf Rayet (WR) stars are direct supernova progenitors undergoing vigorous mass-loss. Understanding the dense and fast outflows of such WR stars is thus crucial for understanding advanced stages of stellar evolution, the dynamical feedback of massive stars on their environments, and characterizing the distribution of black hole masses. In this paper, we develop first time-dependent, multi-dimensional, radiation-hydrodynamical models of the extended optically thick atmospheres and wind outflows of hydrogen-free classical WR stars. A flux limiting radiation hydrodynamics approach is used on a finite volume mesh to model WR outflows. The opacities are described using a combination of tabulated Rosseland mean opacities and the enhanced line opacities expected within a supersonic flow. For high-luminosity models, a radiation-driven, dense, supersonic wind is launched from deep sub-surface regions associated with peaks in the Rosseland mean opacity. For a model with lower luminosity, on the other hand, the Rosseland mean opacity is not sufficient to sustain a net-radial outflow in the sub-surface regions. Rather, what develops in this case is a "standard" line-driven wind launched from the optically thin regions above an extended and highly turbulent atmosphere. We thus find here a natural transition from optically thick outflows of classical WR stars to optically thin winds of hot, compact sub-dwarfs; in our simulations this transition occurs approximately at a luminosity that is about 40% of the Eddington luminosity. Because of the changing character of the wind-launching mechanism, this transition is also accompanied by a large drop (on the low-luminosity end) in average mass-loss rate. △ Less

Submitted 8 September, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 22 pages, 22 figures, Accepted and in publication by A&A

Journal ref: A&A 665, A42 (2022)

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

Authors: C. A. Gottlieb, L. Decin, A. M. S. Richards, F. De Ceuster, W. Homan, S. H. J. Wallstrom, T. Danilovich, T. J. Millar, M. Montarges, K. T. Wong, I. McDonald, A. Baudry, J. Bolte, E. Cannon, E. De Beck, A. de Koter, I. El Mellah, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, P. Kervella, T. Khouri, E. Lagadec , et al. (11 additional authors not shown)

Abstract: This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical e… ▽ More This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field. △ Less

Submitted 13 December, 2021; v1 submitted 8 December, 2021; originally announced December 2021.

Journal ref: A&A 660, A94 (2022)

Spinning black holes magnetically connected to a Keplerian disk -- Magnetosphere, reconnection sheet, particle acceleration and coronal heating

Authors: I. El Mellah, B. Cerutti, B. Crinquand, K. Parfrey

Abstract: Context: Accreting black holes (BHs) may be surrounded by a highly magnetized plasma threaded by a poloidal magnetic field. Non-thermal flares and high energy components could originate from a hot, collisionless and nearly force-free corona. The jets we often observe from these systems are believed to be rotation-powered and magnetically-driven. Aims: We study axisymmetric BH magnetospheres wher… ▽ More Context: Accreting black holes (BHs) may be surrounded by a highly magnetized plasma threaded by a poloidal magnetic field. Non-thermal flares and high energy components could originate from a hot, collisionless and nearly force-free corona. The jets we often observe from these systems are believed to be rotation-powered and magnetically-driven. Aims: We study axisymmetric BH magnetospheres where some magnetic field lines anchored in a surrounding disk can connect to the event horizon of a rotating BH. We identify the sites of magnetic reconnection within 30 gravitational radii depending on the BH spin. Methods: With the fully general relativistic particle-in-cell code GRZeltron, we solve the time-dependent dynamics of the electron-positron pair plasma and of the electromagnetic fields around the BH. The disk is represented by a steady plasma in Keplerian rotation, threaded by a frozen dipolar field. Results: For prograde disks, twisted open magnetic field lines crossing the horizon power a Blandford-Znajek jet while beyond a critical distance, open field lines on the disk are open. In the innermost regions, coupling field lines ensure the transfer of significant amounts of angular momentum and energy between the BH and the disk. From the Y-point at the intersection, a current sheet forms where particle acceleration via magnetic reconnection takes place. We compute the synchrotron images of the current sheet emission. Conclusions: Our estimates for jet power and BH-disk exchanges match those derived from purely force-free models. Dissipation at the Y-point heats the corona and provides a physically motivated source of hard X-rays above the disk for reflection models. Episodic plasmoid ejection might explain millisecond flares observed in Cyg X-1. Particles flowing from the Y-point down to the disk could produce a hot spot at the footpoint of the outermost closed field line. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: Submitted to A&A, 27 pages

Journal ref: A&A 663, A169 (2022)

Binary-object spectral-synthesis in 3D (BOSS-3D) -- Modelling H-alpha emission in the enigmatic multiple system LB-1

Authors: L. Hennicker, N. D. Kee, T. Shenar, J. Bodensteiner, M. Abdul-Masih, I. El Mellah, H. Sana, J. O. Sundqvist

Abstract: Context: To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perf… ▽ More Context: To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perfect example of such a complex multi-D problem. Aims: To improve our understanding of the LB-1 system, we directly modelled the phase-dependent H-alpha line profiles of this system. To this end, we developed a multi-purpose binary-object spectral-synthesis code in 3D (BOSS-3D). Methods: BOSS-3D calculates synthetic line profiles for a given state of the circumstellar material. The standard pz-geometry commonly used for single stars is extended by defining individual coordinate systems for each involved object and by accounting for the appropriate coordinate transformations. The code is then applied to the LB-1 system, considering two main hypotheses, a binary containing a stripped star and Be star, or a B star and a black hole with a disc. Results: Comparing these two scenarios, neither model can reproduce the detailed phase-dependent shape of the H-alpha line profiles. A satisfactory match with the observations, however, is obtained by invoking a disc around the primary object in addition to the Be-star disc or the black-hole accretion disc. Conclusions: The developed code can be used to model synthetic line profiles for a wide variety of binary systems, ranging from transit spectra of planetary atmospheres, to post-asymptotic giant branch binaries including circumstellar and circumbinary discs and massive-star binaries with stellar winds and disc systems. For the LB-1 system, our modelling provides strong evidence that each object in the system contains a disc-like structure. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Journal ref: A&A 660, A17 (2022)

X-ray emission from BH+O star binaries expected to descend from the observed galactic WR+O binaries

Authors: K. Sen, X. -T. Xu, N. Langer, I. El Mellah, C. Schurmann, M. Quast

Abstract: In the Milky Way, $\sim$18 Wolf-Rayet+O (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O (BH+O) binaries are thought to evolve from the former, only one such system is known in the Milky Way. To resolve this disparity, it was suggested that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted. We rea… ▽ More In the Milky Way, $\sim$18 Wolf-Rayet+O (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O (BH+O) binaries are thought to evolve from the former, only one such system is known in the Milky Way. To resolve this disparity, it was suggested that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted. We reassess this issue, with emphasis on the uncertainty in the formation of an accretion disk around wind-accreting BHs in BH+O binaries, which is key to identifying such systems. We follow the methodology of previous work and apply an improved analytic criterion for the formation of an accretion disk around wind accreting BHs. We then use stellar models to predict the properties of the BH+O binaries which are expected to descend from the observed WR+O binaries, if the WR stars would form BHs without a natal kick. We find that disk formation depends sensitively on the O stars' wind velocity, the specific angular momentum carried by the wind, the efficiency of angular momentum accretion by the BH, and the spin of the BH. We show that the assumption of a low wind velocity may lead to predicting that most of the BH+O star binaries will have an extended X-ray bright period. However, this is not the case when typical wind velocities of O stars are considered. We find that a high spin of the BH can boost the duration of the X-ray active phase as well as the X-ray brightness during this phase, producing a strong bias for detecting high mass BH binaries in X-rays with high BH spin parameters. We conclude that large BH formation kicks are not required to understand the sparsity of X-ray bright BH+O stars in the Milky Way. Probing for a population of X-ray silent BH+O systems with alternative methods can inform us about BH kicks and the conditions for high energy emission from high mass BH binaries. (Abridged) △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 14 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 652, A138 (2021)

Revisiting the archetypical wind accretor Vela X-1 in depth -- A case study of a well-known X-ray binary and the limits of our knowledge

Authors: Peter Kretschmar, Ileyk El Mellah, Silvia Martínez-Núñez, Felix Fürst, Victoria Grinberg, Andreas A. C. Sander, Jakob van den Eijnden, Nathalie Degenaar, Jesús Maíz-Apellániz, Francisco Jiménez Esteban, Mercedes Ramos-Lerate, Enrique Utrilla

Abstract: Context: Vela X-1 is one of the best studied X-ray binaries. Frequently though, specific values for its parameters have been used in subsequent studies without considering alternatives. Aims: We aim to provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies and identify specific avenues of possible future research. Methods: We e… ▽ More Context: Vela X-1 is one of the best studied X-ray binaries. Frequently though, specific values for its parameters have been used in subsequent studies without considering alternatives. Aims: We aim to provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies and identify specific avenues of possible future research. Methods: We explore the literature for Vela X-1 and on modelling efforts, describing the evolution of the system knowledge. We also add information derived from public data, especially the Gaia EDR3 release. Results: We update the distance to Vela X-1, the spectral classification for HD 77518 and find that the supergiant may be very close to filling its Roche lobe. Constraints on the clumpiness of the stellar wind have improved. The orbit is very well determined, but the uncertain inclination limits information on the neutron star mass. Estimates for the stellar wind have evolved towards lower velocities, supporting the idea of transient wind-captured disks around the neutron star. Hydrodynamic models and observations are consistent with an accretion wake trailing the neutron star. Conclusions: Vela X-1 is an excellent laboratory, but a lot of room remains to improve. Well-coordinated multi-wavelength observations and campaigns addressing the intrinsic variability are required. New opportunities will arise through new instrumentation. Models of the stellar wind should account for the orbital eccentricity and the non-spherical shape of HD 77581. Realistic multi-dimensional models of radiative transfer in the UV and X-rays are needed, but remain very challenging. Improved MHD models covering a wide range of scales would be required to improve understanding of the plasma-magnetosphere coupling. A full characterization of the accretion column remains another open challenge. (Abbreviated for arXiv) △ Less

Submitted 11 May, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 50 pages, 20 figures, 12 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 652, A95 (2021)

Atomium: The astounding complexity of the near circumstellar environment of the M-type AGB star R Hydrae. I. Morpho-kinematical interpretation of CO and SiO emission

Authors: Ward Homan, Bannawit Pimpanuwat, Fabrice Herpin, Taissa Danilovich, Iain McDonald, Sofia H. J. Wallström, Anita M. S. Richards, Alain Baudry, Raghvendra Sahai, Tom J. Millar, Alex de Koter, C. A. Gottlieb, Pierre Kervella, Miguel Montargès, Marie Van de Sande, Leen Decin, Albert Zijlstra, Sandra Etoka, Manali Jeste, Holger S. P. Müller, Silke Maes, Jolien Malfait, Karl Menten, John Plane, Kelvin Lee , et al. (14 additional authors not shown)

Abstract: Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we… ▽ More Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 21 pages, 23 figures

Journal ref: A&A 651, A82 (2021)

Radiation-Hydrodynamics with MPI-AMRVAC: Flux-Limited Diffusion

Authors: N. Moens, J. O. Sundqvist, I. El Mellah, L. Poniatowski, J. Teunissen, R. Keppens

Abstract: Radiation controls the dynamics and energetics of many astrophysical environments. To capture the coupling between the radiation and matter, however, is often a physically complex and computationally expensive endeavour. We develop a numerical tool to perform radiation-hydrodynamics simulations in various configurations at an affordable cost. We build upon the finite volume code MPI-AMRVAC to solv… ▽ More Radiation controls the dynamics and energetics of many astrophysical environments. To capture the coupling between the radiation and matter, however, is often a physically complex and computationally expensive endeavour. We develop a numerical tool to perform radiation-hydrodynamics simulations in various configurations at an affordable cost. We build upon the finite volume code MPI-AMRVAC to solve the equations of hydrodynamics on multi-dimensional adaptive meshes and introduce a new module to handle the coupling with radiation. A non-equilibrium, flux-limiting diffusion approximation is used to close the radiation momentum and energy equations. The time-dependent radiation energy equation is then solved within a flexible framework, accounting fully for radiation forces and work terms and further allowing the user to adopt a variety of descriptions for the radiation-matter interaction terms (the 'opacities'). We validate the radiation module on a set of standard testcases for which different terms of the radiative energy equation predominate. As a preliminary application to a scientific case, we calculate spherically symmetric models of the radiation-driven and optically thick supersonic outflows from massive Wolf-Rayet stars. This also demonstrates our code's flexibility, as the illustrated simulation combines opacities typically used in static stellar structure models with a parametrised form for the enhanced line-opacity expected in supersonic flows. This new module provides a convenient and versatile tool to perform multi-dimensional and high resolution radiative-hydrodynamics simulations in optically thick environments with the MPI-AMRVAC code. The code is ready to be used for a variety of astrophysical applications, where a first target for us will be multi-dimensional simulations of stellar outflows from Wolf-Rayet stars. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 14 pages, 10 figures Submitted to A&A

Journal ref: A&A 657, A81 (2022)

Atomium: A high-resolution view on the highly asymmetric wind of the AGB star Pi1 Gruis. I. First detection of a new companion and its effect on the inner wind

Authors: Ward Homan, Miguel Montarges, Bannawit Pimpanuwat, Anita M. S. Richards, Sofia H. J. Wallstrom, Pierre Kervella, Leen Decin, Albert Zijlstra, Taissa Danilovich, Alex de Koter, Karl Menten, Raghvendra Sahai, John Plane, Kelvin Lee, Rens Waters, Alain Baudry, Ka Tat Wong, Tom J. Millar, Marie Van de Sande, Eric Lagadec, David Gobrecht, Jeremy Yates, Daniel Price, Emily Cannon, Jan Bolte , et al. (13 additional authors not shown)

Abstract: The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is tho… ▽ More The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is thought to harbour a second, closer-by (<10 au) companion, was observed with the Atacama Large Millimeter/submillimeter Array as part of the ATOMIUM Large programme. In this work, the brightest CO, SiO, and HCN molecular line transitions are analysed. The continuum map shows two maxima, separated by 0.04'' (6 au). The CO data unambiguously reveal that Pi1 Gru's circumstellar environment harbours an inclined, radially outflowing, equatorial density enhancement. It contains a spiral structure at an angle of 38+/-3 deg with the line-of-sight. The HCN emission in the inner wind reveals a clockwise spiral, with a dynamical crossing time of the spiral arms consistent with a companion at a distance of 0.04'' from the asymptotic giant branch star, which is in agreement with the position of the secondary continuum peak. The inner wind dynamics imply a large acceleration region, consistent with a beta-law power of ~6. The CO emission suggests that the spiral is approximately Archimedean within 5'', beyond which this trend breaks down as the succession of the spiral arms becomes less periodic. The SiO emission at scales smaller than 0.5'' exhibits signatures of gas in rotation, which is found to fit the expected behaviour of gas in the wind-companion interaction zone. An investigation of SiO maser emission reveals what could be a stream of gas accelerating from the surface of the AGB star to the companion. Using these dynamics, we have tentatively derived an upper limit on the companion mass to be ~1.1 Msol. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 19 pages, 20 figures

Journal ref: A&A 644, A61 (2020)

(Sub)stellar companions shape the winds of evolved stars

Authors: L. Decin, M. Montargès, A. M. S. Richards, C. A. Gottlieb, W. Homan, I. McDonald, I. El Mellah, T. Danilovich, S. H. J. Wallström, A. Zijlstra, A. Baudry, J. Bolte, E. Cannon, E. De Beck, F. De Ceuster, A. de Koter, J. De Ridder, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, E. Lagadec, P. Kervella, T. Khouri , et al. (10 additional authors not shown)

Abstract: Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their… ▽ More Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe. △ Less

Submitted 28 September, 2020; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: 19 pages main journal, 97 pages Supplementary Information

Journal ref: Science, 2000, Vol 369, Issue 6510, p. 1497

Radiography in high mass X-ray binaries -- Micro-structure of the stellar wind through variability of the column density

Authors: I. El Mellah, V. Grinberg, J. O. Sundqvist, F. A. Driessen, M. A. Leutenegger

Abstract: In high mass X-ray binaries (HMXBs), an accreting compact object orbits a high mass star which loses mass through a dense and inhomogeneous wind. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass loss rates for high mass stars. W… ▽ More In high mass X-ray binaries (HMXBs), an accreting compact object orbits a high mass star which loses mass through a dense and inhomogeneous wind. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass loss rates for high mass stars. We explore the impact of clumpiness on the variability of the column density with a simplified wind model. In particular, we focus on the standard deviation of the column density and the characteristic duration of enhanced absorption episodes, and compare them with analytical predictions based on the porosity length. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a clump in front of the compact object. We provide a recipe to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. The coherence time scale grants direct access to the size of the clumps while their mass can be deduced separately from the amplitude of the variability. If it is due to unaccreted passing-by clumps, the high column density variations in some HMXBs requires high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are hardly compatible with the ones derived from first principles. Alternatively, other components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism still to be identified, or a dense environment in the immediate vicinity of the accretor such as an accretion disk, an outflow or a spherical shell around the magnetosphere of the accreting neutron star. △ Less

Submitted 10 September, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Journal ref: A&A 644, C1 (2020)

Wind morphology around cool evolved stars in binaries: the case of slowly accelerating oxygen-rich outflows

Authors: I. El Mellah, J. Bolte, L. Decin, W. Homan, R. Keppens

Abstract: The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run… ▽ More The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Journal ref: A&A 637, A91 (2020)

Vela X-1 as a laboratory for accretion in High-Mass X-ray Binaries

Authors: P. Kretschmar, S. Martínez-Núñez, F. Fürst, V. Grinberg, M. Lomaeva, I. El Mellah, A. Manousakis, A. A. C. Sander, N. Degenaar, J. van den Eijnden

Abstract: Vela X-1 is an eclipsing high mass X-ray binary (HMXB) consisting of a 283s accreting X-ray pulsar in a close orbit of 8.964 days around the B0.5Ib supergiant HD77581 at a distance of just 2.4 kpc. The system is considered a prototype of wind-accreting HMXB and it has been used as a baseline in different theoretical or modelling studies. We discuss the observational properties of the system and… ▽ More Vela X-1 is an eclipsing high mass X-ray binary (HMXB) consisting of a 283s accreting X-ray pulsar in a close orbit of 8.964 days around the B0.5Ib supergiant HD77581 at a distance of just 2.4 kpc. The system is considered a prototype of wind-accreting HMXB and it has been used as a baseline in different theoretical or modelling studies. We discuss the observational properties of the system and the use of the observational data as laboratory to test recent developments in modelling the accretion process in High-Mass X-ray Binaries (e.g., Sander et al. 2018; El Mellah et al. 2018), which range from detailed descriptions of the wind acceleration to modelling of the structure of the flow of matter close to the neutron star and its variations. △ Less

Submitted 21 May, 2019; originally announced May 2019.

Comments: 4 pages, 2 figures, proceedings of the 12th INTEGRAL conference "INTEGRAL looks AHEAD to Multimessenger astronomy" in Geneva (Switzerland) on 11-15 February 2019

Variability in high-mass X-ray binaries

Authors: F. Fuerst, P. Kretschmar, V. Grinberg, K. Pottschmidt, J. Wilms, M. Kuehnel, I. El Mellah, S. Martinez-Nunez

Abstract: Strongly magnetized, accreting neutron stars show periodic and aperiodic variability over a wide range of time scales. By obtaining spectral and timing information on these different time scales, we can have a closer look into the physics of accretion close to the neutron star and the properties of the accreted material. One of the most prominent time scales is the strong pulsation, i.e., the rota… ▽ More Strongly magnetized, accreting neutron stars show periodic and aperiodic variability over a wide range of time scales. By obtaining spectral and timing information on these different time scales, we can have a closer look into the physics of accretion close to the neutron star and the properties of the accreted material. One of the most prominent time scales is the strong pulsation, i.e., the rotation period of the neutron star itself. Over one rotation, our view of the accretion column and the X-ray producing region changes significantly. This allows us to sample different physical conditions within the column but at the same time requires that we have viewing-angle-resolved models to properly describe them. In wind-fed high-mass X-ray binaries, the main source of aperiodic variability is the clumpy stellar wind, which leads to changes in the accretion rate (i.e., luminosity) as well as absorption column. This variability allows us to study the behavior of the accretion column as a function of luminosity, as well as to investigate the structure and physical properties of the wind, which we can compare to winds in isolated stars. △ Less

Submitted 7 April, 2019; originally announced April 2019.

Comments: 6 pages, 4 figures, accepted for publication in Astronomische Nachrichten (proceedings of the XMM-Newton Workshop 2019)

Reduction of the maximum mass-loss rate of OH/IR stars due to unnoticed binary interaction

Authors: L. Decin, W. Homan, T. Danilovich, A. de Koter, D. Engels, L. B. F. M. Waters, S. Muller, C. Gielen, D. A. García-Hernández, R. J. Stancliffe, M. Vande Sande, G. Molenberghs, F. Kerschbaum, A. A. Zijlstra, I. El Mellah

Abstract: In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind… ▽ More In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind phase is too short for these stars to become white dwarfs. Here, we report on the detection of spiral structures around two cornerstone extreme OH/IR-stars, OH26.5+0.6 and OH30.1-0.7, identifying them as wide binary systems. Hydrodynamical simulations show that the companion's gravitational attraction creates an equatorial density enhancement mimicking a short extreme superwind phase, thereby solving the decades-old conundrum. This discovery restricts the maximum mass-loss rate of AGB stars around the single-scattering radiation-pressure limit of a few 10^(-5) Msun/yr. This brings about crucial implications for nucleosynthetic yields, planet survival, and the wind-driving mechanism. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: Publication date: 25 February 2019 at 16:00h GMT

Journal ref: Nature Astronomy, 2019

Wind Roche lobe overflow in high mass X-ray binaries : a possible mass transfer mechanism for Ultraluminous X-ray sources

Authors: I. El Mellah, J. O. Sundqvist, R. Keppens

Abstract: Ultra-luminous X-ray sources (ULX) have so high X-ray luminosities that they were long thought to be accreting intermediate mass black holes. Yet, some ULX have been shown to display periodic modulations and coherent pulsations, suggestive of a neutron star in orbit around a stellar companion and accreting at super-Eddington rates. In this letter, we propose that the mass transfer in ULX could be… ▽ More Ultra-luminous X-ray sources (ULX) have so high X-ray luminosities that they were long thought to be accreting intermediate mass black holes. Yet, some ULX have been shown to display periodic modulations and coherent pulsations, suggestive of a neutron star in orbit around a stellar companion and accreting at super-Eddington rates. In this letter, we propose that the mass transfer in ULX could be qualitatively the same as in Supergiant X-ray binaries (SgXB), with a wind from the donor star highly beamed towards the compact object. Since the star does not fill its Roche lobe, this mass transfer mechanism known as "wind Roche lobe overflow" can remain stable even for large mass ratios. Based on realistic acceleration profiles derived from spectral observations and modeling of the stellar wind, we compute the bulk motion of the wind to evaluate the fraction of the stellar mass outflow captured by the compact object. We identify the orbital and stellar conditions for a SgXB to transfer mass at rates matching the expectations for ULX and show that the transition from SgXB to ULX luminosity levels is progressive. These results indicate that a high stellar Roche lobe filling factor is not necessary to funnel large quantities of material into the Roche lobe of the accretor. Large stellar mass loss rates such as the ones from the Wolf-Rayet star in M101 ULX-1 or the late B9 Supergiant in NGC 7793 P13 are enough to lead to a highly beamed wind and a significantly enhanced mass transfer rate. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Journal ref: A&A 622, L3 (2019)

Formation of wind-captured discs in Supergiant X-ray binaries : consequences for Vela X-1 and Cygnus X-1

Authors: I. El Mellah, A. A. C. Sander, J. O. Sundqvist, R. Keppens

Abstract: In Supergiant X-ray binaries (SgXB), a compact object captures a fraction of the wind of an O/B supergiant on a close orbit. Proxies exist to evaluate the efficiency of mass and angular momentum accretion but they depend so dramatically on the wind speed that given the current uncertainties, they only set loose constrains. Furthermore, they often bypass the impact of orbital and shock effects on t… ▽ More In Supergiant X-ray binaries (SgXB), a compact object captures a fraction of the wind of an O/B supergiant on a close orbit. Proxies exist to evaluate the efficiency of mass and angular momentum accretion but they depend so dramatically on the wind speed that given the current uncertainties, they only set loose constrains. Furthermore, they often bypass the impact of orbital and shock effects on the flow structure. We study the wind dynamics and the angular momentum gained as the flow is accreted. We identify the conditions for the formation of a disc-like structure around the accretor and the observational consequences for SgXB. We use recent results on the wind launching mechanism to compute 3D streamlines, accounting for the gravitational and X-ray ionizing influence of the compact companion on the wind. Once the flow enters the Roche lobe of the accretor, we solve the hydrodynamics equations with cooling. A shocked region forms around the accretor as the flow is beamed. For wind speeds of the order of the orbital speed, the shock is highly asymmetric compared to the axisymmetric bow shock obtained for a purely planar homogeneous flow. With net radiative cooling, the flow always circularizes for wind speeds low enough. Although the donor star does not fill its Roche lobe, the wind can be significantly beamed and bent by the orbital effects. The net angular momentum of the accreted flow is then sufficient to form a persistent disc-like structure. This mechanism could explain the proposed limited outer extension of the accretion disc in Cygnus X-1 and suggests the presence of a disc at the outer rim of the neutron star magnetosphere in Vela X-1, with dramatic consequences on the spinning up of the accretor. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Journal ref: A&A 622, A189 (2019)

Accretion form a clumpy massive-star wind in Supergiant X-ray binaries

Authors: I. El Mellah, J. O. Sundqvist, R. Keppens

Abstract: Supergiant X-ray Binaries (SgXB) host a compact object, often a neutron star (NS), orbiting an evolved O/B star. Mass transfer proceeds through the intense line-driven wind of the stellar donor, a fraction of which is captured by the gravitational field of the NS. The subsequent accretion process onto the NS is responsible for the abundant X-ray emission from SgXB. They also display peak-to-peak v… ▽ More Supergiant X-ray Binaries (SgXB) host a compact object, often a neutron star (NS), orbiting an evolved O/B star. Mass transfer proceeds through the intense line-driven wind of the stellar donor, a fraction of which is captured by the gravitational field of the NS. The subsequent accretion process onto the NS is responsible for the abundant X-ray emission from SgXB. They also display peak-to-peak variability of the X-ray flux by a factor of a few 10 to 100, along with changes in the hardness ratios possibly due to varying absorption along the line-of-sight. We use recent radiation-hydrodynamic simulations of inhomogeneities (aka clumps) in the non-stationary wind of massive hot stars to evaluate their impact on the time-variable accretion process. For this, we run 3D hydrodynamic simulations of the wind in the vicinity of the accretor to investigate the formation of the bow shock and follow the inhomogeneous flow over several spatial orders of magnitude, down to the NS magnetosphere. In particular, we show that the impact of the wind clumps on the time-variability of the intrinsic mass accretion rate is severely tempered by the crossing of the shock, compared to the purely ballistic Bondi-Hoyle-Lyttleton estimation. We also account for the variable absorption due to clumps passing by the line-of-sight and estimate the final effective variability of the column density and mass accretion rate for different orbital separations. Finally, we compare our results to the most recent analysis of the X-ray flux and the hardness ratio in Vela X-1. △ Less

Submitted 23 November, 2017; originally announced November 2017.

The clumpy absorber in the high-mass X-ray binary Vela X-1

Authors: V. Grinberg, N. Hell, I. El Mellah, J. Neilsen, A. A. C. Sander, M. Leutenegger, F. Fürst, D. P. Huenemoerder, P. Kretschmar, M. Kühnel, S. Martínez-Núñez, S. Niu, K. Pottschmidt, N. S. Schulz, J. Wilms, M. A. Nowak

Abstract: Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion… ▽ More Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium and neon whose strengths and presence depend on the overall level of absorption. They imply a co-existence of cool and hot gas phases in the system that we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries. △ Less

Submitted 17 November, 2017; originally announced November 2017.

Comments: 17 pages, 14 figures; A&A accepted

Journal ref: A&A 608, A143 (2017)

MPI-AMRVAC 2.0 for Solar and Astrophysical Applications

Authors: C. Xia, J. Teunissen, I. El Mellah, E. Chane, R. Keppens

Abstract: We report on the development of MPI-AMRVAC version 2.0, which is an open-source framework for parallel, grid-adaptive simulations of hydrodynamic and magnetohydrodynamic (MHD) astrophysical applications. The framework now supports radial grid stretching in combination with adaptive mesh refinement (AMR). The advantages of this combined approach are demonstrated with one-dimensional, two-dimensiona… ▽ More We report on the development of MPI-AMRVAC version 2.0, which is an open-source framework for parallel, grid-adaptive simulations of hydrodynamic and magnetohydrodynamic (MHD) astrophysical applications. The framework now supports radial grid stretching in combination with adaptive mesh refinement (AMR). The advantages of this combined approach are demonstrated with one-dimensional, two-dimensional and three-dimensional examples of spherically symmetric Bondi accretion, steady planar Bondi-Hoyle-Lyttleton flows, and wind accretion in Supergiant X-ray binaries. Another improvement is support for the generic splitting of any background magnetic field. We present several tests relevant for solar physics applications to demonstrate the advantages of field splitting on accuracy and robustness in extremely low plasma $β$ environments: a static magnetic flux rope, a magnetic null-point, and magnetic reconnection in a current sheet with either uniform or anomalous resistivity. Our implementation for treating anisotropic thermal conduction in multi-dimensional MHD applications is also described, which generalizes the original slope limited symmetric scheme from 2D to 3D. We perform ring diffusion tests that demonstrate its accuracy and robustness, and show that it prevents the unphysical thermal flux present in traditional schemes. The improved parallel scaling of the code is demonstrated with 3D AMR simulations of solar coronal rain, which show satisfactory strong scaling up to 2000 cores. Other framework improvements are also reported: the modernization and reorganization into a library, the handling of automatic regression tests, the use of inline/online Doxygen documentation, and a new future-proof data format for input/output △ Less

Submitted 17 October, 2017; originally announced October 2017.

Comments: 33 pages, 18 figures

Wind accretion onto compact objects

Authors: Ileyk El Mellah

Abstract: X-ray emission associated to accretion onto compact objects displays important levels of photometric and spectroscopic time-variability. When the accretor orbits a Supergiant star, it captures a fraction of the supersonic radiatively-driven wind which forms shocks in its vicinity. The amplitude and stability of this gravitational beaming of the flow conditions the mass accretion rate responsible,… ▽ More X-ray emission associated to accretion onto compact objects displays important levels of photometric and spectroscopic time-variability. When the accretor orbits a Supergiant star, it captures a fraction of the supersonic radiatively-driven wind which forms shocks in its vicinity. The amplitude and stability of this gravitational beaming of the flow conditions the mass accretion rate responsible, in fine, for the X-ray luminosity of those Supergiant X-ray Binaries. The capacity of this low angular momentum inflow to form a disc-like structure susceptible to be the stage of well-known instabilities remains at stake. Using state-of-the-art numerical setups, we characterized the structure of a Bondi-Hoyle-Lyttleton flow onto a compact object, from the shock down to the vicinity of the accretor, typically five orders of magnitude smaller. The evolution of the mass accretion rate and of the bow shock which forms around the accretor (transverse structure, opening angle, stability, temperature profile) with the Mach number of the incoming flow is described in detail. The robustness of those simulations based on the High Performance Computing MPI-AMRVAC code is supported by the topology of the inner sonic surface, in agreement with theoretical expectations. We developed a synthetic model of mass transfer in Supergiant X-ray Binaries which couples the launching of the wind accordingly to the stellar parameters, the orbital evolution of the streamlines in a modified Roche potential and the accretion process. We show that the shape of the permanent flow is entirely determined by the mass ratio, the filling factor, the Eddington factor and the alpha-force multiplier. Provided scales such as the orbital period are known, we can trace back the observables to evaluate the mass accretion rates, the accretion mechanism (stream or wind-dominated) and the shearing of the inflow. △ Less

Submitted 28 July, 2017; originally announced July 2017.

Comments: PhD manuscript

Report number: 2016USPCC120 MSC Class: 85

A numerical investigation of wind accretion in persistent Supergiant X-ray Binaries I - Structure of the flow at the orbital scale

Authors: I. El Mellah, F. Casse

Abstract: Classical Supergiant X-ray Binaries host a neutron star orbiting a supergiant OB star and display persistent X-ray luminosities of 10$^{35}$ to 10$^{37}$ erg/s. The stellar wind from the massive companion is believed to be the main source of matter accreted by the compact object. With this first paper, we introduce a ballistic model to characterize the structure of the wind at the orbital scale as… ▽ More Classical Supergiant X-ray Binaries host a neutron star orbiting a supergiant OB star and display persistent X-ray luminosities of 10$^{35}$ to 10$^{37}$ erg/s. The stellar wind from the massive companion is believed to be the main source of matter accreted by the compact object. With this first paper, we introduce a ballistic model to characterize the structure of the wind at the orbital scale as it accelerates, from the stellar surface to the vicinity of the accretor. Thanks to the parametrization we retained and the numerical pipeline we designed, we can investigate the supersonic flow and the subsequent observables as a function of a reduced set of characteristic numbers and scales. We show that the shape of the permanent flow is entirely determined by the mass ratio, the filling factor, the Eddington factor and the $α$-force multiplier which drives the stellar wind acceleration. Provided scales such as the orbital period are known, we can trace back the observables to evaluate the mass accretion rates, the accretion mechanism, the shearing of the inflow and the stellar parameters. We confront our model to three persistent Supergiant X-ray Binaries (Vela X-1, IGR J18027-2016, XTE J1855-026) and discuss the likelihood of wind-formed accretion discs around the accretors in each case, further investigated in a following paper. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Numerical simulations of axisymmetric hydrodynamical Bondi-Hoyle accretion onto a compact object

Authors: Ileyk El Mellah, Fabien Casse

Abstract: Bondi-Hoyle accretion configurations occur as soon as a gravitating body is immersed in an ambient medium with a supersonic relative velocity. From wind-accreting X-ray binaries to runaway neutron stars, such a regime has been witnessed many times and is believed to account for shock formation, the properties of which can be only marginally derived analytically. In this paper, we present the first… ▽ More Bondi-Hoyle accretion configurations occur as soon as a gravitating body is immersed in an ambient medium with a supersonic relative velocity. From wind-accreting X-ray binaries to runaway neutron stars, such a regime has been witnessed many times and is believed to account for shock formation, the properties of which can be only marginally derived analytically. In this paper, we present the first results of the numerical characterization of the stationary flow structure of Bondi-Hoyle accretion onto a compact object, from the large scale accretion radius down to the vicinity of the compact body. For different Mach numbers, we study the associated bow shock. It turns out that those simulations confirm the analytical prediction by Foglizzo & Ruffert (1996) concerning the topology of the inner sonic surface with an adiabatic index of 5/3. They also enable us to derive the related mass accretion rates, the position and the temperature of the bow shock, as function of the flow parameters, along with the transverse density and temperature profiles in the wake. △ Less

Submitted 25 September, 2015; originally announced September 2015.

Comments: 12 pages, 9 figures

Report number: MN-15-2576-MJ.R1

A Study of the Shortest-Period Planets Found With Kepler

Authors: Roberto Sanchis-Ojeda, Saul Rappaport, Joshua N. Winn, Michael C. Kotson, Alan M. Levine, Ileyk El Mellah

Abstract: We present the results of a survey aimed at discovering and studying transiting planets with orbital periods shorter than one day (ultra--short-period, or USP, planets), using data from the {\em Kepler} spacecraft. We computed Fourier transforms of the photometric time series for all 200,000 target stars, and detected transit signals based on the presence of regularly spaced sharp peaks in the Fou… ▽ More We present the results of a survey aimed at discovering and studying transiting planets with orbital periods shorter than one day (ultra--short-period, or USP, planets), using data from the {\em Kepler} spacecraft. We computed Fourier transforms of the photometric time series for all 200,000 target stars, and detected transit signals based on the presence of regularly spaced sharp peaks in the Fourier spectrum. We present a list of 106 USP candidates, of which 18 have not previously been described in the literature. In addition, among the objects we studied, there are 26 USP candidates that had been previously reported in the literature which do not pass our various tests. All 106 of our candidates have passed several standard tests to rule out false positives due to eclipsing stellar systems. A low false positive rate is also implied by the relatively high fraction of candidates for which more than one transiting planet signal was detected. By assuming these multi-transit candidates represent coplanar multi-planet systems, we are able to infer that the USP planets are typically accompanied by other planets with periods in the range 1-50 days, in contrast with hot Jupiters which very rarely have companions in that same period range. Another clear pattern is that almost all USP planets are smaller than 2 $R_\oplus$, possibly because gas giants in very tight orbits would lose their atmospheres by photoevaporation when subject to extremely strong stellar irradiation. Based on our survey statistics, USP planets exist around approximately $(0.51\pm 0.07)\%$ of G-dwarf stars, and $(0.83\pm 0.18)\%$ of K-dwarf stars. △ Less

Submitted 10 March, 2014; originally announced March 2014.

Comments: 20 pages, 10 figures. Submitted to ApJ. This version has been reviewed by a referee

Triple-Star Candidates Among the Kepler Binaries

Authors: S. Rappaport, K. Deck, A. Levine, T. Borkovits, J. Carter, I. El Mellah, R. Sanchis-Ojeda, B. Kalomeni

Abstract: We present the results of a search through the photometric database of eclipsing Kepler binaries (Prsa et al. 2011; Slawson et al. 2011) looking for evidence of hierarchical triple star systems. The presence of a third star orbiting the binary can be inferred from eclipse timing variations. We apply a simple algorithm in an automated determination of the eclipse times for all 2157 binaries. The "c… ▽ More We present the results of a search through the photometric database of eclipsing Kepler binaries (Prsa et al. 2011; Slawson et al. 2011) looking for evidence of hierarchical triple star systems. The presence of a third star orbiting the binary can be inferred from eclipse timing variations. We apply a simple algorithm in an automated determination of the eclipse times for all 2157 binaries. The "calculated" eclipse times, based on a constant period model, are subtracted from those observed. The resulting O-C (observed minus calculated times) curves are then visually inspected for periodicities in order to find triple-star candidates. After eliminating false positives due to the beat frequency between the ~1/2-hour Kepler cadence and the binary period, 39 candidate triple systems were identified. The periodic O-C curves for these candidates were then fit for contributions from both the classical Roemer delay and so-called "physical" delay, in an attempt to extract a number of the system parameters of the triple. We discuss the limitations of the information that can be inferred from these O-C curves without further supplemental input, e.g., ground-based spectroscopy. Based on the limited range of orbital periods for the triple star systems to which this search is sensitive, we can extrapolate to estimate that at least 20% of all close binaries have tertiary companions. △ Less

Submitted 15 September, 2013; v1 submitted 3 February, 2013; originally announced February 2013.

Comments: 19 pages, 13 figures, 3 tables; ApJ, 2013, 768, 33; corrected Fig. 7, updated references, minor fixes to text

Possible Disintegrating Short-Period Super-Mercury Orbiting KIC 12557548

Authors: S. Rappaport, A. Levine, E. Chiang, I. El Mellah, J. Jenkins, B. Kalomeni, E. S. Kite, M. Kotson, L. Nelson, L. Rousseau-Nepton, K. Tran

Abstract: We report here on the discovery of stellar occultations, observed with Kepler, that recur periodically at 15.685 hour intervals, but which vary in depth from a maximum of 1.3% to a minimum that can be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a K dwarf with T_eff = 4400 K and V = 16. Because the eclipse depths are highly variable, they cannot be due solely to tran… ▽ More We report here on the discovery of stellar occultations, observed with Kepler, that recur periodically at 15.685 hour intervals, but which vary in depth from a maximum of 1.3% to a minimum that can be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a K dwarf with T_eff = 4400 K and V = 16. Because the eclipse depths are highly variable, they cannot be due solely to transits of a single planet with a fixed size. We discuss but dismiss a scenario involving a binary giant planet whose mutual orbit plane precesses, bringing one of the planets into and out of a grazing transit. We also briefly consider an eclipsing binary, that either orbits KIC 12557548 in a hierarchical triple configuration or is nearby on the sky, but we find such a scenario inadequate to reproduce the observations. We come down in favor of an explanation that involves macroscopic particles escaping the atmosphere of a slowly disintegrating planet not much larger than Mercury. The particles could take the form of micron-sized pyroxene or aluminum oxide dust grains. The planetary surface is hot enough to sublimate and create a high-Z atmosphere; this atmosphere may be loaded with dust via cloud condensation or explosive volcanism. Atmospheric gas escapes the planet via a Parker-type thermal wind, dragging dust grains with it. We infer a mass loss rate from the observations of order 1 M_E/Gyr, with a dust-to-gas ratio possibly of order unity. For our fiducial 0.1 M_E planet, the evaporation timescale may be ~0.2 Gyr. Smaller mass planets are disfavored because they evaporate still more quickly, as are larger mass planets because they have surface gravities too strong to sustain outflows with the requisite mass-loss rates. The occultation profile evinces an ingress-egress asymmetry that could reflect a comet-like dust tail trailing the planet; we present simulations of such a tail. △ Less

Submitted 24 March, 2012; v1 submitted 12 January, 2012; originally announced January 2012.

Comments: 14 pages, 7 figures; submitted to ApJ, January 10, 2012; accepted March 21, 2012