title: Eclipsing Binary Stars in Open Clusters

abstract: The study of detached eclipsing binaries in open clusters can provide stringent tests of theoretical stellar evolutionary models, which must simultaneously fit the masses, radii, and luminosities of the eclipsing stars and the radiative properties of every other star in the cluster. We review recent progress in such studies and discuss two unusually interesting objects currently under analysis. GV Carinae is an A0m + A8m binary in the Southern open cluster NGC 3532; its eclipse depths have changed by 0.1 mag between 1990 and 2001, suggesting that its orbit is being perturbed by a relatively close third body. DW Carinae is a high-mass unevolved B1V + B1V binary in the very young open cluster Collinder 228, and displays double-peaked emission in the centre of the Halpha line which is characteristic of Be stars. We conclude by pointing out that the great promise of eclipsing binaries in open clusters can only be satisfied when both the binaries and their parent clusters are well-observed, a situation which is less common than we would like.


title: Revised WMAP constraints on neutrino masses and other extensions of the minimal $\Lambda$CDM model

abstract: Recently, two issues concerning the three-year WMAP likelihood code were pointed out. On large angular scales ($l \lesssim 30$), a sub-optimal likelihood approximation resulted in a small power excess. On small angular scales ($l \gtrsim 300$), over-subtraction of unresolved point sources produced a small power deficit. For a minimal six-parameter cosmological model, these two effects conspired to decrease the value of $n_s$ by $\sim 0.7 \sigma$. In this paper, we study the change in preferred parameter ranges for more extensive cosmological models, including running of $n_s$, massive neutrinos, curvature, and the equation of state for dark energy. We also include large-scale structure and supernova data in our analysis. We find that the parameter ranges for $\alpha_s$, $\Omega_k$ and $w$ are not much altered by the modified analysis. For massive neutrinos the upper limit on the sum of the neutrino masses decreases from $M_\nu < 1.90$eV to $M_\nu < 1.57$eV when using the modified WMAP code and WMAP data only. We also find that the shift of $n_s$ to higher values is quite robust to these extensions of the minimal cosmological model.


title: Accurate stellar masses in the multiple system T Tau

abstract: The goal of this study is to obtain accurate estimates for the individual masses of the components of the tight binary system T Tau S in order to settle the ongoing debate on the nature of T Tau Sa, a so-called infrared companion. We take advantage of the fact that T Tau S belongs to a triple system composed of two hierarchical orbits to simultaneously analyze the motion of T Tau Sb in the rest frames of T Tau Sa and T Tau N. With this method, it is possible to pinpoint the location of the center of mass of T Tau S and, thereby, to determine individual masses for T Tau Sa and T Tau Sb with no prior assumption about the mass/flux ratio of the system. This improvement over previous studies of the system results in much better constraints on orbital parameters. We find individual masses of 2.73+/-0.31 Msun for T Tau Sa and of 0.61+/-0.17 Msun for T Tau Sb (in agreement with its early-M spectral type), including the uncertainty on the distance to the system. These are among the most precise estimates of the mass of any Pre-Main Sequence star, a remarkable result since this is the first system in which individual masses of T Tauri stars can be determined from astrometry only. This model-independent analysis confirms that T Tau Sa is an intermediate-mass star, presumably a very young Herbig Ae star, that may possess an almost edge-on disk.


title: Halo and Galaxy Formation Histories from the Millennium Simulation: Public release of a VO-oriented and SQL-queryable database for studying the evolution of galaxies in the LambdaCDM cosmogony

abstract: The Millennium Run is the largest simulation of the formation of structure within the $\Lambda$CDM cosmogony so far carried out. It uses $10^{10}$ particles to follow the dark matter distribution in a cubic region 500$h^{-1}$Mpc on a side, and has a spatial resolution of 5 $h^{-1}$kpc. Application of simplified modelling techniques to the stored output of this calculation allows the formation and evolution of the $\sim 10^7$ galaxies more luminous than the Small Magellanic Cloud to be simulated for a variety of assumptions about the detailed physics involved. As part of the activities of the German Astrophysical Virtual Observatory we have used a relational database to store the detailed assembly histories both of all the haloes and subhaloes resolved by the simulation, and of all the galaxies that form within these structures for two independent models of the galaxy formation physics. We have created web applications that allow users to query these databases remotely using the standard Structured Query Language (SQL). This allows easy access to all properties of the galaxies and halos, as well as to the spatial and temporal relations between them and their environment. Information is output in table format compatible with standard Virtual Observatory tools and protocols. With this announcement we are making these structures fully accessible to all users. Interested scientists can learn SQL, gain familiarity with the database design and test queries on a small, openly accessible version of the Millennium Run (with volume 1/512 that of the full simulation). They can then request accounts to run similar queries on the databases for the full simulations.


title: The formation of an eccentric gap in a gas disk by a planet in an eccentric orbit

abstract: We investigate the effect of a planet on an eccentric orbit on a two dimensional low mass gaseous disk. At a planet eccentricity above the planet's Hill radius divided by its semi-major axis, we find that the disk morphology differs from that exhibited by a disk containing a planet in a circular orbit. An eccentric gap is created with eccentricity that can exceed the planet's eccentricity and precesses with respect to the planet's orbit. We find that a more massive planet is required to open a gap when the planet is on an eccentric orbit. We attribute this behavior to spiral density waves excited at corotation resonances that increase the disk's eccentricity and exert a torque opposite in sign to that exerted by the Lindblad resonances. The reduced torque makes it more difficult for waves driven by the planet to overcome viscous inflow in the disk. Spectral energy distributions of transitional disks may reveal the presence of an eccentric planet if they are matched by an inner edge comprising a range of temperatures set by the semi-major axis and eccentricity of the hole.


title: Hybrid protoneutron stars with the MIT bag model

abstract: We study the hadron-quark phase transition in the interior of protoneutron stars. For the hadronic sector, we use a microscopic equation of state involving nucleons and hyperons derived within the finite-temperature Brueckner-Bethe-Goldstone many-body theory, with realistic two-body and three-body forces. For the description of quark matter, we employ the MIT bag model both with a constant and a density-dependent bag parameter. We calculate the structure of protostars with the equation of state comprising both phases and find maximum masses below 1.6 solar masses. Metastable heavy hybrid protostars are not found.


title: Methods for improving open cluster fundamental parameters applied to M 52 and NGC 3960

abstract: We derive accurate parameters related to the CMD, structure and dynamical state of M52 and NGC3960, whose fields are affected by differential reddening. Compared to open clusters in different dynamical states studied with similar methods, the core and overall parameters of M52 are consistent with an open cluster more massive than 1000Mo and ~60Myr old, with some mass segregation in the inner region. The core of NGC3960 is in an advanced dynamical state with strong mass segregation in the core/halo region, while the somewhat flat overall MF (x~1.07) suggests low-mass star evaporation. The excess stellar density in the core may suggest post-core collapse. The dynamical evolution of NGC3960 may have been accelerated by the tidal Galactic field, since it lies \~0.5kpc inside the Solar circle.


title: The bright spiky emission of pulsar B0656+14

abstract: We present a detailed study of the single radio pulses of PSR B0656+14, a pulsar also known to be a strong pulsed source of high-energy emission. Despite the extensive studies at high-energy wavelengths, there is little or no published work on its single-pulse behaviour in the radio band. In this report we rectify this omission. We have found that the shape of the pulse profile of PSR B0656+14 requires an unusually long timescale to achieve stability (over 25,000 pulses at 327 MHz). This instability is caused by very bright and narrow pulses with widths and luminosities comparable to those observed for the RRATs. Many pulses are bright enough to qualify as "giant pulses", but are broader than those usually meant by this term. At 327 MHz the brightest pulse was about 116 times brighter than the average pulse. Although the most powerful pulses peak near the centre of the profile, occasional sudden strong pulses are also found on the extreme leading edge of the profile. One of them has a peak flux of about 2000 times the average flux at that pulse longitude. No "break" in the pulse-energy distributions is observed, but nevertheless there is evidence of two separate populations of pulses: bright pulses have a narrow "spiky" appearance consisting of short quasi-periodic bursts of emission with microstructure, in contrast to the underlying weaker broad pulses. Furthermore, the spiky pulses tend to appear in clusters which arise and dissipate over about 10 periods. We demonstrate that the spiky emission builds a narrow and peaked profile, whereas the weak emission produces a broad hump, which is largely responsible for the shoulders in the total emission profiles at both high and low frequencies.


title: A Comparison of Ultraluminous X-ray Sources in NGC 1399 and the Antennae Galaxies (NGC 4038/4039)

abstract: The temporal and spectral properties of ultraluminous X-ray sources (ULXs, L>2x10^39 ergs/s) and bright X-ray sources (L>3x10^38 ergs/s) are examined and compared in two extremely different host environments: the old elliptical galaxy NGC 1399 and the young, starforming Antennae galaxies (NGC 4038/4039). ULXs in NGC 1399 show little variability on either long or short time scales. Only 1 of 8 ULXs and 10 of 63 bright sources in NGC 1399 are variable at a confidence level of 90%. On long timescales, the NGC 1399 sources are steadier than most Galactic black hole X-ray binaries, but similar to GRS 1915+105. The outburst duration of the NGC 1399 sources is about 20 yrs, again, similar to that of GRS 1915+105. The bright X-ray sources in NGC 1399 may be black hole X-ray binaries with giant star companions similar to GRS 1915+105. The brightest ULX (PSX-1) in NGC 1399 is coincident with a globular cluster, shows a hard spectrum with a photon index around 1.5, and has a nearly constant luminosity around 5x10^39 erg/s. It may be an intermediate-mass black hole (IMBH) in a hard spectral state. In contrast to NGC 1399, the ULXs in the Antennae are all variable and a large fraction of the bright sources (9 of 15) are also variable. The variability and luminosity of ULXs in the Antennae suggest they are black hole high mass X-ray binaries accreting via Roche-lobe overflow. A flare with a duration of about 5 ks is found from Antennae X-42. The most luminous ULX, X-16, with a very hard spectrum (Gamma=1.0~1.3) and a luminosity which varies by a factor of 10, could be an IMBH candidate.


title: Studies of Neutron Stars at Optical/IR Wavelengths

abstract: In the last years, optical studies of Isolated Neutron Stars (INSs) have expanded from the more classical rotation-powered ones to other categories, like the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma-ray Repeaters (SGRs), which make up the class of the magnetars, the radio-quiet INSs with X-ray thermal emission and, more recently, the enigmatic Compact Central Objects (CCOs) in supernova remnants. Apart from 10 rotation-powered pulsars, so far optical/IR counterparts have been found for 5 magnetars and for 4 INSs. In this work we present some of the latest observational results obtained from optical/IR observations of different types of INSs.


title: The Structure of Jupiter, Saturn, and Exoplanets: Key Questions for High-Pressure Experiments

abstract: We give an overview of our current understanding of the structure of gas giant planets, from Jupiter and Saturn to extrasolar giant planets. We focus on addressing what high-pressure laboratory experiments on hydrogen and helium can help to elucidate about the structure of these planets.


title: Galaxy Luminosities in 2PIGG groups

abstract: We derive composite luminosity functions (LF) for galaxies in groups and examine the behaviour of the LF as a function of group luminosity (used as an indicator of group or halo mass). We consider both the entire galaxy population and split galaxies into red and blue (quiescent and star forming) samples, in order to examine possible mechanisms behind observed variations of galaxy properties with environment. We find evidence that $M^*$ brightens and $\alpha$ steepens with group luminosity, until a threshold value where the LF parameters stabilize at those found in rich clusters. The effect is seen in the total LF and for the blue and red galaxies separately. The behaviour of the quiescent and star-forming samples is qualitatively consistent with variations resulting from interactions and mergers, where mergers build the bright end of the luminosity function at the same time as dwarf irregulars have their star formation quenched and evolve into dwarf ellipticals. These processes appear to take place preferentially in low luminosity groups and to be complete at a group luminosity of -22.5 in $B$, corresponding to a halo mass of order $10^{13.5} M_{\odot}$.


title: Collapsars in Three Dimensions

abstract: We present the results of 3-dimensional simulations of the direct collapse to a black hole of a rotating, 60-solar-mass, zero metallicity, population III star. Because the structure of this star (angular momentum, density and temperature profiles) is similar to many collapsar gamma-ray burst progenitors, these calculations have implications beyond the fate of population III stars. These simulations provide a first 3-dimensional look at a realistic collapsar progenitor, and the results are very different from any previous 2-dimensional calculations. If the angular momentum of the progenitor is high, non-axisymmetric instabilities in the collapsing core cause spiral structures to form, and these structures shape later outflows. These outflows are driven by the imbalance between viscous heating and inefficient neutrino cooling and ultimately develop into a 1e52 erg explosion. Without magnetic fields, this collapse will not produce relativistic jets, but the explosion is indeed a hypernova. We conclude with a discussion of the implications of such calculations on the explosions, nucleosynthesis, neutrino flux and gravitational wave emission from the collapse of massive stars.


title: On the progenitor of SN 2005gl and the nature of Type IIn supernovae

abstract: We present a study of the type IIn supernova (SN) 2005gl, in the relatively nearby (d~66 Mpc) galaxy NGC 266. Photometry and spectroscopy of the SN indicate it is a typical member of its class. Pre-explosion Hubble Space Telescope (HST) imaging of the location of the SN, along with a precise localization of this event using the Laser-Guide-Star assisted Adaptive Optics (LGS-AO) system at Keck Observatory, are combined to identify a luminous (M_V=-10.3) point source as the possible progenitor of SN 2005gl. If the source is indeed a single star, it was likely a member of the class of luminous blue variable stars (LBVs). This finding leads us to consider the possible general association of SNe IIn with LBV progenitors. We find this is indeed supported by observations of other SNe, and the known properties of LBV stars. For example, we argue that should the prototypical Galactic LBV eta Carina explode in a phase similar to its current state, it will likely produce a type IIn SN. We discuss our findings in the context of current ideas about the evolution of massive stars, and review the census of SNe with identified progenitors. We introduce the concept of the progenitor-SN map as a convenient means to discuss the present status and future prospects of direct searches for SN progenitors. We conclude that this field has matured considerably in recent years, and the transition from anecdotal information about rare single events to robust associations of progenitor classes with specific SN types has already begun.


title: Measuring M dwarf Winds with DAZ White Dwarfs

abstract: Hydrogen atmosphere white dwarfs with metal lines, so-called DAZs, show evidence for ongoing accretion of material onto their surfaces. Some DAZs are known to have unresolved M dwarf companions, which could account for the observed accretion through a stellar wind. I combine observed Ca abundances of the DAZs with information on the orbital separation of their M dwarf companions to infer the mass loss rate of the M dwarfs. I find that for three of the six known DAZs with M dwarf companions, a stellar wind can plausibly explain the observed accretion on the white dwarfs assuming Bondi-Hoyle accretion of solar abundance stellar winds on the order of 10$^{-14}-10^{-16}\Msun$ yr$^{-1}$. The rest of the sample have companions with orbits $\gtorder$ 1~AU, and require companion mass loss rates of $> 10^{-11}\Msun$ yr$^{-1}$. I conclude that there must be an alternative explanation for accretion of material onto DAZs with widely separated companions. The inferred winds for two of the close binaries are orders of magnitude smaller than typically assumed for the angular momentum loss of red dwarf-white dwarf pairs due to magnetic braking from a stellar wind and may seriously affect predictions for the formation rate of CVs with low mass companions.


title: Energy Conditions and Supernovae Observations

abstract: In general relativity, the energy conditions are invoked to restrict general energy-momentum tensors on physical grounds. We show that in the standard Friedmann-Lemaitre-Robertson-Walker approach to cosmological modelling where the equation of state of the cosmological fluid is unknown, the energy conditions provide model-independent bounds on the behavior of the distance modulus of cosmic sources as a function of the redshift. We use both the gold and the legacy samples of current type Ia supenovae to carry out a model-independent analysis of the energy conditions violation in the context of standard cosmology.


title: Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe

abstract: Observations of the high-redshift Universe with the 21 cm hyperfine line of neutral hydrogen promise to open an entirely new window onto the early phases of cosmic structure formation. Here we review the physics of the 21 cm transition, focusing on processes relevant at high redshifts, and describe the insights to be gained from such observations. These include measuring the matter power spectrum at z~50, observing the formation of the cosmic web and the first luminous sources, and mapping the reionization of the intergalactic medium. The epoch of reionization is of particular interest, because large HII regions will seed substantial fluctuations in the 21 cm background. We also discuss the experimental challenges involved in detecting this signal, with an emphasis on the Galactic and extragalactic foregrounds. These increase rapidly toward low frequencies and are especially severe for the highest redshift applications. Assuming that these difficulties can be overcome, the redshifted 21 cm line will offer unique insight into the high-redshift Universe, complementing other probes but providing the only direct, three-dimensional view of structure formation from z~200 to z~6.


title: Nonthermal gamma-ray and X-ray flashes from shock breakout in gamma-ray bursts/supernovae

abstract: Thermal X-ray emission which is simultaneous with the prompt gamma-rays has been detected for the first time from a supernova connected with a gamma-ray burst (GRB), namely GRB060218/SN2006aj. It has been interpreted as arising from the breakout of a mildly relativistic, radiation-dominated shock from a dense stellar wind surrounding the progenitor star. There is also evidence for the presence of a mildly relativistic ejecta in GRB980425/SN1998bw, based on its X-ray and radio afterglow. Here we study the process of repeated bulk Compton scatterings of shock breakout thermal photons by the mildly relativistic ejecta. During the shock breakout process, a fraction of the thermal photons would be repeatedly scattered between the pre-shock material and the shocked material as well as the mildly relativistic ejecta and, as a result, the thermal photons get boosted to increasingly higher energies. This bulk motion Comptonization mechanism will produce nonthermal gamma-ray and X-ray flashes, which could account for the prompt gamma-ray burst emission in low-luminosity supernova-connected GRBs, such as GRB060218. A Monte Carlo code has been developed to simulate this repeated scattering process, which confirms that a significant fraction of the thermal photons get "accelerated" to form a nonthermal component, with a dominant luminosity. This interpretation for the prompt nonthermal emission of GRB060218 may imply that either the usual internal shock emission from highly relativistic jets in these low-luminosity GRBs is weak, or alternatively, that there are no highly relativistic jets in this peculiar class of bursts.


title: The accelerating universe and a limiting curvature proposal

abstract: We consider the hypothesis of a limiting minimal curvature in gravity as a way to construct a class of theories exhibiting late-time cosmic acceleration. Guided by the minimal curvature conjecture (MCC) we are naturally lead to a set of scalar tensor theories in which the scalar is non-minimally coupled both to gravity and to the matter Lagrangian. The model is compared to the Lambda Cold Dark Matter concordance model and to the observational data using the gold SNeIa sample of Riess et. al. (2004). An excellent fit to the data is achieved. We present a toy model designed to demonstrate that such a new, possibly fundamental, principle may be responsible for the recent period of cosmological acceleration. Observational constraints remain to be imposed on these models.


title: An Eccentric Hot Jupiter Orbiting the Subgiant HD 185269

abstract: We report the detection of a Jupiter-mass planet in a 6.838 day orbit around the 1.28 solar mass subgiant HD 185269. The eccentricity of HD 185269b (e = 0.30) is unusually large compared to other planets within 0.1 AU of their stars. Photometric observations demonstrate that the star is constant to +/-0.0001 mag on the radial velocity period, strengthening our interpretation of a planetary companion. This planet was detected as part of our radial velocity survey of evolved stars located on the subgiant branch of the H-R diagram--also known as the Hertzsprung Gap. These stars, which have masses between 1.2 and 2.5 solar masses, play an important role in the investigation of the frequency of extrasolar planets as a function of stellar mass.


title: First Constraints on Rings in the Pluto System

abstract: Simple theoretical calculations have suggested that small body impacts onto Pluto's newly discovered small satellites, Nix and Hydra, are capable of generating time-variable rings or dust sheets in the Pluto system. Using HST/ACS data obtained on 2006 February 15 and 2006 March 2, we find no observational evidence for such a ring system and present the first constraints on the present-day I/F and optical depth of a putative ring system. At the 1500-km radial resolution of our search, we place a 3-sigma upper limit on the azimuthally-averaged normal I/F of ring particles of 5.1x10^-7 at a distance of 42,000 km from the Pluto-Charon barycenter, the minimum distance for a dynamically stable ring (Stern et al., 1994; Nagy et al., 2006); 4.4x10^-7 at the orbit of Nix; and 2.5x10^-7 at the orbit of Hydra. For an assumed ring particle albedo of 0.04 (0.38), these I/F limits translate into 3-sigma upper limits on the normal optical depth of macroscopic ring particles of 1.3x10^-5 (1.4x10^-6), 1.1x10^-5 (1.2x10^-6), 6.4x10^-6 (6.7x10^-7), respectively. Were the New Horizons spacecraft to fly through a ring system with optical depth of 1.3x10^-5, it would collide with a significant number of potentially damaging ring particles. We therefore recommend that unless tighter constraints can be obtained, New Horizons cross the putative ring plane within 42,000 km of the Pluto-Charon barycenter, where rings are dynamically unstable. We derive a crude estimate of the lifetime of putative ring paritcles of 900 years.


title: X-Ray Scattering by Nonspherical Grains. I. Oblate Spheroids

abstract: We calculate the scattering of X-rays by interstellar dust, for a dust model that reproduces the observed wavelength-dependent extinction and polarization of starlight. On interstellar sightlines that produce appreciable starlight polarization, we predict that the dust-scattered X-ray halo around point sources will have measurable azimuthal asymmetry due to scattering by partially-aligned nonspherical grains. We calculate the expected halo asymmetry. X-ray halo asymmetry provides a new test of interstellar dust models.


title: A Survey and Analysis of Spitzer Infrared Spectrograph Spectra of T Tauri Stars in Taurus

abstract: We present mid-infrared spectra of T Tauri stars in the Taurus star-forming region obtained with the Spitzer Infrared Spectrograph (IRS). For the first time, the 5-36 micron spectra of a large sample of T Tauri stars belonging to the same star-forming region is studied, revealing details of the mid-infrared excess due to dust in circumstellar disks. We analyze common features and differences in the mid-IR spectra based on disk structure, dust grain properties, and the presence of companions. Our analysis encompasses spectral energy distributions from the optical to the far-infrared, a morphological sequence based on the IRS spectra, and spectral indices in IRS wave bands representative of continuum emission. By comparing the observed spectra to a grid of accretion disk models, we infer some basic disk properties for our sample of T Tauri stars, and find additional evidence for dust settling.


title: Characterizing Bars at z~0 in the optical and NIR: Implications for the Evolution of Barred Disks with Redshift

abstract: Critical insights on galaxy evolution stem from the study of bars. With the advent of HST surveys that trace bars in the rest-frame optical out to z~1, it is critical to provide a reference baseline for bars at z~0 in the optical band. We present results on bars at z~0 in the optical and NIR bands based on 180 spirals from OSUBSGS. (1) The deprojected bar fraction at z~0 is ~60% +/-6% in the NIR H-band and ~44% +/-6% in the optical B-band. (2) The results before and after deprojection are similar, which is encouraging for high-redshift studies that forego deprojection. (3) Studies of bars at z~0.2-1.0 (lookback time of 3-8 Gyr) have reported an optical bar fraction of ~30% +/-6%, after applying cutoffs in absolute magnitude (M_V < -19.3), bar size (a_bar >= 1.5 kpc), and bar ellipticity (e_bar >= 0.4). Applying these exact cutoffs to the OSUBSGS data yields a comparable optical B-band bar fraction at z~0 of ~ 34%+/-6%. This rules out scenarios where the optical bar fraction in bright disks declines strongly with redshift. (4) Most (~70%) bars have moderate to high strentgh or ellipticity (0.50 <= e_bar <= 0.75). There is no bimodality in the distribution of e_bar. The H-band bar fraction and e_bar show no substantial variation across RC3 Hubble types Sa to Scd. (5) RC3 bar types should be used with caution. Many galaxies with RC3 types "AB" turn out to be unbarred and RC3 bar classes "B" and "AB" have a significant overlap in e_bar. (6) Most bars have sizes below 5 kpc. Bar and disk sizes correlate, and most bars have a_bar/R_25~0.1-0.5. This suggests that the growths of bars and disks are intimately tied.


title: Searching for Low Surface-brightness Galaxies in the Hubble Ultra Deep Field: Implications for the Star Formation Efficiency in Neutral Gas at z ~ 3

abstract: (Abridged) The Kennicutt--Schmidt law relates the face-on star formation rate (SFR) per unit area with the face-on gaseous column density in nearby galaxies. Applying this relation to damped Lya absorption systems (DLAs) of neutral hydrogen column density N > 1.6x10^21 cm^-2 leads to an estimate that three percent of the sky should be covered with extended sources brighter than mu_V~28.4 mag arcsec^-2, if DLAs at redshift z=[2.5,3.5] undergo in situ star formation. We test this hypothesis by searching the Hubble Ultra Deep Field (UDF) F606W image for low surface-brightness features of angular sizes, ranging between theta_dla=0.25" and 4.0". Our search yields upper limits on the comoving SFR densities that are between factors of 30 and 100 lower than predictions, suggesting a reduction by more than a factor of 10 in star formation efficiency at z~3. We consider several mechanisms that could reduce star formation efficiency at high redshift. We find that the cosmological increase with redshift of the critical surface density for the Toomre instability may be sufficient to suppress star formation to the levels implied by the UDF observations. However, the uncertainties are such that Toomre instabilities may still exist. In that case star formation at column densities less than 10^22 cm^-2 may be suppressed by the low molecular content of the DLA gas. The upper limits on in situ star formation reduce the predicted metallicities at z~3 to be significantly lower than observed, and reduce the heat input in the gas to be substantially lower than the inferred cooling rates. In contrast, the radiative output from compact Lyman Break Galaxies (LBGs) with R<27 is sufficient to balance the comoving cooling rate.


title: Understanding the saturation of proton-driven Weibel instabilities and implications for astrophysics

abstract: The linear growth rate and saturation level of magnetic fields for Weibel instabilities driven by ion temperature anisotropy, defined as $\alpha=(T_\perp/T_\parallel)-1$ where $T_\perp$ and $T_\parallel$ are ion temperatures perpendicular and parallel to the wave vector, are derived in the small $\alpha$-limit. It is shown that the ratio of the saturated magnetic energy to the initial ion energy scales as the fourth power of the electron to ion mass ratio, $m/M$, for an initially unmagnetized plasma with $\alpha\leq M/m$. Particle-in-cell simulations confirm the mass scaling and also show that the electron energy gain is of the same order of magnitude as the magnetic field energy. This implies that the Weibel instabilities cannot provide a faster-than-Coulomb collisionless mechanism to equilibrate ion-electron plasmas with ions initially much hotter than electrons, a key component in low-luminosity astrophysical accretion flows. The results here also show that the large $\alpha$-limit formulae used in the study of magnetic field generation in collisionless shocks are only valid if $\alpha\geq M/m$.


title: TeV Astrophysics, A Review

abstract: In this paper I briefly review recent progress in the field of ground-based gamma ray astrophysics.


title: New Self-Similar Solutions of Polytropic Gas Dynamics

abstract: We explore semi-complete self-similar solutions for the polytropic gas dynamics involving self-gravity under spherical symmetry, examine behaviours of the sonic critical curve, and present new asymptotic collapse solutions that describe `quasi-static' asymptotic behaviours at small radii and large times. These new `quasi-static' solutions with divergent mass density approaching the core can have self-similar oscillations. Earlier known solutions are summarized. Various semi-complete self-similar solutions involving such novel asymptotic solutions are constructed, either with or without a shock. In contexts of stellar core collapse and supernova explosion, a hydrodynamic model of a rebound shock initiated around the stellar degenerate core of a massive progenitor star is presented. With this dynamic model framework, we attempt to relate progenitor stars and the corresponding remnant compact stars: neutron stars, black holes, and white dwarfs.


title: Visible spectroscopy of 2003 UB313: Evidence for N2 ice on the surface of the largest TNO?

abstract: The recent discovery of two large trans-Neptunian objects (TNOs) 2003 UB313 and 2005 FY9, with surface properties similar to those of Pluto, provides an exciting new laboratory for the study of processes considered for Pluto and Triton: volatile mixing and transport; atmospheric freeze-out and escape, ice chemistry, and nitrogen phase transitions. We studied the surface composition of TNO 2003 UB313, the first known TNO larger than Pluto. We report a visible spectrum covering the 0.35-0.95 microns spectral range, obtained with the 4.2m William Herschel Telescope at "El Roque de los Muchachos" Observatory (La Palma, Spain). The visible spectrum of this TNO presents very prominent absorpt- ions bands formed in solid CH4. At wavelengths shorter than 0.6 micron the spectrum is almost featureless and slightly red (S'=4%). The icy-CH4 bands are significantly stronger than those of Pluto and slightly weaker than those observed in the spectrum of another giant TNO, 2005 FY9, implying that methane is more abundant on its surface than in Pluto's and close to that of the surface of 2005 FY9. A shift of 15 +/-3 A relative to the position of the bands of the spectrum of laboratory CH4 ice is observed in the bands at larger wavelengths (e.g. around 0.89 microns), but not at shorter wavelengths (the band around 0.73 microns is not shifted) this may be evidence for a vertical compositional gradient. Purer methane could have condensed first while 2003 UB313 moved towards aphelion during the last 200 years, and as the atmosphere gradually collapsed, the composition became more nitrogen-rich as the last, most volatile components condensed, and CH4 diluted in N2 is present in the outer surface layers.


title: Enhancement and suppression of heat transfer by MHD turbulence

abstract: We study of the effect of turbulence on heat transfer within magnetized plasmas for energy injection velocities both larger and smaller that the Alfven speed. We find that in the latter regime the heat transfer is partially suppressed, while in the former regime the effects of turbulence depend on the intensity of driving. In fact, the scale l at which the turbulent velocity is equal the Alfven velocity is a new important parameter. When the electron mean free path is larger than l, the stronger the the turbulence, the lower thermal conductivity by electrons is. The turbulent motions, however, induces their own advective heat transport, which, for the parameters of intracluster medium (ICM) provides effective heat diffusivity that exceeds the classical Spitzer value.


title: SINS of Viscosity Damped Turbulence

abstract: The problems with explaining the Small Ionized and Neutral Structures (SINS) appealing to turbulence stem from inefficiency of the Kolmogorov cascade in creating large fluctuations at sufficiently small scales. However, other types of cascades are possible. When magnetic turbulence in a fluid with viscosity that is much larger than resistivity gets to a viscous damping scale, the turbulence does not vanish. Instead, it gets into a different new regime. Viscosity-damped turbulence produces fluctuations on the small scales. Magnetic fields sheared by turbulent motions by eddies not damped by turbulence create small scale filaments that are confined by the external plasma pressure. This creates small scale density fluctuations. In addition, extended current sheets create even stronger density gradients that accompany field reversals in the plane perpendicular to mean magnetic field. Those can be responsible for the SINS formation. This scenario is applicable to partially ionized gas. More studies of reconnection in the viscosity dominated regime are necessary to understand better the extend to which the magnetic reversals can compress the gas.


title: Theoretical approaches to particle propagation and acceleration in turbulent intergalactic medium

abstract: Intercluster medium is expected to be turbulent with turbulence being superAlfvenic at large scales. Magnetic fields substantially modify the turbulent cascade when the turbulence reaches the scales at which the fluctuation velocity gets less than the Alfven velocity. At those scales it is possible to consider three cascades, of fast, slow and Alfven modes with little energy exchange between them. As Alfvenic and slow modes are anisotropic they marginally scatter and accelerate cosmic rays, while fast modes dominate the processes. However, in the presence of cosmic rays the turbulence is modified as cosmic rays transfer the energy of compressible motions (i.e. slow and fast modes) from large scales to the scale of cosmic ray Larmor radius. This results in generation of a new small-scale Alfvenic component which is not a part of the ordinary MHD cascade. This component does scatter and accelerate cosmic rays. In addition, magnetic reconnection in turbulent medium accelerates cosmic rays. The complexity of the intracluster turbulence calls for observational studies. A new technique Velocity Coordinate Spectrum (VCS) is particularly promising for studies of velocity fluctuations with a new generation of X-ray observatories.


title: Intermittency of Magnetohydrodynamic Turbulence: Astrophysical Perspective

abstract: Intermittency is an essential property of astrophysical fluids, which demonstrate an extended inertial range. As intermittency violates self-similarity of motions, it gets impossible to naively extrapolate the properties of fluid obtained computationally with relatively low resolution to the actual astrophysical situations. In terms of Astrophysics, intermittency affects turbulent heating, momentum transfer, interaction with cosmic rays, radio waves and many more essential processes. Because of its significance, studies of intermittency call for coordinated efforts from both theorists and observers. In terms of theoretical understanding we are still just scratching a surface of a very rich subject. We have some theoretically well justified models that are poorly supported by experiments, we also have She-Leveque model, which could be vulnerable on theoretical grounds, but, nevertheless, is well supported by experimental and laboratory data. I briefly discuss a rather mysterious property of turbulence called ``extended self-similarity'' and the possibilities that it opens for the intermittency research. Then I analyze simulations of MHD intermittency performed by different groups and show that their results do not contradict to each other. Finally, I discuss the intermittency of density, intermittency of turbulence in the viscosity-dominated regime as well as the intermittency of polarization of Alfvenic modes. The latter provides an attractive solution to account for a slower cascading rate that is observed in some of the numerical experiments. I conclude by claiming that a substantial progress in the field may be achieved by studies of the turbulence intermittency via observations.


title: Grain Alignment, Polarization and Magnetic Fields

abstract: Aligned non-spherical dust particles polarize starlight passing through the dust cloud. They also emit polarized far infrared and sub-mm radiation. Substantial progress in understanding of grain alignment theory makes the interpretation of the polarized radiation in terms of underlying magnetic fields much more reliable. I discuss a number of fundamental processes that affect grain alignment. In particular, I shall discuss how subtle effects related to nuclear spins of the atoms alter the dynamics of dust grains. I shall discuss how the theory explains the existing observational data and demonstrate when the polarization can and cannot be interpreted in terms of the underlying magnetic fields.


title: Coherent Development of Neutrino Flavor in the Supernova Environment

abstract: We calculate coherent neutrino and antineutrino flavor transformation in the supernova environment, for the first time including a self-consistent treatment of forward scattering-induced coupling and entanglement of intersecting neutrino/antineutrino trajectories. For the atmospheric neutrino mass-squared difference we find that in the normal (inverted) mass hierarchy the more tangentially-propagating (radially-propagating) neutrinos and antineutrinos can initiate collective, simultaneous medium-enhanced flavor conversion of these particles across broad ranges of energy and propagation direction. Accompanying alterations in neutrino/antineutrino energy spectra and/or fluxes could affect supernova nucleosynthesis and the expected neutrino signal.


title: Density Fluctuations in MHD Turbulence: Spectra, Intermittency and Topology

abstract: We perform three-dimensional (3D) compressible MHD simulations over many dynamical times for an extended range of sonic and Alfven Mach numbers and analyze the statistics of 3D density and 2D column density, which include probability distribution functions, spectra, skewness, kurtosis, She-Leveque exponents, and genus. In order to establish the relation between the statistics of the observables, i.e. column densities, and the underlying 3D statistics of density, we analyze the effects of cloud boundaries. We define the parameter space for 3D measures to be recovered from column densities. In addition, we show that for subsonic turbulence the spectra of density fluctuations are consistent with k^-7/3 in the case of a strong magnetic field and k^-5/3 in the case of a weak magnetic field. For supersonic turbulence we confirm the earlier findings of the shallow spectra of density and Kolmogorov spectra of the logarithm of density. We find that the intermittencies of the density and velocity are very different.


title: A pulsational distance to Omega Centauri based on Near-Infrared Period-Luminosity relations of RR Lyrae stars

abstract: We present new Near-Infrared (J,K) magnitudes for 114 RR Lyrae stars in the globular cluster Omega Cen (NGC 5139) which we combine with data from the literature to construct a sample of 180 RR Lyrae stars with J and K mean magnitudes on a common photometric system. This is presently the largest such sample in any stellar system. We also present updated predictions for J,K-band Period-Luminosity relations for both fundamental and first-overtone RR Lyrae stars, based on synthetic horizontal branch models with metal abundance ranging from Z=0.0001 to Z=0.004. By adopting for the Omega Cen variables with measured metal abundances an alpha-element enhancement of a factor of 3 (about 0.5 dex) with respect to iron we find a true distance modulus of 13.70 (with a random error of 0.06 and a systematic error of 0.06), corresponding to a distance d=5.5 Kpc (with both random and systematic errors equal to 0.03 Kpc). Our estimate is in excellent agreement with the distance inferred for the eclipsing binary OGLEGC-17, but differ significantly from the recent distance estimates based on cluster dynamics and on high amplitude Delta Scuti stars.


title: Sensitivity of an underwater Cerenkov km3 telescope to TeV neutrinos from Galactic Microquasars

abstract: In this paper are presented the results of Monte Carlo simulations on the capability of the proposed NEMO-km$^3$ telescope to detect TeV muon neutrinos from Galactic microquasars. For each known microquasar we compute the number of detectable events, together with the atmospheric neutrino and muon background events. We also discuss the detector sensitivity to neutrino fluxes expected from known microquasars, optimizing the event selection also to reject the background; the number of events surviving the event selection are given. The best candidates are the steady microquasars SS433 and GX339-4 for which we estimate a sensitivity of about $5\cdot10^{-11}$ erg/cm$^2$ s; the predicted fluxes are expected to be well above this sensitivity. For bursting microquasars the most interesting candidates are Cygnus X-3, GRO J1655-40 and XTE J1118+480: their analyses are more complicated because of the stochastic nature of the bursts.


title: Survival of a brown dwarf after engulfment by a red giant star

abstract: (Summary of paper) High resolution optical spectroscopy of the white dwarf WD0137-349 reveals a weak Halpha line in emission, due to a low mass companion in a close orbit. Using this emission feature, and the narrow NLTE core of the white dwarf's Halpha absorption line, we measure the orbital period at 116 minutes, the separation at 0.65Rsun, and the mass ratio m2/m1 = K1/K2 = 0.134+/-0.006. From the mass of the white dwarf (0.39+/-0.035Msun) we derive a mass for the companion of 0.053+/-0.006Msun. This is well below the limit of about 0.075Msun commonly used to distinguish stars from brown dwarfs. The observed infrared flux distribution of WD0137-349 is also consistent with a model of an old brown dwarf spectral type mid T or slightly earlier. The system is detached, and the brown dwarf must have survived a previous phase of common envelope evolution, providing a key observational test for models of this interaction.


title: Pulsar Braking Indices, Glitches and Energy Dissipation In Neutron Stars

abstract: Almost all pulsars with anomalous positive $\ddot \Omega $ measurements (corresponding to anomalous braking indices in the range 5$<n<$100), including all the pulsars with observed large glitches ($\Delta\Omega/\Omega$ $>$ 10$^{-7}$) as well as post glitch or interglitch $\ddot \Omega$ measurements obey the scaling between $\ddot \Omega$ and glitch parameters originally noted in the Vela pulsar. Negative second derivative values can be understood in terms of glitches that were missed or remained unresolved. We discuss the glitch rates and a priori probabilities of positive and negative braking indices according to the model developed for the Vela pulsar. This behavior supports the universal occurrence of a nonlinear dynamical coupling between the neutron star crust and an interior superfluid component. The implied lower limit to dynamical energy dissipation in a neutron star with spindown rate $\dot \Omega$ is $\dot E_{diss}> 1.7 \times 10 ^{-6} \dot E_{rot}$. Thermal luminosities and surface temperatures due to dynamical energy dissipation are estimated for old neutron stars which are spinning down as rotating magnetic dipoles beyond the pulsar death line.


title: An optical search for Low Surface Brightness Galaxies in the Arecibo HI Strip Survey

abstract: In order to estimate the contribution of low surface brightness (LSB) galaxies to the local (z<0.1) galaxy number density, we performed an optical search for LSB candidates in a 15.5 deg^2 part of the region covered by the 65 deg^2 blind Arecibo HI Strip Survey (AHISS). Object detection and galaxy profile fitting were done with analytical algorithms. The detection efficiency and the selection effects were evaluated using large samples of artificial galaxies. Our final catalogue is diameter-limited and contains 306 galaxies with diameters >18" at the limiting surface brightness of mu_{B,lim}=25.2+-0.31 mag/arcsec^2. Of these 306 galaxies, 148 were not catalogued previously. Our results indicate that low surface brightness galaxies contribute at least to 30 % to the local galaxy number density. Without additional distance information, choosing the limiting diameter and the surface brightness at which the diameter is measured is crucial. Depending on these choices, diameter-limited optical catalogues are either biased against LSB galaxies, or contaminated with cosmologically dimmed high surface brightness galaxies, which affects the implied surface brightness distribution. The comparison to the AHISS showed that although optical surveys detect more galaxies per deg^2 than HI surveys, their drawback is the need for spectroscopic follow up observations to derive distances. Blind HI surveys have no diameter limits, but tend to miss gas-poor galaxies and all galaxies which lie outside their redshift limits. HI and optical surveys thus provide complementary information and sample different parts of the LSB galaxy population.


title: The Infrared Emission from the Narrow Line Region

abstract: We present models for the mid- and far- infrared emission from the Narrow Line Region (NLR) of Active Galactic Nuclei (AGN). Using the MAPPINGS III code we explore the effect of typical NLR parameters on the spectral characteristics of the IR emission. These include useful IR emission line ratio diagnostic diagrams for the determination of these parameters, as well as Star formation--AGN mixing diagnostics. We also examine emission line to continuum correlations which would assist in separating the IR emission arising from the NLR from that coming from the inner torus. We find for AGN like NGC 1068 and NGC 4151 that the NLR only contributes ~10% to the total IRAS 25 mum flux, and that other components such as a dusty torus are necessary to explain the total AGN IR emission.


title: Particle Accelerator in Pulsar Magnetospheres: Super Goldreich-Julian Current with Ion Emission from the Neutron Star Surface

abstract: We investigate the self-consistent electrodynamic structure of a particle accelerator in the Crab pulsar magnetosphere on the two-dimensional poloidal plane, solving the Poisson equation for the electrostatic potential together with the Boltzmann equations for electrons, positrons and gamma-rays. If the trans-field thickness of the gap is thin, the created current density becomes sub-Goldreich-Julian, giving the traditional outer-gap solution but with negligible gamma-ray luminosity. As the thickness increases, the created current increases to become super-Goldreich-Julian, giving a new gap solution with substantially screened acceleration electric field in the inner part. In this case, the gap extends towards the neutron star with a small-amplitude positive acceleration field, extracting ions from the stellar surface as a space-charge-limited flow. The acceleration field is highly unscreened in the outer magnetosphere, resulting in a gamma-ray spectral shape which is consistent with the observations.


title: Chaotic zone boundary for low free eccentricity particles near an eccentric planet

abstract: We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Via collisionless particle integration we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semi-major axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7 law has been previously explained by estimating the semi-major at which the first order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first order corotation resonances that have strength proportional to the planet's eccentricity. However, we find the corotation resonance width at low free eccentricity is small. Also the first order resonance width at zero free eccentricity is the same as that for a zero eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.


title: Infrared Properties of Close Pairs of Galaxies

abstract: We discuss spectroscopy and infrared photometry for a complete sample of ~ 800 galaxies in close pairs objectively selected from the CfA2 redshift survey. We use 2MASS to compare near infrared color-color diagrams for our sample with the Nearby Field Galaxy Sample and with a set of IRAS flux-limited pairs from Surace et al. We construct a basic statistical model to explore the physical sources of the substantial differences among these samples. The model explains the spread of near infrared colors and is consistent with a picture where central star formation is triggered by the galaxy-galaxy interaction before a merger occurs. For 160 galaxies we report new, deep JHK photometry within our spectroscopic aperture and we use the combined spectroscopic and photometric data to explore the physical conditions in the central bursts. We find a set of objects with H-K >= 0.45 and with a large F(FIR)/F(H). We interpret the very red H-K colors as evidence for 600-1000 K dust within compact star-forming regions, perhaps similar to super-star clusters identified in individual well-studied interacting galaxies. The galaxies in our sample are candidate ``hidden'' bursts or, possibly, ``hidden'' AGN. Over the entire pair sample, both spectroscopic and photometric data show that the specific star formation rate decreases with the projected separation of the pair. The data suggest that the near infrared color-color diagram is also a function of the projected separation; all of the objects with central near infrared colors indicative of bursts of star formation lie at small projected separation.


title: X-ray Properties of Intermediate-mass Black Holes in Active Galaxies

abstract: We present a pilot study of the X-ray properties of intermediate-mass (~10^5-10^6 M_sun) black holes in active galaxies using the Chandra X-ray telescope. Eight of the 10 active galaxies are detected with a significance of at least 3 sigma, with X-ray luminosities in the range L_(0.5-2 keV) ~ 10^41-10^43 ergs/s. The optical-to-X-ray flux ratios are consistent with expectations, given the known correlations between alpha_ox and ultraviolet luminosity, while a couple of objects appear to be anomalously X-ray weak. The range of 0.5--2 keV photon indices we measure, 1 < Gamma_s < 2.7, is entirely consistent with values found in samples of more luminous sources with more massive black holes. Black hole mass evidently is not a primary driver of soft X-ray spectral index. On the other hand, we do find evidence for a correlation between X-ray power-law slope and both X-ray luminosity and Eddington ratio, which may suggest that X-ray emission mechanisms weaken at high Eddington ratio. Such a weakening may explain the X-ray weakness of one of our most optically luminous objects.


title: Modelling the Pan-Spectral Energy Distribution of Starburst Galaxies: III. Emission Line Diagnostics of Ensembles of Evolving HII Regions

abstract: We build, as far as theory will permit, self consistent model HII regions around central clusters of aging stars. These produce strong emission line diagnostics applicable to either individual HII regions in galaxies, or to the integrated emission line spectra of disk or starburst galaxies. The models assume that the expansion and internal pressure of individual HII regions is driven by the net input of mechanical energy from the central cluster, be it through winds or supernova events. This eliminates the ionization parameter as a free variable, replacing it with a parameter which depends on the ratio of the cluster mass to the pressure in the surrounding interstellar medium. These models explain why HII regions with low abundances have high excitation, and demonstrate that at least part of the warm ionized medium is the result of overlapping faint, old, large, and low pressure HII regions. We present line ratios (at both optical and IR wavelengths) which provide reliable abundance diagnostics for both single HII regions or for integrated galaxy spectra, and we find a number that can be used to estimate the mean age of the cluster stars exciting individual HII regions.


title: Shock-cloud interaction in the Vela SNR II. Hydrodynamic model

abstract: In the framework of the study of the X-ray and optical emission in supernova remnants we focus on an isolated X-ray knot in the northern rim of the Vela SNR (Vela FilD), whose X-ray emission has been studied and discussed in Paper I. We aim at understanding the physical origin of the X-ray and optical emission in FilD, at understanding the role of the different physical processes at work, and at obtaining a key for the interpretation of future X-ray observations of SNRs. To this end we have pursued an accurate ``forward'' modeling of the interaction of the Vela SNR shock with an ISM cloud. We perform hydrodynamic simulations and we directly compare the observables synthesized from the simulations with the data. We explore four different model setups, choosing the values of the physical parameters on the basis of our preliminary analysis of the X-ray data. We synthesize X-ray emission maps and spectra filtered through the XMM-Newton EPIC-MOS instrumental response. The impact of a shock front at 6 million Kelvin on an elliptical cloud 30 times denser than the ambient medium describes well the shock-cloud interaction processes in the Vela FilD region in terms of spectral properties and morphology of the X-ray and optical emission. The bulk of the X-ray emission in the FilD knot originates in the cloud material heated by the transmitted shock front, but significant X-ray emission is also associated to the cloud material, which evaporates, as an effect of thermal conduction, in the intercloud medium. The physical origin of the FilD optical emission is associated to thermal instabilities. In the FilD knot the X-ray emission associated to the reflected shock front is negligible.


title: Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix

abstract: AIMS. The maximum-likelihood method is the standard approach to obtain model fits to observational data and the corresponding confidence regions. We investigate possible sources of bias in the log-likelihood function and its subsequent analysis, focusing on estimators of the inverse covariance matrix. Furthermore, we study under which circumstances the estimated covariance matrix is invertible. METHODS. We perform Monte-Carlo simulations to investigate the behaviour of estimators for the inverse covariance matrix, depending on the number of independent data sets and the number of variables of the data vectors. RESULTS. We find that the inverse of the maximum-likelihood estimator of the covariance is biased, the amount of bias depending on the ratio of the number of bins (data vector variables), P, to the number of data sets, N. This bias inevitably leads to an -- in extreme cases catastrophic -- underestimation of the size of confidence regions. We report on a method to remove this bias for the idealised case of Gaussian noise and statistically independent data vectors. Moreover, we demonstrate that marginalisation over parameters introduces a bias into the marginalised log-likelihood function. Measures of the sizes of confidence regions suffer from the same problem. Furthermore, we give an analytic proof for the fact that the estimated covariance matrix is singular if P>N.


title: Model Atmosphere Analysis of the Weakly Magnetic DZ White Dwarf G165-7

abstract: A reanalysis of the strongly metal-blanketed DZ white dwarf G165-7 is presented. An improved grid of model atmospheres and synthetic spectra is used to analyze BVRI, JHK, and ugriz photometric observations as well as a high quality Sloan Digital Sky Survey spectrum covering the energy distribution from 3600 to 9000 A. The detection of splitting in several lines of Ca, Na, and Fe, suggesting a magnetic field of Bs~650 kG, is confirmed by spectropolarimetric observations that reveal as much as +/- 7.5% circular polarization in many of the absorption lines, most notably Na, Mg, and Fe. Our combined photometric and spectroscopic fit yields Teff = 6440 K, log g=7.99, log H/He=-3.0 and log Ca/He=-8.1. The other heavy elements have solar ratios with respect to calcium, with the exception of Na and Cr that had to be reduced by a factor of two and three, respectively. A crude polarization model based upon the observed local spectral flux gradient yields a longitudinal field of 165 kG, consistent with the mean surface field inferred from the Zeeman splitting. The inclusion of this weak magnetic field in our synthetic spectrum calculations, even in an approximate fashion, is shown to improve our fit significantly.


title: Spectral State Transitions of the Ultraluminous X-ray Sources X-1 and X-2 in NGC 1313

abstract: X-ray spectral state transitions are a key signature of black hole X-ray binaries and follow a well-defined pattern. We examined 12 XMM-Newton observations of the nearby spiral galaxy NGC 1313, which harbors two compact ultraluminous X-ray sources (ULXs), X-1 and X-2, in order to determine if the state transitions in ULXs follow the same pattern. For both sources, the spectra were adequately fitted by an absorbed power-law with the addition of a low temperature (kT=0.1~0.3 keV) disk blackbody component required in 6 of the 12 observations. As the X-ray luminosity of X-1 increases to a maximum at 3x10^40 ergs/s, the power-law photon index softens to 2.5-3.0. This behavior is similar to the canonical spectral state transitions in Galactic black hole binaries, but the source never enters the high/soft or thermal dominant state and instead enters the steep power-law state at high luminosities. X-2 has the opposite behavior and appears to be in the hard state, with a photon index of Gamma=1.7-2.0 at high luminosity, but can soften to Gamma=2.5 at the lower luminosities.


title: Spin Exchange Rates in Electron-Hydrogen Collisions

abstract: The spin temperature of neutral hydrogen, which determines the 21 cm optical depth and brightness temperature, is set by the competition between radiative and collisional processes. In the high-redshift intergalactic medium, the dominant collisions are typically those between hydrogen atoms. However, collisions with electrons couple much more efficiently to the spin state of hydrogen than do collisions with other hydrogen atoms and thus become important once the ionized fraction exceeds ~1%. Here we compute the rate at which electron-hydrogen collisions change the hydrogen spin. Previous calculations included only S-wave scattering and ignored resonances near the n=2 threshold. We provide accurate results, including all partial wave terms through the F-wave, for the de-excitation rate at temperatures T_K < 15,000 K; beyond that point, excitation to n>=2 hydrogen levels becomes significant. Accurate electron-hydrogen collision rates at higher temperatures are not necessary, because collisional excitation in this regime inevitably produces Lyman-alpha photons, which in turn dominate spin exchange when T_K > 6200 K even in the absence of radiative sources. Our rates differ from previous calculations by several percent over the temperature range of interest. We also consider some simple astrophysical examples where our spin de-excitation rates are useful.


title: On the Physics of Type I X-ray Bursts on Accreting Neutron Stars at High Accretion Rates

abstract: We investigate the effect of the hot CNO cycle breakout reaction 15O(alpha,gamma)19Ne on the occurrence of type I X-ray bursts on accreting neutron stars. For f_rp <~ 0.1, where f_rp is a dimensionless factor by which we multiply the 15O(alpha,gamma)19Ne reaction rate of Caughlan & Fowler (1988), our model predicts that bursts should occur only for accretion rates below a critical value of approximately 0.3 times the Eddington limit. This agrees with observations. For larger values of f_rp, including the standard choice f_rp = 1, the model switches to a new regime in which bursts occur all the way up to roughly the Eddington limit. Since the latter regime disagrees with observations, we suggest that the true 15O(alpha,gamma)19Ne reaction rate is lower than usually assumed.


title: Dynamical Evolution of Globular Clusters in Hierarchical Cosmology

abstract: We test the hypothesis that metal-poor globular clusters form within disk galaxies at redshifts z>3. Numerical simulations demonstrate that giant gas clouds, which are cold and dense enough to produce massive star clusters, assemble naturally in hierarchical models of galaxy formation at high redshift. Do model clusters evolve into observed globular clusters or are they disrupted before present as a result of the dynamical evolution? To address this question, we calculate the orbits of model clusters in the time-variable gravitational potential of a Milky Way-sized galaxy, using the outputs of a cosmological N-body simulation. We find that at present the orbits are isotropic in the inner 50 kpc of the Galaxy and preferentially radial at larger distances. All clusters located outside 10 kpc from the center formed in satellite galaxies, some of which are now tidally disrupted and some of which survive as dwarf galaxies. The spatial distribution of model clusters is spheroidal and the fit to the density profile has a power-law slope of 2.7. The combination of two-body relaxation, tidal shocks, and stellar evolution drives the evolution of the cluster mass function from an initial power law to a peaked distribution, in agreement with observations. However, not all initial conditions and not all evolution scenarios are consistent with the observed mass function of the Galactic globular clusters. The successful models require the average cluster density to be constant initially for clusters of all mass and to remain constant with time. Synchronous formation of all clusters at a single epoch (z=4) and continuous formation over a span of 1.6 Gyr (between z=9 and z=3) are both consistent with the data. For both formation scenarios, we provide online catalogs of the main physical properties of model clusters.


title: Stellar Encounters with Massive Star-Disk Systems

abstract: The dense, clustered environment in which massive stars form can lead to interactions with neighboring stars. It has been hypothesized that collisions and mergers may contribute to the growth of the most massive stars. In this paper we extend the study of star-disk interactions to explore encounters between a massive protostar and a less massive cluster sibling using the publicly available SPH code GADGET-2. Collisions do not occur in the parameter space studied, but the end state of many encounters is an eccentric binary with a semi-major axis ~ 100 AU. Disk material is sometimes captured by the impactor. Most encounters result in disruption and destruction of the initial disk, and periodic torquing of the remnant disk. We consider the effect of the changing orientation of the disk on an accretion driven jet, and the evolution of the systems in the presence of on-going accretion from the parent core.


title: Prospects for the Characterization and Confirmation of Transiting Exoplanets via the Rossiter-McLaughlin Effect

abstract: The Rossiter-McLaughlin (RM) effect is the distortion of stellar spectral lines that occurs during eclipses or transits, due to stellar rotation. We assess the future prospects for using the RM effect to measure the alignment of planetary orbits with the spin axes of their parent stars, and to confirm exoplanetary transits. We compute the achievable accuracy for the parameters of interest, in general and for the 5 known cases of transiting exoplanets with bright host stars. We determine the requirements for detecting the effects of differential rotation. For transiting planets with small masses or long periods (as will be detected by forthcoming satellite missions), the velocity anomaly produced by the RM effect can be much larger than the orbital velocity of the star. For a terrestrial planet in the habitable zone of a Sun-like star found by the Kepler mission, it will be difficult to use the RM effect to confirm transits with current instruments, but it still may be easier than measuring the spectroscopic orbit.


title: Search for Oxygen Emission from Warm-Hot Intergalactic Medium around A2218 with Suzaku

abstract: We searched for redshifted O emission lines from the possible warm-hot intergalactic medium (WHIM) surrounding the cluster of galaxies A2218 at z=0.1756 using the XIS instrument on Suzaku. This cluster is thought to have an elongated structure along the line of sight based on previous studies. We studied systematic uncertainties in the spectrum of the Galactic emission and in the soft X-ray response of the detectors due to the contamination building up on the XIS filters. We detected no significant redshifted O lines, and set a tight constraint on the intensity with upper limits for the surface brightness of OVII and OVIII lines of 1.1 x 10^-7 and 3.0 x 10^-7 photons cm^-2 s^-1 arcmin^-2, respectively. These upper limits are significantly lower than the previously reported fluxes from the WHIM around other clusters of galaxies. We also discuss the prospect for the detection of the WHIM lines with Suzaku XIS in the future.


title: A deep mosaic of [O III] 5007 A CCD images of the environment of the LBV star P Cygni

abstract: A mosaic of six, deep, CCD images in the light of the [O III] 5007 A nebular emission line has been obtained with the 1.3-m Skinakas (Crete) telescope of the filamentary nebulosity surrounding P Cygni. The [O III] 5007 A line discriminates against confusing galactic H II regions along the same sight-lines and the new mosaic did not include the 4.8 mag. central star; a source of artifacts in the previous lower angular resolution observations. New giant `lobes' and `shells' are found to be clustered around P Cygni which must be the relics of historic eruptions between 2400 and up to ~10^5 yr ago.


title: The Fluorescence Detector of the Pierre Auger Observatory - A Calorimeter for UHECR

abstract: The Pierre Auger Observatory is a hybrid detector for ultrahigh energy cosmic rays (UHECR) with energies above 10$^{18.5}$ eV. Currently the first part of the Observatory nears completion in the southern hemisphere in Argentina. One detection technique uses over 1600 water Cherenkov tanks at ground where samples of secondary particles of extensive air showers (EAS) are detected. The second technique is a calorimetric measurement of the energy deposited by EAS in the atmosphere. Charged secondary particles of EAS lose part of their energy in the atmosphere via ionization. The deposited energy is converted into excitation of molecules of the air and afterwards partly emitted as fluorescence light mainly from nitrogen in the wavelength region between 300 and 400 nm. This light is observed with 24 fluorescence telescopes in 4 stations placed at the boundary of the surface array. This setup provides a combined measurement of the longitudinal shower development and the lateral particle distribution at ground of the same event. Details on the fluorescence technique and the necessary atmospheric monitoring will be presented, as well as first physics results on UHECR.


title: 3D Lyman-alpha radiation transfer. I. Understanding Lyman-alpha line profile morphologies

abstract: Using a Monte Carlo technique, we have developed a 3D lyman-alpha radiation transfer code allowing for prescribed arbitrary hydrogen density, ionisation, temperature structures, and dust distribution, and arbitrary velocity fields and UV photon sources. We have examined the lyman-alpha line profiles predicted for several simple geometrical configurations and their dependence on the main input parameters. Overall, we find line profiles reaching from doubly peaked symmetric emission to symmetric Voigt (absorption) in static configurations with increasing dust content, and asymmetric red-(blue-) shifted emission lines with a blue (red) counterpart ranging from absorption to emission (with increasing line/continuum strength) in expanding (infalling) media. The following results are of interest for the interpretation of lya profiles from galaxies. 1/ Standard lya absorption line fitting of global spectra of galaxies may lead to an underestimate of the true hydrogen column density in certain geometrical conditions. 2/ Normal (inverted) P-Cygni like lya profiles can be obtained in expanding (infalling) media from objects without any intrinsic lya emission, as a natural consequence of radiation transfer effects. 3/ The formation and the detailed shape of lya profiles resulting from expanding shells has been thoroughly revised: for sufficiently large column densities, the position of the main lya emission peak is redshifted by twice the expansion velocity.This is in excellent agreement with the observations of z~3 LBGs. This finding indicates also that large scale, fairly symmetric shell structures must be a good description for the outflows in LBGs.(shortened abstract)


title: Wolf-Rayet Populations at High Metallicity

abstract: Observed properties of Wolf-Rayet stars at high metallicity are reviewed. Wolf-Rayet stars are more common at higher metallicity, as a result of stronger mass-loss during earlier evolutionary phases with late WC subtypes signatures of solar metallicity or higher. Similar numbers of early (WC4-7) and late (WC8-9) stars are observed in the Solar neighbourhood, whilst late subtypes dominate at higher metallicities, such as Westerlund 1 in the inner Milky Way and in M83. The observed trend to later WC subtype within metal-rich environments is intimately linked to a metallicity dependence of WR stars, in the sense that strong winds preferentially favour late subtypes. This has relevance to (a) the upper mass limit in metal-rich galaxies such as NGC 3049, due to softer ionizing fluxes from WR stars at high metallicity; (b) evolutionary models including a WR metallicity dependence provide a better match to the observed N(WC)/N(WN) ratio. The latter conclusion partially rests upon the assumption of constant line luminosities for WR stars, yet observations and theoretical atmospheric models reveal higher line fluxes at high metallicity.


title: 2.3 micron CO emission and absorption from young high-mass stars in M17

abstract: We are studying the extremely young cluster of M17 to investigate the birth of high-mass stars and the initial mass function. Deep JHKL imaging and K-band spectroscopy from the VLT of 201 stars toward the cluster is presented. The majority of 104 stars show the CO band-head in absorption. Half of them emit X-rays and/or have infrared excess, indicative of very young objects. Their intrinsic IR luminosity is compatible with intermediate and high-mass pre-main sequence stars. Nine additional stars have the CO feature in emission, while sixty sources are lacking any stellar spectral feature due to veiling by circumstellar dust. We suggest that CO absorption is - as in the case of low-mass stars - also a common feature during the early evolution of stars with higher masses. According to model calculations the observed CO absorption is most likely a sign of heavily accreting protostars with mass accretion rates above 10^-5 solar masses/yr.


title: The Galactic WN stars: Spectral analyses with line-blanketed model atmospheres versus stellar evolution models with and without rotation

abstract: CONTEXT: Very massive stars pass through the Wolf-Rayet (WR) stage before they finally explode. Details of their evolution have not yet been safely established, and their physics are not well understood. Their spectral analysis requires adequate model atmospheres, which have been developed step by step during the past decades and account in their recent version for line blanketing by the millions of lines from iron and iron-group elements. However, only very few WN stars have been re-analyzed by means of line-blanketed models yet. AIMS: The quantitative spectral analysis of a large sample of Galactic WN stars with the most advanced generation of model atmospheres should provide an empirical basis for various studies about the origin, evolution, and physics of the Wolf-Rayet stars and their powerful winds. METHODS: We analyze a large sample of Galactic WN stars by means of the Potsdam Wolf-Rayet (PoWR) model atmospheres, which account for iron line blanketing and clumping. The results are compared with a synthetic population, generated from the Geneva tracks for massive star evolution. RESULTS: We obtain a homogeneous set of stellar and atmospheric parameters for the Galactic WN stars, partly revising earlier results. CONCLUSIONS: Comparing the results of our spectral analyses of the Galactic WN stars with the predictions of the Geneva evolutionary calculations, we conclude that there is rough qualitative agreement. However, the quantitative discrepancies are still severe, and there is no preference for the tracks that account for the effects of rotation. It seems that the evolution of massive stars is still not satisfactorily understood.


title: Radio emission from shell-type supernova remnants

abstract: The evolution of the radio emission of shell-type Supernova remnants (SNRs) is modeled within the framework of the simple and commonly used assumptions that the mechanism of diffusive shock acceleration (DSA) is responsible for generating radio emitting electrons and that the magnetic field is the typical interstellar field compressed at the shock. It is considered that electrons are injected into the mechanism in test-particle regime directly from the high energy tail of the downstream Maxwellian distribution function. The model can be applied to most of the observed SNRs. It is shown that the model successfully explains the many averaged observational properties of evolved shell-type SNRs. In particular, the radio surface brightness ($\Sigma$) evolves with diameter as $\sim D^{-(0.3 \div 0.5)}$, while the bounding shock is strong (Mach number is ${\mathcal M} \geq10$), followed by steep decrease (steeper than $\sim D^{-4.5}$) for ${\cal M} <10$. Such evolution of the surface brightness with diameter and its strong dependence on the environmental parameters strongly reduce the usefulness of $\Sigma - D$ relations as a tool for determining the distances to SNRs. The model predicts no radio emission from SNRs in the late radiative stage of evolution and the existence of radio-quiet but relatively active SNRs is possible. Our model easily explains very large-diameter radio sources such as the Galactic Loops and the candidates for Hypernova radio remnants. The model predicts that most of the observed SNRs are located in a tenuous phase of the ISM. From the comparison of the model results with the statistics of evolved shell-type SNRs, we were able to estimate the fraction of electrons accelerated from the thermal pool in the range $(3\div 11) \times 10^{- 4}$.


title: A spectroscopic atlas of post-AGB stars and planetary nebulae selected from the IRAS Point Source Catalogue

abstract: Aims: We study the optical spectral properties of a sample of stars showing far infrared colours similar to those of well-known planetary nebulae. The large majority of them were unidentified sources or poorly known in the literature at the time when this spectroscopic survey started, some 15 years ago. Methods: We present low-resolution optical spectroscopy, finding charts and improved astrometric coordinates of a sample of 253 IRAS sources. Results: We have identified 103 sources as post-AGB stars, 21 as ``transition sources'', and 36 as planetary nebulae, some of them strongly reddened. Among the rest of sources in the sample, we were also able to identify 38 young stellar objects, 5 peculiar stars, and 2 Seyfert galaxies. Up to 49 sources in our spectroscopic sample do not show any optical counterpart, and most of them are suggested to be heavily obscured post-AGB stars, rapidly evolving on their way to becoming planetary nebulae. Conclusions: An analysis of the galactic distribution of the sources identified as evolved stars in the sample is presented together with a study of the distribution of these stars in the IRAS two-colour diagram. Finally, the spectral type distribution and other properties of the sources identified as post-AGB in this spectroscopic survey are discussed in the framework of stellar evolution.


title: CG Draconis, a particularly active dwarf nova

abstract: The dwarf nova CG Dra was intensively monitored during 2005 as part of the Variable Star Sections Recurrent Objects Programme and seven outbursts were detected. These observations were combined with others from the BAA and AAVSO databases which confirmed that the star shows frequent outbursts with a period of around 11 days. Two types of outburst have been detected. Short outbursts lasting about 4 days and long outbursts lasting about 8 days.


title: Determination of the superhump period of the dwarf nova V701 Tau during the 2005 December superoutburst

abstract: We report new measurements of the superhump period of the UGSU-type dwarf nova V701 Tau during the 2005 December superoutburst. Using unfiltered time series CCD observations on 3 nights, we determine a probable superhump period Psh = 0.0690, plus or minus 0.0002d, but note that our data also permit a possible shorter period of 0.0663, plus or minus 0.0002d. The longer period agrees with the value measured during the first recorded superoutburst in 1995.


title: Spitzer Observations of Var Her 04: Possible Detection of Dust Formation in a Super-Outbursting TOAD

abstract: We present four MIPS (24 \micron) and two IRAC (3.6, 4.5, 5.8, and 8.0 \micron) Spitzer observations of the newly discovered Tremendous Outburst Amplitude Dwarf nova (TOAD) Var Her 04 during decline from super-outburst. The four MIPS observations span 271 days and the two IRAC observations span 211 days. Along the line-of-sight to Var Her 04, there is a foreground M-star within 1\arcsec of the variable; as a result, all of the Spitzer photometry presented in this paper is a blend of the foreground M-star and Var Her 04. We estimate the quiescent level of the TOAD to be $\Delta V=4-5$ magnitudes below that of the M-star. Based upon the spectral energy distribution and the 2MASS colors, we find the M-star to be an M3.5V dwarf at a distance of 80-130 pc. Based upon its outburst amplitude and quiescent apparent magnitude, we estimate the distance to Var Her 04 to be 200-400 pc, suggesting that the line-of-sight foreground star is physically unrelated to the cataclysmic variable. All of the Spitzer photometry is consistent with the photospheric emission of the line-of-sight M3.5V star, except for one 24 \micron observation obtained after the variable re-brightened. This 24 \micron flux density is 75 $\mu$Jy ($4\sigma$) above the preceding and following MIPS observations. We tentatively suggest that the mid-infrared brightening of 75 $\mu$Jy may be associated with a dust formation event in the super-outburst ejecta. Assuming a dust temperature of 100-400 K, we have estimated the amount of dust required. We find $10^{-13}-10^{-11}$ M$_\odot$ of dust is needed, consistent with amounts of mass ejection in TOADs expected during super-outburst, and possibly making TOADs important contributors to the recycling of the interstellar medium.


title: X-ray Spectral Variations of U Gem from Quiescence to Outburst

abstract: In this paper we report the discovery of a high energy component of the X-ray spectra of U Gem, which can be observed while the source is in outburst. We used Chandra and XMM-Newton observations to compare the quiescence and outburst X-ray spectra of the source. The additional component may be the result of the reflection of X-rays emitted from an optically thin plasma close to the white dwarf, from the optically thick boundary layer during the outburst. Another possible explanation is that some magnetically channeled accretion may occur onto the equatorial belt of the primary causing shocks similar to the ones in the intermediate polars as it was suggested by \citep{w2002}. We have also found a timing structure at about 73 mHz ($\sim$13.7 s.) in the RXTE observation, resembling dwarf novae oscillations (DNOs).


title: The intergalactic propagation of ultra-high energy cosmic ray nuclei

abstract: We investigate the propagation of ultra-high energy cosmic ray nuclei (A = 1-56) from cosmologically distant sources through the cosmic radiation backgrounds. Various models for the injected composition and spectrum and of the cosmic infra-red background are studied using updated photo-disintegration cross-sections. The observational data on the spectrum and the composition of ultra-high energy cosmic rays are jointly consistent with a model where all of the injected primary cosmic rays are iron nuclei (or a mixture of heavy and light nuclei).


title: The Spiral Host Galaxy of the Double Radio Source 0313-192

abstract: We present new Hubble, Gemini-S, and Chandra observations of the radio galaxy 0313-192, which hosts a 350-kpc double source and jets, even though previous data have suggested that it is a spiral galaxy. We measure the bulge scale and luminosity, radial and vertical profiles of disk starlight, and consider the distributions of H II regions and absorbing dust. In each case, the HST data confirm its classification as an edge-on spiral galaxy, the only such system known to produce such an extended radio source of this kind. The Gemini near-IR images and Chandra spectral fit reveal a strongly obscured central AGN, seen through the entire ISM path length of the disk and showing X-ray evidence of additional absorption from warm or dense material close to the central object. We consider several possible mechanisms for producing such a rare combination of AGN and host properties, some combination of which may be at work. These include an unusually luminous bulge (suggesting a black hole of mass 0.5-0.9 billion solar masses), orientation of the jets near the pole of the gas-rich disk, and some evidence of a weak gravitational interaction which has warped the disk and could have enhanced fuelling of the central engine. An X-ray counterpart of the kiloparsec-scale radio jet emerges to the south; jet/counterjet limits in both radio and X-ray allow them to be symmetric if seen more than 15 degrees from the plane of the sky, still consistent with the jet axes being within ~30 degrees of the poles of the gas-rich galaxy disk. A linear or disklike emission-line structure is seen around thenucleus, inclined by ~20 degrees to the stellar disk but nearly perpendicular to the jets; this may represent the aftermath of a galaxy encounter, where gas is photoionized by the nuclear continuum.


title: The Energy Transfer Process in Planetary Flybys

abstract: We illustrate the energy transfer during planetary flybys as a function of time using a number of flight mission examples. The energy transfer process is rather more complicated than a monotonic increase (or decrease) of energy with time. It exhibits temporary maxima and minima with time which then partially moderate before the asymptotic condition is obtained. The energy transfer to angular momentum is exhibited by an approximate Jacobi constant for the system. We demonstrate this with flybys that have shown unexplained behaviors: i) the possible onset of the "Pioneer anomaly" with the gravity assist of Pioneer 11 by Saturn to hyperbolic orbit (as well as the Pioneer 10 hyperbolic gravity assist by Jupiter) and ii) the Earth flyby anomalies of small increases in energy {\it in the geocentric system} (Galileo-I, NEAR, and Rosetta, in additioon discussing the Cassini and Messenger flybys). Perhaps some small, as yet unrecognized effect in the energy-transfer process can shed light on these anomalies.


title: A relativistic model of the radio jets in 3C 296

abstract: We present new, deep 8.5-GHz VLA observations of the nearby, low-luminosity radio galaxy 3C 296 at resolutions from 0.25 to 5.5 arcsec. These show the intensity and polarization structures of the twin radio jets in detail. We also derive the spectral-index and rotation-measure distributions. We apply our intrinsically symmetrical, decelerating relativistic jet model to the new observations. An optimized model accurately fits the data in both total intensity and linear polarization. We infer that the jets are inclined by 58 deg to the line of sight. On-axis, they decelerate from a (poorly-constrained) initial velocity of 0.8c to approximately 0.4c around 5 kpc from the nucleus, the velocity thereafter remaining constant. The speed at the edge of the jet is low everywhere. The longitudinal profile of proper emissivity has three principal power-law sections and the emission is centre-brightened. Our observations rule out a globally-ordered, helical magnetic-field configuration. Instead, we model the field as random on small scales but anisotropic, with toroidal and longitudinal components only. The ratio of longitudinal to toroidal field falls with distance along the jet, qualitatively but not quantitatively as expected from flux freezing, so that the field is predominantly toroidal far from the nucleus. The toroidal component is relatively stronger at the edges of the jet. A simple adiabatic model fits the emissivity evolution only in the outer region after the jets have decelerated and recollimated; closer to the nucleus, it predicts far too steep an emissivity decline with distance. We also interpret the morphological differences between brightness enhancements ("arcs") in the main and counter-jets as an effect of relativistic aberration. (Abridged)


title: Discovery of twin kHz QPOs in the peculiar X-ray binary Circinus X-1

abstract: We report the discovery with the RXTE/PCA of twin kHz QPOs in Cir X-1. Eleven cases of simultaneous double QPOs occurred, with significances of up to 6.3 and 5.5 sigma and centroid frequencies ranging between approximately 56-225 and 230-500 Hz for the two QPO peaks, respectively, i.e., for the most part at frequencies well below those of other sources. The QPO properties clearly indicate that these double peaks are the kHz QPOs known from low magnetic field neutron stars, and not black-hole high-frequency QPOs, confirming that Cir X-1 is a neutron star. The kHz QPO peak separation varies over a wide range, \~175-340 Hz, and increases with QPO frequency. This is contrary to what is seen in other sources but agrees with predictions of the relativistic precession model and Alfven wave models; beat-frequency models require modification to accommodate this. In other observations single kHz QPOs can be seen down to frequencies as low as ~12 Hz, as well as a strong low-frequency (LF) QPO between 1 and 30 Hz. The relations between the frequencies of the kHz QPOs and the LF QPO are in good agreement with those found previously in Z sources, confirming that Cir X-1 may be a peculiar Z source. We suggest that the low frequencies of the kHz QPOs in Cir X-1 and to a lesser extent in (other) Z sources might be due to a relatively stronger radial inflow to the neutron star than in other kHz QPO sources.


title: Upper Bound on the Dark Matter Total Annihilation Cross Section

abstract: We consider dark matter annihilation into Standard Model particles and show that the least detectable final states, namely neutrinos, define an upper bound on the total cross section. Calculating the cosmic diffuse neutrino signal, and comparing it to the measured terrestrial atmospheric neutrino background, we derive a strong and general bound. This can be evaded if the annihilation products are dominantly new and truly invisible particles. Our bound is much stronger than the unitarity bound at the most interesting masses, shows that dark matter halos cannot be significantly modified by annihilations, and can be improved by a factor of 10--100 with existing neutrino experiments.


title: An Upper Limit to the Degree of Evolution Between Supermassive Black Holes and their Host Galaxies

abstract: We describe a model-independent integral constraint which defines an upper limit to the allowed degree of evolution in the ratio of black hole (BH) mass to host galaxy luminosity or mass, as a function of redshift. Essentially, if the BH/host ratio is excessive at redshift z, then it would imply that the total mass density in BHs above some M_BH(min) is larger at that redshift than at z=0, which is impossible. This argument requires no knowledge of host or BH properties, only a lower limit to the observed luminosity density in the brightest galaxies at some z. We calculate this upper limit from a variety of luminosity and mass functions in different bands from redshifts z=0-2. We show that it is consistent with passive evolution of spheroid populations (with a fixed M_BH/M_host relation) in all cases, and provides tighter constraints than have generally been obtained previously, ruling out at >6 sigma observational and theoretical estimates suggesting that M_BH/M_host was significantly larger at high redshifts than locally, although relatively weak (factor ~2 by z=2) evolution is still allowed. We discuss a variety of possible 'loopholes' or changes in the BH/host populations and correlations, and show that they typically lower the upper limits and strengthen our conclusions.


title: Two-Photon Absorption by Hydrogen Molecules: Origin of the 2175 Angstrom Astronomical Band?

abstract: The near UV spectra of OB stars in our galaxy are often dominated by a very broad extinction band peaking at ~2175 Angstroms. Forty years after its discovery, the origin of this band remains unknown, although interstellar dust particles are generally assumed to be the carriers. Here we report that two-photon absorption by hydrogen molecules in gaseous clouds enveloping OB stars can lead to a band peaking at ~2175 Angstroms. We present astronomical spectral evidence supporting our proposal that this nonlinear absorption mechanism accounts for the lambda2175 feature.


title: Accurate M Dwarf Metallicities from Spectral Synthesis: A Critical Test of Model Atmospheres

abstract: We describe a method for accurately determining M dwarf metallicities with spectral synthesis based on abundance analyses of visual binary stars. We obtained high resolution, high signal-to-noise spectra of each component of five visual binary pairs at McDonald Observatory. The spectral types of the components range from F7 to K3 V for the primaries and M0.5 to M3.5 V for the secondaries. We have determined the metallicities of the primaries differentially with respect to the Sun by fitting synthetic spectra to Fe \textsc{i} line profiles in the observed spectra. In the course of our analysis of the M dwarf secondaries we have made significant improvements to the PHOENIX cool-star model atmospheres and the spectrum analysis code MOOG. Our analysis yields a RMS deviation of 0.11 dex in metallicity values between the binary pairs. We estimate the uncertainties in the derived stellar parameters for the M dwarfs to be 48 K, 0.10 dex, 0.12 dex, 0.15 km s$^{-1}$, and 0.20 km s$^{-1}$ for {T}$_{eff}$, log \textsl{g}, [M/H], $\xi$, and $\eta$ respectively. Accurate stellar evolutionary models are needed to progress further in the analysis of cool-star spectra; the new model atmospheres warrant recalculation of the evolutionary models.


title: Atmospheres and Winds of PN Central Stars

abstract: The progress over the last years in modelling the atmospheres and winds of PN central stars is reviewed. We discuss the effect of the inclusion of the blanketing by millions of metal lines in NLTE on the diagnostics of photospheric and stellar wind lines, which can be used to determine stellar parameters such as effective temperature, gravity, radius, mass loss rate and distance. We also refer to recent work on the winds of massive O-type stars, which indicates that their winds are possibly inhomogeneous and clumped. We investigate implications from this work on the spectral diagnostics of PN central stars and introduce a method to determine wind clumping factors from the relative strengths of Halpha and HeII 4686. Based on new results we discuss the wind properties of CSPN.


title: Tidal Tails Test the Equivalence Principle in the Dark Sector

abstract: Satellite galaxies currently undergoing tidal disruption offer a unique opportunity to constrain an effective violation of the equivalence principle in the dark sector. Theories in which cold dark matter (CDM) couples to a light scalar field naturally lead to a long-range force between dark matter particles. An inverse-square-law force of this kind would manifest itself as a violation of the equivalence principle in the dynamics of CDM compared to baryons in the form of gas or stars. In a previous paper, we showed that an attractive force would displace stars outwards from the bottom of the satellite's gravitational potential well, leading to a higher fraction of stars being disrupted from the tidal bulge further from the Galactic center. Since stars disrupted from the far (near) side of the satellite go on to form the trailing (leading) tidal stream, an attractive dark-matter force will produce a relative enhancement of the trailing stream compared to the leading stream. This distinctive signature of a dark-matter force might be detected through detailed observations of the tidal tails of a disrupting satellite, such as those recently performed by the Two-Micron All-Sky Survey (2MASS) and Sloan Digital Sky Survey (SDSS) on the Sagittarius (Sgr) dwarf galaxy. Here we show that this signature is robust to changes in our models for both the satellite and Milky Way, suggesting that we might hope to search for a dark-matter force in the tidal features of other recently discovered satellite galaxies in addition to the Sgr dwarf.


title: Detecting sterile dark matter in space

abstract: Space-based instruments provide new and, in some cases, unique opportunities to search for dark matter. In particular, if dark matter comprises sterile neutrinos, the x ray detection of their decay line is the most promising strategy for discovery. Sterile neutrinos with masses in the keV range could solve several long-standing astrophysical puzzles, from supernova asymmetries and the pulsar kicks to star formation, reionization, and baryogenesis. The best current limits on sterile neutrinos come from Chandra and XMM-Newton. Future advances can be achieved with a high-resolution x-ray spectrometry in space.