
Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in
the magnitude range 3 to 21 based on observations collected by the European
Space Agency Gaia satellite during the first 22 months of its operational
phase. We describe the input data, models, and processing used for the
astrometric content of Gaia DR2, and the validation of these results performed
within the astrometry task. Some 320 billion centroid positions from the
preprocessed astrometric CCD observations were used to estimate the five
astrometric parameters (positions, parallaxes, and proper motions) for 1332
million sources, and approximate positions at the reference epoch J2015.5 for
an additional 361 million mostly faint sources. Special validation solutions
were used to characterise the random and systematic errors in parallax and
proper motion. For the sources with fiveparameter astrometric solutions, the
median uncertainty in parallax and position at the reference epoch J2015.5 is
about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas
at G=20 mag. In the proper motion components the corresponding uncertainties
are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame
defined by Gaia DR2 is aligned with ICRS and is nonrotating with respect to
the quasars to within 0.15 mas/yr. From the quasars and validation solutions we
estimate that systematics in the parallaxes depending on position, magnitude,
and colour are generally below 0.1 mas, but the parallaxes are on the whole too
small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in
parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and
intermediate (20 deg) angular scales. Important statistics and information for
the users of the Gaia DR2 astrometry are given in the appendices.

We highlight the power of the Gaia DR2 in studying many fine structures of
the HertzsprungRussell 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, postAGB stars, planetary nebulae, and white dwarfs are also
analysed, as well as lowmass 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
HRDbased studies.

We revisit earlier calculations of leptonic decay constants of vector
charmonia and present and illustrate our decomposition of the corresponding
covariant BetheSalpeter amplitudes in terms of orbital angular momentum as
interpreted in the meson's rest frame. Our results confirm our previous
conclusions drawn from the magnitudes of vectormeson leptonic decay constants,
identifying the Psi(3770) and Psi(4160) as Dwave states in our setup.

We investigate the light quarkonium spectrum using a covariant
DysonSchwingerBetheSalpeterequation approach to QCD. We discuss splittings
among as well as orbital angular momentum properties of various states in
detail and analyze common features of mass splittings with regard to properties
of the effective interaction. In particular, we predict the mass of sbars
exotic 1+ states, and identify orbital angular momentum content in the
excitations of the rho meson. Comparing our covariant model results, the rho
and its second excitation being predominantly Swave, the first excitation
being predominantly Dwave, to corresponding conflicting latticeQCD studies,
we investigate the pionmass dependence of the orbitalangularmomentum
assignment and find a crossing at a scale of $m_\pi$ ~ 1.4 GeV. If this
crossing turns out to be a feature of the spectrum generated by latticeQCD
studies as well, it may reconcile the different results, since they have been
obtained at different values of $m_\pi$.

Openflavour meson studies are the necessary completion to any comprehensive
investigation of quarkonia. We extend recent studies of quarkonia in the
DysonSchwingerBetheSalpeterequation approach to explore their results for
all possible flavour combinations. Within the inherent limitations of the
setup, we present the most comprehensive results for meson masses and leptonic
decay constants currently available and put them in perspective with respect to
experiment and other approaches.

Meson states with exotic quantum numbers arise naturally in a covariant
boundstate framework in QCD. We investigate the consequences of shifting quark
masses such that the states are no longer restricted to certain Cparities, but
only by J^P. Then, a priori, one can no longer distinguish exotic or
conventional states. In order to identify signatures of the different states to
look for experimentally, we provide the behavior of masses, leptonic decay
constants, and orbitalangularmomentum decomposition of such mesons, as well
as the constellations in which they could be found. Most prominently, we
consider the case of charged quasiexotic excitations of the pion.

We give a short and basic introduction to our covariant
DysonSchwingerBetheSalpeterequation approach using a rainbowladder
truncated model of QCD, in which we investigate the leptonic decay properties
of heavy quarkonium states in the pseudoscalar and vector channels. Comparing
the magnitudes of decay constants, we identify radial 1 excitations in our
calculation with experimental excitations of J/\Psi and \Upsilon. Particular
attention is paid to those states regarded as Dwave states in the quark model.
We predict e+e decay width of the \Upsilon(1^3D_1) and \Upsilon(2^3D_1) states
of the order of ca. 15 eV or more. We also provide a set of predictions for
decay constants of pseudoscalar radial excitations in heavy quarkonia.

We extend earlier investigations of heavylight pseudoscalar mesons to the
vector case, using a simple model in the context of the
DysonSchwingerBetheSalpeter approach. We investigate the effects of a
dressedquarkgluon vertex in a systematic fashion and illustrate and attempt
to quantify corrections beyond the phenomenologically very useful and
successful rainbowladder truncation. In particular we investigate dressed
quark photon vertex in such a setup and make a prediction for the
experimentally as yet unknown mass of the B_c*, which we obtain at 6.334 GeV
well in line with predictions from other approaches. Furthermore, we combine a
comprehensive set of results from the theory literature. The theory average for
the mass of the B_c* meson is 6.336 + 0.002 GeV.

Using a simple model in the context of the DysonSchwingerBetheSalpeter
approach, we investigate the effects of a dressedquarkgluon vertex on
pseudoscalar meson masses. In particular, we focus on the unequalmass case and
investigate heavylight meson masses; in addition, we study the premise of the
effective treatment of heavy quarks in our approach.

We investigate the properties of mesons with the exotic J^PC = 1^+ quantum
numbers. Starting out from the lightquark domain, where the \pi_1 states are
used as references, we predict the masses of analogous quarkonia for c\bar{c}
and b\bar{b} configurations. We employ a covariant
DysonSchwingerBetheSalpeterequation approach with a rainbowladder
truncated model of quantum chromodynamics.

QCD sum rules for $D$ and $B$ mesons in a strongly interacting environment
require the decomposition of Gibbs averaged operators related to tensor
structures. We present a procedure to decompose these operators into vacuum and
mediumspecific parts, thus defining plain mediumspecific condensates with
coefficients vanishing in vacuum. Our decomposition allows for an
identification of potential elements of order parameters for chiral
restoration, in particular for higher massdimension quarkcondensates which
would be masked otherwise if operators with uncontracted Lorentz indices would
be linked to DIS amplitudes.

The mass spectrum of heavy pseudoscalar mesons, described as quarkantiquark
bound systems, is considered within the BetheSalpeter formalism with momentum
dependent masses of the constituents. This dependence is found by solving the
SchwingerDyson equation for quark propagators in rainbowladder approximation.
Such an approximation is known to provide both a fast convergence of numerical
methods and accurate results for lightest mesons. However, as the meson mass
increases, the method becomes less stable and special attention must be devoted
to details of numerical means of solving the corresponding equations. We focus
on the pseudoscalar sector and show that our numerical scheme describes fairly
accurately the $\pi$, $K$, $D$, $D_s$ and $\eta_c$ ground states. Excited
states are considered as well. Our calculations are directly related to the
future physics programme at FAIR.

The mass spectrum of heavy pseudoscalar mesons, described as quarkantiquark
bound systems, is considered within the BetheSalpeter formalism with
momentumdependent masses of the constituents. This dependence is found by
solving the SchwingerDyson equation for quark propagators in rainbowladder
approximation. Such an approximation is known to provide both a fast
convergence of numerical methods and accurate results for lightest mesons.
However, as the meson mass increases, the method becomes less stable and
special attention must be devoted to details of numerical means of solving the
corresponding equations. We focus on the pseudoscalar sector and show that our
numerical scheme describes fairly accurately the $\pi$, $K$, $D$, $D_s$ and
$\eta_c$ ground states. Excited states are considered as well. Our calculations
are directly related to the future physics programme at FAIR.

We survey some prospects of identifying further inmedium modifications of
hadrons in a strongly interacting medium with respect to ongoing experiment
series of the HADES Collaboration and planned experiments of the CBM
Collaboration at FAIR. Dielectrons, strange and charm mesons are considered
and their potential for signaling imprints of chiral restoration is
highlighted.

Open charm mesons (pseudoscalar and scalar as well as axialvector and
vector) propagating or resting in nuclear matter display an enhanced
sensitivity to the chiral condensate. This offers new prospects to seek for
signals of chiral restoration, in particular in pA and pbarA reactions as
envisaged in firstround experiments by the CBM and PANDA collaborations at
FAIR. Weinberg type sum rules for charming chiral partners are presented, and
the distinct inmedium modifications of opencharm mesons are discussed. We
also address the gluon condensates near Tc and their impact on QCD sum rules.

Pseudoscalar and scalar $D$ mesons are considered within the QCD sum rule
approach. We present an analysis of the mass splitting of the pseudoscalar D 
Dbar mesons and the relation to QCD condensates. Weinberg type sum rules are
derived for chiral partners which highlight the role of the chiral condensate.

QCD sum rules for D and B mesons embedded in cold nuclear matter are
evaluated. We quantify the mass splitting of D  Dbar and B  Bbar mesons as
a function of the nuclear matter density; extrapolated to saturation density it
is in the order of 60 and 130 MeV driven essentially by the condensates
<qadjoint q>, <qadjoint g sigma G q> and <qbar q>. The genuine chiral
condensate <qbar q>, amplified by heavyquark masses, enters the Borel
transformed sum rules for the mass splitting beyond linear density dependence.
Including strange quark condensates reveals a umerically smaller and opposite
effect for the Ds  Dsbar mass splitting.

The QCD sum rule approach to the inmedium behavior of hadrons is discussed
for omega meson, nucleon and D meson. Emphasis is devoted to the impact of
fourquark condensates and to order parameters of spontaneous symmetry
breaking.

The inmedium behavior of the nucleon spectral density including
selfenergies is revisited within the framework of QCD sum rules. Special
emphasis is given to the density dependence of fourquark condensates. A
complete catalog of fourquark condensates is presented and relations among
them are derived. Generic differences of such fourquark condensates occurring
in QCD sum rules for light baryons and light vector mesons are discussed.

Spectral properties of hadrons in nuclear matter are treated in the framework
of QCD sum rules. The influence of the ambient strongly interacting medium is
encoded in various condensates. Especially, the structure of different
fourquark condensates and their density dependencies in light quark systems
are exemplified for the omega meson and the nucleon.

With QCD sum rule evaluations, spectral changes of hadrons inside nuclear
matter are considered, which shed light on QCD condensates and thus on the
nonperturbative structure of the QCD ground state. For some light quark
configurations, omega meson and nucleon, the relevance of fourquark
condensates is compared; the D meson as a representative of heavylight quark
systems is also briefly discussed.