We report the results of XMM-Newton and BeppoSAX observations of the radio and X-ray emitting star LS I +61 303, likely associated with the gamma-ray source 2CG 135+01 and recently detected also at TeV energies. The data include a long XMM-Newton pointing carried out in January 2005, which provides the deepest look ever obtained for this object in the 0.3-12 keV range. During this observation the source flux decreased from a high level of 13E12 erg/cm2/s to 4E12 erg/cm2/s within 2-3 hours.This flux range is the same seen in shorter and less sensitive observations carried out in the past, but the new data show for the first time that transitions between the two levels can occur on short time scales. The flux decrease was accompanied by a significant softening of the spectrum, which is well described by a power law with photon index changing from 1.62+/-0.1 to 1.83+/-0.1. A correlation between hardness and intensity is also found when comparing different short observations spanning almost 10 years and covering various orbital phases.LS I +61 303 was detected in the 15-70 keV range with the PDS instrument in one of the BeppoSAX observations, providing evidence for variability also in the hard X-ray range. The X-ray spectra, discussed in the context of multiwavelength observations, place some interesting constraints on the properties and location of the high-energy emitting region.

The discovery of deep spectral features in the X-ray spectrum of 1E1207.4-5209 focussed the attention of the astronomical community on this radio-quiet NS, making it the most intensively observed INS ever. The harvest of X-ray photons, collected mainly by XMM-Newton, unveiled a number of unique characteristics, raising questions on this source very nature.

In August 2002, XMM-Newton devoted two full orbits to the observation of 1E 1207.4-5209, making this isolated neutron star the most deeply scrutinized galactic target of the mission. Thanks to the high throughput of the EPIC instrument, ~360,000 photons were collected from the source, allowing for a very sensitive study of the temporal and spectral behaviour of this object. The spectral data, both time-averaged and phase-resolved, yield one compelling interpretation of the observed features: cyclotron absorption from one fundamental (~0.7 keV) and three harmonics, at ~1.4, ~2.1 and ~2.8 keV. Possible physical consequences are discussed, also on the basis of the obvious phase variations of the features' shapes and depths. We also present deep VLT optical data which we have used to search for a counterpart, with negative results down to ~27.

The discovery of deep spectral features in the X-ray spectrum of 1E1207.4-5209 has pushed this Isolated Neutron Star (INS) out of the chorus line, since no other INS has shown significant features in its X-ray continuum. On August 2002, XMM-Newton devoted a two-orbit TOO observation to this target with the aim to better understand the nature of such spectral features, using much improved statistics. Indeed, the 260 ksec observation yielded 360,000 photons from 1E1207.4-5209, allowing for a very sensitive study of the temporal and spectral behaviour of this object.