title: Evolution of gravitational orbits in the expanding universe

abstract: The gravitational action of the smooth energy-matter components filling in the universe can affect the orbit of a planetary system. Changes are related to the acceleration of the cosmological scale size R. In a universe with significant dark matter, a gravitational system expands or contracts according to the amount and equation of state of the dark energy. At present time, the Solar system, according to the LambdaCDM scenario emerging from observational cosmology, should be expanding if we consider only the effect of the cosmological background. Its fate is determined by the equation of state of the dark energy alone. The mean motion and periastron precession of a planet are directly sensitive to (d^2 R/d t^2)/R, whereas variations with time in the semi-major axis and eccentricity are related to its time variation. Actual bounds on the cosmological deceleration parameters q_0 from accurate astrometric data of perihelion precession and changes in the third Kepler's law in the Solar system fall short of ten orders of magnitude with respect to estimates from observational cosmology. Future radio-ranging measurements of outer planets could improve actual bounds by five orders of magnitude.


title: The s-process in stellar population synthesis: a new approach to understanding AGB stars

abstract: Thermally pulsating asymptotic giant branch (AGB) stars are the main producers of slow neutron capture (s-) process elements, but there are still large uncertainties associated with the formation of the main neutron source, 13C, and with the physics of these stars in general. Observations of s-process element enhancements in stars can be used as constraints on theoretical models. For the first time we apply stellar population synthesis to the problem of s-process nucleosynthesis in AGB stars, in order to derive constraints on free parameters describing the physics behind the third dredge-up and the properties of the neutron source. We utilize a rapid evolution and nucleosynthesis code to synthesize different populations of s-enhanced stars, and compare them to their observational counterparts to find out for which values of the free parameters in the code the synthetic populations fit best to the observed populations. These free parameters are the amount of third dredge-up, the minimum core mass for third dredge-up, the effectiveness of 13C as a source of neutrons and the size in mass of the 13C pocket. We find that galactic disk objects are reproduced by a spread of a factor of two in the effectiveness of the 13C neutron source. Lower metallicity objects can be reproduced only by lowering by at least a factor of 3 the average value of the effectiveness of the 13C neutron source needed for the galactic disk objects. Using observations of s-process elements in post-AGB stars as constraints we find that dredge-up has to start at a lower core mass than predicted by current theoretical models, that it has to be substantial ($\lambda$ >~ 0.2) in stars with mass M <~ 1.5 M_sun and that the mass of the 13C pocket must be about 1/40 that of the intershell region.


title: The Oblate Schwarzschild Approximation for Light Curves of Rapidly Rotating Neutron Stars

abstract: We present a simple method for including the oblateness of a rapidly rotating neutron star when fitting X-ray light curves. In previous work we showed that the oblateness induced by rotation at frequencies above 300 Hz produces a geometric effect which needs to be accounted for when modelling light curves to extract constraints on the neutron star's mass and radius. In our model X-rays are emitted from the surface of an oblate neutron star and propagate to the observer along geodesics of the Schwarzschild metric for a spherical neutron star. Doppler effects due to rotation are added in the same manner as in the case of a spherical neutron star. We show that this model captures the most important effects due to the neutron star's rotation. We also explain how the geometric oblateness effect can rival the Doppler effect for some emission geometries.


title: The initial mass function of simple and composite stellar populations

abstract: (abridged) The case can be made for a rather universal stellar IMF form that can be approximated by a two-part power-law function in the stellar regime. However, there exists a possible hint for a systematic variation with metallicity. A picture is emerging according to which the binary properties of very-low-mass stars (VLMSs) and BDs may be fundamentally different from those of late-type stars implying the probable existence of a discontinuity in the IMF, but the surveys also appear to suggest the number of BDs per star to be independent of the physical conditions of current Galactic star formation. Star-burst clusters and thus globular cluster may, however, have a much larger abundance of BDs. Very recent advances have allowed the measurement of the physical upper stellar mass limit, which also appears to be disconcertingly robust to variations in metallicity. Furthermore, it now appears that star clusters may be formed such that the most-massive stars just forming terminate further star-formation within the particular cluster. Populations formed from many star clusters, composite populations, would then have steeper IMFs (fewer massive stars per low-mass star) than the simple populations in the constituent clusters. A near invariant star-cluster mass function implies the maximal cluster mass to correlate with the galaxy-wide star-formation rate. This then leads to the result that the composite-stellar IMFs vary in dependence of galaxy type, with potentially dramatic implications for theories of galaxy formation and evolution.


title: Discovery of a Very Bright, Nearby Gravitational Microlensing Event

abstract: We report the serendipitous detection of a very bright, very nearby microlensing event. In late October 2006, an otherwise unremarkable A0 star at a distance ~1 kpc (GSC 3656-1328) brightened achromatically by a factor of nearly 40 over the span of several days and then decayed in an apparently symmetrical way. We present a light curve of the event based on optical photometry from the Center for Backyard Astrophysics and the All Sky Automated Survey, as well as near-infrared photometry from the Peters Automated Infrared Imaging Telescope. This light curve is well-fit by a generic microlensing model. We also report optical spectra, and Swift X-ray and UV observations that are consistent with the microlensing interpretation. We discuss and reject alternative explanations for this variability. The lens star is probably a low-mass star or brown dwarf, with a relatively high proper motion of >20 mas/yr, and may be visible using precise optical/infrared imaging taken several years from now. A modest, all-sky survey telescope could detect ~10 such events per year, which would enable searches for very low-mass planetary companions to relatively nearby stars.


title: The Observed Concentration-Mass Relation for Galaxy Clusters

abstract: The properties of clusters of galaxies offer key insights into the assembly process of structure in the universe. Numerical simulations of cosmic structure formation in a hierarchical, dark matter dominated universe suggest that galaxy cluster concentrations, which are a measure of a halo's central density, decrease gradually with virial mass. However, cluster observations have yet to confirm this correlation. The slopes of the run of measured concentrations with virial mass are often either steeper or flatter than predicted by simulations. In this work, we present the most complete sample of observed cluster concentrations and masses yet assembled, including new measurements for 10 strong lensing clusters, thereby more than doubling the existing number of strong lensing concentration estimates. We fit a power law to the observed concentrations as a function of virial mass, and find that the slope is consistent with the slopes found in simulations, though our normalization factor is higher. Observed lensing concentrations appear to be systematically larger than X-ray concentrations, a more pronounced effect than found in simulations. We also find that at fixed mass, the bulk of observed cluster concentrations are distributed log-normally, with the exception of a few anomalously high concentration clusters. We examine the physical processes likely responsible for the discrepancy between lensing and X-ray concentrations, and for the anomalously high concentrations in particular. The forthcoming Millennium simulation results will offer the most comprehensive comparison set to our findings of an observed concentration-mass power law relation.


title: Comparison of different spectral population models

abstract: We have compared simple stellar populations (SSPs) generated with different population synthesis tools: BC03, Vazdekis and Pegase.HR and different stellar libraries: ELODIE3.1, SteLib and MILES. We find that BC03/SteLib SSPs are biased toward solar metallicity, however Pegase.HR/ELODIE3.1 and Vazdekis/MILES are extremely consistent. The extensive coverage of the space of atmospheric parameters in the large stellar libraries allows precise synthesis for a large range of ages (0.1 .. 10 Gyr) and metallicities (-2 .. +0.4 dex) limited by the quality of the determination of stellar parameters (like temperature scale of the giants).


title: Spectral modelling of the high energy emission of the magnetar 4U 0142+614

abstract: We present an empirical spectral modelling of the high energy emission of the anomalous X-ray pulsar 4U 0142+614, based on simultaneous Swift and INTEGRAL observations from X to gamma-ray energies. We adopted models contained in the XSPEC analysis package, as well as models based on recent theoretical studies, and restricted ourselves to those combinations of up to three components which produce a good fit while requiring the lowest number of free parameters. Only three models were found to fit satisfactorily the 0.5-250keV spectrum of 4U 0142+614: i) a ~0.4keV blackbody and two power-laws, ii) a resonant cyclotron scattering model plus a power-law and iii) two log-parabolic functions. We found that only the latter two models do not over-predict the infrared/optical emission observed simultaneously from this AXP, and only the log-parabolic functions can naturally account for the upper limits set by COMPTEL in the gamma-ray range. A possible interpretation of the two log-parabolae in terms of inverse Compton scattering of soft X-ray photons by very energetic particles is discussed.


title: Some Theoretical Implications of Short-Hard Gamma-Ray Burst observations

abstract: Short-hard and long-soft gamma-ray bursts (GRBs) are two distinct phenomena, but their prompt and afterglow emission show many similarities. This suggests that two different progenitor systems lead to similar physical processes and that the prompt and afterglow observations of short-hard GRBs (SHBs) can be examined using models of long GRBs. Here, I discuss three conclusions that can be drawn from SHB observations. I show that the lower limit on the Lorentz factor of SHBs is typically ``only'' 10-50, significantly lower than that of long GRBs. SHBs with observed X-ray afterglow after 1 day are found to be roughly as efficient as long GRBs in converting the outflow energy into prompt gamma-rays. Finally, I examine the origin of SHBs with X-ray dark afterglows and find that the most plausible explanation is that these SHBs exploded in extremely low density environment ($n \lesssim 10^{-5} cm^{-3}$)


title: Analysis of stellar populations with large empirical libraries at high spectral resolution

abstract: The stellar population models dramatically progressed with the arrival of large and complete libraries, ELODIE, CFLIB (=Indo-US) and MILES at a relatively high resolution. We show that the quality of the fits is not anymore limited by the size of the stellar libraries in a large range of ages (0.1 to 10 Gyrs) and metallicities (-2 to +0.4 dex). The main limitations of the empirical stellar libraries are (i) the coverage of the parameters space (lack of hot stars of low metallicity), (ii) the precision and homogeneity of the atmospheric parameters and (iii) the non-resolution of individual element abundances (in particular [$\alpha$/Fe]). Detailed abundance measurements in the large libraries, and usage of theoretical libraries are probably the next steps, and we show that a combination between an empirical (ELODIE) and a theoretical library (Coelho et al. 2005) immediately improves the modeling of ($\alpha$-enhanced) globular clusters.


title: Black Hole Spin Evolution: Implications for Short-hard Gamma Ray Bursts and Gravitational Wave Detection

abstract: The evolution of the spin and tilt of black holes in compact black hole - neutron star and black hole - black hole binary systems is investigated within the framework of the coalescing compact star binary model for short gamma ray bursts via the population synthesis method. Based on recent results on accretion at super critical rates in slim disk models, estimates of natal kicks, and the results regarding fallback in supernova models, we obtain the black hole spin and misalignment. It is found that the spin parameter, a_spin}, is less than 0.5 for initially non rotating black holes and the tilt angle, i_tilt, is less than 45 deg for 50% of the systems in black hole - neutron star binaries. Upon comparison with the results of black hole - neutron star merger calculations we estimate that only a small fraction (~ 0.01) of these systems can lead to the formation of a torus surrounding the coalesced binary potentially producing a short-hard gamma ray burst. On the other hand, for high initial black hole spin parameters (a_spin>0.6) this fraction can be significant (~ 0.4). It is found that the predicted gravitational radiation signal for our simulated population does not significantly differ from that for non rotating black holes. Due to the (i) insensitivity of signal detection techniques to the black hole spin and the (ii) predicted overall low contribution of black hole binaries to the signal we find that the detection of gravitational waves are not greatly inhibited by current searches with non spinning templates. It is pointed out that the detection of a black hole - black hole binary inspiral system with LIGO or VIRGO may provide a direct measurement of the initial spin of a black hole.


title: Measurement of the pressure dependence of air fluorescence emission induced by electrons

abstract: The fluorescence detection of ultra high energy (> 10^18 eV) cosmic rays requires a detailed knowledge of the fluorescence light emission from nitrogen molecules, which are excited by the cosmic ray shower particles along their path in the atmosphere. We have made a precise measurement of the fluorescence light spectrum excited by MeV electrons in dry air. We measured the relative intensities of 34 fluorescence bands in the wavelength range from 284 to 429 nm with a high resolution spectrograph. The pressure dependence of the fluorescence spectrum was also measured from a few hPa up to atmospheric pressure. Relative intensities and collisional quenching reference pressures for bands due to transitions from a common upper level were found in agreement with theoretical expectations. The presence of argon in air was found to have a negligible effect on the fluorescence yield. We estimated that the systematic uncertainty on the cosmic ray shower energy due to the pressure dependence of the fluorescence spectrum is reduced to a level of 1% by the AIRFLY results presented in this paper.


title: The Solar Neighborhood. XVIII. Discovery of new proper motion stars with 0.40 arcsec yr^-1 > mu >= 0.18 arcsec yr^-1 between declinations -90 degrees and -47 degrees

abstract: We report 1606 new proper motion systems in the southern sky (declinations -90 degrees to -47 degrees with 0.40 arcsec yr^-1 > mu >= 0.18 yr^-1. This effort is a continuation of the SuperCOSMOS-RECONS (SCR) proper motion search to lower proper motions than reported in Papers VIII, X, XII, and XV in this series. Distance estimates are presented for the new systems, assuming that all stars are on the main sequence. We find that 31 systems are within 25 pc, including two systems -- SCR 0838-5855 and SCR 1826-6542 -- we anticipate to be within 10 pc. These new discoveries constitute a more than ten-fold increase in new systems found in the same region of sky searched for systems with mu >= 0.40 arcsec yr^-1, suggesting a happy hunting ground for new nearby slower proper motion systems in the region just north (declinations -47 degrees to 0 degrees, much of which has not been rigorously searched during previous efforts.


title: Binaries in the Kuiper Belt

abstract: Binaries have played a crucial role many times in the history of modern astronomy and are doing so again in the rapidly evolving exploration of the Kuiper Belt. The large fraction of transneptunian objects that are binary or multiple, 48 such systems are now known, has been an unanticipated windfall. Separations and relative magnitudes measured in discovery images give important information on the statistical properties of the binary population that can be related to competing models of binary formation. Orbits, derived for 13 systems, provide a determination of the system mass. Masses can be used to derive densities and albedos when an independent size measurement is available. Angular momenta and relative sizes of the majority of binaries are consistent with formation by dynamical capture. The small satellites of the largest transneptunian objects, in contrast, are more likely formed from collisions. Correlations of the fraction of binaries with different dynamical populations or with other physical variables have the potential to constrain models of the origin and evolution of the transneptunian population as a whole. Other means of studying binaries have only begun to be exploited, including lightcurve, color, and spectral data. Because of the several channels for obtaining unique physical information, it is already clear that binaries will emerge as one of the most useful tools for unraveling the many complexities of transneptunian space.


title: A Deficit of Faint Red Galaxies in the Possible Large-Scale Structures around the RDCS J1252.9-2927 Cluster at z=1.24

abstract: (Abridged) We report a discovery of possible large-scale structures around the RDCS J1252.9-2927 cluster at z=1.24 based on photometric redshifts. We carried out multi-band wide-field imaging with Suprime-Cam on the Subaru Telescope and WFCAM on the United Kingdom Infra-Red Telescope (UKIRT). The distribution of photo-z selected galaxies reveals clumpy structures surrounding the central cluster. We compare the observed structure with an X-ray map and find that two of the four plausible clumps show significant X-ray emissions and one with a marginal detection, which strongly suggest that they are dynamically bound systems. Following the discovery of the possible large-scale structure, we carried out deeper SOFI K_s-band imaging with New Technology Telescope on the four plausible clumps. We construct the optical-to-near-infrared colour-magnitude diagrams of the galaxies in the clumps, and find that the colour-magnitude relation (CMR) of the red galaxies in the clumps is sharply truncated below K_s=22. Interestingly, the main cluster shows a clear relation down to K_s=23 (Lidman et al. 2004). We suggest that galaxies follow the 'environment-dependent down-sizing' evolution. Massive galaxies in high density environments first stop forming stars and become red. Less massive galaxies in less dense environments become red at later times. Based on a few assumptions, we predict that the brightest tip of the CMR appears at z~2.5.


title: High Efficiency of Gamma-Ray Bursts Revisited

abstract: Using the conservation of energy and momentum during collisions of any two shells, we consider the efficiency of gamma-ray bursts by assuming that the ejecta from the central engine are equally massive and have the same Lorentz factors. We calculate the efficiency and the final Lorentz factor of the merged whole shell for different initial diversities of Lorentz factors and for different microscopic radiative efficiency. As a result, a common high efficiency in the range of 0.1 to 0.9 is considerable, and a very high value near 100% is also reachable if the diversity of the Lorentz factors is large enough.


title: Eclipsing binaries in the MACHO database: New periods and classifications for 3031 systems in the Large Magellanic Cloud

abstract: Eclipsing binaries offer a unique opportunity to determine fundamental physical parameters of stars using the constraints on the geometry of the systems. Here we present a reanalysis of publicly available two-color observations of about 6800 stars in the Large Magellanic Cloud, obtained by the MACHO project between 1992 and 2000 and classified as eclipsing variable stars. Of these, less than half are genuine eclipsing binaries. We determined new periods and classified the stars, 3031 in total, using the Fourier parameters of the phased light curves. The period distribution is clearly bimodal, reflecting refer to the separate groups of more massive blue main sequence objects and low mass red giants. The latter resemble contact binaries and obey a period-luminosity relation. Using evolutionary models, we identified foreground stars. The presented database has been cleaned of artifacts and misclassified variables, thus allowing searches for apsidal motion, tertiary components, pulsating stars in binary systems and secular variations with time-scales of several years.


title: Populating the Galaxy with close Be binaries

abstract: Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of how the evolution changes when the star is in close proximity to a companion. Within this work we complete a full population synthesis study of Be binaries for the Galaxy. The results for the first time match aspects of the observational data, most notably the observed upper limit to the period distribution. We conclude that greater detailed studies on the evolution of Be stars within X-ray binaries needs to be completed, so that rapid binary evolution population synthesis packages may best evolve these systems.


title: Evolution of asymptotic giant branch stars I. Updated synthetic TP-AGB models and their basic calibration

abstract: We present new synthetic models of the TP-AGB evolution. They are computed for 7 values of initial metal content (Z from 0.0001 to 0.03) and for initial masses between 0.5 and 5.0 Msun, thus extending the low- and intermediate-mass tracks of Girardi et al. (2000) until the beginning of the post-AGB phase. The calculations are performed by means of a synthetic code that incorporates many recent improvements, among which we mention: (1) the use of detailed and revised analytical relations to describe the evolution of quiescent luminosity, inter-pulse period, third dredge-up, hot bottom burning, pulse cycle luminosity variations, etc.; (2) the use of variable molecular opacities -- i.e. opacities consistent with the changing photospheric chemical composition -- in the integration of a complete envelope model, instead of the standard choice of scaled-solar opacities; (3) the use of formalisms for the mass-loss rates derived from pulsating dust-driven wind models of C- and O-rich AGB stars; and (4) the switching of pulsation modes between the first overtone and the fundamental one along the evolution, which has consequences in terms of the history of mass loss. It follows that, in addition to the time evolution on the HR diagram, the new models predict in a consistent fashion also variations in surface chemical compositions, pulsation modes and periods, and mass-loss rates. The onset and efficiency of the third dredge-up process are calibrated in order to reproduce basic observables like the carbon star luminosity functions in the Magellanic Clouds, and TP-AGB lifetimes (star counts) in Magellanic Cloud clusters. Forthcoming papers will present the theoretical isochrones and chemical yields derived from these tracks, and additional tests performed with the aid of a complete population synthesis code.


title: Kinematics of gas and stars in the circumnuclear starforming ring of NGC 3351

abstract: We have measured gas and stellar velocity dispersions in 5 circumnuclear starforming regions (CNSFRs) and the nucleus of the barred spiral galaxy NGC 3351. The stellar dispersions have been obtained from high resolution spectra of the CaT lines at $\lambda\lambda$ 8494, 8542, 8662 \AA, while the gas velocity dispersions have been measured by Gaussian fits to the H$\beta$ $\lambda$ 4861 \AA line on high dispersion spectra. The CNSFRs, with sizes of about 100 to 150 pc in diameter, are seen to be composed of several individual star clusters with sizes between 1.7 and 4.9 pc on an HST image. Using the stellar velocity dispersions, we have derived dynamical masses for the entire starforming complexes and for the individual star clusters. Values of the stellar velocity dispersions are between 39 and 67 km s$^{-1}$. Dynamical masses for the whole CNSFRs are between 4.9 $\times$ 10$^6$ and 4.3 $\times$ 10$^7$ M$_\odot$ and between 1.8 and 8.7 $\times$ 10$^6$ M$_\odot$ for the individual star clusters. Stellar and gas velocity dispersions are found to differ by about 20 km s$^{-1}$ with the H$\beta$ lines being narrower than both the stellar lines and the [O{\sc iii}] $\lambda$ 5007 \AA lines. We have found indications for the presence of two different kinematical components in the ionised gas of the regions. The radial velocity curve shows deviation from circular motions for the ionised hydrogen consistent with its infall towards the central regions of the galaxy at a velocity of about 25 km s$^{-1}$. To disentangle the origin of these two components it will be necessary to map these regions with high spectral and spatial resolution and much better S/N in particular for the O$^{2+}$ lines.


title: Simulations of partially coherent focal plane imaging arrays: Fisher matrix approach to performance evaluation

abstract: Focal plane arrays of bolometers are increasingly employed in astronomy at far--infrared to millimetre wavelengths. The focal plane fields and the detectors are both partially coherent in these systems, but no account has previously been taken of the effect of partial coherence on array performance. In this paper, we use our recently developed coupled--mode theory of detection together with Fisher information matrix techniques from signal processing to characterize the behaviour of partially coherent imaging arrays. We investigate the effects of the size and coherence length of both the source and the detectors, and the packing density of the array, on the amount of information that can be extracted from observations with such arrays.


title: The main- interpulse interaction of PSR B1702-19

abstract: This paper reports on single-pulse radio observations of PSR B1702-19 and their implications for pulsar emission theories. These observations were made with the Westerbork Synthesis Radio Telescope at 1380 and 328 MHz. The PA-swing is used to constrain possible geometries of the pulsar and the single-pulse data is analysed for subpulse modulation correlations between the main pulse and interpulse. We confirm earlier conclusions that the dipole axis of this pulsar is almost perpendicular to its rotation axis, and report that both its main pulse and interpulse are modulated with a periodicity around 10.4 times the pulsar's rotation. Allowing for the half-period delay between main pulse and interpulse the modulation is found to be precisely in phase. Despite small secular variations in the periodicity, the phase-locking continues over all timescales ranging up to several years. The precision of the phase locking is difficult for current emission theories to explain if the main pulse and interpulse originate from opposing magnetic poles. We therefore also explore the possibility of a bidirectional model, in which all the modulated emission comes from one pole, but is seen from two sides and slightly displaced by aberration and time-delay. In this model the unmodulated emission is directed to us from the opposite pole, requiring the emission of the main pulse to originate from two different poles. This is difficult to reconcile with the observed smooth PA-swing. Whichever model turns out to be correct, the answer will have important implications for emission theories.


title: On the difference between Herbig Ae and Herbig Be stars

abstract: We present linear spectropolarimetric data for eight Herbig Be and four Herbig Ae stars at H alpha, H beta and H gamma. Changes in the linear polarisation are detected across all Balmer lines for a large fraction of the observed objects, confirming that the small-scale regions surrounding these objects are flattened (i.e. disk-like). Furthermore, all objects with detections show similar characteristics at the three spectral lines, despite differences in transition probability and optical depth going from H alpha to H gamma. A large fraction of early Herbig Be stars (B0-B3) observed show line depolarisation effects. However the early Herbig Ae stars (A0-A2), observed for comparison, show intrinsic line polarisation signatures. Our data suggest that the popular magnetic accretion scenario for T Tauri objects may be extended to Herbig Ae stars, but that it may not be extended to early Herbig Be stars, for which the available data are consistent with disc accretion.


title: Gamma-Ray Burst Precursors as the Remnant of the Thermal Radiation Initially Trapped in the Fireball

abstract: In the standard fireball model of gamma-ray bursts (GRBs), the fireball starts with an optically thick phase. As it expands, the fireball becomes optically thin at some stage. The thermal radiation trapped in the originally opaque fireball then leaks out, producing a transient event. The appearance of the event is investigated in the framework of a homogeneous, spherically symmetric, and freely expanding fireball produced instantly by an explosive process without continuous injection of mass and energy. We find that, generally, the event has a time-duration shorter than that of the main burst, which is presumably produced by the internal shock after the fireball becomes optically thin. The event is separated from the main burst by a quiescent time-interval, and is weaker than the main burst at least in a high energy band. Hence, the event corresponds to a GRB precursor. The precursor event predicted by our model has a smooth and FRED (Fast Rise and Exponential Decay) shape lightcurve, and a quasi-thermal spectrum. Typically, the characteristic blackbody photon energy is in the X-ray band. However, if the distortion of the blackbody spectrum by electron scattering is considered, the characteristic photon energy could be boosted to the gamma-ray band. Our model may explain a class of observed GRB precursors--those having smooth and FRED-shape lightcurves and quasi-thermal spectra.


title: Accurate Ritz wavelengths of parity-forbidden [Fe II], [Ti II] and [Cr II] infrared lines of astrophysical interest

abstract: With new astronomical infrared spectrographs the demands of accurate atomic data in the infrared have increased. In this region there is a large amount of parity-forbidden lines, which are of importance in diagnostics of low-density astrophysical plasmas. We present improved, experimentally determined, energy levels for the lowest even LS terms of Fe II, Ti II and Cr II, along with accurate Ritz wavelengths for parity-forbidden transitions between and within these terms. Spectra of Fe II, Ti II and Cr II have been produced in a hollow cathode discharge lamp and acquired using high-resolution Fourier Transform (FT) spectrometry. The energy levels have been determined by using observed allowed ultraviolet transitions connecting the even terms with upper odd terms. Ritz wavelengths of parity-forbidden lines have then been determined. Energy levels of the four lowest Fe II terms (a$^{6}$D, a$^{4}$F, a$^{4}$D and a$^{4}$P) have been determined, resulting in 97 different parity-forbidden transitions with wavelengths between 0.74 and 87 micron. For Ti II the energy levels of the two lowest terms (a$^{4}$F and b$^{4}$F) have been determined, resulting in 24 different parity-forbidden transitions with wavelengths between 8.9 and 130 micron. Also for Cr II the energy levels of the two lowest terms (a$^{6}$S and a$^{6}$D) have been determined, in this case resulting in 12 different parity-forbidden transitions with wavelengths between 0.80 and 140 micron.


title: Magnetic acceleration of relativistic AGN jets

abstract: We present numerical simulations of axisymmetric, magnetically driven relativistic jets. To eliminate the dissipative effects induced by a free boundary with an ambient medium we assume that the flow is confined by a rigid wall of a prescribed shape, which we take to be $z\propto r^a$ (in cylindrical coordinates, with $a$ ranging from 1 to 3). The outflows are initially cold, sub-Alfv\'enic and Poynting flux-dominated, with a total--to--rest-mass energy flux ratio $\mu \sim 15$. We find that in all cases they converge to a steady state characterized by a spatially extended acceleration region. The acceleration process is very efficient: on the outermost scale of the simulation as much as $\sim 77%$ of the Poynting flux has been converted into kinetic energy flux, and the terminal Lorentz factor approaches its maximum possible value ($\Gamma_\infty \simeq \mu$). We also find a high collimation efficiency: all our simulated jets develop a cylindrical core. We argue that this could be the rule for current-carrying outflows that start with a low initial Lorentz factor ($\Gamma_0 \sim 1$). Our conclusions on the high acceleration and collimation efficiencies are not sensitive to the particular shape of the confining boundary or to the details of the injected current distribution, and they are qualitatively consistent with the semi-analytic self-similar solutions derived by Vlahakis & K\"onigl. We apply our results to the interpretation of relativistic jets in AGNs: we argue that they naturally account for the spatially extended accelerations inferred in these sources ($\Gamma_\infty \ga 10$ attained on radial scales $R\ga 10^{17} {\rm cm}$) and are consistent with the transition to the matter-dominated regime occurring already at $R\ga 10^{16} {\rm cm}$.


title: Reconstruction of solar irradiance using the Group sunspot number

abstract: We present a reconstruction of total solar irradiance since 1610 to the present based on variations of the surface distribution of the solar magnetic field. The latter is calculated from the historical record of the Group sunspot number using a simple but consistent physical model. Our model successfully reproduces three independent data sets: total solar irradiance measurements available since 1978, total photospheric magnetic flux from 1974 and the open magnetic flux since 1868 (as empirically reconstructed from the geomagnetic aa-index). The model predicts an increase in the total solar irradiance since the Maunder Minimum of about 1.3 \rm{Wm$^{-2}$}.


title: A New Numerical Method for Solving Radiation Driven Winds from Hot Stars

abstract: We present a general method for solving the non--linear differential equation of monotonically increasing steady--state radiation driven winds. We graphically identify all the singular points before transforming the momentum equation to a system of differential equations with all the gradients explicitly give. This permits a topological classification of all singular points and to calculate the maximum and minimum mass--loss of the wind. We use our method to analyse for the first time the topology of the non--rotating frozen in ionisation m--CAK wind, with the inclusion of the finite disk correction factor and find up to 4 singular points, three of the x--type and one attractor--type. The only singular point (and solution passing through) that satisfies the boundary condition at the stellar surface is the standard m--CAK singular point.


title: Supermassive black holes in the Sbc spiral galaxies NGC 3310, NGC 4303 and NGC 4258

abstract: We present new Space Telescope Imaging Spectrograph (STIS) observations of three spiral galaxies, NGC 4303, NGC 3310 and NGC 4258. The bright optical emission lines H$\alpha$ $\lambda$ $6564 \AA$, [NII] $\lambda$$\lambda$ $6549,6585 \AA$ and [SII] $\lambda$$\lambda$ $ 6718,6732 \AA$ were used to study the kinematics of the ionized gas in the nuclear region of each galaxy with a $\sim 0.07\arcsec$ spatial resolution. In NGC 3310, the observed gas kinematics is well matched by a circularly rotating disk model but we are only able to set an upper limit to the BH mass which, taking into account the allowed disk inclinations, varies in the range $5.0 \times 10^{6} - 4.2 \times 10^{7} M_{\odot}$ at the 95% confidence level. In NGC 4303 the kinematical data require the presence of a BH with mass $M_{BH}=(5.0)^{+0.87}_{-2.26}\times 10^{6}M_{\odot}$ (for a disk inclination $i=70$ deg).In NGC 4258, the observed kinematics require the presence of a black hole with $M_{BH}= (7.9)^{+6.2}_{-3.5} \times 10^{7}M_{\odot}$ ($i=60$ deg). This result is in good agreement with the published value $(3.9 \pm 0.1) \times 10^{7} M_{\odot}$, derived from $H_{2}O$-maser observations. Our attempt at measuring BH masses in these 3 late type Sbc spiral galaxies has shown that these measurements are very challenging and at the limit of the highest spatial resolution currently available. Nonetheless our estimates are in good agreement with the scaling relations between black holes and their host spheroids suggesting that (i) they are reliable and (ii) black holes in spiral galaxies follows the same scaling relations as those in more massive early-type galaxies. A crucial test for the gas kinematical method, the correct recovery of the known BH mass in NGC 4258, has been successful. [abridged]


title: Swift detection of all previously undetected blazars in a micro-wave flux-limited sample of WMAP foreground sources

abstract: Almost the totality of the bright foreground sources in the WMAP CMB maps are blazars, a class of sources that show usually also X-ray emission. However, 23 objects in a flux-limited sample of 140 blazars of the WMAP catalog (first year) were never reported before as X-ray sources. We present here the results of 41 Swift observations which led to the detection of all these 23 blazars in the 0.3-10 keV band. We conclude that all micro-wave selected blazars are X-ray emitters and that the distribution of the micro-wave to X-ray spectral slope $\alpha_{mu x}$ of LBL blazars is very narrow, confirming that the X-ray flux of most blazars is a very good estimator of their micro-wave emission. The X-ray spectral shape of all the objects that were observed long enough to allow spectral analysis is flat and consistent with inverse Compton emission within the commonly accepted view where the radiation from blazars is emitted in a Sychrotron-Inverse-Compton scenario. We predict that all blazars and most radio galaxies above the sensitivity limit of the WMAP and of the Planck CMB missions are X-ray sources detectable by the present generation of X-ray satellites. An hypothetical all-sky soft X-ray survey with sensitivity of approximately $10^{-15}$ erg/s would be crucial to locate and remove over 100,000 blazars from CMB temperature and polarization maps and therefore accurately clean the primordial CMB signal from the largest population of extragalactic foreground contaminants.


title: Cavity opening by a giant planet in a protoplanetary disc and effects on planetary migration

abstract: We study the effect of a Jovian planet on the gas distribution of a protoplanetary disc, using a new numerical scheme that allows us to take into consideration the global evolution of the disc, down to an arbitrarily small inner physical radius. We find that Jovian planets do not open cavities in the inner part of the disc (i.e. interior to their orbits) unless (a) the inner physical edge of the disc is close to the planet's location or (b) the planet is much more massive than the disc. In all other cases the planet simply opens a gap in the gas density distribution, whose global profile is essentially unchanged relative to the one that it would have if the planet were absent. We recognize, though, that the dust distribution can be significantly different from the gas distribution and that dust cavities might be opened in some situations, even if the gas is still present in the inner part of the disc. Concerning the migration of the planet, we find that classical Type-II migration (with speed proportional to the viscosity of the disc) occurs only if the gap opened by the planet is deep and clean. If there is still a significant amount of gas in the gap, the migration of the planet is generally slower than the theoretical Type-II migration rate. In some situations, migration can be stopped or even reversed. We develop a simple model that reproduces satisfactorily the migration rate observed in the simulations, for a wide range of disc viscosities and planet masses and locations relative to the inner disc edge. Our results are relevant for extra-solar planetary systems, as they explain (a) why some hot Jupiters did not migrate all the way down to their parent stars and (b) why the outermost of a pair of resonant planets is typically the most massive one.


title: Protostellar Turbulence Driven by Collimated Outflows

abstract: The majority of stars are thought to form in clusters. Cluster formation in dense clumps of molecular clouds is strongly influenced, perhaps controlled, by supersonic turbulence. We have previously shown that the turbulence in regions of active cluster formation is quickly transformed by the forming stars through protostellar outflows, and that the outflow-driven protostellar turbulence is the environment in which most of the cluster members form. Here, we take initial steps in quantifying the global properties of the protostellar turbulence through 3D MHD simulations. We find that collimated outflows are more efficient in driving turbulence than spherical outflows that carry the same amounts of momentum. Gravity plays an important role in shaping the turbulence, generating infall motions in the cluster forming region that balance the outward motions driven by outflows. The resulting quasi-equilibrium state is maintained through a slow rate of star formation, with a fraction of the total mass converted into stars per free fall time as low as a few percent. Magnetic fields are dynamically important even in magnetically supercritical clumps, provided that their initial strengths are not far below the critical value for static cloud support. We find that the mass weighted PDF of the volume density is often, although not always, approximately lognormal. The PDFs of the column density deviate more strongly from lognormal distributions. There is a prominent break in the velocity power spectrum of the protostellar turbulence, which may provide a way to distinguish it from other types of turbulence.


title: VLT and ACS observations of RDCS J1252.9-2927: dynamical structure and galaxy populations in a massive cluster at z=1.237

abstract: We present results from an extensive spectroscopic survey, carried out with VLT FORS, and from an extensive multiwavelength imaging data set from the HST Advanced Camera for Surveys and ground based facilities, of the cluster of galaxies RDCS J1252.9-2927. We have spectroscopically confirmed 38 cluster members in the redshift range 1.22 < z < 1.25. A cluster median redshift of z=1.237 and a rest-frame velocity dispersion of 747^{+74}_{-84} km/s are obtained. Using the 38 confirmed redshifts, we were able to resolve, for the first time at z > 1, kinematic structure. The velocity distribution, which is not Gaussian at the 95% confidence level, is consistent with two groups that are also responsible for the projected east-west elongation of the cluster. The groups are composed of 26 and 12 galaxies with velocity dispersions of 486^{+47}_{-85} km/s and 426^{+57}_{-105} km/s, respectively. The elongation is also seen in the intracluster gas and the dark matter distribution. This leads us to conclude that RDCS J1252.9-2927 has not yet reached a final virial state. We extend the analysis of the color-magnitude diagram of spectroscopic members to more than 1 Mpc from the cluster center. The scatter and slope of non-[OII]-emitting cluster members in the near-IR red sequence is similar to that seen in clusters at lower redshift. Furthermore, most of the galaxies with luminosities greater than ~ K_s*+1.5 do not show any [OII], indicating that these more luminous, redder galaxies have stopped forming stars earlier than the fainter, bluer galaxies. Our observations provide detailed dynamical and spectrophotometric information on galaxies in this exceptional high-redshift cluster, delivering an in-depth view of structure formation at this epoch only 5 Gyr after the Big Bang.


title: Cosmic-ray positron fraction measurement from 1 to 30 GeV with AMS-01

abstract: A measurement of the cosmic ray positron fraction e+/(e+ + e-) in the energy range of 1-30 GeV is presented. The measurement is based on data taken by the AMS-01 experiment during its 10 day Space Shuttle flight in June 1998. A proton background suppression on the order of 10^6 is reached by identifying converted bremsstrahlung photons emitted from positrons.


title: Proper Motions of the Jets in the Region of HH 30 and HL/XZ Tau. Evidence for a Binary Exciting Source of the HH 30 Jet

abstract: We present [SII] images of the HH 30 and HL/XZ Tau region obtained at two epochs, as well as long-slit optical spectroscopy of the HH 30 jet. We measured proper motions of about 100-300 km/s for the HH 30 jet and counterjet, and of about 120 km/s for the HL Tau jet. Inclination angles with respect to the plane of the sky are 0-40 deg for the HH 30 jet and 60 deg for the HL Tau jet. Comparison with previous observations suggests that most of the jet knots consist of persisting structures. Also, we corroborate that the HH 30-N knots correspond to the head of the HH 30 jet. The overall HH 30 jet structure can be well described by a wiggling ballistic jet, arising either by the orbital motion of the jet source around a primary or by precession of the jet axis because of the tidal effects of a companion. In the first scenario, the orbital period would be 53 yr and the total mass 0.25-2 solar masses. In the precession scenario, the mass of the jet source would be 0.1-1 solar masses, the orbital period <1 yr, and the mass of the companion less than a few times 0.01 solar masses, thus being a substellar object or a giant exoplanet. In both scenarios a binary system with a separation <18 AU (<0.13 arcsec) is required. Since the radius of the flared disk observed with the HST is about 250 AU, we conclude that this disk appears to be circumbinary rather than circumstellar, suggesting that the search for the collimating agent of the HH 30 jet should be carried out at much smaller scales.


title: Images, structural properties and metal abundances of galaxy clusters observed with Chandra ACIS-I at 0.1<z<1.3

abstract: We have assembled a sample of 115 galaxy clusters at 0.1<z<1.3 with archived Chandra ACIS-I observations. We present X-ray images of the clusters and make available region files containing contours of the smoothed X-ray emission. The structural properties of the clusters were investigated and we found a significant absence of relaxed clusters (as determined by centroid shift measurements) at z>0.5. The slope of the surface brightness profiles at large radii were steeper on average by 15% than the slope obtained by fitting a simple beta-model to the emission. This slope was also found to be correlated with cluster temperature, with some indication that the correlation is weaker for the clusters at z>0.5. We measured the mean metal abundance of the cluster gas as a function of redshift and found significant evolution, with the abundances dropping by 50% between z=0.1 and z~1. This evolution was still present (although less significant) when the cluster cores were excluded from the abundance measurements, indicating that the evolution is not solely due to the disappearance of relaxed, cool core clusters (which are known to have enhanced core metal abundances) from the population at z>0.5.


title: SDSS Standard Star Catalog for Stripe 82: the Dawn of Industrial 1% Optical Photometry

abstract: We describe a standard star catalog constructed using multiple SDSS photometric observations (at least four per band, with a median of ten) in the $ugriz$ system. The catalog includes 1.01 million non-variable unresolved objects from the equatorial stripe 82 ($|\delta_{J2000}|<$ 1.266$^\circ$) in the RA range 20h 34m to 4h 00m, and with the corresponding $r$ band (approximately Johnson V band) magnitudes in the range 14--22. The distributions of measurements for individual sources demonstrate that the photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in $ugriz$, respectively (about twice as good as for individual SDSS runs). Several independent tests of the internal consistency suggest that the spatial variation of photometric zeropoints is not larger than $\sim$0.01 mag (rms). In addition to being the largest available dataset with optical photometry internally consistent at the $\sim$1% level, this catalog provides practical definition of the SDSS photometric system. Using this catalog, we show that photometric zeropoints for SDSS observing runs can be calibrated within nominal uncertainty of 2% even for data obtained through 1 mag thick clouds, and demonstrate the existence of He and H white dwarf sequences using photometric data alone. Based on the properties of this catalog, we conclude that upcoming large-scale optical surveys such as the Large Synoptic Survey Telescope will be capable of delivering robust 1% photometry for billions of sources.


title: Dust Echoes from the Ambient Medium of Gamma-Ray Bursts

abstract: Long gamma-ray bursts (GRBs) are likely associated with the collapse of massive stars, which produce dust and are born in dusty environments. Absorption and scattering of ultraviolet/X-ray photons from the prompt, optical flash and afterglow emission of the GRB produce dust echoes. We perform time-dependent calculations of these echoes, accounting for the evolution of the dust grain distribution due to selective grain destruction by the GRB radiation, and for off-axis beaming. We explore cloud configurations of differing density and size -- the echo light curve and spectrum depend on the cloud radius, with larger clouds peaking at longer wavelengths. For a region ~ 3 pc in size with hydrogen density ~ 1000 per cubic centimeter, the echo spectrum peaks at ~ 3.6 microns and ~ 8.8 eV for thermal and scattered components, respectively. Dust echoes should be detectable with the Very Large Telescope up to z ~ 0.1, IRAC onboard the Spitzer Space Telescope up to z ~ 0.2, and NICMOS onboard the Hubble Space Telescope up to z ~ 0.3. Furthermore, the shape of the echo light curve allows one to infer: the jet opening angle; the inclination of the jet axis with respect to the line of sight; the size of the dust-emitting region. For sources with symmetric, bipolar jets, dust echoes exhibit two bumps in the light curve, making them easily distinguishable from the rebrightening due to an underlying supernova.


title: Ground-layer wavefront reconstruction from multiple natural guide stars

abstract: Observational tests of ground layer wavefront recovery have been made in open loop using a constellation of four natural guide stars at the 1.55 m Kuiper telescope in Arizona. Such tests explore the effectiveness of wide-field seeing improvement by correction of low-lying atmospheric turbulence with ground-layer adaptive optics (GLAO). The wavefronts from the four stars were measured simultaneously on a Shack-Hartmann wavefront sensor (WFS). The WFS placed a 5 x 5 array of square subapertures across the pupil of the telescope, allowing for wavefront reconstruction up to the fifth radial Zernike order. We find that the wavefront aberration in each star can be roughly halved by subtracting the average of the wavefronts from the other three stars. Wavefront correction on this basis leads to a reduction in width of the seeing-limited stellar image by up to a factor of 3, with image sharpening effective from the visible to near infrared wavelengths over a field of at least 2 arc minutes. We conclude that GLAO correction will be a valuable tool that can increase resolution and spectrographic throughput across a broad range of seeing-limited observations.


title: Dynamical Origin of Extrasolar Planet Eccentricity Distribution

abstract: We explore the possibility that the observed eccentricity distribution of extrasolar planets arose through planet-planet interactions, after the initial stage of planet formation was complete. Our results are based on ~3250 numerical integrations of ensembles of randomly constructed planetary systems, each lasting 100 Myr. We find that for a remarkably wide range of initial conditions the eccentricity distributions of dynamically active planetary systems relax towards a common final equilibrium distribution, well described by the fitting formula dn ~ e exp[-1/2 (e/0.3)^2] de. This distribution agrees well with the observed eccentricity distribution for e > 0.2, but predicts too few planets at lower eccentricities, even when we exclude planets subject to tidal circularization. These findings suggest that a period of large-scale dynamical instability has occurred in a significant fraction of newly formed planetary systems, lasting 1--2 orders of magnitude longer than the ~1 Myr interval in which gas-giant planets are assembled. This mechanism predicts no (or weak) correlations between semimajor axis, eccentricity, inclination, and mass in dynamically relaxed planetary systems. An additional observational consequence of dynamical relaxation is a significant population of planets (>10%) that are highly inclined (>25deg) with respect to the initial symmetry plane of the protoplanetary disk; this population may be detectable in transiting planets through the Rossiter-McLaughlin effect.


title: Signatures of Planet Formation in Gravitationally Unstable Disks

abstract: In this paper, we calculate simulated scattered light images of a circumstellar disk in which a planet is forming by gravitational instability. The simulated images bear no correlation to the vertically integrated surface density of the disk, but rather trace the density structure in the tenuous upper layers of the disk. Although the density at high altitudes does not bear a direct relation to activity at the midplane, the very existence of structure at high altitudes along with high time variability is an indicator of gravitational instability within the disk. The timescale for variations is much shorter than the orbital period of the planet, which facilitates observation of the phenomenon. Scattered light images may not necessarily be able to tell us where exactly a planet might be forming in a disk, but can still be a useful probe of active planet formation within a circumstellar disk. Although these phenomena are unlikely to be observable by current telescopes, future large telescopes, such as the Giant Magellan Telescope, may be able to detect them.


title: Axisymmetric black hole accretion in the Kerr metric as an autonomous dynamical system

abstract: In a stationary, general relativistic, axisymmetric, inviscid and rotational accretion flow, described within the Kerr geometric framework, transonicity has been examined by setting up the governing equations of the flow as a first-order autonomous dynamical system. The consequent linearised analysis of the critical points of the flow leads to a comprehensive mathematical prescription for classifying these points, showing that the only possibilities are saddle points and centre-type points for all ranges of values of the fixed flow parameters. The spin parameter of the black hole influences the multitransonic character of the flow, as well as some of its specific critical properties. The special case of a flow in the space-time of a non-rotating black hole, characterised by the Schwarzschild metric, has also been studied for comparison and the conclusions are compatible with what has been seen for the Kerr geometric case.


title: Origins of Eccentric Extrasolar Planets: Testing the Planet-Planet Scattering Model

abstract: (Abridged) In planetary systems with two or more giant planets, dynamical instabilities can lead to collisions or ejections through strong planet--planet scattering. Previous studies for simple initial configurations with two equal-mass planets revealed some discrepancies between the results of numerical simulations and the observed orbital elements of extrasolar planets. Here, we show that simulations with two unequal mass planets starting on nearly circular orbits predict a reduced frequency of collisions and a broader range of final eccentricities. The two-planet scattering model can easily reproduce the observed eccentricities with a plausible distribution of planet mass ratios. Further, the two-planet scattering model predicts a maximum eccentricity of about 0.8, independent of the distribution of planet mass ratios, provided that both planets are initially place on nearly circular orbits. This compares favorably with current observations and will be tested by future planet discoveries. The combination of planet--planet scattering and tidal circularization may be able to explain the existence of some giant planets with very short period orbits. Orbital migration due to planet scattering could play an important role in explaining the increased rate of giant planets with very short period orbits. We also re-examine and discuss various possible correlations between eccentricities and other properties of observed extrasolar planets. We demonstrate that the observed distribution of planet masses, orbital periods, and eccentricities can provide constraints for models of planet formation and evolution.


title: Shape, spin and baryon fraction of clusters in the MareNostrum Universe

abstract: The MareNostrum Universe is one of the largest cosmological SPH simulation done so far. It consists of $1024^3$ dark and $1024^3$ gas particles in a box of 500 $h^{-1}$ Mpc on a side. Here we study the shapes and spins of the dark matter and gas components of the 10,000 most massive objects extracted from the simulation as well as the gas fraction in those objects. We find that the shapes of objects tend to be prolate both in the dark matter and gas. There is a clear dependence of shape on halo mass, the more massive ones being less spherical than the less massive objects. The gas distribution is nevertheless much more spherical than the dark matter, although the triaxiality parameters of gas and dark matter differ only by a few percent and it increases with cluster mass. The spin parameters of gas and dark matter can be well fitted by a lognormal distribution function. On average, the spin of gas is 1.4 larger than the spin of dark matter. We find a similar behavior for the spins at higher redshifts, with a slightly decrease of the spin ratios to 1.16 at $z=1.$ The cosmic normalized baryon fraction in the entire cluster sample ranges from $Y_b = 0.94$, at $z=1$ to $Y_b = 0.92$ at $z=0$. At both redshifts we find a slightly, but statistically significant decrease of $Y_b$ with cluster mass.


title: 3D Spectroscopy of Blue Compact Galaxies. Diagnostic Diagrams

abstract: Here we present the analysis of 3D spectroscopic data of three Blue Compact Galaxies (Mrk324, Mrk370, and IIIZw102). Each of the more than 22500 spectra obtained for each galaxy has been fitted by a single gaussian from which we have inferred the velocity dispersion (sigma), the peak intensity (Ipeak), and the central wavelength (lambda_c). The analysis shows that the sigma vs Ipeak diagrams look remarkably similar to those obtained for giant extragalactic HII regions. They all present a supersonic narrow horizontal band that extends across all the range of intensities and that result from the massive nuclear star-forming regions of every galaxy. The sigma vs Ipeak diagrams present also several inclined bands of lower intensity and an even larger sigma, arising from the large galactic volumes that surround the main central emitting knots. Here we also show that the sigma vs lambda_c and lambda_c vs Ipeak diagrams, are powerful tools able to unveil the presence of high and low mass stellar clusters, and thus allow for the possibility of inferring the star formation activity of distant galaxies, even if these are not spatially resolved.


title: Dynamical Outcomes of Planet-Planet Scattering

abstract: Observations in the past decade have revealed extrasolar planets with a wide range of orbital semimajor axes and eccentricities. Based on the present understanding of planet formation via core accretion and oligarchic growth, we expect that giant planets often form in closely packed configurations. While the protoplanets are embedded in a protoplanetary gas disk, dissipation can prevent eccentricity growth and suppress instabilities from becoming manifest. However, once the disk dissipates, eccentricities can grow rapidly, leading to close encounters between planets. Strong planet--planet gravitational scattering could produce both high eccentricities and, after tidal circularization, very short-period planets, as observed in the exoplanet population. We present new results for this scenario based on extensive dynamical integrations of systems containing three giant planets, both with and without residual gas disks. We assign the initial planetary masses and orbits in a realistic manner following the core accretion model of planet formation. We show that, with realistic initial conditions, planet--planet scattering can reproduce quite well the observed eccentricity distribution. Our results also make testable predictions for the orbital inclinations of short-period giant planets formed via strong planet scattering followed by tidal circularization.


title: BVI Photometry and the Luminosity Function of the Globular Cluster M92

abstract: We present new BVI ground-based photometry and VI space-based photometry for the globular cluster M92 (NGC 6341) and examine luminosity functions in B, V, and I containing over 50,000 stars ranging from the tip of the red giant branch to several magnitudes below the main sequence turn off. Once corrected for completeness, the observed luminosity functions agree very well with theoretical models and do not show stellar excesses in any region of the luminosity function. Using reduced chi squared fitting, the new M92 luminosity function is shown to be an excellent match to the previously published luminosity function for M30. These points combine to establish that the "subgiant excess" found in previously published luminosity functions of Galactic globular clusters are due to deficiencies in the stellar models used at that time. Using up to date stellar models results in good agreement between observations and theory. Several statistical methods are presented to best determine the age of M92. These methods prove to be insensitive to the exact choice of metallicity within the published range. Using [Fe/H]=-2.17 to match recent studies we find an age of 14.2 plus or minus 1.2 Gyr for the cluster.


title: Smoothing of an All-sky Survey Map with a Fisher-von Mises Function

abstract: A convolution of the all-sky survey data with a smoothing function is crucial in calculating the smooth surface brightness of the sky survey data. The convolution is usually performed using a spherical version of the convolution theorem. However, a Gaussian function, applicable only in the flat-sky approximation, has been usually adopted as a smoothing kernel. In this paper, we present an exact analytic solution of the spherical-harmonics transformation of a Fisher-von Mises function, the mathematical version of a Gaussian function in spherical space. We also obtain the approximate solutions exp[-l(l + 1)/2k]. The exact and the approximate solutions may be useful when an astrophysical survey map is convolved with a smoothing function of k>1.


title: The origin of steep vertical stellar distribution in the Galactic disk

abstract: Over the past two decades observations have revealed that the vertical density distribution of stars in galaxies near the mid-plane is substantially steeper than the sech^2 function that is expected from an isothermal approximation. However, the physical origin for this has not been explained so far. Here we show that such steep profiles result naturally even within the isothermal regime, on taking into account the gravitational force due to the gas. Due to its low velocity dispersion the gas is concentrated closer to the galactic mid-plane than the stars, and hence it strongly affects the vertical stellar distribution even though its contribution to the total surface density is small. We apply a three-component galactic disk model consisting of gravitationally coupled stars and the HI and H_2 gas, embedded in the dark matter halo, and calculate the vertical density distribution of stars for the Galaxy. The resulting vertical density distribution of stars is shown to be steeper than the sech^2 function, and lies between the sech and an exponential function, in good agreement with observations of galaxies. We also show that a multi-component stellar disk consisting of coupled dwarfs and two populations of giants does not explain the observed steep stellar profiles.


title: Coupling between the 45 Hz Horizontal-Branch Oscillation and the Normal Branch Oscillation in Scorpius X-1

abstract: The observations of the bright persistent neutron star low-mass X-ray binary (LMXB) Sco X-1 performed with the {\it Rossi X-ray Timing Explorer} (RXTE) show a $\sim$ 6 Hz normal-branch oscillation (NBO), a $\sim$ 45 Hz horizontal-branch oscillation (HBO), and twin kHz quasi-periodic oscillations (QPOs) on its normal branch simultaneously. We have found that the fractional amplitude of the HBO corresponding to the NBO phase of high flux is 1.1%, while that of the NBO phase of low flux is undetectable, with a 3$\sigma$ upper limit of 0.4%, implying that the HBO strength varies with the NBO phase in an opposite way to that of the lower kHz QPO previously found, and suggests that the condition for the generation of the HBO is met when the NBO flux is high. The 6 Hz NBO in Sco X-1 connects the 45 Hz HBO and the twin kHz QPO together, showing a unique picture indicating a coupling between the QPOs, which has never been observed in other neutron star LMXBs. We discuss the implications for current models of the 45 Hz HBO, the 6 Hz NBO, and the twin kHz QPOs.


title: High resolution science with high redshift galaxies

abstract: We summarize the high-resolution science that has been done on high redshift galaxies with Adaptive Optics (AO) on the world's largest ground-based facilities and with the Hubble Space Telescope (HST). These facilities complement each other. Ground-based AO provides better light gathering power and in principle better resolution than HST, giving it the edge in high spatial resolution imaging and high resolution spectroscopy. HST produces higher quality, more stable PSF's over larger field-of-view's in a much darker sky-background than ground-based AO, and yields deeper wide-field images and low-resolution spectra than the ground. Faint galaxies have steadily decreasing sizes at fainter fluxes and higher redshifts, reflecting the hierarchical formation of galaxies over cosmic time. HST has imaged this process in great structural detail to z<~6, and ground-based AO and spectroscopy has provided measurements of their masses and other physical properties with cosmic time. Last, we review how the 6.5 meter James Webb Space Telescope (JWST) will measure First Light, reionization, and galaxy assembly in the near--mid-IR after 2013.


title: Rapid magnetic flux variability on the flare star CN Leonis

abstract: We present UVES/VLT observations of the nearby flare star CN Leo covering the Wing-Ford FeH band near 1mu with high spectral resolution. Some of the FeH absorption lines in this band are magnetically sensitive and allow a measurement of the mean magnetic flux on CN Leo. Our observations, covering three nights separated by 48 hours each, allow a clear detection of a mean magnetic field of Bf~2.2kG. The differential flux measurements show a night-to-night variability with extremely high significance. Finally, our data strongly suggest magnetic flux variability on time scales as low as 6 hours in line with chromospheric variability.


title: Occultation measurement of the size of the X-ray emitting region in the Active Galactic Nucleus of NGC 1365

abstract: We present an occultation of the central X-ray emitting region in the Seyfert Galaxy NGC 1365. This extreme spectral variation (from Compton-thin to reflection-dominated and back to Compton-thin in four days) has been caught in a ten days Chandra monitoring campaign consisting of six short (15 ks) observations performed every two days. We discuss the implications of this occultation within the scenario of a Compton-thick cloud crossing the line of sight of the X-ray source. We estimate a source size R<10^14 cm and a distance of the cloud from the source D<10^16 cm. This direct measurement confirms the theoretical expectations of an extremely compact X-ray source, and shows that the Compton-thick circumnuclear gas is located at a distance from the center on the scale of the Broad Line Region.


title: A method to measure the mirror reflectivity of a prime focus telescope

abstract: We have developed a method to measure the mirror reflectivity of telescopes. While it is relatively easy to measure the local reflectivity of the mirror material, it is not so straightforward to measure the amount of light that it focuses in a spot of a given diameter. Our method is based on the use of a CCD camera that is fixed on the mirror dish structure and observes simultaneously part of the telescope's focal plane and the sky region around its optical axis. A white diffuse reflecting disk of known reflectivity is fixed in the telescopes focal plane. During a typical reflectivity measurement the telescope is directed to a selected star. The CCD camera can see two images of the selected star, one directly and another one as a spot focused by the mirror on the white disk. The ratio of the reflected starlight integrated by the CCD from the white disk to the directly measured one provides a precise result of the product of (mirror area x mirror reflectivity).


title: Discovery of hot supergiant stars near the Galactic center

abstract: We report new results of a campaign to find Wolf-Rayet and O (WR/O) stars and high-mass X-ray binaries (HMXBs) in the Galactic center. We searched for candidates by cross-correlating the 2MASS catalog with a deep Chandra catalog of X-ray point sources in the Radio Arches region. Following up with K-band spectroscopy, we found two massive stellar counterparts to CXOGC J174555.3-285126 and CXOGC J174617.0-285131, which we classify as a broad-lined WR star of sub-type WN6b and an O Ia supergiant, respectively. Their X-ray properties are most consistent with those of known colliding-wind binaries in the Galaxy and the Large Magellanic Cloud, although a scenario involving low-rate accretion onto a compact object is also possible. The O Ia star lies 4.4 pc in projection from the Quintuplet cluster, and has a radial velocity consistent with that of the Quintuplet, suggesting that this star might have escaped from the cluster. We also present the discovery of a B2 Ia supergiant, which we identified as a candidate massive star using 8 micron Spitzer maps of the Galactic center in a region near the known massive X-ray-emitting star CXOGC J174516.1-290315. We discuss the origin of these stars in the context of evolving stellar clusters in the Galactic center.


title: Effect of nearby supernova remnants on local Cosmic-Rays

abstract: We study in detail the effect of different particle release times from sources on the cosmic-ray (CR) spectrum below $10^{15}eV$ in the Galaxy. We discuss different possible forms of particle injection such as burst-like injection, continuous injection for a finite time, injection from a stationary source and energy dependent injection. When applied to the nearby known supernova remnants, we find that the observed CR anisotropy data favour the burst-like particle injection model for the CR diffusion coefficient $D(E)\propto E^a$ with $a=0.3-0.6$ in the local region. In the study we have also found that the contribution of the sources G114.3+0.3 and Monogem dominate if the observed anisotropy is a result of the effect of the nearby sources. Further study shows that we should not neglect the contribution of the undetected old sources to the local CR anisotropy.


title: High-resolution simulations of clump-clump collisions using SPH with Particle Splitting

abstract: We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using an SPH code which satisfies the Jeans condition by invoking On-the-Fly Particle Splitting. Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor-Ebert, spheres. Collisions are characterised by M_0 (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius), and M (Mach Number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling, and (iii) adiabatic heating in optically thick protostellar fragments. The efficiency of star formation is found to vary between 10% and 30% in the different collisions studied and it appears to increase with decreasing M_0, and/or decreasing b, and/or increasing M. For b<0.5 collisions produce shock compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high, 1 to 5 x 10^{-5} M_sun yr^{-1}, for the first 1 to 3 x 10^4 yrs. The densities in the filaments, n >~ 5 x 10^5 cm^{-3}, are sufficient that they could be mapped in NH_3 or CS line radiation, in nearby star formation regions.


title: The structure of the nuclear stellar cluster of the Milky Way

abstract: We present high-resolution seeing limited and AO NIR imaging observations of the stellar cluster within about one parsec of Sgr A*, the massive black hole at the centre of the Milky Way. Stellar number counts and the diffuse background light density were extracted from these observations in order to examine the structure of the nuclear stellar cluster.Our findings are as follows: (a) A broken-power law provides an excellent fit to the overall structure of the GC nuclear cluster. The power-law slope of the cusp is $\Gamma=0.19\pm0.05$, the break radius is $R_{\rm break} = 6.0'' \pm 1.0''$ or $0.22\pm0.04$ pc, and the cluster density decreases with a power-law index of $\Gamma=0.75\pm0.1$ outside of $R_{\rm break}$. (b) Using the best velocity dispersion measurements from the literature, we derive higher mass estimates for the central parsec than assumed until now. The inferred density of the cluster at the break radius is $2.8\pm1.3\times 10^{6} {\rm M_{\odot} pc^{-3}}$. This high density agrees well with the small extent and flat slope of the cusp. Possibly, the mass of the stars makes up only about 50% of the total cluster mass. (c) Possible indications of mass segregation in the cusp are found (d) The cluster appears not entirely homogeneous. Several density clumps are detected that are concentrated at projected distances of $R=3''$ and $R=7''$ from Sgr A*.(e) There appears to exist an under-density of horizontal branch/red clump stars near $R=5''$, or an over-density of stars of similar brightness at $R=3''$ and $R=7''$. (f) The extinction map in combination with cometary-like features in an L'-band image may provide support for the assumption of an outflow from Sgr A*.


title: Tensor Microwave Background Fluctuations for Large Multipole Order

abstract: We present approximate formulas for the tensor BB, EE, TT, and TE multipole coefficients for large multipole order l. The error in using the approximate formula for the BB multipole coefficients is less than cosmic variance for l>10. These approximate formulas make various qualitative properties of the calculated multipole coefficients transparent: specifically, they show that, whatever values are chosen for cosmological parameters, the tensor EE multipole coefficients will always be larger than the BB coefficients for all l>15, and that these coefficients will approach each other for l<<100. These approximations also make clear how these multipole coefficients depend on cosmological parameters.


title: Gravitational recoil: signatures on the massive black hole population

abstract: In the last stages of a black hole merger, the binary can experience a recoil due to asymmetric emission of gravitational radiation. Recent numerical relativity simulations suggest that the recoil velocity can be as high as a few thousands kilometers per second for particular configurations. We consider here the effect of a worst case scenario for orbital and phase configurations on the hierarchical evolution of the massive black hole (MBH) population. The orbital configuration and spin orientation in the plane is chosen to be the one yielding the highest possible kick. Masses and spin magnitudes are instead derived self-consistently from the MBH evolutionary models. If seeds form early, e.g. as remnants of the first stars, almost the totality of the first few generation of binaries are ejected. The fraction of lost binaries decreases at later times due to a combination of the binary mass ratio distribution becoming shallower, and the deepening of the hosts potential wells. If seeds form at later times, in more massive halos, the retention rate is much higher. We show that the gravitational recoil does not pose a threat to the evolution of the MBH population that we observe locally in either case, although high mass seeds seem to be favored. The gravitational recoil is instead a real hazard for (i) MBHs in biased halos at high-redshift, where mergers are more common, and the potential wells still relatively shallow. Similarly, it is very challenging to retain (ii) MBHs merging in star clusters.


title: Pulsed X-ray Emission from Pulsar A in the Double Pulsar System J0737-3039

abstract: The double pulsar system J0737-3039 is not only a test bed for General Relativity and theories of gravity, but also provides a unique laboratory for probing the relativistic winds of neutron stars. Recent X-ray observations have revealed a point source at the position of the J0737-3039 system, but have failed to detect pulsations or orbital modulation. Here we report on Chandra X-ray Observatory High Resolution Camera observations of the double pulsar. We detect deeply modulated, double-peaked X-ray pulses at the period of PSR J0737-3039A, similar in appearance to the observed radio pulses. The pulsed fraction is ~70%. Purely non-thermal emission from pulsar A plausibly accounts for our observations. However, the X-ray pulse morphology of A, in combination with previously reported spectral properties of the X-ray emission, allows the existence of both non-thermal magnetospheric emission and a broad sinusoidal thermal emission component from the neutron star surface. No pulsations are detected from pulsar B, and there is no evidence for orbital modulation or extended nebular structure. The absence of orbital modulation is consistent with theoretical expectations of a Poynting-dominated relativistic wind at the termination shock between the magnetosphere of B and the wind from A, and with the small fraction of the energy outflow from A intercepted by the termination shock.


title: PGPG: An Automatic Generator of Pipeline Design for Programmable GRAPE Systems

abstract: We have developed PGPG (Pipeline Generator for Programmable GRAPE), a software which generates the low-level design of the pipeline processor and communication software for FPGA-based computing engines (FBCEs). An FBCE typically consists of one or multiple FPGA (Field-Programmable Gate Array) chips and local memory. Here, the term "Field-Programmable" means that one can rewrite the logic implemented to the chip after the hardware is completed, and therefore a single FBCE can be used for calculation of various functions, for example pipeline processors for gravity, SPH interaction, or image processing. The main problem with FBCEs is that the user need to develop the detailed hardware design for the processor to be implemented to FPGA chips. In addition, she or he has to write the control logic for the processor, communication and data conversion library on the host processor, and application program which uses the developed processor. These require detailed knowledge of hardware design, a hardware description language such as VHDL, the operating system and the application, and amount of human work is huge. A relatively simple design would require 1 person-year or more. The PGPG software generates all necessary design descriptions, except for the application software itself, from a high-level design description of the pipeline processor in the PGPG language. The PGPG language is a simple language, specialized to the description of pipeline processors. Thus, the design of pipeline processor in PGPG language is much easier than the traditional design. For real applications such as the pipeline for gravitational interaction, the pipeline processor generated by PGPG achieved the performance similar to that of hand-written code. In this paper we present a detailed description of PGPG version 1.0.


title: The XENON10 WIMP Search Experiment at the Gran Sasso Underground Laboratory

abstract: XENON10 is a new direct dark matter detection experiment using liquid xenon as target for weakly interacting, massive particles (WIMPs). A two-phase (liquid/gas) time projection chamber with 15 kg fiducial mass has been installed in a low-background shield at the Gran Sasso Underground Laboratory in July 2006. After initial performance tests with various calibration sources, the science data run started on August 24, 2006. The detector has been running stably since then, and a full analysis of more than 75 live days of WIMP search data is now in progress. We present first results on gamma and neutron calibration runs, as well as a preliminary analysis of a subset of the WIMP search data.


title: A bar in the inner halo of barred galaxies I. Structure and kinematics of a representative model

abstract: N-body simulations argue that the inner haloes of barred galaxies should not be spherical, nor even axisymmetric, but triaxial. The departure from sphericity is strongest near the centre and decreases outwards; typical axial ratios for the innermost parts are of the order of 0.8. The halo shape is prolate-like in the inner parts up to a certain radius and then turns to oblate-like. I call this inner halo structure the `halo bar' and analyse here in depth its structure and kinematics in a representative model. It is always considerably shorter than the disc bar. It lags the disc bar by only a few degrees at all radii and the difference between the two bar phases increases with distance from the centre. The two bars turn with roughly the same pattern speed. This means that the halo bar is a slow bar, since its corotation radius is much larger than its length. The bisymmetric component in the halo continues well outside the halo bar in the form of an open spiral, trailing behind the disc bar. The inner parts of the halo display some mean rotation in the same sense as the disc rotation. This is more important for particles nearer to the equatorial plane and decreases with increasing distance from it, but is always much smaller than the disc rotation.


title: Supersoft X-ray Light Curve of RS Ophiuchi (2006)

abstract: One of the candidates for Type Ia supernova progenitors, the recurrent nova RS Ophiuchi underwent the sixth recorded outburst in February 2006, and for the first time a complete light curve of supersoft X-ray has been obtained. It shows the much earlier emergence and longer duration of a supersoft X-ray phase than expected before. These characteristics can be naturally understood when a significant amount of helium layer piles up beneath the hydrogen burning zone during the outburst, suggesting that the white dwarf (WD) is effectively growing up in mass. We have estimated the WD mass in RS Oph to be 1.35 \pm 0.01 M_sun and the growth rate of the WD mass to be at an average rate of about 1 \times 10^{-7} M_sun yr^{-1}. The white dwarf will probably reach the critical mass for Type Ia explosion if the present accretion continues further for a few to several times 10^{5} years.


title: Comparisons for Esta-Task3: Cles and Cesam

abstract: We present the results of comparing three different implementations of the microscopic diffusion process in the stellar evolution codes CESAM and CLES. For each of these implementations we computed models of 1.0, 1.2 and 1.3 M$_{\odot}$. We analyse the differences in their internal structure at three selected evolutionary stages, as well as the variations of helium abundance and depth of the stellar convective envelope. The origin of these differences and their effects on the seismic properties of the models are also considered.


title: Chemical clocks for early-type galaxies

abstract: We present a detailed stellar population analysis of 27 massive elliptical galaxies within 4 very rich clusters at redshift z~0.2. We obtained accurate estimates of the mean luminosity-weighted ages and relative abundances of CN, Mg and Fe as functions of the galaxy velocity dispersion, sigma. Our results are compatible with a scenario in which the stellar populations of massive elliptical galaxies, independently of their environment and mass, had formation timescales shorter than ~1 Gyr. This result implies that massive elliptical galaxies have evolved passively since, at least, as long ago as z~2. For a given galaxy mass the duration of star formation is shorter in those galaxies belonging to more dense environments. Finally, we show that the abundance ratios [CN/Fe] and [Mg/Fe] are the key "chemical clocks" to infer the star formation history timescales in ellipticals. In particular, [Mg/Fe] provides an upper limit for those formation timescales, while [CN/Fe] apperars to be the most suitable parameter to resolve them in elliptical galaxies with sigma<300 km/s.


title: Non-thermal desorption from interstellar dust grains via exothermic surface reactions

abstract: Aims: The gas-phase abundance of methanol in dark quiescent cores in the interstellar medium cannot be explained by gas-phase chemistry. In fact, the only possible synthesis of this species appears to be production on the surfaces of dust grains followed by desorption into the gas. Yet, evaporation is inefficient for heavy molecules such as methanol at the typical temperature of 10 K. It is necessary then to consider non-thermal mechanisms for desorption. But, if such mechanisms are considered for the production of methanol, they must be considered for all surface species. Methods: Our gas-grain network of reactions has been altered by the inclusion of a non-thermal desorption mechanism in which the exothermicity of surface addition reactions is utilized to break the bond between the product species and the surface. Our estimated rate for this process derives from a simple version of classical unimolecular rate theory with a variable parameter only loosely contrained by theoretical work. Results: Our results show that the chemistry of dark clouds is altered slightly at times up to 10^6 yr, mainly by the enhancement in the gas-phase abundances of hydrogen-rich species such as methanol that are formed on grain surfaces. At later times, however, there is a rather strong change. Instead of the continuing accretion of most gas-phase species onto dust particles, a steady-state is reached for both gas-phase and grain-surface species, with significant abundances for the former. Nevertheless, most of the carbon is contained in an undetermined assortment of heavy surface hydrocarbons. Conclusions: The desorption mechanism discussed here will be better constrained by observational data on pre-stellar cores, where a significant accretion of species such as CO has already occurred.


title: Suzaku Observations of the Soft X-ray Background

abstract: We have analyzed a pair of Suzaku XIS1 spectra of the soft X-ray background, obtained by observing towards and to the side of a nearby (d = 230 pc) absorbing filament in the southern Galactic hemisphere. We fit multicomponent spectral models to the spectra in order to separate the foreground emission due to the Local Bubble (LB) from the background emission due to the Galactic halo and unresolved AGN. We obtain LB and halo parameters that are different from those obtained from our analysis of XMM-Newton spectra from these same directions. The LB temperature is lower (log T = 5.93 versus 6.06), and the flux due to the LB in the Suzaku band is an order of magnitude less than is expected from our XMM-Newton analysis. The halo components, meanwhile, are hotter than previously determined, implying our Suzaku spectra are harder than our XMM-Newton spectra.


title: 3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

abstract: Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor $\gamma = 2.5$ jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are $\lesssim c/\sqrt 3$ in the weakly magnetized simulations and $\lesssim 0.3c$ in the strongly magnetized simulations. The Alfven wave speed is $\lesssim 0.07c$ in the weakly magnetized simulations and $\lesssim 0.56c$ in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.


title: On model selection forecasting, Dark Energy and modified gravity

abstract: The Fisher matrix approach (Fisher 1935) allows one to calculate in advance how well a given experiment will be able to estimate model parameters, and has been an invaluable tool in experimental design. In the same spirit, we present here a method to predict how well a given experiment can distinguish between different models, regardless of their parameters. From a Bayesian viewpoint, this involves computation of the Bayesian evidence. In this paper, we generalise the Fisher matrix approach from the context of parameter fitting to that of model testing, and show how the expected evidence can be computed under the same simplifying assumption of a gaussian likelihood as the Fisher matrix approach for parameter estimation. With this `Laplace approximation' all that is needed to compute the expected evidence is the Fisher matrix itself. We illustrate the method with a study of how well upcoming and planned experiments should perform at distinguishing between Dark Energy models and modified gravity theories. In particular we consider the combination of 3D weak lensing, for which planned and proposed wide-field multi-band imaging surveys will provide suitable data, and probes of the expansion history of the Universe, such as proposed supernova and baryonic acoustic oscillations surveys. We find that proposed large-scale weak lensing surveys from space should be able readily to distinguish General Relativity from modified gravity models.


title: GRB 060505: A possible short-duration gamma-ray burst in a star forming region at redshift of 0.09

abstract: On 2006 May 5, a four second duration, low energy, ~10^49 erg, Gamma-Ray Burst (GRB) was observed, spatially associated with a z=0.0894 galaxy. Here, we report the discovery of the GRB optical afterglow and observations of its environment using Gemini-south, Hubble Space Telescope (HST), Chandra, Swift and the Very Large Array. The optical afterglow of this GRB is spatially associated with a prominent star forming region in the Sc-type galaxy 2dFGRS S173Z112. Its proximity to a star forming region suggests that the progenitor delay time, from birth to explosion, is smaller than about 10 Myr. Our HST deep imaging rules out the presence of a supernova brighter than an absolute magnitude of about -11 (or -12.6 in case of ``maximal'' extinction) at about two weeks after the burst, and limits the ejected mass of radioactive Nickel 56 to be less than about 2x10^-4 solar mass (assuming no extinction). Although it was suggested that GRB 060505 may belong to a new class of long-duration GRBs with no supernova, we argue that the simplest interpretation is that the physical mechanism responsible for this burst is the same as for short-duration GRBs.


title: Astrometric-spectroscopic determination of the absolute masses of the HgMn binary star Phi Herculis

abstract: The Mercury-Manganese star Phi Her is a well known spectroscopic binary that has been the subject of a recent study by Zavala et al. (2006), in which they resolved the companion using long-baseline interferometry. The total mass of the binary is now fairly well established, but the combination of the spectroscopy with the astrometry has not resulted in individual masses consistent with the spectral types of the components. The motion of the center of light of Phi Her was clearly detected by the Hipparcos satellite. Here we make use of the Hipparcos intermediate data (`abscissa residuals') and show that by combining them in an optimal fashion with the interferometry the individual masses can be obtained reliably using only astrometry. We re-examine and then incorporate existing radial-velocity measurements into the orbital solution, obtaining improved masses of 3.05 +/- 0.24 M_Sun and 1.614 +/- 0.066 M_Sun that are consistent with the theoretical mass-luminosity relation from recent stellar evolution models. These mass determinations provide important information for the understanding of the nature of this peculiar class of stars.


title: Systematic errors and combination of individual CRF solutions in the framework of the international pilot project for the next ICRF

abstract: A new international Pilot Project for the re-determination of the ICRF was initiated by the International VLBI Service for Geodesy and Astrometry (IVS) in January 2005. The purpose of this project is to compare the individual CRF solutions and to analyze their systematic and random errors with focus on the selection of the optimal strategy for the next ICRF realization. Eight radio source catalogues provided by the IVS Analysis Centers GA, SHAO, DGFI, GIUB-BKG, JPL, MAO NANU, GSFC, USNO were analyzed. In present study, four analytical models were used to investigate the systematic differences between solutions: solid rotation, rotation and deformation (IERS method), and expansion in orthogonal functions: Legendre-Fourier polynomials and spherical functions. It was found that expansions by orthogonal function describe the differences between individual catalogues better than the two former models. Finally, the combined CRF was generated. Using the radio source positions from this combined catalogue for estimation of EOP has shown improvement of the uncertainty of the celestial pole offset time series.


title: Do Mergers Spin up Dark Matter Halos?

abstract: We use a large cosmological N-body simulation to study the origin of possible correlations between the merging history and spin of cold dark matter halos. In particular, we examine claims that remnants of major mergers tend to have higher-than-average spins, and find that the effect is driven largely by unrelaxed systems: equilibrium dark matter halos show no significant correlation between spin and merger history. Out-of-equilibrium halos have, on average, higher spins than relaxed systems, suggesting that the virialization process leads to a net decrease in the value of the spin parameter. We find that this decrease is driven by the internal redistribution of mass and angular momentum that occurs during virialization, a process that is especially efficient during major mergers, when high angular momentum material is pushed beyond the virial radius of the remnant. Since such redistribution likely affects the angular momentum of baryons and dark matter unevenly, our findings question the common practice of identifying the specific angular momentum content of a halo with that of its embedded luminous component. Further work is needed to elucidate the true relation between the angular momentum content of baryons and dark matter in galaxy systems assembled hierarchically.