We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the non-ideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern, the distortion of the electron diffusion size and shape, and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread-functions and the sub-linear variance in flats. The method is implemented in the Photon Simulator (PhoSim) and the code is publically available.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

Given the current limited knowledge of meteor plasma micro-physics and its interaction with the surrounding atmosphere and ionosphere, meteors are a highly interesting observational target for high-resolution wide-field astronomical surveys. Such surveys are capable of resolving the physical size of meteor plasma heads, but they produce large volumes of images that need to be automatically inspected for possible existence of long linear features produced by meteors. Here we show how big aperture sky survey telescopes detect meteors as defocused tracks with a central brightness depression. We derive an analytic expression for a defocused point source meteor track and use it to calculate brightness profiles of meteors modeled as uniform brightness disks. We apply our modeling to meteor images as seen by the SDSS and LSST telescopes. The expression is validated by Monte Carlo ray-tracing simulations of photons traveling through the atmosphere and the LSST telescope optics. We show that estimates of the meteor distance and size can be extracted from the measured FWHM and the strength of the central dip in the observed brightness profile. However, this extraction becomes difficult when the defocused meteor track is distorted by the atmospheric seeing or contaminated by a long lasting glowing meteor trail. The FWHM of satellites tracks is distinctly narrower than meteor values, which enables removal of a possible confusion between satellites and meteors.

We employ electrostatic conversion drift calculations to match CCD pixel signal covariances observed in flat field exposures acquired using candidate sensor devices for the LSST Camera. We thus constrain pixel geometry distortions present at the end of integration, based on signal images recorded. We use available data from several operational voltage parameter settings to validate our understanding. Our primary goal is to optimize flux point-spread function (FPSF) estimation quantitatively, and thereby minimize sensor-induced errors which may limit performance in precision astronomy applications. We consider alternative compensation scenarios that will take maximum advantage of our understanding of this underlying mechanism in data processing pipelines currently under development. To quantitatively capture the pixel response in high-contrast/high dynamic range operational extrema, we propose herein some straightforward laboratory tests that involve altering the time order of source illumination on sensors, within individual test exposures. Hence the word {\it hysteretic} in the title of this paper.

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

We build on previous efforts to model CCD sensors, during illumination and collection of conversions. We use a finite summation of simple, electrostatic field models. The upgraded functionality of our framework provides specific predictions for perturbations in pixel boundary enclosures (e.g., at the backside window) and the bookkeeping capability to stack those perturbations so that they may be utilized as Greens functions -- portable calculation results that may be generically applied to a range of precision astronomy related problems that naturally including astrometric, photometric and shape transfer issues. We approach the topic of using ancillary pixel data, derived from the Greens function and the registered image, to analyze sky data and constrain object parameters of astronomical targets.

Near-future astronomical survey experiments, such as LSST, possess system requirements of unprecedented fidelity that span photometry, astrometry and shape transfer. Some of these requirements flow directly to the array of science imaging sensors at the focal plane. Availability of high quality characterization data acquired in the course of our sensor development program has given us an opportunity to develop and test a framework for simulation and modeling that is based on a limited set of physical and geometric effects. In this paper we describe those models, provide quantitative comparisons between data and modeled response, and extrapolate the response model to predict imaging array response to astronomical exposure. The emergent picture departs from the notion of a fixed, rectilinear grid that maps photo-conversions to the potential well of the channel. In place of that, we have a situation where structures from device fabrication, local silicon bulk resistivity variations and photo-converted carrier patterns still accumulating at the channel, together influence and distort positions within the photosensitive volume that map to pixel boundaries. Strategies for efficient extraction of modeling parameters from routinely acquired characterization data are described. Methods for high fidelity illumination/image distribution parameter retrieval, in the presence of such distortions, are also discussed.

We describe the drift field in thick depleted silicon sensors as a superposition of a one-dimensional backdrop field and various three-dimensional perturbative contributions that are physically motivated. We compute trajectories for the conversions along the field lines toward the channel and into volumes where conversions are confined by the perturbative fields. We validate this approach by comparing predictions against measured response distributions seen in five types of fixed pattern distortion features. We derive a quantitative connection between "tree ring" flat field distortions to astrometric and shape transfer errors with connections to measurable wavelength dependence - as ancillary pixel data that may be used in pipeline analysis for catalog population. Such corrections may be tested on DECam data, where correlations between tree ring flat field distortions and astrometric errors - together with their band dependence - are already under study. Dynamic effects, including the brighter-fatter phenomenon for point sources and the flux dependence of flat field fixed pattern features are approached using perturbations similar in form to those giving rise to the fixed pattern features. These in turn provide drift coefficient predictions that can be validated in a straightforward manner. Once the three parameters of the model are constrained using available data, the model is readily used to provide predictions for arbitrary photo-distributions with internally consistent wavelength dependence provided for free.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

Understanding the origins and distribution of matter in the Universe is one of the most important quests in physics and astronomy. Themes range from astro-particle physics to chemical evolution in the Galaxy to cosmic nucleosynthesis and chemistry in an anticipation of a full account of matter in the Universe. Studies of chemical evolution in the early Universe will answer questions about when and where the majority of metals were formed, how they spread and why they appar today as they are. The evolution of matter in our Universe cannot be characterized as a simple path of development. In fact the state of matter today tells us that mass and matter is under constant reformation through on-going star formation, nucleosynthesis and mass loss on stellar and galactic scales. X-ray absorption studies have evolved in recent years into powerful means to probe the various phases of interstellar and intergalactic media. Future observatories such as IXO and Gen-X will provide vast new opportunities to study structure and distribution of matter with high resolution X-ray spectra. Specifically the capabilities of the soft energy gratings with a resolution of R=3000 onboard IXO will provide ground breaking determinations of element abundance, ionization structure, and dispersion velocities of the interstellar and intergalactic media of our Galaxy and the Local Group

We report on X-ray observations of the supernova remnant 0509-67.5 in the Large Magellanic Cloud with XMM-Newton X-ray observatory. We use the imaging spectroscopy (EPIC) and Reflective Grating Spectrometer (RGS) data to investigate properties of the remnant and its environment. The X-ray spectra were analyzed with SPEX software package. In addition to this we performed a numerical hydrodynamic simulation of the remnant. The EPIC data show prominent Fe K line emission, but the deduced overall amount of iron in the shocked ejecta is low. The data also show that the remnant has an asymmetric ejecta structure: the bright southwest region of the remnant shows an overabundance of metals. The analysis of the RGS spectrum shows that the remnant has a high lines velocity broadening of ~5000 km/s. We found a hydrodynamical model for the remnant with basic hydrodynamical and spectral parameters similar to the observed ones. The data analysis show that the reverse shock just recently reached iron layers of the ejecta. The brightness enhancement in the southwest region could be a sign of an asymmetric explosion or it could be the result of a density enhancement of the interstellar medium. We constructed numerical models which are in good agreement with the observations, with circumstellar density of 3e-25 g/cm^3, age of ~400 years, velocities of ~5000 km/s and an electron to ion temperature ratio of 0.01.

We present 3D gasdynamic simulations of the Herbig Haro object HH 555. HH 555 is a bipolar jet emerging from the tip of an elephant trunk entering the Pelican Nebula from the adjacent molecular cloud. Both beams of HH 555 are curved away from the center of the H II region. This indicates that they are being deflected by a side-wind probably coming from a star located inside the nebula or by the expansion of the nebula itself. HH 555 is most likely an irradiated jet emerging from a highly embedded protostar, which has not yet been detected. In our simulations we vary the incident photon flux, which in one of our models is equal to the flux coming from a star 1 pc away emitting 5x10^48 ionizing (i. e., with energies above the H Lyman limit) photons per second. An external, plane-parallel flow (a ``side-wind'') is coming from the same direction as the photoionizing flux. We have made four simulations, decreasing the photon flux by a factor of 10 in each simulation. We discuss the properties of the flow and we compute Halpha emission maps (integrated along lines of sight). We show that the level of the incident photon flux has an important influence on the shape and visibility of the jet. If the flux is very high, it causes a strong evaporation of the neutral clump, producing a photoevaporated wind traveling in the direction opposite to the incident flow. The interaction of the two flows creates a double shock ``working surface'' around the clump protecting it and the jet from the external flow. The jet only starts to curve when it penetrates through the working surface.

In a series of papers, Nicastro et al. have reported the detection of z>0 OVII absorption features in the spectrum of Mrk421 obtained with the Chandra Low Energy Transmission Grating Spectrometer (LETGS). We evaluate this result in the context of a high quality spectrum of the same source obtained with the Reflection Grating Spectrometer (RGS) on XMM-Newton. The data comprise over 955ks of usable exposure time and more than 26000 counts per 50 milliAngstrom at 21.6 Angstroms. We concentrate on the spectrally clean region (21.3 < lambda < 22.5 Angstroms) where sharp features due to the astrophysically abundant OVII may reveal an intervening, warm--hot intergalactic medium (WHIM). We do not confirm detection of any of the intervening systems claimed to date. Rather, we detect only three unsurprising, astrophysically expected features down to the Log(N_i)~14.6 (3 sigma) sensitivity level. Each of the two purported WHIM features is rejected with a statistical confidence that exceeds that reported for its initial detection. While we can not rule out the existence of fainter, WHIM related features in these spectra, we suggest that previous discovery claims were premature. A more recent paper by Williams et al. claims to have demonstrated that the RGS data we analyze here do not have the resolution or statistical quality required to confirm or deny the LETGS detections. We show that our careful analysis resolves the issues encountered by Williams et al. and recovers the full resolution and statistical quality of the RGS data. We highlight the differences between our analysis and those published by Williams et al. as this may explain our disparate conclusions.

Recently the first detections of highly ionised gas associated with two Warm-Hot Intergalactic Medium (WHIM) filaments have been reported. The evidence is based on X-ray absorption lines due to O VII and other ions observed by Chandra towards the bright blazar Mrk 421. We investigate the robustness of this detection by a re-analysis of the original Chandra LETGS spectra, the analysis of a large set of XMM-Newton RGS spectra of Mrk 421, and additional Chandra observations. We address the reliability of individual spectral features belonging to the absorption components, and assess the significance of the detection of these components. We also use Monte Carlo simulations of spectra. We confirm the apparent strength of several features in the Chandra spectra, but demonstrate that they are statistically not significant. This decreased significance is due to the number of redshift trials that are made and that are not taken into account in the original discovery paper. Therefore these features must be attributed to statistical fluctuations. This is confirmed by the RGS spectra, which have a higher signal to noise ratio than the Chandra spectra, but do not show features at the same wavelengths. Finally, we show that the possible association with a Ly-alpha absorption system also lacks sufficient statistical evidence. We conclude that there is insufficient observational proof for the existence of the two proposed WHIM filaments towards Mrk 421, the brightest X-ray blazar on the sky. Therefore, the highly ionised component of the WHIM still remains to be discovered.

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0<z < 1.5. Such an X-ray survey will discover all collapsed structures with mass above 3.5e14 solar masss at redshifts z<2 (i.e. the full range where such objects are expected) in the high Galactic latitude sky. Above this mass threshold the tight correlations between X-ray observables and mass allow direct interpretation of the data. We describe the constraints on Dark Energy that can be inferred from such a survey, using powerful self-calibration techniques to relate the X-ray observables (luminosity and temperature) to the underlying mass.

We present a non-LTE model atmosphere analysis of Chandra HRC-S/LETG and XMM-Newton RGS spectroscopy of the prototypical supersoft source CAL 83 in the Large Magellanic Cloud. Taken with a 16-month interval, the Chandra and XMM-Newton spectra are very similar. They reveal a very rich absorption line spectrum from the hot white dwarf photosphere, but no spectral signatures of a wind. We also report a third X-ray off-state during a later Chandra observation, demonstrating the recurrent nature of CAL 83. Moreover, we found evidence of short-timescale variability in the soft X-ray spectrum. We completed the analysis of the LETG and RGS spectra of CAL 83 with new NLTE line-blanketed model atmospheres that explicitly include 74 ions of the 11 most abundant species. We successfully matched the Chandra and XMM-Newton spectra assuming a model composition with LMC metallicity. We derived the basic stellar parameters of the hot white dwarf, but the current state of atomic data in the soft X-ray domain precludes a detailed chemical analysis. We have obtained the first direct spectroscopic evidence that the white dwarf is massive (Mwd > 1 Msun). The short timescale of the X-ray off-states is consistent with a high white dwarf mass. Our analysis thus provides direct support for supersoft sources as likely progenitors of SN Ia.

We present a preliminary analysis of XMM-Newton observations of the oxygen-rich supernova remnant G292.0+1.8 (MSH 11-54). Although the spatial extent of the remnant is 8 arcmin the bright central bar is narrow (1'-2') resulting in RGS spectra of a high spectral quality. This allows us to spectroscopically identify a cool, Te = 0.3 keV, and underionized component, resolve details of the Fe-L complex, and resolve the forbidden and resonant lines of the O VII triplet. We are also able to constrain the kinematics of the remnant using Ne IX as observed in the second order spectrum, and O VIII in the first order spectrum. We do not find evidence for O VII line shifts or Doppler broadening (sigma_v < 731 km/s), but line broadening of the Ne X Ly-alpha line seems to be present, corresponding to sigma_v ~ 1500 km/s.

A long 280 ks observation of the Seyfert 1 galaxy NGC 3783 with XMM-Newton is reported. We focus on the oxygen line complex between 17 and 24 A as measured with the RGS spectrometers. Accurate absorption column densities and emission line fluxes are obtained. We explore several options for the geometry and physical form of the emitting and absorbing gas. The lack of change in ionization in the absorber despite an increase in continuum flux during the observation restricts the high-ionization (O-K) and the low-ionization (Fe-M) gas to distances of at least 0.5 pc and 2.8 pc, respectively, away from the central source. Given the P-Cygni type profiles in the resonance spectral lines and the similar velocity widths, column densities, and ionization structure inferred separately from the emission and absorption lines, it is tempting to relate the X-ray narrow-line emitting plasma with the X-ray absorbing gas. Under this assumption, the scenario of dense clumped clouds can be ruled out. Conversely, extended ionization cones (r > 10 pc) are consistent with the observation independent of this assumption. These findings are in stark contrast with the picture of numerous clumpy (n_e > 10^9 cm^-3) clouds drawn recently from UV spectra, but it is consistent with the extended X-ray emission cones observed directly in Seyfert 2 galaxies.

High resolution X-ray spectra of the Reflection Grating Spectrometer (RGS) on board the XMM satellite are used to resolve the oxygen K absorption edge. By combining spectra of low and high extinction sources, the observed absorption edge can be split in the true interstellar (ISM) extinction and the instrumental absorption. The detailed ISM edge structure closely follows the edge structure of neutral oxygen as derived by theoretical R-matrix calculations. However, the position of the theoretical edge requires a wavelength shift. In addition the detailed instrumental RGS absorption edge structure is presented. All results are verified by comparing to a subset of Chandra LETG-HRC observations.

We report on the observation of OVII Doppler line broadening in a compact knot at the edge of SN 1006 detected with the Reflective Grating Spectrometer on board XMM-Newton. The observed line width of sigma = 3.4+/-0.5 eV at a line energy of 574 eV indicates an oxygen temperature of kT = 528+/-150 keV. Combined with the observed electron temperature of 1.5 keV the observed broadening is direct evidence for temperature non-equilibration in high Mach number shocks, and slow subsequent equilibration. The OVII line emission allows an accurate determination of the ionization state of the plasma, which is characterized by a relatively high forbidden line contribution, indicating log(n_e t) ~ 9.2.

Recent observations with the dispersive X-ray spectrometers aboard Chandra and Newton Observatory have begun to probe the properties of the X-ray intergalactic medium (IGM) at small redshifts. Using large quantities (~950 ksec) of spectroscopic data acquired using the RGS aboard Newton Observatory, we investigated the intervening material toward three low redshift, high Galactic latitude AGNs with nominally featureless spectra: Mrk421, PKS2155-304 and 3C273. Each spectrum provides clear evidence for what appears to be a local (z~0), highly ionized absorbing medium betrayed by the OVII 1s-2p resonance transition feature seen at 21.6A (N[OVII] ~ 1E16 cm-2). Measurements are also made for the Lyman alpha transition of the adjacent ionization state, (OVIII; 18.97A), which potentially constrains the absorber's temperature. Finally, in a collisional equilibrium approximation, upper limits to diffuse emission intensities place upper limits on the electron density (ne < 2E-4 cm-3), lower limits on the scale length of the absorber (L > 140 kpc) and lower limits on its mass (M > 5E10 M[sun]). Limits on the absorber's scale length and its velocity distribution lead us to identify it with the Local Group.

The X-ray emission from the young supernova remnant (SNR) N103B is measured and analysed using the high-resolution cameras and spectrometers on board XMMM and CHANDRA. The spectrum from the entire remnant is reproduced very well with three plasma components of kT=0.55, 0.65, and 3.5 keV, corresponding roughly to line emission by the O-K, Fe-L, and Fe-K species, respectively. Narrow band images reveal different morphologies for each component. The kT=0.65keV component, which dominates the emission measure (4.5E65m^{-3}), is in ionisation equilibrium. This provides a lower limit of 3000 yrs to the age of the remnant, which is considerably older than the previously assumed age of the remnant (1500 yrs). Based on the measured energy of the Fe-K feature at 6.5 keV, the hot (3.5 keV) component is found to be recently shocked (~200 yrs) and still ionising. The high elemental abundances of O and Ne and the low abundance of Fe could imply that SNR N103B originated from a type II supernova (SN) rather than a type Ia SN as previously thought.

We present an analysis of the X-ray emission of the supernova remnant MSH14-63, which was partially covered by four observations with XMM-Newton. The detection of Fe K emission at 6.4 keV, and the lack of spatial correlation between hard X-ray and radio emission is evidence against a dominant X-ray synchrotron component. We argue that the hard X-ray continuum is best explained by non-thermal bremsstrahlung from a supra-thermal tail to an otherwise cool electron gas. The existence of low electron temperatures, required to explain the absence of line emission, is supported by low temperatures found in other parts of the remnant, which are as low as 0.2 keV in some regions.

First results from high-resolution coronal spectroscopy of flares with the Reflection Grating Spectrometers on board the \textit{XMM-Newton} satellite are reviewed. Rotational modulation in the X-ray light curve of HR 1099 is discussed. Results from time-dependent spectroscopy of flares in the active stars HR 1099, AB Dor, YY Gem are also presented. Variations in the shape of the emission measure distributions, in the abundances and in the average density of the cool plasma are discussed.

First results from high-resolution coronal spectroscopy with the {\it XMM-Newton} Reflection Grating Spectrometers (RGS) are reviewed. Five stellar systems (HR 1099, Capella, Procyon, YY Gem, AB Dor) have been observed. The emphasis of the present paper is on overall coronal structure. Elemental abundances in {\it active stars} are found to be `anomalous' in the sense that they tend to increase with increasing First Ionization Potential (FIP - i.e., signifying an inverse FIP effect). Coronal densities are measured at levels of a few times $10^{10}$ cm$^{-3}$ for cooler plasma, although there are indications for very high densities in the hotter plasma components.

The observation of the supernova remnant N132D by the scientific instruments on board the XMM-Newton satellite is presented. The X-rays from N132D are dispersed into a detailed line-rich spectrum using the Reflection Grating Spectrometers. Spectral lines of C, N, O, Ne, Mg, Si, S, and Fe are identified. Images of the remnant, in narrow wavelength bands, produced by the European Photon Imaging Cameras reveal a complex spatial structure of the ionic distribution. While K-shell Fe emission seems to originate near the centre, all of the other ions are observed along the shell. A high O VII / O VIII emission ratio is detected on the northeastern edge of the remnant. This can be a sign of hot ionising conditions, or it can reflect relatively cool gas. Spectral fitting of the CCD spectrum suggests high temperatures in this region, but a detailed analysis of the atomic processes involved in producing the O VII spectral lines leads to the conclusion that the intensities of these lines alone cannot provide a conclusive distinction between the two scenarios.

We present the first high resolution photospheric X-ray spectrum of a Supersoft X-ray Source, the famous CAL~83 in the Large Magellanic Cloud. The spectrum was obtained with the Reflection Grating Spectrometer on XMM-Newton during the Calibration/Performance Verification phase of the observatory. The spectrum covers the range 20-40 A at an approximately constant resolution of 0.05 A, and shows very significant, intricate detail, that is very sensitive to the physical properties of the object. We present the results of an initial investigation of the spectrum, from which we draw the conclusion that the spectral structure is probably dominated by numerous absorption features due to transitions in the L-shells of the mid-$Z$ elements and the M-shell of Fe, in addition to a few strong K-shell features due to CNO.

The RS CVn binary system HR 1099 was extensively observed by the XMM-Newton observatory in February 2000 as its first-light target. A total of 570 ks of exposure time was accumulated with the Reflection Grating Spectrometers (RGS). The integrated X-ray spectrum between 5-35 Angstrom is of unprecedented quality and shows numerous features attributed to transitions of the elements C, N, O, Ne, Mg, Si, S, Fe, and Ni. We perform an in-depth study of the elemental composition of the average corona of this system, and find that the elemental abundances strongly depend on the first ionisation potential (FIP) of the elements. But different from the solar coronal case, we find an inverse FIP effect, i.e., the abundances (relative to oxygen) increase with increasing FIP. Possible scenarios, e.g., selective enrichment due to Ne-rich flare-like events, are discussed.