title: A statistical analysis of X-ray variability in pre-main sequence objects of the Taurus Molecular Cloud

abstract: This work is part of a systematic X-ray survey of the Taurus star forming complex with XMM-Newton. We study the time series of all X-ray sources associated with Taurus members, to statistically characterize their X-ray variability, and compare the results to those for pre-main sequence stars in the Orion Nebula Cluster and to expectations arising from a model where all the X-ray emission is the result of a large number of stochastically occurring flares. We find that roughly half of the detected X-ray sources show variability above our sensitivity limit, and in ~ 26 % of the cases this variability is recognized as flares. Variability is more frequently detected at hard than at soft energies. The variability statistics of cTTS and wTTS are undistinguishable, suggesting a common (coronal) origin for their X-ray emission. We have for the first time applied a rigorous maximum likelihood method in the analysis of the number distribution of flare energies on pre-main sequence stars. In its differential form this distribution follows a power-law with index alpha = 2.4 +- 0.5, in the range typically observed on late-type stars and the Sun. The flare energy distribution is probably steep enough to explain the heating of stellar coronae by nano-flares (alpha > 2), albeit associated with a rather large uncertainty that leaves some doubt on this conclusion.


title: Procyon-A and Eta-Bootis: Observational Frequencies Analyzed by the Local-Wave Formalism

abstract: In the present analysis of Procyon-A and Eta-Bootis, we use the local-wave formalism which, despite its lack of precision inherent to any semi-analytical method, uses directly the model profile without any modification when calculating the acoustic mode eigenfrequencies. These two solar-like stars present steep variations toward the center due to the convective core stratification, and toward the surface due to the very thin convective zone. Based on different boundary conditions, the frequencies obtained with this formalism are different from that of the classical numerical calculation. We point out that (1) the frequencies calculated with the local-wave formalism seem to agree better with observational ones. All the frequencies detected with a good confident level including those classified as 'noise' find an identification, (2) some frequencies can be clearly identified here as indications of the core limit.


title: Timing Studies on RXTE Observations of SAX J2103.5+4545

abstract: SAX J2103.5+4545 has been continuously monitored for $\sim $ 900 days by Rossi X-ray Timing Explorer (RXTE) since its outburst in July 2002. Using these observations and previous archival RXTE observations of SAX J2103.5+4545, we refined the binary orbital parameters and find the new orbital period as P= (12.66536 $\pm $ 0.00088) days and the eccentricity as 0.4055$\pm$ 0.0032. With these new orbital parameters, we constructed the pulse frequency and pulse frequency derivative histories of the pulsar and confirmed the correlation between X-ray flux and pulse frequency derivative presented by Baykal, Stark and Swank (2002). We constructed the power spectra for the fluctuations of pulse frequency derivatives and found that the power law index of the noise spectra is 2.13 $\pm$ 0.6. The power law index is consistent with random walk in pulse frequency derivative and is the steepest among the HMXRBs. X-ray spectra analysis confirmed the inverse correlation trend between power-law index and X-ray flux found by Baykal, Stark and Swank (2002).


title: Recovering the star formation rate in the solar neighborhood

abstract: This paper develops a method for obtaining the star formation histories of a mixed, resolved population through the use of color-magnitude diagrams (CMDs). The method provides insight into the local star formation rate, analyzing the observations of the Hipparcos satellite through a comparison with synthetic CMDs computed for different histories with an updated stellar evolution library. Parallax and photometric uncertainties are included explicitly and corrected using the Bayesian Richardson-Lucy algorithm. We first describe our verification studies using artificial data sets. From this sensitivity study, the critical factors determining the success of a recovery for a known star formation rate are a partial knowledge of the IMF and the age-metallicity relation, and sample contamination by clusters and moving groups (special populations whose histories are different than that of the whole sample). Unresolved binaries are less important impediments. We highlight how these limit the method. For the real field sample, complete to Mv < 3.5, we find that the solar neighborhood star formation rate has a characteristic timescale for variation of about 6 Gyr, with a maximum activity close to 3 Gyr ago. The similarity of this finding with column integrated star formation rates may indicate a global origin, possibly a collision with a satellite galaxy. We also discuss applications of this technique to general photometric surveys of other complex systems (e.g. Local Group dwarf galaxies) where the distances are well known.


title: Dust in High-Velocity Clouds : relevance for Planck

abstract: The recent detection of dust emission in Complex C, the largest High-Velocity Cloud (HVC) on the sky, opens a very interesting perspective for Planck. The HVC dust temperature determined using IRAS and Spitzer observations is low (T~10.7 K) in accordance with its great distance from the Galaxy (> 5 kpc). Peak column density in 30 arcmin beam is NH ~10^20 cm-2 which is typical of HVCs and similar to cirrus column density in diffuse regions. On the other hand HVCs appear to be very clumpy at smaller angular scales; several observations at the arcminute scale resolution show significant structure and higher brightness contrasts than in typical cirrus emission. In this contribution we show that, even with their moderate column density, the cold temperature, high emissivity and high column density contrast of HVCs should lead to significant and detectable emission in the Planck-HFI frequency range. In order to separate the HVC emission from the Galactic cirrus emission, the use of 21 cm observations will be mandatory.


title: 0.5 Mpc-scale extended X-ray emission in 4C 23.56

abstract: We present an XMM-Newton observation of the radio galaxy 4C 23.56 at z=2.48 which reveals extended X-ray emission coincident with the radio lobes spanning \~0.5 Mpc. These are the largest X-ray-bright lobes known at z>2. Under the assumption that these features result from inverse-Compton scattering of cosmic microwave background photons by relativistic electrons in the radio source lobes, the measured luminosity of L_0.5-8 keV=7.5x10^44 erg s^-1 implies a minimum energy stored within the lobes of ~10^59 erg. We discuss the potential of the large-scale radio/X-ray lobes to influence evolution of the host galaxy and proto-cluster environment at high redshift.


title: Extremely compact massive galaxies at z~1.4

abstract: The optical rest-frame sizes of 10 of the most massive (~5x10^{11}h_{70}^{-2}M_sun) galaxies found in the near-infrared MUNICS survey at 1.2<z<1.7 are analysed. Sizes were estimated both in the J and K' filters. These massive galaxies are at least a factor of 4_{-1.0}^{+1.9} (+-1 sigma) smaller in the rest-frame V-band than local counterparts of the same stellar mass. Consequently, the stellar mass density of these objects is (at least) 60 times larger than massive ellipticals today. Although the stellar populations of these objects are passively fading, their structural properties are rapidly changing since that redshift. This observational fact disagrees with a scenario where the more massive and passive galaxies are fully assembled at z~1.4 (i.e. a monolithic scenario) and points towards a dry merger scenario as the responsible mechanism for the subsequent evolution of these galaxies.


title: Halting Type I planet migration in non-isothermal disks

abstract: Aims: We investigate the effect of including a proper energy balance on the interaction of a low-mass planet with a protoplanetary disk. Methods: We use a three-dimensional version of the RODEO method to perform hydrodynamical simulations including the energy equation. Radiation is included in the flux-limited diffusion approach. Results: The sign of the torque depends sensitively on the ability of the disk to radiate away the energy generated in the immediate surroundings of the planet. In the case of high opacity, corresponding to the dense inner regions of protoplanetary disks, migration is directed \emph{outward}, instead of the usual inward migration that was found in locally isothermal disks. For low values of the opacity we recover inward migration, and we show that torques originating in the coorbital region are responsible for the change in migration direction.


title: Are isolated planetary-mass objects really isolated? A brown dwarf-exoplanet system candidate in the sigma Orionis cluster

abstract: The recent detection by direct imaging of three giant planets at wide separation (50-250 AU) from their primaries has raised the question about the ``true isolation'' of planetary-mass objects in clusters. Our goal was to test the possibility that some free-floating planetary-mass object could in fact be part of wide planetary systems. We searched in the literature for stellar and brown-dwarf member candidates of the sigma Orionis cluster at small angular separations from published candidate planetary-mass objects. We found one candidate planetary system, SE 70, composed of an X-ray source and a planetary-mass object, namely S Ori 68, separated by only 4.6 arcsec. In order to assess the cluster membership of the X-ray source, we obtained mid-resolution optical spectroscopy using ISIS on the William Herschel Telescope. We also compiled additional data on the target from available astronomical catalogues. We have found that SE 70 follows the spectrophotometric sequence of the cluster and displays spectroscopic features of youth, such as lithium in absorption and chromospheric Halpha emission. The radial velocity is consistent with cluster membership. Hence, SE 70 is very probably a sigma Orionis cluster member. The projected physical separation between SE 70 and S Ori 68 is 1700+/-300 AU at the distance of the cluster. If the common proper motion is confirmed in the near future, the system would be composed of an M5-6 brown dwarf with an estimated mass of ~45 M_Jup and an L5+/-2 giant planet with an estimated mass of ~5 M_Jup. It would be the widest and one of the lowest-mass planetary systems known so far.


title: Simulating planet migration in globally evolving disks

abstract: Numerical simulations of planet-disk interactions are usually performed with hydro-codes that -- because they consider only an annulus of the disk, over a 2D grid -- can not take into account the global evolution of the disk. However, the latter governs planetary migration of type II, so that the accuracy of the planetary evolution can be questioned. To develop an algorithm that models the local planet-disk interactions together with the global viscous evolution of the disk, we surround the usual 2D grid with a 1D grid ranging over the real extension of the disk. The 1D and 2D grids are coupled at their common boundaries via ghost rings, paying particular attention to the fluxes at the interface, especially the flux of angular momentum carried by waves. The computation is done in the frame centered on the center of mass to ensure angular momentum conservation. The global evolution of the disk and the local planet-disk interactions are both well described and the feedback of one on the other can be studied with this algorithm, for a negligible additional computing cost with respect to usual algorithms.


title: The Globular Cluster Kinematics and Galaxy Dark Matter Content of NGC 4649 (M60)

abstract: From observations with the GMOS multi-slit spectrograph on the Gemini North telescope, we have obtained spectra for 39 globular cluster candidates in the Virgo giant elliptical galaxy NGC 4649 (M60), of which 38 are confirmed globular clusters. The clusters extend out to a radius of 260" (3.5 effective radii). We find no rotation of the globular cluster system, with an upper limit of v/sigma < 0.6 at a confidence level of 95%. The globular cluster velocity dispersion is constant with radius, within the uncertainties. We fit isotropic models to the globular cluster and stellar kinematics; these models yield a M/L_V around 16 at 200" radius (16 kpc), an increase of a factor of two from the central M/L. We also use the mass profile as derived from X-rays to determine the orbital structure. Using axisymmetric orbit-based models and the X-ray mass profile, we find the orbital distribution is close to isotropic within 100", and becomes tangentially biased beyond. Furthermore, when using the X-ray profile, we find a better fit to the kinematics compared to using a constant M/L model. Thus, both isotropic and axisymmetric orbit-based models give support for the presence of a dark matter halo in NGC 4649.


title: Pure pseudo-C_l estimators for CMB B-modes

abstract: Fast heuristically weighted, or pseudo-C_l, estimators are a frequently used method for estimating power spectra in CMB surveys with large numbers of pixels. Recently, Challinor & Chon showed that the E-B mixing in these estimators can become a dominant contaminant at low noise levels, ultimately limiting the gravity wave signal which can be detected on a finite patch of sky. We define a modified version of the estimators which eliminates E-B mixing and is near-optimal at all noise levels.


title: Pulse profiles of millisecond pulsars and their Fourier amplitudes

abstract: Approximate analytical formulae are derived for the pulse profile produced by small hot spots on a rapidly rotating neutron star. Its Fourier amplitudes and phases are calculated. The proposed formalism takes into account gravitational bending of light, Doppler effect, anisotropy of emission, and time delays. Its accuracy is checked with exact numerical calculations.


title: The Luminosity Function of high-redshift QSOs - A combined analysis of GOODS and SDSS

abstract: Aims: In this work the luminosity function of QSOs is measured in the redshift range 3.5<z<5.2 for the absolute magnitude interval -21<M_{145}<-28. The determination of the faint end of the luminosity function at these redshifts provides important constraints on models of joint evolution of galaxies and AGNs. Methods: We have defined suitable criteria to select faint QSOs in the GOODS fields, checking in detail their effectiveness and completeness. Spectroscopic follow-up of the resulting QSO candidates has been carried out. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. We have used a Monte-Carlo technique to estimate the properties of the luminosity function, checking various parameterizations for its shape and evolution. Results: Models based on Pure Density Evolution show better agreement with observation than models based on Pure Luminosity Evolution. However a different break magnitude with respect to z~2.1 is required at 3.5<z<5.2. Models with a steeper faint end score a higher probability. We do not find any evidence for a flattening of the bright end at redshift z>3.5. Conclusions: The estimated space density evolution of QSOs indicates a suppression of the formation and/or feeding of Supermassive Black Holes at these redshifts. The QSO contribution to the UV background is insufficient to ionize the IGM at 3.5<z<5.2.


title: Dynamics of Barred Galaxies

abstract: Some 30% of disc galaxies have a pronounced central bar feature in the disc plane and many more have weaker features of a similar kind. Kinematic data indicate that the bar constitutes a major non-axisymmetric component of the mass distribution and that the bar pattern tumbles rapidly about the axis normal to the disc plane. The observed motions are consistent with material within the bar streaming along highly elongated orbits aligned with the rotating major axis. A barred galaxy may also contain a spheroidal bulge at its centre, spirals in the outer disc and, less commonly, other features such as a ring or lens. Mild asymmetries in both the light and kinematics are quite common. We review the main problems presented by these complicated dynamical systems and summarize the effort so far made towards their solution, emphasizing results which appear secure. (Truncated)


title: A Search for Nitrogen-Enhanced Metal-Poor Stars

abstract: Theoretical models of very metal-poor intermediate-mass Asymptotic Giant Branch (AGB) stars predict a large overabundance of primary nitrogen. The very metal-poor, carbon-enhanced, s-process-rich stars, which are thought to be the polluted companions of now-extinct AGB stars, provide direct tests of the predictions of these models. Recent studies of the carbon and nitrogen abundances in metal-poor stars have focused on the most carbon-rich stars, leading to a potential selection bias against stars that have been polluted by AGB stars that produced large amounts of nitrogen, and hence have small [C/N] ratios. We call these stars Nitrogen-Enhanced Metal-Poor (NEMP) stars, and define them as having [N/Fe] > +0.5 and [C/N] < -0.5. In this paper, we report on the [C/N] abundances of a sample of 21 carbon-enhanced stars, all but three of which have [C/Fe] < +2.0. If NEMP stars were made as easily as Carbon-Enhanced Metal-Poor (CEMP) stars, then we expected to find between two and seven NEMP stars. Instead, we found no NEMP stars in our sample. Therefore, this observational bias is not an important contributor to the apparent dearth of N-rich stars. Our [C/N] values are in the same range as values reported previously in the literature (-0.5 to +2.0), and all stars are in disagreement with the predicted [C/N] ratios for both low-mass and high-mass AGB stars. We suggest that the decrease in [C/N] from the low-mass AGB models is due to enhanced extra-mixing, while the lack of NEMP stars may be caused by unfavorable mass ratios in binaries or the difficulty of mass transfer in binary systems with large mass ratios.


title: High-resolution observations of water masers in Bok globules

abstract: We present Very Large Array observations at 1.3 cm of several water maser detections obtained by previous single-dish studies of Bok globules in the Clemens & Barvainis (1988; CB) catalog. We report water maser emission in CB 3 (CB3-mm), CB 54 (IRAS 07020-1618), CB 101 (IRAS 17503-0833), and CB 232 (IRAS 21352+4307), and non-detection towards CB 65 (IRAS 16277-2332) and CB 205 (IRAS 19433+2743). These are the first reported interferometric observations of water masers in Bok globules of the CB catalog. We also present single-dish observations of millimeter and centimeter spectral lines towards CB 101 (IRAS 17503-0833) and CB 65 (IRAS 16277-2332). All the maser emission seems to be associated with star forming regions hosting bipolar molecular outflows, except IRAS 17503-0833 in CB 101, which we suggest to be a possible Mira evolved star, and IRAS 16277-2332 in CB 65, of unknown nature. We have used the precise position of the maser emission to derive information about the powering source of the masers. By analyzing the spatio-kinematical distribution of the water masers, we confirm the millimeter source CB 3-mm as the most likely powering source of the CB 3 masers. We propose the near-IR source CB 232 YC1-I as the best candidate for pumping the maser emission observed in CB 232, while in CB 54, we suggest that the pumping source of the masers could be located at the position of an elongated feature observed in near-infrared maps.


title: The MODEST questions: challenges and future directions in stellar cluster research

abstract: We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)


title: Kinematics of planet-host stars and their relation with dynamical streams in the solar neighbourhood

abstract: We present a detailed study on the kinematics of metal-rich stars with and without planets, and their relation with the Hyades, Sirius and Hercules dynamical streams in the solar neighbourhood. Accurate kinematics have been derived for all the stars belonging to the CORALIE planet search survey. We used precise radial velocity measurements and CCF parameters from the CORALIE database, and parallaxes, photometry and proper motions from the HIPPARCOS and Tycho-2 catalogues. The location of stars with planets in the thin or thick discs has been analysed using both kinematic and chemical constraints. We compare the kinematic behaviour of known planet-host stars to the remaining targets belonging to the volume-limited sample, in particular to its metal-rich population. The high average metallicity of the Hyades stream is confirmed. The planet-host targets show a kinematic behaviour similar to that of the metal-rich comparison subsample, rather than to that of the comparison sample as a whole, thus supporting a primordial origin for the metal excess observed in stars with known planetary companions. According to the scenarios proposed as an explanation for the dynamical streams, systems with giant planets could have formed more easily in metal-rich inner Galactic regions and then been brought into the solar neighbourhood by dynamical streams.


title: The Hybrid Activities of the Pierre Auger Observatory

abstract: The Pierre Auger Observatory detects ultra-high energy cosmic rays by implementing two complementary air-shower techniques. The combination of a large ground array and fluorescence detectors, known as the hybrid concept, means that a rich variety of measurements can be made on a single shower, providing much improved information over what is possible with either detector alone. In this paper I describe the hybrid reconstruction approach and the latest hybrid measurements.


title: On the origin of the ultramassive white dwarf GD50

abstract: We argue on the basis of astrometric and spectroscopic data that the ultramassive white dwarf GD50 is associated with the star formation event that created the Pleiades and is possibly a former member of this cluster. Its cooling age (~60Myrs) is consistent with it having evolved essentially as a single star from a progenitor with a mass M>6Msun so we find no need to invoke a white dwarf-white dwarf binary merger scenario to account for its existence. This result may represent the first direct observational evidence that single star evolution can produce white dwarfs with M>1.1Msun, as predicted by some stellar evolutionary theories. On the basis of its tangential velocity we also provisionally identify the ultramassive (M~1.2Msun) white dwarf PG0136+251 as being related to the Pleiades. These findings may help to alleviate the difficulties in reconciling the observed number of hot nearby ultramassive white dwarfs with the smaller number predicted by binary evolution models under the assumption that they are the products of white dwarf mergers.


title: Enhanced mergers of galaxies in low-redshift clusters

abstract: An ensemble cluster has been formed from a data set comprising a complete magnitude-limited sample of 680 giant galaxies (M(B) < -19) in 8 low-redshift clusters, normalised by the velocity dispersions and virial radii for the early-type cluster populations. Distinct galaxy populations have been identified, including an infall population. A majority (50-70% or greater) of the infall population are found to be in interacting or merging systems characterised by slow gravitational encounters. The observed enhancement of galaxy-galaxy encounters in the infall population compared to the field can be explained by gravitational shocking. It is shown that disc galaxy mergers in the infall population integrated over the estimated lifetime of the cluster (~ 10 Gyr) can readily account for the present cluster S0 population.


title: Test of Thick Pixelated Orbotech Detectors with and without Steering Grids

abstract: We report here on the optimization of 0.5 cm thick pixelated Orbotech CZT detectors with regards to the best contacting materials and the use of steering grids. We evaluated the performance of different contacting materials. Our study differs from earlier ones in that we investigated the performance of different anode and cathode materials separately. We obtain the best performance with Au cathodes. For different anode materials Ti and In give the best energy resolutions. The detector (2.0 x 2.0 x 0.5 cm, 8 x 8 pixels) shows excellent 59 keV, 122 keV and 662 keV energy resolutions of 1.4 keV, 1.9 keV, and 7.4 keV, respectively. Furthermore, we report on using steering grids to improve on the performance of the pixelated detectors. Previously, the benefit of steering grids had been limited by additional electronic noise associated with currents between the negatively biased steering grids and the anode pixels. We are currently exploring the possibility to isolate the steering grid from the CZT substrates by a thin layer of Al$_2$O$_3$. We performed a series of measurements to determine by how much the isolation layer reduces the grid-pixel currents. Comparing the currents between two Au contacts before and after isolating one of the two contacts from the CZT with a 700 nm thick layer of Al2O3, we measure that the isolation layer reduces the currents by a factor of about 10 at 500 V. We present some results from a detector before and after deposition of an isolated steering grid. The grid indeed improves on the detectors energy resolution and detection efficiency. We show that simulations can be used to model the anode to cathode charge correlation in excellent agreement with the experimental results.


title: A VLT/NACO Survey for Triple and Quadruple Systems among Visual Pre-Main Sequence Binaries

abstract: This paper describes a systematic search for high-order multiplicity among wide visual Pre-Main Sequence (PMS) binaries. We conducted an Adaptive Optics survey of a sample of 58 PMS wide binaries from various star-forming regions, which include 52 T Tauri systems with mostly K- and M-type primaries, with the NIR instrument NACO at the VLT. Of these 52 systems, 7 are found to be triple (2 new) and 7 quadruple (1 new). The new close companions are most likely physically bound based on their probability of chance projection and, for some of them, on their position on a color-color diagram. The corresponding degree of multiplicity among wide binaries (number of triples and quadruples divided by the number of systems) is 26.9 +/- 7.2% in the projected separation range 0.07-12 arcsec, with the largest contribution from the Taurus-Auriga cloud. We also found that this degree of multiplicity is twice in Taurus compared to Ophiuchus and Chamaeleon for which the same number of sources are present in our sample. Considering a restricted sample composed of systems at distance 140-190pc, the degree of multiplicity is 26.8 +/- 8.1%, in the separation range 10/14 AU - 1700/2300 AU (30 binaries, 5 triples, 6 quadruples). The observed frequency agrees with results from previous multiplicity surveys within the uncertainties, although a significant overabundance of quadruple systems compared to triple systems is apparent. Tentatively including the spectroscopic pairs in our restricted sample and comparing the multiplicity fractions to those measured for solar-type main-sequence stars in the solar neighborhood leads to the conclusion that both the ratio of triples to binaries and the ratio of quadruples to triples seems to be in excess among young stars. [...]


title: Gravitational Lensing in Modified Gravity and the Lensing of Merging Clusters without Dark Matter

abstract: Gravitational lensing in a modified gravity (MOG) is derived and shown to describe lensing without postulating dark matter. The recent data for merging clusters identified with the interacting cluster 1E0657-56 is shown to be consistent with a weak lensing construction based on MOG without exotic dark matter.


title: Spectropolarimetry of the Deep Impact target comet 9P/Tempel 1 with HiVIS

abstract: High resolution spectropolarimetry of the Deep Impact target, comet 9P/ Tempel 1, was performed during the impact event on July 4th, 2005 with the HiVIS Spectropolarimeter and the AEOS 3.67m telescope on Haleakala, Maui. We observed atypical polarization spectra that changed significantly in the few hours after the impact. The polarization of scattered light as a function of wavelength is very sensitive to the size and composition (complex refractive index) of the scattering particles as well as the scattering geometry. As opposed to most observations of cometary dust, which show an increase in the linear polarization with the wavelength (at least in the visible domain and for phase angles greater than about 30%, a red polarization spectrum) observations of 9P/Tempel 1 at a phase angle of 41 degrees beginning 8 minutes after impact and centered at 6:30UT showed a polarization of 4% at 650 nm falling to 3% at 950 nm. The next observation, centered an hour later showed a polarization of 7% at 650 nm falling to 2% at 950nm. This corresponds to a spectropolarimetric gradient, or slope, of -0.9% per 1000 Angstroms 40 minutes after impact, decreasing to a slope of -2.3% per 1000 Angstroms an hour and a half after impact. This is an atypical blue polarization slope, which became more blue 1 hour after impact. The polarization values of 4% and 7% at 650nm are typical for comets at this scattering angle, whereas the low polarization of 2% and 3% at 950nm is not. We compare observations of comet 9P/Tempel 1 to that of a typical comet, C/2004 Machholz, at a phase angle of 30 degrees which showed a typical red slope, rising from 2% at 650nm to 3% at 950nm in two different observations (+1.0 and +0.9% per 1000 Angstroms).


title: The Distance to the Virgo Cluster from a Recalibrated Tully-Fisher Relation Based on HST Cepheids and a Demonstrated Teerikorpi Cluster Incompleteness Bias

abstract: The importance of the distance of the Virgo cluster in the ongoing debate on the value of the Hubble constant is reviewed. A new calibration of the Tully-Fisher 21-cm line width-absolute magnitude relation is made using Cepheid distances to 25 galaxies determined in various HST programs and reduced with the new Cepheid P-L relations that vary from galaxy-to galaxy. The calibration is applied to a complete sample of Virgo cluster spirals for the purpose of demonstrating the Teerikorpi cluster incompleteness bias. A diagnostic test is shown that should be useful in identifying the presence of bias in incompletely sampled data for distant clusters. The bias-free TF distance modulus for the Virgo cluster is m - M = 31.67 (D = 21.6 Mpc). A systematic correction of 0.07 mag is made because the cluster members are redder in B I on average than the calibrators at a given line width, giving a final adopted modulus for the Virgo cluster core of 31.60 + 0.09. If we assign a generous range of systematic error of ~ 0.3 mag, the distance D = 20.9 Mpc (m - M = 31.60) has a range from 24.0 Mpc to 18.2 Mpc (m - M between 31.9 and 31.3), and a Hubble constant of Ho = 56 between the limits of 49 and 65 when used with a cosmic expansion velocity of 1175 km s-1 determined by the method of distance ratios of remote clusters to Virgo. This range overlaps our preferred value of Ho = 62 from the HST Cepheid calibration of type Ia supernovae recently determined. The TF modulus of Virgo determined here cannot be reconciled with the recent high value of Ho = 72 from Freedman et al.


title: Could GRB 060614 and its presumed host galaxy be a chance superposition?

abstract: The lack of an observed supernova associated with GRB 060614 appears to require a new paradigm for the formation of (a subset of) long-duration GRBs. This requirement is based on the presumed low redshift of the burst, which was inferred from the spatial coincidence of the afterglow with a z=0.125 galaxy. We explore the possibility that this low-redshift galaxy is a chance superposition along the line of sight to GRB 060614. We examine the galaxy distribution of the field of GRB 060614 and find that the probability of a chance association with a galaxy at least as bright as the putative host is only ~ 0.5 - 1.9%. However, for the current ensemble of ~180 Swift GRBs it is likely that several such coincidences have occurred, and given the "non-standard" nature of GRB 060614 it is not implausible that this is one such occurrence. Thus the conclusion that GRB 060614 requires a revision to the formation paradigm for long-duration GRBs should be approached with caution.


title: Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo

abstract: Traditional Markov Chain Monte Carlo methods suffer from low acceptance rate, slow mixing and low efficiency in high dimensions. Hamiltonian Monte Carlo resolves this issue by avoiding the random walk. Hamiltonian Monte Carlo (HMC) is a Markov chain Monte Carlo (MCMC) technique built upon the basic principle of Hamiltonian mechanics. Hamiltonian dynamics allows the chain to move along trajectories of constant energy, taking large jumps in the parameter space with relatively inexpensive computations. This new technique improves the acceptance rate by a factor of 4 and boosts up the efficiency by at least a factor of D in a D-dimensional parameter space. Therefor shorter chains will be needed for a reliable parameter estimation comparing to a traditional MCMC chain yielding the same performance. Besides that, the HMC is well suited for sampling from non-Gaussian and curved distributions which are very hard to sample from using the traditional MCMC methods. The method is very simple to code and can be easily plugged into standard parameter estimation codes such as CosmoMC. In this paper we demonstrate how the HMC can be efficiently used in cosmological parameter estimation.


title: A High-Resolution Spectrum of the Extremely Metal-Rich Bulge G-Dwarf OGLE-2006-BLG-265

abstract: We present an R=45,000 Keck spectrum of the microlensed Galactic bulge G-dwarf OGLE-2006-BLG-265, which has a high (~60) signal-to-noise ratio despite its short (15 min) exposure time because the source was magnified by A~135. While it is very metal-rich ([Fe/H]=0.56), the higher temperature of this star compared with the luminous red giants usually measured in the bulge gives its spectrum many unblended atomic lines. We measure the abundances of 17 elements, including the first abundances for S and Cu in a bulge star. The [alpha/Fe] ratios are subsolar, while the odd-Z elements are slightly supersolar, trends that are also seen in the more metal-rich stars in the bulge and the local Galactic disk. Because the star is a dwarf, the [O/Fe], [Na/Fe], and [Al/Fe] ratios cannot be attributed to internal mixing, as is sometimes claimed for giants. Similar high-resolution spectra could be obtained for about a dozen bulge dwarf stars per year by means of well-designed target-of-opportunity observations.


title: Separating Dark Physics from Physical Darkness: Minimalist Modified Gravity vs. Dark Energy

abstract: The acceleration of the cosmic expansion may be due to a new component of physical energy density or a modification of physics itself. Mapping the expansion of cosmic scales and the growth of large scale structure in tandem can provide insights to distinguish between the two origins. Using Minimal Modified Gravity (MMG) - a single parameter gravitational growth index formalism to parameterize modified gravity theories - we examine the constraints that cosmological data can place on the nature of the new physics. For next generation measurements combining weak lensing, supernovae distances, and the cosmic microwave background we can extend the reach of physics to allow for fitting gravity simultaneously with the expansion equation of state, diluting the equation of state estimation by less than 25% relative to when general relativity is assumed, and determining the growth index to 8%. For weak lensing we examine the level of understanding needed of quasi- and nonlinear structure formation in modified gravity theories, and the trade off between stronger precision but greater susceptibility to bias as progressively more nonlinear information is used.


title: Proto-Neutron Star Winds with Magnetic Fields and Rotation

abstract: We solve the 1D neutrino-heated non-relativistic MHD wind problem for conditions that range from slowly rotating (spin period P > 10 ms) protoneutron stars (PNSs) with surface field strengths typical of radio pulsars (B < 10^13 G), to "proto-magnetars" with B ~ 10^14-10^15 G in their hypothesized rapidly rotating initial states (P ~ 1 ms). We use the simulations of Bucciantini et al. (2006) to map our monopole results onto a more physical dipole geometry and to estimate the spindown of PNSs when their winds are relativistic. We then quantify the effects of rotation and magnetic fields on the mass loss, energy loss, and r-process nucleosynthesis in PNS winds. We describe the evolution of PNS winds through the Kelvin-Helmholtz cooling epoch, emphasizing the transition between (1) thermal neutrino-driven, (2) non-relativistic magnetically-dominated, and (3) relativistic magnetically-dominated outflows. We find that proto-magnetars with P ~ 1 ms and B > 10^15 G drive relativistic winds with luminosities, energies, and Lorentz factors (magnetization sigma ~ 0.1-1000) consistent with those required to produce long duration gamma-ray bursts and hyper-energetic supernovae (SNe). A significant fraction of the rotational energy may be extracted in only a few seconds, sufficiently rapidly to alter the asymptotic energy of the SN remnant, its morphology, and, potentially, its nucleosynthetic yield. Winds from PNSs with more modest rotation periods (2 - 10 ms) and with magnetar-strength fields produce conditions significantly more favorable for the r-process than winds from slowly rotating PNSs. Lastly, we show that energy and momentum deposition by convectively-excited waves further increase the likelihood of successful r-process in PNS winds.


title: CO line emission from Lyman Break Galaxies. Cosmological simulations and predictions for ALMA

abstract: [Abridged] The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA will change the situation dramatically, resulting in the detection of numerous normal galaxies at high redshifts. Maps and spectra of rotational CO line emission were calculated from a cosmological N-body/hydrodynamical TreeSPH simulation of a z ~ 3 "Lyman break galaxy'' of UV luminosity about one order of magnitude below L*. To simulate a typical observation of our system with ALMA, we imposed characteristic noise, angular, and spectral resolution constraints. The CO line properties predicted by our simulation are in good agreement with the two Lyman break systems detected in CO to date. We find that while supernovae explosions from the ongoing star formation carve out large cavities in the molecular ISM, they do not generate large enough gas outflows to make a substantial imprint on the CO line profile. This implies that for most proto-galaxies - except possibly the most extreme cases - stellar feedback effects do not affect CO as a dynamical mass tracer. Detecting CO in sub-L* galaxies at z ~3 will push ALMA to the limits of its cababilities, and whether a source is detected or not may depend critically on its inclination angle. Both these effects (sensitivity and inclination) will severely impair the ability of ALMA to infer the gas kinematics and dynamical masses using line observations.


title: Probing Cool and Warm Infrared Galaxies using Photometric and Structural Measures

abstract: We have analyzed a sample of nearby cool and warm infrared (IR) galaxies using photometric and structural parameters. The set of measures include far-infrared color ($C = \log_{10}[S_{60\mu m}/S_{100\mu m}]$), total IR luminosity ($L_{TIR}$), radio surface brightness as well as radio, near-infrared, and optical sizes. In a given luminosity range cool and warm galaxies are considered as those sources that are found approximately $1 \sigma$ below and above the mean color in the far-infrared $C - L_{TIR}$ diagram. We find that galaxy radio surface brightness is well correlated with color whereas size is less well correlated with color. Our analysis indicates that IR galaxies that are dominated by cool dust are large, massive spirals that are not strongly interacting or merging and presumably the ones with the least active star formation. Dust in these cool objects is less centrally concentrated than in the more typical luminous and ultra-luminous IR galaxies that are dominated by warm dust. Our study also shows that low luminosity early type unbarred and transitional spirals are responsible for the large scatter in the $C - L_{TIR}$ diagram. Among highly luminous galaxies, late type unbarred spirals are predominately warm, and early type unbarred and barred are systematically cooler. We highlight the significance of $C - L_{TIR}$ diagram in terms of local and high redshifts sub-millimeter galaxies.


title: Cluster-Supercluster Alignments

abstract: We study correlations in spatial orientation between galaxy clusters and their host superclusters using a Hubble Volume N-body realization of a concordance cosmology and an analytic model for tidally-induced alignments. We derive an analytic form for distributions of the alignment angle as functions of halo mass (M), ellipticity (epsilon), distance (r) and velocity (v) and show that the model, after tuning of three parameters, provides a good fit to the numerical results. The parameters indicate a high degree of alignment along anisotropic, collapsed filaments. The degree of alignment increases with M and epsilon while it decreases with r and is independent of v. We note the possibility of using the cluster-supercluster alignment effect as a cosmological probe to constrain the slope of the initial power spectrum.


title: New Constraints on the Variable Equation of State Parameter from X-Ray Gas Mass Fractions and SNe Ia

abstract: Recent measurements are suggesting that we live in a flat Universe and that its present accelerating stage is driven by a dark energy component whose equation of state may evolve in time. Assuming two different parameterizations for the function $\omega(z)$, we constrain their free parameters from a joint analysis involving measurements from X-Ray luminosity of galaxy clusters and SNe type Ia data.


title: Multiple dynamical components in Local Group dwarf spheroidals

abstract: The dwarf spheroidal (dSph) satellites of the Local Group have long been thought to be simple spheroids of stars embedded within extended dark matter halos. Recently, however, evidence for the presence of spatially and kinematically distinct stellar populations has been accumulating. Here, we examine the influence of such components on dynamical models of dwarf galaxies embedded in cold dark matter halos. We begin by constructing a model of Andromeda II, a dSph satellite of M31 which shows evidence for spatially distinct stellar components. We find that the two-component model predicts an overall velocity dispersion profile that remains approximately constant at $\sim 10 - 11$ km s$^{-1}$ out to $\sim 1$ kpc from the center; this is despite wide kinematic and spatial differences between the two individual components. The presence of two components may also help to explain oddities in the velocity dispersion profiles of other dSphs; we show that velocity dispersion profiles which appear to rise from the center outwards before leveling off--such as those of Leo I, Draco, and Fornax--can result from the gradual transition from a dynamically cold, concentrated component to a second, hotter, and more spatially extended one, both in equilibrium within the same dark halo. Dwarf galaxies with two stellar components generally have a leptokurtic line-of-sight velocity distribution which is well described by a double Maxwellian. Interestingly, we find that multiple equilibrium components could also provide a potential alternative origin for ``extra-tidal'' stars (normally ascribed to tidal effects) in situations where corroborating evidence for tides may be lacking.


title: Galactic chemical evolution: Carbon through Zinc

abstract: We calculate the evolution of heavy element abundances from C to Zn in the solar neighborhood adopting our new nucleosynthesis yields. Our yields are calculated for wide ranges of metallicity (Z=0-Z_\odot) and the explosion energy (normal supernovae and hypernovae), based on the light curve and spectra fitting of individual supernovae. The elemental abundance ratios are in good agreement with observations. Among the alpha-elements, O, Mg, Si, S, and Ca show a plateau at [Fe/H] < -1, while Ti is underabundant overall. The observed abundance of Zn ([Zn/Fe] ~ 0) can be explained only by the high energy explosion models, which requires a large contribution of hypernovae. The observed decrease in the odd-Z elements (Na, Al, and Cu) toward low [Fe/H] is reproduced by the metallicity effect on nucleosynthesis. The iron-peak elements (Cr, Mn, Co, and Ni) are consistent with the observed mean values at -2.5 < [Fe/H] < -1$, and the observed trend at the lower metallicity can be explained by the energy effect. We also show the abundance ratios and the metallicity distribution functions of the Galactic bulge, halo, and thick disk. Our results suggest that the formation timescale of the thick disk is ~ 1-3 Gyr.


title: On the role of the current loss in radio pulsar evolution

abstract: The aim of this article is to draw attention to the importance of the electric current loss in the energy output of radio pulsars. We remind that even the losses attributed to the magneto-dipole radiation of a pulsar in vacuum can be written as a result of an Ampere force action of the electric currens flowing over the neutron star surface (Michel, 1991, Beskin et al., 1993). It is this force that is responsible for the transfer of angular momentum of a neutron star to an outgoing magneto-dipole wave. If a pulsar is surrounded by plasma, and there is no longitudinal current in its magnetosphere, there is no energy loss (Beskin et al., 1993, Mestel et al., 1999). It is the longitudinal current closing within the pulsar polar cap that exerts the retardation torque acting on the neutron star. This torque can be determined if the structure of longitudinal current is known. Here we remind of the solution by Beskin, Gurevitch & Istomin (1993) and discuss the validity of such an assumption. The behavior of the recently observed "part-time job" pulsar B1931+24 can be naturally explained within the model of current loss while the magneto-dipole model faces difficulties.


title: Environmental Dependence of Cold Dark Matter Halo Formation

abstract: We use a high-resolution $N$-body simulation to study how the formation of cold dark matter (CDM) halos is affected by their environments, and how such environmental effects produce the age-dependence of halo clustering observed in recent $N$-body simulations. We estimate, for each halo selected at redshift $z=0$, an `initial' mass $M_{\rm i}$ defined to be the mass enclosed by the largest sphere which contains the initial barycenter of the halo particles and within which the mean linear density is equal to the critical value for spherical collapse at $z=0$. For halos of a given final mass, $M_{\rm h}$, the ratio $M_{\rm i}/M_{\rm h}$ has large scatter, and the scatter is larger for halos of lower final masses. Halos that form earlier on average have larger $M_{\rm i}/M_{\rm h}$, and so correspond to higher peaks in the initial density field than their final masses imply. Old halos are more strongly clustered than younger ones of the same mass because their initial masses are larger. The age-dependence of clustering for low-mass halos is entirely due to the difference in the initial/final mass ratio. Low-mass old halos are almost always located in the vicinity of big structures, and their old ages are largely due to the fact that their mass accretions are suppressed by the hot environments produced by the tidal fields of the larger structure. The age-dependence of clustering is weaker for more massive halos because the heating by large-scale tidal fields is less important.


title: Non-Gaussianity of the density distribution in accelerating universes II:N-body simulations

abstract: We explore the possibility of putting constraints on dark energy models with statistical property of large scale structure in the non-linear region. In particular, we investigate the $w$ dependence of non-Gaussianity of the smoothed density distribution generated by the nonlinear dynamics. In order to follow the non-linear evolution of the density fluctuations, we apply N-body simulations based on $P^3 M$ scheme. We show that the relative difference between non-Gaussianity of $w=-0.8$ model and that of $w=-1.0$ model is $0.67 %$ (skewness) and $1.2 %$ (kurtosis) for $R=8h^{-1}$ Mpc. We also calculate the correspondent quantities for $R=2h^{-1}$ Mpc, $3.0 %$ (skewness) and $4.5 %$ (kurtosis), and the difference turn out to be greater, even though non-linearity in this scale is so strong that the complex physical processes about galaxy formation affect the galaxy distribution. From this, we can expect that the difference can be tested by all sky galaxy surveys with the help of mock catalogs created by selection functions, which suggests that non-Gaussianity of the density distribution potentially plays an important role for extracting information on dark energy.


title: XMM search for $X$-ray emission from the black hole candidate MACHO-96-BLG-5

abstract: MACHO-96-BLG-5 is a microlensing event observed towards the bulge of the Galaxy with exceptionally long duration of $\sim 970$ days. The microlensing parallax fit parameters were used to estimate the lens mass $M=6^{+10}_{-3}$ M$_{\odot}$ and its distance $d$ which results to be in the range 0.5 kpc - 2 kpc. The upper limit on the absolute brightness for main-sequence stars of the same mass is less than 1 $L_{\odot}$ so that the lens is a good black hole candidate. If it is so, the black hole would accrete by interstellar medium thereby emitting in the $X$-ray band. Here, the analysis of an {\it XMM} deep observation towards MACHO-96-BLG-5 lens position is reported. Only an upper limit (at 99.8% confidence level) to the $X$-ray flux from the lens position of $9.10\times 10^{-15}$-- $1.45\times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the energy band $0.2-10$ keV has been obtained from that deep observation allowing to constrain the putative black hole accretion parameters.


title: On the X-Ray and Mass Distribution in the Merging Galaxy Cluster 1E 0657-56: Ram Pressure-Stripping in Substructures with an NFW Density Profile

abstract: We investigate the X-ray and mass distribution in the merging galaxy cluster 1E 0657-56. We study head-on collisions of two virialized clusters with an NFW density profile in the $\Lambda$CDM universe using an N-body + hydrodynamical code. A clear off-set of an X-ray peak from a mass peak, which is like what is reported in Clow et al. (2004), is first reproduced in the N-body + hydrodynamical simulations. We estimate the ram pressure-stripping conditions of the substructure in mergers of two NFW dark halos using a simple analytical model. We find that the ram pressure dominates the gravity of the substructure when the smaller cluster's mass is less than approximately one tenth of the larger cluster's mass. The characteristic X-ray and mass structures found in 1E 0657-56 suggest that neither the ram pressure nor the gravitational bound force overwhelms the other and that the mass ratio between the progenitors is near the critical value mentioned above.


title: One-armed oscillations in Be star discs

abstract: In this paper we study the effect of the quadrupole-term in the gravitational potential of a rotationally deformed central Be star on one armed density waves in the circumstellar disc. The aim is to explain the observed long-term violet over red (V/R) intensity variations of the double peaked Balmer emission-lines, not only in cool Be star systems, but also in the hot systems like gamma Cas. We have carried out semi-analytic and numerical studies of low-frequency one armed global oscillations in near Keplerian discs around Be stars. In these we have investigated surface density profiles for the circumstellar disc which have inner narrow low surface density or gap regions, just interior to global maxima close to the rapidly rotating star, as well as the mode inner boundary conditions. Our results indicate that it is not necessary to invoke extra forces such as caused by line absorption from the stellar flux in order to explain the long-term V/R variations in the discs around massive Be stars. When there exists a narrow gap between the star and its circumstellar disc, with the result that the radial velocity perturbation is non-zero at the inner disc boundary, we find oscillation (and V/R) periods in the observed range for plausible magnitudes for the rotational quadrupole term.


title: Self-induced conversion in dense neutrino gases: Pendulum in flavour space

abstract: Neutrino-neutrino interactions can lead to collective flavour conversion effects in supernovae and in the early universe. We demonstrate that the case of "bipolar" oscillations, where a dense gas of neutrinos and antineutrinos in equal numbers completely converts from one flavour to another even if the mixing angle is small, is equivalent to a pendulum in flavour space. Bipolar flavour conversion corresponds to the swinging of the pendulum, which begins in an unstable upright position (the initial flavour), and passes through momentarily the vertically downward position (the other flavour) in the course of its motion. The time scale to complete one cycle of oscillation depends logarithmically on the vacuum mixing angle. Likewise, the presence of an ordinary medium can be shown analytically to contribute to a logarithmic increase in the bipolar conversion period. We further find that a more complex (and realistic) system of unequal numbers of neutrinos and antineutrinos is analogous to a spinning top subject to a torque. This analogy easily explains how such a system can oscillate in both the bipolar and the synchronised mode, depending on the neutrino density and the size of the neutrino-antineutrino asymmetry. Our simple model applies strictly only to isotropic neutrino gasses. In more general cases, and especially for neutrinos streaming from a supernova core, different modes couple to each other with unequal strength, an effect that can lead to kinematical decoherence in flavour space rather than collective oscillations. The exact circumstances under which collective oscillations occur in non-isotropic media remain to be understood.


title: X-ray emission from the young brown dwarfs of the Taurus Molecular Cloud

abstract: The XMM-Newton Extended Survey of the TMC (XEST) is a large program designed to systematically investigate the X-ray properties of young stellar/substellar objects in the TMC. In particular, the area surveyed by 15 XMM-Newton pointings (of which three are archival observations), supplemented with one archival Chandra observation, allows us to study 17 BDs with M spectral types. Half of this sample (9 out of 17 BDs) is detected; 7 BDs are detected here for the first time in X-rays. We observed a flare from one BD. We confirm several previous findings on BD X-ray activity: a log-log relation between X-ray and bolometric luminosity for stars (with L*<10 Lsun) and BDs detected in X-rays; a shallow log-log relation between X-ray fractional luminosity and mass; a log-log relation between X-ray fractional luminosity and effective temperature; a log-log relation between X-ray surface flux and effective temperature. We find no significant log-log correlation between the X-ray fractional luminosity and EW(Halpha). Accreting and nonaccreting BDs have a similar X-ray fractional luminosity. The median X-ray fractional luminosity of nonaccreting BDs is about 4 times lower than the mean saturation value for rapidly rotating low-mass field stars. Our TMC BDs have higher X-ray fractional luminosity than BDs in the Chandra Orion Ultradeep Project. The X-ray fractional luminosity declines from low-mass stars to M-type BDs, and as a sample, the BDs are less efficient X-ray emitters than low-mass stars. We thus conclude that while the BD atmospheres observed here are mostly warm enough to sustain coronal activity, a trend is seen that may indicate its gradual decline due to the drop in photospheric ionization degree (abridged).


title: PAMELA - A Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics

abstract: The PAMELA experiment is a satellite-borne apparatus designed to study charged particles in the cosmic radiation with a particular focus on antiparticles. PAMELA is mounted on the Resurs DK1 satellite that was launched from the Baikonur cosmodrome on June 15th 2006. The PAMELA apparatus comprises a time-of-flight system, a magnetic spectrometer, a silicon-tungsten electromagnetic calorimeter, an anticoincidence system, a shower tail catcher scintillator and a neutron detector. The combination of these devices allows antiparticles to be reliably identified from a large background of other charged particles. This paper reviews the design, space qualification and on-ground performance of PAMELA. The in-orbit performance will be discussed in future publications.


title: Spectral Synthesis of Massive Stars in Clusters

abstract: Stellar clusters are thought to be the simplest stellar systems and the closest observational counterparts to theoretical models for single stellar populations. Progress in our understanding of the atmospheres and evolution of massive stars has led to generally reliable synthesis models. The future release of new evolution models with rotation, however, will require non-trivial updates to previously published synthesis models, in particular for all Wolf-Rayet and red supergiant related quantities. Cluster synthesis work is currently progressing from a purely stellar approach to a more comprehensive stellar+cluster perspective. The photometric evolution of stars and the dynamical evolution of clusters are delicately interwoven. Recent work attempts to combine these seemingly related fields.


title: Assessing The Starburst Contribution to the Gamma-Ray and Neutrino Backgrounds

abstract: If cosmic ray protons interact with gas at roughly the mean density of the interstellar medium in starburst galaxies, then pion decay in starbursts is likely to contribute significantly to the diffuse extra-galactic background in both gamma-rays and high energy neutrinos. We describe the assumptions that lead to this conclusion and clarify the difference between our estimates and those of Stecker (2006). Detection of a single starburst by GLAST would confirm the significant contribution of starburst galaxies to the extra-galactic neutrino and gamma-ray backgrounds.


title: An improved deprojection and PSF-deconvolution technique for galaxy-cluster X-ray surface-brightness profiles

abstract: We have developed a regularisation procedure for the direct deprojection and PSF-deconvolution of X-ray surface brightness profiles of clusters of galaxies. This procedure allows us to obtain accurate density profiles in a straightforward manner from X-ray observations (in particular data from XMM-Newton, where the PSF correction is important), while retaining information about substructure in the gas distribution, in contrast to analytic modelling of the profiles. In addition to describing our procedure, we present here a detailed investigation of the accuracy of the method and its error calculations over a wide range of input profile characteristics and data quality using Monte Carlo simulations. We also make comparisons with gas density profiles obtained from Chandra observations, where the PSF correction is small, and with profiles obtained using analytic modelling, which demonstrate that our procedure is a useful improvement over standard techniques. This type of method will be especially valuable in the ongoing analysis of unbiased and complete samples of X-ray clusters, both local and distant, helping to improve the quality of their results.


title: The Cosmic Battery Revisited

abstract: We reinvestigate the generation and accumulation of magnetic flux in optically thin accretion flows around active gravitating objects. The source of the magnetic field is the azimuthal electric current associated with the Poynting-Robertson drag on the electrons of the accreting plasma. This current generates magnetic field loops which open up because of the differential rotation of the flow. We show through simple numerical simulations that what regulates the generation and accumulation of magnetic flux near the center is the value of the plasma conductivity. Although the conductivity is usually considered to be effectively infinite for the fully ionized plasmas expected near the inner edge of accretion disks, the turbulence of those plasmas may actually render them much less conducting due to the presence of anomalous resistivity. We have discovered that if the resistivity is sufficiently high throughout the turbulent disk while it is suppressed interior to its inner edge, an interesting steady-state process is established: accretion carries and accumulates magnetic flux of one polarity inside the inner edge of the disk, whereas magnetic diffusion releases magnetic flux of the opposite polarity to large distances. In this scenario, magnetic flux of one polarity grows and accumulates at a steady rate in the region inside the inner edge and up to the point of equipartition when it becomes dynamically important. We argue that this inward growth and outward expulsion of oppositely-directed magnetic fields that we propose may account for the ~30 min cyclic variability observed in the galactic microquasar GRS1915+105.


title: Searching for massive pre--stellar cores through observations of N2H+ and N2D+

abstract: We have measured the deuterium fractionation (through the column density ratio N(N2D+)/N(N2H+)) and the CO depletion factor (ratio between expected and observed CO abundance) in a sample of 10 high-mass protostellar candidates, in order to understand whether the earliest evolutionary stages of high-mass stars have chemical characteristics similar to those of low-mass ones. The observations were carried out with the IRAM-30m telescope and the JCMT. We have detected N2D+ emission in 7 of the 10 sources of our sample, and found an average value N(N2D+)/N(N2H+)~0.015. This value is 3 orders of magnitude larger than the interstellar D/H ratio, indicating the presence of cold and dense gas, in which the physical-chemical conditions are similar to those observed in low-mass pre-stellar cores. Also, the integrated CO depletion factors show that in the majority of the sources the expected CO abundances are larger than the observed values, with a median ratio of 3.2. In principle, the cold gas that gives origin to the N2D+ emission can be the remnant of the massive molecular core in which the high-mass (proto-)star was born, not yet heated up by the central object. If so, our results indicate that the chemical properties of the clouds in which high-mass stars are born are similar to their low-mass counterparts. Alternatively, this cold gas can be located into one (or more) starless core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot decide which is the correct scenario.


title: Nonlinear simulations of magnetic Taylor-Couette flow with current-free helical magnetic fields

abstract: The magnetorotational instability (MRI) in cylindrical Taylor-Couette flow with external helical magnetic field is simulated for infinite and finite aspect ratios. We solve the MHD equations in their small Prandtl number limit and confirm with time-dependent nonlinear simulations that the additional toroidal component of the magnetic field reduces the critical Reynolds number from $O(10^6)$ (axial field only) to $O(10^3)$ for liquid metals with their small magnetic Prandtl number. Computing the saturated state we obtain velocity amplitudes which help designing proper experimental setups. Experiments with liquid gallium require axial field $\sim 50$ Gauss and axial current $\sim 4$ kA for the toroidal field. It is sufficient that the vertical velocity $u_z$ of the flow can be measured with a precision of 0.1mm/s. We also show that the endplates enclosing the cylinders do not destroy the traveling wave instability which can be observed as presented in earlier studies. For TC containers without and with endplates the angular momentum transport of the MRI instability is shown as to be outwards.


title: K-band spectroscopy of deeply embedded, young OB stars

abstract: We have obtained high resolution (R = 10,000) K-band spectra of candidate young massive stars deeply embedded in high-mass star-forming regions. These objects were selected from a near-infrared survey of 44 regions of high-mass star-formation (Kaper et al, 2006). In these clusters, 38 OB stars are identified whose K-band spectra are dominated by photospheric emission. In almost all those stars, the K-band spectra are indistinguishable from field stars. However, in some stars the profile of the Bracket gamma line is different (less deep, or absent) from those of the O field stars. One of the explanations of these profiles might be an enhanced mass-loss.


title: The population in the background of open clusters: Tracer of the Norma-Cygnus arm

abstract: We present colour-magnitude diagrams of open clusters, located in the range $112^\circ < l < 252^\circ$, manifesting stellar populations in the background of clusters. Some of the populations are found to be located beyond the Perseus arm and may be the part of Norma-Cygnus (outer) arm. The outer arm seems to be continued from $l\sim120^\circ$ to $l\sim235^\circ$. The background populations follow the downward warp of the Galactic plane around $l\sim240^\circ$.


title: Dark Matter: the Connection with Gamma-Ray Astrophysics

abstract: We review the status of indirect Dark Matter searches, focusing in particular on the connection with gamma-ray Astrophysics, and on the prospects for detection with the upcoming space telescope GLAST and Air Cherenkov Telescopes such as CANGAROO, HESS, MAGIC and VERITAS. After a brief introduction where we review the fundamental motivations for indirect searches, we tackle the question of whether it is possible to obtain strong enough evidence from astrophysical observations, to claim discovery of Dark Matter particles. To this purpose, we discuss some recent conflicting claims that have generated some confusion in the field, and present new strategies that may provide the long-awaited smoking-gun for Dark Matter.


title: High-energy gamma-ray emission from the inner jet of LS I+61 303: the hadronic contribution revisited

abstract: LS I+61 303 has been detected by the Cherenkov telescope MAGIC at very high energies, presenting a variable flux along the orbital motion with a maximum clearly separated from the periastron passage. In the light of the new observational constraints, we revisit the discussion of the production of high-energy gamma rays from particle interactions in the inner jet of this system. The hadronic contribution could represent a major fraction of the TeV emission detected from this source. The spectral energy distribution resulting from p-p interactions is recalculated. Opacity effects introduced by the photon fields of the primary star and the stellar decretion disk are shown to be essential in shaping the high-energy gamma-ray light curve at energies close to 200 GeV. We also present results of Monte Carlo simulations of the electromagnetic cascades developed very close to the periastron passage. We conclude that a hadronic microquasar model for the gamma-ray emission in LS I +61 303 can reproduce the main features of its observed high-energy gamma-ray flux.


title: Ultraviolet Spectroscopy of Comet 9P/Tempel 1 with Alice/Rosetta during the Deep Impact Encounter

abstract: We report on spectroscopic observations of periodic comet 9P/Tempel 1 by the Alice ultraviolet spectrograph on the Rosetta spacecraft in conjunction with NASA's Deep Impact mission. Our objectives were to measure an increase in atomic and molecular emissions produced by the excavation of volatile sub-surface material. We unambiguously detected atomic oxygen emission from the quiescent coma but no enhancement at the 10% (1-sigma) level following the impact. We derive a quiescent water production rate of 9 x 10^27 molecules per second with an estimated uncertainty of 30%. Our upper limits to the volatiles produced by the impact are consistent with other estimates.


title: A multi-wavelength census of star formation activity in the young embedded cluster around Serpens/G3-G6

abstract: Aims. The aim of this paper is to characterise the star formation activity in the poorly studied embedded cluster Serpens/G3-G6, located ~ 45' (3 pc) to the south of the Serpens Cloud Core, and to determine the luminosity and mass functions of its population of Young Stellar Objects (YSOs). Methods. Multi-wavelength broadband photometry was obtained to sample the near and mid-IR spectral energy distributions to separate YSOs from field stars and classify the YSO evolutionary stage. ISOCAM mapping in the two filters LW2 (5-8.5 um) and LW3 (12-18 um) of a 19' x 16' field was combined with JHKs data from 2MASS, Ks data from Arnica/NOT, and L' data from SIRCA/NOT. Continuum emission at 1.3 mm (IRAM) and 3.6 cm (VLA) was mapped to study the cloud structure and the coldest/youngest sources. Deep narrow band imaging at the 2.12 um S(1) line of H2 from NOTCam/NOT was obtained to search for signs of bipolar outflows. Results. We have strong evidence for a stellar population of 31 Class II sources, 5 flat-spectrum sources, 5 Class I sources, and two Class 0 sources. Our method does not sample the Class III sources. The cloud is composed of two main dense clumps aligned along a ridge over ~ 0.5 pc plus a starless core coinciding with absorption features seen in the ISOCAM maps. We find two S-shaped bipolar collimated flows embedded in the NE clump, and propose the two driving sources to be a Class 0 candidate (MMS3) and a double Class I (MMS2). For the Class II population we find a best age of ~ 2 Myr and compatibility with recent Initial Mass Functions (IMFs) by comparing the observed Class II luminosity function (LF), which is complete to 0.08 L_sun, to various model LFs with different star formation scenarios and input IMFs.


title: The Origin of Diffuse X-ray and $\gamma$-ray Emission from the Galactic Center Region: Cosmic Ray Particles

abstract: (Abridged) The inner couple hundred pcs of our Galaxy is characterized by significant amount of synchrotron-emitting gas, which appears to co-exist with a large reservoir of molecular gas. The spatial correlation between fluorescent Fe K-alpha line emission at 6.4 keV and molecular line emission from Galactic center molecular clouds has been explained as reflected X-rays from a past outburst of Sgr A*. Here we present multi-wavelength study of a representative Galactic center cloud Sgr C using Chandra, VLA and FCRAO. We note a correlation between the nonthermal radio filaments in Sgr C and the X-ray features, suggesting that the two are related. We account for this distribution in terms of the impact of the relativistic particles from local (nonthermal filaments) and extended sources with diffuse neutral gas producing both a nonthermal bremsstrahlung X-ray continuum emission, as well as diffuse 6.4 keV line emission. The production rate of Fe Kalpha photons associated with the injection of electrons into a cloud as a function of column density is calculated. The required energy density of low-energy cosmic rays associated with the synchrotron emitting radio filaments or extended features is estimated to be in the range between 20 and $\sim10^3$ eV cm^-3 for Sgr C, Sgr B1, Sgr B2, and ``the 45 and -30 \kms'' clouds. We also generalize this idea to explain the pervasive production of diffuse Kalpha line and TeV emission from the Galactic center molecular clouds. In particular, we suggest that IC scattering of the sub-millimeter radiation from dust by relativistic electrons may contribute substantially to the large-scale diffuse TeV emission observed towards the central regions of the Galaxy.


title: CZT imaging detectors for ProtoEXIST

abstract: We describe the detector development for a balloon-borne wide-field hard X-ray (20 - 600 keV) telescope, ProtoEXIST. ProtoEXIST is a pathfinder for both technology and science of the proposed implementation of the Black Hole Finder Probe, Energetic X-ray Imaging Survey telescope (EXIST). The principal technology challenge is the development of large area, close-tiled modules of imaging CZT detectors (1000 cm2 for ProtoEXIST1). We review the updates of the detector design and package concept for ProtoEXIST1 and report the current development status of the CZT detectors, using calibration results of our basic detector unit - 2 x 2 x 0.5 cm CZT crystals with 2.5 mm pixels (8 x 8 array). The current prototype (Rev1) of our detector crystal unit (DCU) shows ~4.5 keV electronics noise (FWHM), and the radiation measurements show the energy resolution (FWHM) of the units is 4.7 keV (7.9%) at 59.5 keV, 5.6 keV (4.6%) at 122 keV, and 7.6 keV (2.1%) at 356 keV. The new (Rev2) DCU with revised design is expected to improve the resolution by ~30%.


title: Galactic structure studies from BATC survey

abstract: We present an analysis of the photometric parallaxes of stars in 21 BATC fields carried out with the National Astronomical Observatories (NAOC) 60/90 cm Schmidt Telescope in 15 intermediate-band filters from 3000 to 10000 {\AA}. In this study, we have adopted a three-component (thin disk, thick disk and halo) model to analyze star counts information. By calculating the stellar space density as a function of distance from the Galactic plane, we determine that the range of scale height for the thin disk varies from 220 to 320 pc. Although 220 pc seems an extreme value, it is close to the lower limit in the literature. The range of scale height for the thick disk is from 600 to 1100 pc, and the corresponding space number density normalization is 7.0-1.0% of the thin disk. We find that the scale height of the disk may be variable with observed direction, which cannot simply be attributed to statistical errors. Possibly the main reasons can be attributed to the disk (mainly the thick disk) is flared, with a scale height increasing with radius. The structure is consistent with merger origin for the thick disk formation. Adopting a de Vaucouleurs $r^{1/4}$ law halo, we also find that the axis ratio towards the Galactic center is somewhat flatter ($\sim 0.4$), while the shape of the halo in the anticentre and antirotation direction is rounder with $c/a> 0.4$. Our results show that star counts in different lines of sight can be used directly to obtain a rough estimate of the shape of the stellar halo. Our solutions support the Galactic models with a flattened inner halo, possibly it is formed by a merger early in the Galaxy's history.


title: Modeling the Emission Processes in Blazars

abstract: Blazars are the most violent steady/recurrent sources of high-energy gamma-ray emission in the known Universe. They are prominent emitters of electromagnetic radiation throughout the entire electromagnetic spectrum. The observable radiation most likely originates in a relativistic jet oriented at a small angle with respect to the line of sight. This review starts out with a general overview of the phenomenology of blazars, including results from a recent multiwavelength observing campaign on 3C279. Subsequently, issues of modeling broadband spectra will be discussed. Spectral information alone is not sufficient to distinguish between competing models and to constrain essential parameters, in particular related to the primary particle acceleration and radiation mechanisms in the jet. Short-term spectral variability information may help to break such model degeneracies, which will require snap-shot spectral information on intraday time scales, which may soon be achievable for many blazars even in the gamma-ray regime with the upcoming GLAST mission and current advances in Atmospheric Cherenkov Telescope technology. In addition to pure leptonic and hadronic models of gamma-ray emission from blazars, leptonic/hadronic hybrid models are reviewed, and the recently developed hadronic synchrotron mirror model for TeV gamma-ray flares which are not accompanied by simultaneous X-ray flares (``orphan TeV flares'') is revisited.


title: Towards a Solution of the Polycyclic Aromatic Hydrocarbon - Diffuse Interstellar Band Hypothesis

abstract: A novel theoretical method is developed to study the polycyclic aromatic hydrocarbon - diffuse interstellar band (PAH-DIB) hypothesis. In this method, a computer program is used to enumerate all PAH molecules with up to a specific number of fused benzene rings. Fast quantum chemical calculations are then employed to calculate the electronic transition energies, oscillator strengths, and rotational constants of these molecules. An electronic database of all PAHs with up to any specific number of benzene rings can be constructed this way. Comparison of the electronic transition energies, oscillator strengths, and rotational band contours of all PAHs in the database with astronomical spectra allows one to identify possible individual PAH carriers of some of the intense narrow DIBs. Using the current database containing up to 10 benzene rings we have selected 8 closed-shell PAHs as possible carriers of the intense lambda6614 DIB.


title: HST Stellar Standards with 1% Accuracy in Absolute Flux

abstract: Free of any atmospheric contamination, HST provides the best available spectrophotometry from the far-UV to the near-IR for stars as faint as V~16. The HST CALSPEC standard star network is based on 3 standard candles: the hot, pure hydrogen white dwarf (WD) stars G191B2B, GD153, and GD71, which have Hubeny NLTE model flux calculations that require the atomic physics for only one atom. These model flux distributions are normalized to the absolute flux for Vega of 3.46x10^{-9} erg cm^{-2} s^{-1} \AA^{-1} at 5556\AA using precise Landolt V band photometry and the V bandpass function corrected for atmospheric transmission by M. Cohen. The 3 primary WD standards provide absolute flux calibrations for FOS, STIS, and NICMOS spectrophotometry from these instruments on the HST. About 32 stellar spectral energy distributions (SEDs) have been constructed with a primary pedigree from the STIS data, which extends from 1150 \AA for the hot stars to a long wavelength limit of 1\mu m. NICMOS grism spectrophotometry provides an extension to 1.9\mu m in the IR for 17 of the HST standards and longward to 2.5\mu m for a few of the brighter stars. Included among these HST standards are Vega, the Sloan standard BD+17$^{\circ}$4708, 3 bright solar analog candidates, 3 cool stars of type M or later, and 5 hot WDs. In addition, 4 K giants and 4 main sequence A-stars have NICMOS spectrophotometry from 0.8-2.5\mu m. The WD fluxes are compared to their modeled SEDs and demonstrate an internal precision of 1-2%, while the A-stars agree with the Cohen IR fluxes to ~2%. Three solar analog candidate stars differ from the solar spectrum by up to 10% in the region of heavy line blanketing from 3000-4000 \AA and show differences in shape of ~5% in the IR around 1.8\mu m.


title: Using the Galactic Dynamics of M7 Dwarfs to Infer the Evolution of Their Magnetic Activity

abstract: We present a spectroscopic study and dynamical analysis of ~2600 M7 dwarfs. We confirm our previous finding that the fraction of magnetically active stars decreases with vertical distance from the Galactic plane. We also show that the mean luminosity of the H-alpha emission has a small but statistically significant decrease with distance. Using space motions for ~1300 stars and a simple one-dimensional dynamical simulation, we demonstrate that the drop in the activity fraction of M7 dwarfs can be explained by thin disk dynamical heating and a rapid decrease of magnetic activity at a mean stellar age of ~6-7 Gyr.


title: Does the AGN Unified Model Evolve with Redshift? Utilizing the X-ray Background to Predict the Mid-Infrared Emission of AGN

abstract: Deep X-ray surveys by Chandra and XMM-Newton have resolved about 80% of the 2-10 keV cosmic extragalactic X-ray background (CXRB) into point sources, the majority of which are obscured AGN. The obscuration might be connected to processes within the host galaxy, possibly the star-formation rate. Here, we use the results of CXRB synthesis calculations as input to detailed Cloudy simulations in order to predict the evolution of AGN properties at several mid-IR wavelengths. Computations were performed for three different evolutions of the AGN type 2/type 1 ratio between z=0 and 1: where the ratio increased as (1+z)^{0.9}, as (1+z)^{0.3} and one with no redshift evolution. Models were calculated with the inner radius of the absorbing gas and dust at 1 pc or at 10 pc. Comparing the results of the calculations to combined X-ray and Spitzer data of AGN shows that the predicted spectral energy distributions are a good description of average AGNs found in the deep surveys. The existing data indicates that the mid-IR emission from an average AGN is best described by models where the attenuating material is ~10 pc from the central engine. We present the expected Spitzer cumulative number count distributions and the evolution of the total AGN (type 1 + type 2) luminosity function (LF) between z=0 and 1 at rest-frame 8 microns and 30 microns for the three evolutionary scenarios. The mid-IR AGN LF will be an excellent tool to measure the evolution of the covering factor of the gas and dust from z~0 to 1.


title: Slow Roll Reconstruction: Constraints on Inflation from the 3 Year WMAP Dataset

abstract: We study the constraints on the inflationary parameter space derived from the 3 year WMAP dataset using ``slow roll reconstruction'', using the SDSS galaxy power spectrum to gain further leverage where appropriate. This approach inserts the inflationary slow roll parameters directly into a Monte Carlo Markov chain estimate of the cosmological parameters, and uses the inflationary flow hierarchy to compute the parameters' scale-dependence. We work with the first three parameters (epsilon, eta and xi) and pay close attention to the possibility that the 3 year WMAP dataset contains evidence for a ``running'' spectral index, which is dominated by the xi term. Mirroring the WMAP team's analysis we find that the permitted distribution of xi is broad, and centered away from zero. However, when we require that inflationary parameters yield at least 30 additional e-folds of inflation after the largest observable scales leave the horizon, the bounds on xi tighten dramatically. We make use of the absence of an explicit pivot scale in the slow roll reconstruction formalism to determine the dependence of the computed parameter distributions on the pivot. We show that the choice of pivot has a significant effect on the inferred constraints on the inflationary variables, and the spectral index and running derived from them. Finally, we argue that the next round of cosmological data can be expected to place very stringent constraints on the region of parameter space open to single field models of slow roll inflation.


title: RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

abstract: We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.


title: The cosmic star formation rate evolution from z=5 to z=0 from the VIMOS VLT Deep Survey

abstract: We present the evolution of the comoving SFR density in the redshift range 0 < z < 5 using the first epoch data release of the VVDS, that is 11564 spectra selected at I_AB=24 over 2200 arcmin^2 in two fields of view, the VVDS-0226-04 and the VVDS-CDFS-0332-27, and the cosmological parameters (Omega_M, Omega_L, h)=(0.3, 0.7, 0.7). We study the multi-wavelength non dust-corrected luminosity densities at 0 < z < 2 from the rest-frame FUV to the optical passbands, and the rest-frame 1500A luminosity functions and densities at 2.7 < z < 5. They evolve from z=1.2 to 0.05 according to (1+z)^{x} with x = 2.05, 1.94, 1.92, 1.14, 0.73, 0.42, 0.30 in the FUV-1500, NUV-2800, U-3600, B-4400, V-5500, R-6500, and I-7900 passbands, respectively. From z = 1.2 to 0.2 the B-band density for the irregular-like galaxies decreases markedly by a factor 3.5 while it increases by a factor 1.7 for the elliptical-like galaxies. We identify several SFR periods; from z = 5 to 3.4 the FUV-band density increases by at most 0.5dex, from z=3.4 to 1.2 it decreases by 0.08dex, from z=1.2 to 0.05 it declines steadily by 0.6dex. For the most luminous M_AB(1500) < -21 galaxies the FUV-band density drops by 2dex from z = 3.9 to 1.2, and for the intermediate -21 < M_AB(1500) < -20 galaxies it drops by 2dex from z = 0.2 to 0. Comparing with dust corrected surveys, at 0.4 < z < 2 the FUV seems obscured by a constant factor of ~1.8-2 mag, while at z < 0.5 it seems progressively less obscured by up to ~0.9-1 mag when the dust-deficient early-type population is increasingly dominating the B-band density. The VVDS results agree with a downsizing picture where the most luminous sources cease to efficiently produce new stars 12 Gyrs ago (at z~4), while intermediate luminosity sources keep producing stars until 2.5 Gyrs ago (at z~0.2).(abridged)


title: The Electron Temperature Gradient in the Galactic Disk

abstract: We derive the electron temperature gradient in the Galactic disk using a sample of HII regions that spans Galactocentric distances 0--17 kpc. The electron temperature was calculated using high precision radio recombination line and continuum observations for more than 100 HII regions. Nebular Galactocentric distances were calculated in a consistent manner using the radial velocities measured by our radio recombination line survey. The large number of nebulae widely distributed over the Galactic disk together with the uniformity of our data provide a secure estimate of the present electron temperature gradient in the Milky Way. Because metals are the main coolants in the photoionized gas, the electron temperature along the Galactic disk should be directly related to the distribution of heavy elements in the Milky Way. Our best estimate of the electron temperature gradient is derived from a sample of 76 sources for which we have the highest quality data. The present gradient in electron temperature has a minimum at the Galactic Center and rises at a rate of 287 +/- 46 K/kpc. There are no significant variations in the value of the gradient as a function of Galactocentric radius or azimuth. The scatter we find in the HII region electron temperatures at a given Galactocentric radius is not due to observational error, but rather to intrinsic fluctuations in these temperatures which are almost certainly due to fluctuations in the nebular heavy element abundances. Comparing the HII region gradient with the much steeper gradient found for planetary nebulae suggests that the electron temperature gradient evolves with time, becoming flatter as a consequence of the chemical evolution of the Milky Way's disk.


title: Stability Properties of Magnetic Tower Jets

abstract: Stability properties of ``magnetic tower'' jets propagating in the gravitationally stratified background have been examined by performing three-dimensional magnetohydrodynamic simulations. The current-carrying, Poynting flux-dominated magnetic tower jet, which possesses a highly wound helical magnetic field, is subject to the current-driven instability (CDI). We find that, under general physical conditions including small perturbations in the initial background profiles, the propagating magnetic tower jets develop the non-axisymmetric, $m=1$ kink mode of the CDI. The kink mode grows on the local Alfv\'en crossing time scale. In addition, two types of kink modes appear in the system. At the central region where external thermal pressure confinement is strong, only the internal kink mode is excited and will grow. A large distance away from the central region where the external thermal pressure becomes low, the external kink mode is observed. As a result, the exterior of magnetic tower jets will be deformed into a large-scale wiggled structure. We also discuss extensively the different physical processes that contribute to the overall stability properties of the magnetic tower jets. Specifically, when the jet propagates in an initially unperturbed background, we find that they can survive the kink mode beyond the point predicted by the well-known Kruskal-Shafranov (K-S) criterion. The stabilization in this case comes mainly from the dynamical relaxation of magnetic twists during the propagation of magnetic towers; the magnetic pitch is reduced and the corresponding K-S critical wavelength becomes longer as the tower jet proceeds. Furthermore, we show that the pressure-driven and Kelvin-Helmholtz instabilities do not occur in the magnetic tower jets.


title: New Galactic Wolf-Rayet stars, and candidates. An annex to The VIIth Catalogue of Galactic Wolf-Rayet Stars

abstract: This paper gathers, from the literature and private communication, 72 new Galactic Population I Wolf-Rayet stars and 17 candidate WCLd stars, recognized and/or discovered after the publication of The VIIth Catalogue of Galactic Wolf-Rayet Stars. This brings the total number of known Galactic Wolf-Rayet stars to 298, of which 24 (8%) are in open cluster Westerlund 1, and 60 (20%) are in open clusters near the Galactic Center.


title: The optical spectral slope variability of 17 blazars

abstract: Many quasi-simultaneous optical observations of 17 blazars are obtained from previous papers published over the last 19 years in order to investigate the spectral slope variability and understand the radiation mechanism of blazars. The long-period dereddened optical spectral slopes are calculated. We analyse the average spectral slope distribution, which suggests that the spectra of flat spectrum radio quasars (FSRQs) and high energy peaked BL Lac objects (HBLs) are probably deformed by other emission components. The average spectral slopes of low energy peaked BL Lac objects(LBLs), which scatter around 1.5, show a good accordance with the synchrotron self-Compton (SSC) loss-dominated model. We present and discuss the variability between the spectral slope and optical luminosity. The spectra of all HBLs and LBLs get flatter when they turn brighter, while for FSRQs this trend does not exist or may even be reversed. This phenomenon may imply that there is a thermal contribution to the optical spectrum for FSRQs. For the FSRQ 1156+295, there is a hint that the slope gets flatter at both the brightest and faintest states. Our result shows that three subclasses locate in different regions in the pattern of slope variability indicator versus average spectral slope. The relativistic jet mechanism is supported by the significant correlation between the optical Doppler factor and the average spectral slope.


title: The Southern Galactic Plane Survey: Polarized Radio Continuum Observations and Analysis

abstract: The Southern Galactic Plane Survey (SGPS) is a radio survey in the 21 cm H I line and in 1.4 GHz full-polarization continuum, observed with the Australia Telescope Compact Array and the Parkes 64m single dish telescope. The survey spans a Galactic longitude of 253 deg < l < 358 deg and a latitude of |b| < 1 deg at a resolution of 100 arcsec and a sensitivity below 1 mJy/beam. This paper presents interferometer only polarized continuum survey data and describes the data taking, analysis processes and data products. The primary data products are the four Stokes parameters I, Q, U, and V in 25 overlapping fields of 5.5 deg by 2 deg, from which polarized intensity, polarization angle and rotation measure are calculated. We describe the effects of missing short spacings, and discuss the importance of the polarized continuum data in the SGPS for studies of fluctuations and turbulence in the ionized interstellar medium and for studying the strength and structure of the Galactic magnetic field.


title: Implication of Existence of Hybrid stars and Theoretical Expectation of Submillisecond Pulsars

abstract: We derive the bulk viscous damping timescale of hybrid stars, neutron stars with quark matter core. The r-mode instability windows of the stars show that the theoretical results are consistent with the rapid rotation pulsar data, which may give an indication for the existence of quark matter in the interior of neutron stars. Hybrid stars instead of neutron or strange stars may lead to submillisecond pulsars.


title: Comparison between Windowed FFT and Hilbert-Huang Transform for Analyzing Time Series with Poissonian Fluctuations: A Case Study

abstract: Hilbert-Huang Transform (HHT) is a novel data analysis technique for nonlinear and non-stationary data. We present a time-frequency analysis of both simulated light curves and an X-ray burst from the X-ray burster 4U 1702-429 with both the HHT and the Windowed Fast Fourier Transform (WFFT) methods. Our results show that the HHT method has failed in all cases for light curves with Poissonian fluctuations which are typical for all photon counting instruments used in astronomy, whereas the WFFT method can sensitively detect the periodic signals in the presence of Poissonian fluctuations; the only drawback of the WFFT method is that it cannot detect sharp frequency variations accurately.


title: Upper bounds on the low-frequency stochastic gravitational wave background from pulsar timing observations: current limits and future prospects

abstract: Using a statistically rigorous analysis method, we place limits on the existence of an isotropic stochastic gravitational wave background using pulsar timing observations. We consider backgrounds whose characteristic strain spectra may be described as a power-law dependence with frequency. Such backgrounds include an astrophysical background produced by coalescing supermassive black-hole binary systems and cosmological backgrounds due to relic gravitational waves and cosmic strings. Using the best available data, we obtain an upper limit on the energy density per unit logarithmic frequency interval of \Omega^{\rm SMBH}_g(1/8yr) h^2 <= 1.9 x 10^{-8} for an astrophysical background which is five times more stringent than the earlier Kaspi et al. (1994) limit of 1.1 x 10^{-7}. We also provide limits on a background due to relic gravitational waves and cosmic strings of \Omega^{\rm relic}_g(1/8yr) h^2 <= 2.0 x 10^{-8} and \Omega^{\rm cs}_g(1/8yr) h^2 <= 1.9 x 10^{-8} respectively. All of the quoted upper limits correspond to a 0.1% false alarm rate together with a 95% detection rate. We discuss the physical implications of these results and highlight the future possibilities of the Parkes Pulsar Timing Array project. We find that our current results can 1) constrain the merger rate of supermassive binary black hole systems at high red shift, 2) rule out some relationships between the black hole mass and the galactic halo mass, 3) constrain the rate of expansion in the inflationary era and 4) provide an upper bound on the dimensionless tension of a cosmic string background.


title: Large Misalignment between Stellar Bar and Dust Pattern in NGC 3488 Revealed by Spitzer and SDSS

abstract: A large position angle misalignment between the stellar bar and the distribution of dust in the late-type barred spiral NGC 3488 was discovered, using mid-infrared images from the Spitzer Space Telescope and optical images from the Sloan Digital Sky Survey (SDSS). The angle between the bar and dust patterns was measured to be 25+-2deg, larger than most of the misalignments found previously in barred systems based on Ha or HI/CO observations. The stellar bar is bright at optical and 3.6um, while the dust pattern is more prominent in the 8um band but also shows up in the SDSS u and g-band images, suggesting a rich interstellar medium environment harboring ongoing star formation. This angular misalignment is unlikely to have been caused by spontaneous bar formation. We suggest that the stellar bar and the dust pattern may have different formation histories, and that the large misalignment was triggered by a tidal interaction with a small companion. A statistical analysis of a large sample of nearby galaxies with archival Spitzer data indicates that bar structure such as that seen in NGC 3488 is quite rare in the local Universe.