Large-amplitude Sgr A* near-infrared flares result from energy injection into electrons near the black hole event horizon. Astrometry data show continuous rotation of the emission region during bright flares, and corresponding rotation of the linear polarization angle. One broad class of physical flare models invokes magnetic reconnection. Here we show that such a scenario can arise in a general relativistic magnetohydrodynamic simulation of a magnetically arrested disc. Saturation of magnetic flux triggers eruption events, where magnetically dominated plasma is expelled from near the horizon and forms a rotating, spiral structure. Dissipation occurs via reconnection at the interface of the magnetically dominated plasma and surrounding fluid. This dissipation is associated with large increases in near-infrared emission in models of Sgr A*, with durations and amplitudes consistent with the observed flares. Such events occur at roughly the timescale to re-accumulate the magnetic flux from the inner accretion disc, 10h for Sgr A*. We study near-infrared observables from one sample event to show that the emission morphology tracks the boundary of the magnetically dominated region. As the region rotates around the black hole, the near-infrared centroid and linear polarization angle both undergo continuous rotation, similar to the behavior seen in Sgr A* flares.

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.

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.

The Galactic center black hole candidate Sgr A* is the best target for studies of low-luminosity accretion physics, including with near-infrared and submillimeter wavelength long baseline interferometry experiments. Here we compare images and spectra generated from a parameter survey of general relativistic MHD simulations to a set of radio to near-infrared observations of Sgr A*. Our models span the limits of weak and strong magnetization and use a range of sub-grid prescriptions for electron heating. We find two classes of scenarios can explain the broad shape of the submillimeter spectral peak and the highly variable near-infrared flaring emission. Weakly magnetized "disk-jet" models where most of the emission is produced near the jet wall, consistent with past work, as well as strongly magnetized (magnetically arrested disk) models where hot electrons are present everywhere. Disk-jet models are strongly depolarized at submillimeter wavelengths as a result of strong Faraday rotation, inconsistent with observations of Sgr A*. We instead favor the strongly magnetized models, which provide a good description of the median and highly variable linear polarization signal. The same models can also explain the observed mean Faraday rotation measure and potentially the polarization signals seen recently in Sgr A* near-infrared flares.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

In this paper, we explore a scenario where the dark matter is a mixture of interacting and non interacting species. Assuming dark matter-photon interactions for the interacting species, we find that the suppression of the matter power spectrum in this scenario can mimic that expected in the case of massive neutrinos. Our numerical studies include present limits from Planck Cosmic Microwave Background data, which render the strength of the dark matter photon interaction unconstrained when the fraction of interacting dark matter is small. Despite the large entangling between mixed dark matter and neutrino masses, we show that future measurements from the Dark Energy Instrument (DESI) could help in establishing the dark matter and the neutrino properties simultaneously, provided that the interaction rate is very close to its current limits and the fraction of interacting dark matter is at least of $\mathcal{O}\left(10\%\right)$. However, for that region of parameter space where a small fraction of interacting DM coincides with a comparatively large interaction rate, our analysis highlights a considerable degeneracy between the mixed dark matter parameters and the neutrino mass scale.

Mixed damping is a physical effect that occurs when a heavy species is coupled to a relativistic fluid which is itself free streaming. As a cross-case between collisional damping and free-streaming, it is crucial in the context of neutrino-dark matter interactions. In this work, we establish the parameter space relevant for mixed damping, and we derive an analytical approximation for the evolution of dark matter perturbations in the mixed damping regime to illustrate the physical processes responsible for the suppression of cosmological perturbations. Although extended Boltzmann codes implementing neutrino-dark matter scattering terms automatically include mixed damping, this effect has not been systematically studied. In order to obtain reliable numerical results, it is mandatory to reconsider several aspects of neutrino-dark matter interactions, such as the initial conditions, the ultra-relativistic fluid approximation and high order multiple moments in the neutrino distribution. Such a precise treatment ensures the correct assessment of the relevance of mixed damping in neutrino-dark matter interactions.

We study the small scale structure of axion dark matter in the post-inflationary scenario, which predicts the formation of low-mass, high-density clumps of gravitationally bound axions called axion miniclusters. To this end, we follow numerically the cosmological evolution of the axion field and the network of strings and domain walls until the density contrast is frozen. Our simulations, comprising up to $8192^3$ points, are the largest studies of the axion field evolution in the non-linear regime presented so far. Axitons, pseudo-breathers of the Klein-Gordon equation, are observed to form in our simulation at late times. Studying their properties analytically and numerically, we observe that in particular, the earliest axitons contribute to density perturbations at the typical length scale of miniclusters. We analyse the small scale structure of the density field, giving the correlation length, power spectrum and the distribution of high-density regions that will collapse into axion miniclusters. The final density field of our simulations can be used to calculate the minicluster mass fraction in simulations including gravity. In particular, we find that typical minicluster progenitors are smaller than previously thought and only of moderate, $\mathcal{O}(1)$ overdensity. We expect these miniclusters to have a rich sub-structure, emerging from small-scale fluctuations produced in the collapse of the string-wall network and from axitons.

The Planck experiment has measured the Cosmic Microwave Background (CMB) angular spectrum with unprecedented accuracy, and these results are likely to remain the state-of-the art for a decade or more. Since these measurements are going to be used to forecast Euclid's ability to test alternative dark matter models, we aim to set the final CMB constraint on the cross section of dark matter-photon scattering in the early universe. As the formalism has not been published in the literature yet, we present the key steps to derive the full set of linearized Einstein and Boltzmann equations and elucidate necessary modifications to the line of sight and the tight coupling approximation. These approximation schemes are commonly used in Boltzmann codes to reduce computational costs, but the tight coupling regime has been neglected in previous works related to dark matter-photon interactions. Our computations account for a dark matter sound speed, which was neglected in previous work. We find that its effect is negligible on the CMB spectra but fairly important on the linear matter power spectrum. We then analyze the Planck data, with and without the 2015 polarization data at high multipoles. In both cases our constraint on the ratio of the elastic scattering cross section to the dark matter mass becomes more stringent with respect to previous constraints. The change is about 35\% when the full Planck dataset is added (which includes E-mode polarization at high multipoles and the lensing potential power spectrum) and about 20\%. when we only account for the polarization data at low multipoles and the full temperature spectrum. Our most conservative limit reads as $\sigma_{\mathrm{DM}-\gamma} \leq 2.25\times 10^{-6}\,\sigma_\mathrm{Th}\left(m_\mathrm{DM}/\mathrm{GeV}\right)$ at 95\% confidence. We do not expect this limit to change by adding the recently published SPT data.