Showing 1–50 of 253 results for author: Nowak, M

Multistructured accretion flow of Sgr A* II: Signatures of a Cool Accretion Disk in Hydrodynamic Simulations of Stellar Winds

Authors: Mayura Balakrishnan, Christopher M. P. Russell, Lia Corrales, Diego Calderón, Jorge Cuadra, Daryl Haggard, Sera Markoff, Joey Neilsen, Michael Nowak, Q. Daniel Wang, Fred Baganoff

Abstract: Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulat… ▽ More Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulation runtime of 3500 years. This result is consistent with previous grid-based simulations, demonstrating that a cold disk can form regardless of numerical method. We examine the plasma scenarios arising from an environment with and without this cold disk, by generating synthetic spectra for comparison to the quiescent Fe K alpha Sgr A* spectrum from Chandra HETG-S, taken through the Chandra X-ray Visionary Program. We find that current and future X-ray missions are unlikely to distinguish between the kinematic signatures in the plasma in these two scenarios. Nonetheless, the stellar wind plasma model presents a good fit to the dispersed Chandra spectra within 1.5" of Sgr A*. We compare our results to the Radiatively Inefficient Accretion Flow (RIAF) model fit to the HETG-S spectrum presented in Paper I and find that the Bayesian model evidence does not strongly favor either model. With 9" angular resolution and high spectral resolution of the X-IFU, NewAthena will offer a clearer differentiation between the RIAF plasma model and hydrodynamic simulations, but only a future X-ray mission with arcsecond resolution will significantly advance our understanding of Sgr A*'s accretion flow in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ

Multistructured accretion flow of Sgr A* I: Examination of a RIAF model

Authors: Mayura Balakrishnan, Lia Corrales, Sera Markoff, Michael Nowak, Daryl Haggard, Q. Daniel Wang, Joey Neilsen, Christopher M. P. Russell, Diego Calderón, Jorge Cuadra, Fred Baganoff

Abstract: The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Progr… ▽ More The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution gratings spectrum indicate an inflow balanced by an outflow ($s \sim 1$) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following $1/r$. The data require that the abundance of Iron relative to solar is $Z_{Fe} < 0.32 Z_\odot$ (90\% credible interval), much lower than the $2~Z_\odot$ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra's unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic Center in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures, 1 table. Accepted for publication in ApJ

Exoplanets in reflected starlight with dual-field interferometry: A case for shorter wavelengths and a fifth Unit Telescope at VLTI/Paranal

Authors: S. Lacour, Ó. Carrión-González, M. Nowak

Abstract: The direct observation of cold and temperate planets within 1 to 10 AU would be extremely valuable for uncovering their atmospheric compositions but remains a formidable challenge with current astronomical methods. Ground-based optical interferometry, capable of high angular-resolution imaging, offers a promising avenue for studying these exoplanets, complementing space-based observations. Our obj… ▽ More The direct observation of cold and temperate planets within 1 to 10 AU would be extremely valuable for uncovering their atmospheric compositions but remains a formidable challenge with current astronomical methods. Ground-based optical interferometry, capable of high angular-resolution imaging, offers a promising avenue for studying these exoplanets, complementing space-based observations. Our objective is to explore the fundamental limits of dual-field interferometry and assess its potential for characterizing exoplanets in reflected light using the Very Large Telescope Interferometer (VLTI). We developed analytical expressions to describe the performance of dual-field interferometry and integrated these with simulations of atmospheric wavefronts corrected by extreme Adaptive Optics. An analytical solution for optimal phase apodization was formulated to enhance starlight rejection when injected into a single-mode fibre. This framework was applied to determine the detectability of known exoplanets in reflected light across various wavelength bands for both the current VLTI and a proposed extended version. Our results indicate that employing shorter wavelengths improves detectability, enabling at least seven Jupiter-mass exoplanets to be observed in the J band with current VLTI's baselines. Adding new baselines with lengths beyond 200 meters significantly enhances VLTI's capabilities, increasing the number of detectable exoplanets and revealing potential habitable zone candidates such as $τ$ Ceti e and Proxima Centauri b. To substantially improve the VLTI's exoplanet characterization capabilities, we recommend developing instrumentation at wavelengths shorter than 1$\,μ$m, as well as the addition of a fifth Unit Telescope (UT5). △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: To be discussed with the community. Please contact us if you wish to send feedback. Submitted to A&A. 15 pages plus appendix

The $β$ Pictoris b Hill sphere transit campaign. Paper II: Searching for the signatures of the $β$ Pictoris exoplanets through time delay analysis of the $δ$ Scuti pulsations

Authors: Sebastian Zieba, Konstanze Zwintz, Matthew Kenworthy, Daniel Hey, Simon J. Murphy, Rainer Kuschnig, Lyu Abe, Abdelkrim Agabi, Djamel Mekarnia, Tristan Guillot, François-Xavier Schmider, Philippe Stee, Yuri De Pra, Marco Buttu, Nicolas Crouzet, Samuel Mellon, Jeb Bailey III, Remko Stuik, Patrick Dorval, Geert-Jan J. Talens, Steven Crawford, Eric Mamajek, Iva Laginja, Michael Ireland, Blaine Lomberg , et al. (12 additional authors not shown)

Abstract: The $β$ Pictoris system is the closest known stellar system with directly detected gas giant planets, an edge-on circumstellar disc, and evidence of falling sublimating bodies and transiting exocomets. The inner planet, $β$ Pictoris c, has also been indirectly detected with radial velocity (RV) measurements. The star is a known $δ$ Scuti pulsator, and the long-term stability of these pulsations op… ▽ More The $β$ Pictoris system is the closest known stellar system with directly detected gas giant planets, an edge-on circumstellar disc, and evidence of falling sublimating bodies and transiting exocomets. The inner planet, $β$ Pictoris c, has also been indirectly detected with radial velocity (RV) measurements. The star is a known $δ$ Scuti pulsator, and the long-term stability of these pulsations opens up the possibility of indirectly detecting the gas giant planets through time delays of the pulsations due to a varying light travel time. We search for phase shifts in the $δ$ Scuti pulsations consistent with the known planets $β$ Pictoris b and c and carry out an analysis of the stellar pulsations of $β$ Pictoris over a multi-year timescale. We used photometric data collected by the BRITE-Constellation, bRing, ASTEP, and TESS to derive a list of the strongest and most significant $δ$ Scuti pulsations. We carried out an analysis with the open-source python package maelstrom to study the stability of the pulsation modes of $β$ Pictoris in order to determine the long-term trends in the observed pulsations. We did not detect the expected signal for $β$ Pictoris b or $β$ Pictoris c. The expected time delay is 6 seconds for $β$ Pictoris c and 24 seconds for $β$ Pictoris b. With simulations, we determined that the photometric noise in all the combined data sets cannot reach the sensitivity needed to detect the expected timing drifts. An analysis of the pulsational modes of $β$ Pictoris using maelstrom showed that the modes themselves drift on the timescale of a year, fundamentally limiting our ability to detect exoplanets around $β$ Pictoris via pulsation timing. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 16 pages, 16 figures, 4 tables, accepted for publication in A&A

High contrast at short separation with VLTI/GRAVITY: Bringing Gaia companions to light

Authors: N. Pourré, T. O. Winterhalder, J. -B. Le Bouquin, S. Lacour, A. Bidot, M. Nowak, A. -L. Maire, D. Mouillet, C. Babusiaux, J. Woillez, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube , et al. (151 additional authors not shown)

Abstract: Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working… ▽ More Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($Δ\mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($Δ\mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 16 pages, 14 figures. Submitted to A&A

GRAVITY for MATISSE -- Improving the MATISSE performance with the GRAVITY fringe tracker

Authors: J. Woillez, R. Petrov, R. Abuter, F. Allouche, P. Berio, R. Dembet, F. Eisenhauer, R. Frahm, F. Gonté, X. Haubois, M. Houllé, W. Jaffe, S. Lacour, S. Lagarde, J. Leftley, B. Lopez, A. Matter, A. Meilland, F. Millour, M. Nowak, C. Paladini, T. Rivinius, D. Salabert, N. Schuhler, J. Varga , et al. (1 additional authors not shown)

Abstract: Context: MATISSE, the mid-infrared spectro-imaging instrument of VLTI, was designed to deliver its advertised performance when paired with an external second generation fringe tracker. Science observation started in 2019, demonstrating imaging capabilities and faint science target observations. Now, The GRAVITY fringe tracker stabilizes the MATISSE fringes which allows using all spectroscopic mode… ▽ More Context: MATISSE, the mid-infrared spectro-imaging instrument of VLTI, was designed to deliver its advertised performance when paired with an external second generation fringe tracker. Science observation started in 2019, demonstrating imaging capabilities and faint science target observations. Now, The GRAVITY fringe tracker stabilizes the MATISSE fringes which allows using all spectroscopic modes and improves sensitivity and data accuracy. Aims: We present how the MATISSE and GRAVITY instruments were adapted to make the GRAVITY fringe tracker work with MATISSE, under the umbrella of the aptly-named GRA4MAT project, led by ESO in collaboration with the two instrument consortia. Methods: We detail the software modifications needed to implement an acquisition and observing sequence specific to GRA4MAT, including simultaneous fringe tracking and chopping and a narrow off-axis capability inspired by the galactic center and exoplanet capability of GRAVITY. We explain the modified data collection and reduction processes. We show how we leveraged the recent fringe tracker upgrade to implement features specific to its use with MATISSE, e.g. fringe jumps mitigation with an improved group delay control and simultaneous fringe tracking and chopping with a new state machine. Results: We successfully demonstrate significant improvements to the MATISSE instrument. Observations can now be performed at higher spectral resolutions of up to $R\sim3300$ and across the full LM bands at once. Long detector integration times, made possible with stabilized fringes, have improved the LM-bands sensitivity by a factor of 10. Low flux biases in coherently-reduced N-band data have been eliminated. The L-band transfer function is now higher and more stable. We finally illustrate the scientific potential of GRA4MAT with a preview of the first exoplanet observation made by MATISSE on $β$ Pictoris b. △ Less

Submitted 30 June, 2024; v1 submitted 31 May, 2024; originally announced May 2024.

Comments: Revised after language editing

Long term variability of Cygnus X-1. VIII. A spectral-timing look at low energies with NICER

Authors: Ole König, Guglielmo Mastroserio, Thomas Dauser, Mariano Méndez, Jingyi Wang, Javier A. García, James F. Steiner, Katja Pottschmidt, Ralf Ballhausen, Riley M. Connors, Federico García, Victoria Grinberg, David Horn, Adam Ingram, Erin Kara, Timothy R. Kallman, Matteo Lucchini, Edward Nathan, Michael A. Nowak, Philipp Thalhammer, Michiel van der Klis, Jörn Wilms

Abstract: The Neutron Star Interior Composition Explorer (NICER) monitoring campaign of Cyg X-1 allows us to study its spectral-timing behavior at energies ${<}1$ keV across all states. The hard state power spectrum can be decomposed into two main broad Lorentzians with a transition at around 1 Hz. The lower-frequency Lorentzian is the dominant component at low energies. The higher-frequency Lorentzian begi… ▽ More The Neutron Star Interior Composition Explorer (NICER) monitoring campaign of Cyg X-1 allows us to study its spectral-timing behavior at energies ${<}1$ keV across all states. The hard state power spectrum can be decomposed into two main broad Lorentzians with a transition at around 1 Hz. The lower-frequency Lorentzian is the dominant component at low energies. The higher-frequency Lorentzian begins to contribute significantly to the variability above 1.5 keV and dominates at high energies. We show that the low- and high-frequency Lorentzians likely represent individual physical processes. The lower-frequency Lorentzian can be associated with a (possibly Comptonized) disk component, while the higher-frequency Lorentzian is clearly associated with the Comptonizing plasma. At the transition of these components, we discover a low-energy timing phenomenon characterized by an abrupt lag change of hard (${\gtrsim}2$ keV) with respect to soft (${\lesssim}1.5$ keV) photons, accompanied by a drop in coherence, and a reduction in amplitude of the second broad Lorentzian. The frequency of the phenomenon increases with the frequencies of the Lorentzians as the source softens and cannot be seen when the power spectrum is single-humped. A comparison to transient low-mass X-ray binaries shows that this feature does not only appear in Cyg X-1, but that it is a general property of accreting black hole binaries. In Cyg X-1, we find that the variability at low and high energies is overall highly coherent in the hard and intermediate states. The high coherence shows that there is a process at work which links the variability, suggesting a physical connection between the accretion disk and Comptonizing plasma. This process fundamentally changes in the soft state, where strong red noise at high energies is incoherent to the variability at low energies. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 18+29 pages, 17+54 figures, accepted for publication in Astronomy & Astrophysics

SPHERE RefPlanets: Search for epsilon Eridani b and warm dust

Authors: C. Tschudi, H. M. Schmid, M. Nowak, H. Le Coroller, S. Hunziker, R. G. van Holstein, C. Perrot, D. Mouillet, J. -C. Augereau, A. Bazzon, J. L. Beuzit, A. Boccaletti, M. J. Bonse, G. Chauvin, S. Desidera, K. Dohlen, C. Dominik, N. Engler, M. Feldt, J. H. Girard, R. Gratton, Th. Henning, M. Kasper, P. Kervella, A. -M. Lagrange , et al. (13 additional authors not shown)

Abstract: We carried out very deep VLT/SPHERE imaging polarimetry of the nearby system Eps Eri based on 38.5 hours of integration time with a 600 - 900 nm broadband filter to search for polarized scattered light from a planet or from circumstellar dust using AO, coronagraphy, high precision differential polarimetry, and angular differential imaging. We have improved several data reduction and post-processin… ▽ More We carried out very deep VLT/SPHERE imaging polarimetry of the nearby system Eps Eri based on 38.5 hours of integration time with a 600 - 900 nm broadband filter to search for polarized scattered light from a planet or from circumstellar dust using AO, coronagraphy, high precision differential polarimetry, and angular differential imaging. We have improved several data reduction and post-processing techniques and also developed new ones to further increase the sensitivity of SPHERE/ZIMPOL. The data provide unprecedented contrast limits, but no significant detection of a point source or an extended signal from circumstellar dust. For each observing epoch, we obtained a point source contrast for the polarized intensity between $2\cdot 10^{-8}$ and $4\cdot 10^{-8}$ at the expected separation of the planet Eps Eri b of 1'' near quadrature phase. The polarimetric contrast limits are about six to 50 times better than the intensity limits because polarimetric imaging is much more efficient in speckle suppression. Combining the entire 14-month data set to the search for a planet moving on a Keplerian orbit with the K-Stacker software further improves the contrast limits by a factor of about two, to about $8 \cdot 10^{-9}$ at 1''. This would allow the detection of a planet with a radius of about 2.5 Jupiter radii. The surface brightness contrast limits achieved for the polarized intensity from an extended scattering region are about 15 mag arcsec$^{-2}$ at 1'', or up to 3 mag arcsec$^{-2}$ deeper than previous limits. For Eps Eri, these limits exclude the presence of a narrow dust ring and they constrain the dust properties. This study shows that the polarimetric contrast limits for reflecting planets with SPHERE/ZIMPOL can be improved to a level $<10^{-8}$ simply by collecting more data over many nights and using the K-Stacker software. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 26 pages, 20 figures, 5 tables. Accepted for publication in A&A

Chandra Study of the Proper Motion of HST-1 in the Jet of M87

Authors: Rameshan Thimmappa, Joey Neilsen, Daryl Haggard, Michael A. Nowak, Sera Markoff

Abstract: The radio galaxy M87 is well known for its jet, which features a series of bright knots observable from radio to X-ray wavelengths. We analyze the X-ray image and flux variability of the knot HST-1 in the jet. Our analysis includes all 112 available Chandra ACIS-S observations from 2000-2021, with a total exposure time of $\sim$887 ks. We use de-convolved images to study the brightness profile of… ▽ More The radio galaxy M87 is well known for its jet, which features a series of bright knots observable from radio to X-ray wavelengths. We analyze the X-ray image and flux variability of the knot HST-1 in the jet. Our analysis includes all 112 available Chandra ACIS-S observations from 2000-2021, with a total exposure time of $\sim$887 ks. We use de-convolved images to study the brightness profile of the X-ray jet and measure the relative separation between the core and HST-1. From 2003-2005 (which coincides with a bright flare from HST-1), we find a correlation between the flux of HST-1 and its offset from the core. In subsequent data, we find a steady increase in this offset, which implies a bulk superluminal motion for HST-1 of 6.6$\pm$0.9 c (2.0$\pm$0.3 pc yr$^{-1}$), in keeping with prior results. We discuss models for the flux-offset correlation that feature either two or four emission regions separated by tens of parsecs. We attribute these results to moving shocks in the jet, that allow us to measure the internal structure of the jet. △ Less

Submitted 1 May, 2024; v1 submitted 30 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ. 9 pages, 13 figures, 1 table

Astrometric detection of a Neptune-mass candidate planet in the nearest M-dwarf binary system GJ65 with VLTI/GRAVITY

Authors: GRAVITY Collaboration, R. Abuter, A. Amorim, M. Benisty, J-P. Berger, H. Bonnet, G. Bourdarot, P. Bourget, W. Brandner, Y. Clénet, R. Davies, F. Delplancke-Ströbele, R. Dembet, A. Drescher, A. Eckart, F. Eisenhauer, H. Feuchtgruber, G. Finger, N. M. Förster-Schreiber, P. Garcia, R. Garcia-Lopez, F. Gao, E. Gendron, R. Genzel, S. Gillessen , et al. (43 additional authors not shown)

Abstract: The detection of low-mass planets orbiting the nearest stars is a central stake of exoplanetary science, as they can be directly characterized much more easily than their distant counterparts. Here, we present the results of our long-term astrometric observations of the nearest binary M-dwarf Gliese 65 AB (GJ65), located at a distance of only 2.67 pc. We monitored the relative astrometry of the tw… ▽ More The detection of low-mass planets orbiting the nearest stars is a central stake of exoplanetary science, as they can be directly characterized much more easily than their distant counterparts. Here, we present the results of our long-term astrometric observations of the nearest binary M-dwarf Gliese 65 AB (GJ65), located at a distance of only 2.67 pc. We monitored the relative astrometry of the two components from 2016 to 2023 with the VLTI/GRAVITY interferometric instrument. We derived highly accurate orbital parameters for the stellar system, along with the dynamical masses of the two red dwarfs. The GRAVITY measurements exhibit a mean accuracy per epoch of 50-60 microarcseconds in 1.5h of observing time using the 1.8m Auxiliary Telescopes. The residuals of the two-body orbital fit enable us to search for the presence of companions orbiting one of the two stars (S-type orbit) through the reflex motion they imprint on the differential A-B astrometry. We detected a Neptune-mass candidate companion with an orbital period of p = 156 +/- 1 d and a mass of m = 36 +/- 7 Mearth. The best-fit orbit is within the dynamical stability region of the stellar pair. It has a low eccentricity, e = 0.1 - 0.3, and the planetary orbit plane has a moderate-to-high inclination of i > 30° with respect to the stellar pair, with further observations required to confirm these values. These observations demonstrate the capability of interferometric astrometry to reach microarcsecond accuracy in the narrow-angle regime for planet detection by reflex motion from the ground. This capability offers new perspectives and potential synergies with Gaia in the pursuit of low-mass exoplanets in the solar neighborhood. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: Corresponding authors: G.Bourdarot, P.Kervella, O.Pfuhl. Accepted in A&A Letters

Four-of-a-kind? Comprehensive atmospheric characterisation of the HR 8799 planets with VLTI/GRAVITY

Authors: E. Nasedkin, P. Mollière, S. Lacour, M. Nowak, L. Kreidberg, T. Stolker, J. J. Wang, W. O. Balmer, J. Kammerer, J. Shangguan, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli , et al. (73 additional authors not shown)

Abstract: With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels… ▽ More With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5σ$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 45 pages, 25 figures, submitted to Astronomy & Astrophysics

Combining Gaia and GRAVITY: Characterising five new Directly Detected Substellar Companions

Authors: T. O. Winterhalder, S. Lacour, A. Mérand, A. -L. Maire, J. Kammerer, T. Stolker, N. Pourré, C. Babusiaux, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (74 additional authors not shown)

Abstract: Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observati… ▽ More Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+. △ Less

Submitted 24 June, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

Rapid Variability of Mrk 421 During Extreme Flaring as Seen Through the Eyes of XMM-Newton

Authors: A. Gokus, J. Wilms, M. Kadler, D. Dorner, M. A. Nowak, A. Kreikenbohm, K. Leiter, T. Bretz, B. Schleicher, A. G. Markowitz, K. Pottschmidt, K. Mannheim, I. Kreykenbohm, M. Langejahn, F. McBride, T. Beuchert, T. Dauser, M. Kreter, J. Abhir, D. Baack, M. Balbo, A. Biland, K. Brand, J. Buss, L. Eisenberger , et al. (9 additional authors not shown)

Abstract: By studying the variability of blazars across the electromagnetic spectrum, it is possible to resolve the underlying processes responsible for rapid flux increases, so-called flares. We report on an extremely bright X-ray flare in the high-peaked BL Lacertae object Mrk 421 that occurred simultaneously with enhanced $γ$-ray activity detected at very high energies (VHE) by FACT on 2019 June 9. We tr… ▽ More By studying the variability of blazars across the electromagnetic spectrum, it is possible to resolve the underlying processes responsible for rapid flux increases, so-called flares. We report on an extremely bright X-ray flare in the high-peaked BL Lacertae object Mrk 421 that occurred simultaneously with enhanced $γ$-ray activity detected at very high energies (VHE) by FACT on 2019 June 9. We triggered an observation with XMM-Newton, which observed the source quasi-continuously for 25 hours. We find that the source was in the brightest state ever observed using XMM-Newton, reaching a flux of $2.8\times10^{-9}$ erg cm$^{-2}$ s$^{-1}$ over an energy range of 0.3 - 10 keV. We perform a spectral and timing analysis to reveal the mechanisms of particle acceleration and to search for the shortest source-intrinsic time scales. Mrk 421 exhibits the typical harder-when-brighter behaviour throughout the observation and shows a clock-wise hysteresis pattern, which indicates that the cooling dominates over the acceleration process. While the X-ray emission in different sub-bands is highly correlated, we can exclude large time lags as the computed zDCFs are consistent with a zero lag. We find rapid variability on time scales of 1 ks for the 0.3 - 10 keV band and down to 300s in the hard X-ray band (4 - 10 keV). Taking these time scales into account, we discuss different models to explain the observed X-ray flare, and find that a plasmoid-dominated magnetic reconnection process is able to describe our observation best. △ Less

Submitted 29 February, 2024; originally announced February 2024.

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

A catalogue of dual-field interferometric binary calibrators

Authors: M. Nowak, S. Lacour, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, B. Charnay, G. Chauvin, A. Chavez, E. Choquet, V. Christiaens, Y. Clénet, V. Coudé du Foresto, A. Cridland , et al. (75 additional authors not shown)

Abstract: Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to with… ▽ More Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages, 6 figures. Submitted to A&A

The orbit of HD 142527 B is too compact to explain many of the disc features

Authors: M. Nowak, S. Rowther, S. Lacour, F. Meru, R. Nealon, D. J. Price

Abstract: HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13… ▽ More HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13 $\pm$ 0.03 $M_\odot{}$), which is thought to be responsible for most of the structures observed in the surrounding disc. We gathered VLTI/GRAVITY observations of HD 142527, either from our own programmes or from the ESO archive. We used this inhomogeneous set of data to extract a total of seven high-precision measurements of the relative astrometry between HD 142527 A and B, spread from mid-2017 to early 2021. Combined with what is available in the literature, we now have 9 yr of astrometric monitoring on HD 142527. We used orbit fitting tools to determine the orbital parameters of HD 142527 B, and used them as inputs for a 3D hydrodynamical model of the disc to determine whether or not the binary is able to create the structures observed in the disc. Our VLTI/GRAVITY astrometry gives excellent constraints on the orbit of HD 142527 B. We show that the secondary is following an orbit of semi-major axis a = 10.80 $\pm$ 0.22 au, with moderate eccentricity (e = 0.47 $\pm$ 0.01). With such a compact orbit, we show that HD 142527 B can only generate a gap and spiral arms of $\sim$30 au in the disc, which is much smaller than what is revealed by observations. Even from a theoretical standpoint, the observed cavity size of $\sim$100 au far exceeds even the most generous predictions for a companion like HD 142527 B on such a compact orbit. Thus, we conclude that the low-mass companion cannot be solely responsible for the observed morphology of the disc surrounding the system. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 11 pages, 6 figures. Accepted for publication in A&A

Upgrading the GRAVITY fringe tracker for GRAVITY+: Tracking the white light fringe in the non-observable Optical Path Length state-space

Authors: M. Nowak, S. Lacour, R. Abuter, J. Woillez, R. Dembet, M. S. Bordoni, G. Bourdarot, B. Courtney-Barrer, D. Defrère, A. Drescher, F. Eisenhauer, M. Fabricius, H. Feuchtgruber, R. Frahm, P. Garcia, S. Gillessen, V. Gopinath, J. Graf, S. Hoenig, L. Kreidberg, R. Laugier, J. B. Le Bouquin, D. Lutz, F. Mang, F. Millour , et al. (13 additional authors not shown)

Abstract: Aims. As part of the ongoing GRAVITY+ upgrade of the Very Large Telescope Interferometer infrastructure, we aim to improve the performance of the GRAVITY Fringe-Tracker, and to enable its use by other instruments. Methods. We modify the group delay controller to consistently maintain tracking in the white light fringe, characterised by a minimum group delay. Additionally, we introduce a novel appr… ▽ More Aims. As part of the ongoing GRAVITY+ upgrade of the Very Large Telescope Interferometer infrastructure, we aim to improve the performance of the GRAVITY Fringe-Tracker, and to enable its use by other instruments. Methods. We modify the group delay controller to consistently maintain tracking in the white light fringe, characterised by a minimum group delay. Additionally, we introduce a novel approach in which fringe-tracking is performed in the non-observable Optical Path Length state-space, using a covariance-weighted Kalman filter and an auto-regressive model of the disturbance. We outline this new state-space representation, and the formalism we use to propagate the state-vector and generate the control signal. While our approach is presented specifically in the context of GRAVITY/GRAVITY+, it can easily be adapted to other instruments or interferometric facilities. Results. We successfully demonstrate phase delay tracking within a single fringe, with any spurious phase jumps detected and corrected in less than 100 ms. We also report a significant performance improvement, as evidenced by a reduction of about 30 to 40% in phase residuals, and a much better behaviour under sub-optimal atmospheric conditions. Compared to what was observed in 2019, the median residuals have decreased from 150 nm to 100 nm on the Auxiliary Telescopes and from 250 nm to 150 nm on the Unit Telescopes. Conclusions. The improved phase-delay tracking combined with whit light fringe tracking means that from now-on, the GRAVITY Fringe-Tracker can be used by other instruments operating in different wavebands. The only limitation remains the need for an optical path dispersion adjustment. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 16 pages, 8 figures. Accepted for publication in A&A

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

A dynamical measure of the black hole mass in a quasar 11 billion years ago

Authors: R. Abuter, F. Allouche, A. Amorim, C. Bailet, A. Berdeu, J. -P. Berger, P. Berio, A. Bigioli, O. Boebion, M. -L. Bolzer, H. Bonnet, G. Bourdarot, P. Bourget, W. Brandner, Y. Cao, R. Conzelmann, M. Comin, Y. Clénet, B. Courtney-Barrer, R. Davies, D. Defrère, A. Delboulbé, F. Delplancke-Ströbele, R. Dembet, J. Dexter , et al. (102 additional authors not shown)

Abstract: Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves… ▽ More Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves with time; a key epoch to probe this relationship is at the peaks of star formation and black hole growth 8-12 billion years ago (redshifts 1-3). Here we report a dynamical measurement of the mass of the black hole in a luminous quasar at a redshift of 2, with a look back time of 11 billion years, by spatially resolving the broad line region. We detect a 40 micro-arcsecond (0.31 pc) spatial offset between the red and blue photocenters of the H$α$ line that traces the velocity gradient of a rotating broad line region. The flux and differential phase spectra are well reproduced by a thick, moderately inclined disk of gas clouds within the sphere of influence of a central black hole with a mass of 3.2x10$^{8}$ solar masses. Molecular gas data reveal a dynamical mass for the host galaxy of 6x10$^{11}$ solar masses, which indicates an under-massive black hole accreting at a super-Eddington rate. This suggests a host galaxy that grew faster than the supermassive black hole, indicating a delay between galaxy and black hole formation for some systems. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 5 pages Main text, 8 figures, 2 tables, to be published in Nature, under embargo until 29 January 2024 16:00 (London)

VLTI/GRAVITY Provides Evidence the Young, Substellar Companion HD 136164 Ab formed like a "Failed Star"

Authors: William O. Balmer, L. Pueyo, S. Lacour, J. J. Wang, T. Stolker, J. Kammerer, N. Pourré, M. Nowak, E. Rickman, S. Blunt, A. Sivaramakrishnan, D. Sing, K. Wagner, G. -D. Marleau, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, M. S. Bordoni , et al. (71 additional authors not shown)

Abstract: Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or c… ▽ More Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or core fragmentation)? Additionally, their orbits can yield model-independent dynamical masses, and when paired with wide wavelength coverage and accurate system age estimates, can constrain evolutionary models in a regime where the models have a wide dispersion depending on initial conditions. We present new interferometric observations of the $16\,\mathrm{Myr}$ substellar companion HD~136164~Ab (HIP~75056~Ab) with VLTI/GRAVITY and an updated orbit fit including proper motion measurements from the Hipparcos-Gaia Catalogue of Accelerations. We estimate a dynamical mass of $35\pm10\,\mathrm{M_J}$ ($q\sim0.02$), making HD~136164~Ab the youngest substellar companion with a dynamical mass estimate. The new mass and newly constrained orbital eccentricity ($e=0.44\pm0.03$) and separation ($22.5\pm1\,\mathrm{au}$) could indicate that the companion formed via the low-mass tail of the Initial Mass Function. Our atmospheric fit to the \texttt{SPHINX} M-dwarf model grid suggests a sub-solar C/O ratio of $0.45$, and $3\times$ solar metallicity, which could indicate formation in the circumstellar disk via disk fragmentation. Either way, the revised mass estimate likely excludes ``bottom-up" formation via core accretion in the circumstellar disk. HD~136164~Ab joins a select group of young substellar objects with dynamical mass estimates; epoch astrometry from future \textit{Gaia} data releases will constrain the dynamical mass of this crucial object further. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astronomical Journal. 9 figures, 3 tables

The High Energy X-ray Probe (HEX-P): Probing Accretion onto Stellar Mass Black Holes

Authors: Riley Connors, John Tomsick, Paul Draghis, Benjamin Coughenour, Aarran Shaw, Javier Garcia, Dominic Walton, Kristin Madsen, Daniel Stern, Nicole Cavero Rodriguez, Thomas Dauser, Melania Del Santo, Jiachen Jiang, Henric Krawczynski, Honghui Liu, Joseph Neilsen, Michael Nowak, Sean Pike, Andrea Santangelo, Navin Sridhar, Andrew West, Joern Wilms, the HEX-P Team

Abstract: Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity -- and most importantly, to measure the angular momentum, or spin, of black holes, and its r… ▽ More Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity -- and most importantly, to measure the angular momentum, or spin, of black holes, and its role as a powering mechanism for relativistic astrophysical phenomena. Comprehensive characterization of the disk-corona system of accreting black holes, and their co-evolution, is fundamental to measurements of black hole spin. Here, we use simulated data to demonstrate how key unanswered questions in the study of accreting stellar-mass black holes will be addressed by the {\it High Energy X-ray Probe} (\hexp). \hexp\ is a probe-class mission concept that will combine high spatial resolution X-ray imaging and broad spectral coverage ($0.2\mbox{--}80$keV) with a sensitivity superior to current facilities (including \xmm\ and \nustar) to enable revolutionary new insights into a variety of important astrophysical problems. We illustrate the capability of \hexp\ to: 1) measure the evolving structures of black hole binary accretion flows down to low ($\lesssim0.1\%$) Eddington-scaled luminosities via detailed X-ray reflection spectroscopy; 2) provide unprecedented spectral observations of the coronal plasma, probing its elusive geometry and energetics; 3) perform detailed broadband studies of stellar mass black holes in nearby galaxies, thus expanding the repertoire of sources we can use to study accretion physics and determine the fundamental nature of black holes; and 4) act as a complementary observatory to a range of future ground and space-based astronomical observatories, thus providing key spectral measurements of the multi-component emission from the inner accretion flows of BH-XRBs. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 7 figures

Polarization analysis of the VLTI and GRAVITY

Authors: GRAVITY Collaboration, F. Widmann, X. Haubois N. Schuhler, O. Pfuhl, F. Eisenhauer, S. Gillessen, N. Aimar, A. Amorim, M. Bauböck, J. B. Berger, H. Bonnet, G. Bourdarot, W. Brandner, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, A. Drescher, A. Eckart, H. Feuchtgruber, N. M. Förster Schreiber, P. Garcia, E. Gendron, R. Genzel, M. Hartl , et al. (37 additional authors not shown)

Abstract: The goal of this work is to characterize the polarization effects of the VLTI and GRAVITY. This is needed to calibrate polarimetric observations with GRAVITY for instrumental effects and to understand the systematic error introduced to the astrometry due to birefringence when observing targets with a significant intrinsic polarization. By combining a model of the VLTI light path and its mirrors an… ▽ More The goal of this work is to characterize the polarization effects of the VLTI and GRAVITY. This is needed to calibrate polarimetric observations with GRAVITY for instrumental effects and to understand the systematic error introduced to the astrometry due to birefringence when observing targets with a significant intrinsic polarization. By combining a model of the VLTI light path and its mirrors and dedicated experimental data, we construct a full polarization model of the VLTI UTs and the GRAVITY instrument. We first characterize all telescopes together to construct a UT calibration model for polarized targets. We then expand the model to include the differential birefringence. With this, we can constrain the systematic errors for highly polarized targets. Together with this paper, we publish a standalone Python package to calibrate the instrumental effects on polarimetric observations. This enables the community to use GRAVITY to observe targets in a polarimetric observing mode. We demonstrate the calibration model with the galactic center star IRS 16C. For this source, we can constrain the polarization degree to within 0.4 % and the polarization angle within 5 deg while being consistent with the literature. Furthermore, we show that there is no significant contrast loss, even if the science and fringe-tracker targets have significantly different polarization, and we determine that the phase error in such an observation is smaller than 1 deg, corresponding to an astrometric error of 10 μas. With this work, we enable the use of the polarimetric mode with GRAVITY/UTs and outline the steps necessary to observe and calibrate polarized targets. We demonstrate that it is possible to measure the intrinsic polarization of astrophysical sources with high precision and that polarization effects do not limit astrometric observations of polarized targets. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: Accepted by A&A

IXPE observation confirms a high spin in the accreting black hole 4U 1957+115

Authors: L. Marra, M. Brigitte, N. Rodriguez Cavero, S. Chun, J. F. Steiner, M. Dovčiak, M. Nowak, S. Bianchi, F. Capitanio, A. Ingram, G. Matt, F. Muleri, J. Podgorný, J. Poutanen, J. Svoboda, R. Taverna, F. Ursini, A. Veledina, A. De Rosa, J. A. Garcia, A. A. Lutovinov, I. A. Mereminskiy, R. Farinelli, S. Gunji, P. Kaaret , et al. (91 additional authors not shown)

Abstract: We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer in May 2023. The binary system has been in a high-soft spectral state since its discovery and is thought to host a black hole. The $\sim$571 ks observation reveals a linear polarisation degree of $1.9\% \pm 0.6\%$ and a polarisation angl… ▽ More We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer in May 2023. The binary system has been in a high-soft spectral state since its discovery and is thought to host a black hole. The $\sim$571 ks observation reveals a linear polarisation degree of $1.9\% \pm 0.6\%$ and a polarisation angle of $-41^\circ.8 \pm 7^\circ.9$ in the 2-8 keV energy range. Spectral modelling is consistent with the dominant contribution coming from the standard accretion disc, while polarimetric data suggest a significant role of returning radiation: photons that are bent by strong gravity effects and forced to return to the disc surface, where they can be reflected before eventually reaching the observer. In this setting, we find that models with a black hole spin lower than 0.96 and an inclination lower than $50^\circ$ are disfavoured. △ Less

Submitted 8 February, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 12 pages, 10 figures, 2 tables, accepted for publication in A&A

First VLTI/GRAVITY Observations of HIP 65426 b: Evidence for a Low or Moderate Orbital Eccentricity

Authors: S. Blunt, W. O. Balmer, J. J. Wang, S. Lacour, S. Petrus, G. Bourdarot, J. Kammerer, N. Pourré, E. Rickman, J. Shangguan, T. Winterhalder, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (73 additional authors not shown)

Abstract: Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging… ▽ More Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior-dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best fit ExoREM model with solar metallicity and C/O=0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available at https://github.com/sblunt/hip65426. △ Less

Submitted 6 October, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 28 pages, 14 figures. Revised and resubmitted to AJ

VLTI/GRAVITY Observations and Characterization of the Brown Dwarf Companion HD 72946 B

Authors: W. O. Balmer, L. Pueyo, T. Stolker, H. Reggiani, S. Lacour, A. -L. Maire, P. Mollière, M. Nowak, D. Sing, N. Pourré, S. Blunt, J. J. Wang, E. Rickman, Th. Henning, K. Ward-Duong, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet , et al. (74 additional authors not shown)

Abstract: Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. W… ▽ More Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint $\mathrm{M_B}=69.5\pm0.5~\mathrm{M_{Jup}}$ assuming a strong prior on the host star mass $\mathrm{M_A}=0.97\pm0.01~\mathrm{M_\odot}$ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution on effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure-temperature or cloud parameters within our current retrieval setup. Independent of cloud prescription and prior assumptions, we find agreement within $\pm2\,σ$ between the C/O ratio of the host ($0.52\pm0.05)$ and brown dwarf ($0.43$ to $0.63$), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high ($0.6-0.8$) in light of an excellent agreement of the data with the solar abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface-gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis. △ Less

Submitted 15 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astrophysical Journal. 32 pages, 14 figures, 7 tables. v2&3 correct errors in co-author's affiliations, figure rendering, and some grant acknowledgements

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

NICER/NuSTAR Characterization of 4U 1957+11: A Near Maximally Spinning Black Hole Potentially in the Mass Gap

Authors: Erin Barillier, Victoria Grinberg, David Horn, Michael A. Nowak, Ronald A. Remillard, James F. Steiner, Dominic J. Walton, Jörn Wilms

Abstract: 4U 1957+11 is a black hole candidate system that has been in a soft X-ray spectral state since its discovery. We present analyses of recent joint NICER and NuSTAR spectra, which are extremely well-described by a highly inclined disk accreting into a near maximally spinning black hole. Owing to the broad X-ray coverage of NuSTAR the fitted spin and inclination are strongly constrained for our hypot… ▽ More 4U 1957+11 is a black hole candidate system that has been in a soft X-ray spectral state since its discovery. We present analyses of recent joint NICER and NuSTAR spectra, which are extremely well-described by a highly inclined disk accreting into a near maximally spinning black hole. Owing to the broad X-ray coverage of NuSTAR the fitted spin and inclination are strongly constrained for our hypothesized disk models. The faintest spectra are observed out to 20 keV, even though their hard tail components are almost absent when described with a simple corona. The hard tail increases with luminosity, but shows clear two track behavior with one track having appreciably stronger tails. The disk spectrum color-correction factor is anti-correlated with the strength of the hard tail (e.g., as measured by the Compton $y$ parameter). Although the spin and inclination parameters are strongly constrained for our chosen model, the mass and distance are degenerate parameters. We use our spectral fits, along with a theoretical prior on color-correction, an observational prior on likely fractional Eddington luminosity, and an observational prior on distance obtained from Gaia studies, to present mass and distance contours for this system. The most likely parameters, given our presumed disk model, suggest a 4.6 $\mathrm{M_\odot}$ black hole at 7.8 kpc observed at luminosities ranging from $\approx 1.7\%$--$9\%$ of Eddington. This would place 4U 1957+11 as one of the few actively accreting sources within the `mass gap' of ${\approx} 2$--$5\,\mathrm{M_\odot}$ where there are few known massive neutron stars or low mass black holes. Higher mass and distance, however, remain viable. △ Less

Submitted 22 January, 2023; originally announced January 2023.

Comments: 21 pages. To be published in the Astrophysical Journal

The GRAVITY+ Project: Towards All-sky, Faint-Science, High-Contrast Near-Infrared Interferometry at the VLTI

Authors: GRAVITY+ Collaboration, :, Roberto Abuter, Patricio Alarcon, Fatme Allouche, Antonio Amorim, Christophe Bailet, Helen Bedigan, Anthony Berdeu, Jean-Philippe Berger, Philippe Berio, Azzurra Bigioli, Richard Blaho, Olivier Boebion, Marie-Lena Bolzer, Henri Bonnet, Guillaume Bourdarot, Pierre Bourget, Wolfgang Brandner, Cesar Cardenas, Ralf Conzelmann, Mauro Comin, Yann Clénet, Benjamin Courtney-Barrer, Yigit Dallilar , et al. (112 additional authors not shown)

Abstract: The GRAVITY instrument has been revolutionary for near-infrared interferometry by pushing sensitivity and precision to previously unknown limits. With the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) in GRAVITY+, these limits will be pushed even further, with vastly improved sky coverage, as well as faint-science and high-contrast capabilities. This upgrade includes the im… ▽ More The GRAVITY instrument has been revolutionary for near-infrared interferometry by pushing sensitivity and precision to previously unknown limits. With the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) in GRAVITY+, these limits will be pushed even further, with vastly improved sky coverage, as well as faint-science and high-contrast capabilities. This upgrade includes the implementation of wide-field off-axis fringe-tracking, new adaptive optics systems on all Unit Telescopes, and laser guide stars in an upgraded facility. GRAVITY+ will open up the sky to the measurement of black hole masses across cosmic time in hundreds of active galactic nuclei, use the faint stars in the Galactic centre to probe General Relativity, and enable the characterisation of dozens of young exoplanets to study their formation, bearing the promise of another scientific revolution to come at the VLTI. △ Less

Submitted 19 January, 2023; originally announced January 2023.

Comments: Published in the ESO Messenger

The spectral-timing analysis of Cygnus X-1 with Insight-HXMT

Authors: M. Zhou, V. Grinberg, Q. -C. Bu, A. Santangelo, F. Cangemi, C. M. Diez, O. König, L. Ji, M. A. Nowak, K. Pottschmidt, J. Rodriguez, J. Wilms, S. Zhang, J. -L. Qu, S. -N. Zhang

Abstract: Cygnus X-1, as the first discovered black hole binary, is a key source for understanding the mechanisms of state transitions, and the scenarios of accretion in extreme gravity fields. We present a spectral-timing analysis of observations taken with the Insight-HXMT mission, focusing on the spectral-state dependent timing properties in the broad energy range of 1--150 keV, thus extending previous R… ▽ More Cygnus X-1, as the first discovered black hole binary, is a key source for understanding the mechanisms of state transitions, and the scenarios of accretion in extreme gravity fields. We present a spectral-timing analysis of observations taken with the Insight-HXMT mission, focusing on the spectral-state dependent timing properties in the broad energy range of 1--150 keV, thus extending previous RXTE-based studies to both lower and higher energies. Our main results are the following: a) We successfully use a simple empirical model to fit all spectra, confirming that the reflection component is stronger in the soft state than in the hard state; b) The evolution of the total fractional root mean square (rms) depends on the selected energy band and the spectral shape, which is a direct result of the evolution of the power spectral densities (PSDs); c) In the hard/intermediate state, we see clear short-term variability features and a positive correlation between central frequencies of the variability components and the soft photon index $Γ_1$, also at energies above 15 keV. The power spectrum is dominated by red noise in the soft state instead. These behaviors can be traced to at least 90 keV; d) The coherence and the phase-lag spectra show different behaviors dependent on different spectral shapes. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 18 pages, 16 figures, 1 table. accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 666, A172 (2022)

Direct discovery of the inner exoplanet in the HD206893 system. Evidence for deuterium burning in a planetary-mass companion

Authors: S. Hinkley, S. Lacour, G. -D. Marleau, A. M. Lagrange, J. J. Wang, J. Kammerer, A. Cumming, M. Nowak, L. Rodet, T. Stolker, W. -O. Balmer, S. Ray, M. Bonnefoy, P. Mollière, C. Lazzoni, G. Kennedy, C. Mordasini, R. Abuter, S. Aigrain, A. Amorim, R. Asensio-Torres, C. Babusiaux, M. Benisty, J. -P. Berger, H. Beust , et al. (89 additional authors not shown)

Abstract: Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoi… ▽ More Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $μ$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au. △ Less

Submitted 3 April, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted to A&A

Multiwavelength astrophysics of the blazar OJ 287 and the project MOMO

Authors: S. Komossa, A. Kraus, D. Grupe, M. L. Parker, A. Gonzalez, L. C. Gallo, M. A. Gurwell, S. Laine, S. Yao, S. Chandra, L. Dey, J. L. Gomez, A. Gopakumar. K. Hada, D. Haggard, A. R. Hollett, H. Jermak, S. Jorstad, T. P. Krichbaum, S. Markoff, C. McCall, J. Neilsen, M. Nowak

Abstract: We are carrying out the densest and longest multiyear, multiwavelength monitoring project of OJ 287 ever done. The project MOMO (Multiwavelength Observations and Modelling of OJ 287) covers wavelengths from the radio to the high-energy regime. A few selected observations are simultaneous with those of the Event Horizon Telescope (EHT). MOMO aims at understanding disk-jet physics and at testing pre… ▽ More We are carrying out the densest and longest multiyear, multiwavelength monitoring project of OJ 287 ever done. The project MOMO (Multiwavelength Observations and Modelling of OJ 287) covers wavelengths from the radio to the high-energy regime. A few selected observations are simultaneous with those of the Event Horizon Telescope (EHT). MOMO aims at understanding disk-jet physics and at testing predictions of the binary black hole scenario of OJ 287. Here, we present a discussion of extreme outburst and minima states in context, and then focus on the recent flux and spectral evolution between 2021 and May 2022, including an ongoing bright radio flare. Further, we show that there is no evidence for precursor flare activity in our optical-UV-X-ray light curves that would be associated with any secondary supermassive black hole (SMBH) disk impact and that was predicted to start as thermal flare on 2021 December 23. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: Submitted to AN. Contribution to the XMM-Newton Workshop 2022 (Madrid, Spain). 5 figures

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

Authors: Maciek Wielgus, Nicola Marchili, Ivan Marti-Vidal, Garrett K. Keating, Venkatessh Ramakrishnan, Paul Tiede, Ed Fomalont, Sara Issaoun, Joey Neilsen, Michael A. Nowak, Lindy Blackburn, Charles F. Gammie, Ciriaco Goddi, Daryl Haggard, Daeyoung Lee, Monika Moscibrodzka, Alexandra J. Tetarenko, Geoffrey C. Bower, Chi-Kwan Chan, Koushik Chatterjee, Paul M. Chesler, Jason Dexter, Sheperd S. Doeleman, Boris Georgiev, Mark Gurwell , et al. (6 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data s… ▽ More The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Journal ref: The Astrophysical Journal Letters, 930:L19 (2022)

Efficiently combining Alpha CenA multi-epoch high-contrast imaging data. Application of K-Stacker to the 80 hrs NEAR campaign

Authors: Hervé Le Coroller, Mathias Nowak, Kevin Wagner, Markus Kasper, Gael Chauvin, Celia Desgrange, Simon Conseil

Abstract: Keplerian-Stacker is an algorithm able to combine multiple observations acquired at different epochs taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. In 2019, a total of 100 hours of observation were allocated to VLT VISIR-NEAR, a collaboration between ESO and Breakthrough Initiatives, to search for low mass planets in the ha… ▽ More Keplerian-Stacker is an algorithm able to combine multiple observations acquired at different epochs taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. In 2019, a total of 100 hours of observation were allocated to VLT VISIR-NEAR, a collaboration between ESO and Breakthrough Initiatives, to search for low mass planets in the habitable zone of the Alpha Cen AB binary system. A weak signal (S/N = 3) was reported around Alpha Cen A, at a separation of 1.1 a.u. which corresponds to the habitable zone. We have re-analysed the NEAR data using K-Stacker. This algorithm is a brute-force method able to find planets in time series of observations and to constrain their orbital parameters, even if they remain undetected in a single epoch. We scanned a total of about 3.5e+5 independent orbits, among which about 15 % correspond to fast moving orbits on which planets cannot be detected without taking into account the orbital motion. We find only a single planet candidate, which matches the C1 detection reported in Wagner et al. 2021. Despite the significant amount of time spent on this target, the orbit of this candidate remains poorly constrained due to these observations being closely distributed in 34 days. We argue that future single-target deep surveys would benefit from a K-Stacker based strategy, where the observations would be split over a significant part of the expected orbital period to better constrain the orbital parameters. This application of K-Stacker on high contrast imaging data in the mid-infrared demonstrates the capability of this algorithm to aid in the search for Earth-like planets in the habitable zone of the nearest stars with future instruments of the E-ELT such as METIS. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 9 pages, 11 figures, K-Stacker github link

Journal ref: A&A 667, A142 (2022)

First Light for GRAVITY Wide: Large Separation Fringe Tracking for the Very Large Telescope Interferometer

Authors: GRAVITY+ Collaboration, :, R. Abuter, F. Allouche, A. Amorim, C. Bailet, M. Bauböck, J. -P. Berger, P. Berio, A. Bigioli, O. Boebion, M. L. Bolzer, H. Bonnet, G. Bourdarot, P. Bourget, W. Brandner, Y. Clénet, B. Courtney-Barrer, Y. Dallilar, R. Davies, D. Defrère, A. Delboulbé, F. Delplancke, R. Dembet, P. T. de Zeeuw , et al. (92 additional authors not shown)

Abstract: GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines t… ▽ More GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines the dual-beam capabilities of the VLTI and the GRAVITY instrument to increase the maximum separation between the science target and the reference star from 2 arcseconds with the 8 m UTs up to several 10 arcseconds, limited only by the Earth's turbulent atmosphere. This increases the sky-coverage of GRAVITY by two orders of magnitude, opening up milliarcsecond resolution observations of faint objects, and in particular the extragalactic sky. The first observations in 2019 - 2022 include first infrared interferometry of two redshift $z\sim2$ quasars, interferometric imaging on the binary system HD 105913A, and repeated observations of multiple star systems in the Orion Trapezium Cluster. We find the coherence loss between the science object and fringe-tracking reference star well described by the turbulence of the Earth's atmosphere. We confirm that the larger apertures of the UTs result in higher visibilities for a given separation due to larger overlap of the projected pupils on sky and give predictions for visibility loss as a function of separation to be used for future planning. △ Less

Submitted 23 August, 2022; v1 submitted 1 June, 2022; originally announced June 2022.

Comments: 15 pages, 12 figures, 5 tables. Accepted by A&A

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

In-depth direct imaging and spectroscopic characterization of the young Solar System analog HD 95086

Authors: C. Desgrange, G. Chauvin, V. Christiaens, F. Cantalloube, L. -X. Lefranc, H. Le Coroller, P. Rubini, G. P. P. L. Otten, H. Beust, M. Bonavita, P. Delorme, M. Devinat, R. Gratton, A. -M. Lagrange, M. Langlois, D. Mesa, J. Milli, J. Szulágyi, M. Nowak, L. Rodet, P. Rojo, S. Petrus, M. Janson, T. Henning, Q. Kral , et al. (26 additional authors not shown)

Abstract: Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmosphe… ▽ More Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmospheric properties of HD 95086 b, and to search for additional planets in this system. Methods. SPHERE observations, spread over ten epochs from 2015 to 2019 and including five new datasets, were used. Combined with archival observations, from VLT/NaCo (2012-2013) and Gemini/GPI (2013-2016), the extended set of astrometric measurements allowed us to refine the orbital properties of HD 95086 b. We also investigated the spectral properties and the presence of a circumplanetary disk around HD 95086 b by using the special fitting tool exploring the diversity of several atmospheric models. In addition, we improved our detection limits in order to search for a putative planet c via the K-Stacker algorithm. Results. We extracted for the first time the JH low-resolution spectrum of HD 95086 b by stacking the six best epochs, and confirm its very red spectral energy distribution. Combined with additional datasets from GPI and NaCo, our analysis indicates that this very red color can be explained by the presence of a circumplanetary disk around planet b, with a range of high-temperature solutions (1400-1600 K) and significant extinction (Av > 10 mag), or by a super-solar metallicity atmosphere with lower temperatures (800-1300 K), and small to medium amount of extinction (Av < 10 mag). We do not find any robust candidates for planet c, but give updated constraints on its potential mass and location. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 29 pages, 20 figures, A&A, accepted

Journal ref: A&A 664, A139 (2022)

Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions

Authors: Paul Mollière, Tamara Molyarova, Bertram Bitsch, Thomas Henning, Aaron Schneider, Laura Kreidberg, Christian Eistrup, Remo Burn, Evert Nasedkin, Dmitry Semenov, Christoph Mordasini, Martin Schlecker, Kamber R. Schwarz, Sylvestre Lacour, Mathias Nowak, Matthäus Schulik

Abstract: Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric comp… ▽ More Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ, base formation inversion code is available at: https://gitlab.com/mauricemolli/formation-inversion

Multi-wavelength Variability of Sagittarius A* in July 2019

Authors: H. Boyce, D. Haggard, G. Witzel, S. von Fellenberg, S. P. Willner, E. E. Becklin, T. Do, A. Eckart, G. G. Fazio, M. A. Gurwell, J. L. Hora, S. Markoff, M. R. Morris, J. Neilsen, M. Nowak, H. A. Smith, S. Zhang

Abstract: We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Tw… ▽ More We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Two dates show moderate variability with no significant lags between the submm and the infrared at 99% confidence. July 18 captured a moderately bright NIR flare (F_K ~ 15 mJy) simultaneous with an X-ray flare (F ~ 0.1 cts/s) that most likely preceded bright submm flux (F ~ 5.5 Jy) by about +34 (+14 -33) minutes at 99% confidence. The uncertainty in this lag is dominated by the fact that we did not observe the peak of the submm emission. A synchrotron source cooled through adiabatic expansion can describe a rise in the submm once the synchrotron-self-Compton NIR and X-ray peaks have faded. This model predicts high GHz and THz fluxes at the time of the NIR/X-ray peak and electron densities well above those implied from average accretion rates for Sgr A*. However, the higher electron density postulated in this scenario would be in agreement with the idea that 2019 was an extraordinary epoch with a heightened accretion rate. Since the NIR and X-ray peaks can also be fit by a non-thermal synchrotron source with lower electron densities, we cannot rule out an unrelated chance coincidence of this bright submm flare with the NIR/X-ray emission. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: Accepted for publication in The Astrophysical Journal

Deep Images of the Galactic Center with GRAVITY

Authors: GRAVITY Collaboration, R. Abuter, N. Aimar, A. Amorim, P. Arras, M. Bauböck, J. P. Berger, H. Bonnet, W. Brandner, G. Bourdarot, V. Cardoso, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, Y. Dallilar, A. Drescher, F. Eisenhauer, T. Enßlin, N. M. Förster Schreiber, P. Garcia, F. Gao, E. Gendron, R. Genzel, S. Gillessen , et al. (43 additional authors not shown)

Abstract: Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole s… ▽ More Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole spin. Furthermore, it provides the starting point for astrometric fitting to derive highly accurate stellar positions. Here, we present $\mathrm{G^R}$, a new imaging tool specifically designed for Galactic Center observations with GRAVITY. The algorithm is based on a Bayesian interpretation of the imaging problem, formulated in the framework of information field theory and building upon existing works in radio-interferometric imaging. Its application to GRAVITY observations from 2021 yields the deepest images to date of the Galactic Center on scales of a few milliarcseconds. The images reveal the complicated source structure within the central $100\,\mathrm{mas}$ around Sgr A*, where we detected the stars S29 and S55 and confirm S62 on its trajectory, slowly approaching Sgr A*. Furthermore, we were able to detect S38, S42, S60, and S63 in a series of exposures for which we offset the fiber from Sgr A*. We provide an update on the orbits of all aforementioned stars. In addition to these known sources, the images also reveal a faint star moving to the west at a high angular velocity. We cannot find any coincidence with any known source and, thus, we refer to the new star as S300. From the flux ratio with S29, we estimate its K-band magnitude as $m_\mathrm{K}\left(\mathrm{S300}\right)\simeq 19.0 - 19.3$. Images obtained with CLEAN confirm the detection. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 24 pages, 14 figures

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

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

The mass of Beta Pictoris c from Beta Pictoris b orbital motion

Authors: S. Lacour, J. J. Wang, L. Rodet, M. Nowak, J. Shangguan, H. Beust, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, S. Blunt, A. Boccaletti, A. Bohn, M. -L. Bolzer, M. Bonnefoy, H. Bonnet, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay, G. Chauvin, E. Choquet , et al. (74 additional authors not shown)

Abstract: We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was… ▽ More We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was also reobserved at a separation of 96\,mas, confirming the previous orbital estimations. From the astrometry of planet b only, we can (i) detect the presence of $β$ Pictoris c and (ii) constrain its mass to $10.04^{+4.53}_{-3.10}\,M_{\rm Jup}$. If one adds the astrometry of $β$ Pictoris c, the mass is narrowed down to $9.15^{+1.08}_{-1.06}\,M_{\rm Jup}$. The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to $8.89^{+0.75}_{-0.75}\,M_{\rm Jup}$. With a semimajor axis of $2.68\pm0.02$\,au, a period of $1221\pm15$ days, and an eccentricity of $0.32\pm0.02$, the orbital parameters of $β$ Pictoris c are now constrained as precisely as those of $β$ Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets' dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted in A&A Letter

Journal ref: A&A 654, L2 (2021)

Constraining particle acceleration in Sgr A* with simultaneous GRAVITY, Spitzer, NuSTAR and Chandra observations

Authors: R. Abuter, A. Amorim, M. Bauböck, F. Baganoff, J. P. Berge, H. Boyce, H. Bonnet, W. Brandner, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, Y. Dallilar, A. Drescher, A. Eckart, F. Eisenhauer, G. G. Fazio, N. M. Förster Schreiber, K. Foster, C. Gammie, P. Garcia, F. Gao, E. Gendron, R. Genzel, G. Ghisellini , et al. (59 additional authors not shown)

Abstract: We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $νL_ν\propto ν^{0.5\pm0.2}$; the spectral slope observed in the X-ray ban… ▽ More We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $νL_ν\propto ν^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $νL_ν\propto ν^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $γ_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $γ_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $γ_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $γ_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

Journal ref: A&A 654, A22 (2021)

GRAVITY K-band spectroscopy of HD 206893 B: brown dwarf or exoplanet

Authors: J. Kammerer, S. Lacour, T. Stolker, P. Mollière, D. K. Sing, E. Nasedkin, P. Kervella, J. J. Wang, K. Ward-Duong, M. Nowak, R. Abuter, A. Amorim, R. Asensio-Torres, M. Bauböck, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, A. Bohn, M. -L. Bolzer, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube , et al. (72 additional authors not shown)

Abstract: We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectra… ▽ More We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectral features, we include additional extinction by high-altitude dust clouds made of enstatite grains in the atmospheric model fits. We also infer the orbital parameters of HD 206893 B by combining the $\sim 100~μ\text{as}$ precision astrometry from GRAVITY with data from the literature and constrain the mass and position of HD 206893 C based on the Gaia proper motion anomaly of the system. The extremely red color and the very shallow $1.4~μ\text{m}$ water absorption feature of HD 206893 B can be fit well with the adapted atmospheric models and spectral retrievals. Altogether, our analysis suggests an age of $\sim 3$-$300~\text{Myr}$ and a mass of $\sim 5$-$30~\text{M}_\text{Jup}$ for HD 206893 B, which is consistent with previous estimates but extends the parameter space to younger and lower-mass objects. The GRAVITY astrometry points to an eccentric orbit ($e = 0.29^{+0.06}_{-0.11}$) with a mutual inclination of $< 34.4~\text{deg}$ with respect to the debris disk of the system. While HD 206893 B could in principle be a planetary-mass companion, this possibility hinges on the unknown influence of the inner companion on the mass estimate of $10^{+5}_{-4}~\text{M}_\text{Jup}$ from radial velocity and Gaia as well as a relatively small but significant Argus moving group membership probability of $\sim 61\%$. However, we find that if the mass of HD 206893 B is $< 30~\text{M}_\text{Jup}$, then the inner companion HD 206893 C should have a mass between $\sim 8$-$15~\text{M}_\text{Jup}$. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 28 pages, 19 figures, accepted for publication in A&A

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

MOLsphere and pulsations of the Galactic Center's red supergiant GCIRS 7 from VLTI/GRAVITY

Authors: GRAVITY Collaboration, G. Rodríguez-Coira, T. Paumard, G. Perrin, F. Vincent, R. Abuter, A. Amorim, M. Bauböck, J. P. Berger, H. Bonnet, W. Brandner, Y. Clénet, P. T. de Zeeuw, J. Dexter, A. Drescher, A. Eckart, F. Eisenhauer, N. M. Förster Schreiber, F. Gao, P. Garcia, E. Gendron, R. Genzel, S. Gillessen, M. Habibi, X. Haubois , et al. (33 additional authors not shown)

Abstract: GCIRS 7, the brightest star in the Galactic central parsec, formed $6\pm2$ Myr ago together with dozens of massive stars in a disk orbiting the central black-hole. It has been argued that GCIRS 7 is a pulsating body, on the basis of photometric variability. We present the first medium-resolution ($R=500$), K-band spectro-interferometric observations of GCIRS 7, using the GRAVITY instrument with th… ▽ More GCIRS 7, the brightest star in the Galactic central parsec, formed $6\pm2$ Myr ago together with dozens of massive stars in a disk orbiting the central black-hole. It has been argued that GCIRS 7 is a pulsating body, on the basis of photometric variability. We present the first medium-resolution ($R=500$), K-band spectro-interferometric observations of GCIRS 7, using the GRAVITY instrument with the four auxiliary telescopes of the ESO VLTI. We looked for variations using two epochs, namely 2017 and 2019. We find GCIRS 7 to be moderately resolved with a uniform-disk photospheric diameter of $θ^*_\text{UD}=1.55 \pm 0.03$ mas ($R^*_\text{UD}=1368 \pm 26$ $R_\odot$) in the K-band continuum. The narrow-band uniform-disk diameter increases above 2.3 $μ$m, with a clear correlation with the CO band heads in the spectrum. This correlation is aptly modeled by a hot ($T_\text{L}=2368\pm37$ K), geometrically thin molecular shell with a diameter of $θ_\text{L}=1.74\pm0.03$ mas, as measured in 2017. The shell diameter increased ($θ_\text{L}=1.89\pm0.03$ mas), while its temperature decreased ($T_\text{L}=2140\pm42$ K) in 2019. In contrast, the photospheric diameter $θ^*_\text{UD}$ and the extinction up to the photosphere of GCIRS 7 ($A_{\mathrm{K}_\mathrm{S}}=3.18 \pm 0.16$) have the same value within uncertainties at the two epochs. In the context of previous interferometric and photo-spectrometric measurements, the GRAVITY data allow for an interpretation in terms of photospheric pulsations. The photospheric diameter measured in 2017 and 2019 is significantly larger than previously reported using the PIONIER instrument ($θ_*=1.076 \pm 0.093$ mas in 2013 in the H band). The parameters of the photosphere and molecular shell of GCIRS 7 are comparable to those of other red supergiants that have previously been studied using interferometry. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 12 pages, 11 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics (A&A)

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

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

X-ray spectral components of the blazar and binary black hole candidate OJ 287 (2005-2020)

Authors: S. Komossa, D. Grupe, M. L. Parker, J. L. Gómez, M. J. Valtonen, M. A. Nowak, S. G. Jorstad, D. Haggard, S. Chandra, S. Ciprini, L. Dey, A. Gopakumar, K. Hada, S. Markoff, J. Neilsen

Abstract: We present a comprehensive analysis of all XMM-Newton spectra of OJ 287 spanning 15 years of X-ray spectroscopy of this bright blazar. We also report the latest results from our dedicated Swift UVOT and XRT monitoring of OJ 287 which started in 2015, along with all earlier public Swift data since 2005. During this time interval, OJ 287 was caught in extreme minima and outburst states. Its X-ray sp… ▽ More We present a comprehensive analysis of all XMM-Newton spectra of OJ 287 spanning 15 years of X-ray spectroscopy of this bright blazar. We also report the latest results from our dedicated Swift UVOT and XRT monitoring of OJ 287 which started in 2015, along with all earlier public Swift data since 2005. During this time interval, OJ 287 was caught in extreme minima and outburst states. Its X-ray spectrum is highly variable and encompasses all states seen in blazars from very flat to exceptionally steep. The spectrum can be decomposed into three spectral components: Inverse Compton (IC) emission dominant at low-states, super-soft synchrotron emission which becomes increasingly dominant as OJ 287 brightens, and an intermediately-soft (Gamma_x=2.2) additional component seen at outburst. This last component extends beyond 10 keV and plausibly represents either a second synchrotron/IC component and/or a temporary disk corona of the primary supermassive black hole (SMBH). Our 2018 XMM-Newton observation, quasi-simultaneous with the Event Horizon Telescope observation of OJ 287, is well described by a two-component model with a hard IC component of Gamma_x=1.5 and a soft synchrotron component. Low-state spectra limit any long-lived accretion disk/corona contribution in X-rays to a very low value of L_x/L_Edd < 5.6 times 10^(-4) (for M_(BH, primary) = 1.8 times 10^10 M_sun). Some implications for the binary SMBH model of OJ 287 are discussed. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 14 pages. MNRAS, in press

Supermassive binary black holes and the case of OJ 287

Authors: S. Komossa, S. Ciprini, L. Dey, L. C. Gallo, J. L. Gomez, A. Gonzalez, D. Grupe, A. Kraus, S. J. Laine, M. L. Parker, M. J. Valtonen, S. Chandra, A. Gopakumar, D. Haggard, M. A. Nowak

Abstract: Supermassive binary black holes (SMBBHs) are laboratories par excellence for relativistic effects, including precession effects in the Kerr metric and the emission of gravitational waves. Binaries form in the course of galaxy mergers, and are a key component in our understanding of galaxy evolution. Dedicated searches for SMBBHs in all stages of their evolution are therefore ongoing and many syste… ▽ More Supermassive binary black holes (SMBBHs) are laboratories par excellence for relativistic effects, including precession effects in the Kerr metric and the emission of gravitational waves. Binaries form in the course of galaxy mergers, and are a key component in our understanding of galaxy evolution. Dedicated searches for SMBBHs in all stages of their evolution are therefore ongoing and many systems have been discovered in recent years. Here we provide a review of the status of observations with a focus on the multiwavelength detection methods and the underlying physics. Finally, we highlight our ongoing, dedicated multiwavelength program MOMO (for Multiwavelength Observations and Modelling of OJ 287). OJ 287 is one of the best candidates to date for hosting a sub-parsec SMBBH. The MOMO program carries out a dense monitoring at >13 frequencies from radio to X-rays and especially with Swift since 2015. Results so far included: (1) The detection of two major UV-X-ray outbursts with Swift in 2016/17 and 2020; exhibiting softer-when-brighter behaviour. The non-thermal nature of the outbursts was clearly established and shown to be synchrotron radiation. (2) Swift multi-band dense coverage and XMM-Newton spectroscopy during EHT campaigns caught OJ 287 at an intermediate flux level with synchrotron and IC spectral components. (3) Discovery of a remarkable, giant soft X-ray excess with XMM and NuSTAR during the 2020 outburst. (4) Spectral evidence (at 2sigma) for a relativistically shifted iron absorption line in 2020. (5) The non-thermal 2020 outburst is consistent with an after-flare predicted by the SMBBH model of OJ 287. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Invited review at the COST Special Session on "The gravitational-wave Universe" during the XIX Serbian Astronomical Conference (Oct. 2020). 14 pages incl. 7 figures

Broadband Multi-wavelength Properties of M87 during the 2017 Event Horizon Telescope Campaign

Authors: J. C. Algaba, J. Anczarski, K. Asada, M. Balokovic, S. Chandra, Y. -Z. Cui, A. D. Falcone, M. Giroletti, C. Goddi, K. Hada, D. Haggard, S. Jorstad, A. Kaur, T. Kawashima, G. Keating, J. -Y. Kim, M. Kino, S. Komossa, E. V. Kravchenko, T. P. Krichbaum, S. -S. Lee, R. -S. Lu, M. Lucchini, S. Markoff, J. Neilsen , et al. (14 additional authors not shown)

Abstract: In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground,… ▽ More In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87's spectrum. We can exclude that the simultaneous gamma-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the gamma-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 43 pages, 19 figures, 10 tables; complete author list available in manuscript; The Astrophysical Journal Letters, 2021, 911, L11; publication doi: 3847/2041-8213/abef71, data doi: 10.25739/mhh2-cw46

High-precision astrometric studies in direct imaging with SPHERE

Authors: A. -L. Maire, G. Chauvin, A. Vigan, R. Gratton, M. Langlois, J. H. Girard, M. A. Kenworthy, J. -U. Pott, T. Henning, P. Kervella, S. Lacour, E. L. Rickman, A. Boccaletti, P. Delorme, M. R. Meyer, M. Nowak, S. P. Quanz, A. Zurlo

Abstract: Orbital monitoring of exoplanetary and stellar systems is fundamental for analysing their architecture, dynamical stability and evolution, and mechanisms of formation. Current high-contrast extreme-adaptive optics imagers like SPHERE, GPI, and SCExAO+CHARIS explore the population of giant exoplanets and brown dwarf and stellar companions beyond typically 10 au, covering generally a small fraction… ▽ More Orbital monitoring of exoplanetary and stellar systems is fundamental for analysing their architecture, dynamical stability and evolution, and mechanisms of formation. Current high-contrast extreme-adaptive optics imagers like SPHERE, GPI, and SCExAO+CHARIS explore the population of giant exoplanets and brown dwarf and stellar companions beyond typically 10 au, covering generally a small fraction of the orbit (<20%) leading to degeneracies and biases in the orbital parameters. Precise and robust measurements over time of the position of the companions are critical, which require good knowledge of the instrumental limitations and dedicated observing strategies. The homogeneous dedicated calibration strategy for astrometry implemented for SPHERE has facilitated high-precision studies by its users since its start of operation in 2014. As the precision of exoplanet imaging instruments is now reaching milliarcseconds and is expected to improve with the upcoming facilities, we initiated a community effort, triggered by the SPHERE experience, to share lessons learned for high-precision astrometry in direct imaging. A homogeneous strategy would strongly benefit the VLT community, in synergy with VLTI instruments like GRAVITY/GRAVITY+, future instruments like ERIS and MAVIS, and in preparation for the exploitation of the ELT's first instruments MICADO, HARMONI, and METIS. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Comments: In press for the ESO Messenger. 10 pages, 5 figures, 1 table

The $β$ Pictoris b Hill sphere transit campaign. Paper I: Photometric limits to dust and rings

Authors: M. A. Kenworthy, S. N. Mellon, J. I. Bailey III, R. Stuik, P. Dorval, G. -J. J. Talens, S. R. Crawford, E. E. Mamajek, I. Laginja, M. Ireland, B. Lomberg, R. B. Kuhn, I. Snellen, K. Zwintz, R. Kuschnig, G. M. Kennedy, L. Abe, A. Agabi, D. Mekarnia, T. Guillot, F. Schmider, P. Stee, Y. de Pra, M. Buttu, N. Crouzet , et al. (11 additional authors not shown)

Abstract: Photometric monitoring of Beta Pictoris in 1981 showed anomalous fluctuations of up to 4% over several days, consistent with foreground material transiting the stellar disk. The subsequent discovery of the gas giant planet Beta Pictoris b and the predicted transit of its Hill sphere to within 0.1 au projected distance of the planet provided an opportunity to search for the transit of a circumplane… ▽ More Photometric monitoring of Beta Pictoris in 1981 showed anomalous fluctuations of up to 4% over several days, consistent with foreground material transiting the stellar disk. The subsequent discovery of the gas giant planet Beta Pictoris b and the predicted transit of its Hill sphere to within 0.1 au projected distance of the planet provided an opportunity to search for the transit of a circumplanetary disk in this $21\pm 4$ Myr-old planetary system. Continuous broadband photometric monitoring of Beta Pictoris requires ground-based observatories at multiple longitudes to provide redundancy and to provide triggers for rapid spectroscopic followup. These observatories include the dedicated Beta Pictoris monitoring observatory bRing at Sutherland and Siding Springs, the ASTEP400 telescope at Concordia, and observations from the space observatories BRITE and Hubble Space Telescope. We search the combined light curves for evidence of short period transient events caused by rings and for longer term photometric variability due to diffuse circumplanetary material. We find no photometric event that matches with the event seen in November 1981, and there is no systematic photometric dimming of the star as a function of the Hill sphere radius. We conclude that the 1981 event was not caused by the transit of a circumplanetary disk around Beta Pictoris b. The upper limit on the long term variability of Beta Pictoris places an upper limit of $1.8\times 10^{22}$ g of dust within the Hill sphere. Circumplanetary material is either condensed into a non-transiting disk, is condensed into a disk with moons that has a small obliquity, or is below our detection threshold. This is the first time that a dedicated international campaign has mapped the Hill sphere transit of a gas giant extrasolar planet at 10 au. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 12 pages, 9 figures, 1 table, accepted for publication in A&A. Reduced data and reduction scripts on GitHub at https://github.com/mkenworthy/beta_pic_b_hill_sphere_transit

Journal ref: A&A 648, A15 (2021)