Source author record

M. Nowak

M. Nowak appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

Direct observations of exoplanet and brown dwarf companions with near-infrared interferometry, first enabled by the dual-field mode of VLTI/GRAVITY, provide unique measurements of the objects' orbital motions and atmospheric compositions. Here, we compile a homogeneous library of all exoplanet and brown dwarf K-band spectra observed by GRAVITY thus far. We re-reduced all the available GRAVITY dual-field high-contrast data and, where companions are detected, extract their ~2.0-2.4 $μ$m K-band contrast spectra. We then derived stellar model atmospheres for all employed flux references, which we used to convert the companion contrast into companion flux spectra. Solely from the resulting GRAVITY spectra, we extracted spectral types, spectral indices, and bulk physical properties for all companions. Finally, and with the help of age constraints from the literature, we also derived isochronal masses for most companions using evolutionary models. The resulting library contains R ~ 500 GRAVITY spectra of 39 substellar companions from late M to late T spectral types, including the entire L-T transition. Throughout this transition, a shift from CO-dominated late M- and L-type dwarfs to CH4-dominated T-type dwarfs can be observed in the K-band. The GRAVITY spectra alone constrain the objects' bolometric luminosity to typically within $\pm$0.15 dex. The derived isochronal masses agree with dynamical masses from the literature where available, except for HD 4113 c for which we confirm its previously reported potential underluminosity. Medium-resolution spectroscopy of substellar companions with GRAVITY provides insight into the carbon chemistry and the cloudiness of these objects' atmospheres. It also constrains these objects' bolometric luminosities which can yield measurements of their formation entropy if combined with dynamical masses, for instance from Gaia and GRAVITY astrometry.

preprint2025arXiv

Multi-modal atmospheric characterization of $β$ Pictoris b: Adding high-resolution continuum spectra from GRAVITY

We present the first VLTI/GRAVITY observations at R$_λ\sim 4000$ of $β$ Pic b. These four high S/N ($\sim$20) K-band spectra conserve both the pseudo-continuum and molecular absorption patterns. We analyze them with four self-consistent forward model grids (Exo-REM, ATMO, BT-Settl, Sonora) exploring $T_{\mathrm{eff}}$, log(g), metallicity, C/O, and $^{12}$CO/$^{13}$CO ratio. We also upgrade our forward modeling code ForMoSA to account for the data multi-modality and combine the GRAVITY epochs with published 1-5 $μ$m photometry, low- to medium-resolution spectra (0.9-7 $μ$m), and high-resolution echelle spectra (2.1-5.2 $μ$m). Sonora and Exo-REM are statistically preferred. Exo-REM yields $T_{\mathrm{eff}}$ $=1607.45^{+4.85}_{-6.20}$ K and log(g) $=4.46^{+0.02}_{-0.04}$ dex from GRAVITY alone, and $T_{\mathrm{eff}}$ $=1502.74^{+2.32}_{-2.14}$ K and log(g) $=4.00\pm0.01$ dex when including all datasets. Archival data significantly affect the retrieved parameters. C/O remains solar ($0.552^{+0.003}_{-0.002}$) while [M/H] reaches super-solar values (0.50$\pm$0.01). We report the first tentative constraint on log($^{12}$CO/$^{13}$CO) $\sim$1.12, though this remains inconclusive due to telluric residuals. Additionally, we estimate the luminosity to be log(L/L$_\odot$) $=-4.01^{+0.04}_{-0.05}$, implying a heavy-element content of up to $\sim$5% (20-80 M$_\oplus$) given the system age and dynamical mass measurements. Access to both continuum and molecular lines at K-band significantly impacts the metallicity, possibly owing to collision-induced absorption shaping the continuum. Echelle spectra do not dominate the final fit with respect to lower resolution data. Future multi-modal frameworks should include weighting schemes reflecting bandwidth and central wavelength coverage.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Multi-wavelength Variability of Sagittarius A* in July 2019

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.

preprint2021arXiv

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $μ$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b.

preprint2021arXiv

Deep Images of the Galactic Center with GRAVITY

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.

preprint2021arXiv

Detection of faint stars near SgrA* with GRAVITY

The spin of the supermassive black hole that resides at the Galactic Centre can in principle be measured by accurate measurements of the orbits of stars that are much closer to SgrA* than S2, the orbit of which recently provided the measurement of the gravitational redshift and the Schwarzschild precession. The GRAVITY near-infrared interferometric instrument combining the four 8m telescopes of the VLT provides a spatial resolution of 2-4 mas, breaking the confusion barrier for adaptive-optics-assisted imaging with a single 8-10m telescope. We used GRAVITY to observe SgrA* over a period of six months in 2019 and employed interferometric reconstruction methods developed in radio astronomy to search for faint objects near SgrA*. This revealed a slowly moving star of magnitude 18.9 in K band within 30mas of SgrA*. The position and proper motion of the star are consistent with the previously known star S62, which is at a substantially larger physical distance, but in projection passes close to SgrA*. Observations in August and September 2019 easily detected S29, with K-magnitude of 16.6, at approximately 130 mas from SgrA*. The planned upgrades of GRAVITY, and further improvements in the calibration, hold the promise of finding stars fainter than magnitude 19 at K.

preprint2021arXiv

Improved GRAVITY astrometric accuracy from modeling of optical aberrations

The GRAVITY instrument on the ESO VLTI pioneers the field of high-precision near-infrared interferometry by providing astrometry at the $10 - 100\,μ$as level. Measurements at such high precision crucially depend on the control of systematic effects. Here, we investigate how aberrations introduced by small optical imperfections along the path from the telescope to the detector affect the astrometry. We develop an analytical model that describes the impact of such aberrations on the measurement of complex visibilities. Our formalism accounts for pupil-plane and focal-plane aberrations, as well as for the interplay between static and turbulent aberrations, and successfully reproduces calibration measurements of a binary star. The Galactic Center observations with GRAVITY in 2017 and 2018, when both Sgr A* and the star S2 were targeted in a single fiber pointing, are affected by these aberrations at a level of less than 0.5 mas. Removal of these effects brings the measurement in harmony with the dual beam observations of 2019 and 2020, which are not affected by these aberrations. This also resolves the small systematic discrepancies between the derived distance $R_0$ to the Galactic Center reported previously.

preprint2021arXiv

The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances.

preprint2020arXiv

Detection of the Schwarzschild precession in the orbit of the star S2 near the Galactic centre massive black hole

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

preprint2020arXiv

Dynamically important magnetic fields near the event horizon of Sgr A*

We study the time-variable linear polarisation of Sgr A* during a bright NIR flare observed with the GRAVITY instrument on July 28, 2018. Motivated by the time evolution of both the observed astrometric and polarimetric signatures, we interpret the data in terms of the polarised emission of a compact region ('hotspot') orbiting a black hole in a fixed, background magnetic field geometry. We calculated a grid of general relativistic ray-tracing models, created mock observations by simulating the instrumental response, and compared predicted polarimetric quantities directly to the measurements. We take into account an improved instrument calibration that now includes the instrument's response as a function of time, and we explore a variety of idealised magnetic field configurations. We find that the linear polarisation angle rotates during the flare, which is consistent with previous results. The hotspot model can explain the observed evolution of the linear polarisation. In order to match the astrometric period of this flare, the near horizon magnetic field is required to have a significant poloidal component, which is associated with strong and dynamically important fields. The observed linear polarisation fraction of $\simeq 30\%$ is smaller than the one predicted by our model ($\simeq 50\%$). The emission is likely beam depolarised, indicating that the flaring emission region resolves the magnetic field structure close to the black hole.

preprint2020arXiv

Effect of Na doping on photoluminescence and laser stimulated nonlinear T optical features of ZnO nanostructures

In this work, nonlinear optical properties of Na: ZnO thin films (Na: ZnO) have been experimentally elaborated. The principal possibility to operate the nonlinear optical features using external laser beams is shown. The Na: ZnO films were synthesized by spray pyrolysis technique at a deposition temperature equal to about 400 C. XRD graph reveals that the grown films were polycrystalline in nature with a dominant peak corresponding to (0 0 2) plane. Despite the difference in the ionic radii of the Na and Zn, an angle shift in the XRD peak was not observed, whereas there was a significant change in peak intensity. The photoluminescence (PL) spectra resulted in three emission centres spectrally situated in violet, blue and green colour region due to the presence of native defect states in the forbidden energy gap. Second and third harmonic generation (SHG, THG) ex-periments stimulated by external coherent light beams show the existence of THG maxima for 15% Na doped Zn:O which was quite different with respect to SHG maximum for pure ZnO. The main innovation of this work is the possibility to change nonlinear optical susceptibility varying Na doping concentration and by coherent laser treatment, contrary to the previous works where these parameters have not been explored.

preprint2020arXiv

High resolution X-ray spectroscopy of Supergiant HMXB 4U1700$-$37 during the compact object eclipse

We present an analysis of the first observation of the iconic High Mass X-ray Binary \so with the \chandra High Energy Transmission Gratings during an X-ray eclipse. The goal of the observation was to study the structure/physical conditions in the clumpy stellar wind through high resolution spectroscopy. We find that: a) emission line brightness from K shell transitions, corresponding to near neutral species, directly correlates with continuum illumination. However, these lines do not greatly diminish during eclipse. This is readily explained if fluorescence K$α$ emission comes from the bulk of the wind. b) The highly ionised Fexxv and Fexxvi Ly$α$ diminish during eclipse. Thus, they must be produced in the vicinity of the compact object where $\log ξ>3$. c) to describe the emission line spectrum, the sum of two self consistent photo ionisation models with low ionisation ($\log ξ\sim -1$) and high ionisation ($\log ξ\sim 2.4$) is required. From their emission measures, the clump-to-interclump density ratio can be estimated to be $n_c/n_i\sim 300$. To fit the complex He-like \ion{Si}{xiii}{} profile, the plasma requires a broadening with $v_{\rm bulk}\sim 840$ km s$^{-1}$. Reproducing the observed $r\approx f$ line fluxes requires the addition of a third collisionally ionised plasma. d) Emission lines widths appear unresolved at the \textsc{hetg} gratings resolution with exception of Silicon. There is no clear radial segregation between (quasi)neutral and ionised species, consistent with cold wind clumps interspersed in a hot rarefied interclump medium.

preprint2020arXiv

K-Stacker, an algorithm to hack the orbital parameters of planets hidden in high-contrast imaging. First applications to VLT SPHERE multi-epoch observations

Recent high-contrast imaging surveys, looking for planets in young, nearby systems showed evidence of a small number of giant planets at relatively large separation beyond typically 20 au where those surveys are the most sensitive. Access to smaller physical separations between 5 and 20 au is the next step for future planet imagers on 10 m telescopes and ELTs in order to bridge the gap with indirect techniques (radial velocity, transit, astrometry with Gaia). In that context, we recently proposed a new algorithm, Keplerian-Stacker, combining multiple observations acquired at different epochs and taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of SPHERE even when a planet remains undetected at one epoch. Here, we validate the K-Stacker algorithm performances on real SPHERE datasets, to demonstrate its resilience to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns in high-contrast imaging to search for planets in emitted and reflected light. Results. We show that K-Stacker achieves high success rate when the SNR of the planet in the stacked image reaches 7. The improvement of the SNR ratio goes as the square root of the total exposure time. During the blind test and the redetection of HD 95086 b, and betaPic b, we highlight the ability of K-Stacker to find orbital solutions consistent with the ones derived by the state of the art MCMC orbital fitting techniques, confirming that in addition to the detection gain, K-Stacker offers the opportunity to characterize the most probable orbital solutions of the exoplanets recovered at low signal to noise.

preprint2020arXiv

Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e

Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

preprint2020arXiv

The Chandra High Resolution X-ray Spectrum of Quiescent Emission from Sgr A*

In quiescence, Sgr A* is surprisingly dim, shining 100,000 times less than expected for its environment. This problem has motivated a host of theoretical models to explain radiatively inefficient accretion flows (RIAFs). The Chandra Galactic Center (GC) X-ray Visionary Program obtained approximately 3 Ms (one month) of Chandra HETG data, offering the only opportunity to examine the quiescent X-ray emission of Sgr A* with high resolution spectroscopy. Utilizing custom background regions and filters for removing overlapping point sources, this work provides the first ever look at stacked HETG spectra of Sgr A*. We model the background datasets with a cubic spline and fit the unbinned Sgr A* spectra with a simple parametric model of a power law plus Gaussian lines under the effects of interstellar extinction. We detect a strong 6.7 keV iron emission line in the HEG spectra and a 3.1 keV emission line in the MEG spectra. In all cases, the line centroids and equivalent widths are consistent with those measured from low-resolution CCD spectra. An examination of the unbinned, stacked HEG+/-1 spectrum reveals fine structure in the iron line complex. In addition to resolving the resonant and forbidden lines from He-like iron, there are apparent emission features arising with higher statistical significance at lower energy, potentially associated with FeXX-XXIV ions in a ~1 keV plasma arising near the Bondi radius of Sgr A*. With this work, we release the cleaned and stacked Sgr A* and background HETG spectra to the public as a special legacy dataset.

preprint2020arXiv

The Flux Distribution of Sgr A*

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

preprint2020arXiv

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

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

preprint2019arXiv

Peering into the formation history of beta Pictoris b with VLTI/GRAVITY long baseline interferometry

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $β$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $β$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $β$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

preprint2016arXiv

An Ultra-Fast X-ray Disk Wind in the Neutron Star Binary GX 340+0

We present a spectral analysis of a brief Chandra/HETG observation of the neutron star low-mass X-ray binary GX~340+0. The high-resolution spectrum reveals evidence of ionized absorption in the Fe K band. The strongest feature, an absorption line at approximately 6.9 keV, is required at the 5 sigma level of confidence via an F-test. Photoionization modeling with XSTAR grids suggests that the line is the most prominent part of a disk wind with an apparent outflow speed of v = 0.04c. This interpretation is preferred at the 4 sigma level over a scenario in which the line is H-like Fe XXVI at a modest red-shift. The wind may achieve this speed owing to its relatively low ionization, enabling driving by radiation pressure on lines; in this sense, the wind in GX 340+0 may be the stellar-mass equivalent of the flows in broad absorption line quasars (BALQSOs). If the gas has a unity volume filling factor, the mass ouflow rate in the wind is over 10^-5 Msun/year, and the kinetic power is nearly 10^39 erg/s (or, 5-6 times the radiative Eddington limit for a neutron star). However, geometrical considerations - including a small volume filling factor and low covering factor - likely greatly reduce these values.

preprint2016arXiv

Nulling Data Reduction and On-Sky Performance of the Large Binocular Telescope Interferometer

The Large Binocular Telescope Interferometer (LBTI) is a versatile instrument designed for high-angular resolution and high-contrast infrared imaging (1.5-13 microns). In this paper, we focus on the mid-infrared (8-13 microns) nulling mode and present its theory of operation, data reduction, and on-sky performance as of the end of the commissioning phase in March 2015. With an interferometric baseline of 14.4 meters, the LBTI nuller is specifically tuned to resolve the habitable zone of nearby main-sequence stars, where warm exozodiacal dust emission peaks. Measuring the exozodi luminosity function of nearby main-sequence stars is a key milestone to prepare for future exoEarth direct imaging instruments. Thanks to recent progress in wavefront control and phase stabilization, as well as in data reduction techniques, the LBTI demonstrated in February 2015 a calibrated null accuracy of 0.05% over a three-hour long observing sequence on the bright nearby A3V star beta Leo. This is equivalent to an exozodiacal disk density of 15 to 30 zodi for a Sun-like star located at 10pc, depending on the adopted disk model. This result sets a new record for high-contrast mid-infrared interferometric imaging and opens a new window on the study of planetary systems.

preprint2016arXiv

Reaching sub-milimag photometric precision on Beta Pictoris with a nanosat: the PicSat mission

PicSat is a nanosatellite currently being developed to observe the transit of the giant planet \b{eta} Pictoris, expected some time between July 2017 and June 2018. The mission is based on a Cubesat architecture, with a small but ambitious 2 kg opto-mechanical payload specifically designed for high precision photometry. The satellite will be launched in early 2017, probably on a 600 km Sun synchronous orbit. The main objective of the mission is the constant monitoring of the brightness of Pic at an unprecedented combination of reliability and precision (200 ppm per hour, with interruptions of at most 30 minutes) to finely characterize the transiting exoplanet and detect exocomets in the Pictoris system. To achieve this difficult objective, the payload is designed with a 3.5 cm effective aperture telescope which injects the light in a single-mode optical fiber linked to an avalanche photodioode. A two-axis piezoelectric actuation system, driven by a tailor-made feedback loop control algorithm, is used to lock the fiber on the center of the star in the focal plane. These actuators complement the attitude determination and control system of the satellite to maintain the sub-arcsecond pointing accuracy required to reach the excellent level of photometric precision. Overall, the mission raises multiple very difficult challenges: high temperature stability of the avalanche detector (achieved with a thermoelectric colling device), high pointing accuracy and stability, and short timeframe for the development.

preprint2015arXiv

K-Stacker, a new way of detecting and characterizing exoplanets with high contrast imaging instruments

This year, a second generation of coronagraphs dedicated to high-contrast direct imaging of exoplanets is starting operations. Among them, SPHERE, installed at the focus of the UT3 Very Large Telescope, reaches unprecedented contrast ratios up to $10^{-6}$ -$ 10^{-7}$, using eXtreme Adaptive Optics and the Angular Differential Imaging (ADI) techniques. In this paper, we present a new method called Keplerian-Stacker that improves the detection limit of high contrast instruments like SPHERE, by up to a factor of 10. It consists of observing a star on a long enough period to let a hypothetical planet around that star move along its orbit. Even if in each individual observation taken during one night, we do not detect anything, we show that it is possible, using an optimization algorithm, to re-center the images according to keplerian motions (ex: 10-100 images taken over a long period of typically 1-10 years) and detect planets otherwise unreachable. This method can be used in combination with the ADI technics (or possibly any other high contrast data reduction method) to improve the Signal to Noise Ratio in each individual image, and to further improve the global detection limit. It also directly provides orbital parameters of the detected planets, as a by-product of the optimization algorithm.

preprint2015arXiv

NuSTAR and Suzaku observations of the hard state in Cygnus X-1: locating the inner accretion disk

We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. We use the reflection models of Garcia et al. (2014) to simultaneously measure the black hole spin, disk inner radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3σ confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection model.

preprint2014arXiv

Shell-Shocked: The Interstellar Medium Near Cygnus X-1

We conduct a detailed case-study of the interstellar shell near the high-mass X-ray binary, Cygnus X-1. We present new WIYN optical spectroscopic and Chandra X-ray observations of this region, which we compare with detailed MAPPINGS III shock models, to investigate the outflow powering the shell. Our analysis places improved, physically motivated constraints on the nature of the shockwave and the interstellar medium (ISM) it is plowing through. We find that the shock is traveling at less than a few hundred km/s through a low-density ISM (< 5 cm^-3). We calculate a robust, 3 sigma upper limit to the total, time-averaged power needed to drive the shockwave and inflate the bubble, < 2 x 10^38 erg/s. We then review possible origins of the shockwave. We find that a supernova origin to the shockwave is unlikely and that the black hole jet and/or O-star wind can both be central drivers of the shockwave. We conclude that the source of the Cygnus X-1 shockwave is far from solved.

preprint2014arXiv

Spitzer/infrared spectrograph investigation of MIPSGAL 24 μm compact bubbles : Low resolution observations

We present Spitzer/IRS low resolution observations of 11 compact circumstellar bubbles from the MIPSGAL 24 μm Galactic Plane Survey. We find that this set of MIPSGAL bubbles (MBs) is divided into two categories, and that this distinction correlates with the morphologies of the MBs in the mid- IR. The four MBs with central sources in the mid-IR exhibit dust-rich, low excitation spectra, and their 24 μm emission is accounted for by the dust continuum. The seven MBs without central sources in the mid-IR have spectra dominated by high excitation gas lines (e.g., [O IV] 26.0 μm, [Ne V] 14.3 and 24.3 μm, [Ne III] 15.5 μm), and the [O IV] line accounts for 50 to almost 100% of the 24 μm emission in five of them. In the dust-poor MBs, the [Ne V] and [Ne III] line ratios correspond to high excitation conditions. Based on comparisons with published IRS spectra, we suggest that the dust-poor MBs are highly excited planetary nebulae with peculiar white dwarfs (e.g., [WR], novae) at their centers. The central stars of the four dust-rich MBs are all massive star candidates. Dust temperatures range from 40 to 100 K in the outer shells. We constrain the extinction along the lines of sight from the IRS spectra. We then derive distance, dust masses, and dust production rate estimates for these objects. These estimates are all consistent with the nature of the central stars. We summarize the identifications of MBs made to date and discuss the correlation between their mid-IR morphologies and natures. Candidate Be/B[e]/LBV and WR stars are mainly "rings" with mid-IR central sources, whereas PNe are mostly "disks" without mid-IR central sources. Therefore we expect that most of the 300 remaining unidentified MBs will be classified as PNe.

preprint2014arXiv

The Large Area Detector of LOFT: the Large Observatory for X-ray Timing

LOFT (Large Observatory for X-ray Timing) is one of the five candidates that were considered by ESA as an M3 mission (with launch in 2022-2024) and has been studied during an extensive assessment phase. It is specifically designed to perform fast X-ray timing and probe the status of the matter near black holes and neutron stars. Its pointed instrument is the Large Area Detector (LAD), a 10 m 2 -class instrument operating in the 2-30keV range, which holds the capability to revolutionise studies of variability from X-ray sources on the millisecond time scales. The LAD instrument has now completed the assessment phase but was not down-selected for launch. However, during the assessment, most of the trade-offs have been closed leading to a robust and well documented design that will be re- proposed in future ESA calls. In this talk, we will summarize the characteristics of the LAD design and give an overview of the expectations for the instrument capabilities.

preprint2014arXiv

The Large Observatory For x-ray Timing

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

preprint2013arXiv

A Census of X-ray gas in NGC 1068: Results from 450ks of Chandra HETG Observation

We present models for the X-ray spectrum of the Seyfert 2 galaxy NGC 1068. These are fitted to data obtained using the High Energy Transmission Grating (HETG) on the Chandra X-ray observatory. The data show line and radiative recombination continuum (RRC) emission from a broad range of ions and elements. The models explore the importance of excitation processes for these lines including photoionization followed by recombination, radiative excitation by absorption of continuum radiation and inner shell fluorescence. The models show that the relative importance of these processes depends on the conditions in the emitting gas, and that no single emitting component can fit the entire spectrum. In particular, the relative importance of radiative excitation and photoionization/recombination differs according to the element and ion stage emitting the line. This in turn implies a diversity of values for the ionization parameter of the various components of gas responsible for the emission, ranging from log(xi)=1 -- 3. Using this, we obtain an estimate for the total amount of gas responsible for the observed emission. The mass flux through the region included in the HETG extraction region is approximately 0.3 Msun/yr assuming ordered flow at the speed characterizing the line widths. This can be compared with what is known about this object from other techniques.

preprint2012arXiv

A Large Area Detector proposed for the Large Observatory for X-ray Timing (LOFT)

The Large Observatory for X-ray Timing (LOFT) is one of the four candidate ESA M3 missions considered for launch in the 2022 time-frame. It is specifically designed to perform fast X-ray timing and probe the status of the matter near black holes and neutron stars. The LOFT scientific payload is composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a 10 m2-class pointed instrument with 20 times the collecting area of the best past timing missions (such as RXTE) over the 2-30 keV range, which holds the capability to revolutionize studies of X-ray variability down to the millisecond time scales. Its ground-breaking characteristic is a low mass per unit surface, enabling an effective area of ~10 m^2 (@10 keV) at a reasonable weight. The development of such large but light experiment, with low mass and power per unit area, is now made possible by the recent advancements in the field of large-area silicon detectors - able to time tag an X-ray photon with an accuracy <10 μs and an energy resolution of ~260 eV at 6 keV - and capillary-plate X-ray collimators. In this paper, we will summarize the characteristics of the LAD instrument and give an overview of its capabilities.

preprint2012arXiv

LOFT: the Large Observatory For X-ray Timing

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

preprint2010arXiv

Exploring Accretion and Disk-Jet Connections in the LLAGN M81*

We report on a year-long effort to monitor the central supermassive black hole in M81 in the X-ray and radio bands. Using Chandra and the VLA, we obtained quasi-simultaneous observations of M81* on seven occasions during 2006. The X-ray and radio luminosity of M81* are not strongly correlated on the approximately 20-day sampling timescale of our observations, which is commensurate with viscous timescales in the inner flow and orbital timecales in a radially-truncated disk. This suggests that short-term variations in black hole activity may not be rigidly governed by the "fundamental plane", but rather adhere to the plane in a time-averaged sense. Fits to the X-ray spectra of M81* with bremsstrahlung models give temperatures that are inconsistent with the outer regions of very simple advection-dominated inflows. However, our results are consistent with the X-ray emission originating in a transition region where a truncated disk and advective flow may overlap. We discuss our results in the context of models for black holes accreting at small fractions of their Eddington limit, and the fundamental plane of black hole accretion.

preprint2009arXiv

Constraining Jet/Disk Geometry and Radiative Processes in Stellar Black Holes XTE J1118+480 and GX 339-4

We present results from modeling of quasi-simultaneous broad band (radio through X-ray) observations of the galactic stellar black hole (BH) transient X-ray binary (XRB) systems XTE J1118+480 and GX 339-4 using an irradiated disc + compact jet model. In addition to quantifying the physical properties of the jet, we have developed a new irradiated disc model which also constrains the geometry and temperature of the outer accretion disc by assuming a disc heated by viscous energy release and X-ray irradiation from the inner regions. For the source XTE J1118+480, which has better spectral coverage of the two in optical and near-IR (OIR) wavelengths, we show that the entire broad band continuum can be well described by an outflow-dominated model + an irradiated disc. The best-fit radius of the outer edge of the disc is consistent with the Roche lobe geometry of the system, and the temperature of the outer edge of the accretion disc is similar to those found for other XRBs. Irradiation of the disc by the jet is found to be negligible for this source. For GX 339-4, the entire continuum is well described by the jet-dominated model only, with no disc component required. For the two XRBs, which have very different physical and orbital parameters and were in different accretion states during the observations, the sizes of the jet base are similar and both seem to prefer a high fraction of non-thermal electrons in the acceleration/shock region and a magnetically dominated plasma in the jet. These results, along with recent similar results from modeling other galactic XRBs and AGNs, may suggest an inherent unity in diversity in the geometric and radiative properties of compact jets from accreting black holes.