title: Missing Thermal Energy of the Intracluster Medium

abstract: The Sunyaev-Zel'dovich (SZ) effect is a direct probe of thermal energy content of the Universe, induced in the cosmic microwave background (CMB) sky through scattering of CMB photons off hot electrons in the intracluster medium (ICM). We report a 9-sigma detection of the SZ signal in the CMB maps of Wilkinson Microwave Anisotropy Probe (WMAP) 3yr data, through study of a sample of 193 massive galaxy clusters with observed X-ray temperatures greater than 3 keV. For the first time, we make a model-independent measurement of the pressure profile in the outskirts of the ICM, and show that it closely follows the profiles obtained by X-ray observations and numerical simulations. We find that our measurements of the SZ effect would account for only half of the thermal energy of the cluster, if all the cluster baryons were in the hot ICM phase. Our measurements indicate that a significant fraction (35 +/- 8 %) of baryonic mass is missing from the hot ICM, and thus must have cooled to form galaxies, intracluster stars, or an unknown cold phase of the ICM. There does not seem to be enough mass in the form of stars or cold gas in the cluster galaxies or intracluster space, signaling the need for a yet-unknown baryonic component (at 3-sigma level), or otherwise new astrophysical processes in the ICM.


title: Testing neutrino spectra formation in collapsing stars with the diffuse supernova neutrino flux

abstract: I address the question of what can be learned from the observation of the diffuse supernova neutrino flux in the precision phase, at next generation detectors of Megaton scale. An analytical study of the spectrum of the diffuse flux shows that, above realistic detection thresholds of 10 MeV or higher, the spectrum essentially reflects the exponential-times-polynomial structure of the original neutrino spectrum at the emission point. There is only a weak (tens of per cent) dependence on the power \beta describing the growth of the supernova rate with the redshift. Different original neutrino spectra correspond to large differences in the observed spectrum of events at a water Cerenkov detector: for typical supernova rates, the ratio of the numbers of events in the first and second energy bins (of 5 MeV width) varies in the interval 1.5 - 4.3 for pure water (energy threshold 18 MeV) and in the range 1 - 2.5 for water with Gadolinium (10 MeV threshold). In the first case discrimination would be difficult due to the large errors associated with background. With Gadolinium, instead, the reduction of the total error down to 10-20 % level would allow spectral sensitivity, with a dramatic improvement of precision with respect to the SN1987A data. Even in this latter case, for typical neutrino luminosity the dependence on \beta is below sensitivity, so that it can be safely neglected in data analysis.


title: Spitzer Observations of 3C Quasars and Radio Galaxies: Mid-Infrared Properties of Powerful Radio Sources

abstract: We have measured mid-infrared radiation from an orientation-unbiased sample of 3CRR galaxies and quasars at redshifts 0.4 < z < 1.2 with the IRS and MIPS instruments on the Spitzer Space Telescope. Powerful emission (L_24micron > 10^22.4 W/Hz/sr) was detected from all but one of the sources. We fit the Spitzer data as well as other measurements from the literature with synchrotron and dust components. The IRS data provide powerful constraints on the fits. At 15 microns, quasars are typically four times brighter than radio galaxies with the same isotropic radio power. Based on our fits, half of this difference can be attributed to the presence of non-thermal emission in the quasars but not the radio galaxies. The other half is consistent with dust absorption in the radio galaxies but not the quasars. Fitted optical depths are anti-correlated with core dominance, from which we infer an equatorial distribution of dust around the central engine. The median optical depth at 9.7 microns for objects with core-dominance factor R > 10^-2 is approximately 0.4; for objects with R < 10^-2, it is 1.1. We have thus addressed a long-standing question in the unification of FR II quasars and galaxies: quasars are more luminous in the mid-infrared than galaxies because of a combination of Doppler-boosted synchrotron emission in quasars and extinction in galaxies, both orientation-dependent effects.


title: Possible Solutions to the Radius Anomalies of Transiting Giant Planets

abstract: We calculate the theoretical evolution of the radii of all fourteen of the known transiting extrasolar giant planets (EGPs) for a variety of assumptions concerning atmospheric opacity, dense inner core masses, and possible internal power sources. We incorporate the effects of stellar irradiation and customize such effects for each EGP and star. Looking collectively at the family as a whole, we find that there are in fact two radius anomalies to be explained. Not only are the radii of a subset of the known transiting EGPs larger than expected from previous theory, but many of the other objects are smaller than the default theory would allow. We suggest that the larger EGPs can be explained by invoking enhanced atmospheric opacities that naturally retain internal heat. This explanation might obviate the necessity for an extra internal power source. We explain the smaller radii by the presence in perhaps all the known transiting EGPs of dense cores, such as have been inferred for Saturn and Jupiter. Importantly, we derive a rough correlation between the masses of our "best-fit" cores and the stellar metallicity that seems to buttress the core-accretion model of their formation. Though many caveats and uncertainties remain, the resulting comprehensive theory that incorporates enhanced-opacity atmospheres and dense cores is in reasonable accord with all the current structural data for the known transiting giant planets.


title: The Angstrom Project Alert System: real-time detection of extragalactic microlensing

abstract: The Angstrom Project is undertaking an optical survey of stellar microlensing events across the bulge region of the Andromeda Galaxy (M31) using a distributed network of two-meter class telescopes. The Angstrom Project Alert System (APAS) has been developed to identify in real time candidate microlensing and transient events using data from the Liverpool and Faulkes North robotic telescopes. This is the first time that real-time microlensing discovery has been attempted outside of the Milky Way and its satellite galaxies. The APAS is designed to enable follow-up studies of M31 microlensing systems, including searches for gas giant planets in M31. Here we describe the APAS and we present a few example light curves obtained during its commissioning phase which clearly demonstrate its real-time capability to identify microlensing candidates as well as other transient sources.


title: Proper Motions of Dwarf Spheroidal Galaxies from Hubble Space Telescope Imaging. V: Final Measurement for Fornax

abstract: The measured proper motion of Fornax, expressed in the equatorial coordinate system, is $(\mu_{\alpha},\mu_{\delta})=(47.6\pm 4.6,-36.0\pm 4.1)$ mas century$^{-1}$. This proper motion is a weighted mean of four independent measurements for three distinct fields. Each measurement uses a quasi-stellar object as a reference point. Removing the contribution of the motion of the Sun and of the Local Standard of Rest to the measured proper motion produces a Galactic rest-frame proper motion of $(\mu_{\alpha}^{\mbox{\tiny{Grf}}}, \mu_{\delta}^{\mbox{\tiny{Grf}}}) = (24.4\pm 4.6,-14.3\pm 4.1)$ mas century$^{-1}$. The implied space velocity with respect to the Galactic center has a radial component of $V_{r}=-31.8 \pm 1.7$ km s$^{-1}$ and a tangential component of $V_{t}=196 \pm 29$ km s$^{-1}$. Integrating the motion of Fornax in a realistic potential for the Milky Way produces orbital elements. The perigalacticon and apogalacticon are 118 (66, 137) kpc and 152 (144, 242) kpc, respectively, where the values in the parentheses represent the 95% confidence intervals derived from Monte Carlo experiments. The eccentricity of the orbit is 0.13 (0.11, 0.38), and the orbital period is 3.2 (2.5, 4.6) Gyr. The orbit is retrograde and inclined by $101^{\circ}$ ($94^{\circ}$, $107^{\circ}$) to the Galactic plane. Fornax could be a member of a proposed ``stream'' of galaxies and globular clusters, however the membership of another proposed galaxy in the stream, Sculptor, has been previously ruled out. Fornax is in the Kroupa-Theis-Boily plane that contains eleven of the Galactic satellite galaxies, but its orbit will take it out of that plane.


title: The Validity of the Cosmic String Pattern Search with the Cosmic Microwave Background

abstract: We introduce a new technique to detect the discrete temperature steps that cosmic strings might have left in the cosmic microwave background (CMB) anisotropy map. The technique provides a validity test on the pattern search of cosmic strings that could serve as the groundwork for future pattern searches. The detecting power of the technique is only constrained by two unavoidable features of CMB data: (1) the finite pixelization of the sky map and (2) the Gaussian fluctuation from instrumental noise and primordial anisotropy. We set the upper limit on the cosmic string parameter as $G\mu\lesssim 3.7\times 10^{-6}$ at the 95% confidence level (CL) and find that the amplitude of the temperature step has to be greater than $44\mu K$ in order to be detected for the {\it{Wilkinson Microwave Anisotropy Probe (WMAP)}} 3 year data.


title: Object detection in multi-epoch data

abstract: In astronomy multiple images are frequently obtained at the same position of the sky for follow-up co-addition as it helps one go deeper and look for fainter objects. With large scale panchromatic synoptic surveys becoming more common, image co-addition has become even more necessary as new observations start to get compared with co-added fiducial sky in real time. The standard co-addition techniques have included straight averages, variance weighted averages, medians etc. A more sophisticated nonlinear response chi-square method is also used when it is known that the data are background noise limited and the point spread function is homogenized in all channels. A more robust object detection technique capable of detecting faint sources, even those not seen at all epochs which will normally be smoothed out in traditional methods, is described. The analysis at each pixel level is based on a formula similar to Mahalanobis distance. The method does not depend on the point spread function.


title: A study of the Sunyaev-Zel'dovich increment using archival SCUBA data

abstract: In a search for evidence of the short wavelength increment in the Sunyaev-Zel'dovich (SZ) effect, we have analyzed archival galaxy cluster data from the Sub-millimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope, resulting in the most complete pointed survey of clusters at 850 microns to date. SCUBA's 850 microns passband overlaps the peak of the SZ increment. The sample consists of 44 galaxy clusters in the range 0 < z < 1.3. Maps of each of the clusters have been made and sources have been extracted; as an ancillary product we generate the most thorough galaxy cluster point source list yet from SCUBA. Seventeen of these clusters are free of obvious AGN and have data deep enough to provide interesting measurements of the expected SZ signal. Specialized analysis techniques are employed to extract the SZ effect signal from these SCUBA data, including using SCUBA's short wavelength band as an atmospheric monitor and fitting the long wavelength channel to a model of the spatial distribution of each cluster's SZ effect. By explicitly excising the exact cluster centre from our analysis we demonstrate that emission from galaxies within the cluster does not contaminate our measurement. The SZ amplitudes from our measurements are consistently higher than the amplitudes inferred from low frequency measurements of the SZ decrement.


title: A low density of the Extragalactic Background Light revealed by the H.E.S.S. spectra of the BLLac objects 1ES 1101-232 and H 2356-309

abstract: The unexpectedly hard spectra measured by HESS for the BLLacs 1ES 1101-232 and H 2356-309 has allowed an upper limit on the Extragalactic Background Light (EBL) to be derived in the optical/near-infrared range, which is very close to the lower limit given by the resolved galaxy counts. This result seems to exclude a large contribution to the EBL from other sources (e.g. Population III stars) and indicates that the intergalactic space is more transparent to gamma-rays than previously thought. A brief discussion of EBL absorption effects on blazar spectra and further observational tests to check this conclusion are presented, including the selection of new candidates for observations with Cherenkov telescopes.


title: Beryllium abundance in turn-off stars of NGC 6752

abstract: Aims: To measure the beryllium abundance in two TO stars of the Globular Cluster NGC 6752, one oxygen rich and sodium poor, the other presumably oxygen poor and sodium rich. Be abundances in these stars are used to put on firmer grounds the hypothesis of Be as cosmochronometer and to investigate the formation of Globular Clusters. Method:We present near UV spectra with resolution R$\sim 45000$ obtained with the UVES spectrograph on the 8.2m VLT Kueyen telescope, analysed with spectrum synthesis based on plane parallel LTE model atmospheres. Results:Be is detected in the O rich star with log(Be/H)=-12.04 $\pm$0.15, while Be is not detected in the other star for which we obtain the upper limit log(Be/H)$<$-12.2. A large difference in nitrogen abundance (1.6 dex) is found between the two stars. Conclusions:The Be measurement is compatible with what found in field stars with the same [Fe/H] and [O/H]. The 'Be age' of the cluster is found to be 13.3 Gyrs, in excellent agreement with the results from main sequence fitting and stellar evolution. The presence of Be confirms the results previously obtained for the cluster NGC 6397 and supports the hypothesis that Be can be used as a clock for the early formation of the Galaxy. Since only an upper limit is found for the star with low oxygen abundance, we cannot decide between competing scenarios of Globular Cluster formation, but we can exclude that 'polluted' stars are substantially younger than 'unpolluted' ones. We stress that the Be test might be the only measurement capable of distinguishing between these scenarios.


title: SN 2006jc: A Wolf-Rayet Star Exploding in a Dense He-Rich Circumstellar Medium

abstract: We present optical photometry and spectra of the peculiar Type Ib supernova (SN) 2006jc. Strong and relatively narrow He I emission lines indicate the progenitor star exploded inside a dense circumstellar medium (CSM) rich in He. An exceptionally blue apparent continuum persists from our first spectrum obtained 15 days after discovery through our last spectrum ~1 month later. Based on the presence of isolated Fe II emission lines, we interpret the blue "continuum" as blended, perhaps fluorescent, Fe emission. One or two of the reddest He I line profiles in our spectra are double peaked, suggesting that the CSM has an aspherical geometry. The He I lines that are superposed on the blue continuum show P-Cygni profiles, while the redder He I lines do not, implying that the blue continuum also originates from an asymmetric mass distribution. The He-rich CSM, aspherical geometry, and line velocities indicate that the progenitor star was a WNE Wolf-Rayet (WR) star. A recent (2 years before the SN), coincident, luminous outburst similar to those seen in luminous blue variables (LBVs) is the leading candidate for the dense CSM. Such an eruption associated with a WR star has not been seen before, indicating that the progenitor star may have recently transitioned from the LBV phase. We also present unpublished spectral and photometric data on SN 2002ao which, along with SN 1999cq, is very similar to SN 2006jc. We propose that these three objects may represent a new and distinct class of SNe arising from WR progenitors surrounded by a dense CSM.


title: The Broad Line Region of AGN: Kinematics and Physics

abstract: In this paper a discussion of kinematics and physics of the Broad Line Region (BLR) is given. The possible physical conditions in the BLR and problems in determination of the physical parameters (electron temperature and density) are considered. Moreover, one analyses the geometry of the BLR and the probability that (at least) a fraction of the radiation in the Broad Emission Lines (BELs) originates from a relativistic accretion disk.


title: Primordial non-Gaussianity and the CMB bispectrum

abstract: We present a new formalism, together with efficient numerical methods, to directly calculate the CMB bispectrum today from a given primordial bispectrum using the full linear radiation transfer functions. Unlike previous analyses which have assumed simple separable ansatze for the bispectrum, this work applies to a primordial bispectrum of almost arbitrary functional form, for which there may have been both horizon-crossing and superhorizon contributions. We employ adaptive methods on a hierarchical triangular grid and we establish their accuracy by direct comparison with an exact analytic solution, valid on large angular scales. We demonstrate that we can calculate the full CMB bispectrum to greater than 1% precision out to multipoles l<1800 on reasonable computational timescales. We plot the bispectrum for both the superhorizon ('local') and horizon-crossing ('equilateral') asymptotic limits, illustrating its oscillatory nature which is analogous to the CMB power spectrum.


title: Transport of exotic anti-nuclei: II- Antiproton and Antideuteron astrophysical uncertainties

abstract: We use a 1D propagation model to study the dependence of the pbar and dbar exotic fluxes on the transport parameters. The simple analytical solutions allow us i) to clarify the origin of the astrophysical uncertainties, and ii) to compare two models used for {\em signal} predictions, namely the constant and the linear Galactic wind models. We also study how these uncertainties should be reduced using forthcoming nuclear cosmic ray data. We confirm that the degeneracy of the transport parameters for a given propagation model leads to very different fluxes for primary antinuclei (~10^2). However, we show that with forthcoming data, these uncertainties could be greatly reduced (~2). As the precision will increase, the astrophysical uncertainty could then be dominated by our ignorance of the correct spatial dependence for some of the transport parameters: for instance, the constant and the linear wind models do not predict the same amount of exotic pbar at low energy.


title: Hemispherical transport equation: modeling of quasiparallel collisionless shocks

abstract: Using a so-called hemispherical model we derive a general transport equation for cosmic ray and thermal particles scattered in pitch angle by magnetic inhomogeneities in a moving collisionless plasma. The weak scattering through 90 degrees results in isotropic particle distributions in each hemisphere. The consideration is not limited by small anisotropies and by the condition that particle velocities are higher than characteristic flow velocity differences. For high velocities and small anisotropies the standard cosmic ray transport equation is recovered. We apply the equations derived for investigation of injection and acceleration of particles by collisionless shocks.


title: Signal Search and Reconstruction by a Trend Filtering Algorithm

abstract: We present additional tests of our algorithm aimed at filtering out systematics due to data reduction and instrumental imperfections in time series obtained by ensemble photometry. Signal detection efficiency is demonstrated, and a method of decreasing the false alarm probability is presented. Including the recently discovered transiting extrasolar planet HAT-P-1, we show various examples on the signal reconstruction capability of the method.


title: Analytical solutions for energy spectra of electrons accelerated by nonrelativistic shock-waves in shell type supernova remnants

abstract: %context {Recent observations of hard X-rays and very high energy gamma-rays from a number of young shell type supernova remnants indicate the importance of detailed quantitative studies of energy spectra of relativistic electrons formed via diffusive shock acceleration accompanied by intense nonthermal emission through synchrotron radiation and inverse Compton scattering.} %aim {The aim of this work was derivation of exact asymptotic solutions of the kinetic equation which describes the energy distribution of shock-accelerated electrons for an arbitrary energy-dependence of the diffusion coefficient.} %method {The asymptotic solutions at low and very high energy domains coupled with numerical calculations in the intermediate energy range allow analytical presentations of energy spectra of electrons for the entire energy region.} %results {Under the assumption that the energy losses of electrons are dominated by synchrotron cooling, we derived the exact asymptotic spectra of electrons without any restriction on the diffusion coefficient. We also obtained simple analytical approximations which describe, with accuracy better than ten percent, the energy spectra of nonthermal emission of shock-accelerated electrons due to the synchrotron radiation and inverse Compton scattering.} %conclusions {The results can be applied for interpretation of X-ray and gamma-ray observations of shell type supernova remnants, as well as other nonthermal high energy source populations like microquasars and large scale synchrotron jets of active galactic nuclei.


title: Wide field aplanatic two-mirror telescopes for ground-based gamma-ray astronomy

abstract: Paper describes novel wide field Atmospheric Cherenkov Telescopes based on a modified Ritchey-Chretien design. The performance characteristics of aplanatic telescopes are compared with traditional Davies-Cotton reflectors.


title: Luminosity and stellar mass functions of disks and spheroids in the SDSS and the supermassive black hole mass function

abstract: Using the Galactica code of Benson et al., we obtain quantitative measurements of spheroid-to-disk ratios for a sample of 8839 galaxies observed in the Sloan Digital Sky Survey. We carry out extensive tests of this code and of Gim2D, finding that they perform similarly in all respects. From the spheroid and disk luminosities, we construct luminosity and stellar mass functions for each component and estimate the relative luminosity and stellar mass densities of disks and spheroids in the local Universe. Assuming a simple one-to-one mapping between between spheroid mass and the mass of a central supermassive black hole, we provide the most accurate determination so far of the black hole mass function in the local universe. From this, we infer a cosmological mass density of black holes of rho_BH=3.77 +/- 0.97 x 10^5 h M_sun Mpc^-3. We compare our results to predictions from current hierarchical models of galaxy formation and these are found to fare well in predicting the qualitative trends observed. We find that stars in disks contribute 35-51% of the local stellar mass density.


title: Cosmology from Planck

abstract: I briefly review some of the main scientific outputs expected from the upcoming Planck mission. Planck will map the CMB sky with 5' resolution and $\mu$K sensitivity, with minimal foreground contribution and superb control on systematics. It will collect the entire information enclosed in the temperature primary anisotropy signal and will also get a good measurement of the polarized component of the CMB. This will have profound implications on our knowledge of the physics of the early universe and on the determination of cosmological parameters.


title: North-South Asymmetry of Solar Activity during Cycle 23

abstract: In this paper, we have made a statistical analysis of solar H$\alpha$ flares that occurred during the period 1996 to 2005 to investigate their spatial distribution with respect to northern and southern hemispheres of the Sun. The analysis includes a total of 21608 single events. The study shows a significant N$-$S asymmetry which is persistent with the evolution of the solar cycle. The flare activity favors the northern hemisphere in general during the rising and maximum phase of the solar cycle (i.e., in 1997, 1999, and 2000), while the declining phase (i.e., from 2001 to 2005) shows a southern dominance. Further, the monthly N$-$S asymmetry index for flares, sunspot numbers and sunspot areas suggests similar variations for these phenomena with the progress of solar cycle. We also find that in terms of asymmetric behavior of solar flares, cycle 23 seems to act quite differently from cycle 22 but comparably to cycle 21.


title: Reverberation Mapping of High-Luminosity Quasars: First Results

abstract: Reverberation mapping of nearby active galactic nuclei has led to estimates of broad-line-region (BLR) sizes and central-object masses for some 37 objects to date. However, successful reverberation mapping has yet to be performed for quasars of either high luminosity (above L_opt~10^{46} erg/s) or high redshift (z>0.3). Over the past six years, we have carried out, at the Hobby-Eberly Telescope, rest-frame-ultraviolet spectrophotometric monitoring of a sample of six quasars at redshifts z=2.2--3.2, with luminosities of L_opt~10^{46.4}--10^{47.6} erg/s, an order of magnitude greater than those of previously mapped quasars. The six quasars, together with an additional five having similar redshift and luminosity properties, were monitored photometrically at the Wise Observatory during the past decade. All 11 quasars monitored show significant continuum variations of order 10%--70%. This is about a factor of two smaller variability than for lower luminosity quasars monitored over the same rest-frame period. In the six objects which have been spectrophotometrically monitored, significant variability is detected in the CIV1550 broad emission line. In several cases the variations track the continuum variations in the same quasar, with amplitudes comparable to, or even greater than, those of the corresponding continua. In contrast, no significant Ly\alpha variability is detected in any of the four objects in which it was observed. Thus, UV lines may have different variability trends in high-luminosity and low-luminosity AGNs. For one quasar, S5~0836+71 at z=2.172, we measure a tentative delay of 595 days between CIV and UV-continuum variations, corresponding to a rest-frame delay of 188 days and a central black-hole mass of 2.6\times10^9 M_\odot.


title: The Angular Diameter of $\lambda$ Bo\"{o}tis}

abstract: Using the CHARA Array and the Palomar Testbed Interferometer, the chemically peculiar star $\lambda$ Bo\"{o}tis has been spatially resolved. We have measured the limb darkened angular diameter to be $\theta_{LD} = 0.533\pm0.029$ mas, corresponding to a linear radius of $R_{\star} = 1.70 \pm 0.10 R_\odot$. The measured angular diameter yields an effective temperature for $\lambda$ Boo of $T_{eff} = 8887 \pm 242$ K. Based upon literature surface gravity estimates spanning $\log{(g)} = 4.0-4.2$ $[\rm{cm s}^{-\rm{2}}]$, we have derived a stellar mass range of $M_{\star} = 1.1 - 1.7$ $M_\odot$. For a given surface gravity, the linear radius uncertainty contributes approximately $\sigma(M_\star) = 0.1-0.2 M_\odot$ to the total mass uncertainty. The uncertainty in the mass (i.e., the range of derived masses) is primarily a result of the uncertainty in the surface gravity. The upper bound of our derived mass range ($\log(g)=4.2, M_\star = 1.7\pm0.2 M_\odot$) is consistent with 100-300 MYr solar-metallicity evolutionary models. The mid-range of our derived masses ($\log(g)=4.1, M_\star = 1.3\pm0.2 M_\odot$) is consistent with 2-3 GYr metal-poor evolutionary models. A more definitive surface gravity determination is required to determine a more precise mass for $\lambda$ Boo.


title: Multi-wavelength Study of Young Massive Star Clusters in the Interacting Galaxy ARP 24

abstract: We made a multi-wavelength study of young massive star clusters (YSCs) in the interacting galaxy ARP 24, using the optical and ultraviolet images from Hubble Space Telescope (HST), Sloan Digital Sky Survey, and Galaxy Evolution Explorer; the mid-infrared images from Spitzer Space Telescope; and the narrow-band Ha image and optical spectra from the NAOC 2.16m telescope. Based on the HST images, we found that the brightest infrared knot in ARP 24 is associated with a complex of five young massive star clusters, within a region of ~ 0.95" radius (127pc) in size. The ages and masses of the star clusters in this complex and other regions were estimated using HST broadband photometries and the Starburst99 synthesis models. The star clusters in this complex are very young (within ages of ~ 3-5 Myr) and massive (masses of ~ 10^5 Msun). The ionization parameter and metallicity of the complex were estimated using the emission line ratios, and the star formation rates were calculated using monochromatic 24um, FUV, and Ha line luminosities. We speculate that ARP 24 may formed by a retrograde fly-by encounter indicated by its one-armed appearance and fan-like structure, and the formation of the YSCs in this galaxy is triggered by the interaction. The clusters in the YSC complex may formed in a single giant molecular cloud simultaneously. From the ultraviolet to mid-infrared spectral energy distributions, we found that the region of the YSC complex is relatively bluer in optical and has higher 24um dust emission relative to the starlight and 8um emission. This warm infrared color may due to strong UV radiation field or other mechanisms (e.g., shocks) within this region which may destroy the Polycyclic Aromatic Hydrocarbons and enhance the small grain emission at 24um.


title: On the theory of magnetar QPOs

abstract: We consider torsional oscillations of magnetars. This problem features rich dynamics due to the strong interaction between the normal modes of a magnetar's crust and a continuum of Magneto-Hydro-Dynamic (MHD) modes in its fluid core. We study the dynamics using a simple model of a magnetar possessing a uniform magnetic field and a thin spherical crust. Firstly, we show that global torsional modes only exist when one introduces unphysically large dissipative terms into the equations of motion; thus global modes are not helpful for understanding the magnetar Quasi-Periodic Oscillations (QPOs). Secondly, we solve the initial-value problem by simulating the sudden release of an initially strained crust and monitoring the subsequent crustal motion. We find that the crustal torsional modes quickly exchange their energy with the MHD continuum in the core, and decay by several orders of magnitude over the course of about 10 oscillation periods. After the initial rapid decay, the crustal motion is stabilized and several time-varying QPOs are observed. The dynamical spectrum of the simulated crustal motion is in qualitative agreement with that of the x-ray light-curve in the tail of a giant magnetar flare. The asymptotic frequencies of some of the QPOs are associated with the special spectral points--the turning points or edges--of the MHD continuum, and are not related to those of the crust. The observed steady low-frequency QPO at 18 Hz is almost certainly associated with the lowest frequency of the MHD continuum, or its first overtone. We also find that drifting QPOs get amplified when they come near the frequencies of the crustal modes. This explains why some of the observed QPOs have frequencies close to the expected crustal frequencies, and why these QPOs are highly variable with time.


title: Evidence Contrary to the Existing Exo-Planet Migration Concept

abstract: Exo-planet migration is assumed to have occurred to explain close-to-star gas giant exo-planets within the context of the so-called standard model of solar system formation, rather than giving cause to question the validity of that particular model. I present evidence against the concept of planet migration, evidence that is historical, interdisciplinary, and model-independent. First, I demonstrate a flaw in the standard model of solar system formation that would lead to the contradiction of terrestrial planets having insufficiently massive cores. Then, I discuss the evidence that points to the Earth previously having been a Jupiter-like close-to-Sun gas giant and the consequences that arise there from. Observations of close-to-star gas giant exo-planets orbiting stars other than our own Sun, rather than being evidence for planet migration, I submit, are evidence for differing degrees of violence associated with the thermonuclear ignition of their particular stars. As observational resolution improves, one might expect to find other exo-planetary systems, like our own, with their innermost planets stripped of gaseous envelopes. Such systems may be a necessary requirement for the existence of life.


title: BBN, ADTree and MLP Comparison in Separating Quasars from Large Survey Catalogues

abstract: We compare the performance of Bayesian Belief Networks (BBN), Multilayer Perceptron (MLP) networks and Alternating Decision Trees (ADtree) on separating quasars from stars with the database from the 2MASS and FIRST survey catalogs. Having a training sample of sources of known object types, the classifiers are trained to separate quasars from stars. By the statistical properties of the sample, the features important for classification are selected. We compare the classification results with and without feature selection. Experiments show that the results with feature selection are better than those without feature selection. From the high accuracy, it is concluded that these automated methods are robust and effective to classify point sources, moreover they all may be applied for large survey projects (e.g. selecting input catalogs) and for other astronomical issues, such as the parameter measurement of stars and the redshift estimation of galaxies and quasars.


title: Probing for dynamics of dark energy and curvature of universe with latest cosmological observations

abstract: We use the newly released 182 Type Ia supernova data combined with the third-year Wilkinson Microwave Anisotropic Probe data (WMAP3) and large scale structure (LSS) information including SDSS and 2dFGRS to constrain the dark energy equation of state (EoS) as well as the curvature of universe $\Omega_K$. Using the full dataset of Cosmic Microwave Background (CMB) and LSS rather than the shift parameter and linear growth factor, we make a Markov Chain Monte Carlo (MCMC) global fit, while paying particular attention to the dark energy perturbation. Parameterizing the EoS as $w_{DE}(a) = w_{0} + w_{1}(1-a)$, we find the best fit of ($w_0,w_1$) is ($-1.053,0.944$) and for $w_{DE}(a) = w_{0} + w_{1}\sin({3/2}\pi \ln(a))$, the best fit for ($w_0,w_1$) is ($-1.614,-1.046$). We find that a flat universe is a good approximation, namely, $|\Omega_K|>0.06$ has been excluded by 2$\sigma$ yet the inclusion of $\Omega_K$ can affect the measurement of DE parameters owing to their correlation and the present systematic effects of SNIa measurements.


title: Roche lobe effects on the atmospheric loss of "Hot Jupiters"

abstract: Observational evidence of a hydrodynamically evaporating upper atmosphere of HD209458b (Vidal-Madjar et al. 2003; 2004) and recent theoretical studies on evaporation scenarios of ``Hot Jupiters'' in orbits around solar-like stars with the age of the Sun indicate that the upper atmospheres of short-periodic exoplanets experience hydrodynamic blow-off conditions resulting in loss rates of the order of about 10^10 - 10^12 g s^-1 (Lammer et al. 2003; Yelle 2004; Baraffe et al. 2004; Lecavlier des Etangs et al. 2004; Jaritz et al. 2005, Tian et al. 2005; Penz et al. 2007). By studying the effect of the Roche lobe on the atmospheric loss from short-periodic gas giants we found, that the effect of the Roche lobe can enhance the hydrodynamic evaporation from HD209458b by about 2 and from OGLE-TR-56b by about 2.5 times. For similar exoplanets which are closer to their host star than OGLE-TR-56b, the enhancement of the mass loss can be even larger. Moreover, we show that the effect of the Roche lobe raises the possibility that ``Hot Jupiters'' can reach blow-off conditions at temperatures which are less than expected (< 10000 K) due to the stellar X-ray and EUV (XUV) heating.


title: Epicyclic Frequencies and Resonant Phenomena Near Black Holes: The Current Status

abstract: We briefly discuss the theory of epicyclic motion near rotating black holes and its relation to the origin of high frequency Quasi-Periodic Oscillations (QPOs) seen in many cases of accreting black hole binaries. We also point out some new frequency ratios predicted by the theory.


title: Physics Results of the Pierre Auger Observatory

abstract: The Southern Pierre Auger Observatory is now nearing completion and the accumulating data sample has already allowed to extract scientific results of relevance for ultra-high energy cosmic-ray (UHECR) physics. This lecture focuses on the description and outcomes of the studies on the UHECR arrival directions (and in particular the anisotropy searches around the Galactic Center) and composition (with the extraction of a limit on the flux of UHE photons and its implications). Perspectives concerning the study of inclined showers and the potential of detection of UHE neutrinos are also briefly discussed.


title: Bright and Dark Matter in Elliptical Galaxies: Mass and Velocity Distributions from Self-consistent Hydrodynamical Simulations

abstract: [Abridged] We have analysed the mass and velocity distributions of two samples of relaxed elliptical-like-objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. Our analysis shows that they are embedded in extended, massive dark matter haloes, and they also have an extended corona of hot diffuse gas. Dark matter haloes have experienced adiabatic contraction along their assembly process. The relative ELO dark- to bright-mass content and space distributions show broken homology, and they are consistent with observational results on the dark matter fraction at the central regions, as well as on the gradients of the mass-to-light ratio profiles for boxy ellipticals, as a function of their stellar masses. These results indicate that massive ellipticals miss stars at their central regions, as compared with less massive ones. Our simulations indicate that these missing baryons could be found beyond the virial radii as a hot, diffuse plasma. The projected stellar mass profiles of our ELOs can be well fit by the S\'ersic (1968) formula. Parameters of these fits reproduce observational correlations of the shape parameter with size, luminosity and velocity dispersion. The total mass density profiles show a power-law behaviour over a large r/r_{vir} interval, consistent with data on massive lens ellipticals at shorter radii. The velocity dispersion profiles show kinematical segregation, with no systematic mass dependence and a positive anisotropy, roughly independent of the radial distance outside the central regions. The consistency with observations strongly suggests that they could also describe important intrinsic characteristics of real ellipticals, as well as some of their properties recently inferred from observational data.


title: Anisotropy of dark matter annihilation with respect to the Galactic plane

abstract: We describe the anisotropy of dark matter clump distribution caused by tidal destruction of clumps in the Galactic disk. A tidal destruction of clumps with orbit planes near the disk plane occurs more efficiently as compared with destruction of clumps at near-polar orbits. A corresponding annihilation of dark matter particles in small-scale clumps produces the anisotropic gamma-ray signal with respect to the Galactic disk. This anisotropy is rather small, 9%, and superimposed on that due to off-centering position of the Sun in the Galaxy. The anisotropy of annihilation signal with respect to the Galactic disk provides a possibility to discriminate dark matter annihilation from the diffuse gamma-ray backgrounds of other origin.


title: Circumstellar Interaction Around Type Ib/c Supernovae and the GRB Connection

abstract: Radio observations of Type Ib/c supernovae suggest that circumstellar interaction takes place with a wide range of wind densities, comparable to that seen in Galactic Wolf-Rayet stars. Efficient production of magnetic field in the shocked region is needed. The X-ray emission observed from some Type Ib/c supernovae is higher than would be expected by the thermal or inverse Compton mechanisms; a synchrotron interpretation requires a flattening of the electron energy spectrum at high energies, as might occur in a cosmic ray dominated shock wave. The wind density variations that are indicated in two supernovae may be due to a binary companion, although variable mass loss from a single star remains a possibility. Other than the optical supernova radiation, the emission from the nearby SN 2006aj/GRB 060218 appears to be powered by a central engine, while that from SN 1998bw/GRB 980425 could be powered by either a central engine or the outer supernova ejecta.


title: Rotating Braneworld Black Holes

abstract: We present a Kerr-Newman type stationary and axisymmetric solution that describes rotating black holes with a tidal charge in the Randall-Sundrum braneworld. The tidal charge appears as an imprint of nonlocal gravitational effects from the bulk space. We also discuss the physical properties of these black holes and their possible astrophysical appearance.


title: Insights into Galaxy Evolution from Mid-infrared Wavelengths

abstract: In this paper, I have attempted to highlight key results from deep extragalactic surveys at mid-infrared wavelengths. I discuss advances in our understanding of dust enshrouded star-formation and AGN activity at 0<z<3 from IRAS, ISO and Spitzer. The data seem to indicate that about 70% of the co-moving star-formation rate density at 0.5<z<3 is obscured by dust and that AGN, including obscured sources, account for <20% of the co-moving bolometric luminosity density. There is tentative evidence that the mode of star-formation changes as a function of redshift; star-formation at z~2 is preferentially in massive, ultraluminous infrared galaxies (ULIRGs) while z~1 sources are luminous infrared galaxies (LIRGs) which are about 1 mag fainter than ULIRGs in the near-infrared. This evolution of the star-formation mode, is similar to the evolution of the redshift distribution of X-ray sources as a function of X-ray luminosity and would suggest an extension of the downsizing hypothesis to include both AGN and star-forming galaxies. Measuring dust-enshrouded star-formation at z>3 will become possible only with future facilities like ALMA. Currently, the presence of dust can only be assessed in a small fraction of the youngest starbursts at z>5 by looking for redshifted large equivalent width H_alpha emission in broadband filters like the IRAC 4.5 micron passband. H_alpha to UV ratios in these objects are a tracer of dust extinction and measuring this ratio in GOODS galaxies indicate dust in ~20% of star-forming galaxies at z>5. Finally, implications for reionization based on the measured stellar mass density and star-formation rates of galaxies at these redshifts are discussed.


title: Arp102B: An ADAF and a Torus ?

abstract: Arp102B is a nearby radio galaxy which displays the presence of double peaked Balmer emission lines. Sub-arcsec Keck mid-infrared imaging and Spitzer spectroscopy reveal a spatially compact mid-infrared source which displays tentative evidence for variability. The F$_{\nu}\propto\nu^{-1.2}$ spectral energy distribution is suggestive of an advection dominated accretion flow. The absence of dust features over the 5-40 micron range make it unlikely that thermal dust emission dominates the mid-infrared luminosity. We also detect the presence of molecular hydrogen in emission which is asymmetrically redshifted by ~500-1000 km/s from the systemic velocity of the galaxy. Since the forbidden, low ionization lines in this galaxy are at the systemic velocity, we suggest that the molecular hydrogen emission arises from a rotating molecular gas structure surrounding the nuclear black hole at a distance of ~1 pc.


title: Dynamic Scalar Torsion and an Oscillating Universe

abstract: For the Poincare gauge theory of gravity we consider the dynamical scalar torsion mode in a cosmological context. We explore in particular the possibility of using dynamical torsion to explain the current state of the accelerating Universe. With certain suitable sets of chosen parameters, this model can give a (qualitatively) proper description of the current universe without a cosmological constant, and the universe described is oscillating with a period of the Hubble time.


title: Dark Energy Parametrizations and the Curvature of the Universe

abstract: We investigate observational constraints on the curvature of the universe not restricting ourselves to a cosmological constant as dark energy, in particular allowing a dark energy equation of state to evolve with time in several ways. We use type Ia supernovae (SNeIa) data from the latest gold data set which includes 182 SNeIa, along with the CMB shift parameter and the baryon acoustic oscillation peak. We show quantitatively that the constraint on the curvature of the universe depends on dark energy model: some popular parametrizations give constraints closely around the flat universe at 5% level (2 sigma C.L.) whereas some parametrizations allow the universe to be as open as Omega_k ~ 0.2.


title: Strangelets, Nuclearites, Q-balls--A Brief Overview

abstract: Astrophysical bounds on the properties and abundances of primordial quark nuggets and cosmic ray strangelets are reviewed. New experiments to search for cosmic ray strangelets in lunar soil and from the International Space Station are described. Analogies with baryonic and supersymmetric Q-balls are briefly mentioned, as are prospects for strangelets as ultra-high energy cosmic rays.


title: The morphology of minor axis gaseous outflows in edge-on Seyfert galaxies

abstract: (Abridged) Spiral galaxies often have extended outflows that permeate beyond the region of the disk. Such outflows have been seen both in starburst galaxies, actively star forming galaxies and galaxies with an AGN. In the latter galaxies it is unknown whether the large-scale outflows are driven by star formation activity or purely by the active nucleus. The aim of our investigation is to study the frequency of extended minor-axis outflows in edge-on Seyfert galaxies to investigate the role of the AGN, the circumnuclear environment and star formation activity within the disk regions, and their importance for IGM enrichment on large scales. Narrowband imaging in two different ionizational stages (H-alpha and [OIII]) was performed to attempt a discrimination between processes associated with the active nucleus and those connected to star forming activity within the disk. The H-alpha morphology of the Seyfert galaxies is usually complex, but only in three out of 14 galaxies did we find evidence for minor axis disk outflows. At the sensitivity of our observations [OIII] emission is generally detected only in the nuclear region. Overall, our results show that extraplanar emission of similar brightness and extent as in the previously known cases of NGC3079 and NGC4388 is not common in Seyfert galaxies of otherwise similar properties. Comparison with our previous results shows that for nearby edge-on spiral galaxies star formation may be a more powerful mechanism for producing DIG than AGN activity. While in general AGN activity undoubtedly plays some role in driving minor-axis outflows, this probably requires higher AGN luminosities than are encountered in our small distance-limited sample.


title: Spectral Methods for Time-Dependent Studies of Accretion Flows. III. Three-Dimensional, Self-gravitating, Magnetohydrodynamic Disks

abstract: Accretion disks are three-dimensional, turbulent, often self-gravitating, magnetohydrodynamic flows, which can be modeled in detail with numerical simulations. In this paper, we present a new algorithm that is based on a spectral decomposition method to simulate such flows. Because of the high order of the method, we can solve the induction equation in terms of the magnetic potential and, therefore, ensure trivially that the magnetic fields in the numerical solution are divergence free. The spectral method also suffers minimally from numerical dissipation and allows for an easy implementation of models for sub-grid physics. Both properties make our method ideal for studying MHD turbulent flows such as those found in accretion disks around compact objects. We verify our algorithm with a series of standard tests and use it to show the development of MHD turbulnce in a simulation of an accretion disk. Finally, we study the evolution and saturation of the power spectrum of MHD turbulence driven by the magnetorotational instability.


title: Cosmological Parameters 2006

abstract: Recent measurements of the cosmic microwave background radiation (CMB), particularly when combined with other datasets, have revolutionised our knowledge of the values of the basic cosmological parameters. Here we summarize the state of play at the end of 2006, focusing on the combination of CMB measurements with the power spectrum of galaxy clustering. We compare the constraints derived from the extant CMB data circa 2005 and the final 2dFGRS galaxy power spectrum, with the results obtained when the WMAP 1-year data is replaced by the 3-year measurements (hereafter WMAP1 and WMAP3). Remarkably, the picture has changed relatively little with the arrival of WMAP3, though some aspects have been brought into much sharper focus. One notable example of this is the index of primordial scalar fluctuations, n_s. Prior to WMAP3, Sanchez et al. (2006) found that the scale invariant value of n_s = 1 was excluded at the 95% level. With WMAP3, this becomes a 3sigma result, with implications for models of inflation. We find some disagreement between the constraints on certain parameters when the 2dFGRS P(k) is replaced by the SDSS measurement. This suggests that more work is needed to understand the relation between the clustering of different types of galaxies and the linear perturbation theory prediction for the power spectrum of matter fluctuations.


title: Exoplanets and the Rossiter-McLaughlin Effect

abstract: A transiting planet eclipses part of the rotating stellar surface, thereby producing an anomalous Doppler shift of the stellar spectrum. Here I review how this "Rossiter-McLaughlin Effect" can be used to characterize exoplanetary systems. In particular, one can measure the angle on the sky between the orbital axis and the stellar rotation axis. This may help to discriminate among migration theories. Measurements have been made for 4 exoplanets, and in all cases the spin and orbital axes are fairly well-aligned. In the future, the Rossiter-McLaughlin effect may also be important as an alternative means of probing exoplanetary atmospheres, and for confirming the transits of objects identified by the satellite missions Corot and Kepler.


title: Splitting neutrino masses and showering into Sky

abstract: Neutrino masses might be as light as a few time the atmospheric neutrino mass splitting. High Energy ZeV cosmic neutrinos (in Z-Showering model) might hit relic ones at each mass in different resonance energies in our nearby Universe. This non-degenerated density and energy must split UHE Z-boson secondaries (in Z-Burst model) leading to multi injection of UHECR nucleons within future extreme AUGER energy. Secondaries of Z-Burst as neutral gamma, below a few tens EeV are better surviving local GZK cut-off and they might explain recent Hires BL-Lac UHECR correlations at small angles. A different high energy resonance must lead to Glashow's anti-neutrino showers while hitting electrons in matter. In air, Glashow's anti-neutrino showers lead to collimated and directional air-showers offering a new Neutrino Astronomy. At greater energy around PeV, Tau escaping mountains and Earth and decaying in flight are effectively showering in air sky. These Horizontal showering is splitting by geomagnetic field in forked shapes. Such air-showers secondaries release amplified and beamed gamma bursts (like observed TGF), made also by muon and electron pair bundles, with their accompanying rich Cherenkov flashes. Also planet' s largest (Saturn, Jupiter) atmosphere limbs offer an ideal screen for UHE GZK and Z-burst tau neutrino, because their largest sizes. Titan thick atmosphere and small radius are optimal for discovering up-going resonant Glashow resonant showers. Earth detection of Neutrino showering by twin Magic Telescopes on top mountains, or by balloons and satellites arrays facing the limbs are the simplest and cheapest way toward UHE Neutrino Astronomy .


title: Reconstruction of Hessence Dark Energy and the Latest Type Ia Supernovae Gold Dataset

abstract: Recently, many efforts have been made to build dark energy models whose equation-of-state parameter can cross the so-called phantom divide $w_{de}=-1$. One of them is the so-called hessence dark energy model in which the role of dark energy is played by a non-canonical complex scalar field. In this work, we develop a simple method based on Hubble parameter $H(z)$ to reconstruct the hessence dark energy. As examples, we use two familiar parameterizations for $H(z)$ and fit them to the latest 182 type Ia supernovae Gold dataset. In the reconstruction, measurement errors are fully considered.


title: New Analytical Formulae for Optically-Thin Accretion Flows

abstract: In a previous paper, we described new analytic formulae for optically-thick supercritical accretion flows (Watarai 2006, hereafter paper 1). Here we present analytic formulae for optically-thin one-temperature accretion flows including the advection-dominated regime, using the ``semi-iterative'' method described in paper 1. Our analytic formulae have two real solutions. The first solution corresponds to the advection-dominated accretion flow (ADAF), and the second solution corresponds to the radiation-dominated accretion flow described by Shapiro, Lightman, & Eardley (the so-called SLE model). Both solutions are given by a cubic equation for the advection parameter $f$, which is the ratio of the advection cooling rate $Q_{\rm adv}$ to the viscous heating rate $Q_{\rm vis}$, i.e., $f=Q_{\rm adv}/Q_{\rm vis}$. Most previous studies assume that $f$ is constant ($f \sim 1$ for the ADAF). However, it is clear that $f$ should be a function of the physical parameters of the radiative-cooling dominated regime. We found that the ratio $f$ can be written as a function of the radius, mass accretion rate, and viscous parameter $\alpha$. Using this formula, we can estimate the transition radius from the inner optically-thin ADAF to the outer optically-thick standard disk, which can be measured using observations of the quiescent state in black hole X-ray binaries.


title: Infrared Molecular Starburst Fingerprints in Deeply Obscured (U)LIRG Nuclei

abstract: High resolution spectra of the Spitzer Space Telescope show vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense (n>~ 10^7 cm^-3), warm (T_ex = 200-700 K) molecular gas with high column densities of these molecules ranging from a few 10^15 - 10^17 cm^-2. Abundances relative to H2, inferred from the silicate optical depth, range from ~10^-7 to 10^-6 and show no correlation with temperature. Theoretical studies show that the high abundances of both C2H2 and HCN exclude a X-ray dominated region (XDR) associated with the toroid surrounding an AGN as the origin of this dense warm molecular gas. Galactic massive protostars in the so-called Hot Core phase have similar physical characteristics with comparable high abundances of C2H2, HCN, and CO2 in the hot phase. However, the abundances of C2H2 and HCN and the C2H2/CO2 and HCN/CO2 ratios are much higher toward the (U)LIRGs in the cooler (T_ex <= 400 K) phase. We suggest that the warm dense molecular gas revealed by the mid-IR absorption lines is associated with a phase of deeply embedded star formation where the extreme pressures and densities of the nuclear starburst environment have inhibited the expansion of HII regions and the global disruption of the star forming molecular cloud cores, and `trapped' the star formation process in an `extended' Hot Core phase.


title: Multi-parameter estimating photometric redshifts with artificial neural networks

abstract: We calculate photometric redshifts from the Sloan Digital Sky Survey Data Release 2 Galaxy Sample using artificial neural networks (ANNs). Different input patterns based on various parameters (e.g. magnitude, color index, flux information) are explored and their performances for redshift prediction are compared. For ANN technique, any parameter may be easily incorporated as input, but our results indicate that using dereddening magnitude produces photometric redshift accuracies often better than the Petrosian magnitude or model magnitude. Similarly, the model magnitude is also superior to Petrosian magnitude. In addition, ANNs also show better performance when the more effective parameters increase in the training set. Finally, the method is tested on a sample of 79, 346 galaxies from the SDSS DR2. When using 19 parameters based on the dereddening magnitude, the rms error in redshift estimation is sigma(z)=0.020184. The ANN is highly competitive tool when compared with traditional template-fitting methods where a large and representative training set is available.


title: Probing Cold Dark Matter Cusps by Gravitational Lensing

abstract: I elaborate on my prediction that an indirect detection of cold dark matter (CDM) may be possible by observing the gravitational lensing effects of the CDM cusp caustics at cosmological distances. Cusps in the distribution of CDM are plentiful once density perturbations enter the nonlinear regime of structure formation. Caustic ring model of galactic halo formation provides a well defined density profile and geometry near the cusps of the caustic rings. I calculate the gravitational lensing effects of the cusps in this model. As a pointlike background source passes behind a cusp of a cosmological foreground halo, the magnification in its image may be detected by present instruments. Depending on the strength of detected effect and the time scale of brightness change, it may even be possible to discriminate between the CDM candidates: axions and weakly interacting massive particles.


title: Time-dependent models of two-phase accretion discs around black holes

abstract: We present time-dependent simulations of a two-phase accretion flow around a black hole. The accretion flow initially is composed of an optically thick and cool disc close to the midplane, while on top and below the disc there is a hot and optically thin corona. We consider several interaction mechanisms as heating of the disc by the corona and Compton cooling of the corona by the soft photons of the disc. Mass and energy can be exchanged between the disc and the corona due to thermal conduction. For the course of this more exploratory work, we limit ourselves to one particular model for a stellar mass black hole accreting at a low accretion rate. We confirm earlier both theoretical and observational results which show that at low accretion rates the disc close to the black hole cannot survive and is evaporated. Given the framework of this model, we now can follow through this phase of disc evaporation time dependently.


title: Flickering in Black Hole Accretion discs

abstract: We present an extension of the King et al. (2004) model for the flickering of black hole accretion discs by taking proper account for the thermal properties of the disc. First we develop a one-dimensional, vertically averaged, one-zone model for an optically thick accretion disc and study the temporal evolution. This limits the current model to the so-called high/soft state, where the X-Ray spectrum is dominated by a thermal black-body component. Then we couple this disc model to the flickering process as described in King et al. (2004). Thus we consider the evolution of a poloidal magnetic field subject to a magnetic dynamo. By comparing to observations of X-Ray binaries in the high-soft state, we can constrain the strength of the energy density of the poloidal magnetic field to a few percent of the energy density of the intrinsic disc magnetic field.


title: Hard X-Ray Properties of Groups of Galaxies as Observed with ASCA

abstract: X-ray spectra of groups of galaxies, obtained with the GIS instrument onboard ASCA, were investigated for diffuse hard X-rays in excess of the soft thermal emission from their inter-galactic medium (IGM). In total, 18 objects with the IGM temperature of 0.7--1.7 keV were studied, including HCG 62 in particular. Non X-ray backgrounds in the GIS spectra were carefully estimated and subtracted. The IGM emission was represented by up to two temperature thermal models, which was determined in a soft energy band below 2.5 keV mainly by the SIS data. When extrapolated to a higher energy range of 4--8 keV, this thermal model under-predicted the background-subtracted GIS counts in HCG 62 and RGH 80 by > 2 sigma significance, even though the background uncertainties and the IGM modeling errors are carefully accounted. A hard excess could be also present in NGC 1399. The excess was successfully explained by a power-law model with a photon index ~ 2, or a thermal emission with a temperature exceeding ~ 3 keV. In HCG 62, the 2--10 keV luminosity of the excess hard component was found to be 5.5E41 erg/s at 2--10 keV, which is ~ 30% of the thermal IGM luminosity in 0.7--2.5 keV. Non-thermal and thermal interpretations of this excess components are discussed.


title: Probing the Disk-jet Connection of the Radio Galaxy 3C120 Observed with Suzaku

abstract: Broad line radio galaxies (BLRGs) are a rare type of radio-loud AGN, in which the broad optical permitted emission lines have been detected in addition to the extended jet emission. Here we report on deep (40ksec x4) observations of the bright BLRG 3C~120 using Suzaku. The observations were spaced a week apart, and sample a range of continuum fluxes. An excellent broadband spectrum was obtained over two decades of frequency (0.6 to 50 keV) within each 40 ksec exposure. We clearly resolved the iron K emission line complex, finding that it consists of a narrow K_a core (sigma ~ 110 eV or an EW of 60 eV), a 6.9 keV line, and an underlying broad iron line. Our confirmation of the broad line contrasts with the XMM-Newton observation in 2003, where the broad line was not required. The most natural interpretation of the broad line is iron K line emission from a face-on accretion disk which is truncated at ~10 r_g. Above 10 keV, a relatively weak Compton hump was detected (reflection fraction of R ~ 0.6), superposed on the primary X-ray continuum of Gamma ~ 1.75. Thanks to the good photon statistics and low background of the Suzaku data, we clearly confirm the spectral evolution of 3C120, whereby the variability amplitude decreases with increasing energy. More strikingly, we discovered that the variability is caused by a steep power-law component of Gamma ~2.7, possibly related to the non-thermal jet emission. We discuss our findings in the context of similarities and differences between radio-loud/quiet objects.


title: A WENO Algorithm of the Temperature and Ionization Profiles around a Point Source

abstract: We develop a numerical solver for radiative transfer problems based on the weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive flux corrections, in order to solve the temperature and ionization profiles around a point source of photons in the reionization epoch. Algorithms for such simulation must be able to handle the following two features: 1. the sharp profiles of ionization and temperature at the ionizing front (I-front) and the heating front (T-front), and 2. the fraction of neutral hydrogen within the ionized sphere is extremely small due to the stiffness of the rate equations of atom processes. The WENO scheme can properly handle these two features, as it has been shown to have high order of accuracy and good convergence in capturing discontinuities and complicated structures in fluid as well as to be significantly superior over piecewise smooth solutions containing discontinuities. With this algorithm, we show the time-dependence of the preheated shell around a UV photon source. In the first stage the I-front and T-front are coincident, and propagate with almost the speed of light. In later stage, when the frequency spectrum of UV photons is hardened, the speeds of propagation of the ionizing and heating fronts are both significantly less than the speed of light, and the heating front is always beyond the ionizing front. In the spherical shell between the I- and T-fronts, the IGM is heated, while atoms keep almost neutral. The time scale of the preheated shell evolution is dependent on the intensity of the photon source. We also find that the details of the pre-heated shell and the distribution of neutral hydrogen remained in the ionized sphere are actually sensitive to the parameters used. The WENO algorithm can provide stable and robust solutions to study these details.


title: Galactic dynamos supported by magnetic helicity fluxes

abstract: We present a simple semi-analytical model of nonlinear, mean-field galactic dynamos and use it to study the effects of various magnetic helicity fluxes. The dynamo equations are reduced using the `no-$z$' approximation to a nonlinear system of ordinary differential equations in time; we demonstrate that the model reproduces accurately earlier results, including those where nonlinear behaviour is driven by a magnetic helicity flux. We discuss the implications and interplay of two types of magnetic helicity flux, one produced by advection (e.g., due to the galactic fountain or wind) and the other, arising from anisotropy of turbulence as suggested by Vishniac & Cho(2001). We argue that the latter is significant if the galactic differential rotation is strong enough: in our model, for $\Rw\la-10$ in terms of the corresponding turbulent magnetic Reynolds number. We confirm that the intensity of gas outflow from the galactic disc optimal for the dynamo action is close to that expected for normal spiral galaxies. The steady-state strength of the large-scale magnetic field supported by the helicity advection is still weaker than that corresponding to equipartition with the turbulent energy. However, the Vishniac-Cho helicity flux can boost magnetic field further to achieve energy equipartition with turbulence. For stronger outflows that may occur in starburst galaxies, the Vishniac-Cho flux can be essential for the dynamo action. However, this mechanism requires a large-scale magnetic field of at least $\simeq1\mkG$ to be launched, so that it has to be preceded by a conventional dynamo assisted by the advection of magnetic helicity by the fountain or wind.


title: Dynamics of Planetary Systems in Star Clusters

abstract: At least 10-15% of nearby sun-like stars have known Jupiter-mass planets. In contrast, very few planets are found in mature open and globular clusters such as the Hyades and 47 Tuc. We explore here the possibility that this dichotomy is due to the post-formation disruption of planetary systems associated with the stellar encounters in long-lived clusters. One supporting piece of evidence for this scenario is the discovery of freely floating low-mass objects in star forming regions. We use two independent numerical approaches, a hybrid Monte Carlo and a direct $N$-body method, to simulate the impact of the encounters. We show that the results of numerical simulations are in reasonable agreement with analytical determinations in the adiabatic and impulsive limits. They indicate that distant stellar encounters generally do not significantly modify the compact and nearly circular orbits. However, moderately close stellar encounters, which are likely to occur in dense clusters, can excite planets' orbital eccentricity and induce dynamical instability in systems which are closely packed with multiple planets. We compute effective cross sections for the dissolution of planetary systems and show that, for all initial eccentricities, dissolution occurs on time scales which are longer than the dispersion of small stellar associations, but shorter than the age of typical open and globular clusters. Although it is much more difficult to disrupt short-period planets, close encounters can excite modest eccentricity among them, such that subsequent tidal dissipation leads to orbital decay, tidal inflation, and even disruption of the close-in planets.


title: Spectroscopic Binary Orbits from Photoelectric Radial Velocities - Paper 191: HD 17310, HD 70645 and HD 80731

abstract: The three objects have been identified as members of the recently recognized class of Gamma Doradus stars, which exhibit multi-periodic photometric variations that are thought to arise from non-radial pulsation. The particular objects treated here also prove to be spectroscopic binaries, for which we provide reliable orbits. The radial velocities exhibit unusually large residuals, in which some of the photometric periodicities can be traced. Some of the same periodicities are also demonstrated by the observed variations in the line profiles, which are quantified here simply in terms of the line-widths.


title: A retrospective view of Miriad

abstract: Miriad is a radio interferometry data-reduction package, designed for taking raw data through to the image analysis stage. The Miriad project, begun in 1988, is now middle-aged. With the wisdom of hindsight, we review design decisions and some of Miriad's characteristics.


title: Inflating Fat Bubbles in Clusters of Galaxies by Wide Jets

abstract: We conduct two-dimensional hydrodynamical simulations of jets expanding in the intra-cluster medium (ICM). We find that for a fat, i.e. more or less spherical, bubble attached to the center to be formed the jet should have high momentum flux and a large opening angle. Typically, the half opening angle should be >50 degrees, and the large momentum flux requires a jet speed of \~10,000 km/sec. The inflation process involves vortices and local instabilities which mix some ICM with the hot bubble. These results predict that most of the gas inside the bubble has a temperature of 3x10^8<T<3x10^9 K, and that large quantities of the cooling gas in cooling flow clusters are expelled back to the intra-cluster medium, and heated up. The magnetic fields and relativistic electrons that produce the synchrotron radio emission might be formed in the shock wave of the jet.


title: Line profile variations in rapidly oscillating Ap stars: resolution of the enigma

abstract: We have carried out the first survey of the pulsational line profile variability in rapidly oscillating Ap (roAp) stars. We analysed high signal-to-noise time-series observations of ten sharp-lined roAp stars obtained with the high-resolution spectrographs attached to the VLT and CFHT telescopes. We investigated in detail the variations of Pr III, Nd II, Nd III and Tb III lines and discovered a prominent change of the profile variability pattern with height in the atmospheres of all studied roAp stars. In every investigated star profile variability of at least one rare-earth ion is characterized by unusual blue-to-red moving features - a behaviour inexplicable in the framework of the standard oblique pulsator model of slowly rotating roAp stars. Using analysis of the line profile moments and spectrum synthesis calculations, we demonstrate that unusual oscillations in spectral lines of roAp stars arise from the pulsational modulation of line widths. This variation occurs approximately in quadrature with the radial velocity changes, and its amplitude rapidly increases with height in stellar atmosphere. We propose that the line width modulation is a consequence of the periodic expansion and compression of turbulent layers in the upper atmospheres of roAp stars. Thus, the line profile changes observed in slowly rotating magnetic pulsators should be interpreted as a superposition of two types of variability: the usual time-dependent velocity field due to an oblique low-order pulsation mode and an additional line width modulation, synchronized with the changes of stellar radius. Our explanation of the line profile variations of roAp stars solves the long-standing observational puzzle and opens new possibilities for constraining geometric and physical properties of the stellar magnetoacoustic pulsations.


title: A Review of Early-Time Optical Follow-ups with 2-m Robotic Telescopes

abstract: We summarise recent deep, rapid GRB follow-up observations using the RoboNet-1.0 network which comprises three fully-robotic 2-m telescopes, the Liverpool Telescope and the Faulkes Telescopes North and South. Observations begin automatically within minutes of receipt of a GRB alert and may continue for hours or days to provide well-sampled multi-colour light curves or deep upper limits. Our light curves show a variety of early afterglow behaviour, from smooth, simple or broken power laws to 'bumpy', for a wide range of optical brightness (from the unprecedented faint detections of GRB 060108 and GRB 060510B to classical bright ones). We discuss GRB 051111 as an example of how the combination of optical and X-ray light curves can provide insight into the circumburst environment, in particular the role played by intrinsic extinction soon after the burst.


title: Transiting Planet Simulations from the All Sky Extrasolar Planets Survey

abstract: Many of the planets discovered via the radial velocity technique are hot Jupiters in 3-5 day orbits with ~10$% chance of transiting their parent star. However, radial velocity surveys for extra-solar planets generally require substantial amounts of large telescope time in order to monitor a sufficient number of stars due to the single-object capabilities of the spectrograph. A multi-object Doppler survey instrument has been developed which is based on the dispersed fixed-delay interferometer design. We present simulations of the expected results from the Sloan Doppler survey based on calculated noise models and sensitivity for the instrument and the known distribution of exoplanetary system parameters. We have developed code for automatically sifting and fitting the planet candidates produced by the survey to allow for fast follow-up observations to be conducted. A transit ephemeris is automatically calculated by the code for each candidate and updated when new data becomes available. The techniques presented here may be applied to a wide range of multi-object planet surveys.


title: Quantum inflaton, primordial metric perturbations and CMB fluctuations

abstract: We compute the primordial scalar, vector and tensor metric perturbations arising from quantum field inflation. Quantum field inflation takes into account the nonperturbative quantum dynamics of the inflaton consistently coupled to the dynamics of the (classical) cosmological metric. For chaotic inflation, the quantum treatment avoids the unnatural requirements of an initial state with all the energy in the zero mode. For new inflation it allows a consistent treatment of the explosive particle production due to spinodal instabilities. Quantum field inflation (under conditions that are the quantum analog of slow roll) leads, upon evolution, to the formation of a condensate starting a regime of effective classical inflation. We compute the primordial perturbations taking the dominant quantum effects into account. The results for the scalar, vector and tensor primordial perturbations are expressed in terms of the classical inflation results. For a N-component field in a O(N) symmetric model, adiabatic fluctuations dominate while isocurvature or entropy fluctuations are negligible. The results agree with the current WMAP observations and predict corrections to the power spectrum in classical inflation. Such corrections are estimated to be of the order of m^2/H^2 where m is the inflaton mass and H the Hubble constant at horizon crossing. This turns to be about 4% for the cosmologically relevant scales. This quantum field treatment of inflation provides the foundations to the classical inflation and permits to compute quantum corrections to it.


title: Spectroscopy of Optical Counterparts of Ultraluminous X-ray Sources

abstract: Here we present the results of panoramic and long-slit observations of eight ULX nebular counterparts held with the 6m SAO telescope. In two ULXNe we detected for the first time signatures of high excitation ([OIII]5007 / H\beta > 5). Two of the ULXs were identified with young (T ~ 5-10 Myr) massive star clusters. Four of the eight ULX Nebulae (ULXNe) show bright high-excitation lines. This requires existence of luminous (~ 10^{38} .. 10^{40} erg/s) UV/EUV sources coinciding with the X-ray sources. Other 4 ULXNe require shock excitation of the gas with shock velocities of 20-100km/s. However, all the studied ULXN spectra show signatures of shock excitation, but even those ULXNe where the shocks are prevailing show presence of a hard ionizing source with the luminosity at least ~10^{38} erg/s. Most likely shock waves, X-ray and EUV ionization act simultaneously in all the ULXNe, but they may be roughly separated in two groups, shock-dominated and photoionization-dominated ULXNe. The ULXs have to produce strong winds and/or jets powering their nebulae with \~10^{39} erg/s. Both the wind/jet activity and the EUV source needed are consistent with the suggestion that ULXs are high-mass X-ray binaries with the supercritical accretion disks of the SS433 type.


title: The Aptly Named Phoenix Dwarf Galaxy

abstract: The Local Group galaxy Phoenix has the properties of a dwarf spheroidal galaxy, but an adjacent HI cloud has been found to be at the same radial velocity as the stars. The proximity suggests that this cloud is associated with the most recent ($\le 100$ Myr) star formation in Phoenix. We have obtained relatively high sensitivity and high resolution HI imaging with the VLA with the goal of distinguishing between different processes for displacing the gas from the galaxy. Due to the outer curvature of the HI cloud, it appears that expulsion from the galaxy by winds from supernovae is more likely than ram-pressure stripping. The isolation of the galaxy makes tidal stripping highly unlikely. Using a star formation history constructed from HST imaging, we construct a simple kinematic model which implies that the HI cloud is still gravitationally bound to the galaxy. Gas which is expelled from the centers of dwarf galaxies but which remains gravitationally bound may explain the episodic star formation observed in several dwarfs. In the specific case of Phoenix, there may be future star formation in this currently dSph-like galaxy.


title: Heating of Heavy Ions by Interplanetary Coronal Mass Ejection (ICME) Driven Collisionless Shocks

abstract: Shock heating and particle acceleration processes are some of the most fundamental physical phenomena of plasma physics with countless applications in laboratory physics, space physics, and astrophysics. This study is motivated by previous observations of non-thermal heating of heavy ions in astrophysical shocks (Korreck et al. 2004). Here, we focus on shocks driven by Interplanetary Coronal Mass Ejections (ICMEs) which heat the solar wind and accelerate particles. This study focuses specifically on the heating of heavy ions caused by these shocks. Previous studies have focused only on the two dynamically dominant species, H+ and He2+ . This study utilizes thermal properties measured by the Solar Wind Ion Composition Spectrometer (SWICS) aboard the Advanced Composition Explorer (ACE) spacecraft to examine heavy ion heating. This instrument provides data for many heavy ions not previously available for detailed study, such as Oxygen (O6+, O7+), Carbon (C5+, C6+), and Iron (Fe10+). The ion heating is found to depend critically on the upstream plasma


title: Continuous statistics of the Lyman-alpha forest at 0 < z < 1.6: the mean flux, flux distribution, and autocorrelation from HST FOS spectra

abstract: We measure the amount of absorption in the Lyman-alpha forest at 0 < z < 1.6 in HST FOS spectra of 74 QSOs. At 0 < z < 1.6 we find that 79% of the absorption is from the low density intergalactic medium, 12% from metals and 9% from the strong H I lines, nearly identical to the percentages (78, 15 and 7) that we measured independently at z=2 from spectra taken with the Kast spectrograph on the Lick 3-m. At z=1 the low density intergalactic medium absorbs 0.037 +/- 0.004 of the flux. The error includes some but not all of the uncertainty in the continuum level. The remaining part gives relative errors of approximately 0.21 when we report the mean absorption in eight independent redshift intervals, and 0.047 when we average over all redshifts. We find 1.46 times more absorption from the low density intergalactic medium than comes from Ly-alpha lines that Bechtold et al. 2002 listed in the same spectra. The amount of absorption increases with z and can be fit by a power law in (1+z) with index 1.01. This corresponds to no change in the number of lines, of fixed rest frame equivalent widths, per unit redshift, consistent with the Janknecht et al. 2006 results on the distribution of lines. When we include similar measurements from higher redshifts, we need more degrees of freedom to fit the amount of absorption at 0 < z < 3.2. A power law with a break in slope, changing from index 1.5 at low z to 3.0 above z ~ 1.1 is a better but only marginally acceptable fit. We also calculate two other continuous statistics, the flux probability distribution function and the flux autocorrelation function that is non zero out to v ~ 500 km/sec at 0.5 < z < 1.5.


title: Longitude-resolved imaging of Jupiter at lambda = 2 cm

abstract: We present a technique for creating a longitude-resolved image of Jupiter's thermal radio emission. The technique has been applied to VLA data taken on 25 January 1996 at a wavelength of 2 cm. A comparison with infrared data shows a good correlation between radio hot spots and the 5 micron hot spots seen on IRTF images. The brightest spot on the radio image is most likely the hot spot through which the Galileo probe entered Jupiter's atmosphere. We derived the ammonia abundance (= volume mixing ratio) in the hot spot, which is ~ 3 x 10^{-5}, about half that seen in longitude-averaged images of the NEB, or less than 1/3 of the longitude-averaged ammonia abundance in the EZ. This low ammonia abundance probably extends down to at least the 4 bar level.


title: The Role of Ambipolar Diffusion in the Formation Process of Moderately Magnetized Diffuse Clouds

abstract: We investigate the dynamical condensation process in a magnetized thermally bistable medium. We perform one-dimensional two fluid numerical simulations that describe the neutral and ionized components in the interstellar medium with purely transverse magnetic fields. We find that the clouds that are formed as a consequence of the thermal instability always have a magnetic field strength on the order of a few microgauss irrespective of the initial strength. This shows good agreement with the measurements of the magnetic field strength in diffuse clouds. We show analytically that the final magnetic field strength in the clouds is determined by the balance between ambipolar diffusion and the accumulation of the magnetic field due to condensation.


title: G315.1+2.7: a new Galactic SNR from the AAO/UKST HAlpha survey

abstract: New narrow-band HAlpha imaging and subsequent optical spectra confirm G315.1+2.7, a previously identified candidate supernova remnant (SNR), as a bona-fide Galactic SNR. Present observations are based on independent discovery of filamentary optical emission nebulosity on images of the AAO/UKST HAlpha survey of the southern Galactic plane which were found to co-incide with existing multi-frequency radio detections. Separate medium and high dispersion spectra were taken across two locations of this 11 arcmin N-S aligned optical filament. The resulting spectral signatures were found to strongly confirm the SNR identification based on standard emission line ratio discriminators which characterise emission from shock heated gas. The average observed ratios of [SII]/HAlpha=1.13, [NII]/HAlpha=1.43 and [SII] 6717/6731=1.46, together with the simultaneous detection of [OII] at 3727A, [OIII] 5007A and [OI] 6300A, all point to an SNR origin of the observed optical emission. Careful scrutiny of the low-resolution but high sensitivity SHASSA HAlpha survey also revealed a low-level but distinct optical emission arc. This arc precisely correlates with the large, 2.5 degree, north-south angular extent of the proposed new SNR also seen as a fractured structure in the extant radio data. G315.1+2.7 was detected previously at 2400 and 4800 MHz and at 408 and 1420 MHz. We also identified associated radio emission at 843 MHz from the now publicly available SUMSS survey. On the basis of optical imaging and spectra and radio observations at five frequencies, we identify G315.1+2.7 as a new Galactic SNR. The large projected angular extent of the new remnant, together with the distance estimate of ~1.7 kpc and diameter of ~80pc, make G315.1+2.7 one of the largest remnants known.


title: NUclei of GAlaxies: V. Radio emission in 7 NUGA sources

abstract: We present high angular resolution radio snap-shot observations of seven nearby low-luminosity active galaxies (LLAGN) from the NUclei of GAlaxies (NUGA) survey. The observations were conducted with MERLIN and EVN/VLBI at 18cm and 6cm. At all observed angular resolutions and frequencies, we find indications for extended emission in about ~40% of the sources, consistent with the decrease of flux with increasing angular resolution. The extended components resemble jet emission in a majority of cases, consistent with the optically thin synchrotron emission implied by their steep spectra. We consider the compact 6cm EVN/VLBI radio emission of our sources in the context of the ``fundamental plane'' that previous LLAGN studies identified within the three-dimensional parameter space of radio luminosity, X-ray luminosity, and black hole mass. We demonstrate, using NGC7217 and NGC1068 as particular examples, that high-resolution, multi-epoch radio observations offer useful information about the origin of offsets from the fundamental plane.


title: The matter power spectrum as a tool to discriminate dark matterr-dark energy interactions

abstract: The coincidence problem of late cosmic acceleration gets significantly alleviated when a suitable interaction between matter and dark energy, either of phantom type or not, enters the picture. We show that a class of models featuring this interaction fares rather well when contrasted with the anisotropies of the CMBR and the matter power spectrum. The latter test is very sensible to the interaction and may be used to discriminate between different models.


title: Reinterpreting dark energy through backreaction: the minimally coupled morphon field

abstract: In the context of averaged cosmologies, the effective equations can be written in the form of "regional" Friedmannian equations with additional sources arising from the so-called backreaction of inhomogeneities. We propose a mean field description of this backreaction in terms of a regionally homogeneous scalar field: this provides a physical motivation to the phenomenological scalar fields generically called quintessence fields. We explicitly reconstruct the potential of the scalar field for a one-parameter family of scaling solutions to the backreaction problem, showing that it entails most of the standard scalar fields including e.g. standard and phantom quintessence scenarii.


title: Static spherically symmetric perfect fluid solutions in $f(R)$ theories of gravity

abstract: Static spherically symmetric perfect fluid solutions are studied in metric $f(R)$ theories of gravity. We show that pressure and density do not uniquely determine $f(R)$ ie. given a matter distribution and an equation state, one cannot determine the functional form of $f(R)$. However, we also show that matching the outside Schwarzschild-de Sitter-metric to the metric inside the mass distribution leads to additional constraints that severely limit the allowed fluid configurations.


title: Asymmetric neutrino emission in quark matter and pulsar kicks

abstract: The puzzling phenomenon of pulsar kicks, i.e. the observed large escape velocities of pulsars out of supernova remnants, is examined for compact stars with a strange quark matter core. The direct Urca process in quark matter is studied in the presence of a strong magnetic field. Conditions for an asymmetric emission of the produced neutrinos are worked out in detail, giving constraints on the temperature, the strength of the magnetic field and the electron chemical potential in the quark matter core. In addition, the neutrino mean free paths for quark matter and a possible hadronic mantle are considered.


title: Neutron stars and quark stars: Two coexisting families of compact stars?

abstract: The mass-radius relation of compact stars is discussed with relation to the presence of quark matter in the core. The existence of a new family of compact stars with quark matter besides white dwarfs and ordinary neutron stars is outlined.


title: On the structure of the turbulent interstellar atomic hydrogen. I- Physical characteristics

abstract: {We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}


title: On the structure of the turbulent interstellar atomic hydrogen. II- First comparison between observation and theory

abstract: {It is necessary to understand the dynamics of the atomic gas to use complex modeling and to carry out detailed comparisons between theoretical models and observations.}{In a companion paper, we present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen. Here, we further characterize these simulations and we compare our results with various observations.}{We give statistics of the column density and velocity along the line of sight and show that they compare favorably with observations of high-latitude lines of sight. We compute synthetic HI spectra and qualitatively discuss the information that could be inferred if these spectra were observed. Finally, we extract CNM clouds and study their physical properties finding strong similarities with real clouds. In particular, we find that the clouds follow Larson-type relations, i.e. $M \propto L^\gamma$, where $\gamma \simeq 1.7$ (we speculate that in 3D, $\gamma \simeq 2.5$) and $\sqrt{<\delta v^2>} \propto L^{0.4}$. We also find that the distribution, ${\cal N}(N)$, of the column density, $N$, of the CNM structures formed in the simulation follows ${\cal N}(N) \propto N^{-1.2}$ which is marginally compatible with the observational result obtained by Heiles & Troland (2005). From the mass-size relation and the mass spectrum, we derive an exponent for the column density distribution close to the value obtained in the numerical simulation.} {We conclude that the simulations reproduce various observational features reasonably well....}


title: Neutron stars in Microquasars

abstract: We discuss the "basic" condition for an accreting neutron star to become a microquasar, i.e. ejecting relativistic particles orthogonal to the accretion disk instead of confining disk-material down to the magnetic poles and creating the two emitting caps typical for a X-ray pulsar. Jet creation is prevented for B >/= 10^12 G independent of the accretion rate. This excludes the possibility for a classic X-ray pulsar to develop a "microquasar-phase" and is consistent with the lack of radio emission from such pulsar systems. Millisecond accretion-powered pulsars, on the contrary, may show a "microquasar-phase", where B < 10^7.5 G is valid, whereas the limit for Z sources is B < 10^8.2 G. The implication of our analysis is that the jet might be the suitable agent of angular momentum sink for millisecond accretion-powered pulsars.


title: Dynamic confinement of jets by magneto-torsional oscillations

abstract: Many quasars and active galactic nuclei (AGN) appear in radio, optical, and X-ray maps, as a bright nuclear sources from which emerge single or double long, thin jets (Thomson et al., 1993). When observed with high angular resolution these jets show structure with bright knots separated by relatively dark regions. High percentages of polarization, sometimes more then 50% in some objects, indicates the nonthermal nature of the radiation which is well explained as the synchrotron radiation of the relativistic electrons in an ordered magnetic field. A strong collimation of jets is most probably connected with ordered magnetic fields. The mechanism of magnetic collimation, first suggested by Bisnovatyi-Kogan et al. (1969), was based on the initial charge separation, leading to creation of oscillating electrical current, which produces azimuthal magnetic field, preventing jet expansion and disappearance. Here we consider magnetic collimation, connected with torsional oscillations of a cylinder with elongated magnetic field. Instead of initial blobs with charge separation, we consider a cylinder with a periodically distributed initial rotation around the cylinder axis. The stabilizing azimuthal magnetic field is created here by torsional oscillations, where charge separation is not necessary. Approximate simplified model is developed. Ordinary differential equation is derived, and solved numerically, what gives a possibility to estimate quantitatively the range of parameters where jets may be stabilized by torsional oscillations.


title: Forced oscillations in relativistic accretion disks and QPOs

abstract: In this work we explore the idea that the high frequency QPOs observed in LMXBs may be explained as a resonant coupling between the neutron star spin and epicyclic modes of accretion disk oscillations. We propose a new model for these QPOs based on forced oscillations induced in the accretion disk due to a stellar asymmetric rotating gravitational or magnetic field. It is shown that particles evolving in a rotating non-axisymmetric field are subject to three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driving force, and a parametric resonance due to the time varying epicyclic frequencies. These results are extends by means of 2D numerical simulations of a simplified version of the accretion disk. The simulations are performed for the Newtonian gravitational potential, as well as for a pseudo-general relativistic potential, which enables us to explore the behavior of the resonances around both rotating neutron stars and black holes. Density perturbations are only significant in the region located close to the inner edge of the disk near the ISCO where the gravitational or magnetic perturbation is maximal. It is argued that the nearly periodic motion induced in the disk will produce high quality factor QPOs. Finally, applying this model to a typical neutron star, we found that the strongest response occurs when the frequency difference of the two modes equals either the spin frequency (for "slow rotators") or half of it (for "fast rotators"). The two main excited modes may both be connected to vertical oscillations of the disk. We emphasize that strong gravity is not needed to excite the modes.


title: Compact star constraints on the high-density EoS

abstract: A new scheme for testing the nuclear matter (NM) equation of state (EoS) at high densities using constraints from compact star (CS) phenomenology is applied to neutron stars with a core of deconfined quark matter (QM). An acceptable EoS shall not to be in conflict with the mass measurement of 2.1 +/- 0.2 solar mass (1 sigma level) for PSR J0751+1807 and the mass radius relation deduced from the thermal emission of RX J1856-3754. Further constraints for the state of matter in CS interiors come from temperature-age data for young, nearby objects. The CS cooling theory shall agree not only with these data, but also with the mass distribution inferred via population synthesis models as well as with LogN-LogS data. The scheme is applied to a set of hybrid EsoS with a phase transition to stiff, color superconducting QM which fulfills all above constraints and is constrained otherwise from NM saturation properties and flow data of heavy-ion collisions. We extrapolate our description to low temperatures and draw conclusions for the QCD phase diagram to be explored in heavy-ion collision experiments.