Showing 1–50 of 76 results for author: Heida, M

Voronoi Graph -- Improved raycasting and integration schemes for high dimensional Voronoi diagrams

Authors: Alexander Sikorski, Martin Heida

Abstract: The computation of Voronoi Diagrams, or their dual Delauney triangulations is difficult in high dimensions. In a recent publication Polianskii and Pokorny propose an iterative randomized algorithm facilitating the approximation of Voronoi tesselations in high dimensions. In this paper, we provide an improved vertex search method that is not only exact but even faster than the bisection method that… ▽ More The computation of Voronoi Diagrams, or their dual Delauney triangulations is difficult in high dimensions. In a recent publication Polianskii and Pokorny propose an iterative randomized algorithm facilitating the approximation of Voronoi tesselations in high dimensions. In this paper, we provide an improved vertex search method that is not only exact but even faster than the bisection method that was previously recommended. Building on this we also provide a depth-first graph-traversal algorithm which allows us to compute the entire Voronoi diagram. This enables us to compare the outcomes with those of classical algorithms like qHull, which we either match or marginally beat in terms of computation time. We furthermore show how the raycasting algorithm naturally lends to a Monte Carlo approximation for the volume and boundary integrals of the Voronoi cells, both of which are of importance for finite Volume methods. We compare the Monte-Carlo methods to the exact polygonal integration, as well as a hybrid approximation scheme. △ Less

Submitted 16 May, 2024; originally announced May 2024.

The NuSTAR Serendipitous Survey: the 80-month catalog and source properties of the high-energy emitting AGN and quasar population

Authors: Claire L. Greenwell, Lizelke Klindt, George B. Lansbury, David J. Rosario, David M. Alexander, James Aird, Daniel Stern, Karl Forster, Michael J. Koss, Franz E. Bauer, Claudio Ricci, John Tomsick, William N. Brandt, Thomas Connor, Peter G. Boorman, Adlyka Annuar, David R. Ballantyne, Chien-Ting Chen, Francesca Civano, Andrea Comastri, Victoria A. Fawcett, Francesca M. Fornasini, Poshak Gandhi, Fiona Harrison, Marianne Heida , et al. (10 additional authors not shown)

Abstract: We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80-month (NSS80) catalog has an unprecedented $\sim$ 62 Ms of effective exposure time over 894 unique fields (a factor of three increase over the 40-month catalog), with an areal coverage of $\sim $36 deg… ▽ More We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80-month (NSS80) catalog has an unprecedented $\sim$ 62 Ms of effective exposure time over 894 unique fields (a factor of three increase over the 40-month catalog), with an areal coverage of $\sim $36 deg$^2$, larger than all NuSTAR extragalactic surveys. NSS80 provides 1274 hard X-ray sources in the $3-24$ keV band (822 new detections compared to the previous 40-month catalog). Approximately 76% of the NuSTAR sources have lower-energy ($<10$ keV) X-ray counterparts from Chandra, XMM-Newton, and Swift-XRT. We have undertaken an extensive campaign of ground-based spectroscopic follow-up to obtain new source redshifts and classifications for 427 sources. Combining these with existing archival spectroscopy provides redshifts for 550 NSS80 sources, of which 547 are classified. The sample is primarily composed of active galactic nuclei (AGN), detected over a large range in redshift ($z$ = 0.012-3.43), but also includes 58 spectroscopically confirmed Galactic sources. In addition, five AGN/galaxy pairs, one dual AGN system, one BL Lac candidate, and a hotspot of 4C 74.26 (radio quasar) have been identified. The median rest-frame $10-40$ keV luminosity and redshift of the NSS80 are $\langle{L_\mathrm{10-40 keV}}\rangle$ = 1.2 $\times$ 10$^{44}$ erg s$^{-1}$ and $\langle z \rangle = 0.56$. We investigate the optical properties and construct composite optical spectra to search for subtle signatures not present in the individual spectra, finding an excess of redder BL AGN compared to optical quasar surveys predominantly due to the presence of the host-galaxy and, at least in part, due to dust obscuration. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ:S. 57 pages, 32 figures

An ergodic and isotropic zero-conductance model with arbitrarily strong local connectivity

Authors: Martin Heida, Benedikt Jahnel, Anh Duc Vu

Abstract: We exhibit a percolating ergodic and isotropic lattice model in all but at least two dimensions that has zero effective conductivity in all spatial directions and for all non-trivial choices of the connectivity parameter. The model is based on the so-called randomly stretched lattice where we additionally elongate layers containing few open edges. We exhibit a percolating ergodic and isotropic lattice model in all but at least two dimensions that has zero effective conductivity in all spatial directions and for all non-trivial choices of the connectivity parameter. The model is based on the so-called randomly stretched lattice where we additionally elongate layers containing few open edges. △ Less

Submitted 11 March, 2024; originally announced March 2024.

MSC Class: 60K35 60K37; secondary: 60G55 90B18

Digging a little deeper: characterising three new extreme ULX candidates

Authors: T. P. Roberts, D. J. Walton, A. D. A. Mackenzie, M. Heida, S. Scaringi

Abstract: A prime motivation for compiling catalogues of any celestial X-ray source is to increase our numbers of rare sub-classes. In this work we take a recent multi-mission catalogue of ultraluminous X-ray sources (ULXs) and look for hitherto poorly-studied ULX candidates that are luminous ($L_{\rm X} \geq 10^{40} \rm ~erg~s^{-1}$), bright ($f_{\rm X} \geq 5 \times 10^{-13} \rm ~erg~cm~s^{-1}$) and have… ▽ More A prime motivation for compiling catalogues of any celestial X-ray source is to increase our numbers of rare sub-classes. In this work we take a recent multi-mission catalogue of ultraluminous X-ray sources (ULXs) and look for hitherto poorly-studied ULX candidates that are luminous ($L_{\rm X} \geq 10^{40} \rm ~erg~s^{-1}$), bright ($f_{\rm X} \geq 5 \times 10^{-13} \rm ~erg~cm~s^{-1}$) and have archival XMM-Newton data. We speculate that this luminosity regime may be ideal for identifying new pulsating ULXs (PULXs), given that the majority of known PULXs reach similar high luminosities. We find three sources that match our criteria, and study them using archival data. We find 4XMM J165251.5-591503 to possess a bright and variable Galactic optical/IR counterpart, and so conclude it is very likely to be a foreground interloper. 4XMM J091948.8-121429 does appear an excellent ULX candidate associated with the dwarf irregular galaxy PGC 26378, but has only one detection to date with low data quality. The best dataset belongs to 4XMM J112054.3+531040 which we find to be a moderately variable, spectrally hard ($Γ\approx 1.4$) X-ray source located in a spiral arm of NGC 3631. Its spectral hardness is similar to known PULXs, but no pulsations are detected by accelerated pulsation searches in the available data. We discuss whether other missions provide objects for similar studies, and compare this method to others suggested for identifying good PULX candidates. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 15 pages, 9 figures, accepted for publication in MNRAS

A new sample of transient ultraluminous X-ray sources serendipitously discovered by Swift/XRT

Authors: Murray Brightman, Jean-Marie Hameury, Jean-Pierre Lasota, Ranieri D. Baldi, Gabriele Bruni, Jenna M. Cann, Hannah Earnshaw, Felix Fürst, Marianne Heida, Amruta Jaodand, Margaret Lazzarini, Matthew J. Middleton, Dominic J. Walton, Kimberly A. Weaver

Abstract: Ultraluminous X-ray sources (ULXs) are our best laboratories for studying extreme super-Eddington accretion. Most studies of these objects are of relatively persistent sources, however there is growing evidence to suggest a large fraction of these sources are transient. Here we present a sample of five newly reported transient ULXs in the galaxies NGC 4945, NGC 7793 and M81 serendipitously discove… ▽ More Ultraluminous X-ray sources (ULXs) are our best laboratories for studying extreme super-Eddington accretion. Most studies of these objects are of relatively persistent sources, however there is growing evidence to suggest a large fraction of these sources are transient. Here we present a sample of five newly reported transient ULXs in the galaxies NGC 4945, NGC 7793 and M81 serendipitously discovered in Swift/XRT observations. Swift monitoring of these sources have provided well sampled lightcurves, allowing for us to model the lightcurves with the disk instability model of Hameury & Lasota (2020) which implies durations of 60-400 days and that the mass accretion rate through the disk is close to or greater than the Eddington rate. Of the three source regions with prior HST imaging, color magnitude diagrams of the potential stellar counterparts show varying ages of the possible stellar counterparts. Our estimation of the rates of these sources in these three galaxies is 0.4-1.3 year$^{-1}$. We find that while persistent ULXs dominate the high end of galaxy luminosity functions, the number of systems that produce ULX luminosities are likely dominated by transient sources. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ

The Enhanced Resolution Imager and Spectrograph for the VLT

Authors: R. Davies, O. Absil, G. Agapito, A. Agudo Berbel, A. Baruffolo, V. Biliotti, M. Bonaglia, M. Bonse, R. Briguglio, P. Campana, Y. Cao, L. Carbonaro, A. Cortes, G. Cresci, Y. Dallilar, F. Dannert, R. J. De Rosa, M. Deysenroth, I. Di Antonio, A. Di Cianno, G. Di Rico, D. Doelman, M. Dolci, R. Dorn, F. Eisenhauer , et al. (59 additional authors not shown)

Abstract: ERIS, the Enhanced Resolution Imager and Spectrograph, is an instrument that both extends and enhances the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It replaces two instruments that were being maintained beyond their operational lifetimes, combines their functionality on a single focus, provides a new wavefront sensing module for natural and laser guide stars… ▽ More ERIS, the Enhanced Resolution Imager and Spectrograph, is an instrument that both extends and enhances the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It replaces two instruments that were being maintained beyond their operational lifetimes, combines their functionality on a single focus, provides a new wavefront sensing module for natural and laser guide stars that makes use of the Adaptive Optics Facility, and considerably improves on their performance. The observational modes ERIS provides are integral field spectroscopy at 1-2.5 μm, imaging at 1-5 μm with several options for high contrast imaging, and longslit spectroscopy at 3-4 μm, The instrument is installed at the Cassegrain focus of UT4 at the VLT and, following its commissioning during 2022, has been made available to the community. △ Less

Submitted 26 April, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 20 pages with 29 figures; accepted for A&A (minor changes)

Probing the nature of the low state in the extreme ultraluminous X-ray pulsar NGC 5907 ULX1

Authors: F. Fuerst, D. J. Walton, G. L. Israel, M. Bachetti, D. Barret, M. Brightman, H. P. Earnshaw, A. Fabian, M. Heida, M. Imbrogno, M. J. Middleton, C. Pinto, R. Salvaterra, T. P. Roberts, G. A. Rodríguez Castillo, N. Webb

Abstract: NGC 5907 ULX1 is the most luminous ultra-luminous X-ray pulsar (ULXP) known to date, reaching luminosities in excess of 1e41 erg/s. The pulsar is known for its fast spin-up during the on-state. Here, we present a long-term monitoring of the X-ray flux and the pulse period between 2003-2022. We find that the source was in an off- or low-state between mid-2017 to mid-2020. During this state, our pul… ▽ More NGC 5907 ULX1 is the most luminous ultra-luminous X-ray pulsar (ULXP) known to date, reaching luminosities in excess of 1e41 erg/s. The pulsar is known for its fast spin-up during the on-state. Here, we present a long-term monitoring of the X-ray flux and the pulse period between 2003-2022. We find that the source was in an off- or low-state between mid-2017 to mid-2020. During this state, our pulse period monitoring shows that the source had spun down considerably. We interpret this spin-down as likely being due to the propeller effect, whereby accretion onto the neutron star surface is inhibited. Using state-of-the-art accretion and torque models, we use the spin-up and spin-down episodes to constrain the magnetic field. For the spin-up episode, we find solutions for magnetic field strengths of either around 1e12G or 1e13G, however, the strong spin-down during the off-state seems only to be consistent with a very high magnetic field, namely, >1e13G. This is the first time a strong spin-down is seen during a low flux state in a ULXP. Based on the assumption that the source entered the propeller regime, this gives us the best estimate so far for the magnetic field of NGC 5907 ULX1. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 10 pages, 1 figure, accepted for publication in A&A

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

First on-sky results of ERIS at VLT

Authors: Kateryna Kravchenko, Yigit Dallilar, Olivier Absil, Alex Agudo Berbel, Andrea Baruffolo, Markus J. Bonse, Alexander Buron, Yixian Cao, Angela Cortes, Felix Dannert, Richard Davies, Robert J. De Rosa, Matthias Deysenroth, David S. Doelman, Frank Eisenhauer, Simone Esposito, Helmut Feuchtgruber, Natascha Förster Schreiber, Xiaofeng Gao, Hans Gemperlein, Reinhard Genzel, Stefan Gillessen, Christian Ginski, Adrian M. Glauser, Andreas Glindemann , et al. (24 additional authors not shown)

Abstract: ERIS (Enhanced Resolution Imager and Spectrograph) is a new adaptive optics instrument installed at the Cassegrain focus of the VLT-UT4 telescope at the Paranal Observatory in Chile. ERIS consists of two near-infrared instruments: SPIFFIER, an integral field unit (IFU) spectrograph covering J to K bands, and NIX, an imager covering J to M bands. ERIS has an adaptive optics system able to work with… ▽ More ERIS (Enhanced Resolution Imager and Spectrograph) is a new adaptive optics instrument installed at the Cassegrain focus of the VLT-UT4 telescope at the Paranal Observatory in Chile. ERIS consists of two near-infrared instruments: SPIFFIER, an integral field unit (IFU) spectrograph covering J to K bands, and NIX, an imager covering J to M bands. ERIS has an adaptive optics system able to work with both LGS and NGS. The Assembly Integration Verification (AIV) phase of ERIS at the Paranal Observatory was carried out starting in December 2021, followed by several commissioning runs in 2022. This contribution will describe the first preliminary results of the on-sky performance of ERIS during its commissioning and the future perspectives based on the preliminary scientific results. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: Proceeding of SPIE Astronomical Telescopes + Instrumentation 2022

BASS XXII: The BASS DR2 AGN Catalog and Data

Authors: Michael J. Koss, Claudio Ricci, Benny Trakhtenbrot, Kyuseok Oh, Jakob S. den Brok, Julian E. Mejia-Restrepo, Daniel Stern, George C. Privon, Ezequiel Treister, Meredith C. Powell, Richard Mushotzky, Franz E. Bauer, Tonima T. Ananna, Mislav Balokovic, Rudolf E. Bar, George Becker, Patricia Bessiere, Leonard Burtscher, Turgay Caglar, Enrico Congiu, Phil Evans, Fiona Harrison, Marianne Heida, Kohei Ichikawa, Nikita Kamraj , et al. (10 additional authors not shown)

Abstract: We present the AGN catalog and optical spectroscopy for the second data release of the Swift BAT AGN Spectroscopic Survey (BASS DR2). With this DR2 release we provide 1425 optical spectra, of which 1181 are released for the first time, for the 858 hard X-ray selected AGN in the Swift BAT 70-month sample. The majority of the spectra (813/1425, 57%) are newly obtained from VLT/Xshooter or Palomar/Do… ▽ More We present the AGN catalog and optical spectroscopy for the second data release of the Swift BAT AGN Spectroscopic Survey (BASS DR2). With this DR2 release we provide 1425 optical spectra, of which 1181 are released for the first time, for the 858 hard X-ray selected AGN in the Swift BAT 70-month sample. The majority of the spectra (813/1425, 57%) are newly obtained from VLT/Xshooter or Palomar/Doublespec. Many of the spectra have both higher resolution (R>2500, N~450) and/or very wide wavelength coverage (3200-10000 A, N~600) that are important for a variety of AGN and host galaxy studies. We include newly revised AGN counterparts for the full sample and review important issues for population studies, with 44 AGN redshifts determined for the first time and 780 black hole mass and accretion rate estimates. This release is spectroscopically complete for all AGN (100%, 858/858) with 99.8% having redshift measurements (857/858) and 96% completion in black hole mass estimates of unbeamed AGN (outside the Galactic plane). This AGN sample represents a unique census of the brightest hard X-ray selected AGN in the sky, spanning many orders of magnitude in Eddington ratio (Ledd=10^-5-100), black hole mass (MBH=10^5-10^10 Msun), and AGN bolometric luminosity (Lbol=10^40-10^47 ergs/s). △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 34 pages, 11 figures, 15 tables, published in ApJS as part of BASS DR2 special issue

Journal ref: ApJS, 261, 2 (2022)

Detecting Stripped Stars While Searching for Quiescent Black Holes

Authors: J. Bodensteiner, M. Heida, M. Abdul-Masih, D. Baade, G. Banyard, D. M. Bowman, M. Fabry, A. Frost, L. Mahy, P. Marchant, A. Mérand, M. Reggiani, Th. Rivinius, H. Sana, F. Selman, T. Shenar

Abstract: While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent bl… ▽ More While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent black hole candidates have been proposed, the majority have been challenged by follow-up investigations. A confusion that emerged recently concerns binary systems that appear to contain a normal B-type star with an unseen companion, believed to be a black hole. On closer inspection, some of these seemingly normal B-type stars instead turn out to be stars stripped of most of their mass through an interaction with their binary companion, which in at least two cases is a rapidly rotating star rather than a compact object. These contaminants in the search for quiescent black holes are themselves extremely interesting objects as they represent a rare phase of binary evolution, and should be given special attention when searching for binaries hosting black holes in large spectroscopic studies. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 7 pages, 5 figures

Journal ref: The Messenger, 2022, vol. 186, p. 3-9

The variability behavior of NGC 925 ULX-3

Authors: Hannah P. Earnshaw, Murray Brightman, Fiona A. Harrison, Marianne Heida, Amruta Jaodand, Matthew J. Middleton, Timothy P. Roberts, Dominic J. Walton

Abstract: We report the results of a 2019-2021 monitoring campaign with Swift and associated target-of-opportunity observations with XMM-Newton and NuSTAR, examining the spectral and timing behavior of the highly variable ultraluminous X-ray source (ULX) NGC 925 ULX-3. We find that the source exhibits a 127-128 day periodicity, with fluxes typically ranging from 1e-13 to 8e-13 ergs/s/cm2. We do not find str… ▽ More We report the results of a 2019-2021 monitoring campaign with Swift and associated target-of-opportunity observations with XMM-Newton and NuSTAR, examining the spectral and timing behavior of the highly variable ultraluminous X-ray source (ULX) NGC 925 ULX-3. We find that the source exhibits a 127-128 day periodicity, with fluxes typically ranging from 1e-13 to 8e-13 ergs/s/cm2. We do not find strong evidence for a change in period over the time that NGC 925 ULX-3 has been observed, although the source may have been in a much lower flux state when first observed with Chandra in 2005. We do not detect pulsations, and we place an upper limit on the pulsed fraction of ~40% in the XMM-Newton band, consistent with some previous pulsation detections at low energies in other ULXs. The source exhibits a typical ULX spectrum that turns over in the NuSTAR band and can be fitted using two thermal components. These components have a high temperature ratio that may indicate the lack of extreme inner disk truncation by a magnetar-level magnetic field. We examine the implications for a number of different models for superorbital periods in ULXs, finding that a neutron star with a magnetic field of ~10^12 G may be plausible for this source. The future detection of pulsations from this source would allow for the further testing and constraining of such models. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in ApJ

The Redshift Evolution of Ultraluminous X-Ray Sources out to z~0.5: Comparison with X-Ray Binary Populations and Contribution to the Cosmic X-Ray Background

Authors: R. Scott Barrows, Julia M. Comerford, Daniel Stern, Marianne Heida

Abstract: Ultraluminous X-ray sources (ULXs) are thought to be powerful X-ray binaries (XRBs) and may contribute significantly to the redshift-dependent X-ray emission from star forming galaxies. We have assembled a uniform sample of 259 ULXs over the redshift range z=0.002-0.51 to constrain their physical nature and their contribution to the Cosmic X-Ray Background (CXB). The sample is constructed by cross… ▽ More Ultraluminous X-ray sources (ULXs) are thought to be powerful X-ray binaries (XRBs) and may contribute significantly to the redshift-dependent X-ray emission from star forming galaxies. We have assembled a uniform sample of 259 ULXs over the redshift range z=0.002-0.51 to constrain their physical nature and their contribution to the Cosmic X-Ray Background (CXB). The sample is constructed by crossmatching galaxies from the Sloan Digital Sky Survey with the Chandra Source Catalog and selecting off-nuclear X-ray sources after applying astrometric corrections. The fraction of contaminants is ~30% and shows no evolution with redshift. The host galaxy star formation rates (SFRs) are systematically elevated relative to the parent sample when matched in host stellar mass. The specific SFRs suggest a slight preference for high-mass XRBs, and the X-ray luminosity scaling relations with host galaxy stellar mass and SFR indicate that the highest redshift sources represent relatively luminous XRB populations that dominate their host galaxy X-ray emission. The fraction of galaxies hosting at least one ULX of a given luminosity increases with redshift over the full range of our sample, as expected if ULXs are preferentially found in galaxies with high SFRs and low metallicities. At z~0.5, the ULX X-ray flux is consistent with the X-ray emission from star-forming galaxies. Moreover, ULXs may account for up to ~40% of the integrated flux from XRBs in the normal galaxy population out to z~0.5, suggesting they may contribute significantly to the overall ionizing radiation from galaxies. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: 21 pages, 17 figures. Accepted for publication in the Astrophysical Journal

An 8.56 keV absorption line in the hyperluminous X-ray source in NGC 4045: ultra-fast outflow or cyclotron line?

Authors: Murray Brightman, Peter Kosec, Felix Fürst, Hannah Earnshaw, Marianne Heida, Matthew J Middleton, Daniel Stern, Dominic J Walton

Abstract: We report on the discovery of an absorption line at $E=8.56^{+0.05}_{-0.11}$ keV detected with a significance of $>3.3σ$ in the NuSTAR and XMM-Newton spectra of a newly discovered hyperluminous X-ray source (HLX, $L_{\rm X}>10^{41}$ ergs$^{-1}$) in the galaxy NGC 4045 at a distance of 32 Mpc. The source was first discovered serendipitously in a Swift/XRT observation of the galaxy, and Swift monito… ▽ More We report on the discovery of an absorption line at $E=8.56^{+0.05}_{-0.11}$ keV detected with a significance of $>3.3σ$ in the NuSTAR and XMM-Newton spectra of a newly discovered hyperluminous X-ray source (HLX, $L_{\rm X}>10^{41}$ ergs$^{-1}$) in the galaxy NGC 4045 at a distance of 32 Mpc. The source was first discovered serendipitously in a Swift/XRT observation of the galaxy, and Swift monitoring reveals a highly variable source changing by over an order of magnitude from maximum to minimum. The origin of the absorption line appears likely to be by highly ionized iron with a blue shift of 0.19$c$, indicating an ultrafast outflow (UFO). However, the large equivalent width of the line (EW$=-0.22^{+0.08}_{-0.09}$ keV) paired with the lack of other absorption lines detected are difficult to reconcile with models. An alternative explanation is that the line is due to a cyclotron resonance scattering feature (CRSF), produced by the interaction of X-ray photons with the powerful magnetic field of a neutron star. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

HR 6819 is a binary system with no black hole -- revisiting the source with infrared interferometry and optical integral field spectroscopy

Authors: A. J. Frost, J. Bodensteiner, Th. Rivinius, D. Baade, A. Merand, F. Selman, M. Abdul-Masih, G. Banyard, E. Bordier, K. Dsilva, C. Hawcroft, L. Mahy, M. Reggiani, T. Shenar, M. Cabezas, P. Hadrava, M. Heida, R. Klement, H. Sana

Abstract: Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stella… ▽ More Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 16 pages, 11 figures, subject of ESO press release eso2204

Journal ref: A&A, Volume 659, March 2022, L3, 12, Letters to the Editor

Orbital decay in M82 X-2

Authors: Matteo Bachetti, Marianne Heida, Thomas Maccarone, Daniela Huppenkothen, Gian Luca Israel, Didier Barret, Murray Brightman, McKinley Brumback, Hannah P. Earnshaw, Karl Forster, Felix Fürst, Brian W. Grefenstette, Fiona A. Harrison, Amruta D. Jaodand, Kristin K. Madsen, Matthew Middleton, Sean N. Pike, Maura Pilia, Juri Poutanen, Daniel Stern, John A. Tomsick, Dominic J. Walton, Natalie Webb, Jörn Wilms

Abstract: M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenar… ▽ More M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenarios. In this Paper, we follow the orbit of the neutron star for seven years, measure the decay of the orbit ($\dot{P}_{orb}/{P}_{orb}\approx-8\cdot10^{-6}\mathrm{yr}^{-1}$), and argue that this orbital decay is driven by extreme mass transfer of more than 150 times the mass transfer limit set by the Eddington luminosity. If this is true, the mass available to the accretor is more than enough to justify its luminosity, with no need for beaming. This also strongly favors models where the accretor is a highly-magnetized neutron star. △ Less

Submitted 29 August, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 18 pages, 6 figures, accepted for publication in ApJ

Journal ref: ApJ 937, 125 (2022)

Science with the Ultraviolet Explorer (UVEX)

Authors: S. R. Kulkarni, Fiona A. Harrison, Brian W. Grefenstette, Hannah P. Earnshaw, Igor Andreoni, Danielle A. Berg, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Jessie L. Christiansen, Michael W. Coughlin, Alexander Wuollet Criswell, Behnam Darvish, Kaustav K. Das, Kishalay De, Luc Dessart, Don Dixon, Bas Dorsman, Kareem El-Badry, Christopher Evans, K. E. Saavik Ford, Christoffer Fremling, Boris T. Gansicke, Suvi Gezari, Y. Goetberg , et al. (31 additional authors not shown)

Abstract: UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky… ▽ More UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program. △ Less

Submitted 17 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 69 pages, 43 figures

Evolution of the spin, spectrum and super-orbital period of the ultraluminous X-ray pulsar M51 ULX7

Authors: Murray Brightman, Matteo Bachetti, Hannah Earnshaw, Felix Fürst, Marianne Heida, Gian Luca Israel, Sean Pike, Daniel Stern, Dominic J Walton

Abstract: M51 ULX7 is among a small group of known ultraluminous X-ray pulsars (ULXP). The neutron star powering the source has a spin period of 2.8s, orbits its companion star with a period of 2 days, and a super-orbital period of 38 days is evident in its X-ray lightcurve. Here we present NuSTAR and XMM-Newton data on the source from 2019 obtained when the source was near its peak brightness. We detect th… ▽ More M51 ULX7 is among a small group of known ultraluminous X-ray pulsars (ULXP). The neutron star powering the source has a spin period of 2.8s, orbits its companion star with a period of 2 days, and a super-orbital period of 38 days is evident in its X-ray lightcurve. Here we present NuSTAR and XMM-Newton data on the source from 2019 obtained when the source was near its peak brightness. We detect the pulsations, having spun up at a rate of 3$\pm0.5\times10^{-10}$ s s$^{-1}$ since they were previously detected in 2018. The data also provide the first high-quality broadband spectrum of the source. We find it to be very similar to that of other ULXPs, with two disk-like components, and a high energy tail. When combined with XMM-Newton data obtained in 2018, we explore the evolution of the spectral components with super-orbital phase, finding that the luminosity of the hotter component drives the super-orbital flux modulation. The inclination the disk components appear to change with phase, which may support the idea that these super-orbital periods are caused by disk precession. We also reexamine the super-orbital period with 3 years of Swift/XRT monitoring, finding that the period is variable, increasing from 38.2$\pm0.5$ days in 2018--2019 to 44.2$\pm0.9$ days in 2020--2021, which rules out alternative explanations for the super-orbital period. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJ

Stochastic Homogenization on Irregularly Perforated Domains

Authors: Martin Heida, Benedikt Jahnel, Anh Duc Vu

Abstract: We study stochastic homogenization of a quasilinear parabolic PDE with nonlinear microscopic Robin conditions on a perforated domain. The focus of our work lies on the underlying geometry that does not allow standard homogenization techniques to be applied directly. Instead we prove homogenization on a regularized geometry and demonstrate afterwards that the form of the homogenized equation is ind… ▽ More We study stochastic homogenization of a quasilinear parabolic PDE with nonlinear microscopic Robin conditions on a perforated domain. The focus of our work lies on the underlying geometry that does not allow standard homogenization techniques to be applied directly. Instead we prove homogenization on a regularized geometry and demonstrate afterwards that the form of the homogenized equation is independent from the regularization. Then we pass to the regularization limit to obtain the anticipated limit equation. Furthermore, we show that Boolean models of Poisson point processes are covered by our approach. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 41 pages, 7 figures

MSC Class: 74Q05; 47J30; 60K35

$ν$ Gem: a hierarchical triple system with an outer Be star

Authors: Robert Klement, Petr Hadrava, Thomas Rivinius, Dietrich Baade, Mauricio Cabezas, Marianne Heida, Gail H. Schaefer, Tyler Gardner, Douglas R. Gies, Narsireddy Anugu, Cyprien Lanthermann, Claire L. Davies, Matthew D. Anderson, John D. Monnier, Jacob Ennis, Aaron Labdon, Benjamin R. Setterholm, Stefan Kraus, Theo A. ten Brummelaar, Jean-Baptiste le Bouquin

Abstract: Time series of spectroscopic, speckle-interferometric, and optical long-baseline-interferometric observations confirm that $ν$ Gem is a hierarchical triple system. It consists of an inner binary composed of two B-type stars and an outer classical Be star. Several photospheric spectral lines of the inner components were disentangled, revealing two stars with very different rotational broadening (… ▽ More Time series of spectroscopic, speckle-interferometric, and optical long-baseline-interferometric observations confirm that $ν$ Gem is a hierarchical triple system. It consists of an inner binary composed of two B-type stars and an outer classical Be star. Several photospheric spectral lines of the inner components were disentangled, revealing two stars with very different rotational broadening ($\sim$260 and $\sim$140 kms$^{-1}$, respectively), while the photospheric lines of the Be star remain undetected. From the combined spectroscopic and astrometric orbital solution it is not possible to unambiguously cross-identify the inner astrometric components with the spectroscopic components. In the preferred solution based on modeling of the disentangled line profiles, the inner binary is composed of two stars with nearly identical masses of 3.3 M$_\odot$ and the more rapidly rotating star is the fainter one. These two stars are in a marginally elliptical orbit ($e$ = 0.06) about each other with a period of 53.8 d. The third star also has a mass of 3.3 M$_\odot$ and follows a more eccentric ($e$ = 0.24) orbit with a period of 19.1 yr. The two orbits are co-directional and, at inclinations of 79$^{\circ}$ and 76$^{\circ}$ of the inner and the outer orbit, respectively, about coplanar. No astrometric or spectroscopic evidence could be found that the Be star itself is double. The system appears dynamically stable and not subject to eccentric Lidov-Kozai oscillations. After disentangling, the spectra of the components of the inner binary do not exhibit peculiarities that would be indicative of past interactions. Motivations for a wide range of follow-up studies are suggested. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: accepted to ApJ

Stochastic two-scale convergence and Young measures

Authors: Martin Heida, Stefan Neukamm, Mario Varga

Abstract: In this paper we compare the notion of stochastic two-scale convergence in the mean (by Bourgeat, Mikelić and Wright), the notion of stochastic unfolding (recently introduced by the authors), and the quenched notion of stochastic two-scale convergence (by Zhikov and Pyatnitskii). In particular, we introduce stochastic two-scale Young measures as a tool to compare mean and quenched limits. Moreover… ▽ More In this paper we compare the notion of stochastic two-scale convergence in the mean (by Bourgeat, Mikelić and Wright), the notion of stochastic unfolding (recently introduced by the authors), and the quenched notion of stochastic two-scale convergence (by Zhikov and Pyatnitskii). In particular, we introduce stochastic two-scale Young measures as a tool to compare mean and quenched limits. Moreover, we discuss two examples, which can be naturally analyzed via stochastic unfolding, but which cannot be treated via quenched stochastic two-scale convergence. △ Less

Submitted 26 May, 2021; originally announced May 2021.

Comments: 30 pages. arXiv admin note: substantial text overlap with arXiv:1805.09546

MSC Class: 74Q05; 47J30

Stochastic homogenization on perforated domains I -- Extension operators

Authors: Martin Heida

Abstract: We study the existence of uniformly bounded extension and trace operators for W^{1,p}-functions on randomly perforated domains, where the geometry is assumed to be stationary ergodic. Such extension and trace operators are important for compactness in stochastic homogenization. In contrast to former approaches and results, we use very weak assumptions on the geometry which we call local (δ,M)-regu… ▽ More We study the existence of uniformly bounded extension and trace operators for W^{1,p}-functions on randomly perforated domains, where the geometry is assumed to be stationary ergodic. Such extension and trace operators are important for compactness in stochastic homogenization. In contrast to former approaches and results, we use very weak assumptions on the geometry which we call local (δ,M)-regularity, isotropic cone mixing and bounded average connectivity. The first concept measures local Lipschitz regularity of the domain while the second measures the mesoscopic distribution of void space. The third is the most tricky part and measures the "mesoscopic" connectivity of the geometry. In contrast to former approaches we do not require a minimal distance between the inclusions and we allow for globally unbounded Lipschitz constants and percolating holes. We will illustrate our method by applying it to the Boolean model based on a Poisson point process and to a Delaunay pipe process. This is Part I of a series of 3 papers that will be uploaded during the next weeks. It is part of a major revision of the original ArXiv 2001.10373. The original work was corrected, improved and extended significantly. In particular, the theory is now able to also deal with extensions that preserve the norm of the symmetrized gradient. △ Less

Submitted 9 December, 2021; v1 submitted 23 May, 2021; originally announced May 2021.

Comments: arXiv admin note: text overlap with arXiv:2001.10373

Long-term pulse period evolution of the ultra-luminous X-ray pulsar NGC 7793 P13

Authors: F. Fuerst, D. J. Walton, M. Heida, M. Bachetti, C. Pinto, M. J. Middleton, M. Brightman, H. P. Earnshaw, D. Barret, A. C. Fabian, P. Kretschmar, K. Pottschmidt, A. Ptak, T. Roberts, D. Stern, N. Webb, J. Wilms

Abstract: Ultra-luminous X-ray pulsars (ULXPs) provide a unique opportunity to study super-Eddington accretion. We present the results of a monitoring campaign of ULXP NGC 7793 P13. Over our four-year monitoring campaign with Swift, XMM-Newton, and NuSTAR, we measured a continuous spin-up with $\dot P$ ~ -3.8e-11 s/s. The strength of the spin-up is independent of the observed X-ray flux, indicating that des… ▽ More Ultra-luminous X-ray pulsars (ULXPs) provide a unique opportunity to study super-Eddington accretion. We present the results of a monitoring campaign of ULXP NGC 7793 P13. Over our four-year monitoring campaign with Swift, XMM-Newton, and NuSTAR, we measured a continuous spin-up with $\dot P$ ~ -3.8e-11 s/s. The strength of the spin-up is independent of the observed X-ray flux, indicating that despite a drop in observed flux in 2019, accretion onto the source has continued at largely similar rates. The source entered an apparent off-state in early 2020, which might have resulted in a change in the accretion geometry as no pulsations were found in observations in July and August 2020. We used the long-term monitoring to update the orbital ephemeris and the periodicities seen in both the observed optical/UV and X-ray fluxes. We find that the optical/UV period is very stable over the years, with $P_\text{UV}$ = 63.75 (+0.17, -0.12) d. The best-fit orbital period determined from our X-ray timing results is 64.86 +/- 0.19 d, which is almost a day longer than previously implied, and the X-ray flux period is 65.21+/- 0.15 d, which is slightly shorter than previously measured. The physical origin of these different flux periods is currently unknown. We study the hardness ratio to search for indications of spectral changes. We find that the hardness ratios at high energies are very stable and not directly correlated with the observed flux. At lower energies we observe a small hardening with increased flux, which might indicate increased obscuration through outflows at higher luminosities. We find that the pulsed fraction is significantly higher at low fluxes. This seems to imply that the accretion geometry already changed before the source entered the deep off-state. We discuss possible scenarios to explain this behavior, which is likely driven by a precessing accretion disk. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 10 pages, 6 figures, accepted for publication in A&A

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

A New Transient Ultraluminous X-ray Source in NGC 7090

Authors: D. J. Walton, M. Heida, M. Bachetti, F. Furst, M. Brightman, H. Earnshaw, P. A. Evans, A. C. Fabian, B. W. Grefenstette, F. A. Harrison, G. L. Israel, G. B. Lansbury, M. J. Middleton, S. Pike, V. Rana, T. P. Roberts, G. A. Rodriguez Castillo, R. Salvaterra, X. Song, D. Stern

Abstract: We report on the discovery of a new, transient ultraluminous X-ray source (ULX) in the galaxy NGC 7090. This new ULX, which we refer to as NGC 7090 ULX3, was discovered via monitoring with $Swift$ during 2019-20, and to date has exhibited a peak luminosity of $L_{\rm{X}} \sim 6 \times 10^{39}$ erg s$^{-1}$. Archival searches show that, prior to its recent transition into the ULX regime, ULX3 appea… ▽ More We report on the discovery of a new, transient ultraluminous X-ray source (ULX) in the galaxy NGC 7090. This new ULX, which we refer to as NGC 7090 ULX3, was discovered via monitoring with $Swift$ during 2019-20, and to date has exhibited a peak luminosity of $L_{\rm{X}} \sim 6 \times 10^{39}$ erg s$^{-1}$. Archival searches show that, prior to its recent transition into the ULX regime, ULX3 appeared to exhibit a fairly stable luminosity of $L_{\rm{X}} \sim 10^{38}$ erg s$^{-1}$. Such strong long-timescale variability may be reminiscent of the small population of known ULX pulsars, although deep follow-up observations with $XMM$-$Newton$ and $NuSTAR$ do not reveal any robust X-ray pulsation signals. Pulsations similar to those seen from known ULX pulsars cannot be completely excluded, however, as the limit on the pulsed fraction of any signal that remains undetected in these data is $\lesssim$20\%. The broadband spectrum from these observations is well modelled with a simple thin disc model, consistent with sub-Eddington accretion, which may instead imply a moderately large black hole accretor ($M_{\rm{BH}} \sim 40 ~ M_{\odot}$). Similarly, though, more complex models consistent with the super-Eddington spectra seen in other ULXs (and the known ULX pulsars) cannot be excluded given the limited signal-to-noise of the available broadband data. The nature of the accretor powering this new ULX therefore remains uncertain. △ Less

Submitted 19 November, 2020; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 11 pages, 7 figures; accepted for publication in MNRAS

NIR Counterparts to ULXs (III): Completing the photometric survey and selected spectroscopic results

Authors: K. M. López, M. Heida, P. G. Jonker, M. A. P. Torres, T. P. Roberts, D. J. Walton, D. -S. Moon, F. A. Harrison

Abstract: We present results from the remaining sources in our search for near-infrared (NIR) candidate counterparts to ultraluminous X-ray sources (ULXs) within $\simeq$ 10 Mpc. We observed 23 ULXs in 15 galaxies and detected NIR candidate counterparts to six of them. Two of these have an absolute magnitude consistent with a single red supergiant (RSG). Three counterparts are too bright for a RSG and spati… ▽ More We present results from the remaining sources in our search for near-infrared (NIR) candidate counterparts to ultraluminous X-ray sources (ULXs) within $\simeq$ 10 Mpc. We observed 23 ULXs in 15 galaxies and detected NIR candidate counterparts to six of them. Two of these have an absolute magnitude consistent with a single red supergiant (RSG). Three counterparts are too bright for a RSG and spatially extended, and thus we classify them as stellar clusters. The other candidate is too faint for a RSG. Additionally, we present the results of our NIR spectroscopic follow-up of five sources: four originally classified as RSG and one as a stellar cluster on the basis of previous photometry. The stellar cluster candidate is actually a nebula. Of the four RSGs candidates, one source has a broad H$α$ emission line redshifted by $\sim z = 1$, making it a background AGN. Two other sources show stellar spectra consistent with them being RSGs. The final RSG candidate is too faint to classify, but does not show strong (nebular) emission lines in its spectrum. After our search for NIR counterparts to 113 ULXs, where we detected a candidate counterpart for 38 ULXs, we have spectroscopically confirmed the nature of 12: five sources are nebulae, one source is not classified, one source is an AGN and five are RSGs. These possible five ULX-RSG binary systems would constitute $\simeq$ $(4 \pm 2)\%$ of the observed ULXs, a fraction almost four times larger than what was predicted by binary evolution simulations. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 20 pages, 10 figures, 8 tables, accepted for publication in MNRAS

A naked-eye triple system with a nonaccreting black hole in the inner binary

Authors: Th. Rivinius, D. Baade, P. Hadrava, M. Heida, R. Klement

Abstract: Several dozen optical echelle spectra demonstrate that HR 6819 is a hierarchical triple. A classical Be star is in a wide orbit with an unconstrained period around an inner 40 d binary consisting of a B3 III star and an unseen companion in a circular orbit. The radial-velocity semi-amplitude of 61.3 km/s of the inner star and its minimum (probable) mass of 5.0 Msun (6.3 +- 0.7 Msun) imply a mass o… ▽ More Several dozen optical echelle spectra demonstrate that HR 6819 is a hierarchical triple. A classical Be star is in a wide orbit with an unconstrained period around an inner 40 d binary consisting of a B3 III star and an unseen companion in a circular orbit. The radial-velocity semi-amplitude of 61.3 km/s of the inner star and its minimum (probable) mass of 5.0 Msun (6.3 +- 0.7 Msun) imply a mass of the unseen object of >= 4.2 Msun (>= 5.0 +- 0.4 Msun), that is, a black hole (BH). The spectroscopic time series is stunningly similar to observations of LB-1. A similar triple-star architecture of LB-1 would reduce the mass of the BH in LB-1 from ~70 Msun to a level more typical of Galactic stellar remnant BHs. The BH in HR 6819 probably is the closest known BH to the Sun, and together with LB-1, suggests a population of quiet BHs. Its embedment in a hierarchical triple structure may be of interest for models of merging double BHs or BH + neutron star binaries. Other triple stars with an outer Be star but without BH are identified; through stripping, such systems may become a source of single Be stars. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: Accepted as Letter to the Editor for A&A

The binary mass ratio in the black hole transient MAXI J1820+070

Authors: M. A. P. Torres, J. Casares, F. Jiménez-Ibarra, A. Álvarez-Hernández, T. Muñoz-Darias, M. Armas Padilla, P. G. Jonker, M. Heida

Abstract: We present intermediate resolution spectroscopy of the optical counterpart to the black hole X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained with the OSIRIS spectrograph on the 10.4-m Gran Telescopio Canarias. The observations were performed with the source close to the quiescent state and before the onset of renewed activity in August 2019. We make use of these data and K-type dwarf templa… ▽ More We present intermediate resolution spectroscopy of the optical counterpart to the black hole X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained with the OSIRIS spectrograph on the 10.4-m Gran Telescopio Canarias. The observations were performed with the source close to the quiescent state and before the onset of renewed activity in August 2019. We make use of these data and K-type dwarf templates taken with the same instrumental configuration to measure the projected rotational velocity of the donor star. We find $v_{rot} \sin i = 84 \pm 5$ km s$^{-1}$ ($1\!-\!σ$), which implies a donor to black-hole mass ratio $q = {M_2}/{M_1} = 0.072 \pm 0.012$ for the case of a tidally locked and Roche-lobe filling donor star. The derived dynamical masses for the stellar components are $M_1 = (5.95 \pm 0.22)\sin ^{-3}i$ $M_\odot$ and $M_2 = (0.43 \pm 0.08) \sin^{-3}i$ $M_\odot$. The use of $q$, combined with estimates of the accretion disk size at the time of the optical spectroscopy, allows us to revise our previous orbital inclination constraints to $66^{\circ} < i < 81^{\circ}$. These values lead to 95% confidence level limits on the masses of $5.73 <M_1(M_\odot) < 8.34$ and $0.28 < M_2(M_\odot) < 0.77$. Adopting instead the $63 \pm 3^{\circ}$ orientation angle of the radio jet as the binary inclination leads to $M_1 = 8.48^{+0.79}_{-0.72} M_\odot$ and $M_2 = 0.61^{+0.13}_{-0.12} M_\odot$ ($1\!-\!σ$). △ Less

Submitted 4 April, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Comments: 8 pages, 2 figures. Accepted for publication in ApJL

Consistency and convergence for a family of finite volume discretizations of the Fokker--Planck operator

Authors: Martin Heida, Markus Kantner, Artur Stephan

Abstract: We introduce a family of various finite volume discretization schemes for the Fokker--Planck operator, which are characterized by different weight functions on the edges. This family particularly includes the well-established Scharfetter--Gummel discretization as well as the recently developed square-root approximation (SQRA) scheme. We motivate this family of discretizations both from the numeric… ▽ More We introduce a family of various finite volume discretization schemes for the Fokker--Planck operator, which are characterized by different weight functions on the edges. This family particularly includes the well-established Scharfetter--Gummel discretization as well as the recently developed square-root approximation (SQRA) scheme. We motivate this family of discretizations both from the numerical and the modeling point of view and provide a uniform consistency and error analysis. Our main results state that the convergence order primarily depends on the quality of the mesh and in second place on the quality of the weights. We show by numerical experiments that for small gradients the choice of the optimal representative of the discretization family is highly non-trivial while for large gradients the Scharfetter--Gummel scheme stands out compared to the others. △ Less

Submitted 21 February, 2020; originally announced February 2020.

MSC Class: 35Q84; 49M25; 65N08

From SN2010da to NGC 300 ULX-1: Ten Years of Observations of an Unusual High Mass X-ray Binary in NGC 300

Authors: Breanna A. Binder, Stefania Carpano, Marianne Heida, Ryan Lau

Abstract: In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and… ▽ More In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and attempts to understand the 2010 transient event and, later, the re-classification of this system as NGC~300 ULX-1: a red supergiant + neutron star ultraluminous X-ray source. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: Accepted to special edition of Galaxies: High-Mass X-ray Binaries and Ultraluminous X-ray Sources Hosting Neutron Stars

Stochastic homogenization on randomly perforated domains

Authors: Martin Heida

Abstract: We study the existence of uniformly bounded extension and trace operators for $W^{1,p}$-functions on randomly perforated domains, where the geometry is assumed to be stationary ergodic. Such extension and trace operators are important for compactness in stochastic homogenization. In contrast to former approaches and results, we use very weak assumptions on the geometry which we call local $(δ,M)$-… ▽ More We study the existence of uniformly bounded extension and trace operators for $W^{1,p}$-functions on randomly perforated domains, where the geometry is assumed to be stationary ergodic. Such extension and trace operators are important for compactness in stochastic homogenization. In contrast to former approaches and results, we use very weak assumptions on the geometry which we call local $(δ,M)$-regularity, isotropic cone mixing and bounded average connectivity. The first concept measures local Lipschitz regularity of the domain while the second measures the mesoscopic distribution of void space. The third is the most tricky part and measures the "mesoscopic" connectivity of the geometry. In contrast to former approaches we do not require a minimal distance between the inclusions and we allow for globally unbounded Lipschitz constants and percolating holes. We will illustrate our method by applying it to the Boolean model based on a Poisson point process and to a Delaunay pipe process. We finally introduce suitable Sobolev spaces on $\mathbb{R}^d$ and $Ω$ in order to construct a stochastic two-scale convergence method and apply the resulting theory to the homogenization of a $p$-Laplace problem on a randomly perforated domain. △ Less

Submitted 21 September, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

MSC Class: 80M40; 60D05

The (re)appearance of NGC 925 ULX-3, a new transient ULX

Authors: Hannah P. Earnshaw, Marianne Heida, Murray Brightman, Felix Fürst, Fiona A. Harrison, Amruta Jaodand, Matthew J. Middleton, Timothy P. Roberts, Rajath Sathyaprakash, Daniel Stern, Dominic J. Walton

Abstract: We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$… ▽ More We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$ erg s$^{-1}$. With an additional Chandra non-detection in 2005, this object demonstrates a high dynamic range of flux of factor >26. In its high-luminosity detections, ULX-3 exhibits a hard power-law spectrum with $Γ=1.6\pm0.1$, whereas the XMM-Newton detection is slightly softer, with $Γ=1.8^{+0.2}_{-0.1}$ and also well-fitted with a broadened disc model. The long-term light curve is sparsely covered and could be consistent either with the propeller effect or with a large-amplitude superorbital period, both of which are seen in ULXs, in particular those with neutron star accretors. Further systematic monitoring of ULX-3 will allow us to determine the mechanism by which ULX-3 undergoes its extreme variability and to better understand the accretion processes of ULXs. △ Less

Submitted 19 February, 2020; v1 submitted 2 January, 2020; originally announced January 2020.

Comments: 9 pages, 3 figures, 1 table. Accepted for publication in ApJ; updated to accepted version with minor revisions

Swift monitoring of M51: A 38-day super-orbital period for the pulsar ULX7 and a new transient ULX

Authors: Murray Brightman, Hannah Earnshaw, Felix Fürst, Fiona A. Harrison, Marianne Heida, Gianluca Israel, Sean Pike, Daniel Stern, Dominic J. Walton

Abstract: We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining th… ▽ More We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: Submitted to ApJ

The Unusual Broadband X-ray Spectral Variability of NGC 1313 X-1 seen with XMM$-$Newton, Chandra and NuSTAR

Authors: D. J. Walton, C. Pinto, M. Nowak, M. Bachetti, R. Sathyaprakash, E. Kara, T. P. Roberts, R. Soria, M. Brightman, C. R. Canizares, H. P. Earnshaw, F. Furst, M. Heida, M. J. Middleton, D. Stern, L. Tao, N. Webb, W. N. Alston, D. Barret, A. C. Fabian, F. A. Harrison, P. Kosec

Abstract: We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X… ▽ More We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability. △ Less

Submitted 21 April, 2020; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: 18 pages, 7 figures, accepted for publication in MNRAS

Discovery of a red supergiant donor star in SN2010da/NGC 300 ULX-1

Authors: M. Heida, R. M. Lau, B. Davies, M. Brightman, F. Fürst, B. W. Grefenstette, J. A. Kennea, F. Tramper, D. J. Walton, F. A. Harrison

Abstract: SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS ste… ▽ More SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $α\approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 10 pages, 6 figures, accepted for publication in ApJ Letters

Discovery and analysis of a ULX nebula in NGC 3521

Authors: K. M. López, P. G. Jonker, M. Heida, M. A. P. Torres, T. P. Roberts, D. J. Walton, D. -S. Moon, F. A. Harrison

Abstract: We present Very Large Telescope/X-shooter and Chandra X-ray observatory/ACIS observations of the ULX [SST2011] J110545.62+000016.2 in the galaxy NGC 3521. The source identified as a candidate near-infrared counterpart to the ULX in our previous study shows an emission line spectrum of numerous recombination and forbidden lines in the visible and near-infrared spectral regime. The emission from the… ▽ More We present Very Large Telescope/X-shooter and Chandra X-ray observatory/ACIS observations of the ULX [SST2011] J110545.62+000016.2 in the galaxy NGC 3521. The source identified as a candidate near-infrared counterpart to the ULX in our previous study shows an emission line spectrum of numerous recombination and forbidden lines in the visible and near-infrared spectral regime. The emission from the candidate counterpart is spatially extended ($\sim$ 34 pc) and appears to be connected with an adjacent H II region, located $\sim$ 138 pc to the NE. The measured velocities of the emission lines confirm that both the candidate counterpart and H II region reside in NGC 3521. The intensity ratios of the emission lines from the ULX counterpart show that the line emission originates from the combined effect of shock and photoionisation of low metallicity (12 + log (O/H) = 8.19 $\pm$ 0.11) gas. Unfortunately, there is no identifiable spectral signature directly related to the photosphere of the mass-donor star in our spectrum. From the archival Chandra data, we derive the X-ray luminosity of the source in the 0.3-7 keV range to be (1.9 $\pm$ 0.8) $\times$ 10$^{40}$ erg cm$^{-2}$ s$^{-1}$, almost a factor of four higher than what is previously reported. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 17 pages, 10 figures, accepted for publication in MNRAS

Dynamical confirmation of a black hole in MAXI J1820+070

Authors: M. A. P. Torres, J. Casares, F. Jiménez-Ibarra, T. Muñoz-Darias, M. Armas-Padilla, P. G. Jonker, M. Heida

Abstract: We present time-resolved 10.4-m GTC and 4.2-m WHT intermediate resolution spectroscopy of the X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained during its decline to the quiescent state. Cross-correlation of the 21 individual spectra against late-type templates reveals a sinusoidal velocity modulation with a period of 0.68549 +/- 0.00001 d and semi-amplitude of 417.7 +/- 3.9 km/s. We derive a… ▽ More We present time-resolved 10.4-m GTC and 4.2-m WHT intermediate resolution spectroscopy of the X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained during its decline to the quiescent state. Cross-correlation of the 21 individual spectra against late-type templates reveals a sinusoidal velocity modulation with a period of 0.68549 +/- 0.00001 d and semi-amplitude of 417.7 +/- 3.9 km/s. We derive a mass function f(M) = 5.18 +/- 0.15 Msun, dynamically confirming the black hole nature of the compact object. Our analysis of the stellar absorption features supports a K3-5 spectral classification for the donor star, which contributes 20% of the total flux at 5200-6800 Angs. The photometric 0.703 +/- 0.003 d periodicity observed during outburst is 2.6% longer than the orbital period supporting the presence of a superhump modulation in the outburst lightcurves. In line with this interpretation, we constrain the binary mass ratio to be q=0.12. In addition, we observe a sharp increase in the Halpha emission line equivalent width during inferior conjunction of the donor star that we interpret as a grazing eclipse of the accretion disc and allows us to constrain the binary inclination to > 69 deg. On the other hand, the absence of X-ray eclipses during outburst imply i < 77 deg. These inclination limits, together with our dynamical solution, lead to a black hole mass in the range 7-8 Msun. We also measure a systemic velocity = -21.6 +/- 2.3 km/s which, combined with the Gaia DR2 proper motion and parallax, implies a large peculiar velocity of 100 km/s. △ Less

Submitted 11 August, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: 7 pages, 3 figures, accepted for publcation in ApJ Letters

On a thermodynamic framework for developing boundary conditions for Korteweg fluids

Authors: Ondřej Souček, Martin Heida, Josef Málek

Abstract: We provide a derivation of several classes of boundary conditions for fluids of Korteweg-type using a simple and transparent thermodynamic approach that automatically guarentees that the derived boundary conditions are compatible with the second law of thermodynamics. The starting assumption of our approach is to describe the boundary of the domain as the membrane separating two different continua… ▽ More We provide a derivation of several classes of boundary conditions for fluids of Korteweg-type using a simple and transparent thermodynamic approach that automatically guarentees that the derived boundary conditions are compatible with the second law of thermodynamics. The starting assumption of our approach is to describe the boundary of the domain as the membrane separating two different continua, one inside the domain, and the other outside the domain. With this viewpoint one may employ the framework of continuum thermodynamics involving singular surfaces. This approach allows us to identify, for various classes of surface Helmholtz free energies, the corresponding surface entropy production mechanisms. By establishing the constitutive relations that guarantee that the surface entropy production is non-negative, we identify a new class of boundary conditions, which on one hand generalizes in a nontrivial manner the Navier's slip boundary conditions, and on the other hand describes dynamic and static contact angle conditions. We explore the general model in detail for a particular case of Korteweg fluid where the Helmholtz free energy in the bulk is that of a van der Waals fluid. We perform a series of numerical experiments to document the basic qualitative features of the novel boundary conditions and their practical applicability to model phenomena such as the contact angle hysteresis. △ Less

Submitted 8 June, 2019; originally announced June 2019.

MSC Class: 76T10

Constraining the nature of the accreting binary in CXOGBS J174623.5-310550

Authors: M. A. P. Torres, S. Repetto, T. Wevers, M. Heida, P. G. Jonker, R. I. Hynes, G. Nelemans, Z. Kostrzewa-Rutkowska, L. Wyrzykowski, C. T. Britt, C. O. Heinke, J. Casares, C. B. Johnson, T. J. Maccarone, D. T. H. Steeghs

Abstract: We report optical and infrared observations of the X-ray source CXOGBS J174623.5-310550. This Galactic object was identified as a potential quiescent low-mass X-ray binary accreting from an M-type donor on the basis of optical spectroscopy and the broad Halpha emission line. The analysis of X-shooter spectroscopy covering 3 consecutive nights supports an M2/3-type spectral classification. Neither… ▽ More We report optical and infrared observations of the X-ray source CXOGBS J174623.5-310550. This Galactic object was identified as a potential quiescent low-mass X-ray binary accreting from an M-type donor on the basis of optical spectroscopy and the broad Halpha emission line. The analysis of X-shooter spectroscopy covering 3 consecutive nights supports an M2/3-type spectral classification. Neither radial velocity variations nor rotational broadening is detected in the photospheric lines. No periodic variability is found in I- and r'-band light curves. We derive r' = 20.8, I = 19.2 and Ks = 16.6 for the optical and infrared counterparts with the M-type star contributing 90% to the I-band light. We estimate its distance to be 1.3-1.8 kpc. The lack of radial velocity variations implies that the M-type star is not the donor star in the X-ray binary. This could be an interloper or the outer body in a hierarchical triple. We constrain the accreting binary to be a < 2.2 hr orbital period eclipsing cataclysmic variable or a low-mass X-ray binary lying in the foreground of the Galactic Bulge. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: (9 pages, 5 figures, accepted for publication in MNRAS)

All at once: transient pulsations, spin down and a glitch from the Pulsating Ultraluminous X-ray Source M82 X-2

Authors: Matteo Bachetti, Thomas J. Maccarone, Murray Brightman, McKinley C. Brumback, Felix Fürst, Fiona A. Harrison, Marianne Heida, Gian Luca Israel, Matthew J. Middleton, John A. Tomsick, Natalie A. Webb, Dominic J. Walton

Abstract: M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down… ▽ More M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down rate between 2014 and 2016 of 5x10^-11 Hz/s, which is in contrast to the strong spin up seen during the 2014 observations and represent the first detection of spin down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 dataset. We find a glitch superimposed on the very strong and variable spin-up already reported, the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs. △ Less

Submitted 19 November, 2019; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 15 pages, 7 figures, 5 tables. Revised version

A broadband look at the old and new ULXs of NGC 6946

Authors: Hannah P. Earnshaw, Brian W. Grefenstette, Murray Brightman, Dominic J. Walton, Didier Barret, Felix Fürst, Fiona A. Harrison, Marianne Heida, Sean N. Pike, Daniel Stern, Natalie A. Webb

Abstract: Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $Γ=3.6^{+0.4}_{-0.3}$ in both cas… ▽ More Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $Γ=3.6^{+0.4}_{-0.3}$ in both cases. Their properties are consistent with being super-Eddington accreting sources with the majority of their hard emission obscured and down-scattered. ULX-3 (NGC 6946 X-1) is significantly detected by both XMM-Newton and NuSTAR at $L_{\rm X}=(6.5\pm0.1)\times10^{39}$ erg s$^{-1}$, and has a power-law spectrum with $Γ=2.51\pm0.05$. We are unable to identify a high-energy break in its spectrum like that found in other ULXs, but the soft spectrum likely hinders our ability to detect one. We also characterise the new source, ULX-4, which is only detected in the joint XMM-Newton and NuSTAR observation, at $L_{\rm X}=(2.27\pm0.07)\times10^{39}$ erg s$^{-1}$, and is absent in a Chandra observation ten days later. It has a very hard cut-off power-law spectrum with $Γ=0.7\pm0.1$ and $E_{\rm cut}=11^{+9}_{-4}$ keV. We do not detect pulsations from ULX-4, but its transient nature can be explained either as a neutron star ULX briefly leaving the propeller regime or as a micro-tidal disruption event induced by a stellar-mass compact object. △ Less

Submitted 8 May, 2019; originally announced May 2019.

Comments: 19 pages, 11 figures, 6 tables. Accepted for publication by ApJ

Stochastic homogenization of $Λ$-convex gradient flows

Authors: Martin Heida, Stefan Neukamm, Mario Varga

Abstract: In this paper we present a stochastic homogenization result for a class of Hilbert space evolutionary gradient systems driven by a quadratic dissipation potential and a $Λ$-convex energy functional featuring random and rapidly oscillating coefficients. Specific examples included in the result are Allen-Cahn type equations and evolutionary equations driven by the $p$-Laplace operator with… ▽ More In this paper we present a stochastic homogenization result for a class of Hilbert space evolutionary gradient systems driven by a quadratic dissipation potential and a $Λ$-convex energy functional featuring random and rapidly oscillating coefficients. Specific examples included in the result are Allen-Cahn type equations and evolutionary equations driven by the $p$-Laplace operator with $p\in (1,\infty)$. The homogenization procedure we apply is based on a stochastic two-scale convergence approach. In particular, we define a stochastic unfolding operator which can be considered as a random counterpart of the well-established notion of periodic unfolding. The stochastic unfolding procedure grants a very convenient method for homogenization problems defined in terms of ($Λ$-)convex functionals. △ Less

Submitted 5 May, 2019; originally announced May 2019.

Comments: arXiv admin note: text overlap with arXiv:1805.09546

MSC Class: 49J40; 74Q10; 35K57

A Mathematical model for Alzheimer's disease: An approach via stochastic homogenization of the Smoluchowski equation

Authors: Bruno Franchi, Martin Heida, Silvia Lorenzani

Abstract: In this note, we apply the theory of stochastic homogenization to find the asymptotic behavior of the solution of a set of Smoluchowski's coagulation-diffusion equations with non-homogeneous Neumann boundary conditions. This system is meant to model the aggregation and diffusion of $β$-amyloid peptide (A$β$) in the cerebral tissue, a process associated with the development of Alzheimer's disease.… ▽ More In this note, we apply the theory of stochastic homogenization to find the asymptotic behavior of the solution of a set of Smoluchowski's coagulation-diffusion equations with non-homogeneous Neumann boundary conditions. This system is meant to model the aggregation and diffusion of $β$-amyloid peptide (A$β$) in the cerebral tissue, a process associated with the development of Alzheimer's disease. In contrast to the approach used in our previous works, in the present paper we account for the non-periodicity of the cellular structure of the brain by assuming a stochastic model for the spatial distribution of neurons. Further, we consider non-periodic random diffusion coefficients for the amyloid aggregates and a random production of A$β$ in the monomeric form at the level of neuronal membranes. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Uncovering Red and Dusty Ultraluminous X-ray Sources with Spitzer

Authors: Ryan M. Lau, Marianne Heida, Dominic J. Walton, Mansi M. Kasliwal, Scott M. Adams, Ann Marie Cody, Kishalay De, Robert D. Gehrz, Felix Furst, Jacob E. Jencson, jamie A. Kennea, Frank Masci

Abstract: We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candida… ▽ More We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $α=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 40 pages, 11 figures, accepted for publication in ApJ

The fractional p-Laplacian emerging from homogenization of the random conductance model with degenerate ergodic weights and unbounded-range jumps

Authors: Franziska Flegel, Martin Heida

Abstract: We study a general class of discrete $p$-Laplace operators in the random conductance model with long-range jumps and ergodic weights. Using a variational formulation of the problem, we show that under the assumption of bounded first moments and a suitable lower moment condition on the weights, the homogenized limit operator is a fractional $p$-Laplace operator. Under strengthened lower moment co… ▽ More We study a general class of discrete $p$-Laplace operators in the random conductance model with long-range jumps and ergodic weights. Using a variational formulation of the problem, we show that under the assumption of bounded first moments and a suitable lower moment condition on the weights, the homogenized limit operator is a fractional $p$-Laplace operator. Under strengthened lower moment conditions, we can apply our insights also to the spectral homogenization of the discrete Laplace operator to the continuous fractional Laplace operator. △ Less

Submitted 15 April, 2019; originally announced April 2019.

MSC Class: 80M40; 60H25; 60K37; 35B27; 35R60; 47B80; 47A75

Astro 2020 Science White Paper: Time Domain Studies of Neutron Star and Black Hole Populations: X-ray Identification of Compact Object Types

Authors: N. Vulic, A. E. Hornschemeier, V. Antoniou, A. R. Basu-Zych, B. Binder, F. M. Fornasini, F. Furst, F. Haberl, M. Heida, B. D. Lehmer, T. J. Maccarone, A. F. Ptak, G. R. Sivakoff, P. Tzanavaris, D. R. Wik, B. F. Williams, J. Wilms, M. Yukita, A. Zezas

Abstract: What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, a… ▽ More What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates. △ Less

Submitted 23 March, 2019; originally announced March 2019.

Comments: 9 pages, 2 figures. Submitted to the Astro2020 Decadal Survey

Breaking the limit: Super-Eddington accretion onto black holes and neutron stars

Authors: M. Brightman, M. Bachetti, H. P. Earnshaw, F. Fürst, J. García, B. Grefenstette, M. Heida, E. Kara, K. K. Madsen, M. J. Middleton, D. Stern, F. Tombesi, D. J. Walton

Abstract: With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (V… ▽ More With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (Volonteri & Rees, 2005). Indeed, for any stellar remnant seed, the rate required would need to exceed the Eddington limit. This is the point at which the outward force produced by radiation pressure is equal to the gravitational attraction experienced by the in-falling matter. In principle, this implies that there is a maximum luminosity an object of mass $M$ can emit; assuming spherical accretion and that the opacity is dominated by Thompson scattering, this Eddington luminosity is $L_{\rm{E}} = 1.38 \times 10^{38} (M/M_{\odot})$ erg s$^{-1}$. In reality, it is known that this limit can be violated, due to non-spherical geometry or various kinds of instabilities. Nevertheless, the Eddington limit remains an important reference point, and many of the details of how accretion proceeds above this limit remain unclear. Understanding how this so-called super-Eddington accretion occurs is of clear cosmological importance, since it potentially governs the growth of the first supermassive black holes (SMBHs) and the impact this growth would have had on their host galaxies (`feedback') and the epoch of reionization, as well as improving our understanding of accretion physics more generally. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: A Science White Paper submitted to the Astro2020 Decadal Survey

Searching for the Donor Stars of ULX Pulsars

Authors: M. Heida, F. A. Harrison, M. Brightman, F. Fürst, D. Stern, D. J. Walton

Abstract: We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival HST observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3… ▽ More We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival HST observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3/UVIS F336W image whose photometry matches that of 10-15 M$_\odot$ stars turning off the main sequence. Such stars have densities consistent with the properties of the donor star of M82 X-2 as inferred from X-ray timing analysis, although it is also possible that the donor is a lower mass star below our detection limit or that there is a significant contribution from the accretion disc to the optical emission. We detect three candidate counterparts to NGC 5907 X-1 in the near-infrared. All of these are too bright to be the donor star of the ULX, which based on its orbital period is a red giant. The high background at the location of NGC 5907 X-1 precludes us from detecting this expected donor star. The recently discovered NGC 5907 ULX-2 also falls within the field of view of the near-infrared imaging; we detect four sources in the error circle, with photometry that matches AGB stars. The star suggested to be the counterpart of NGC 5907 ULX-2 by Pintore et al. (2018) falls outside our 2-$σ$ error circle. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: 10 pages, 3 figures, accepted for publication in ApJ

A new class of flares from accreting supermassive black holes

Authors: Benny Trakhtenbrot, Iair Arcavi, Claudio Ricci, Sandro Tacchella, Daniel Stern, Hagai Netzer, Peter G. Jonker, Assaf Horesh, Julián Esteban Mejía-Restrepo, Griffin Hosseinzadeh, Valentina Hallefors, D. Andrew Howell, Curtis McCully, Mislav Baloković, Marianne Heida, Nikita Kamraj, George Benjamin Lansbury, Łukasz Wyrzykowski, Mariusz Gromadzki, Aleksandra Hamanowicz, S. Bradley Cenko, David J. Sand, Eric Y. Hsiao, Mark M. Phillips, Tiara R. Diamond , et al. (4 additional authors not shown)

Abstract: Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed… ▽ More Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: Published in Nature Astronomy

Quiescent NIR and optical counterparts to candidate black hole X-ray binaries

Authors: K. M. López, P. G. Jonker, M. A. P. Torres, M. Heida, A. Rau, D. Steeghs

Abstract: We present near-infrared and optical imaging of fifteen candidate black hole X-ray binaries. In addition to quiescent observations for all sources, we also observed two of these sources (IGR J17451-3022 and XTE J1818-245) in outburst. We detect the quiescent counterpart for twelve out of fifteen sources, and for the remaining three we report limiting magnitudes. The magnitudes of the detected coun… ▽ More We present near-infrared and optical imaging of fifteen candidate black hole X-ray binaries. In addition to quiescent observations for all sources, we also observed two of these sources (IGR J17451-3022 and XTE J1818-245) in outburst. We detect the quiescent counterpart for twelve out of fifteen sources, and for the remaining three we report limiting magnitudes. The magnitudes of the detected counterparts range between $K_s$ = 17.59 and $K_s$ = 22.29 mag. We provide (limits on) the absolute magnitudes and finding charts of all sources. Of these twelve detections in quiescence, seven represent the first quiescent reported values (for MAXI J1543-564, XTE J1726-476, IGR J17451-3022, XTE J1818-245, MAXI J1828-249, MAXI J1836-194, Swift J1910.2-0546) and two detections show fainter counterparts to XTE J1752-223 and XTE J2012+381 than previously reported. We used theoretical arguments and observed trends, for instance between the outburst and quiescent X-ray luminosity and orbital period $P_{orb}$ to derive an expected trend between $ΔK_s$ and $P_{orb}$ of $ΔK_s \propto \log P_{orb}^{0.565}$. Comparing this to observations we find a different behaviour. We discuss possible explanations for this result. △ Less

Submitted 14 October, 2018; originally announced October 2018.

Comments: 18 pages, 6 figures. Accepted for publication in MNRAS

The Fast, Luminous Ultraviolet Transient AT2018cow: Extreme Supernova, or Disruption of a Star by an Intermediate-Mass Black Hole?

Authors: Daniel A. Perley, Paolo A. Mazzali, Lin Yan, S. Bradley Cenko, Suvi Gezari, Kirsty Taggart, Nadia Blagorodnova, Christoffer Fremling, Brenna Mockler, Avinash Singh, Nozomu Tominaga, Masaomi Tanaka, Alan M. Watson, Tomás Ahumada, G. C. Anupama, Chris Ashall, Rosa L. Becerra, David Bersier, Varun Bhalerao, Joshua S. Bloom, Nathaniel R. Butler, Chris Copperwheat, Michael W. Coughlin, Kishalay De, Andrew J. Drake , et al. (38 additional authors not shown)

Abstract: Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwi… ▽ More Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins. △ Less

Submitted 23 November, 2018; v1 submitted 2 August, 2018; originally announced August 2018.

Comments: Corrected Figure 8 / Table 4 to use final fits. Includes machine-readable photometry table (hopefully for real this time)

Stochastic unfolding and homogenization

Authors: Martin Heida, Stefan Neukamm, Mario Varga

Abstract: The notion of periodic two-scale convergence and the method of periodic unfolding are prominent and useful tools in multiscale modeling and analysis of PDEs with rapidly oscillating periodic coefficients. In this paper we are interested in the theory of stochastic homogenization for continuum mechanical models in form of PDEs with random coefficients, describing random heterogeneous materials. The… ▽ More The notion of periodic two-scale convergence and the method of periodic unfolding are prominent and useful tools in multiscale modeling and analysis of PDEs with rapidly oscillating periodic coefficients. In this paper we are interested in the theory of stochastic homogenization for continuum mechanical models in form of PDEs with random coefficients, describing random heterogeneous materials. The notion of periodic two-scale convergence has been extended in different ways to the stochastic case. In this work we introduce a stochastic unfolding method that features many similarities to periodic unfolding. In particular it allows to characterize the notion of stochastic two-scale convergence in the mean by mere weak convergence in an extended space. We illustrate the method on the (classical) example of stochastic homogenization of convex integral functionals, and prove a new result on stochastic homogenization for a non-convex evolution equation of Allen-Cahn type. Moreover, we discuss the relation of stochastic unfolding to previously introduced notions of (quenched and mean) stochastic two-scale convergence. The method described in the present paper extends to the continuum setting the notion of discrete stochastic unfolding, as recently introduced by the second and third author in the context of discrete-to-continuum transition. △ Less

Submitted 27 May, 2021; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: 46 pages. This preprint has been split into two parts. For streamlined and extended versions of these parts see arXiv:2105.12447 and arXiv:1905.02562. Please refer to those

MSC Class: 49J40; 74Q05; 74Q10; 35K57