We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

We present the spectral analysis of a 200~ks observation of the broad-line radio galaxy 3C~120 performed with the high energy transmission grating (HETG) spectrometer on board the \emph{Chandra} X-ray Observatory. We find (i) a neutral absorption component intrinsic to the source with column density of $\text{log}N_H = 20.67\pm0.05$~cm$^{-2}$, (ii) no evidence for a warm absorber with an upper limit on the column density of just $\text{log}N_H < 19.7$~cm$^{-2}$ assuming the typical ionization parameter log$\xi$$\simeq$2.5~erg~s$^{-1}$~cm, the warm absorber may instead be replaced by (iii) a hot emitting gas with temperature $kT \simeq 0.7$~keV observed as soft X-ray emission from ionized Fe L-shell lines which may originate from a kpc scale shocked bubble inflated by the AGN wind or jet with a shock velocity of about 1,000~km~s$^{-1}$ determined by the emission line width, (iv) a neutral Fe K$\alpha$ line and accompanying emission lines indicative of a Compton-thick cold reflector with low reflection fraction $R\simeq0.2$, suggesting a large opening angle of the torus, (v) a highly ionized Fe~XXV emission feature indicative of photoionized gas with ionization parameter log$\xi$$=$$3.75^{+0.27}_{-0.38}$~erg~s$^{-1}$~cm and a column density of $\text{log}N_H > 22$~cm$^{-2}$ localized within $\sim$2~pc from the X-ray source, and (vi) possible signatures for a highly ionized disk wind. Together with previous evidence for intense molecular line emission, these results indicate that 3C~120 is likely a late state merger undergoing strong AGN feedback.

The project Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS (MUCHFUSS) aims at finding hot subdwarf stars with massive compact companions like massive white dwarfs (M>1.0 M$_\odot$), neutron stars, or stellar-mass black holes. We present orbital and atmospheric parameters and put constraints on the nature of the companions of 12 close hot subdwarf B star (sdB) binaries found in the course of the MUCHFUSS project. The systems show periods between 0.14 and 7.4 days. Three systems most likely have white dwarf companions. SDSS J083006.17+475150.3 is likely to be a rare example of a low-mass helium-core white dwarf. SDSS J095101.28+034757.0 shows an excess in the infrared that probably originates from a third companion in a wide orbit. SDSS J113241.58-063652.8 is the first helium deficient sdO star with a confirmed close companion. This study brings to 142 the number of sdB binaries with orbital periods of less than 30 days and with measured mass functions. We present an analysis of the minimum companion mass distribution and show that it is bimodal. One peak around 0.1 M$_\odot$ corresponds to the low-mass main sequence and substellar companions. The other peak around 0.4 M$_\odot$ corresponds to the white dwarf companions. The derived masses for the white dwarf companions are significantly lower than the average mass for single carbon-oxygen white dwarfs. In a T$_{\rm eff}$-log(g) diagram of sdB+dM companions, we find signs that the sdB components are more massive than the rest of the sample. The full sample was compared to the known population of extremely low-mass white dwarf binaries as well as short-period white dwarfs with main sequence companions. Both samples show a significantly different companion mass distribution. We calculate merger timescales and timescales when the companion will fill its Roche Lobe and the system evolves into a cataclysmic variable.

We present solutions for the velocity field and mass-loss rates for 2D axisymmetric outflows, as well as for the case of mass accretion through the use of the Lambert W-function. For the case of a rotating radiation-driven wind the velocity field is obtained analytically using a parameterised description of the line acceleration that only depends on radius r at any given latitude $\theta$. The line acceleration g(r) is obtained from Monte-Carlo multi-line radiative transfer calculations. The critical/sonic point of our equation of motion varies with latitude $\theta$. Furthermore, an approximate analytical solution for the supersonic flow of a rotating wind is derived, which is found to closely resemble the exact solution. For the simultaneous solution of the mass-loss rate and velocity field, we use the iterative method of our 1D method extended to the non-spherical 2D case. We apply the new theoretical expressions with our iterative method to the stellar wind from a differentially rotating 40 $M_{sun}$ O5-V main sequence star as well as to a 60 $M_{sun}$ O-giant star, and we compare our results to previous studies that are extensions of the Castor et al. (1975, ApJ, 195, 157) CAK formalism. Next, we account for the effects of oblateness and gravity darkening. Our numerical results predict an equatorial decrease of the mass-loss rate, which would imply that (surface-averaged) total mass-loss rates are lower than for the spherical 1D case, in contradiction to the Maeder & Meynet (2000, A&A, 361, 159) formalism that is oftentimes employed in stellar evolution calculations for rotating massive stars. To clarify the situation in nature we discuss observational tests to constrain the shapes of large-scale 2D stellar winds.

Observations show that the time of onset of dust formation in classical novae depends strongly on their speed class, with dust typically taking longer to form in slower novae. Using empirical relationships between speed class, luminosity and ejection velocity, it can be shown that dust formation timescale is expected to be essentially independent of speed class. However, following a nova outburst the spectrum of the central hot source evolves, with an increasing proportion of the radiation being emitted short-ward of the Lyman limit. The rate at which the spectrum evolves also depends on the speed class. We have therefore refined the simple model by assuming photons at energies higher than the Lyman limit are absorbed by neutral hydrogen gas internal to the dust formation sites, therefore preventing these photons reaching the nucleation sites. With this refinement the dust formation timescale is theoretically dependent on speed class and the results of our theoretical modification agree well with the observational data. We consider two types of carbon-based dust, graphite and amorphous carbon, with both types producing similar relationships. Our results can be used to predict when dust will form in a nova of a given speed class and hence when observations should optimally be taken to detect the onset of dust formation.

We study the X-ray properties of LINER sources with definite detection of a broad H(alpha) emission line in their optical spectra, LINER 1s from Ho et al. sample. These objects preferentially harbor a low luminosity active nucleus at the center and show small or no intrinsic absorption (<10^(22) cm^(-1)). We analyzed all available X-ray archived XMM-Newton and Chandra observations of 13 LINER 1s satisfying the above criterion in a systematic homogeneous way. We looked for any correlations between the X-ray properties and the intrinsic parameters of our sample of LINER 1s. An absorbed power-law gave a good fit to the spectra of 9 out of the 13 sources. A thermal component and an absorbed power-law were required in the remaining 4 sources. We found a photon index between 1.3\pm0.2 for the hardest source and 2.4^(+0.2)_(-0.3) for the softest one with a mean value of 1.9\pm0.2 and a dispersion sigma=0.3. The thermal component had a mean temperature kT~0.6 keV. Significant short (hours to days) time-scale variability is not common in the present sample and was observed in only 2 sources (NGC 3226 and NGC 4278). Three other sources indicate a possible variability with a low K-S test probability (2%-4%) that the nuclear emission originates from a constant source. Significant variability on months to years time-scales is detected in 7 out of the 9 sources observed more than once. No significant Fe K(alpha) emission line at 6.4 keV was detected and upper limits were derived for the 4 sources with a high enough signal to noise ratio around 6 keV. Finally, we established, for the first time for a sample of LINER 1s, that the photon index is significantly anticorrelated to L_(2-10 keV)/L_(Edd). This anticorrelation is similar to the one seen in XRBs in their low/hard state where a radiatively inefficient accretion flow is thought to be responsible for the X-ray emitted energy.

We present observations of the Rossiter-McLaughlin effect for the transiting exoplanet systems WASP-1, WASP-24, WASP-38 and HAT-P-8, and deduce the orientations of the planetary orbits with respect to the host stars' rotation axes. The planets WASP-24b, WASP-38b and HAT-P-8b appear to move in prograde orbits and be well aligned, having sky-projected spin orbit angles consistent with zero: {\lambda} = -4.7 \pm 4.0{\deg}, {\lambda} = 15 + 33{\deg}/-43{\deg} and {\lambda} = -9.7 +9.0{\deg}/-7.7{\deg}, respectively. The host stars have Teff < 6250 K and conform with the trend of cooler stars having low obliquities. WASP-38b is a massive planet on a moderately long period, eccentric orbit so may be expected to have a misaligned orbit given the high obliquities measured in similar systems. However, we find no evidence for a large spin-orbit angle. By contrast, WASP-1b joins the growing number of misaligned systems and has an almost polar orbit, {\lambda} = -79 +4.5{\deg}/-4.3{\deg}. It is neither very massive, eccentric nor orbiting a hot host star, and therefore does not share the properties of many other misaligned systems.

We present light curves of three classical novae (KT Eridani, V598 Puppis, V1280 Scorpii) and one recurrent nova (RS Ophiuchi) derived from data obtained by the Solar Mass Ejection Imager (SMEI) on board the Coriolis satellite. SMEI provides near complete sky-map coverage with precision visible-light photometry at 102-minute cadence. The light curves derived from these sky maps offer unprecedented temporal resolution around, and especially before, maximum light, a phase of the nova eruption normally not covered by ground-based observations. They allow us to explore fundamental parameters of individual objects including the epoch of the initial explosion, the reality and duration of any pre-maximum halt (found in all three fast novae in our sample), the presence of secondary maxima, speed of decline of the initial light curve, plus precise timing of the onset of dust formation (in V1280 Sco) leading to estimation of the bolometric luminosity, white dwarf mass and object distance. For KT Eri, Liverpool Telescope SkyCamT data confirm important features of the SMEI light curve and overall our results add weight to the proposed similarities of this object to recurrent rather than to classical novae. In RS Oph, comparison with hard X-ray data from the 2006 outburst implies that the onset of the outburst coincides with extensive high velocity mass-loss. It is also noted that two of the four novae we have detected (V598 Pup and KT Eri) were only discovered by ground-based observers weeks or months after maximum light, yet these novae reached peak magnitudes of 3.46 and 5.42 respectively. This emphasizes the fact that many bright novae per year are still overlooked, particularly those of the very fast speed class. Coupled with its ability to observe novae in detail even when relatively close to the Sun in the sky, we estimate that as many as 5 novae per year may be detectable by SMEI.

We present VLT eclipse photometry for the giant planet CoRoT-1b. We observed a transit in the R-band filter and an occultation in a narrow filter centered on 2.09 microns. Our analysis of this new photometry and published radial velocities, in combination with stellar-evolutionary modeling, leads to a planetary mass and radius of 1.07 (+0.13,-0.18) M_Jup and 1.45 (+0.07,-0.13) R_Jup, confirming the very low density previously deduced from CoRoT photometry. The large occultation depth that we measure at 2.09 microns (0.278 (+0.043,-0.066) %) is consistent with thermal emission and is better reproduced by an atmospheric model with no redistribution of the absorbed stellar flux to the night side of the planet.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

Our \textit{Swift} observations of RS Oph form an unprecedented X-ray dataset to undertake investigations of both the central source and the interaction of the outburst ejecta with the circumstellar environment. Over the first month, the XRT data are dominated by emission from rapidly evolving shocks. We discuss the differences in derived parameters from those found for \textit{RXTE} at early times and the evolution of the X-ray emission to much later times. It is apparent that at late times several emission components are present. We find no strong evidence of the proposed shock break-out in our data.

We present new Spitzer Space Telescope observations of stars in the young ~5 Myr gamma Velorum stellar cluster. Combining optical and 2MASS photometry, we have selected 579 stars as candidate members of the cluster. With the addition of the Spitzer mid-infrared data, we have identified 5 debris disks around A-type stars, and 5-6 debris disks around solar-type stars, indicating that the strong radiation field in the cluster does not completely suppress the production of planetesimals in the disks of cluster members. However, we find some evidence that the frequency of circumstellar primordial disks is lower, and the IR flux excesses are smaller than for disks around stellar populations with similar ages. This could be evidence for a relatively fast dissipation of circumstellar dust by the strong radiation field from the highest mass star(s) in the cluster. Another possibility is that gamma Velorum stellar cluster is slightly older than reported ages and the the low frequency of primordial disks reflects the fast disk dissipation observed at ~5 Myr.

We present a Spitzer Space Telescope imaging survey of the most massive Galactic globular cluster, omega Centauri, and investigate stellar mass loss at low metallicity and the intracluster medium (ICM). The survey covers approximately 3.2x the cluster half-mass radius at 3.6, 4.5, 5.8, 8, and 24 microns, resulting in a catalog of over 40,000 point-sources in the cluster. Approximately 140 cluster members ranging 1.5 dex in metallicity show a red excess at 24 microns, indicative of circumstellar dust. If all of the dusty sources are experiencing mass loss, the cumulative rate of loss is estimated at 2.9 - 4.2 x 10^(-7) solar masses per year, 63% -- 66% of which is supplied by three asymptotic giant branch stars at the tip of the Red Giant Branch (RGB). There is little evidence for strong mass loss lower on the RGB. If this material had remained in the cluster center, its dust component (> 1 x 10^(-4) solar masses) would be detectable in our 24 and 70 micron images. While no dust cloud located at the center of omega Cen is apparent, we do see four regions of very faint, diffuse emission beyond two half-mass radii at 24 microns. It is unclear whether these dust clouds are foreground emission or are associated with omega Cen. In the latter case, these clouds may be the ICM in the process of escaping from the cluster.

In ESO period 65 (April-September 2000) the large programme 165.N-0276, led by Roger Cayrel, began making use of UVES at the Kueyen VLT telescope. Known within the Team and outside as ``First Stars'', it was aimed at obtaining high resolution, high signal-to-noise ratio spectra in the range 320 nm -- 1000 nm for a large sample of extremely metal-poor (EMP) stars identified from the HK objective prism survey.In this contribution we highlight the main results of the large programme.

The evolution of the first generations of stars at zero or extremly low metallicity, and especially some crucial properties like the primary N14 production, is charactarized by convective-reactive mixing events that are mostly absent from similar evolution phases at solar-like metallicity. These episodes occur when unprocessed H-rich material is mixed accross a convective boundary into C12 rich He-burning material, as for example in He-shell flashes of extremely-low metallicity AGB stars. In this paper we describe the astrophysical context of such convective-reactive events, including the difficulty of current one-dimensional stellar evolution models to correctly simulate these evolutionary phases. We then describe the requirements and current state of modeling convective-reactive processes in the first stars environment. We demonstrate some of the new concepts that we are applying to this problem, i.e. the highly accurate PPB advection scheme in the framework of PPM hydrodynamic simulations of mixing accross a very stiff convective boundary. We show initial results of such simulations that address the first non-reactive step of this problem, which is the entrainment of H at the top boundary of the He-shell flash convection zone.

We have used the Advanced Camera for Surveys on board the Hubble Space Telescope to image the core of the globular cluster M15 in the far-ultraviolet (FUV) waveband. Based on these observations, we identify the FUV counterpart of the recently discovered low-mass X-ray binary M15 X-2. Our time-resolved FUV photometry shows a modulation with 0.062+/-0.004 mag semi-amplitude and we clearly detect a period of 22.5806+/-0.0002 min. We have carried out extensive Monte Carlo simulations which show that the signal is consistent with being coherent over the entire observational time range of more than 3000 cycles. This strongly suggests that it represents the orbital period of the binary system. M15 X-2 is FUV bright (approx. 17 mag) and is characterized by an extremely blue spectral energy distribution (F_lambda ~ lambda^{-2.0}). We also find evidence for an excess of flux between 1500 and 1600 \AA and probably between 1600 and 2000 \AA, which might be due to CIV 1550 and HeII 1640 emission lines. We also show that M15 X-2's X-ray luminosity can be powered by accretion at the rate expected for gravitational-wave-driven mass transfer at this binary period. The observed FUV emission appears to be dominated by an irradiated accretion disk around the neutron star primary, and the variability can be explained by irradiation of the low-mass white dwarf donor if the inclination of the system is approx. 34 degree. We conclude that all observational characteristics of M15 X-2 are consistent with it being an ultracompact X-ray binary, only the third confirmed such object in a globular cluster.

Understanding post-common-envelope evolution is important in the studies of close-binary systems. The majority of all interacting binaries with a compact object in their system (e.g. cataclysmic variables, X-ray binaries) are thought to have passed through a common-envelope (CE) phase. Current models of interacting binaries all assume, as a starting point, that there is no significant modification of the secondary star compared with a normal star of similar spectral type. The extent to which the CE significantly alters the composition of the secondary star has yet to be determined. We are studying the M-type secondary in the pre-CV system EG UMa, in order to determine its atmospheric parameters by comparison with synthetic spectra generated using PHOENIX. Absorption lines due to several elements have been used in investigating the stellar parameters of effective temperature, gravity and over elemental abundance. In addition, we are searching for anomalies due to s-process elements (such as Ba, Sr, Rb, Y), which may have been deposited in the atmosphere during the CE stage, and prove that CE evolution has occurred. EG UMa displays strong YO absorption bands, which are normally associated with giants and S stars. These bands may have been formed as a result of an Y-abundance enhancement introduced by the CE, which would be consistent with the observed possible Rb and Sr enrichment.

We take a critical look at the lithium depletion boundary (LDB) technique that has recently been used to derive the ages of open clusters. We identify the sources of experimental and systematic error and show that the probable errors are larger by approximately a factor two than presently claimed in the literature. We then use the Pleiades LDB age and photometry in combination with evolutionary models to define empirical colour-T_eff relations that can be applied to younger clusters. We find that these relationships DO NOT produce model isochrones that match the younger cluster data. We propose that this is due either to systematic problems in the evolutionary models or an age (gravity) sensitivity in the colour-T_eff relation which is not present in published atmospheric models.

We present high resolution spectroscopy of a sample of 24 solar-type stars in the young (15-40 Myr), open cluster, NGC 2547. We use our spectra to confirm cluster membership in 23 of these stars, determine projected equatorial velocities and chromospheric activity, and to search for the presence of accretion discs. We have found examples of both fast (vsini>50kms) and slow (vsini<10kms) rotators, but find no evidence for active accretion in any of the sample. The distribution of projected rotation velocities is indistinguishable from the slightly older IC 2391 and IC 2602 clusters, implying similar initial angular momentum distributions and circumstellar disc lifetimes. The presence of very slow rotators indicates that either long (10-40 Myr) disc lifetimes or internal differential rotation are needed, or that NGC 2547 (and IC 2391/2602) were born with more slowly rotating stars than are presently seen in even younger clusters and associations. The solar-type stars in NGC 2547 follow a similar rotation-activity relationship to that seen in older clusters. X-ray activity increases until a saturation level is reached for vsini>15-20kms. We are unable to explain why this saturation level, of log (L_x/L_bol)~-3.3, is a factor of two lower than in other clusters, but rule out anomalously slow rotation rates or uncertainties in X-ray flux calculations.

We provide a review of the current status of several topics on the ages, distances, and mass functions of open clusters, with a particular emphasis on illuminating the areas of uncertainty. Hipparcos has obtained parallaxes for nearby open clusters that have expected accuracies much better than has been previously achievable. By using the lithium depletion boundary method and isochrone fitting based on much improved new theoretical evolutionary models for low mass stars, it is arguable that we will soon have have much better age scales for clusters and star-forming regions. With improved optical and near-IR cameras, we are just now beginning to extend the mass function of open clusters like the Pleiades into the regime below the hydrogen burning mass limit. Meanwhile, observations in star-forming regions are in principle capable of identifying objects down to of order 10 Jupiter masses.

(abridged abstract) We present multi-epoch high-resolution spectroscopy and photoelectric polarimetry of the long-period polar (AM Herculis star) QQ Vul. The blue emission lines show several distinct components, the sharpest of which can unequivocally be assigned to the illuminated hemisphere of the secondary star and used to trace its orbital motion. This narrow emission line can be used in combination with NaI-absorption lines from the photosphere of the companion to build a stable long-term ephemeris for the star: inferior conjunction of the companion occurs at HJD = 2448446.4710(5) + E 0.15452011(11).

First principles electronic structure calculations based on density functional theory have been used to study the thermodynamic, structural and transport properties of solid solutions and liquid alloys of iron and oxygen at Earth's core conditions. Aims of the work are to determine the oxygen concentration needed to account for the inferred density in the outer core, to probe the stability of the liquid against phase separation, to interpret the bonding in the liquid, and to find out whether the viscosity differs significantly from that of pure liquid iron at the same conditions. It is shown that the required concentration of oxygen is in the region 25-30 mol percent, and evidence is presented for phase stability at these conditions. The Fe-O bonding is partly ionic, but with a strong covalent component. The viscosity is lower than that of pure liquid iron at Earth's core conditions. It is shown that earlier first-principles calculations indicating very large enthalpies of formation of solid solutions may need reinterpretation, since the assumed crystal structures are not the most stable at the oxygen concentration of interest.

First principles electronic structure calculations, based upon density functional theory within the generalized gradient approximation and ultra-soft Vanderbilt pseudopotentials, have been used to simulate a liquid alloy of iron and sulfur at Earth's core conditions. We have used a sulfur concentration of $\approx 12 % $wt, in line with the maximum recent estimates of the sulfur abundance in the Earth's outer core. The analysis of the structural, dynamical and electronic structure properties has been used to report on the effect of the sulfur impurities on the behavior of the liquid. Although pure sulfur is known to form chains in the liquid phase, we have not found any tendency towards polymerization in our liquid simulation. Rather, a net S-S repulsion is evident, and we propose an explanation for this effect in terms of the electronic structure. The inspection of the dynamical properties of the system suggests that the sulfur impurities have a negligible effect on the viscosity of Earth's liquid core.

Ab initio molecular dynamics simulation is used to study the structure and electronic properties of the liquid Ga-Se system at the three compositions Ga$_2$Se, GaSe and Ga$_2$Se$_3$, and of the GaSe and Ga$_2$Se$_3$ crystals. The calculated equilibrium structure of GaSe crystal agrees well with available experimental data. The neutron-weighted liquid structure factors calculated from the simulations are in reasonable agreement with recent neutron diffraction measurements. Simulation results for the partial radial distribution functions show that the liquid structure is closely related to that of the crystals. A close similarity between solid and liquid is also found for the electronic density of states and charge density. The calculated electronic conductivity decreases strongly with increasing Se content, in accord with experimental measurements.

First-principles calculations based on density functional theory and the pseudopotential method have been used to investigate the energetics of H$_2$O adsorption on the (110) surface of TiO$_2$ and SnO$_2$. Full relaxation of all atomic positions is performed on slab systems with periodic boundary conditions, and the cases of full and half coverage are studied. Both molecular and dissociative (H$_2$O $\rightarrow$ OH$^-$ + H$^+$) adsorption are treated, and allowance is made for relaxation of the adsorbed species to unsymmetrical configurations. It is found that for both TiO$_2$ and SnO$_2$ an unsymmetrical dissociated configuration is the most stable. The symmetrical molecularly adsorbed configuration is unstable with respect to lowering of symmetry, and is separated from the fully dissociated configuration by at most a very small energy barrier. The calculated dissociative adsorption energies for TiO$_2$ and SnO$_2$ are in reasonable agreement with the results of thermal desorption experiments. Calculated total and local electronic densities of states for dissociatively and molecularly adsorbed configurations are presented and their relation with experimental UPS spectra is discussed.