Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

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.

We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

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.

Aims. The fueling and feedback of Active Galactic Nuclei (AGNs) are important in understanding the co-evolution between black holes and host galaxies. Mergers are thought to have the capability to bring gas inwards and ignite nuclear activity, especially for more powerful AGNs. However, there is still significant ongoing debate on whether mergers can trigger AGNs and, if they do, whether mergers are a significant triggering mechanism. Methods. We select a low-redshift ($0.005<z<0.1$) sample from the Sloan Digital Sky Survey (SDSS) and a high-redshift ($0 < z < 0.6$) sample from the Galaxy And Mass Assembly (GAMA) survey. We take advantage of the convolutional neural network technique to identify mergers. We use mid-infrared (MIR) color cut and optical emission line diagnostics to classify AGNs. We also include Low Excitation Radio Galaxies (LERGs) to investigate the connection between mergers and low accretion rate AGNs. Results. We find that AGNs are more likely to be found in mergers than non-mergers, with an AGN excess up to $1.81\pm{0.16}$, suggesting that mergers can trigger AGNs. We also find the fraction of mergers in AGNs is higher than that in non-AGN controls, for both MIR and optically selected AGNs, as well as LERGs, with values between $16.40\pm{0.5}\%$ and $39.23\pm{2.10}\%$, implying a non-negligible to potentially significant role of mergers in triggering AGNs. This merger fraction in AGNs increases as stellar mass increases which supports the idea that mergers are more important for triggering AGNs in more massive galaxies. In terms of merger fraction as a function of AGN power we find a positive trend for MIR selected AGNs and a complex trend for optically selected AGNs, which we interpret under an evolutionary scenario proposed by previous studies. In addition, obscured MIR selected AGNs are more likely to be hosted in mergers than unobscured MIR selected AGNs.

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.

As the identified number of Jupiter's moons has skyrocketed to 79, some of them have been regrouped. In this work, we continue to identify the potential distributions of the physical characteristics of Jupiter's irregular moons. By using nonparametric Kolmogorov-Smirnov tests, we verified more than 20 commonly used distributions and found that surprisingly, almost all the physical characteristics (i.e., the equatorial radius, equatorial circumference, circumference, volume, mass, surface gravity and escape velocity) of the moons in the Ananke and Carme groups follow log-logistic distributions. Additionally, more than half of the physical characteristics of the moons in the Pasiphae group are theoretically subject to this type of distribution. Combined with strict analytical derivations, it is increasingly clear and possible to anticipate that the physical characteristics of most irregular moons follow log-logistic distributions with the discovery of an increasing number of Jupiter's irregular moons.

From the statistical point of view, this paper mainly emphasizes the orbital distribution laws of Jupiter's moons, most of which are located in Ananke group, Carme group and Pasiphae group. By comparing 19 known continuous distributions, it is verified that there are suitable distribution functions to describe the distribution of these natural satellites. For each distribution type, interval estimation is used to estimate the corresponding parameter values. At a given significant level, one-sample Kolmogorov-Smirnov nonparametric test is applied to verify the specified distribution, and we often select the one with the largest $p$-value. The results show that all the semi-major axis, mean inclination and the orbital period of the moons in Ananke group and Carme group obey the Stable distribution. In addition, according to Kepler's third planetary motion law, and by comparing the theoretically calculated best-fit cumulative distribution function (CDF) with the observed CDF, we demonstrate that the theoretical distribution is in good agreement with the empirical distribution. Therefore, these characteristics of Jupiter's moons are indeed very likely to follow some specific distribution laws, and it will be possible to use these laws to help study certain features of poorly investigated moons or even predict undiscovered ones.

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.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

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.

We present a measurement of the Hubble constant made using geometric distance measurements to megamaser-hosting galaxies. We have applied an improved approach for fitting maser data and obtained better distance estimates for four galaxies previously published by the Megamaser Cosmology Project: UGC 3789, NGC 6264, NGC 6323, and NGC 5765b. Combining these updated distance measurements with those for the maser galaxies CGCG 074-064 and NGC 4258, and assuming a fixed velocity uncertainty of 250 km s$^{-1}$ associated with peculiar motions, we constrain the Hubble constant to be $H_0 = 73.9 \pm 3.0$ km s$^{-1}$ Mpc$^{-1}$ independent of distance ladders and the cosmic microwave background. This best value relies solely on maser-based distance and velocity measurements, and it does not use any peculiar velocity corrections. Different approaches for correcting peculiar velocities do not modify $H_0$ by more than ${\pm}1{\sigma}$, with the full range of best-fit Hubble constant values spanning 71.8-76.9 km s$^{-1}$ Mpc$^{-1}$. We corroborate prior indications that the local value of $H_0$ exceeds the early-Universe value, with a confidence level varying from 95-99% for different treatments of the peculiar velocities.

As part of the survey component of the Megamaser Cosmology Project, we have discovered a disk megamaser system in the galaxy CGCG 074-064. Using the GBT and the VLA, we have obtained spectral monitoring observations of this maser system at a monthly cadence over the course of two years. We find that the systemic maser features display line-of-sight accelerations of ~4.4 km s$^{-1}$ yr$^{-1}$ that are nearly constant with velocity, while the high-velocity maser features show accelerations that are consistent with zero. We have also used the HSA to make a high-sensitivity VLBI map of the maser system in CGCG 074-064, which reveals that the masers reside in a thin, edge-on disk with a diameter of ~1.5 mas (0.6 pc). Fitting a three-dimensional warped disk model to the data, we measure a black hole mass of $2.42^{+0.22}_{-0.20} \times 10^7$ M$_{\odot}$ and a geometric distance to the system of $87.6^{+7.9}_{-7.2}$ Mpc. Assuming a CMB-frame recession velocity of $7308 \pm 150$ km s$^{-1}$, we constrain the Hubble constant to $H_0 = 81.0^{+7.4}_{-6.9}$ (stat.) $\pm 1.4$ (sys.) km s$^{-1}$ Mpc$^{-1}$.

Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.

We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_\chi$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_\chi > 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_\chi$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$.

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 $\beta$ 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 $\beta$ 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 $\beta$ 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.

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

We present VLBI images and kinematics of water maser emission in three active galaxies: NGC 5728, Mrk 1, and IRAS 08452-0011. IRAS 08452-0011 is a triple-peaked H2O megamaser, consistent with a Keplerian rotating disk, indicating a black hole mass of (3.2+/-0.2)x10^7 M_sun. NGC 5728 and Mrk 1 display double-peaked spectra and VLBI imaging reveal complicated gas kinematics, which do not allow for a robust determination of black hole mass. The two systems are either gas disks perturbed by AGN winds or part of outflows. We find that disturbed morphology and kinematics are a ubiquitous feature of all double-peaked maser systems, implying that these maser sources may reside in environments where AGN winds are prominent at ~1 pc scale and have significant impact on the masing gas. Such AGN tend to have black hole mass M_BH < 8x10^6 M_sun and Eddington ratios lambda_Edd >~ 0.1, while the triple-peaked megamasers show an opposite trend.

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

Aims. We aim to study the far-infrared radio correlation (FIRC) at 150 MHz in the local Universe (at a median redshift z~0:05) and improve the use of the rest-frame 150-MHz luminosity, L150, as a star-formation rate (SFR) tracer, which is unaffected by dust extinction. Methods. We cross-match the 60-um selected Revised IRAS Faint Source Survey Redshift (RIFSCz) catalogue and the 150-MHz selected LOFAR value-added source catalogue in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring Field. We estimate L150 for the cross-matched sources and compare it with the total infrared (IR) luminosity, LIR, and various SFR tracers. Results. We find a tight linear correlation between log L150 and log LIR for star-forming galaxies, with a slope of 1.37. The median qIR value (defined as the logarithm of the LIR to L150 ratio) and its rms scatter of our main sample are 2.14 and 0.34, respectively. We also find that log L150 correlates tightly with the logarithm of SFR derived from three different tracers, i.e., SFR_Halpha based on the Halpha line luminosity, SFR_60 based on the rest-frame 60-um luminosity and SFR_IR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (Lsun) = 1.35(0.06) x log SFR_Halpha (Msun/yr) + 3.20(0.06), log L150 (Lsun) = 1.31(0.05) x log SFR_60 (Msun/yr) + 3.14(0.06), and log L150 (Lsun) = 1.37(0.05) x log SFR_IR (Msun/yr) + 3.09(0.05), which show excellent agreement with each other.

Precise measurements of the S-stars orbiting SgrA* have set strong constraints on the nature of the compact object at the centre of the Milky Way. The presence of a black hole in that region is well established, but its neighboring environment is still an open debate. In that respect, the existence of dark matter in that central region may be detectable due to its strong signatures on the orbits of stars: the main effect is a Newtonian precession which will affect the overall pericentre shift of S2, the latter being a target measurement of the GRAVITY instrument. The exact nature of this dark matter (e.g., stellar dark remnants or diffuse dark matter) is unknown. This article assumes it to be an scalar field of toroidal distribution, associated with ultra-light dark matter particles, surrounding the Kerr black hole. Such a field is a form of "hair" expected in the context of superradiance, a mechanism that extracts rotational energy from the black hole. Orbital signatures for the S2 star are computed and shown to be detectable by GRAVITY. The scalar field can be constrained because the variation of orbital elements depends both on the relative mass of the scalar field to the black hole and on the field mass coupling parameter.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

Employing multiple pulsars and using an appropriate algorithm to establish ensemble pulsar timescale can reduce the influences of various noises on the long-term stability of pulsar timescale, compared to a single pulsar. However, due to the low timing precision and the significant red noises of some pulsars, their participation in the construction of ensemble pulsar timescale is often limited. Inspired by the principle of solving non-stationary sequence modeling using co-integration theory, we puts forward an algorithm based on the co-integration theory to establish ensemble pulsar timescale. It is found that this algorithm can effectively suppress some noise sources if a co-integration relationship between different pulsar data exist. Different from the classical weighted average algorithm, the co-integration method provides the chances of the pulsar with significant red noises to attend the establishment of ensemble pulsar timescale. Based on the data from the North American Nanohertz Observatory for Gravitational Waves, we found that the co-integration algorithm can successfully reduce several timing noises and improve the long-term stability of the ensemble pulsar timescale.

In a dynamically relaxed cluster around a massive black hole a dense stellar cusp of old stars is expected to form. Previous observations showed a relative paucity of red giant stars within the central 0.5 pc in the Galactic Center. By co-adding spectroscopic observations taken over a decade, we identify new late-type stars, including the first five warm giants (G2-G8III), within the central 1 arcsec 2 (0.04 {\times} 0.04 pc^2) of the Galaxy. Our findings increase the number of late-type stars to 21, of which we present deep spectra for 16. The updated star count, based on individual spectral classification, is used to reconstruct the surface density profile of giant stars. Our study, for the first time, finds a cusp in the surface number density of the spectroscopically identified old (>3 Gyr) giants population (m K<17) within 0.02-0.4 pc described by a single power law with an exponent {\Gamma}= 0.34 {\pm} 0.04.

The Galactic Center black hole Sgr A* is the archetypical example of an underfed massive black hole. The extremely low accretion rate can be understood in radiatively inefficient accretion flow models. Testing those models has proven to be difficult due to the lack of suitable probes. Radio and submm polarization measurements constrain the flow very close to the event horizon. X-ray observations resolving the Bondi radius yield an estimate roughly four orders of magnitude further out. Here, we present a new, indirect measurement of the accretion flow density at intermediate radii. We use the dynamics of the gas cloud G2 to probe the ambient density. We detect the presence of a drag force slowing down G2 with a statistical significance of approx 9 {\sigma}. This probes the accretion flow density at around 1000 Schwarzschild radii and yields a number density of approx. 4 x 10^3 cm^-3. Self-similar accretion models where the density follows a power law radial profile between the inner zone and the Bondi radius have predicted similar values.

NU Ori is a massive spectroscopic and visual binary in the Orion Nebula Cluster, with 4 components: Aa, Ab, B, and C. The B0.5 primary (Aa) is one of the most massive B-type stars reported to host a magnetic field. We report the detection of a spectroscopic contribution from the C component in high-resolution ESPaDOnS spectra, which is also detected in a Very Large Telescope Interferometer (VLTI) dataset. Radial velocity (RV) measurements of the inner binary (designated Aab) yield an orbital period of 14.3027(7) d. The orbit of the third component (designated C) was constrained using both RVs and interferometry. We find C to be on a mildly eccentric 476(1) d orbit. Thanks to spectral disentangling of mean line profiles obtained via least-squares deconvolution we show that the Zeeman Stokes $V$ signature is clearly associated with C, rather than Aa as previously assumed. The physical parameters of the stars were constrained using both orbital and evolutionary models, yielding $M_{\rm Aa} = 14.9 \pm 0.5 M_\odot$, $M_{\rm Ab} = 3.9 \pm 0.7 M_\odot$, and $M_{\rm C} = 7.8 \pm 0.7 M_\odot$. The rotational period obtained from longitudinal magnetic field $\langle B_z \rangle$ measurements is $P_{\rm rot} = 1.09468(7)$ d, consistent with previous results. Modeling of $\langle B_z \rangle$ indicates a surface dipole magnetic field strength of $\sim 8$ kG. NU Ori C has a magnetic field strength, rotational velocity, and luminosity similar to many other stars exhibiting magnetospheric H$\alpha$ emission, and we find marginal evidence of emission at the expected level ($\sim$1% of the continuum).

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

We compile multi-wavelength data from ultraviolet to infrared (IR) bands as well as redshift and source-type information for a large sample of 178,341 sources in the Hawaii-Hubble Deep Field-North field. A total of 145,635 sources among the full sample are classified/treated as galaxies and have redshift information available. We derive physical properties for these sources utilizing the spectral energy distribution fitting code CIGALE that is based on Bayesian analysis. Through various consistency and robustness check, we find that our stellar-mass and star-formation rate (SFR) estimates are reliable, which is mainly due to two facts. First, we adopt the most updated and accurate redshifts and point spread function-matched photometry; and second, we make sensible parameter choices with the CIGALE code and take into account influences of mid-IR/far-IR data, star-formation history models, and AGN contribution. We release our catalog of galaxy properties publicly (including, e.g., redshift, stellar mass, SFR, age, metallicity, dust attenuation), which is the largest of its kind in this field and should facilitate future relevant studies on formation and evolution of galaxies.

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

We report the first detection of the Galactic Centre massive black hole, Sgr~A*, in the far infrared. Our measurements were obtained with PACS on board the \emph{Herschel} satellite at $100~\mathrm{\mu m}$ and $160~\mathrm{\mu m}$. While the warm dust in the Galactic Centre is too bright to allow for a direct detection of Sgr~A*, we measure a significant and simultaneous variation of its flux of $\Delta F_{\nu\widehat{=}160 ~\mathrm{\mu m}} = (0.27\pm0.06)~\mathrm{Jy}$ and $\Delta F_{\nu\widehat{=}100 ~\mathrm{\mu m}}= (0.16\pm0.10)~\mathrm{Jy}$ during one observation. The significance level of the $160 ~\mathrm{\mu m}$ band variability is $4.5\sigma$ and the corresponding $100 ~\mathrm{\mu m}$ band variability is significant at $1.6\sigma$. We find no example of an equally significant false positive detection. Conservatively assuming a variability of $25\%$ in the FIR, we can provide upper limits to the flux. Comparing the latter with theoretical models we find that 1D RIAF models have difficulties explaining the observed faintness. However, the upper limits are consistent with modern ALMA and VLA observations. Our upper limits provide further evidence for a spectral peak at $\sim 10^{12} ~ \mathrm{Hz}$ and constrain the number density of $\gamma \sim 100$ electrons in the accretion disk and or outflow.

Water megamasers from circumnuclear disks in galaxy centers provide the most accurate measurements of supermassive black hole masses and uniquely probe the sub-parsec accretion processes. At the same time, these systems offer independent crucial constraints of the Hubble Constant in the nearby universe, and thus, the arguably best single constraint on the nature of dark energy. The chances of finding these golden standards are however abysmally low, at an overall =< 3% for any level of water maser emission detected at 22 GHz, and =< 1% for those exhibiting the disk-like configuration. We provide here a thorough summary of the current state of the detection of water megamaser disks, along with a novel investigation of the likelihood of increasing their detection rates based on a multivariate parameter analysis of optical and mid-infrared photometric properties of the largest database of galaxies surveyed for the 22 GHz emission. We find that galaxies with water megamaser emission tend to associate with strong emission in all WISE mid-infrared wavelengths, with the strongest enhancement in the W4 band, at 22 micron, as well as with previously proposed and newly found indicators of AGN strength in the mid-infrared, such as red W1-W2 and W1-W4 colors, and the integrated mid-infrared luminosity of the host galaxy. These trends offer a potential boost of the megamaser detection rates to 6 - 15%, or a factor of 2-8 relative to the current rates, depending on the chosen sample selection criteria, while fostering real chances for discovering >= 20 new megamaser disks.

H2O maser disks with Keplerian rotation in active galactic nuclei offer a clean way to determine accurate black hole mass and the Hubble constant. An important assumption made in using a Keplerian H2O maser disk for measuring the black hole mass and the Hubble constant is that the disk mass is negligible compared to the black hole mass. To test this assumption, a simple and useful model can be found in Hure et al. (2011). In this work, the authors apply a linear disk model to a position-dynamical mass diagram and re-analyze position-velocity data from H2O maser disks associated with active galactic nuclei. They claim that a maser disk with nearly perfect Keplerian rotation could have disk mass comparable to the black hole mass. This would imply that ignoring the effects of disk self-gravity can lead to large systematic errors in the measurement of black hole mass and the Hubble constant. We examine their methods and find that their large estimated disk masses of Keplerian disks are likely the result of their use of projected instead of 3-dimensional position and velocity information. To place better constraints on the disk masses of Keplerian maser systems, we incorporate disk self-gravity into a 3-dimensional Bayesian modelling program for maser disks and also evaluate constraints based on the physical conditions for disks which support water maser emission. We find that there is little evidence that disk masses are dynamically important at the ~<1% level compared to the black holes.

In this research note, we present a set of distortion solutions that may be used to correct geometric optical distortion in images taken with the S13 camera of the NACO adaptive optics imager.

In this paper, we describe the XENON100 data analyses used to assess the target-intrinsic background sources radon ($^{222}$Rn), thoron ($^{220}$Rn) and krypton ($^{85}$Kr). We detail the event selections of high-energy alpha particles and decay-specific delayed coincidences. We derive distributions of the individual radionuclides inside the detector and quantify their abundances during the main three science runs of the experiment over a period of $\sim$ 4 years, from January 2010 to January 2014. We compare our results to external measurements of radon emanation and krypton concentrations where we find good agreement. We report an observed reduction in concentrations of radon daughters that we attribute to the plating-out of charged ions on the negatively biased cathode.

Astrometric and spectroscopic monitoring of individual stars orbiting the supermassive black hole in the Galactic Center offer a promising way to detect general relativistic effects. While low-order effects are expected to be detected following the periastron passage of S2 in Spring 2018, detecting higher-order effects due to black hole spin will require the discovery of closer stars. In this paper, we set out to determine the requirements such a star would have to satisfy to allow the detection of black hole spin. We focus on the instrument GRAVITY, which saw first light in 2016 and which is expected to achieve astrometric accuracies $10-100 \mu$as. For an observing campaign with duration $T$ years, $N_{obs}$ total observations, astrometric precision $\sigma_x$ and normalized black hole spin $\chi$, we find that $a_{orb}(1-e^2)^{3/4} \lesssim 300 R_S \sqrt{\frac{T}{4 \text{years}}} \left(\frac{N_{obs}}{120}\right)^{0.25} \sqrt{\frac{10 \mu as}{\sigma_x}} \sqrt{\frac{\chi}{0.9}}$ is needed. For $\chi=0.9$ and a potential observing campaign with $\sigma_x = 10 \mu$as, 30 observations/year and duration 4-10 years, we expect $\sim 0.1$ star with $K<19$ satisfying this constraint based on the current knowledge about the stellar population in the central 1". We also propose a method through which GRAVITY could potentially measure radial velocities with precision $\sim 50$ km/s. If the astrometric precision can be maintained, adding radial velocity information increases the expected number of stars by roughly a factor of two. While we focus on GRAVITY, the results can also be scaled to parameters relevant for future extremely large telescopes.

We report on the response of liquid xenon to low energy electronic recoils below 15 keV from beta decays of tritium at drift fields of 92 V/cm, 154 V/cm and 366 V/cm using the XENON100 detector. A data-to-simulation fitting method based on Markov Chain Monte Carlo is used to extract the photon yields and recombination fluctuations from the experimental data. The photon yields measured at the two lower fields are in agreement with those from literature; additional measurements at a higher field of 366 V/cm are presented. The electronic and nuclear recoil discrimination as well as its dependence on the drift field and photon detection efficiency are investigated at these low energies. The results provide new measurements in the energy region of interest for dark matter searches using liquid xenon.

We present high resolution (sub-mas) VLBI maps of nuclear H2O megamasers for seven galaxies. In UGC6093, the well-aligned systemic masers and high-velocity masers originate in an edge-on, flat disk and we determine the mass of the central SMBH to be M_SMBH = 2.58*10^7Msun(+-7%). For J1346+5228, the distribution of masers is consistent with a disk, but the faint high-velocity masers are only marginally detected, and we constrain the mass of the SMBH to be in the range 1.5-2.0*10^7Msun. The origin of the masers in Mrk1210 is less clear, as the systemic and high-velocity masers are misaligned and show a disorganized velocity structure. We present one possible model in which the masers originate in a tilted, warped disk, but we do not rule out the possibility of other explanations including outflow masers. In NGC6926, we detect a set of redshifted masers, clustered within a pc of each other, and a single blueshifted maser about 4.4pc away, an offset that would be unusually large for a maser disk system. Nevertheless, if it is a disk system, we estimate the enclosed mass to be M_SMBH<4.8*10^7 Msun . For NGC5793, we detect redshifted masers spaced about 1.4pc from a clustered set of blueshifted features. The orientation of the structure supports a disk scenario as suggested by Hagiwara et al.(2001). We estimate the enclosed mass to be M SMBH<1.3*10^7 Msun. For NGC2824 and J0350-0127, the masers may be associated with pc or sub-pc scale jets or outflows.

We report the first dark matter search results from XENON1T, a $\sim$2000-kg-target-mass dual-phase (liquid-gas) xenon time projection chamber in operation at the Laboratori Nazionali del Gran Sasso in Italy and the first ton-scale detector of this kind. The blinded search used 34.2 live days of data acquired between November 2016 and January 2017. Inside the (1042$\pm$12) kg fiducial mass and in the [5, 40] $\mathrm{keV}_{\mathrm{nr}}$ energy range of interest for WIMP dark matter searches, the electronic recoil background was $(1.93 \pm 0.25) \times 10^{-4}$ events/(kg $\times$ day $\times \mathrm{keV}_{\mathrm{ee}}$), the lowest ever achieved in a dark matter detector. A profile likelihood analysis shows that the data is consistent with the background-only hypothesis. We derive the most stringent exclusion limits on the spin-independent WIMP-nucleon interaction cross section for WIMP masses above 10 GeV/c${}^2$, with a minimum of 7.7 $\times 10^{-47}$ cm${}^2$ for 35-GeV/c${}^2$ WIMPs at 90% confidence level.

We present the first search for dark matter-induced delayed coincidence signals in a dual-phase xenon time projection chamber, using the 224.6 live days of the XENON100 science run II. This very distinct signature is predicted in the framework of magnetic inelastic dark matter which has been proposed to reconcile the modulation signal reported by the DAMA/LIBRA collaboration with the null results from other direct detection experiments. No candidate event has been found in the region of interest and upper limits on the WIMP's magnetic dipole moment are derived. The scenarios proposed to explain the DAMA/LIBRA modulation signal by magnetic inelastic dark matter interactions of WIMPs with masses of 58.0 GeV/c$^2$ and 122.7 GeV/c$^2$ are excluded at 3.3 $\sigma$ and 9.3 $\sigma$, respectively.

We present results of searches for vector and pseudo-scalar bosonic super-WIMPs, which are dark matter candidates with masses at the keV-scale, with the XENON100 experiment. XENON100 is a dual-phase xenon time projection chamber operated at the Laboratori Nazionali del Gran Sasso. A profile likelihood analysis of data with an exposure of 224.6 live days $\times$ 34\,kg showed no evidence for a signal above the expected background. We thus obtain new and stringent upper limits in the $(8-125)$\,keV/c$^2$ mass range, excluding couplings to electrons with coupling constants of $g_{ae} > 3\times10^{-13}$ for pseudo-scalar and $\alpha'/\alpha > 2\times10^{-28}$ for vector super-WIMPs, respectively. These limits are derived under the assumption that super-WIMPs constitute all of the dark matter in our galaxy.

The XENON1T experiment at the Laboratori Nazionali del Gran Sasso (LNGS) is the first WIMP dark matter detector operating with a liquid xenon target mass above the ton-scale. Out of its 3.2t liquid xenon inventory, 2.0t constitute the active target of the dual-phase time projection chamber. The scintillation and ionization signals from particle interactions are detected with low-background photomultipliers. This article describes the XENON1T instrument and its subsystems as well as strategies to achieve an unprecedented low background level. First results on the detector response and the performance of the subsystems are also presented.

We report on WIMP search results in the XENON100 detector using a non-relativistic effective field theory approach. The data from science run II (34 kg $\times$ 224.6 live days) was re-analyzed, with an increased recoil energy interval compared to previous analyses, ranging from $(6.6 - 240)~\mathrm{keV_\mathrm{nr}}$. The data is found to be compatible with the background-only hypothesis. We present 90% confidence level exclusion limits on the coupling constants of WIMP-nucleon effective operators using a binned profile likelihood method. We also consider the case of inelastic WIMP scattering, where incident WIMPs may up-scatter to a higher mass state, and set exclusion limits on this model as well.

We present the first constraints on the spin-dependent, inelastic scattering cross section of Weakly Interacting Massive Particles (WIMPs) on nucleons from XENON100 data with an exposure of 7.64$\times$10$^3$\,kg\,day. XENON100 is a dual-phase xenon time projection chamber with 62\,kg of active mass, operated at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy and designed to search for nuclear recoils from WIMP-nucleus interactions. Here we explore inelastic scattering, where a transition to a low-lying excited nuclear state of $^{129}$Xe is induced. The experimental signature is a nuclear recoil observed together with the prompt de-excitation photon. We see no evidence for such inelastic WIMP-$^{129}$Xe interactions. A profile likelihood analysis allows us to set a 90\% C.L. upper limit on the inelastic, spin-dependent WIMP-nucleon cross section of $3.3 \times 10^{-38}$\,cm$^{2}$ at 100\,GeV/c$^2$. This is the most constraining result to date, and sets the pathway for an analysis of this interaction channel in upcoming, larger dual-phase xenon detectors.

We report on WIMP search results of the XENON100 experiment, combining three runs summing up to 477 live days from January 2010 to January 2014. Data from the first two runs were already published. A blind analysis was applied to the last run recorded between April 2013 and January 2014 prior to combining the results. The ultra-low electromagnetic background of the experiment, ~$5 \times 10^{-3}$ events/(keV$_{\mathrm{ee}}\times$kg$\times$day) before electronic recoil rejection, together with the increased exposure of 48 kg $\times$ yr improves the sensitivity. A profile likelihood analysis using an energy range of (6.6 - 43.3) keV$_{\mathrm{nr}}$ sets a limit on the elastic, spin-independent WIMP-nucleon scattering cross section for WIMP masses above 8 GeV/$c^2$, with a minimum of 1.1 $\times 10^{-45}$ cm$^2$ at 50 GeV/$c^2$ and 90% confidence level. We also report updated constraints on the elastic, spin-dependent WIMP-nucleon cross sections obtained with the same data. We set upper limits on the WIMP-neutron (proton) cross section with a minimum of 2.0 $\times 10^{-40}$ cm$^2$ (52$\times 10^{-40}$ cm$^2$) at a WIMP mass of 50 GeV/$c^2$, at 90% confidence level.

We report on a search for electronic recoil event rate modulation signatures in the XENON100 data accumulated over a period of 4 years, from January 2010 to January 2014. A profile likelihood method, which incorporates the stability of the XENON100 detector and the known electronic recoil background model, is used to quantify the significance of periodicity in the time distribution of events. There is a weak modulation signature at a period of $431^{+16}_{-14}$ days in the low energy region of $(2.0-5.8)$ keV in the single scatter event sample, with a global significance of $1.9\,\sigma$, however no other more significant modulation is observed. The expected annual modulation of a dark matter signal is not compatible with this result. Single scatter events in the low energy region are thus used to exclude the DAMA/LIBRA annual modulation as being due to dark matter electron interactions via axial vector coupling at $5.7\,\sigma$.

We present the first results of searches for axions and axion-like-particles with the XENON100 experiment. The axion-electron coupling constant, $g_{Ae}$, has been probed by exploiting the axio-electric effect in liquid xenon. A profile likelihood analysis of 224.6 live days $\times$ 34 kg exposure has shown no evidence for a signal. By rejecting $g_{Ae}$, larger than $7.7 \times 10^{-12}$ (90\% CL) in the solar axion search, we set the best limit to date on this coupling. In the frame of the DFSZ and KSVZ models, we exclude QCD axions heavier than 0.3 eV/c$^2$ and 80 eV/c$^2$, respectively. For axion-like-particles, under the assumption that they constitute the whole abundance of dark matter in our galaxy, we constrain $g_{Ae}$, to be lower than $1 \times 10^{-12}$ (90\% CL) for mass range from 1 to 40 keV/c$^2$, and set the best limit to date as well.

As part of the Megamaser Cosmology Project (MCP), we present VLBI maps of nuclear water masers toward five galaxies. The masers originate in sub-parsec circumnuclear disks. For three of the galaxies, we fit Keplerian rotation curves to estimate their supermassive black hole (SMBH) masses, and determine (2.9 $\pm$ 0.3) $\times~10^{6}M_\odot$ for J0437+2456, (1.7 $\pm$ 0.1) $\times~10^{7}M_\odot$ for ESO 558$-$G009, and (1.1 $\pm$ 0.2) $\times~10^{7}M_\odot$ for NGC 5495. In the other two galaxies, Mrk 1029 and NGC 1320, the geometry and dynamics are more complicated and preclude robust black hole mass estimates. Including our new results, we compiled a list of 15 VLBI-confirmed disk maser galaxies with robust SMBH mass measurements. With this sample, we confirm the empirical relation of $R_{out} \propto 0.3 M_{SMBH}$ reported in Wardle & Yusef-Zadeh (2012). We also find a tentative correlation between maser disk outer radii and WISE luminosity. We find no correlations of maser disk size with X-ray 2-10 keV luminosity or [O III] luminosity.

We use new precision measurements of black hole masses from water megamaser disks to investigate scaling relations between macroscopic galaxy properties and supermassive black hole (BH) mass. The megamaser-derived BH masses span 10^6-10^8 M_sun, while all the galaxy properties that we examine (including stellar mass, central mass density, central velocity dispersion) lie within a narrow range. Thus, no galaxy property correlates tightly with M_BH in ~L* spiral galaxies. Of them all, stellar velocity dispersion provides the tightest relation, but at fixed sigma* the mean megamaser M_BH are offset by -0.6+/-0.1 dex relative to early-type galaxies. Spiral galaxies with non-maser dynamical BH masses do not show this offset. At low mass, we do not yet know the full distribution of BH mass at fixed galaxy property; the non-maser dynamical measurements may miss the low-mass end of the BH distribution due to inability to resolve the spheres of influence and/or megamasers may preferentially occur in lower-mass BHs.

As part of the Megamaser Cosmology Project (MCP), here we present a new geometric distance measurement to the megamaser galaxy NGC 5765b. Through a series of VLBI observations, we have confirmed the water masers trace a thin, sub-parsec Keplerian disk around the nucleus, implying an enclosed mass of 4.55 $\pm$ 0.40 $\times~10^{7}M_\odot$. Meanwhile, from single dish monitoring of the maser spectra over two years, we measured the secular drifts of maser features near the systemic velocity of the galaxy with rates between 0.5 and 1.2 km s$^{-1}$ yr$^{-1}$. Fitting a warped, thin disk model to these measurements, we determine a Hubble Constant $H_{0}$ of 66.0 $\pm$ 6.0 km s$^{-1}$ Mpc$^{-1}$ with the angular-diameter distance to NGC 5765b of 126.3 $\pm$ 11.6 Mpc. Apart from the distance measurement, we also investigate some physical properties related to the maser disk in NGC 5765b. The high-velocity features are spatially distributed into several clumps, which may indicate the existence of a spiral density wave associated with the accretion disk. For the red-shifted features, the envelope defined by the peak maser intensities increases with radius. The profile of the systemic masers in NGC 5765b is smooth and shows almost no structural changes over the two years of monitoring time, which differs from the more variable case of NGC 4258.

Relic gravitational waves (RGWs) , a background originated during inflation, would give imprints on the pulsar timing residuals. This makes RGWs be one of important sources for detection using the method of pulsar timing. In this paper, we discuss the effects of RGWs on the single pulsar timing, and give quantitively the timing residuals caused by RGWs with different model parameters. In principle, if the RGWs are strong enough today, they can be detected by timing a single millisecond pulsar with high precision after the intrinsic red noise in pulsar timing residuals were understood, even though observing simultaneously multiple millisecond pulsars is a more powerful technique in extracting gravitational wave signals. We corrected the normalization of RGWs using observations of the cosmic microwave background (CMB), which leads to the amplitudes of RGWs being reduced by two orders of magnitude or so compared to our previous works. We made new constraints on RGWs using the recent observations from the Parkes Pulsar Timing Array, employing the tensor-to-scalar ratio $r=0.2$ due to the tensor-type polarization observations of CMB by BICEP2 as a referenced value even though it has been denied. Moreover, the constraints on RGWs from CMB and BBN (Big Bang nucleosynthesis) will also be discussed for comparison.

Laboratory experiments searching for galactic dark matter particles scattering off nuclei have so far not been able to establish a discovery. We use data from the XENON100 experiment to search for dark matter interacting with electrons. With no evidence for a signal above the low background of our experiment, we exclude a variety of representative dark matter models that would induce electronic recoils. For axial-vector couplings to electrons, we exclude cross-sections above 6x10^(-35) cm^2 for particle masses of m_chi = 2 GeV/c^2. Independent of the dark matter halo, we exclude leptophilic models as explanation for the long-standing DAMA/LIBRA signal, such as couplings to electrons through axial-vector interactions at a 4.4 sigma confidence level, mirror dark matter at 3.6 sigma, and luminous dark matter at 4.6 sigma.

We have searched for periodic variations of the electronic recoil event rate in the (2-6) keV energy range recorded between February 2011 and March 2012 with the XENON100 detector, adding up to 224.6 live days in total. Following a detailed study to establish the stability of the detector and its background contributions during this run, we performed an un-binned profile likelihood analysis to identify any periodicity up to 500 days. We find a global significance of less than 1 sigma for all periods suggesting no statistically significant modulation in the data. While the local significance for an annual modulation is 2.8 sigma, the analysis of a multiple-scatter control sample and the phase of the modulation disfavor a dark matter interpretation. The DAMA/LIBRA annual modulation interpreted as a dark matter signature with axial-vector coupling of WIMPs to electrons is excluded at 4.8 sigma.

We use single-dish radio spectra of known 22 GHz H$_2$O megamasers, primarily gathered from the large dataset observed by the Megamaser Cosmology Project, to identify Keplerian accretion disks and to investigate several aspects of the disk physics. We test a mechanism for maser excitation proposed by Maoz & McKee (1998), whereby population inversion arises in gas behind spiral shocks traveling through the disk. Though the flux of redshifted features is larger on average than that of blueshifted features, in support of the model, the high-velocity features show none of the predicted systematic velocity drifts. We find rapid intra-day variability in the maser spectrum of ESO 558-G009 that is likely the result of interstellar scintillation, for which we favor a nearby ($D \approx 70$ pc) scattering screen. In a search for reverberation in six well-sampled sources, we find that any radially-propagating signal must be contributing $\lesssim$10% of the total variability. We also set limits on the magnetic field strengths in seven sources, using strong flaring events to check for the presence of Zeeman splitting. These limits are typically 200--300 mG ($1\sigma$), but our most stringent limits reach down to 73 mG for the galaxy NGC 1194.

The XENON100 dark matter experiment uses liquid xenon in a time projection chamber (TPC) to measure xenon nuclear recoils resulting from the scattering of dark matter Weakly Interacting Massive Particles (WIMPs). In this paper, we report the observation of single-electron charge signals which are not related to WIMP interactions. These signals, which show the excellent sensitivity of the detector to small charge signals, are explained as being due to the photoionization of impurities in the liquid xenon and of the metal components inside the TPC. They are used as a unique calibration source to characterize the detector. We explain how we can infer crucial parameters for the XENON100 experiment: the secondary-scintillation gain, the extraction yield from the liquid to the gas phase and the electron drift velocity.

The XENON100 experiment, situated in the Laboratori Nazionali del Gran Sasso, aims at the direct detection of dark matter in the form of weakly interacting massive particles (WIMPs), based on their interactions with xenon nuclei in an ultra low background dual-phase time projection chamber. This paper describes the general methods developed for the analysis of the XENON100 data. These methods have been used in the 100.9 and 224.6 live days science runs from which results on spin-independent elastic, spin-dependent elastic and inelastic WIMP-nucleon cross-sections have already been reported.

The XENON100 experiment, installed underground at the Laboratori Nazionali del Gran Sasso (LNGS), aims to directly detect dark matter in the form of Weakly Interacting Massive Particles (WIMPs) via their elastic scattering off xenon nuclei. This paper presents a study on the nuclear recoil background of the experiment, taking into account neutron backgrounds from ($\alpha$,n) and spontaneous fission reactions due to natural radioactivity in the detector and shield materials, as well as muon-induced neutrons. Based on Monte Carlo simulations and using measured radioactive contaminations of all detector components, we predict the nuclear recoil backgrounds for the WIMP search results published by the XENON100 experiment in 2011 and 2012, 0.11$^{+0.08}_{-0.04}$ events and 0.17$^{+0.12}_{-0.07}$ events, respectively, and conclude that they do not limit the sensitivity of the experiment.

Results from the nuclear recoil calibration of the XENON100 dark matter detector installed underground at the Laboratori Nazionali del Gran Sasso (LNGS), Italy are presented. Data from measurements with an external 241AmBe neutron source are compared with a detailed Monte Carlo simulation which is used to extract the energy dependent charge-yield Qy and relative scintillation efficiency Leff. A very good level of absolute spectral matching is achieved in both observable signal channels - scintillation S1 and ionization S2 - along with agreement in the 2-dimensional particle discrimination space. The results confirm the validity of the derived signal acceptance in earlier reported dark matter searches of the XENON100 experiment.

We report on a search for particle dark matter with the XENON100 experiment, operated at the Laboratori Nazionali del Gran Sasso (LNGS) for 13 months during 2011 and 2012. XENON100 features an ultra-low electromagnetic background of (5.3 \pm 0.6) \times 10^-3 events (kg day keVee)^-1 in the energy region of interest. A blind analysis of 224.6 live days \times 34 kg exposure has yielded no evidence for dark matter interactions. The two candidate events observed in the pre-defined nuclear recoil energy range of 6.6-30.5 keVnr are consistent with the background expectation of (1.0 \pm 0.2) events. A Profile Likelihood analysis using a 6.6-43.3 keVnr energy range sets the most stringent limit on the spin-independent elastic WIMP-nucleon scattering cross section for WIMP masses above 8 GeV/c^2, with a minimum of 2 \times 10^-45 cm^2 at 55 GeV/c^2 and 90% confidence level.

We present new experimental constraints on the elastic, spin-dependent WIMP-nucleon cross section using recent data from the XENON100 experiment, operated in the Laboratori Nazionali del Gran Sasso in Italy. An analysis of 224.6 live days x 34 kg of exposure acquired during 2011 and 2012 revealed no excess signal due to axial-vector WIMP interactions with 129-Xe and 131-Xe nuclei. This leads to the most stringent upper limits on WIMP-neutron cross sections for WIMP masses above 6 GeV, with a minimum cross section of 3.5 x 10^{-40} cm^2 at a WIMP mass of 45 GeV, at 90% confidence level.

The XENON100 experiment, in operation at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, was designed to search for evidence of dark matter interactions inside a volume of liquid xenon using a dual-phase time projection chamber. This paper describes the Slow Control System (SCS) of the experiment with emphasis on the distributed architecture as well as on its modular and expandable nature. The system software was designed according to the rules of Object-Oriented Programming and coded in Java, thus promoting code reusability and maximum flexibility during commissioning of the experiment. The SCS has been continuously monitoring the XENON100 detector since mid 2008, remotely recording hundreds of parameters on a few dozen instruments in real time, and setting emergency alarms for the most important variables.

In a recent manuscript (arXiv:1208.5046) Peter Sorensen claims that XENON100's upper limits on spin-independent WIMP-nucleon cross sections for WIMP masses below 10 GeV "may be understated by one order of magnitude or more". Having performed a similar, though more detailed analysis prior to the submission of our new result (arXiv:1207.5988), we do not confirm these findings. We point out the rationale for not considering the described effect in our final analysis and list several potential problems with his study.

We use a model of an accretion flow coupled with an emergent flare to interpret the latest 1.3mm VLBI measurements for Sagittarius A*. The visibility data constrained the distances from the flare center to the black hole center as $d_{\rm EW}\lesssim20{\rm R_g}$ and $d_{\rm NS}\lesssim80{\rm R_g}$ in the East-West and North-South directions, respectively. If interpreted by the hot-spot model, the flare was preferred to pass in front of the black hole at a radius much larger than $d_{\rm EW}$. If interpreted by the episodic jet launched from a nearly edge-on hot accretion flow, the flare was preferred to be ejected with $\theta_{\rm j}\gtrsim40^\circ$ off the black hole rotating axis. This method can be generalized to help us understand future sub-millimeter VLBI observations, and study the millimeter/sub-millimeter variabilities in the vicinity of the Galactic Center supermassive black hole.

The XENON100 experiment has recently completed a dark matter run with 100.9 live-days of data, taken from January to June 2010. Events in a 48kg fiducial volume in the energy range between 8.4 and 44.6 keVnr have been analyzed. A total of three events have been found in the predefined signal region, compatible with the background prediction of (1.8 \pm 0.6) events. Based on this analysis we present limits on the WIMP-nucleon cross section for inelastic dark matter. With the present data we are able to rule out the explanation for the observed DAMA/LIBRA modulation as being due to inelastic dark matter scattering off iodine at a 90% confidence level.

We present results from the direct search for dark matter with the XENON100 detector, installed underground at the Laboratori Nazionali del Gran Sasso of INFN, Italy. XENON100 is a two-phase time projection chamber with a 62 kg liquid xenon target. Interaction vertex reconstruction in three dimensions with millimeter precision allows to select only the innermost 48 kg as ultra-low background fiducial target. In 100.9 live days of data, acquired between January and June 2010, no evidence for dark matter is found. Three candidate events were observed in a pre-defined signal region with an expected background of 1.8 +/- 0.6 events. This leads to the most stringent limit on dark matter interactions today, excluding spin-independent elastic WIMP-nucleon scattering cross-sections above 7.0x10^-45 cm^2 for a WIMP mass of 50 GeV/c^2 at 90% confidence level.