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.

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.

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.

Aims. The Sloan Digital Sky Survey (SDSS) offers a unique possibility of not only detecting Colour Induced Displacement (CID) double stars but also confirming these detections. Methods. Successive cuts are applied to the SDSS DR12 database in order to reduce the size of the sample to be considered. The resulting dataset is then screened with a criterion based on the distance and orientation of the photocentres in different photometric bands. Results. About 3 200 distinct objects are classified as CID double stars, 40 of which are confirmed with at least a second detection. A consistency check further validates these detections.

Context. The Mass loss of Evolved StarS (MESS) sample observed with PACS on board the Herschel Space Observatory revealed that several asymptotic giant branch (AGB) stars are surrounded by an asymmetric circumstellar envelope (CSE) whose morphology is most likely caused by the interaction with a stellar companion. The evolution of AGB stars in binary systems plays a crucial role in understanding the formation of asymmetries in planetary nebul{\ae} (PNe), but at present, only a handful of cases are known where the interaction of a companion with the stellar AGB wind is observed. Aims. We probe the environment of the very evolved AGB star $\pi^1$ Gruis on large and small scales to identify the triggers of the observed asymmetries. Methods. Observations made with Herschel/PACS at 70 $\mu$m and 160 $\mu$m picture the large-scale environment of $\pi^1$ Gru. The close surroundings of the star are probed by interferometric observations from the VLTI/AMBER archive. An analysis of the proper motion data of Hipparcos and Tycho-2 together with the Hipparcos Intermediate Astrometric Data help identify the possible cause for the observed asymmetry. Results. The Herschel/PACS images of $\pi^1$ Gru show an elliptical CSE whose properties agree with those derived from a CO map published in the literature. In addition, an arc east of the star is visible at a distance of $38^{\prime\prime}$ from the primary. This arc is most likely part of an Archimedean spiral caused by an already known G0V companion that is orbiting the primary at a projected distance of 460 au with a period of more than 6200 yr. However, the presence of the elliptical CSE, proper motion variations, and geometric modelling of the VLTI/AMBER observations point towards a third component in the system, with an orbital period shorter than 10 yr, orbiting much closer to the primary than the G0V star.

Gaia is a very ambitious mission of the European Space Agency. At the heart of Gaia lie the measurements of the positions, distances, space motions, brightnesses and astrophysical parameters of stars, which represent fundamental pillars of modern astronomical knowledge. We provide a brief description of the Gaia mission with an emphasis on binary stars. In particular, we summarize results of simulations, which estimate the number of binary stars to be processed to several tens of millions. We also report on the catalogue release scenarios. In the current proposal, the first results for binary stars will be available in 2017 (for a launch in 2013).

Two upcoming large scale surveys, the ESA Gaia and LSST projects, will bring a new era in astronomy. The number of binary systems that will be observed and detected by these projects is enormous, estimations range from millions for Gaia to several tens of millions for LSST. We review some tools that should be developed and also what can be gained from these missions on the subject of binaries and exoplanets from the astrometry, photometry, radial velocity and their alert systems.

The Hyades stream has long been thought to be a dispersed vestige of the Hyades cluster. However, recent analyses of the parallax distribution, of the mass function, and of the action-space distribution of stream stars have shown it to be rather composed of orbits trapped at a resonance of a density disturbance. This resonant scenario should leave a clearly different signature in the element abundances of stream stars than the dispersed cluster scenario, since the Hyades cluster is chemically homogeneous. Here, we study the metalllicty as well as the element abundances of Li, Na, Mg, Fe, Zr, Ba, La, Ce, Nd, and Eu for a random sample of stars belonging to the Hyades stream, and compare them with those of stars from the Hyades cluster. From this analysis: (i) we independently confirm that the Hyades stream cannot be solely composed of stars originating in the Hyades cluster; (ii) we show that some stars from the Hyades stream nevertheless have abundances compatible with an origin in the cluster; (iii) we emphasize that the use of Li as a chemical tag of the cluster origin of main-sequence stars is very efficient in the range 5500 - 6200 K, since the Li sequence in the Hyades cluster is very tight; (iv) we show that, while this evaporated population has a metallicity excess of ~0.2 dex w.r.t. the local thin disk population, identical to that of the Hyades cluster, the remainder of the Hyades stream population has still a metallicity excess of ~0.06 to 0.15 dex, consistent with an origin in the inner Galaxy; (v) we show that the Hyades stream can be interpreted as an inner 4:1 resonance of the spiral pattern: this then also reproduces an orbital family compatible with the Sirius stream, and places the origin of the Hyades stream up to 1 kpc inwards from the solar radius, which might explain the observed metallicity excess of the stream population.

We obtained Spitzer/IRAC 3.6-8 micron images of the nearby spiral galaxy NGC 4258 to study possible interactions between dust and the radio jet. In our analysis we also included high-resolution radio continuum, H-alpha, CO, and X-ray data. Our data reveal that the 8 micron emission, believed to originate largely from PAH molecules and hot dust, is an excellent tracer of the normal spiral structure in NGC 4258, and hence it originates from the galactic plane. We investigated the possibility of dust destruction by the radio jet by calculating correlation coefficients between the 8 micron and radio continuum emissions along the jet in two independent ways, namely (i) from wavelet-transformed maps of the original images at different spatial scales, and (ii) from one-dimensional intensity cuts perpendicular to the projected path of the radio jet on the sky. No definitive sign of a correlation (or anticorrelation) was detected on relevant spatial scales with either approach, implying that any dust destruction must take place at spatial scales that are not resolved by our observations.

We report five new measurements of central black hole masses based on STIS and WFPC2 observations with the Hubble Space Telescope and on axisymmetric, three-integral, Schwarzschild orbit-library kinematic models. We selected a sample of galaxies within a narrow range in velocity dispersion that cover a range of galaxy parameters (including Hubble type and core/power-law surface density profile) where we expected to be able to resolve the galaxy's sphere of influence based on the predicted value of the black hole mass from the M-sigma relation. We find masses in units of 10^8 solar masses for the following galaxies: NGC 3585, M_BH = 3.4 (+1.5, -0.6); NGC 3607, M_BH = 1.2 (+0.4, -0.4); NGC 4026, M_BH = 2.1 (+0.7, -0.4); and NGC 5576, M_BH = 1.8 (+0.3, -0.4), all significantly excluding M_BH = 0. For NGC 3945, M_BH = 0.09 (+0.17, -0.21), which is significantly below predictions from M-sigma and M-L relations and consistent with M_BH = 0, though the presence of a double bar in this galaxy may present problems for our axisymmetric code.

We derive improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on 49 M_BH measurements and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. We find log(M_BH / M_sun) = alpha + beta * log(sigma / 200 km/s) with (alpha, beta, epsilon_0) = (8.12 +/- 0.08, 4.24 +/- 0.41, 0.44 +/- 0.06) for all galaxies and (alpha, beta, epsilon_0) = (8.23 +/- 0.08, 3.96 +/- 0.42, 0.31 +/- 0.06) for ellipticals. The results for ellipticals are consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter inferred reinforces the need for its consideration when calculating local black hole mass function based on the M-sigma relation, and further implies that there may be substantial selection bias in studies of the evolution of the M-sigma relation. We estimate the M-L relationship as log(M_BH / M_sun) = alpha + beta * log(L_V / 10^11 L_sun,V) of (alpha, beta, epsilon_0) = (8.95 +/- 0.11, 1.11 +/- 0.18, 0.38 +/- 0.09); using only early-type galaxies. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to larger mean masses, larger slopes, and incorrect intrinsic residuals.

We determine the mass of the black hole at the center of the spiral galaxy NGC 4258 by constructing axisymmetric dynamical models of the galaxy. These models are constrained by high spatial resolution imaging and long-slit spectroscopy of the nuclear region obtained with the {\em Hubble Space Telescope}, complemented by ground-based observations extending to larger radii. Our best mass estimate is $\MBH = (3.3 \pm 0.2) \times 10^7 \MSun $ for a distance of 7.28 Mpc (statistical errors only). This is within 15% of $ (3.82\pm 0.01) \times 10^7 \MSun$, the mass determined from the kinematics of water masers (rescaled to the same distance) assuming they are in Keplerian rotation in a warped disk. The construction of accurate dynamical models of NGC 4258 is somewhat compromised by an unresolved active nucleus and color gradients, the latter caused by variations in the stellar population and/or obscuring dust. These problems are not present in the $\sim 30$ other black hole mass determinations from stellar dynamics that have been published by us and other groups; thus, the relatively close agreement between the stellar dynamical mass and the maser mass in NGC 4258 enhances our confidence in the black hole masses determined in other galaxies from stellar dynamics using similar methods and data of comparable quality.

The largest galaxies, and in particular central galaxies in clusters, offer unique insight into understanding the mechanism for the growth of nuclear black holes. We present Hubble Space Telescope kinematics for NGC1399, the central galaxy in Fornax. We find the best-fit model contains a black hole of 5.1 +-0.7 x 10^8 Msun (at a distance of 21.1 Mpc), a factor of over 2 below the correlation of black hole mass and velocity dispersion. We also find a dramatic signature for central tangential anisotropy. The velocity profiles on adjacent sides 0.5" away from the nucleus show strong bimodality, and the central spectrum shows a large drop in the dispersion. Both of these observations point to an orbital distribution that is tangentially biased. The best-fit orbital model suggests a ratio of the tangential to radial internal velocity dispersions of three. This ratio is the largest seen in any galaxy to date and will provide an important measure for the mode by which the central black hole has grown.

We analyze HST+WFPC2 images of 77 early-type galaxies. Brightness profiles are classed into "core" or "power-law" forms. Cores are typically rounder than power-law galaxies. Nearly all power-laws with central ellipticity >=0.3 have stellar disks, implying that disks are present in power-laws with epsilon <0.3, but are not visible due to unfavorable geometry. A few low-luminosity core galaxies also have disks; these may be transition forms from power-laws. Cores and power-laws both have twisting isophotes at small radii. Core galaxies have somewhat weaker color gradients than power-laws. Nuclei are found in 29% of the cores and 60% of the power-laws. Nuclei are typically bluer than the surrounding galaxy. NGC 4073 and 4382 have central minima in their intrinsic starlight distributions; NGC 4382 resembles the double nucleus of M31. In general, the peak brightness location is coincident with the photocenter of the core to <1 pc. Five galaxies, however, have centers significantly displaced from their cores; these may be unresolved asymmetric double nuclei. Central dust is visible in half of the galaxies. The presence and strength of dust correlates with nuclear emission. The prevalence of dust and its morphology suggests that dust clouds form, settle to the center, and disappear repeatedly on ~10^8 yr timescales. We discuss the hypothesis that cores are created by the decay of a massive black hole binary. Apart from their brightness profiles, there are no strong differences between cores and power-laws that demand this scenario; however, the rounder shapes of cores, their lack of disks, and their reduced color gradients may be consistent with it.

We show that orbit-superposition dynamical models (Schwarzschild's method) provide reliable estimates of nuclear black hole masses and errors when constructed from adequate orbit libraries and kinematic data. We thus rebut two recent papers that argue that BH masses obtained from this method are unreliable. These papers claim to demonstrate that the range of allowable BH masses derived from a given dataset is artificially too narrow as a result of an inadequate number of orbits in the library used to construct dynamical models. This is an elementary error that is easily avoided. We describe a method to estimate the number and nature of orbits needed for the library. We provide an example that shows that this prescription is adequate, in the sense that the range of allowable BH masses is not artificially narrowed by use of too few orbits. A second point raised by critics is that kinematic data are generally obtained with insufficient spatial resolution to obtain a reliable mass. We make the distinction between unreliable determinations and imprecise ones. We show that there are several different properties of a kinematic dataset that can lead to imprecise BH determinations, but none of the attributes we have investigated leads to an unreliable determination. In short, the degree to which the BH radius of influence is resolved by spectroscopic observations is already reflected in the BH-mass error envelope, and is not a hidden source of error. The BH masses published by our group and the Leiden group are reliable.

This paper investigates chaos and chaotic phase mixing in triaxial Dehnen potentials which have been proposed to describe realistic ellipticals. Earlier work is extended by exploring the effects of (1) variable axis ratios, (2) `graininess' associated with stars and bound substructures, idealised as friction and white noise, and (3) large-scale organised motions presumed to induce near-random forces idealised as coloured noise with finite autocorrelation time. Three important conclusions are: (1) not all the chaos can be attributed to the cusp; (2) significant chaos can persist even for axisymmetric systems; and (3) introducing a supermassive black hole can increase both the relative number of chaotic orbits and the size of the largest Lyapunov exponent. Sans perturbations, distribution functions associated with initially localised chaotic ensembles evolve exponentially towards a nearly time-independent form at a rate L that correlates with the finite time Lyapunov exponents associated with the evolving orbits. Perturbations accelerate phase space transport by increasing the rate of phase mixing in a given phase space region and by facilitating diffusion along the Arnold web that connects different phase space regions, thus facilitating an approach towards a true equilibrium. The details of the perturbation appear unimportant. All that matters are the amplitude and the autocorrelation time, upon which there is a weak logarithmic dependence. Even comparatively weak perturbations can increase L by a factor of three or more, a fact that has potentially significant implications for violent relaxation.

Cuspy triaxial potentials admit a large number of chaotic orbits, which moreover exhibit extreme "stickiness" that makes the process of chaotic mixing surprisingly inefficient. Environmental effects, modeled as noise and/or periodic driving, help accelerate phase space transport but probably not as much as in simpler potentials. This could mean that cuspy triaxial ellipticals cannot exist as time-independent systems.

This talk provides a progress report on an extended collaboration which has aimed to address two basic questions, namely: Should one expect to see cuspy, triaxial galaxies in nature? And can one construct realistic cuspy, triaxial equilibrium models that are robust? Three technical results are described: (1) Unperturbed chaotic orbits in cuspy triaxial potentials can be extraordinarily sticky, much more so than orbits in many other three-dimensional potentials. (2) Even very weak perturbations can be important by drastically reducing, albeit not completely eliminating, this stickiness. (3) A simple toy model facilitates a simple understanding of why black holes and cusps can serve as an effective source of chaos. These results suggest that, when constructing models of galaxies using Schwarzschild's method or any analogue thereof, astronomers would be well advised to use orbital building blocks that have been perturbed by `noise' or other weak irregularities, since such building blocks are likely to be more nearly time-independent than orbits evolved in the absence of all perturbations.

(Abridged) This paper studies chaotic orbit ensembles evolved in triaxial generalisations of the Dehnen potential which have been proposed to model ellipticals with a strong density cusp that manifest significant deviations from axisymmetry. Allowance is made for a possible supermassive black hole, as well as low amplitude friction, noise, and periodic driving which can mimic irregularities associated with discreteness effects and/or an external environment. The degree of chaos is quantified by determining how (1) the relative number of chaotic orbits and (2) the size of the largest Lyapunov exponent depend on the steepness of the cusp and the black hole mass, and (3) the extent to which Arnold webs significantly impede phase space transport, both with and without perturbations. In the absence of irregularities, chaotic orbits tend to be extremely `sticky,' so that different pieces of the same chaotic orbit can behave very differently for 10000 dynamical times or longer, but even very low amplitude perturbations can prove efficient in erasing many -- albeit not all -- these differences. The implications thereof are discussed both for the structure and evolution of real galaxies and for the possibility of constructing approximate near-equilibrium models using Schwarzschild's method. Much of the observed qualitative behaviour can be reproduced with a toy potential given as the sum of an anisotropic harmonic oscillator and a spherical Plummer potential, which suggests that the results may be generic.

This paper summarises an investigation of the statistical properties of orbits escaping from three different two-degree-of-freedom Hamiltonian systems which exhibit global stochasticity. Each H=H_{0}+eH', with H_{0} integrable and eH' a nonintegrable correction, not necessarily small. For e below a critical e_{0} escapes are impossible energetically. For somewhat higher values, escape is allowed energetically but many orbits never escape and the escape probability P for a generic orbit ensemble decays exponentially. At or near a critical e_{1}>e_{0}, there is an abrupt qualitative change in behaviour. Above e_{1}, P typically exhibits (1) a rapid evolution towards a nonzero P_{0}(e) followed by (2) a much slower subsequent decay towards zero which, in at least one case, is well fit by a power law P=const x t^{m}, with m=0.35-0.40. In all three cases, P_{0} and the time T required to converge towards P_{0} scales in e-e_{1}, i.e., P_{0}=const x (e-e_{1})^{a} and T=const x (e-e_{1})^{b}, and T also scales in the size r of the region sampled for initial conditions, i.e., T=const x r^{-d}. To within statistical uncertainties, the best fit values of the critical exponents are the same for all three potentials, namely: a=0.5, b=0.4, and d=0.1, and satisfy a-b-d=0. The transitional behaviour observed near e_{1} is attributed to the breakdown of some especially significant KAM tori or cantori. The power law behaviour at late times is interpreted as reflecting intrinsic diffusion of chaotic orbits through cantori surrounding islands of regular orbits.

Schwarzschild models of a Plummer sphere were constructed and compared with analytical results to test the robustness of Schwarzschild's method and its efficiency in probing the degeneracy of the solution space. The method was then used to construct and study the nonuniqueness of triaxial equilibria with central density cusps, for which no analytical solutions are known. The stability of some models was also tested via N-body simulations.

(Abridged:) Schwarzschild's method was used to construct equilibrium solutions to the collisionless Boltzmann equation corresponding to a Plummer sphere, which were compared with analytical results to test the robustness of the numerical method and its efficiency in probing the degeneracy of the solution space. The method was then used to construct triaxial equilibria with no known analytical solutions and to study their nonuniqueness. The model that was studied, the triaxial Dehnen potential, contains a central density cusp and admits both regular and chaotic orbits. It was found that, for a weak density cusp, self-consistent models do not exist if the chaotic orbits are assumed to be completely mixed so as to yield a time-independent building-block: only the innermost 65% of the mass can be mixed. In these inner regions, it is possible to obtain alternative solutions that contain significantly different numbers of chaotic orbits, yet yield (at least approximately) the same mass density distribution. However, it is not clear whether these solutions are truly time-independent, since the unmixed chaotic orbits in the outer regions, which do not sample an invariant measure, can cause secular evolution. When using Schwarzschild's method, one must be very careful to sample accessible phase space as comprehensively and as densely as possible, while ensuring that each orbit is a truly time-independent building block. Some of the numerical equilibria were sampled to generate initial conditions for N-body simulations to test the structural stability of the models. Preliminary work showed that no catastrophic evolution takes place, but there is a weak tendency for the configuration to become more nearly axisymmetric over several dynamical times. It is not yet clear whether this tendency is real or a numerical artifact.