In the era of "big data" and with the advent of web 2.0 technologies, ESASky (http://sky.esa.int) aims at providing a modern and visual way to access astronomical science-ready data products and metadata. The main goal of the application is to simplify the interaction between the scientific community and the ever-growing amount of data collected over the past decades from the most important astronomy missions. The ESASky concept is to offer a complementary scientific application to the more-traditional table-oriented exploitation of astronomical data, by allowing a more natural and visual approach and enabling the exploration of astronomical objects across the entire electromagnetic spectrum. To fulfill this goal, ESASky provides a multiwavelength interface to a set of astronomy data from an increasing number of missions and surveys, with the intention of becoming the single-point of entry to perform visual analysis and cross-matching among different energy ranges. A lot of effort has been invested on the design of a user-friendly, responsive Graphical User Interface (GUI) by the definition and optimisation of algorithms running behind each visual feature offered. In this contribution, we describe in detail the design and solutions adopted for the technical challenges arising during the development. We present the data services and features implemented in the latest version of ESASky (v2.1), including a Mission Planning Tool to support current James Webb Space Telescope (JWST) planning, the possibility to search for observations of Solar System Objects (planets, comets and moons) taken by astronomy missions, the integration of the SAO/NASA Astrophysics Data System (ADS) publication system in the ESASky GUI, and the retrieval of metadata and data products available within a specific region of the sky.

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.

-In the Alvarez-Macovski method [R.E. Alvarez and A. Macovski, Phys. Med. Biol., 1976, 733-744], the attenuation coefficient is approximated as a linear combination of functions of energy multiplied by coefficients that depend on the material composition at points within the object. The method then computes the line integrals of the basis set coefficient from measurements with different x-ray spectra. This paper shows that the transformation from photon counting detector data with pileup to the line integrals can become ill-conditioned under some circumstances leading to highly increased noise. Methods: An idealized model that includes pileup and quantum noise is used. The noise variance of the line integral estimates is computed using the Cram{\`e}r-Rao lower bound (CRLB). The CRLB is computed as a function of object thickness for photon counting detector data with three and four bin pulse height analysis (PHA) and low and high pileup. Results: With four bin PHA data the transformation is well conditioned with either high or low pileup. With three bin PHA and high pileup, the transformation becomes ill-conditioned for specific values of object attenuation. At these values the CRLB variance increases by approximately 10 5 compared with the four bin PHA or low pileup results. The condition number of the forward transformation matrix also shows a spike at those attenuation values. Conclusion: Designers of systems using counting detectors should study the stability of the line integral estimator output with their data.

Experimental efforts to measure neutrinos by radio-frequency (RF) signals resulting from neutrino interactions in-ice have intensified over the last decade. Recent calculations indicate that one may dramatically improve the sensitivity of ultra-high energy ("UHE"; >EeV) neutrino experiments via detection of radio waves trapped along the air-ice surface. Detectors designed to observe the "Askaryan effect" currently search for RF electromagnetic pulses propagating through bulk ice, and could therefore gain sensitivity if signals are confined to the ice-air boundary. To test the feasibilty of this scenario, measurements of the complex radio-frequency properties of several air-dielectric interfaces were performed for a variety of materials. Two-dimensional surfaces of granulated fused silica (sand), both in the lab as well as occurring naturally, water doped with varying concentrations of salt, natural rock salt formations, granulated salt and ice itself were studied, both in North America and also Antarctica. In no experiment do we observe unambiguous surface wave propagation, as would be evidenced by signals travelling with reduced signal loss and/or superluminal velocities, compared to conventional EM wave propagation.

This paper presents statistics of the line-doubling phenomenon in a sample of 81 long-period variable (LPV) stars of various periods, spectral types and brightness ranges. When correlated with a mask mimicking a K0III spectrum, 54% of the sample stars clearly showed a double-peaked cross-correlation profile around maximum light, reflecting double absorption lines. Several pieces of evidence are presented that point towards the double absorption lines as being caused by the propagation of a shock wave through the photosphere. The observation of the Balmer lines appearing in emission around maximum light in these stars corroborates the presence of a shock wave. The observed velocity discontinuities, ranging between 10 and 25 km/s, are not correlated with the brightness ranges. A comparison with the center-of-mass (COM) velocity obtained from submm CO lines originating in the circumstellar envelope reveals that the median velocity between the red and blue peaks is blueshifted with respect to the COM velocity, as expected if the shock moves upwards.

A tomographic method is described that makes it possible to follow the propagation of shock waves across the photosphere of long-period variable stars. The method relies on the correlation of the observed spectrum with numerical masks probing layers of different atmospheric depths. The formation depth of spectral lines is derived from synthetic spectra of non-variable red giant stars. When applied to Mira stars around maximum light, the tomographic masks reveal that the deepest photospheric layers are generally characterized by blueshifted absorption lines (attesting their upward motion), whereas the uppermost layers generally exhibit redshifted absorption lines (indicating their infalling motion). Double absorption lines are found in intermediate layers, where the shock front is located. At later phases, the shock front is seen moving towards upper layers, until it leaves the photosphere.

In this paper HIPPARCOS astrometric and kinematic data are used to calibrate both infrared luminosities and kinematical parameters of Long Period Variable stars (LPVs). Individual absolute K and IRAS 12 and 25 luminosities of 800 LPVs are determined and made available in electronic form. The estimated mean kinematics is analyzed in terms of galactic populations. LPVs are found to belong to galactic populations ranging from the thin disk to the extended disk. An age range and a lower limit of the initial mass is given for stars of each population. A difference of 1.3mag in K for the upper limit of the Asymptotic Giant Branch is found between the disk and old disk galactic populations, confirming its dependence on the mass in the main sequence. LPVs with a thin envelope are distinguished using the estimated mean IRAS luminosities. The level of attraction (in the classification sense) of each group for the usual classifying parameters of LPVs (variability and spectral types) is examined.

This paper is the first one in a series devoted to the study of the dynamics of the atmospheres of long-period variable stars. Results from a two-month-long monitoring of the Mira variables RT Cyg and X Oph around maximum light with the ELODIE spectrograph at the Haute-Provence Observatory are presented. The monitoring covers phases 0.80 to 1.16 for RT Cyg and phases 0.83 to 1.04 for X Oph. The cross-correlation profile of the spectrum of RT Cyg with a K0 III mask confirms that the absorption lines of RT Cyg in the optical domain appear double around maximum light. No line doubling was found in the optical spectrum of X Oph around maximum light, indicating that this feature is not common to all long-period variables. This paper also presents the application to RT Cyg of a new tomographic technique deriving the velocity field across the atmosphere by cross-correlating the optical spectrum with numerical masks constructed from synthetic spectra and probing layers of increasing depths. This technique reveals that both the temporal evolution of the line doubling, and its variation with depth in the atmosphere of RT Cyg, are consistent with the ``Schwarzschild scenario''. This scenario relates the temporal evolution of the red and blue peaks of the double absorption lines to the progression of a shock wave in the atmosphere. The temporal evolution of the Balmer Halpha, H beta, Hgamma and Hdelta emission lines around maximum light is also presented for RT Cyg and X Oph. The velocity variations of Halpha and of the absorption lines are discussed in the framework of two competing models for the formation of Balmer emission lines in long-period variable stars.

Based on a new and large sample of optical and near-infrared spectra obtained at the Mount Stromlo and Siding Spring Observatories (Lancon & Wood 1998; Lancon & Wood, in preparation), spectrophotometric properties of cool oxygen- and carbon-rich Long Period Variables and supergiants are presented. Temperatures of oxygen-rich stars are assigned by comparison with synthetic spectra computed from up-to-date oxygen-rich model atmosphere grids. The existence of reliable optical and near-infrared temperature indicators is investigated. A narrow relation between the bolometric correction BC(I) and the broad-band colour I-J is obtained for oxygen-rich cool stars. The ability of specific near-infrared indices to separate luminosity classes, atmospheric chemistry or variability subtypes is discussed. Some comments are also given on extinction effects, water band strengths in Long Period Variables and the evaluation of 12CO/13CO ratio in red giants.

The absolute K magnitudes and kinematic parameters of about 350 oxygen-rich Long-Period Variable stars are calibrated, by means of an up-to-date maximum-likelihood method, using Hipparcos parallaxes and proper motions together with radial velocities and, as additional data, periods and V-K colour indices. Four groups, differing by their kinematics and mean magnitudes, are found. For each of them, we also obtain the distributions of magnitude, period and de-reddened colour of the base population, as well as de-biased period-luminosity-colour relations and their two-dimensional projections. The SRa semiregulars do not seem to constitute a separate class of LPVs. The SRb appear to belong to two populations of different ages. In a PL diagram, they constitute two evolutionary sequences towards the Mira stage. The Miras of the disk appear to pulsate on a lower-order mode. The slopes of their de-biased PL and PC relations are found to be very different from the ones of the Oxygen Miras of the LMC. This suggests that a significant number of so-called Miras of the LMC are misclassified. This also suggests that the Miras of the LMC do not constitute a homogeneous group, but include a significant proportion of metal-deficient stars, suggesting a relatively smooth star formation history. As a consequence, one may not trivially transpose the LMC period-luminosity relation from one galaxy to the other.

From near-infrared, narrow-band photometry of 256 oxygen-rich Mira variables we obtain evidence about the loops that these stars follow in colour-colour diagrams. We also find a phase lag between indices related to molecular band-strength of titanium oxide and vanadium oxide. We compute colours for normal M-giants and Miras using hydrostatic and hydrodynamic model atmospheres and very extensive up-to-date line lists. Normal M-giants colours are well reproduced, reaching a high quantitative agreement with observations for spectral types earlier than M7. The out-of-phase variations of the various spectral features of Miras are also acceptably reproduced, despite limitations in the modelling. This enables us to confirm that the phase lag phenomenon results from the propagation of perturbations in the extended atmosphere. It opens new perspectives in the spectral modelling of Miras.

HIPPARCOS astrometric and kinematical data of oxygen-rich Mira variables are used to calibrate absolute near-infrared magnitudes and kinematic parameters. Two sets of near-infrared magnitudes compiled from different authors are used: broad-band K and narrow-band photometric measurements at 1.04 micron (104 filter). Three distinct classes of stars with different kinematics and scale height have been identified. The two most significant groups present characteristics close to the ones usually assigned to extended/thick disk-halo population and old disk population respectively, and thus they might differ by their metallicity abundance. They exhibit different period distributions, as expected if these two groups actually correspond to populations of distinct initial masses, ages and metallicities. Two parallel period-luminosity relations are found in K as well as in 104, one for each significant population. The shift between these relations is interpreted as the consequence of the effects of metallicity abundance on the luminosity.