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.

Using data of the Baksan Underground Scintillation Telescope (BUST) we have made a search for muon neutrinos and antineutrinos with energies above 1 GeV coinciding with the gravitational wave event GW170817 that was recorded on August 17, 2017 by the Advanced LIGO and Advanced Virgo observatories. This is a first detection of the new type of events occurring as a result of a merger of two neutron stars in a binary system. A short gamma-ray burst GRB170817A accompanying this event is an evidence of particle acceleration in the source whose precise position was determined by detection of the subsequent optical signal. No neutrino signals were found with the BUST in the interval $\pm 500$ s around the moment of the gravitational wave event GW170817, as well as during the next 14 days. The upper limits on integral fluxes of muon neutrino and antineutrino from the source are derived.

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 highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Astrometric detection of an unseen companion is based on analysis of apparent motion of its host star around the system's barycentre. Systems with orbital period close to one year may escape detection if orbital motion of their host stars are observationally indistinguishable from parallax effect. Additionally, the astrometric solution may produce a biased parallax estimation for such systems. We examine effects of orbital motion of the Earth on astrometric detectability in terms of correlation between the Earth's orbital position and position of the star relative to its system barycentre. The $\chi^2$ statistic for parallax estimation is calculated analytically, leading to expressions that relate the decrease in detectability and accompanying parallax bias to the position correlation function. The impact of the Earth's motion critically depends on the exoplanet's orbital period, diminishing rapidly as the period deviates from one year. Selection effects against one-year period systems is therefore expected. Statistical estimation shows that the corresponding loss of sensitivity results in a typical 10 per cent increase in detection threshold. Consideration of eccentric orbits shows that the Earth's motion has no effect on detectability for $e\ge0.5$. Dependence of detectability on other parameters, such as orbital phases and inclination of the orbit plane to the ecliptic, are smooth and monotonic because they are described by simple trigonometric functions.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Determination of absolute parallaxes by means of a scanning astrometric satellite such as Hipparcos or Gaia relies on the short-term stability of the so-called basic angle between the two viewing directions. Uncalibrated variations of the basic angle may produce systematic errors in the computed parallaxes. We examine the coupling between a global parallax shift and specific variations of the basic angle, namely those related to the satellite attitude with respect to the Sun. The changes in observables produced by small perturbations of the basic angle, attitude, and parallaxes are calculated analytically. We then look for a combination of perturbations that has no net effect on the observables. In the approximation of infinitely small fields of view, it is shown that certain perturbations of the basic angle are observationally indistinguishable from a global shift of the parallaxes. If such perturbations exist, they cannot be calibrated from the astrometric observations but will produce a global parallax bias. Numerical simulations of the astrometric solution, using both direct and iterative methods, confirm this theoretical result. For a given amplitude of the basic angle perturbation, the parallax bias is smaller for a larger basic angle and a larger solar aspect angle. In both these respects Gaia has a more favourable geometry than Hipparcos. In the case of Gaia, internal metrology is used to monitor basic angle variations. Additionally, Gaia has the advantage of detecting numerous quasars, which can be used to verify the parallax zero point.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

As part of the data processing for Gaia Data Release~1 (Gaia DR1) a special astrometric solution was computed, the so-called auxiliary quasar solution. This gives positions for selected extragalactic objects, including radio sources in the second realisation of the International Celestial Reference Frame (ICRF2) that have optical counterparts bright enough to be observed with Gaia. A subset of these positions was used to align the positional reference frame of Gaia DR1 with the ICRF2. We describe the properties of the Gaia auxiliary quasar solution for a subset of sources matched to ICRF2, and compare their optical and radio positions at the sub-mas level. Their formal standard errors are better than 0.76~milliarcsec (mas) for 50% of the sources and better than 3.35~mas for 90%. Optical magnitudes are obtained in Gaia's unfiltered photometric G band. The comparison with the radio positions of the defining sources shows no systematic differences larger than a few tenths of a mas. The fraction of questionable solutions, not readily accounted for by the statistics, is less than 6%. Normalised differences have extended tails requiring case-by-case investigations for around 100 sources, but we have not seen any difference indisputably linked to an optical-radio offset in the sources.

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame.

Gaia's astrometric solution aims to determine at least five parameters for each star, together with appropriate estimates of their uncertainties and correlations. This requires at least five distinct observations per star. In the early data reductions the number of observations may be insufficient for a five-parameter solution, and even after the full mission many stars will remain under-observed, including faint stars at the detection limit and transient objects. In such cases it is reasonable to determine only the two position parameters. Their formal uncertainties would however grossly underestimate the actual errors, due to the neglected parallax and proper motion. We aim to develop a recipe to calculate sensible formal uncertainties that can be used in all cases of under-observed stars. Prior information about the typical ranges of stellar parallaxes and proper motions is incorporated in the astrometric solution by means of Bayes' rule. Numerical simulations based on the Gaia Universe Model Snapshot (GUMS) are used to investigate how the prior influences the actual errors and formal uncertainties when different amounts of Gaia observations are available. We develop a criterion for the optimum choice of priors, apply it to a wide range of cases, and derive a global approximation of the optimum prior as a function of magnitude and galactic coordinates. The feasibility of the Bayesian approach is demonstrated through global astrometric solutions of simulated Gaia observations. With an appropriate prior it is possible to derive sensible positions with realistic error estimates for any number of available observations. Even though this recipe works also for well-observed stars it should not be used where a good five-parameter astrometric solution can be obtained without a prior. Parallaxes and proper motions from a solution using priors are always biased and should not be used.

High-precision astrometric and radial-velocity observations require accurate modelling of stellar motions in order to extrapolate measurements over long time intervals, and to detect deviations from uniform motion caused for example by unseen companions. We aim to explore the simplest possible kinematic model of stellar motions, namely that of uniform rectilinear motion relative to the Solar System Barycentre, in terms of observable quantities including error propagation. The apparent path equation for uniform rectilinear motion is solved analytically in a classical (special-relativistic) framework, leading to rigorous expressions which relate the (apparent) astrometric parameters and radial velocity to the (true) kinematic parameters of the star in the barycentric reference system. We present rigorous and explicit formulae for the transformation of stellar positions, parallaxes, proper motions, and radial velocities from one epoch to another, assuming uniform rectilinear motion and taking into account light-time effects. The Jacobian matrix of the transformation is also given, allowing accurate and reversible propagation of errors over arbitrary time intervals. The light-time effects are generally very small but exceeds 0.1 mas or 0.1 m/s over 100 yr for at least 33 stars in the Hipparcos Catalogue. For high-velocity stars within a few tens of pc from the Sun light-time effects are generally more important than the effects of the curvature of their orbits in the Galactic potential.

This paper summarizes the analysis of the consequences of the violation of the Local Lorentz Invariance (LLI) on astrometric observations. We demonstrate that from the point of view of the LLI astrometric observations represent an experiment of Michelson-Morley type. The future high-accuracy astrometric projects (e.g., Gaia) will be used to test the LLI.

The problem of determination of the orbital velocity of an astrometric satellite from its own observational data is studied. It is well known that data processing of microarcsecond-level astrometric observations imposes very stringent requirements on the accuracy of the orbital velocity of the satellite (a velocity correction of 1.45 mm/s implies an aberrational correction of 1 $\mu$as). Because of a number of degeneracies the orbital velocity cannot be fully restored from observations provided by the satellite. Seven constraints that must be applied on a velocity parameterization are discussed and formulated mathematically. It is shown what part of velocity can be recovered from astrometric data by a combined fit of both velocity parameters and astrometric parameters of the sources. Numerical simulations show that, with the seven constraints applied, the velocity and astrometric parameters can be reliably estimated from observational data. It is also argued that the idea to improve the velocity of an astrometric satellite from its own observational data is only useful if the a priori information on the orbital velocity justifies the applicability of the velocity constraints. The proposed model takes into account only translational motion of the satellite and ignores any satellite-specific parameters. Therefore, the results of this study are equally applicable to both scanning missions similar to Gaia, and pointing ones like SIM, provided that enough sources were observed sufficiently uniformly.

We derive the luminosity-temperature relation for the super-critically accreting black holes (BHs) and compare it to the data on ultraluminous X-ray sources (ULXs). At super-Eddington accretion rates, an outflow forms within the spherization radius. We construct the accretion disc model accounting for the advection and the outflow, and compute characteristic disc temperatures. The bolometric luminosity exceeds the Eddington luminosity L_Edd by a logarithmic factor 1+0.6 ln mdot (where mdot is the accretion rate in Eddington units) and the wind kinetic luminosity is close to L_Edd. The apparent luminosity for the face-on observer is 2-7 times higher because of geometrical beaming. Such an observer has a direct view of the inner hot accretion disc, which has a peak temperature T_max of a few keV in stellar-mass BHs. The emitted spectrum extends as a power-law F_E ~ E**{-1} down to the temperature at the spherization radius T_sp ~ mdot**(-1/2) keV. We associate T_max with a few keV spectral components and T_sp with the soft, 0.1-0.2 keV components observed in ULXs. An edge-on observer sees only the soft emission from the extended envelope, with the photosphere radius exceeding the spherization radius by orders of magnitude. The dependence of the photosphere temperature on luminosity is consistent with that observed in the super-Eddington accreting BHs SS 433 and V4641 Sgr. Strong outflows combined with the large intrinsic X-ray luminosity of the central BH explain naturally the presence of the photoionized nebulae around ULXs. An excellent agreement between the model and the observational data strongly argues in favour of ULXs being super-critically accreting, stellar-mass BHs similar to SS 433, but viewed close to the symmetry axis.