The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Convection is one of the fundamental energy transport processes in physics and astrophysics, and its description is central to allstellar models. In the context of stellar astrophysics, the mixing length theory is the most successful approximation to handle theconvection zones inside the stars because of its simplicity and rapidity. The price to pay is the mixing length parameter that isintroduced to derive the velocity of convective elements, the temperature gradients in the convective regions and finally the energy flux carried by convection. The mixing length is a free parameter that needs to be calibrated on observational data. Pasetto et al. (2014) have proposed a new theory that determines all the properties of convective regions and the convective transport of energy with no need for a free parameter. In this study, we aim to discuss the merits of this new approach and the limits of its applicability in comparison with the mixing length theory. We present an analytical and numerical investigation of the main physical assumptions made by Pasetto et al. (2014) and compare them with the counterparts of the mixing length theory. We also present here the homogeneous isotropic limit of the Pasetto et al. (2014) theory and discuss some numerical examples to address and clarify misconceptions often associated with the new formalism.

(Abridged summary) Responding to ESA's Voyage 2050 call to define the long-term plan for the future space missions that will address the astrophysics science questions during the 2035-2050 cycle, we propose a dedicated, ultra-deep spectroscopic survey in the near infrared (NIR), that will target a mass-limited sample of galaxies during two of the most fundamental epochs of cosmic evolution: the formation of the first galaxies (at z>6; cosmic dawn), and at the peak of galaxy formation activity (between redshift z=1 and 3; cosmic noon). By way of NIR observations, it is possible to study the Lyman-alpha region in the former, and the optical rest-frame in the latter, allowing us to extract fundamental observables such as gas and stellar kinematics, chemical abundances, and ages, providing a unique legacy database covering these two crucial stages of cosmic evolution. A dedicated, space-based facility will overcome the challenges faced by ground-based telescopes, no matter how large the aperture, or the reduced field of view and low multiplex factor of the best space-based instrument in the near future, namely NIRSpec at the JWST. Our project (codename Chronos) aims to produce about 1 million high quality spectra, with a high S/N in the continuum, where information about the underlying stellar populations is encoded.

In this paper we derive a full expression for the propagation of multiplicative and additive shape measurement biases into the cosmic shear power spectrum. In doing so we identify several new terms that are associated with selection effects, as well as cross-correlation terms between the multiplicative and additive biases and the shear field. The computation of the resulting bias in the shear power spectrum scales as the fifth power of the maximum multipole considered. Consequently the calculation is unfeasible for large l-modes, and the only tractable way to assess the full impact of shape measurement biases on cosmic shear power spectrum is through forward modelling of the effects. To linear order in bias parameters the shear power spectrum is only affected by the mean of the multiplicative bias field over a survey and the cross correlation between the additive bias field and the shear field. If the mean multiplicative bias is zero then second order convolutive terms are expected to be orders of magnitude smaller.

Gaia is a revolutionary space mission developed by ESA and is delivering 5 parameter astrometry, photometry and radial velocities over the whole sky with astrometric accuracies down to a few tens of micro-arcseconds. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. This problem can be overcome by switching to the Near Infra-Red (NIR) but this is not possible with silicon CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires the detectors operate in Time Delayed Integration (TDI) mode or similar.

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.

By taking advantage of the superb measurements of position and velocity for an unprecedented large number of stars provided in Gaia DR2, we have generated the first maps of the rotation velocity, $V_{\rm rot}$, and vertical velocity, $V_{\rm z}$, distributions as a function of the Galactocentric radius, $R_{\rm gal}$, across a radial range of $5<R_{\rm gal}<12$~kpc. In the $R-V_{\rm rot}$ map, we have identified many diagonal ridge features, which are compared with the location of the spiral arms and the expected outer Lindblad resonance of the Galactic bar. We have detected also radial wave-like oscillations of the peak of the vertical velocity distribution.

Aims. The Radial Velocity Spectrometer (RVS) on board the ESA satellite mission Gaia has no calibration device. Therefore, the radial velocity zero point needs to be calibrated with stars that are proved to be stable at a level of 300 m/s during the Gaia observations. Methods. We compiled a dataset of ~71000 radial velocity measurements from five high-resolution spectrographs. A catalogue of 4813 stars was built by combining these individual measurements. The zero point was established using asteroids. Results. The resulting catalogue has seven observations per star on average on a typical time baseline of six years, with a median standard deviation of 15 m/s. A subset of the most stable stars fulfilling the RVS requirements was used to establish the radial velocity zero point provided in Gaia Data Release 2. The stars that were not used for calibration are used to validate the RVS data.

This paper presents the specification, design, and development of the Radial Velocity Spectrometer (RVS) on the European Space Agency's Gaia mission. Starting with the rationale for the full six dimensions of phase space in the dynamical modelling of the Galaxy, the scientific goals and derived top-level instrument requirements are discussed, leading to a brief description of the initial concepts for the instrument. The main part of the paper is a description of the flight RVS, considering the optical design, the focal plane, the detection and acquisition chain, and the as-built performance drivers and critical technical areas. After presenting the pre-launch performance predictions, the paper concludes with the post-launch developments and mitigation strategies, together with a summary of the in-flight performance at the end of commissioning.

For Gaia DR2 (GDR2), 280 million spectra, collected by the RVS instrument on-board Gaia, were processed and median radial velocities were derived for 9.8 million sources brighter than Grvs = 12 mag. This paper describes the validation and properties of the median radial velocities published in GDR2. Quality tests and filters are applied to select, from the 9.8 million radial velocities, those with the quality to be published in GDR2. The accuracy of the selected sample is assessed with respect to ground-based catalogues. Its precision is estimated using both ground-based catalogues and the distribution of the Gaia radial velocity uncertainties. GDR2 contains median radial velocities for 7 224 631 stars, with Teff in the range [3550, 6900] K, which passed succesfully the quality tests. The published median radial velocities provide a full sky-coverage and have a completness with respect to the astrometric data of 77.2\% (for $G \leq 12.5$ mag). The median radial velocity residuals with respect to the ground-based surveys vary from one catalogue to another, but do not exceed a few 100s m/s. In addition, the Gaia radial velocities show a positive trend as a function of magnitude, which starts around Grvs $\sim 9$ mag and reaches about $+500$ m/s at Grvs $= 11.75$ mag. The overall precision, estimated from the median of the Gaia radial velocity uncertainties, is 1.05 km/s. The radial velocity precision is function of many parameters, in particular the magnitude and effective temperature. For bright stars, Grvs in [4, 8] mag, the precision is in the range 200-350 m/s, which is about 3 to 5 times more precise than the pre-launch specification of 1 km/s. At the faint end, Grvs = 11.75 mag, the precisions for Teff = 5000 K and 6500 K are respectively 1.4 km/s and 3.7 km/s.

The European Space Agency Gaia satellite was launched into orbit around L2 in December 2013. This ambitious mission has strict requirements on residual systematic errors resulting from instrumental corrections in order to meet a design goal of sub-10 microarcsecond astrometry. During the design and build phase of the science instruments, various critical calibrations were studied in detail to ensure that this goal could be met in orbit. In particular, it was determined that the video-chain offsets on the analogue side of the analogue-to-digital conversion electronics exhibited instabilities that could not be mitigated fully by modifications to the flight hardware. We provide a detailed description of the behaviour of the electronic offset levels on microsecond timescales, identifying various systematic effects that are known collectively as offset non-uniformities. The effects manifest themselves as transient perturbations on the gross zero-point electronic offset level that is routinely monitored as part of the overall calibration process. Using in-orbit special calibration sequences along with simple parametric models, we show how the effects can be calibrated, and how these calibrations are applied to the science data. While the calibration part of the process is relatively straightforward, the application of the calibrations during science data processing requires a detailed on-ground reconstruction of the readout timing of each charge-coupled device (CCD) sample on each device in order to predict correctly the highly time-dependent nature of the corrections. We demonstrate the effectiveness of our offset non-uniformity models in mitigating the effects in Gaia data. We demonstrate for all CCDs and operating instrument and modes on board Gaia that the video-chain noise-limited performance is recovered in the vast majority of science samples.

The Gaia Data Release 2 contains the 1st release of radial velocities complementing the kinematic data of a sample of about 7 million relatively bright, late-type stars. Aims: This paper provides a detailed description of the Gaia spectroscopic data processing pipeline, and of the approach adopted to derive the radial velocities presented in DR2. Methods: The pipeline must perform four main tasks: (i) clean and reduce the spectra observed with the Radial Velocity Spectrometer (RVS); (ii) calibrate the RVS instrument, including wavelength, straylight, line-spread function, bias non-uniformity, and photometric zeropoint; (iii) extract the radial velocities; and (iv) verify the accuracy and precision of the results. The radial velocity of a star is obtained through a fit of the RVS spectrum relative to an appropriate synthetic template spectrum. An additional task of the spectroscopic pipeline was to provide 1st-order estimates of the stellar atmospheric parameters required to select such template spectra. We describe the pipeline features and present the detailed calibration algorithms and software solutions we used to produce the radial velocities published in DR2. Results: The spectroscopic processing pipeline produced median radial velocities for Gaia stars with narrow-band near-IR magnitude Grvs < 12 (i.e. brighter than V~13). Stars identified as double-lined spectroscopic binaries were removed from the pipeline, while variable stars, single-lined, and non-detected double-lined spectroscopic binaries were treated as single stars. The scatter in radial velocity among different observations of a same star, also published in DR2, provides information about radial velocity variability. For the hottest (Teff > 7000 K) and coolest (Teff < 3500 K) stars, the accuracy and precision of the stellar parameter estimates are not sufficient to allow selection of appropriate templates. [Abridged]

We investigate the vertical metallicity gradients of five mono-age stellar populations between 0 and 11 Gyr for a sample of 18 435 dwarf stars selected from the cross-matched Tycho-Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiment (RAVE) Data Release 5. We find a correlation between the vertical metallicity gradients and age, with no vertical metallicity gradient in the youngest population and an increasingly steeper negative vertical metallicity gradient for the older stellar populations. The metallicity at disc plane remains almost constant between 2 and 8 Gyr, and it becomes significantly lower for the $8 < \tau \leqslant 11$ Gyr population. The current analysis also reveals that the intrinsic dispersion in metallicity increases steadily with age. We discuss that our results are consistent with a scenario that (thin) disc stars formed from a flaring (thin) star-forming disc.

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.

Previous studies have found that the Galactic rotation velocity-metallicity (V-[Fe/H]) relations for the thin and thick disk populations show negative and positive slopes, respectively. The first Gaia Data Release includes the Tycho-Gaia Astrometric Solution (TGAS) information, which we use to analyze the V-[Fe/H] relation for a strictly selected sample with high enough astrometric accuracy. We aim to arrive at an explanation for the slopes of the V-[Fe/H] relationship. We measure the V-[Fe/H] relation for thin and thick disk stars classified on the basis of their [$\alpha$/Fe] and [Fe/H] abundances. We find dV/d[Fe/H]= -18 +/- 2 km/s/dex for stars in the thin disk and dV/d[Fe/H]= +23 +/- 10 km/s/dex for thick disk stars, so we confirm the different signs for the slopes. The negative value of dV/d[Fe/H] for thick disk stars is consistent with previous studies, but the combination of TGAS and APOGEE data provide higher precision, even though systematic errors could exceed +/-5 km/s/dex . Our average measurement of dV/d[Fe/H] for local thick disk stars shows a somewhat flatter slope than the previous studies, but we confirm a significant spread and a dependence of the slope on the [alpha/Fe] ratio of the stars. Using a simple N-body model, we demonstrate that the observed trend for the thick and thin disk can be explained by the observed radial metallicity gradients and the correlation between orbital eccentricity and metallicity in the thick disk.

Euclid is a space-based optical/near-infrared survey mission of the European Space Agency (ESA) to investigate the nature of dark energy, dark matter and gravity by observing the geometry of the Universe and on the formation of structures over cosmological timescales. Euclid will use two probes of the signature of dark matter and energy: Weak gravitational Lensing, which requires the measurement of the shape and photometric redshifts of distant galaxies, and Galaxy Clustering, based on the measurement of the 3-dimensional distribution of galaxies through their spectroscopic redshifts. The mission is scheduled for launch in 2020 and is designed for 6 years of nominal survey operations. The Euclid Spacecraft is composed of a Service Module and a Payload Module. The Service Module comprises all the conventional spacecraft subsystems, the instruments warm electronics units, the sun shield and the solar arrays. In particular the Service Module provides the extremely challenging pointing accuracy required by the scientific objectives. The Payload Module consists of a 1.2 m three-mirror Korsch type telescope and of two instruments, the visible imager and the near-infrared spectro-photometer, both covering a large common field-of-view enabling to survey more than 35% of the entire sky. All sensor data are downlinked using K-band transmission and processed by a dedicated ground segment for science data processing. The Euclid data and catalogues will be made available to the public at the ESA Science Data Centre.

Gaia's Radial Velocity Spectrometer (RVS) has been operating in routine phase for over one year since initial commissioning. RVS continues to work well but the higher than expected levels of straylight reduce the limiting magnitude. The end-of-mission radial-velocity (RV) performance requirement for G2V stars was 15 km/s at V = 16.5 mag. Instead, 15 km/s precision is achieved at 15 < V < 16 mag, consistent with simulations that predict a loss of 1.4 mag. Simulations also suggest that changes to Gaia's onboard software could recover ~0.14 mag of this loss. Consequently Gaia's onboard software was upgraded in April 2015. The status of this new commissioning period is presented, as well as the latest scientific performance of the on-ground processing of RVS spectra. We illustrate the implications of the RVS limiting magnitude on Gaia's view of the Milky Way's halo in 6D using the Gaia Universe Model Snapshot (GUMS).

Euclid-VIS is the large format visible imager for the ESA Euclid space mission in their Cosmic Vision program, scheduled for launch in 2020. Together with the near infrared imaging within the NISP instrument, it forms the basis of the weak lensing measurements of Euclid. VIS will image in a single r+i+z band from 550-900 nm over a field of view of ~0.5 deg2. By combining 4 exposures with a total of 2260 sec, VIS will reach to deeper than mAB=24.5 (10sigma) for sources with extent ~0.3 arcsec. The image sampling is 0.1 arcsec. VIS will provide deep imaging with a tightly controlled and stable point spread function (PSF) over a wide survey area of 15000 deg2 to measure the cosmic shear from nearly 1.5 billion galaxies to high levels of accuracy, from which the cosmological parameters will be measured. In addition, VIS will also provide a legacy dataset with an unprecedented combination of spatial resolution, depth and area covering most of the extra-Galactic sky. Here we will present the results of the study carried out by the Euclid Consortium during the period up to the Critical Design Review.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We present here the first stellar models on the Hertzsprung-Russell diagram (HRD), in which convection is treated according to the novel scale-free convection theory (SFC theory) by Pasetto et al. (2014). The aim is to compare the results of the new theory with those from the classical, calibrated mixing-length (ML) theory to examine differences and similarities. We integrate the equations describing the structure of the atmosphere from the stellar surface down to a few percent of the stellar mass using both ML theory and SFC theory. The key temperature over pressure gradients, the energy fluxes, and the extension of the convective zones are compared in both theories. The analysis is first made for the Sun and then extended to other stars of different mass and evolutionary stage. The results are adequate: the SFC theory yields convective zones, temperature gradients of the ambient and of the convective element, and energy fluxes that are very similar to those derived from the "calibrated" MT theory for main sequence stars. We conclude that the old scale dependent ML theory can now be replaced with a self-consistent scale-free theory able to predict correct results, one which is simpler and more physically grounded than the ML theory. Fundamentally, the SFC theory offers a deeper insight of the underlying physics than numerical simulations.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

One of the most powerful techniques to study the dark sector of the Universe is weak gravitational lensing. In practice, to infer the reduced shear, weak lensing measures galaxy shapes, which are the consequence of both the intrinsic ellipticity of the sources and of the integrated gravitational lensing effect along the line of sight. Hence, a very large number of galaxies is required in order to average over their individual properties and to isolate the weak lensing cosmic shear signal. If this `shape noise' can be reduced, significant advances in the power of a weak lensing surveys can be expected. This paper describes a general method for extracting the probability distributions of parameters from catalogues of data using Voronoi cells, which has several applications, and has synergies with Bayesian hierarchical modelling approaches. This allows us to construct a probability distribution for the variance of the intrinsic ellipticity as a function of galaxy property using only photometric data, allowing a reduction of shape noise. As a proof of concept the method is applied to the CFHTLenS survey data. We use this approach to investigate trends of galaxy properties in the data and apply this to the case of weak lensing power spectra.

We investigate the kinematic signatures induced by spiral and bar structure in a set of simulations of Milky Way-sized spiral disc galaxies. The set includes test particle simulations that follow a quasi-stationary density wave-like scenario with rigidly rotating spiral arms, and $N$-body simulations that host a bar and transient, co-rotating spiral arms. From a location similar to that of the Sun, we calculate the radial, tangential and line-of-sight peculiar velocity fields of a patch of the disc and quantify the fluctuations by computing the power spectrum from a two-dimensional Fourier transform. We find that the peculiar velocity power spectrum of the simulation with a bar and transient, co-rotating spiral arms fits very well to that of APOGEE red clump star data, while the quasi-stationary density wave spiral model without a bar does not. We determine that the power spectrum is sensitive to the number of spiral arms, spiral arm pitch angle and position with respect to the spiral arm. However, it is necessary to go beyond the line of sight velocity field in order to distinguish fully between the various spiral models with this method. We compute the power spectrum for different regions of the spiral discs, and discuss the application of this analysis technique to external galaxies.

In this paper we investigate the impact that realistic scale-dependence systematic effects may have on cosmic shear tomography. We model spatially varying residual ellipticity and size variations in weak lensing measurements and propagate these through to predicted changes in the uncertainty and bias of cosmological parameters. We show that the survey strategy - whether it is regular or randomised - is an important factor in determining the impact of a systematic effect: a purely randomised survey strategy produces the smallest biases, at the expense of larger parameter uncertainties, and a very regularised survey strategy produces large biases, but unaffected uncertainties. However, by removing, or modelling, the affected scales (l-modes) in the regular cases the biases are reduced to negligible levels. We find that the integral of the systematic power spectrum is not a good metric for dark energy performance, and we advocate that systematic effects should be modelled accurately in real space, where they enter the measurement process, and their effect subsequently propagated into power spectrum contributions.

Radiation damage to space-based Charge-Coupled Device (CCD) detectors creates defects which result in an increasing Charge Transfer Inefficiency (CTI) that causes spurious image trailing. Most of the trailing can be corrected during post-processing, by modelling the charge trapping and moving electrons back to where they belong. However, such correction is not perfect -- and damage is continuing to accumulate in orbit. To aid future development, we quantify the limitations of current approaches, and determine where imperfect knowledge of model parameters most degrade measurements of photometry and morphology. As a concrete application, we simulate $1.5\times10^{9}$ "worst case" galaxy and $1.5\times10^{8}$ star images to test the performance of the Euclid visual instrument detectors. There are two separable challenges: If the model used to correct CTI is perfectly the same as that used to add CTI, $99.68$ % of spurious ellipticity is corrected in our setup. This is because readout noise is not subject to CTI, but gets over-corrected during correction. Second, if we assume the first issue to be solved, knowledge of the charge trap density within $\Delta\rho/\rho\!=\!(0.0272\pm0.0005)$ %, and the characteristic release time of the dominant species to be known within $\Delta\tau/\tau\!=\!(0.0400\pm0.0004)$ % will be required. This work presents the next level of definition of in-orbit CTI calibration procedures for Euclid.

We compare the stellar motion around a spiral arm created in two different scenarios, transient/co-rotating spiral arms and density-wave-like spiral arms. We generate Gaia mock data from snapshots of the simulations following these two scenarios using our stellar population code, SNAPDRAGONS, which takes into account dust extinction and the expected Gaia errors. We compare the observed rotation velocity around a spiral arm similar in position to the Perseus arm, and find that there is a clear difference in the velocity features around the spiral arm between the co-rotating spiral arm and the density-wave-like spiral arm. Our result demonstrates that the volume and accuracy of the Gaia data are sufficient to clearly distinguish these two scenarios of the spiral arms.

Euclid is the next ESA mission devoted to cosmology. It aims at observing most of the extragalactic sky, studying both gravitational lensing and clustering over $\sim$15,000 square degrees. The mission is expected to be launched in year 2020 and to last six years. The sheer amount of data of different kinds, the variety of (un)known systematic effects and the complexity of measures require efforts both in sophisticated simulations and techniques of data analysis. We review the mission main characteristics, some aspects of the the survey and highlight some of the areas of interest to this meeting

In this paper we present the testing of a back-illuminated development Euclid Visible Instrument (VIS) Charge-Coupled Device (CCD) to measure the intrinsic CCD Point Spread Function (PSF) characteristics using a novel modelling technique. We model the optical spot projection system and the CCD273-84 PSF jointly. We fit a model using Bayesian posterior probability density function, sampling to all available data simultaneously. The generative model fitting is shown, using simulated data, to allow good parameter estimations even when these data are not well sampled. Using available spot data we characterise a CCD273-84 PSF as a function of wavelength and intensity. The CCD PSF kernel size was found to increase with increasing intensity and decreasing wavelength.

We have observed an N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way like barred spiral galaxy. We present a simple method that samples N-body model particles into mock Gaia stellar observations and takes into account stellar populations, dust extinction and Gaia's science performance estimates. We examine the kinematics around a nearby spiral arm at a similar position to the Perseus arm at three lines of sight in the disc plane; (l,b)=(90,0), (120,0) and (150,0) degrees. We find that the structure of the peculiar kinematics around the co-rotating spiral arm, which is found in Kawata et al. (2014b), is still visible in the observational data expected to be produced by Gaia despite the dust extinction and expected observational errors of Gaia. These observable kinematic signatures will enable testing whether the Perseus arm of the Milky Way is similar to the co-rotating spiral arms commonly seen in N-body simulations.

Radial migration is defined as the change in guiding centre radius of stars and gas caused by gains or losses of angular momentum that result from gravitational interaction with non-axisymmetric structure. This has been shown to have significant impact on the metallicity distribution in galactic discs, and therefore affects the interpretation of Galactic archeology. We use a simulation of a Milky Way-sized galaxy to examine the effect of radial migration on the star and gas radial metallicity distribution. We find that both the star and gas component show significant radial migration. The stellar radial metallicity gradient remains almost unchanged but the radial metallicity distribution of the stars is broadened to produce a greater dispersion at all radii. However, the metallicity dispersion of the gas remains narrow. We find that the main drivers of the gas metallicity distribution evolution are metal enrichment and mixing: more efficient metal enrichment in the inner region maintains a negative slope in the radial metallicity distribution, and the metal mixing ensures the tight relationship of the gas metallicity with the radius. The metallicity distribution function reproduces the trend in the age-metallicity relation found from observations for stars younger than 1.0 Gyr in the Milky Way.

Stellar convection is customarily described by Mixing-Length Theory, which makes use of the mixing-length scale to express the convective flux, velocity, and temperature gradients of the convective elements and stellar medium. The mixing-length scale is taken to be proportional to the local pressure scale height, and the proportionality factor (the mixing-length parameter) must be determined by comparing the stellar models to some calibrator, usually the Sun. No strong arguments exist to suggest that the mixing-length parameter is the same in all stars and at all evolutionary phases. The aim of this study is to present a new theory of stellar convection that does not require the mixing length parameter. We present a self-consistent analytical formulation of stellar convection that determines the properties of stellar convection as a function of the physical behaviour of the convective elements themselves and of the surrounding medium. This new theory is formulated starting from a conventional solution of the Navier-Stokes/Euler equations, i.e. the Bernoulli equation for a perfect fluid, but expressed in a non-inertial reference frame co-moving with the convective elements. In our formalism the motion of stellar convective cells inside convectively-unstable layers is fully determined by a new system of equations for convection in a non-local and time-dependent formalism. We obtain an analytical, non-local, time-dependent sub-sonic solution for the convective energy transport that does not depend on any free parameter. The theory is suitable for the outer convective zones of solar type stars and stars of all mass on the main sequence band. The predictions of the new theory are compared with those from the standard mixing-length paradigm for the most accurate calibrator, the Sun, with very satisfactory results.

We develop a simple analytical criterion to investigate the role of the environment on the onset of star formation. We will consider the main external agents that influence the star formation (i.e. ram pressure, tidal interaction, Rayleigh-Taylor and Kelvin-Helmholtz instabilities) in a spherical galaxy moving through an external environment. The theoretical framework developed here has direct applications to the cases of dwarf galaxies in galaxy clusters and dwarf galaxies orbiting our Milky Way system, as well as any primordial gas-rich cluster of stars orbiting within its host galaxy. We develop an analytic formalism to solve the fluid dynamics equations in a non-inertial reference frame mapped with spherical coordinates. The two-fluids instability at the interface between a stellar system and its surrounding hotter and less dense environment is related to the star formation processes through a set of differential equations. The solution presented here is quite general, allowing us to investigate most kinds of orbits allowed in a gravitationally bound system of stars in interaction with a major massive companion. We present an analytical criterion to elucidate the dependence of star formation in a spherical stellar system (as a dwarf galaxy or a globular cluster) on its surrounding environment useful in theoretical interpretations of numerical results as well as observational applications. We show how spherical coordinates naturally enlighten the interpretation of the two-fluids instability in a geometry that directly applies to astrophysical case. This criterion predicts the threshold value for the onset of star formation in a mass vs. size space for any orbit of interest. Moreover, we show for the first time the theoretical dependencies of the different instability phenomena acting on a system in a fully analytical way.

We have observed a snapshot of our N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way-sized barred spiral galaxy in a similar way to how we can observe the Milky Way. The simulated galaxy shows a co-rotating spiral arm, i.e. the spiral arm rotates with the same speed as the circular speed. We observed the rotation and radial velocities of the gas and stars as a function of the distance from our assumed location of the observer at the three lines of sight on the disc plane, (l, b) = (90, 0), (120, 0) and (150,0) deg. We find that the stars tend to rotate slower (faster) behind (at the front of) the spiral arm and move outward (inward), because of the radial migration. However, because of their epicycle motion, we see a variation of rotation and radial velocities around the spiral arm. On the other hand, the cold gas component shows a clearer trend of rotating slower (faster) and moving outward (inward) behind (at the front of) the spiral arm, because of the radial migration. We have compared the results with the velocity of the maser sources from Reid et al. (2014), and find that the observational data show a similar trend in the rotation velocity around the expected position of the spiral arm at l = 120 deg. We also compared the distribution of the radial velocity from the local standard of the rest, V_LSR, with the APOGEE data at l = 90 deg as an example.

(Abridged) We combine NUV, optical and IR imaging of the nearby starburst galaxy M82 to explore the properties of the dust both in the interstellar medium of the galaxy and the dust entrained in the superwind. The three NUV filters of Swift/UVOT enable us to probe in detail the properties of the extinction curve in the region around the 2175A bump. The NUV colour-colour diagram strongly rules out a Calzetti-type law, which can either reflect intrinsic changes in the dust properties or in the star formation history compared to starbursts well represented by such an attenuation law. We emphasize that it is mainly in the NUV region where a standard Milky-Way-type law is preferred over a Calzetti law. The age and dust distribution of the stellar populations is consistent with the scenario of an encounter with M81 in the recent 400 Myr. The radial gradients of the NUV and optical colours in the superwind region support the hypothesis that the emission in the wind cone is driven by scattering from dust grains entrained in the ejecta. The observed wavelength dependence reveals either a grain size distribution $n(a)\propto a^{-2.5}$, where $a$ is the size of the grain, or a flatter distribution with a maximum size cutoff, suggesting that only small grains are entrained in the supernovae-driven wind.

Recent numerical N-body simulations of spiral galaxies have shown that spiral arms in N-body simulations seem to rotate at a similar speed to the local rotation speed of the stellar disc material. This in turn yields winding, transient and recurrent spiral structure, whose co-rotating nature gives rise to changes in the angular momentum (radial migration) of star particles close to the spiral arm at many radii. From high resolution N-body simulations, we highlight the evolution of strongly migrating star particles (migrators) and star particles that do not migrate (non-migrators) around a spiral arm. We investigate the individual orbit histories of migrators and non-migrators and find that there are several types of migrator and non-migrator, each with unique radial evolution. We contrast each type of orbit to establish the reasons for the differences between them. We find that the positive (negative) migrators sustain a position behind (in front of) the spiral arm, and feel continuous tangential force as long as the spiral arm persists. On the other hand, non-migrators stay close to the spiral arm, and pass or are passed by the spiral arm one or two times. Although they gain or lose the angular momentum when they are behind or in front of the spiral arm, their net angular momentum change becomes close to zero. We discuss also the long term effects of radial migration on the radial metallicity distribution and radial angular momentum and mass profiles. [Abridged]

Charge-Coupled Device (CCD) detectors, widely used to obtain digital imaging, can be damaged by high energy radiation. Degraded images appear blurred, because of an effect known as Charge Transfer Inefficiency (CTI), which trails bright objects as the image is read out. It is often possible to correct most of the trailing during post-processing, by moving flux back to where it belongs. We compare several popular algorithms for this: quantifying the effect of their physical assumptions and tradeoffs between speed and accuracy. We combine their best elements to construct a more accurate model of damaged CCDs in the Hubble Space Telescope's Advanced Camera for Surveys/Wide Field Channel, and update it using data up to early 2013. Our algorithm now corrects 98% of CTI trailing in science exposures, a substantial improvement over previous work. Further progress will be fundamentally limited by the presence of read noise. Read noise is added after charge transfer so does not get trailed - but it is incorrectly untrailed during post-processing.

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

This paper describes the definition of a typical next-generation space-based weak gravitational lensing experiment. We first adopt a set of top-level science requirements from the literature, based on the scale and depth of the galaxy sample, and the avoidance of systematic effects in the measurements which would bias the derived shear values. We then identify and categorise the contributing factors to the systematic effects, combining them with the correct weighting, in such a way as to fit within the top-level requirements. We present techniques which permit the performance to be evaluated and explore the limits at which the contributing factors can be managed. Besides the modelling biases resulting from the use of weighted moments, the main contributing factors are the reconstruction of the instrument point spread function (PSF), which is derived from the stellar images on the image, and the correction of the charge transfer inefficiency (CTI) in the CCD detectors caused by radiation damage.

Spiral galaxies are observed to exhibit a range of morphologies, in particular in the shape of spiral arms. A key diagnostic parameter is the pitch angle, which describes how tightly wound the spiral arms are. Observationally and analytically, a correlation between pitch angle and galactic shear rate has been detected. For the first time, we examine whether this effect is detected in N-body simulations by calculating and comparing pitch angles of both individual density waves and overall spiral structure in a suite of N-body simulations. We find that higher galactic shear rates produce more tightly wound spiral arms, both in individual mode patterns (density waves) and in the overall density enhancement. Although the mode pattern pitch angles by construction remain constant with time, the overall logarithmic spiral arm winds over time, which could help to explain the scatter in the relation between pitch angle versus shear seen from observations. The correlation between spiral arm pitch angle and galactic shear rate that we find in N-body simulations may also explain why late Hubble type of spiral galaxies tend to have more open arms.

Many different classes of X-ray sources contribute to the Galactic landscape at high energies. Although the nature of the most luminous X-ray emitters is now fairly well understood, the population of low-to-medium X-ray luminosity (Lx = 10^27-10^34 erg/s) sources remains much less studied, our knowledge being mostly based on the observation of local members. The advent of wide field and high sensitivity X-ray telescopes such as XMM-Newton now offers the opportunity to observe this low-to-medium Lx population at large distances. We report on the results of a Galactic plane survey conducted by the XMM-Newton Survey Science Centre (SSC). Beyond its astrophysical goals, this survey aims at gathering a representative sample of identified X-ray sources at low latitude that can be used later on to statistically identify the rest of the serendipitous sources discovered in the Milky Way. The survey is based on 26 XMM-Newton observations, obtained at |b| < 20 deg, distributed over a large range in Galactic longitudes and covering a summed area of 4 deg2. The flux limit of our survey is 2 x 10-15 erg/cm^2/s in the soft (0.5 - 2 keV) band and 1 x 10^-14 erg/cm^2/s in the hard (2 - 12 keV) band. We detect a total of 1319 individual X-ray sources. Using optical follow-up observations supplemented by cross-correlation with a large range of multi-wavelength archival catalogues we identify 316 X-ray sources. This constitutes the largest group of spectroscopically identified low latitude X-ray sources at this flux level. The majority of the identified X-ray sources are active coronae with spectral types in the range A - M at maximum distances of ~ 1 kpc. The number of identified active stars increases towards late spectral types, reaching a maximum at K. (abridged)

X-ray observations unveiled various types of radio-silent Isolated Neutron Stars (INSs), phenomenologically very diverse, e.g. the Myr old X-ray Dim INS (XDINSs) and the kyr old magnetars. Although their phenomenology is much diverse, the similar periods (P=2--10 s) and magnetic fields (~10^{14} G) suggest that XDINSs are evolved magnetars, possibly born from similar populations of supermassive stars. One way to test this hypothesis is to identify their parental star clusters by extrapolating backward the neutron star velocity vector in the Galactic potential. By using the information on the age and space velocity of the XDINS RX J1856.5-3754, we computed backwards its orbit in the Galactic potential and searched for its parental stellar cluster by means of a closest approach criterion. We found a very likely association with the Upper Scorpius OB association, for a neutron star age of 0.42+/-0.08 Myr, a radial velocity V_r^NS =67+/- 13$ km s^{-1}, and a present-time parallactic distance d_\pi^NS = 123^{+11}_{-15} pc. Our result confirms that the "true" neutron star age is much lower than the spin-down age (tau_{sd}=3.8 Myrs), and is in good agreement with the cooling age, as computed within standard cooling scenarios. The mismatch between the spin-down and the dynamical/cooling age would require either an anomalously large breaking index (n~20) or a decaying magnetic field with initial value B_0 ~ 10^{14} G. Unfortunately, owing to the uncertainty on the age of the Upper Scorpius OB association and the masses of its members we cannot yet draw firm conclusions on the estimated mass of the RX J1856.5-3754 progenitor.

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset.

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

Euclid is an ESA Cosmic-Vision wide-field-space mission which is designed to explain the origin of the acceleration of Universe expansion. The mission will investigate at the same time two primary cosmological probes: Weak gravitational Lensing (WL) and Galaxy Clustering (in particular Baryon Acoustic Oscillations, BAO). The extreme precision requested on primary science objectives can only be achieved by observing a large number of galaxies distributed over the whole sky in order to probe the distribution of dark matter and galaxies at all scales. The extreme accuracy needed requires observation from space to limit all observational biases in the measurements. The definition of the Euclid survey, aiming at detecting billions of galaxies over 15 000 square degrees of the extragalactic sky, is a key parameter of the mission. It drives its scientific potential, its duration and the mass of the spacecraft. The construction of a Reference Survey derives from the high level science requirements for a Wide and a Deep survey. The definition of a main sequence of observations and the associated calibrations were indeed a major achievement of the Definition Phase. Implementation of this sequence demonstrated the feasibility of covering the requested area in less than 6 years while taking into account the overheads of space segment observing and maneuvering sequence. This reference mission will be used for sizing the spacecraft consumables needed for primary science. It will also set the framework for optimizing the time on the sky to fulfill the primary science and maximize the Euclid legacy.

Euclid-VIS is a large format visible imager for the ESA Euclid space mission in their Cosmic Vision program, scheduled for launch in 2019. Together with the near infrared imaging within the NISP instrument it forms the basis of the weak lensing measurements of Euclid. VIS will image in a single r+i+z band from 550-900 nm over a field of view of ~0.5 deg2. By combining 4 exposures with a total of 2240 sec, VIS will reach to V=24.5 (10{\sigma}) for sources with extent ~0.3 arcsec. The image sampling is 0.1 arcsec. VIS will provide deep imaging with a tightly controlled and stable point spread function (PSF) over a wide survey area of 15000 deg2 to measure the cosmic shear from nearly 1.5 billion galaxies to high levels of accuracy, from which the cosmological parameters will be measured. In addition, VIS will also provide a legacy imaging dataset with an unprecedented combination of spatial resolution, depth and area covering most of the extra-Galactic sky. Here we will present the results of the study carried out by the Euclid Consortium during the Euclid Definition phase.

We run N-body smoothed particle hydrodynamics (SPH) simulations of a Milky Way sized galaxy. The code takes into account hydrodynamics, self-gravity, star formation, supernova and stellar wind feedback, radiative cooling and metal enrichment. The simulated galaxy is a barred-spiral galaxy consisting of a stellar and gas disc, enveloped in a static dark matter halo. Similar to what is found in our pure N-body simulation of a non-barred galaxy in Grand et. al. (2012), we find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is similar to the rotation of star particles. Compared to the non-barred case, we find that the spiral arm pattern speed is slightly faster than the rotation speed of star particles: the bar appears to boost the pattern speed ahead of the rotational velocity. We trace particle motion around the spiral arms at different radii, and demonstrate that there are star particles that are drawn towards and join the arm from behind (in front of) the arm and migrate toward the outer (inner) regions of the disc until the arm disappears as a result of their transient nature. We see this migration over the entire radial range analysed, which is a consequence of the spiral arm rotating at similar speeds to star particles at all radii, which is inconsistent with the prediction of classical density wave theory. The bar does not prevent this systematic radial migration, which is shown to largely preserve circular orbits. We also demonstrate that there is no significant offset of different star forming tracers across the spiral arm, which is also inconsistent with the prediction of classical density wave theory.

The XMM-Newton Serendipitous Ultraviolet Source Survey (XMM-SUSS) is a catalogue of ultraviolet (UV) sources detected serendipitously by the Optical Monitor (XMM-OM) on-board the XMM-Newton observatory. The catalogue contains ultraviolet-detected sources collected from 2,417 XMM-OM observations in 1-6 broad band UV and optical filters, made between 24 February 2000 and 29 March 2007. The primary contents of the catalogue are source positions, magnitudes and fluxes in 1 to 6 passbands, and these are accompanied by profile diagnostics and variability statistics. The XMM-SUSS is populated by 753,578 UV source detections above a 3 sigma signal-to-noise threshold limit which relate to 624,049 unique objects. Taking account of substantial overlaps between observations, the net sky area covered is 29-54 square degrees, depending on UV filter. The magnitude distributions peak at 20.2, 20.9 and 21.2 in UVW2, UVM2 and UVW1 respectively. More than 10 per cent of sources have been visited more than once using the same filter during XMM-Newton operation, and > 20 per cent of sources are observed more than once per filter during an individual visit. Consequently, the scope for science based on temporal source variability on timescales of hours to years is broad. By comparison with other astrophysical catalogues we test the accuracy of the source measurements and define the nature of the serendipitous UV XMM-OM source sample. The distributions of source colours in the UV and optical filters are shown together with the expected loci of stars and galaxies, and indicate that sources which are detected in multiple UV bands are predominantly star-forming galaxies and stars of type G or earlier.

We study the star forming regions in the spiral galaxy NGC4321, taking advantage of the spatial resolution (2.5 arcsec FWHM) of the Swift/UVOT camera and the availability of three UV passbands in the region 1600-3000 A, in combination with optical and IR imaging from SDSS, KPNO/Ha and Spitzer/IRAC, to obtain a catalogue of 787 star forming regions out to three disc scale lengths. We determine the properties of the young stellar component and its relationship with the spiral arms. The Ha luminosities of the sources have a strong decreasing radial trend, suggesting more massive star forming regions in the central part of the galaxy. When segregated with respect to NUV-optical colour, blue sources have a significant excess of flux in the IR at 8 micron, revealing the contribution from PAHs, although the overall reddening of these sources stays below E(B-V)=0.2 mag. The distribution of distances to the spiral arms is compared for subsamples selected according to Ha luminosity, NUV-optical colour, or ages derived from a population synthesis model. An offset is expected between these subsamples as a function of radius if the pattern speed of the spiral arm were constant - as predicted by classic density wave theory. No significant offsets are found, favouring instead a mechanism where the pattern speed has a radial dependence.

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in terms of a long-lived wave structure with a constant pattern speed in order to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. We run N-body simulations of a giant disc galaxy consisting of a pure stellar disc and a static dark matter halo, and find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace particle motion around the spiral arms. We show that particles from behind and in front of the spiral arm are drawn towards and join the arm. Particles move along the arm in the radial direction and we find a clear trend that they migrate toward the outer (inner) radii on the trailing (leading) side of the arm. Our simulations demonstrate that because the spiral arm feature is co-rotating, the particles continue to be accelerated (decelerated) by the spiral arm for long periods, which leads to strong and efficient migration, at all radii in the disc.

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in terms of a long-lived wave structure with a constant pattern speed in order to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. We run an N-body/SPH simulation of a Milky Way-sized barred disk, and find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace particle motion around the spiral arms. We show that particles from behind and in front of the spiral arm are drawn towards and join the arm. Particles move along the arm in the radial direction and we find a clear trend that they migrate toward the outer (inner) radii on the trailing (leading) side of the arm. Our simulations demonstrate that at all radii where there is a co-rotating spiral arm the particles continue to be accelerated (decelerated) by the spiral arm for long periods, which leads to strong migration.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We report the discovery of 31 blue, short period, pulsators made using data taken as part of the Rapid Temporal Survey (RATS). We find they have periods between 51-83 mins and full-amplitudes between 0.05-0.65 mag. Using the period-luminosity relationship for short period pulsating stars we determine their distance. Assuming they are pulsating in either the fundamental or first over-tone radial mode the majority are located at a distance greater than 3kpc, with several being more than 20 kpc distant. Most stars are at least 1 kpc from the Galactic plane, with three being more than 10 kpc. One is located in the direction of the Galactic anti-center and has Galactocentric distance of ~30 kpc and is ~20 kpc below the plane: they are therefore potential tracers of Galactic structure. We have obtained low-resolution spectra for a small number our targets and find they have temperatures between 7200--7900K and a metal content less than Solar. The colours of the pulsators and the spectral fits to those stars for which we have spectra indicate that they are either SX Phe or delta Scuti stars. We estimate the number of SX Phe stars in our Galaxy and find significantly fewer per unit mass than reported in massive globular clusters or dwarf spheroidal galaxies.

Cosmic shear has been identified as the method with the most potential to constrain dark energy. To capitalise on this potential it is necessary to measure galaxy shapes with great accuracy, which in turn requires a detailed model for the image blurring, the Point Spread Function (PSF). In general the PSF varies with wavelength and therefore the PSF integrated over an observing filter depends on the spectrum of the object. For a typical galaxy the spectrum varies across the galaxy image, thus the PSF depends on the position within the image. We estimate the bias on the shear due to such colour gradients by modelling galaxies using two co-centered, co-elliptical Sersic profiles, each with a different spectrum. We estimate the effect of ignoring colour gradients and find the shear bias from a single galaxy can be very large depending on the properties of the galaxy. We find that halving the filter width reduces the shear bias by a factor of about 5. We show that, to first order, tomographic cosmic shear two point statistics depend on the mean shear bias over the galaxy population at a given redshift. For a single broad filter, and averaging over a small galaxy catalogue, we find a mean shear bias which is subdominant to the predicted statistical errors for future cosmic shear surveys. However, the true mean shear bias may exceed the statistical errors, depending on how accurately the catalogue represents the observed distribution of galaxies in the cosmic shear survey. We then investigate the bias on the shear for two-filter imaging and find that the bias is reduced by at least an order of magnitude. Lastly, we find that it is possible to calibrate galaxies for which colour gradients were ignored using two-filter imaging of a fair sample of noisy galaxies, if the galaxy model is known. For a S/N of 25 the number of galaxies required in each tomographic redshift bin is of order 1e4.

X-ray observations performed with ROSAT led to the discovery of a group (seven to date) of X-ray dim and radio-silent middle-aged isolated neutron stars (a.k.a. XDINSs), which are characterised by pure blackbody spectra (kT~40-100 eV), long X-ray pulsations (P=3-12 s), and appear to be endowed with relatively high magnetic fields, (B~10d13-14 G). RBS 1774 is one of the few XDINSs with a candidate optical counterpart, which we discovered with the VLT. We performed deep observations of RBS 1774 in the R band with the VLT to disentangle a non-thermal power-law spectrum from a Rayleigh-Jeans, whose contributions are expected to be very much different in the red part of the spectrum. We did not detect the RBS 1774 candidate counterpart down to a 3 sigma limiting magnitude of R~27. The constraint on its colour, (B-R)<0.6, rules out that it is a background object, positionally coincident with the X-ray source. Our R-band upper limit is consistent with the extrapolation of the B-band flux (assuming a 3 sigma uncertainty) for a set of power-laws F_nu ~nu^alpha with spectral indeces alpha<0.07. If the optical spectrum of RBS 1774 were non-thermal, its power-law slope would be very much unlike those of all isolated neutron stars with non-thermal optical emission, suggesting that it is most likely thermal. For instance, a Rayleigh-Jeans with temperature T_O = 11 eV, for an optically emitting radius r_O=15 km and a source distance d=150 pc, would be consistent with the optical measurements. The implied low distance is compatible with the 0.04 X-ray pulsed fraction if either the star spin axis is nearly aligned with the magnetic axis or with the line of sight, or it is slightly misaligned with respect to both the magnetic axis and the line of sight by 5-10 degrees

The lack of radial velocity data in the Hipparcos catalogue was considered a significant deficiency, so when Gaia was conceived, a spectrometer was a core constituent of its payload. The Gaia Radial Velocity Spectrometer faced a number of design challenges, in particular set by the need to balance kinematic and astrophysical capability. We present an overview of the evolution of the instrument to its present form, identifying the competing technical, performance and programmatic factors which have shaped it.

The Galactic population of open clusters provides an insight into star formation in the Galaxy. The open cluster catalogue by Dias et al.(2002b) is a rich source of data, including kinematic information. This large sample made it possible to carry out a systematic analysis of 481 open cluster orbits, using parameters based on orbit eccentricity and separation from the Galactic plane. These two parameters may be indicative of origin, and we find them to be correlated. We also find them to be correlated with metallicity, another parameter suggested elsewhere to be a marker for origin in that high values of any of these two parameters generally indicates a low metallicity ([Fe/H] Solar$<-$0.2 dex). The resulting analysis points to four open clusters in the catalogue being of extra-Galactic origin by impact of high velocity cloud on the disk: Berkeley21, 32, 99, and Melotte66, with a possible further four due to this origin (NGC2158, 2420, 7789, IC1311). A further three may be due to Galactic globular cluster impact on the disk i.e of internal Galactic origin (NGC6791, 1817, and 7044).

We report the results of an optical campaign carried out by the XMM-Newton Survey Science Centre with the specific goal of identifying the brightest X-ray sources in the XMM-Newton Galactic Plane Survey of Hands et al. (2004). In addition to photometric and spectroscopic observations obtained at the ESO-VLT and ESO-3.6m, we used cross-correlations with the 2XMMi, USNO-B1.0, 2MASS and GLIMPSE catalogues to progress the identification process. Active coronae account for 16 of the 30 identified X-ray sources. Many of the identified hard X-ray sources are associated with massive stars emitting at intermediate X-ray luminosities of 10^32-34 erg/s. Among these are a very absorbed likely hyper-luminous star with X-ray/optical spectra and luminosities comparable with those of eta Carina, a new X-ray selected WN8 Wolf-Rayet star, a new Be/X-ray star belonging to the growing class of Gamma-Cas analogs and a possible supergiant X-ray binary of the kind discovered recently by INTEGRAL. One of the sources, XGPS-25 has a counterpart which exhibits HeII 4686 and Bowen CIII-NIII emission lines suggesting a quiescent or X-ray shielded Low Mass X-ray Binary, although its properties might also be consistent with a rare kind of cataclysmic variable (CV). We also report the discovery of three new CVs, one of which is a likely magnetic system. The soft (0.4-2.0 keV) band LogN-LogS curve is completely dominated by active stars in the flux range of 1x10^-13 to 1x10^-14 erg/cm2/s. In total, we are able to identify a large fraction of the hard (2-10 keV) X-ray sources in the flux range of 1x10^-12 to 1x10^-13 erg/cm2/s with Galactic objects at a rate consistent with that expected for the Galactic contribution only. (abridged)

So far, six mechanisms have been proposed to account for the Galactic disk heating. Of these, the most important appear to be a combination of scattering of stars by molecular clouds and by spiral arms. We study a further mechanism, namely, the repeated disk impact of the original Galactic Globular Cluster population up to the present. We find that Globular Clusters could have contributed at most a small fraction of the current vertical energy of the disk, as they could heat the whole disk to {$\sigma $}$_{z}$ = 5.5kms$^{{\rm -} {\rm 1}}$ (c.f. the observed 18 and 39 kms$^{{\rm -} {\rm 1}}$ for the thick and thin disks respectively). We find that the rate of rise of disk heat ($\alpha$=0.22 in \textit{$\sigma $}$_{z}$ $\sim t^{\alpha}$ with \textit t being time), is close to that found for scattering by molecular clouds.

The crossing of the Galactic disk by a Globular Cluster could produce star formation due to gravitational focussing or compression of disk material. We report on simulations of the effect on disk material which reveal that the crossing can sometimes cause local gravitational focussing of disk material. We also present the salient points of a little-known paper by Levy (2000), that shows that strong compression can result from the shock wave generated by GC disk crossing. The main thrust of our paper is a search for remnants of disk crossings by Globular Clusters. Using the gravitational potential of the Galaxy to locate the position of the most recent crossings of a subset of fifty-four Globular Clusters reveals that systematic errors and uncertainties in initial conditions limit the scope for unequivocal identification. From the subset of fifty-four, six possible search sites with the best constraints are retained for further scrutiny. Three of the six potentially promising search areas in the disk are from Globular Clusters NGC 3201, 6397 and NGC 6838, for which we cannot rule out some observed star associations observed nearby as being remnants. The three other of the six areas are too large to provide meaningful identification of remnants. Also, a possible remnant (open cluster NGC6231) is shown not to be due to Globular Cluster impact, contrary to a previous report. In a more wide-ranging screening of one hundred and fifty-five Globular Clusters we identify which Globular Clusters are compatible with being responsible for the formation of any of the Galaxy's five most prominent Star Super Clusters.

The Gaia satellite is a high-precision astrometry, photometry and spectroscopic ESA cornerstone mission, currently scheduled for launch in late 2011. Its primary science drivers are the composition, formation and evolution of the Galaxy. Gaia will achieve its scientific requirements with detailed calibration and correction for radiation damage. Microscopic models of Gaia's CCDs are being developed to simulate the charge trapping effect of radiation damage, which causes charge transfer inefficiency. The key to calculating the probability of a photoelectron being captured by a trap is the 3D electron density within each CCD pixel. However, this has not been physically modelled for Gaia CCD pixels. In this paper, the first of a series, we motivate the need for such specialised 3D device modelling and outline how its future results will fit into Gaia's overall radiation calibration strategy.

DUNE (Dark Universe Explorer) is a proposed mission to measure parameters of dark energy using weak gravitational lensing The particular challenges of both optical and infrared focal planes and the DUNE baseline solution is discussed. The DUNE visible Focal Plane Array (VFP) consists of 36 large format red-sensitive CCDs, arranged in a 9x4 array together with the associated mechanical support structure and electronics processing chains. Four additional CCDs dedicated to attitude control measurements are located at the edge of the array. All CCDs are 4096 pixel red-enhanced e2v CCD203-82 devices with square 12 $\mu$m pixels, operating from 550-920nm. Combining four rows of CCDs provides a total exposure time of 1500s. The VFP will be used in a closed-loop system by the spacecraft, which operates in a drift scan mode, in order to synchronize the scan and readout rates. The Near Infrared (NIR) FPA consists of a 5 x 12 mosaic of 60 Hawaii 2RG detector arrays from Teledyne, NIR bandpass filters for the wavelength bands Y, J, and H, the mechanical support structure, and the detector readout and signal processing electronics. The FPA is operated at a maximum temperature of 140 K for low dark current of 0.02e$-$/s. Each sensor chip assembly has 2048 x 2048 square pixels of 18 $\mu$m size (0.15 arcsec), sensitive in the 0.8 to 1.7 $\mu$m wavelength range. As the spacecraft is scanning the sky, the image motion on the NIR FPA is stabilized by a de-scanning mirror during the integration time of 300 s per detector. The total integration time of 1500 seconds is split among the three NIR wavelengths bands. DUNE has been proposed to ESA's Cosmic Vision program and has been jointly selected with SPACE for an ESA Assessment Phase which has led to the joint Euclid mission concept.

Multiwavelength studies of the seven identified X-ray dim isolated neutron stars (XDINSs) offer a unique opportunity to investigate their surface thermal and magnetic structure and the matter-radiation interaction in presence of strong gravitational and magnetic fields. As a part of an ongoing campaign aimed at a complete identification and spectral characterization of XDINSs in the optical band, we performed deep imaging with the ESO Very Large Telescope (VLT) of the field of the XDINS RBS1774 (1RXS J214303.7 +065419). The recently upgraded FORS1 instrument mounted on the VLT provided the very first detection of a candidate optical counterpart in the B band. The identification is based on a very good positional coincidence with the X-ray source (chance probability ~2E-3). The source has B=27.4 +/- 0.2 (1 sigma confidence level), and the optical flux exceeds the extrapolation of the X-ray blackbody at optical wavelengths by a factor ~35 (+/- 20 at 3sigma confidence level). This is barely compatible with thermal emission from the neutron star surface, unless the source distance is d~200-300 pc, and the star is an almost aligned rotator or its spin axis is nearly aligned with the line of sight. At the same time, such a large optical excess appears difficult to reconcile with rotation-powered magnetospheric emission, unless the source has an extremely large optical emission efficiency. The implications and possible similarities with the optical spectra of other isolated NSs are discussed.

We present the XMM-Newton Medium sensitivity Survey (XMS), including a total of 318 X-ray sources found among the serendipitous content of 25 XMM-Newton target fields. The XMS comprises four largely overlapping source samples selected at soft (0.5-2 keV), intermediate (0.5-4.5 keV), hard (2-10 keV) and ultra-hard (4.5-7.5 keV) bands, the first three of them being flux-limited. We report on the optical identification of the XMS samples, complete to 85-95%. At the intermediate flux levels sampled by the XMS we find that the X-ray sky is largely dominated by Active Galactic Nuclei. The fraction of stars in soft X-ray selected samples is below 10%, and only a few per cent for hard selected samples. We find that the fraction of optically obscured objects in the AGN population stays constant at around 15-20% for soft and intermediate band selected X-ray sources, over 2 decades of flux. The fraction of obscured objects amongst the AGN population is larger (~35-45%) in the hard or ultra-hard selected samples, and constant across a similarly wide flux range. The distribution in X-ray-to-optical flux ratio is a strong function of the selection band, with a larger fraction of sources with high values in hard selected samples. Sources with X-ray-to-optical flux ratios in excess of 10 are dominated by obscured AGN, but with a significant contribution from unobscured AGN.

We present XMM-Newton observations of the eclipsing polar V2301 Oph which cover nearly 2.5 binary orbital cycles and 2 eclipses. This polar is believed to have the lowest magnetic field strength (7 MG) of any known polar. We find evidence for structure in the X-ray eclipse profile which shows a `standstill' feature lasting 26+/-4 sec. This allows us to place an upper limit on the mass of the white dwarf of ~1.2 Msun. We find no evidence for QPOs in the frequency range 0.02-10 Hz. This coupled with the absence of QPOs in RXTE data suggest that, if present, any oscillations in the shock front have a minimal effect on the resultant X-ray flux. We find no evidence for a distinct soft X-ray component in its spectrum - it therefore joins another 7 systems which do not show this component. We suggest that those systems which are asynchronous, have low mass transfer rates, or have accretion occurring over a relatively large fraction of the white dwarf are more likely to show this effect. We find that the specific mass transfer rate has to be close to 0.1 g cm^-2 s^-1 to predict masses which are consistent with that derived from our eclipse analysis. This maybe due to the fact that the low magnetic field strength allows accretion to take place along a wide range of azimuth.

We investigate the hydrodynamics of accretion channelled by a dipolar magnetic field (funnel flows). We consider situations in which the electrons and ions in the flow cannot maintain thermal equilibrium (two-temperature effects) due to strong radiative loss, and determine the effects on the keV X-ray properties of the systems. We apply this model to investigate the accretion shocks of white dwarfs in magnetic cataclysmic variables. We have found that the incorporation of two-temperature effects could harden the keV X-rays. Also, the dipolar model yields harder X-ray spectra than the standard planar model if white dwarf is sufficiently massive (>~1M_sun). When fitting observed keV X-ray spectra of magnetic cataclysmic variables, the inclusion of two-temperature hydrodynamics and a dipolar accretion geometry lowers estimates for white-dwarf masses when compared with masses inferred from models excluding these effects. We find mass reductions <~9% in the most massive cases.

We present a model for the prediction of the optical/infra-red emission from ULXs. In the model, ULXs are binary systems with accretion taking place through Roche lobe overflow. We show that irradiation effects and presence of an accretion disk significantly modify the optical/infrared flux compared to single stars, and also that the system orientation is important. We include additional constraints from the mass transfer rate to constrain the parameters of the donor star, and to a lesser extent the mass of the BH. We apply the model to fit photometric data for several ULX counterparts. We find that most donor stars are of spectral type B and are older and less massive than reported elsewhere, but that no late-type donors are admissable. The degeneracy of the acceptable parameter space will be significantly reduced with observations over a wider spectral range, and if time-resolved data become available.

We have constructed a model which describes the optical emission from ultraluminous X-ray sources (ULXs), and have used it to constrain the parameters of seven ULX systems. Our model assumes a binary nature for ULXs, and accounts for optical emission from an X-ray irradiated companion star and accretion disk. We apply our model to six different ULX optical counterparts observed with HST, and one observed with the ESO VLT, and determine the mass, radius and age of the donor stars in these systems. In addition, we obtained constraints for the black hole (BH) mass in some cases. We use the mass accretion rate implied by the X-ray luminosity of these sources as an additional constraint on the donor star, by assuming the mass transfer is driven by the stellar nuclear evolution. We find that in general the donors are older and less massive than previously thought, and are consistent with being of spectral type B. We discuss how these results affect our understanding of the evolution and history of ULXs. Where we can constrain the BH masses, we find them to be consistent with stellar mass BHs or intermediate mass BHs of order 100 solar masses. We make predictions for future observations of optical/infrared ULX counterparts, calculating binary periods for different BH masses in each of the seven sources.

We report the discovery with XMM-Newton of quasiperiodic variability in the 0.2 - 10 keV X-ray flux from the ultraluminous X-ray source NGC 5408 X-1. The average power spectrum of all EPIC-pn data reveals a strong 20 mHz QPO with an average amplitude (rms) of 9%, and a coherence of Q~6. In a 33 ksec time interval when the 20 mHz QPO is strongest we also find evidence for a 2nd QPO peak at 15 mHz. This is the first indication for a close pair of QPOs in a ULX source. Interestingly, the frequency ratio of this QPO pair is inconsistent with 3:2 at the 3sigma level, but is consistent with a 4:3 ratio. A powerlaw noise component with slope near 1.5 is also present below 0.1 Hz with evidence for a break to a flatter slope at about 3 mHz. The source shows substantial broadband variability, with a total amplitude (rms) of about 30% in the 0.1 - 100 mHz frequency band, and there is strong energy dependence to the variability. We discuss the implications of these findings in the context of models for ULXs, and their implications for the object's mass.

We present time-resolved spectroscopy of the optical counterpart to the proposed ultra-compact binary system V407 Vul (=RX J1914.4+2456). Our Gemini spectra resolve the 9.48 minute periodicity that has previously been reported for this source. We find that the optical counterpart is dominated by a reddened late-type spectrum of type G9V and contains solely unresolved absorption features. No radial velocity signatures exceeding 10 km/s could be detected on periods from minutes to hours. Using interstellar extinction estimates, we derive a distance to the G9 star of 1.1-3.5 kpc. In addition to this stellar spectrum, we detect a blue component that modulates solely on the 9.48 minute period, and peaks ~0.15 in phase ahead of the X-ray peak. This blue component which contributes up to 40% of the light shows no evidence for emission line features that are the usual hallmarks of an interacting binary. Good seeing images obtained with the Magellan telescopes indicate that the variable and the G-star are aligned to better than 0.1". Despite the low probability of a chance alignment of a field star along the line of sight, the G9 light cannot be directly associated with the 9.48 minute variable that powers the luminous (~10^35 ergs/s) and highly variable X-ray source. The outlook for the detection of conclusive radial velocity measurements will remain challenging due to the extinction along the line of sight in conjunction with the contaminating effect of the G9V star.

We present XMM-Newton and Chandra X-ray observations of the middle-aged radio pulsar PSR B0355+54. Our X-ray observations reveal emission not only from the pulsar itself, but also from a compact diffuse component extending ~50'' in the opposite direction to the pulsar's proper motion. There is also evidence for the presence of fainter diffuse emission extending ~5' from the point source. The compact diffuse feature is well-fitted with a power-law, the index of which is consistent with the values found for other pulsar wind nebulae. The morphology of the diffuse component is similar to the ram-pressure confined pulsar wind nebulae detected for other sources. The X-ray emission from the pulsar itself is described well by a thermal plus power-law fit, with the thermal emission most likely originating in a hot polar cap.

The nature of the X-ray source RX J1914+24 has been the subject of much debate. It shows a prominent period of 569 sec in X-rays and the optical/infra-red: in most models this has been interpreted as the binary orbital period. We present our analysis of new XMM-Newton and Chandra data. We find a longer term trend in the XMM-Newton data and power at 556 and 585 sec in 5 sets of data. It is not clear if they are produced as a result of a beat between a longer intrinsic period and the 569 sec modulation or if they are due to secular variations. We obtain a good fit to the XMM-Newton spectrum with a low temperature thermal plasma model with an edge at 0.83keV. This model implies an unabsorbed bolometric X-ray luminosity of 1x10^{33} ergs/s (for a distance of 1kpc) - this is 2 orders of magnitude lower than our previous estimate (derived using a different model). If the distance is much less, as the absorption derived from the X-ray fits suggest, then it is even lower at \~3x10^{31} ergs/s.

We report on the first XMM-Newton observation of the Vela-like pulsar PSR B2334+61. Spectral analysis reveals soft X-ray emission, with the bulk of the photons emitted at energies below ~1.5 keV. We find that the spectrum has a thermal origin and is well-fitted with either a blackbody or a magnetized, pure H atmospheric model. In the latter case, for a neutron star with a radius of 13 km and a magnetic field of 10e13 G, the best-fit gives an hydrogen column density nH = 0.33 x 10^22 cm^-2 and an effective temperature T_eff^infinity = 0.65 x 10^6 K, as measured at Earth. A comparison of the surface temperature of PSR B2334+61 obtained from this fit with cooling curves favor a medium mass neutron star with M ~ 1.45 solar masses or M ~ 1.6 solar masses, depending on two different models of proton superfluidity in the interior. We do not detect any pulsed emission from the source, and determine an upper limit of 5% for the modulation amplitude of the emission on the pulsar's radio frequency.

The Gaia mission is designed as a Galaxy explorer, and will measure simultaneously, in a survey mode, the five or six phase space parameters of all stars brighter than 20th magnitude, as well as providing a description of their astrophysical characteristics. These measurements are obtained by combining an astrometric instrument with micro-arcsecond capabilities, a photometric system giving the magnitudes and colours in 15 bands and a medium resolution spectrograph named the Radial Velocity Spectrometer (RVS). The latter instrument will produce spectra in the 848 to 874 nm wavelength range, with a resolving power R = 11 500, from which radial velocities, rotational velocities, atmospheric parameters and abundances can be derived. A companion paper (Katz et al. 2004) presents the characteristics of the RVS and its performance. This paper details the outstanding scientific impact of this important part of the Gaia satellite on some key open questions in present day astrophysics. The unbiased and simultaneous acquisition of multi-epoch radial velocities and individual abundances of key elements in parallel with the astrometric parameters is essential for the determination of the dynamical state and formation history of our Galaxy. Moreover, for stars brighter than V=15, the resolving power of the RVS will give information about most of the effects which influence the position of a star in the Hertzsprung-Russell diagram, placing unprecedented constraints on the age, internal structure and evolution of stars of all types. Finally, the RVS multi-epoch observations are ideally suited to the identification, classification and characterisation of the many types of double, multiple and variable stars.

We have constructed a model to describe the optical emission from ultra-luminous X-ray sources (ULXs). We assume a binary model with a black hole accreting matter from a Roche lobe filling companion star. We consider the effects of radiative transport and radiative equilibrium in the irradiated surfaces of both the star and a thin accretion disk. We have developed this model as a tool with which to positively identify the optical counterparts of ULXs, and subsequently derive parameters such as the black hole mass and the luminosity class and spectral type of the counterpart. We examine the dependence of the optical emission on these and other variables. We extend our model to examine the magnitude variation at infrared wavelengths, and we find that observations at these wavelengths may have more diagnostic power than in the optical. We apply our model to existing HST observations of the candidates for the optical counterpart of ULX X-7 in NGC 4559. All candidates could be consistent with an irradiated star alone, but we find that a number of them are too faint to fit with an irradiated star and disk together. Were one of these the optical counterpart to X-7, it would display a significant temporal variation.

We present the results of XMM-Newton observations of four AM CVn systems -- AM CVn, CR Boo, HP Lib and GP Com. Their light curves show very different characteristics. The X-ray light curves show no coherent pulsations, suggesting the accreting white dwarfs have relatively low magnetic field strengths. Their spectra were best modelled using a multi-temperature emission model and a strong UV component. We find that CR Boo and HP Lib have X-ray spectra with abundances consistent with relatively low temperature CNO processed material, while AM CVn and GP Com show an enhancement of nitrogen. A large fraction of the accretion luminosity is emitted in the UV. We determine accretion luminosities of ~1.6x10^{33} ergs/s and 1.7x10^{31} ergs/s for AM CVn and GP Com respectively. Comparing the implied mass transfer rates with that derived using model fits to optical and UV spectra, we find evidence that in the case of AM CVn, we do not detect a significant proportion of the accretion energy. This missing component could be lost in the form of a wind.

We use a two-temperature hydrodynamical formulation to determine the temperature and density structures of the post-shock accretion flows in magnetic cataclysmic variables (mCVs) and calculate the corresponding X-ray spectra. The effects of two-temperature flows are significant for systems with a massive white dwarf and a strong white-dwarf magnetic field. Our calculations show that two-temperature flows predict harder keV spectra than one-temperature flows for the same white-dwarf mass and magnetic field. This result is insensitive to whether the electrons and ions have equal temperature at the shock but depends on the electron-ion exchange rate, relative to the rate of radiative loss along the flow. White-dwarf masses obtained by fitting the X-ray spectra of mCVs using hydrodynamic models including the two-temperature effects will be lower than those obtained using single-temperature models. The bias is more severe for systems with a massive white dwarf.

Field-channelled accretion flows occur in a variety of astrophysical objects, including T Tauri stars,magnetic cataclysmic variables and X-ray pulsars. We consider a curvilinear coordinate system and derive a general hydrodynamic formulation for accretion onto stellar objects confined by a stellar dipole magnetic field. The hydrodynamic equations are solved to determine the velocity, density and temperature profiles of the flow. We use accreting magnetic white-dwarf stars as an illustrative example of astrophysical applications. Our calculations show that the compressional heating due to the field geometry is as important as radiative cooling and gravity in determining the structure of the post-shock flow in accreting white-dwarf stars. The generalisation of the formulation to accretion flows channelled by higher-order fields and the applications to other astrophysical systems are discussed.

We report on the results of a deep XMM-Newton observation of RBS 1774, the most recent dim isolated neutron star candidate found in the ROSAT archive data. Spectral and timing analysis of the high-quality PN and MOS data confirm the association of this source with an isolated neutron star. The spectrum is thermal and blackbody-like, and there is evidence at a significance level > 4sigma that the source is an X-ray pulsar, with spin period of 9.437 s. Spectral fitting reveils the presence of an absorption feature at ~0.7 keV, but at this level data do not have enough resolution to allow us to discriminate between an absorption line or an edge. We compare the newly measured properties of RBS 1774 with those of other known dim isolated neutron stars, and discuss possible interpretations for the absorption feature.

We present high time-resolution optical energy-resolved photometry of the eclipsing cataclysmic variable V2301 Oph made using the ESA S-Cam detector, an array of photon counting super-conducting tunnel junction (STJ) devices with intrinsic energy resolution. Three eclipses were observed, revealing considerable variation in the eclipse shape, particularly during ingress. The eclipse shape is shown to be understood in terms of AM Her accretion via a bright stream, with very little contribution from the white dwarf photosphere and/or hotspot. About two thirds of the eclipsed light arises in the threading region. Variation in the extent of the threading region can account for most of the variations observed between cycles. Spectral fits to the data reveal a 10,000K blackbody continuum with strong, time-varying emission lines of hydrogen and helium. This is the first time that stellar emission lines have been detected in the optical band using a non-dispersive photon-counting system.

We present XMM-Newton observations of the 569 sec period system RX J1914+24 (V407 Vul). This period is believed to represent the binary orbital period making it an ultra-compact binary system. By comparing the phase of the rise to maximum X-ray flux at various epochs (this includes observations made using ROSAT, ASCA and Chandra) we find that the system is spinning up at a rate of 3.17+/-0.07x10^{-12} s/s. We find that the spectra softens as the X-ray flux declines towards the off-phase of the 569 sec period. Further, the spectra are best fitted by an absorbed blackbody component together with a broad emission feature around 0.59keV. This emission feature is most prominent at the peak of the on-phase. We speculate on its origin.

The definition and optimisation studies for the Gaia satellite spectrograph, the Radial Velocity Spectrometer (RVS), converged in late 2002 with the adoption of the instrument baseline. This paper reviews the characteristics of the selected configuration and presents its expected performance. The RVS is a 2.0 by 1.6 degree integral field spectrograph, dispersing the light of all sources entering its field of view with a resolving power R=11 500 over the wavelength range [848, 874] nm. The RVS will continuously and repeatedly scan the sky during the 5 years of the Gaia mission. On average, each source will be observed 102 times over this period. The RVS will collect the spectra of about 100-150 million stars up to magnitude V~17-18. At the end of the mission, the RVS will provide radial velocities with precisions of ~2 km/s at V=15 and \~15-20 km/s at V=17, for a solar metallicity G5 dwarf. The RVS will also provide rotational velocities, with precisions (at the end of the mission) for late type stars of sigma_vsini ~5 km/s at V~15 as well as atmospheric parameters up to V~14-15. The individual abundances of elements such as Silicon and Magnesium, vital for the understanding of Galactic evolution, will be obtained up to V~12-13. Finally, the presence of the 862.0 nm Diffuse Interstellar Band (DIB) in the RVS wavelength range will make it possible to derive the three dimensional structure of the interstellar reddening.

We have studied the candidate optical counterparts and the stellar population in the star-forming complex around a bright ULX in NGC4559, using HST/WFPC2, XMM-Newton/OM, and ground-based data. We find that the ULX is located near a small group of OB stars. The brightest point source in the Chandra error circle is consistent with a single blue supergiant of mass ~ 20 M_sun and age ~ 10 Myr. A few other stars are resolved inside the error circle: mostly blue and red supergiants with masses ~ 10-15 M_sun and ages ~ 20 Myr. This is consistent with the interpretation of this ULX as a black hole (BH) accreting from a high-mass donor star in its supergiant phase, via Roche-lobe overflow. The observed optical colors and the blue-to-red supergiant ratio suggest a low metal abundance: 0.2 <~ Z/Z_sun <~ 0.4 (Padua tracks), or 0.05 <~ Z/Z_sun <~ 0.2 (Geneva tracks). The age of the star-forming complex is <~ 30 Myr. H-alpha images show that this region has a ring-like appearance. We propose that it is an expanding wave of star formation, triggered by an initial density perturbation, in a region where the gas was only marginally stable to gravitational collapse. A possible trigger was the collision with a satellite dwarf galaxy, visible a few arcsec north-west of the complex, going through the gas-rich outer disk of NGC4559. The X-ray data favour a BH more massive (M > 50 M_sun) than typical Milky Way BH candidates. The optical data favour a ``young'' BH originating in the recent episode of massive star formation; however, they also rule out an association with young massive star clusters. We speculate that other mechanisms may lead to the formation of relatively massive BHs (~ 50-100 M_sun) from stellar evolution processes in low-metallicity environments, or when star formation is triggered by galactic collisions.

Ultraluminous X-ray sources (ULXs) are likely to include different physical types of objects. We discuss some possible subclasses, reviewing the properties of a sample of ULXs recently observed by Chandra and XMM-Newton. Sources with an isotropic X-ray luminosity up to a few times 10^{39} erg/s are consistent with ``normal'' stellar-mass X-ray binaries (mostly high-mass X-ray binaries in star-forming regions). Higher black hole (BH) masses (~ 50-100 M_sun) may be the end product of massive stellar evolution in peculiar environments: they may explain ULXs with luminosities ~ 1-2 x 10^{40} erg/s. Only a handful of ULXs require a true intermediate-mass BH (M >~ 500 M_sun). Finally, a small subclass of ULXs shows flaring or rapid variability in its power-law spectral component.

We present XMM-Newton observations of the eclipsing polar EP Dra which cover nearly 3 binary orbital cycles. The X-ray and UV data show evidence for a prominent dip before the eclipse which is due to the accretion stream obscuring the accretion region. The dip ingress is rapid in hard X-rays suggesting there is a highly collimated core of absorption. We find that a different level of absorption column density is required to match the observed count rates in different energy bands. We propose that this is due to the fact that different absorption components should be used to model the reprocessed X-rays, the shocked X-ray component and the UV emission and explore the affect that this has on the resulting fits to the spectrum. Further, there is evidence that absorption starts to obscure the softer X-rays shortly after the onset of the bright phase. This suggests that material is threaded by an unusually wide range of magnetic field lines, consistent with the suggestion of Bridge et al. We find that the period is slightly greater than that determined by Schwope & Mengel.

We report on the analysis of new X-ray data obtained with XMM-Newton and Chandra from two ROSAT-discovered X-ray dim isolated neutron stars. RX J0806.4-4123 was observed with XMM-Newton in April 2003, 2.5 years after the first observation. The EPIC-pn data confirm that this object is an X-ray pulsar with 11.371 s neutron star spin period. The X-ray spectrum is consistent with absorbed black-body emission with a temperature kT = 96 eV and NH = 4x10^19 cm^-2 without significant changes between the two observations. Four XMM-Newton observations of RX J0420.0-5022 between December 2002 and July 2003 did not confirm the 22.7 s pulsations originally indicated in ROSAT data, but clearly reveal a 3.453 s period. A fit to the X-ray spectrum using an absorbed black-body model yields kT = 45 eV, the lowest value found from the small group of XDINs and NH = 1.0x10^20 cm^-2. Including a broad absorption line improves the quality of the spectral fits considerably for both objects and may indicate the presence of absorption features similar to those reported from RBS1223, RX J1605.3+3249 and RX J0720.4-3125. For both targets we derive accurate X-ray positions from the Chandra data and present an optical counterpart candidate for RX J0420.0-5022 with B = 26.6+/-0.3 mag from VLT imaging.

We apply an eclipse mapping technique using `genetically modified fire-flies' to the eclipse light curves of HU Aqr and EP Dra. The technique makes as few assumptions as possible about the location of accretion stream material, allowing the emission to be located anywhere within the Roche lobe of the white dwarf. We model two consecutive eclipses in the UBVR_c band for HU Aqr, and four consecutive `white'-light eclipses for EP Dra, to determine the changing brightness distribution of stream material. We find fire-fly distributions which are consistent with accretion through a curtain of material in both HU Aqr and EP Dra, and show that the previously assumed two part ballistic and magnetic trajectory is a good approximation for polars. Model fits to the colour band data of HU Aqr indicate that the material confined to the magnetic field lines is brightest, and most of the emission originates from close to the white dwarf. There is evidence for emission from close to a calculated ballistic stream in both HU Aqr and EP Dra.We propose that a change in the stream density causes a change in the location of the bright material in the accretion stream in EP Dra.

We present a reanalysis of the X-ray data for RX J0720.4-3125 presented in Zane et al. (2002), using more data recently available from XMM-Newton and Chandra. This analysis also corrects the Rosat data used in that paper to the TDB time system, incorporates the revised XMM-Newton barycentric correction available since then, and corrects the definition of the instantaneous period in the maximum likelihood periodogramme search. However, we are now unable to find a single coherent period that is consistent with all Rosat, Chandra and XMM-Newton datasets. From an analysis of the separate datasets, we have derived limits on the period change of Pdot = 1.4+/-0.6x10^-13 s/s at 99% confidence level. This is stronger than the value presented in Zane et al. (2002), but sufficiently similar that their scientific conclusions remain unchanged. We examine the implications in more detail, and find that RX J0720.4-3125 can have been born as a magnetar provided that it has a young age of ~10^4 yr. A more conservative interpretation is that the field strength has remained relatively unchanged at just over 10^13 G, over the ~10^6 yr lifetime of the star.

We have made a series of snap-shot observations of 37 polars using XMM-Newton. We found that 16 of these systems were in a low, or much reduced, accretion state. Of those, 6 were not detected in X-rays. This suggests that in any survey of polars, around half will be in a low accretion state. We tested if there was a bias towards certain orbital periods: this is not the case. Of the 10 systems which were detected at low, but significant rates in X-rays, 8 showed significant variability in their X-ray light curves. This implies that non-uniform accretion still takes place during low accretion epochs. The bolometric luminosity of these systems is ~10^30 ergs s, two orders of magnitude less than for systems in a high accretion state. The X-ray spectra show no evidence of a distinct soft X-ray component. However, the X-ray and UV data imply that such a low temperature component exists: its temperature is low enough for its flux distribution to move outside the bandpass of the X-ray instruments.

We have studied an ultraluminous X-ray source (ULX) in NGC 4559 with XMM-Newton, and its peculiar star-forming environment with HST/WFPC2. The X-ray source is one of the brightest in its class (L_x ~ 2 x 10^{40} erg/s). Luminosity and timing arguments suggest a mass >~ 50 M_sun for the accreting black hole. The ULX is located near the rim of a young (age < 30 Myr), large (diameter ~ 700 pc) ring-like star forming complex possibly triggered by the impact of a dwarf satellite galaxy through the gas-rich outer disk of NGC 4559. We speculate that galaxy interactions (including the infall of high-velocity clouds and satellites on a galactic disk) and low-metallicity environments offer favourable conditions for the formation of compact remnants more massive than "standard" X-ray binaries, and accreting from a massive Roche-lobe filling companion.

We report on a ~37-ks XMM-Newton observation of the intermediate polar FO Aquarii, presenting X-ray and UV data from the EPIC and OM cameras. We find that the system has changed from its previously reported state of disc-overflow accretion to one of purely disc-fed accretion. We detect the previously reported `notch' feature in the X-ray spin pulse, and explain it as a partial occultation of the upper accretion pole. Maximum flux of the quasi-sinusoidal UV pulse coincides with the notch, in keeping with this idea. However, an absorption dip owing to the outer accretion curtains occurs 0.27 later than the expected phase, which implies that the accretion curtains are twisted, trailing the magnetic poles. This result is the opposite of that reported in PQ Gem, where accreting field lines were found to lead the pole. We discuss how such twists relate to the accretion torques and thus the observed period changes of the white dwarfs, but find no simple connection.

We have examined X-ray and optical observations of two ultra-luminous X-ray sources, X7 and X10 in NGC4559, using XMM-Newton, Chandra and HST. The UV/X-ray luminosity of X7 exceeds 2.1e40 erg/s in the XMM-Newton observation, and that of X10 is >1.3e40 erg/s. X7 has both thermal and power-law spectral components. The characteristic temperature of the thermal component is 0.12 keV. The power-law components in the two sources both have slopes with photon index \~2.1. A timing analysis of X7 indicates a break frequency at 28 mHz in the power spectrum, while that for X10 is consistent with an unbroken power law. The luminosity of the blackbody component in the X-ray spectrum of X7 and the nature of its time variability provides evidence that this object is an intermediate mass black hole accreting at sub-Eddington rates, but other scenarios which require high advection efficiencies from a hollowed-out disk might be possible. The emission from X10 can be characterised by a single power-law. This source can be interpreted either as an intermediate mass black hole, or as a stellar-mass black hole with relativistically-beamed Comptonised emission. There are four optical counterparts in the error circle of X7. No counterparts are evident in the error circle for X10.

In the EXOSAT and ROSAT eras a significant number of polars were found to show a soft/hard X-ray ratio much greater than that expected from the standard accretion shock model. This was known as the `soft X-ray excess'. We have made an snapshot survey of polars using XMM-Newton and determined their soft/hard ratios. We find that less than one in five of systems show a significant soft X-ray excess, while the rest show ratios consistent with that predicted by the standard model. We have investigated the discrepancy between this and the previous investigations by re-examining all the available ROSAT PSPC pointed observations of polars using more recent calibrations than in the original studies. We find that these data show an energy balance ratio which is broadly consistent with that of our XMM-Newton results. We conclude that the previous studies were affected by the data being less well calibrated. We discuss which physical mechanisms might give rise to a high soft X-ray excess and whether systems with high ratios show more variation in soft X-rays. Surprisingly, we find that 6 out of 21 systems found in a high accretion state did not show a distinct soft X-ray component. Two systems showed one pole with such a component and one which did not. Based on the ratio of the observed soft X-ray to UV flux measurements (which were obtained simultaneously using the Optical Monitor) we suggest that this is because the reprocessed component in these systems is cool enough to have moved out of the soft X-ray band and into the EUV or UV band.

We report on the XMM-Newton observations of the young, 102 ms pulsar PSR B1706-44. We have found that both a blackbody plus power-law and a magnetized atmospheric model plus power-law provide an excellent fit to the EPIC spectra. The two scenarios are therefore indistinguishable on a statistical basis, although we are inclined to prefer the latter on physical grounds. In this case, assuming a source distance of ~2.3 kpc, the size of the region responsible for the thermal emission is R~13 km, compatible with the surface of a neutron star. A comparison of the surface temperature of PSR B1706-44 obtained from this fit with cooling curves favor a medium mass neutron star with M~1.45 solar masses or M~1.59 solar masses, depending on two different models of proton superfluidity in the interior. The large collecting area of XMM-Newton allows us to resolve a substructure in the broad soft X-ray modulation detected by Chandra, revealing the presence of two separate peaks with pulsed fractions of 7 +/- 4% and 15 +/- 3%, respectively.

We present phase-resolved XMM_Newton data of three short period polars: V347 Pav, GG Leo and EU UMa. All three systems show one dominant accretion region which is seen for approximately half of the orbital cycle. GG Leo shows a strong dip feature in its X-ray and UV light curves which is due to absorption of X-rays from the accretion site by the accretion stream. The emission in the case of EU UMa is dominated by soft X-rays: its soft/hard X-ray ratio is amongst the highest seen in these objects. In contrast, GG Leo and V347 Pav shows a ratio consistent with that predicted by the standard shock model. We infer the mass of the white dwarf and explore the affect of restricting the energy range on the derived parameters.

RX J0806+15 has recently been identified as the binary system with the shortest known orbital period. We present a series of observations of RX J0806+15 including new optical observations taken one month apart. Using these observations and archival data we find that the period of this system is decreasing over time. Our measurements imply f_dot = 6.11x10^-16 Hz/s, which is in agreement with a rate expected from the gravitational radiation for two white dwarfs orbiting at a given period. However, a smaller value of f_dot = 3.14x10^-16 Hz/s cannot be ruled out. Our result supports the idea that the 321.5 s period is the orbital period and that the system is the shortest period binary known so far and that it is one of the strongest sources of constant gravitational radiation in the sky. Furthermore, the decrease of the period strongly favours the unipolar inductor (or electric star) model rather than the accretion models.

We describe a multi-order spectrograph concept suitable for 8m-class telescopes, using the intrinsic spectral resolution of Superconducting Tunneling Junction detectors to sort the spectral orders. The spectrograph works at low orders, 1-5 or 1-6, and provides spectral coverage with a resolving power of R~8000 from the atmospheric cutoff at 320 nm to the long wavelength end of the infrared H or K band at 1800 nm or 2400 nm. We calculate that the spectrograph would provide substantial throughput and wavelength coverage, together with high time resolution and sufficient dynamic range. The concept uses currently available technology, or technologies with short development horizons, restricting the spatial sampling to two linear arrays; however an upgrade path to provide more spatial sampling is identified. All of the other challenging aspects of the concept - the cryogenics, thermal baffling and magnetic field biasing - are identified as being feasible.

We report on the XMM-Newton observations of the 143 ms pulsar PSR J0538+2817. We present evidence for the first detections of pulsed X-rays from the source at a frequency which is consistent with the predicted radio frequency. The pulse profile is broad and asymmetric, with a pulse fraction of 18 +/- 3%. We find that the spectrum of the source is well-fit with a blackbody with T^{infty} = (2.12^{+0.04}_{-0.03}) x 10^6 K and N_{H} = 2.5 x 10^21 cm^{-2}. The radius determined from the model fit of 1.68 +/- 0.05 km suggests that the emission is from a heated polar cap. A fit to the spectrum with an atmospheric model reduces the inferred temperature and hence increases the radius of the emitting region, however the pulsar distance determined from the fit is then smaller than the dispersion distance.

We present phase resolved optical spectroscopy and X-ray light curves of the asynchronous polar RX J2115-58 as they change on the beat period.

We apply a new method of eclipse mapping to the light curves of eclipsing polars. The technique aims to locate the bright emission associated with the accretion stream, using a technique that makes the fewest prior assumptions about the location of the accretion stream material. We have obtained data of EP Dra and HU Aqr with the S-Cam 2 superconducting tunnel junction camera using the William Herschel Telescope. The location of emission regions in both systems show that previously assumed trajectories are consistent with those found using our technique. Most of the emission is located in a region where we expect material to be confined to magnetic field lines, particularly for HU Aqr, while there appears to be less emission from where we conventionally expect material to follow a ballistic trajectory from the L1-point.

We review the current observational status of the ROSAT sources RX J1914.4+2456 and RX J0806.3+1527, and the evidence that these are ultra-short period (<10min) binary systems. We argue that an Intermediate Polar interpretation can be ruled out, that they are indeed compact binaries with a degenerate secondary, and that the period seen in the X-ray and optical is the orbital period. A white dwarf primary is preferred, but a neutron star cannot be excluded. We examine the capability of the three current double-degenerate models (Polar, Direct Accretor and Electric Star) to account for the observational characteristics of these systems. All models have difficulties with some aspects of the observations, but none can be excluded with confidence at present. The Electric Star model provides the best description, but the lifetime of this phase requires further investigation. These ultra-short period binaries will be strong gravitational wave emitters in the LISA bandpass, and because of their known source properties will be important early targets for gravitational wave studies.

A hydrodynamic formulation for accretion flow channeled by a dipolar magnetic field is constructed using a curvi-linear coordinate system natural to the field structure. We solve the hydrodynamic equations and determine the velocity, density and temperature profiles of the post-shock accretion flow. The results are applied to accretion flows in intermediate polars. We have found that for systems with massive white dwarfs (~1Msolar) the temperature profiles in the flow can differ significantly to those obtained from models in which the accretion column is assumed to be cylindrical.

The XMM-Newton-MSSL survey of polars is a snapshot survey of nearly 40 polars. We present the initial results, with an emphasis on the energy balance in these systems and the number of systems which were observed in a low accretion state. In contrast to the ROSAT results, we find that the vast majority of polars show no `soft X-ray excess'. We find that half of our sample were observed in low or off accretion states. We find no evidence that these systems were biased towards certain orbital periods.

We present the first high time resolution light curves for six eclipses of the magnetic cataclysmic variable EP Dra, taken using the superconducting tunnel junction imager S-Cam2. The system shows a varying eclipse profile between consecutive eclipses over the two nights of observation. We attribute the variable stream eclipse after accretion region ingress to a variation in the amount and location of bright material in the accretion stream. This material creates an accretion curtain as it is threaded by many field lines along the accretion stream trajectory. We identify this as the cause of absorption evident in the light curves when the system is in a high accretion state. We do not see direct evidence in the light curves for an accretion spot on the white dwarf; however, the variation of the stream brightness with the brightness of the rapid decline in flux at eclipse ingress indicates the presence of some form of accretion region. This accretion region is most likely located at high colatitude on the white dwarf surface, forming an arc shape at the foot points of the many field lines channeling the accretion curtain.

We report simultaneous UBVRI photo-polarimetric observations of the long period (7.98 h) AM Her binary V1309 Ori. The length and shape of the eclipse ingress and egress varies from night to night. We suggest this is due to the variation in the brightness of the accretion stream. By comparing the phases of circular polarization zero-crossovers with previous observations, we confirm that V1309 Ori is well synchronized, and find an upper limit of 0.002 percent for the difference between the spin and orbital periods. We model the polarimetry data using a model consisting of two cyclotron emission regions at almost diametrically opposite locations, and centered at colatitude 35 (deg) and 145 (deg) on the surface of the white dwarf. We also present archive X-ray observations which show that the negatively polarised accretion region is X-ray bright.

We present XMM-Newton data of three strongly magnetic cataclysmic variables (polars) EV UMa, RX J1002-19 and RX J1007-20. These include the polar with the shortest orbital period (EV UMa) and the polar with one of the highest magnetic field strengths (RX J1007-20). They exhibit a range of X-ray spectral characteristics which are consistent with their known magnetic field strength. We find that two of the systems show evidence for an absorption dip in soft X-rays. Their profiles are well defined, implying that the stream is highly collimated. We determine the mass transfer rate for the two systems with known distances. We determine that the mass of the white dwarf in EV UMa and RX J1007-20 is ~1Msun while in RX J1002-19 it is closer to ~0.5Msun.

We apply an eclipse mapping technique to observations of the eclipsing magnetic cataclysmic variable HU Aqr. The observations were made with the S-Cam2 Superconducting Tunnel Junction detector at the WHT in October 2000, providing high signal-to-noise observations with simultaneous spectral and temporal resolution. HU Aqr was in a bright (high accretion) state (V=14.7) and the stream contributes as much to the overall system brightness as the accretion region on the white dwarf. The stream is modelled assuming accretion is occuring onto only one pole of the white dwarf. We find enhanced brightness towards the accretion region from irradiation and interpret enhanced brightness in the threading region, where the ballistic stream is redirected to follow the magnetic field lines of the white dwarf, as magnetic heating from the stream-field interaction, which is consistent with recent theoretical results. Changes in the stream eclipse profile over one orbital period indicate that the magnetic heating process is unstable.

Using RXTE and ASCA data, we investigate the roles played by occultation and absorption in the X-ray spin pulse profile of the Intermediate Polar PQ Gem. From the X-ray light curves and phase-resolved spectroscopy, we find that the intensity variations are due to a combination of varying degrees of absorption and the accretion regions rotating behind the visible face of the white dwarf. These occultation and absorption effects are consistent with those expected from the accretion structures calculated from optical polarisation data. We can reproduce the changes in absorber covering fraction either from geometrical effects, or by considering that the material in the leading edge of the accretion curtain is more finely fragmented than in other parts of the curtain. We determine a white dwarf mass of ~1.2 using the RXTE data.

We report the detection of the polar CE Gru in X-rays for the first time. We find evidence for a dip seen in the hard X-ray light curve which we attribute to the accretion stream obscuring the accretion region in the lower hemisphere of the white dwarf. The X-ray spectrum can be fitted using only a shock model: there is no distinct soft X-ray component. We suggest that this is because the reprocessed component is cool enough so that it is shifted into the UV. We determine a mass for the white dwarf of ~1.0Msun.

We report observations of the asynchronous polar BY Cam made using XMM-Newton. We find evidence for two accretion regions which have significantly different spectra. In both regions we find evidence for hard X-ray emission from the post-shock flow while we observe a distinct soft X-ray blackbody component in only one. We detect two emission lines in the RGS detectors which we attribute to Nitrogen NVI and NVII: the first time that Nitrogen lines have been detected in an X-ray spectrum of a cataclysmic variable. In the time interval where we observe short timescale variability we find the hard X-ray light curve is correlated with the softest energies in the sense that the hard X-rays trail the soft X-rays and are anti-correlated. We suggest that there are a small number of dense blobs of material which impact the white dwarf without forming a shock and release an amount of optically thick material which obscures hard X-rays for a short duration.

We present a combined analysis of XMM-Newton, Chandra and Rosat observations of the isolated neutron star RXJ0720.4-3125, spanning a total period of \sim 7 years. We develop a maximum likelihood periodogramme for our analysis based on the \Delta C-statistic and the maximum likelihood method, which are appropriate for the treatment of sparse event lists. Our results have been checked "a posteriori" by folding a further BeppoSAX dataset with the period predicted at the time of that observation: the phase is found to be consistent. The study of the spin history and the measure of the spin-down rate is of extreme importance in discriminating between the possible mechanisms suggested for the nature of the X-ray emission. The value of \dot P, here measured for the first time, is \approx 10^{-14} s/s. This value can not be explained in terms of torque from a fossil disk. When interpreted in terms of dipolar losses, it gives a magnetic field of B \approx 10^{13} G, making also implausible that the source is accreting from the underdense surroundings. On the other hand, we also find unlikely that the field decayed from a much larger value (B\approx 10^{15} G, as expected for a magnetar powered by dissipation of a superstrong field) since this scenario predicts a source age of \approx 10^4 yrs, too young to match the observed X-ray luminosity. The observed properties are more compatible with a scenario in which the source is \approx 10^6 yrs old, and its magnetic field has not changed substantially over the lifetime.

We present a combined analysis of XMM-Newton, Chandra and Rosat observations of the isolated neutron star RX J0720.4-3125, spanning a total period of \sim 7 years. We develop a maximum likelihood periodogramme based on \Delta C statistic and maximum likelihood method, which are appropriate for sparse event lists. As an "a posteriori" check, we have folded a further BeppoSAX dataset with the period predicted at the time of that observation, finding that the phase is consistent. The value of the spin down rate, here measured for the first time, is \approx 10^{-14} s/s and can not be explained in terms of propeller or torque from a fossil disk. When interpreted in terms of dipolar losses, it gives a magnetic field of B \approx 10^{13} G, making also implausible that the source is accreting from the underdense surroundings. We discuss the implications of this measure for the different mechanisms that have been suggested to explain the X-ray emission. We conclude that the observed properties are more compatible with a scenario in which the source is \approx 10^6 yrs old, and its magnetic field has not changed substantially over the lifetime.

The X-ray source RX J0806+15 was discovered using ROSAT, and shows an X-ray light curve with a prominent modulation on a period of 321.5 sec. We present optical observations in which we report the detection of its optical counterpart. We find an optical period consistent with the X-ray period. We do not find convincing evidence for a second period in the data: this implies the 321.5 sec period is the orbital period. As such it would be the shortest period stellar binary system yet known. We discuss the nature of this system. We conclude that an isolated neutron star and an intermediate polar interpretation is unlikely and that a double degenerate interpretation is the most likely.

We present observations of the proposed double degenerate polar RX J1914+24. Our optical and infrared spectra show no emission lines. This, coupled with the lack of significant levels of polarisation provide difficulties for a double degenerate polar interpretation. Although we still regard the double degenerate polar model as feasible, we have explored alternative scenarios for RX~J1914+24. These include a double degenerate algol system, a neutron star-white dwarf pair and an electrically powered system. The latter model is particularly attractive since it naturally accounts for the lack of both emission lines and detectable polarisation in RX J1914+24. The observed X-ray luminosity is consistent with the predicted power output. If true, then RX J1914+24 would be the first known stellar binary system radiating largely by electrical energy.

XMM-Newton was used to observe two eclipsing, magnetic cataclysmic variables, DP Leo and WW Hor, continuously for three orbital cycles each. Both systems were in an intermediate state of accretion. For WW Hor we also obtained optical light curves with the XMM-Newton Optical Monitor and from ground-based observations. Our analysis of the X-ray and optical light curves allows us to constrain physical and geometrical parameters of the accretion regions and derive orbital parameters and eclipse ephemerides of the systems. For WW Hor we directly measure horizontal and vertical temperature variations in the accretion column. From comparisons with previous observations we find that changes in the accretion spot longitude are correlated with the accretion rate. For DP Leo the shape of the hard X-ray light curve is not as expected for optically thin emission, showing the importance of optical depth effects in the post-shock region. We find that the spin period of the white dwarf is slightly shorter than the orbital period and that the orbital period is decreasing faster than expected for energy loss by gravitational radiation alone.

We present spectroscopic and high speed photometric data of the eclipsing polar V895 Cen. We find that the eclipsed component is consistent with it being the accretion regions on the white dwarf. This is in contrast to Stobie et al who concluded that the eclipsed component was not the white dwarf. Further, we find no evidence for an accretion disc in our data. From our Doppler tomography results, we find that the white dwarf has M geq 0.7 Msun. Our indirect imaging of the accretion stream suggests that the stream is brightest close to the white dwarf. When we observed V895 Cen in its highest accretion state, emission is concentrated along field lines leading to the upper pole. There is no evidence for enhanced emission at the magnetic coupling region.

We propose a model for stellar binary systems consisting of a magnetic and a non-magnetic white-dwarf pair which is powered principally by electrical energy. In our model the luminosity is caused by resistive heating of the stellar atmospheres due to induced currents driven within the binary. This process is reminiscent of the Jupiter-Io system, but greatly increased in power because of the larger companion and stronger magnetic field of the primary. Electrical power is an alternative stellar luminosity source, following on from nuclear fusion and accretion. We find that this source of heating is sufficient to account for the observed X-ray luminosity of the 9.5-min binary RX J1914+24, and provides an explanation for its puzzling characteristics.

We present an analysis of the X-ray spectra of two strongly magnetic cataclysmic variables, DP Leo and WW Hor, made using XMM-Newton. Both systems were in intermediate levels of accretion. Hard optically thin X-ray emission from the shocked accreting gas was detected from both systems, while a soft blackbody X-ray component from the heated surface was detected only in DP Leo. We suggest that the lack of a soft X-ray component in WW Hor is due to the fact that the accretion area is larger than in previous observations with a resulting lower temperature for the re-processed hard X-rays. Using a multi-temperature model of the post-shock flow, we estimate that the white dwarf in both systems has a mass greater than 1 Msun. The implications of this result are discussed. We demonstrate that the `soft X-ray excess' observed in many magnetic cataclysmic variables can be partially attributed to using an inappropriate model for the hard X-ray emission.

We have investigated the ionisation structure of the post-shock regions of magnetic cataclysmic variables using an analytic density and temperature structure model in which effects due to bremsstrahlung and cyclotron cooling are considered. We find that in the majority of the shock-heated region where H- and He-like lines of the heavy elements are emitted, the collisional-ionisation and corona-condition approximations are justified. We have calculated the line emissivity and ionisation profiles for Iron as a function of height within the post-shock flow. For low-mass white dwarfs, line emission takes place near the shock. For high-mass white dwarfs, most of the line emission takes place in regions well below the shock and hence it is less sensitive to the shock temperature. Thus, the line ratios are useful to determine the white-dwarf masses for the low-mass white dwarfs, but the method is less reliable when the white dwarfs are massive. Line spectra can, however, be used to map the hydrodynamic structure of the post-shock accretion flow.

Using a parameterised function for the mass loss at the base of the post-shock region, we have constructed a formulation for magnetically confined accretion flows which avoids singularities, such as the infinity in density, at the base associated with all previous formulations. With the further inclusion of a term allowing for the heat input into the base from the accreting white dwarf we are able also to obtain the hydrodynamic variables to match the conditions in the stellar atmosphere. (We do not, however, carry out a mutually consistent analysis for the match). Changes to the emitted X-ray spectra are negligible unless the thickness of mass leakage region at the base approaches or exceeds one percent of the height of the post-shock region. In this case the predicted spectra from higher-mass white dwarfs will be harder, and fits to X-ray data will predict lower white-dwarf masses than previous formulations.

We present high time-resolution optical photometry of the eclipsing binary UZ For using a superconducting tunnel junction (STJ) device, a photon-counting array detector with intrinsic energy resolution. Three eclipses of the $\sim$18 mag 126.5 min orbital binary were observed using a $6\times6$ array of Tantalum STJs at the 4.2-m William Herschel Telescope on La Palma. The detector presently provides individual photon arrival time accuracy to about 5 $\mu$s, and a wavelength resolution of about 60 nm at 500 nm, with each array element capable of counting up to $\sim$5000 photons s$^{-1}$. The data allow us to place accurate constraints on the accretion geometry from our time- and spectrally-resolved monitoring, especially of the eclipse ingress and egress. We find that there are two small accretion regions, located close to the poles of the white dwarf. The positions of these are accurately constrained, and show little movement from eclipse to eclipse, even over a number of years. The colour of the emission from the two regions appears similar, although their X-ray properties are known to be significantly different: we argue that the usual accretion shock may be absent at the non-X-ray emitting region, and instead the flow here interacts directly with the white dwarf surface; alternatively, a special grazing occultation of this region is required. There is no evidence for any quasi-periodic oscillations on time-scales of the order of seconds, consistent with relatively stable cyclotron cooling in each accretion region.

The XMM-OM instrument extends the spectral coverage of the XMM-Newton observatory into the ultraviolet and optical range. It provides imaging and time-resolved data on targets simultaneously with observations in the EPIC and RGS. It also has the ability to track stars in its field of view, thus providing an improved post-facto aspect solution for the spacecraft. An overview of the XMM-OM and its operation is given, together with current information on the performance of the instrument.

We present the high resolution spectrum of the isolated neutron star RXJ0720.4-3125, obtained with the Reflection Grating Spectrometer on XMM-Newton, complemented with the broad band spectrum observed with the EPIC PN camera. The spectrum appears smooth, with no evidence for strong photospheric absorption or emission features. We briefly discuss the implications of our failure to detect structure in the spectrum.

We model the spin pulse intensity and hardness ratio profiles of the isolated neutron star RXJ0720.4--3125 using XMM-Newton data. The observed variation is approximately sinusoidal with a peak-to-peak amplitude of 15%, and the hardness ratio is softest slightly before flux maximum. By using polar cap models we are able to derive maximum polar cap sizes and acceptable viewing geometries. The inferred sizes of the caps turn out to be more compatible with a scenario in which the neutron star is heated by accretion, and place limits on the magnetic field strength. The hardness ratio modulation can then be explained in terms of energy-dependent beaming effects, and this constrains the acceptable models of the emerging radiation to cases in which softer photons are more strongly beamed than harder photons. An alternative explanation in terms of spatially variable absorption co-rotating in the magnetosphere is also discussed.

We present simultaneous optical polarimetry and X-ray data of the near synchronous polar RX J2115-5840. We model the polarisation data using the Stokes imaging technique of Potter et al. We find that the data are best modelled using a relatively high binary inclination and a small angle between the magnetic and spin axes. We find that for all spin-orbit beat phases, a significant proportion of the accretion flow is directed onto the lower hemisphere of the white dwarf, producing negative circular polarisation. Only for a small fraction of the beat cycle is a proportion of the flow directed onto the upper hemisphere. However, the accretion flow never occurs near the upper magnetic pole, whatever the orientation of the magnetic poles. This indicates the presence of a non-dipole field with the field strength at the upper pole significantly higher. We find that the brightest parts of the hard X-ray emitting region and the cyclotron region are closely coincident.

We have detected the optical counterpart of the proposed double degenerate polar RX J1914+24. The I band light curve is modulated on the 9.5 min period seen in X-rays. There is no evidence for any other periods. No significant modulation is seen in J. The infrared colours of RX J1914+24 are not consistent with a main sequence dwarf secondary star. Our ASCA spectrum of RX J1914+24 is typical of a heavily absorbed polar and our ASCA light curve also shows only the 9.5 min period. We find that the folded I band and X-ray light curves are out of phase. We attribute the I band flux to the irradiated face of the donor star. The long term X-ray light curve shows a variation in the observed flux of up to an order of magnitude. These observations strengthen the view that RX J1914+24 is indeed the first double degenerate polar to be detected. In this light, we discuss the synchronising mechanisms in such a close binary and other system parameters.

We re-examine the vertical structure of the post-shock flow in the accretion region of mCVs, and the X-ray emission as a function of height. We then predict X-ray light curves and phase-resolved spectra, taking into account the vertical structure, examine the implications and check whether the predicted heights are compatible with observation.

We have calculated the temperature and density structure of the hot postshock plasma in magnetically confined accretion flows, including the gravitational potential. This avoids the inconsistency of previous calculations which assume that the height of the shock is negligible. We assume a stratified accretion column with 1-d flow along the symmetry axis. We find that the calculations predict a lower shock temperature than previous calculations, with a flatter temperature profile with height. We have revised previous determinations of the masses of the white dwarf primary stars and find that for higher mass white dwarfs there is a general reduction in derived masses when the gravitational potential is included. This is because the spectrum from such flows is harder than that of previous prescriptions at intermediate energies.

Following the suggestion of Schwope et al that the magnetic cataclysmic variable RX J2115-5840 maybe a near-synchronous polar, we obtained optical polarimetry of this system over a 2 week period. From a power spectrum of the circular polarimetry data we determine that the spin period of the white dwarf and the binary orbital period which differ by 1.2%. RX J2115-5840 is thus the fourth near synchronous polar and has the shortest spin-orbit beat period: 6.3 days. By folding the data on spin, beat and orbital periods we find evidence that the accretion stream is directed towards opposite magnetic poles as the stream precesses around the white dwarf on the spin-orbit beat period. The phasing requires that the accretion flow must be directed onto the same magnetic field line at all spin-orbit beat phases implying that at some phases the flow must follow a path around the white dwarf before accreting. This is difficult to reconcile with simple views of how the accretion stream attaches onto the magnetic field of the white dwarf.

We apply our technique for indirect imaging of the accretion stream to UBVR eclipse profiles of HU Aqr obtained when the system was in a high state. We perform model fits to the eclipses using two geometries: one where the stream accretes onto one footpoint of the field line, and the other where the stream accretes onto both footpoints. We find that the stream brightness is not uniform, nor is it a simple function of the radial distance from the white dwarf. The model provides estimates of the radius where the stream couples to the magnetic field (1.0-1.3 x10^{10} cm) and of the mass transfer rate (8-76 x10^{16} g/s). We examine the wavelength dependence of the stream brightness, and find that in the case of the two-footpoint model, the sections of the magnetically-constrained stream immediately adjacent to the coupling radius appear to be hotter than the remainder of the stream. This may be due to shock heating of the plasma as it threads onto the magnetic field.

Schwope et al (1997) suggested that the newly discovered Polar RX J2115-5840 is a near-synchronous system. We have obtained circular polarisation observations of RX J2115-5840 which show that the spin and orbital periods differ by 1.2%. We find the first direct evidence of `pole-switching' in a near-synchronous Polar. Further our data requires that the accretion flow must be directed onto the same magnetic field line at all spin-orbit beat phases implying that at some phases the flow must follow a path around the white dwarf before accreting.

Cropper et al 1998 presented X-ray observations of RX J1914+24 and claimed that the most likely interpretation of their data was that of a double degenerate Polar. Here we show the preliminary results of optical and further X-ray observations. We identify the optical counterpart and show that no other period apart from the 569 sec period seen in X-rays is identified. Although the optical counterpart has no significant polarised flux in the I band, our ASCA spectrum (which is typical of a Polar) together with our photometric result adds weight to the conclusion that RX J1914+24 is the first double degenerate Polar to be discovered.

The post-shock flow in the accretion region of MCVs is considered to be a magnetically confined cooling flow. The X-ray spectrum from this hot plasma is emitted from a range of temperatures and densities. Extracting information from such spectra obtained with the next generation of X-ray satellites will rely on models for the structure of this flow which include appropriate treatment of the cooling mechanisms, boundary conditions and of absorptions and reflection components. This paper summarises the current situation and explores some of the difficulties arising in the exploitation of the line emission from these systems. It concludes by looking at improvements to the models now in the pipeline.

Modelling of the polarized cyclotron emission from magnetic cataclysmic variables (MCVs) has been a powerful technique for determining the structure of the accretion zones on the white dwarf. Until now, this has been achieved by constructing emission regions (for example arcs and spots) put in by hand, in order to recover the polarized emission. These models were all inferred indirectly from arguments based on polarization and X-ray light curves. Potter, Hakala & Cropper (1998) presented a technique (Stokes imaging) which objectively and analytically models the polarized emission to recover the structure of the cyclotron emission region(s) in MCVs. We demonstrate this technique with the aid of a test case, then we apply the technique to polarimetric observations of the AM Her system V347 Pav. As the system parameters of V347 Pav (for example its inclination) have not been well determined, we describe an extension to the Stokes imaging technique which also searches the system parameter space (GOSSIP).

The XMM satellite has a scientific payload of 3 EPIC and 2 RGS instruments which provide for both imaging and non-dispersive spectroscopy as well as medium-to-high resolution dispersive spectroscopy, in the soft to medium X-ray bands. The co-aligned OM adds the facility for simultaneous optical and UV studies. In this poster paper we describe the instrumental characteristics of XMM, comparing them to those of ASCA. We also forecast the expected X-ray spectra of a selection of both the Polars and Intermediate Polars by means of simulations based on ASCA model fits.

We review the main observational characteristics of AM Herculis stars (polars) at X-ray, EUV, UV, IR and optical wavelengths. Particular emphasis is given to multi-epoch, multi-wavelength observations of the eclipsing polar HU Aqr (RX J2107.9-0518). In AM Herculis stars the broad-band spectral energy distribution from X-rays to the IR is governed by only very small structures: the hot accretion regions on the footpoints af accreting field lines. The most extended structures in the binary systems on the other hand, the mass-donating secondary stars and the accretion streams, distinctly appear only as Doppler-shifted emission or absorption lines. They can best be studied by investigating selected narrow spectral features in the optical, ultraviolet or the near infrared. In this contribution both aspects will be highlighted, the structure of the accretion regions as inferred from multi-wavelength observations with low or no spectral resolution, as well as the structure of the secondary stars and the accretion streams as inferred from high-resolution spectral observations and Doppler mapping.