Nonrotating, zero metallicity stars with initial masses 140 < M < 260 solar masses are expected to end their lives as pair-production supernovae (PPSNe), in which an electron-positron pair-production instability triggers explosive nuclear burning. Interest in such stars has been rekindled by recent theoretical studies that suggest primordial molecular clouds preferentially form stars with these masses. Since metal enrichment is a local process, the resulting PPSNe could occur over a broad range of redshifts, in pockets of metal-free gas. Using the implicit hydrodynamics code KEPLER, we have calculated a set of PPSN light curves that addresses the theoretical uncertainties and allows us to assess observational strategies for finding these objects at intermediate redshifts. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches may soon be able to limit the contribution of these very massive stars to < 1% of the total star formation rate density out to z=2 which already provides useful constraints for theoretical models. The planned Joint Dark Energy Mission satellite will be able to extend these limits out to z=6.
We present preliminary results of a new quasar survey being undertaken with multi-colour optical imaging from the Canada-France-Hawaii Telescope. The current data consists of 3.83 sq. deg. of imaging in the i' and z' filters to a 10 sigma limit of z'<23.35. Near-infrared photometry of 24 candidate 5.7<z<6.4 quasars confirms them all to be low mass stars including two T dwarfs and four or five L dwarfs. Photometric estimates of the spectral type of the two T dwarfs are T3 and T6. We use the lack of high-redshift quasars in this survey volume to constrain the z=6 quasar luminosity function. For reasonable values of the break absolute magnitude M*_1450 and faint-end slope alpha, we determine that the bright-end slope beta>-3.2 at 95% confidence. We find that the comoving space-density of quasars brighter than M_1450=-23.5 declines by a factor >25 from z=2 to z=6, mirroring the decline observed for high-luminosity quasars. We consider the contribution of the quasar population to the ionizing photon density at z=6 and the implications for reionization. We show that the current constraints on the quasar population give an ionizing photon density <<30% that of the star-forming galaxy population. We conclude that active galactic nuclei make a negligible contribution to the reionization of hydrogen at z~6.
If general relativity is the correct theory of physics on large scales, then there is a differential equation that relates the Hubble expansion function, inferred from measurements of angular diameter distance and luminosity distance, to the growth rate of large scale structure. Deviations from this consistency relationship could be the signature of the breakdown of general relativity on cosmological scales. As an example, we consider a universe described by a recently proposed 5-dimensional modified gravity theory. We demonstrate that this leads to an inconsistency within the dark energy parameter space. We propose a procedure to detect this signature of modified gravity theory using different combinations of cosmological probes.
Using the Chandra Advanced CCD Imaging Spectrometer Imaging array (ACIS-I), we have carried out a deep hard X-ray observation of the Galactic plane region at (l,b) ~ (28.5, 0.0), where no discrete X-ray source had been reported previously. We have detected 274 new point X-ray sources (4 sigma confidence) as well as strong Galactic diffuse emission within two partially overlapping ACIS-I fields (~250 arcmin^2in total). Sum of all the detected point source fluxes accounts for only ~ 10 % of the total X-ray flux in the field of view. Even hypothesizing a new population of much dimmer and numerous Galactic point sources, the total observed X-ray flux cannot be explained. Therefore, we conclude that X-ray emission from the Galactic plane has truly diffuse origin. Only 26 point sources were detected both in the soft and hard bands, indicating that there are two distinct classes of the X-ray sources distinguished by the spectral hardness ratio. Surface number density of the hard sources is only slightly higher than that measured at the high Galactic latitude regions, indicating that majority of the hard sources are background AGNs. Following up the Chandra observation, we have performed a near-infrared (NIR) survey with SOFI at ESO/NTT. Almost all the soft X-ray sources have been identified in NIR and their spectral types are consistent with main-sequence stars, suggesting most of them are nearby X-ray active stars. On the other hand, only 22 % of the hard sources had NIR counterparts, which are presumably Galactic. From X-ray and NIR spectral study, they are most likely to be quiescent cataclysmic variables. We have also carried out a precise spectral study of the Galactic diffuse X-ray emission excluding the point sources.
We apply the recently defined multipole vector framework to the
frequency-specific first-year WMAP sky maps, estimating the low-l multipole
coefficients from the high-latitude sky by means of a power equalization
filter. While most previous analyses of this type have considered only heavily
processed (and foreground-contaminated) full-sky maps, the present approach
allows for greater control of residual foregrounds, and therefore potentially
also for cosmologically important conclusions. The low-l spherical harmonics
coefficients and corresponding multipole vectors are tabulated for easy
reference.
Using this formalism, we re-assess a set of earlier claims of both
cosmological and non-cosmological low-l correlations based on multipole
vectors. First, we show that the apparent l=3 and 8 correlation claimed by Copi
et al. (2004) is present only in the heavily processed map produced by Tegmark
et al. (2003), and must therefore be considered an artifact of that map.
Second, the well-known quadrupole-octopole correlation is confirmed at the 99%
significance level, and shown to be robust with respect to frequency and sky
cut. Previous claims are thus supported by our analysis. Finally, the low-l
alignment with respect to the ecliptic claimed by Schwarz et al. (2004) is
nominally confirmed in this analysis, but also shown to be very dependent on
severe a-posteriori choices. Indeed, we show that given the peculiar
quadrupole-octopole arrangement, finding such a strong alignment with the
ecliptic is not unusual.
We present narrowband images of eta Carinae in the light of [O III] 5007 obtained with HST/WFPC2, as well as a ground-based image in the same emission line with a larger field of view. These images show a thin veil of [O III] emission around eta Car and its ejecta, confirming the existence of an oxygen-bearing ``cocoon'' inferred from spectra. This [O III] veil may be the remnant of the pre-outburst wind of eta Car, and its outer edge probably marks the interface where eta Car's ejecta meet the stellar wind of the nearby O4 V((f)) star HD303308 or other ambient material -- i.e., it marks the ``astropause'' in eta Car's wind. This veil is part of a more extensive [O III] shell that appears to be shaped and ionized by HD303308. A pair of HST images with a 10 yr baseline shows no proper motion, limiting the expansion speed away from eta Car to 12pm13 km/s, or an expansion age of a few times 10^4 yr. Thus, this is probably the decelerated pre-outburst LBV wind of eta Car. The [O III] morphology is very different from that seen in [N II], which traces young knots of CNO-processed material; this represents a dramatic shift in the chemical makeup of material recently ejected by eta Car. This change in the chemical abundances may have resulted from the sudden removal of the star's outer envelope during eta Car's 19th century outburst or an earlier but similar event.
We report the analysis of time-resolved spectroscopy of IP Pegasi in outburst with eclipse mapping techniques to investigate the location and geometry of the observed spiral structures. We were able to obtain an improved view of the spiral structures with the aid of light curves extracted in velocity bins matching the observed range of velocities of the spiral arms combined with a double default map tailored for reconstruction of asymmetric structures. Two-armed spiral structures are clearly seen in all eclipse maps. The arms are located at different distances from the disc centre. The ``blue'' arm is farther out in the disc (R= 0.55 +/- 0.05 R_{L1}) than the ``red'' arm (R= 0.30 +/- 0.05 R_{L1}). There are evidences that the velocity of the emitting gas along the spiral pattern is lower than the Keplerian velocity for the same disc radius. The discrepancy is smaller in the outer arm (measured velocities 10-15 per cent lower than Keplerian) and is more significant in the inner arm (observed velocities up to 40 per cent lower than Keplerian). We measured the opening angle of the spirals from the azimuthal intensity distribution of the eclipse maps to be \phi= 25 +/- 3 degrees. A comparison with similar measurements on data at different outburst stages reveals that the opening angle of the spiral arms in IP Peg decreases while the outbursting accretion disc cools and shrinks, in agreement with the expected evolution of a tidally driven spiral wave. The sub-Keplerian velocities along the spiral pattern and the clear correlation between the opening angle of the spirals and the outburst stage favors the interpretation of these asymmetric structures as tidally-induced spiral shocks.
Observations of the hot gas surrounding Sgr A* and a few other nearby galactic nuclei imply electron and proton mean free paths comparable to the gas capture radius: hot accretion likely proceeds under weakly-collisional conditions in these systems. As a result, thermal conduction, rather than convection, may be important on all scales and affect the global flow properties. The self-similar ADAF solution of Narayan & Yi (1994) is generalized to include a saturated form of thermal conduction, as is appropriate for the weakly-collisional regime of interest. Conduction provides extra heating and yet it reduces the free-free radiative efficiency of the accretion flow (by potentially large factors). These idealized solutions suggest that thermal conduction may be an important physical ingredient to understand hot accretion onto dim accreting black holes. Conduction could also play a role in reducing the rate at which black holes capture ambient gas and in providing an evaporation mechanism for an underlying cold thin disk.
The Cryogenic Dark Matter Search (CDMS-II) employs low-temperature Ge and Si detectors to seek Weakly Interacting Massive Particles (WIMPs) via their elastic scattering interactions with nuclei. Simultaneous measurements of both ionization and phonon energy provide discrimination against interactions of background particles. For recoil energies above 10 keV, events due to background photons are rejected with >99.99% efficiency. Electromagnetic events very near the detector surface can mimic nuclear recoils because of reduced charge collection, but these surface events are rejected with >96% efficiency by using additional information from the phonon pulse shape. Efficient use of active and passive shielding, combined with the the 2090 m.w.e. overburden at the experimental site in the Soudan mine, makes the background from neutrons negligible for this first exposure. All cuts are determined in a blind manner from in situ calibrations with external radioactive sources without any prior knowledge of the event distribution in the signal region. Resulting efficiencies are known to ~10%. A single event with a recoil of 64 keV passes all of the cuts and is consistent with the expected misidentification rate of surface-electron recoils. Under the assumptions for a standard dark matter halo, these data exclude previously unexplored parameter space for both spin-independent and spin-dependent WIMP-nucleon elastic scattering. The resulting limit on the spin-independent WIMP-nucleon elastic-scattering cross-section has a minimum of 4x10^-43 cm^2 at a WIMP mass of 60 GeV/c^2. The minimum of the limit for the spin-dependent WIMP-neutron elastic-scattering cross-section is 2x10^-37 cm^2 at a WIMP mass of 50 GeV/c^2.
We report the results of mapping observations of the bipolar nebula with SiO maser emission, IRAS 19312+1950, in the CO (J=1-0 and J=2-1), 13CO (J=1-0 and J=2-1), C18O (J=1-0), CS (J=2-1), SO (J_K=3_2-2_1) and HCO+ (J=3-2) lines with the Berkeley-Illinois-Maryland Association array. Evolutional status of this source has been evoking a controversy since its discovery, though SiO maser sources are usually identified as late-type stars with active mass loss. In line profiles, two kinematical components are found as reported in previous single-dish observations: a broad pedestal component and a narrow component. Spatio-kinetic properties of a broad component region traced by 12CO lines are roughly explained by a simple spherical outflow model with a typical expanding velocity of an AGB star, though some properties of the broad component region still conflict with properties of a typical AGB spherical outflow. A narrow component region apparently exhibits a bipolar flow. The angular size of the narrow component region is spatially larger than that of a broad component region. Intensity distribution of the CS emission avoids the central region of the source, and that of an SO broad component emission exhibits a small feature peaked exactly at the mapping center. According to the present results, if a broad component really originates in a spherical outflow, an oxygen-rich evolved stellar object seems to be a natural interpretation for the central star of IRAS 19312+1950.
In this paper we consider the use of Gamma Ray Bursts (GRBs) as distance markers to study the unification of dark energy and dark matter in the context of the so-called Generalized Chaplygin Gas (GCG) model. We consider that the GRB luminosity may be estimated from its variability, time-lag, and also use the so-called Ghirlanda relation. We evaluate the improvements one may expect once more GRBs and their redshift become available. We show that although GRBs allow for extending the Hubble diagram to higher redshifts, its use as a dark energy probe is limited when compared to SNe Ia. We find that the information from GRBs can provide some bounds on the amount of dark matter and dark energy independently of the equation of state. This is particularly evident for XCDM-type models, which are, for low-redshifts ($z\leq2$), degenerate with the GCG.
Presented is the exact solution of the problem of finding a potential of an inflaton scalar field for which adiabatic perturbations generated during a de Sitter (inflationary) stage in the early Universe have the exactly flat (or, the Harrison-Zeldovich) initial spectrum. This solution lies outside the scope of the slow-roll approximation and higher-order corrections to it. The potential found depends on two arbitrary physical constants, one of those determines an amplitude of perturbations. For small (zero) values of the other constant, a long (infinite) inflationary stage with slow rolling of the inflaton field exists.
In the framework of test-particle and cosmic-ray modified hydrodynamics, we calculate synchrotron emission radial profiles in young ejecta-dominated supernova remnants (SNRs) evolving in an ambient medium which is uniform in density and magnetic field. We find that, even without any magnetic field amplification by Raleigh-Taylor instabilities, the radio synchrotron emission peaks at the contact discontinuity because the magnetic field is compressed and is larger there than at the forward shock. The X-ray synchrotron emission sharply drops behind the forward shock as the highest energy electrons suffer severe radiative losses.
We discuss the issue of energy nonconservation in the virialzation process of spherical collapse model with homogeneous dark energy. We propose an approximation scheme to find the virialization radius as well as modify a scheme recently proposed by Maor and Lahav. By comparing various schemes and estimating the parameter $q$, we conclude that the problem of energy conservation may have sizable effect in fitting models to observations.
Transient spectral features have been discovered in the X-ray spectra of Active Galactic Nuclei, mostly in the 5--7 keV energy range. Several interpretations were proposed for the origin of these features. We examined a model of Doppler boosted blue horns of the iron line originating from a spot in a black hole accretion disc, taking into account different approximations of general relativistic light rays and the resulting shift of energy of photons. We provide a practical formula for the blue horn energy of an intrinsically narrow line and assess its accuracy by comparing the approximation against an exact value, predicted under the assumption of a planar accretion disc. The most accurate approximation provides excellent agreement with the spot orbital radius down to the marginally stable orbit of a non-rotating black hole.
We present a new fast method for simulating pencil-beam type light-cones, using the MLAPM-code (Multi Level Adaptive Particle Mesh) with light-cone additions. We show that by a careful choice of the light-cone orientation, it is possible to avoid extra periodicities in the light-cone. As an example, we apply the method to simulate a 6 Gpc deep light-cone, create the dark matter halo catalogue for the light-cone and study the evolution of haloes from $z=6$ up to the present time. We determine the spatial density of the haloes, their large-scale correlation function, and study the evolution of the mass function. We find a surprisingly simple relation for the dependence of halo maximum mass on redshift, and apply it to derive redshift limits for bright quasars.
We present the analysis of Rossi X-ray Timing Explorer (RXTE) observations of the turn-on phase of a 35 day cycle of the X-ray binary Her X-1. During the early phases of the turn-on, the energy spectrum is composed of X-rays scattered into the line of sight plus heavily absorbed X-rays. The energy spectra in the 3-17 keV range can be described by a partial covering model, where one of the components is influenced by photoelectric absorption and Thomson scattering in cold material plus an iron emission line at 6.5 keV. In this paper we show the evolution of spectral parameters as well as the evolution of the pulse profile during the turn-on. We describe this evolution using Monte Carlo simulations which self-consistently describe the evolution of the X-ray pulse profile and of the energy spectrum.
It is shown that the existence of radiation-induced solitary waves in hot plasmas of accretion disks depends on the radial temperature profile.
We present the discovery of nine quasars behind the Large Magellanic Cloud, with emission redshifts ranging from 0.07 to 2.0. Six of them were identified as part of the systematic variability-based search for QSOs in the objects from the OGLE-II database. Combination of variability-based selection of candidates with the candidates' colours appears to be a powerful technique for identifying quasars, potentially reaching ca. 50% efficiency. We report an apparent correlation between variability magnitude and variability timescale, which - if confirmed - could put even more constraints on QSO candidate selection. The remaining three quasars were identified via followup spectroscopy of optical counterparts to X-ray sources found serendipitously by the Chandra X-ray Observatory satellite. Even though the locations of the candidates were quite uniformly distributed over the LMC bar, the confirmed QSOs all appear near the bar's outskirts.
Polar magnetospheric gaps consume a fraction of the electric potential that develops accross open field lines. This effect modifies significantly the structure of the axisymmetric pulsar magnetosphere. We present numerical stead-state solutions for various values of the gap potential. We show that a charge starved magnetosphere contains significantly less electric current than one with freely available electric charges. As a result, electromagnetic neutron star braking becomes inefficient. We argue that the magnetosphere may spontaneously rearrange itself to a lower energy configuration through a dramatic release of electromagnetic field energy and magentic flux. Our results might be relevant in understanding the recent December 27, 2004 burst observed in SGR 1806-20.
In the present work we report on several CH stars identified in a sample of
Faint High Latitude Carbon stars from Hamburg survey and discuss their medium
resolution spectra covering a wavelength range 4000 - 6800 \AA . Estimation of
the depths of bands (1,0) $^{12}$C$^{12}$C ${\lambda}$4737 and (1,0)
$^{12}$C$^{13}$C ${\lambda}$4744 in these stars indicate isotopic ratio
$^{12}$C/$^{13}$C ${\sim}$ 3, except for a few exceptions; these ratios are
consistent with existing theories of CH stars evolution. The stars of Hamburg
survey, a total of 403 objects were reported to be carbon star candidates with
strong C$_{2}$ and CN molecular bands.
In the first phase of observation, we have acquired spectra of ninety one
objects. Inspection of the objects spectra show fifty one objects with C$_{2}$
molecular bands in their spectra of which thirteen stars have low flux below
about 4300 \AA . Twenty five objects show weak or moderate CH and CN bands,
twelve objects show weak but detectable CH bands in their spectra and there are
three objects which do not show any molecular bands due to C$_{2}$, CN or CH in
their spectra. Objects with C$_{2}$ molecular bands and with good signals
bluewards of 4300 \AA which show prominent CH bands in their spectra are
potential candidate CH stars. Thirty five such candidates are found in the
present sample of ninty one objects observed so far. The set of CH stars
identified could be the targets of subsequent observation at high resolution
for a detail and comprehensive analysis for understanding their role in early
Galactic chemical evolution.
We investigate the properties of the pairwise velocity dispersion as a function of galaxy luminosity in the context of a halo model. We derive the distribution of velocities of pairs at a given separation taking into account both one-halo and two-halo contributions. We show that pairwise velocity distribution in real space is a complicated mixture of host-satellite, satellite-satellite and two-halo pairs. The peak value is reached at around 1 Mpc/h and does not reflect the velocity dispersion of a typical halo hosting these galaxies, but is instead dominated by the satellite-satellite pairs in high mass clusters. This is true even for cross-correlations between bins separated in luminosity. As a consequence the velocity dispersion at a given separation can decrease with luminosity, even if the underlying typical halo host mass is increasing, in agreement with some recent observations. We compare our findings to numerical simulations and find a good agreement. Numerical simulations also suggest a luminosity dependent velocity bias, which depends on the subhalo mass. We develop models of the auto and cross-correlation function of luminosity subsamples of galaxies in the observable r_proj - pi space and argue that the fingers of god effects extend to large separations and depend both on velocity and density correlations. Finally, we develop a method to extract an average pair-weighted velocity dispersion, dominated by satellite-satellite pairs, from the observations and show that it can estimate it to about 10%.
A detailed analysis is presented of the dominant ionised knot in the halo of the planetary nebula NGC 6543. Observations were made at high spectral and spatial resolution of the [OIII] 5007 line using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. A 20-element multislit was stepped across the field to give almost complete spatial coverage of the knot and surrounding halo. The spectra reveal, for the first time, gas flows around the kinematically inert knot. The gas flows are found to have velocities comparable to the sound speed as gas is photo-evaporated off an ionised surface. No evidence is found of fast wind interaction with the knot and we find it likely that the fast wind is still contained in a pressure-driven bubble in the core of the nebula. This rules out the possibility of the knot having its origin in instabilities at the interface of the fast and AGB winds. We suggest that the knot is embedded in the slowly expanding Red Giant wind and that its surfaces are being continually photoionised by the central star.
At the extreme densities and temperatures typical of the central engine of GRBs, the accreting torus is cooled mainly by advection and by neutrino emission. The latter process is dominated by electron and positron capture onto nucleons ($\beta$ reactions). We calculate the reaction rates and the nuclear composition of matter, assuming that the torus consists of helium, eletron-positron pairs, free neutrons and protons. After determining the equation of state and solving for the disk structure for a given initial accretion rate, we subsequently follow its time evolution. We find that, for accretion rates of the order of $10 M_{\odot}$/s, likely typical for the early stages of the accretion event, the disk becomes unstable, giving rise to variable energy output. This instability may play an important role for producing internal shocks.
We present new optical spectroscopy of the eclipsing binary pulsar Sk 160/SMC X-1. From the He I absorption lines, taking heating corrections into account, we determine the radial velocity semi-amplitude of Sk 160 to be 21.8 +/- 1.8 km/s. Assuming Sk 160 fills its Roche-lobe, the inclination angle of the system is i=65.3 deg +/- 1.3 deg and in this case we obtain upper limits for the mass of the neutron star as Mx = 1.21 +/- 0.10 Msolar and for Sk 160 as Mo= 16.6 +/- 0.4 Msolar. However if we assume that the inclination angle is i=90 deg, then the ratio of the radius of Sk 160 to the radius of its Roche-lobe is beta = 0.79 +/- 0.02, and the lower limits for the masses of the two stars are Mx = 0.91 +/- 0.08 Msolar and Mo = 12.5 +/- 0.1 Msolar. We also show that the HeII 4686A emission line tracks the motion of the neutron star, but with a radial velocity amplitude somewhat less than that of the neutron star itself. We suggest that this emission may arise from a hotspot where material accreting via Roche lobe overflow impacts the outer edge of an accretion disc.
A sample of selected active galactic nuclei (AGN) was observed in 2003 and 2004 with the High Energy Stereoscopic System (H.E.S.S.), an array of imaging atmospheric-Cherenkov telescopes in Namibia. The redshifts of these candidate very-high-energy (VHE, >100 GeV) gamma-ray emitters range from z=0.00183 to z=0.333. Significant detections were already reported for some of these objects, such as PKS 2155-304 and Markarian 421. Marginal evidence (3.1 sigma) for a signal is found from large-zenith-angle observations of Markarian 501, corresponding to an integral flux of I(>1.65 TeV) = (1.5 +/- 0.6 (stat) +/- 0.3 (syst)) x 10^{-12} cm^{-2} s^{-1} or ~15% of the Crab Nebula flux. Integral flux upper limits for 19 other AGN, based on exposures of ~1 to ~8 hrs live time, and with average energy thresholds between 160 GeV and 610 GeV, range from 0.4% to 5.1% of the Crab Nebula flux. All the upper limits are the most constraining ever reported for these objects.
A large population of Intermediate Mass Black Holes (IMBHs) might be produced at early cosmic times as a left over of the evolution of the very massive first stars. Accretion onto IMBHs provides a source of (re)ionizing radiation. We show that the baryon mass fraction locked into IMBHs and their growth is strongly constrained by the observed residual Soft X-ray Background (SXRB) intensity. Thus, unless they are extremely X-ray quiet, miniquasars must be quite rare and/or have a short shining phase. As a byproduct, we find that miniquasars can not be the only source of reionization and that their alleged contribution to the near infrared bands is completely negligible.
Primordial stars are expected to be very massive and hot, producing copious amounts of hard ionizing radiation. The best place to study hard ionizing radiation in the local universe is in very metal-deficient Blue Compact Dwarf (BCD) galaxies. We have carried out a MMT spectroscopic search for [Ne V] 3426 (ionization potential of 7.1 Ryd), [Fe V] 4227 (ionization potential of 4 Ryd) and He II 4686 (ionization potential of 4 Ryd) emission in a sample of 18 BCDs. We have added data from previous work and from the Data Release 3 of the Sloan Digital Sky Survey. In total, we have assembled a BCD high-ionization sample with [Ne V] emission in 4 galaxies, [Fe V] emission in 15 galaxies and He II emission in 465 galaxies. With this large sample, we have reached the following conclusions. There is a general trend of higher [Ne V], [Fe V] and He II emission at lower metallicities. However metallicity is not the only factor which controls the hardness of the radiation. High-mass X-ray binaries and main-sequence stars are probably excluded as the main sources of the very hard ionizing radiation responsible for [Ne V] emission. The most likely source of [Ne V] emission is probably fast radiative shocks moving with velocities > 450 km/s through a dense interstellar medium with an electron number density of several hundreds per cm^-3 and associated with supernova explosions of the most massive stars. These have masses of ~ 50 - 100 Msun and are formed in very compact super-star clusters. The softer ionizing radiation required for He II emission is likely associated with less massive evolved stars and/or radiative shocks moving through a less dense interstellar medium.
We present microlensing calculations for a Galactic model based on Han & Gould (2003), which is empirically normalised by star counts. We find good agreement between this model and data recently published by the MACHO and OGLE collaborations for the optical depth in various Galactic fields, and the trends thereof with Galactic longitude l and latitude b. We produce maps of optical depth and, by adopting simple kinematic models, of average event timescales for microlensing towards the Galactic bulge. We also find that our model predictions are in reasonable agreement with the OGLE data for the expected timescale distribution. We show that the fractions of events with very long and short timescales due to a lens of mass M are weighted by M^2 n(M)dM and M^(-1) n(M)dM respectively, independent of the density and kinematics of the lenses.
The binary star HD 45166 has been observed since 1922 but its orbital period has not yet been found. It is considered a peculiar Wolf-Rayet star, and its assigned classification varied along the years. High-resolution spectroscopic observations show that the spectrum, in emission and in absorption, is quite rich. The emission lines have great diversity of widths and profiles. The Hydrogen and Helium lines are systematically broader than the CNO lines. Assuming that HD 45166 is a double-line spectroscopic binary, it presents an orbital period of P = 1.596 days, with an eccentricity of e = 0.18. In addition, a search for periodicity using standard techniques reveals that the emission lines present at least two other periods, of 5 hours and of 15 hours. The secondary star has a spectral type of B7 V and, therefore, should have a mass of about 4.8 solar masses. Given the radial velocity amplitudes, we determined the mass of the hot (primary) star as being 4.2 solar masses and the inclination angle of the system, i = 0.77 degr. As the eccentricity of the orbit is non zero, the Roche lobes increase and decrease as a function of the orbital phase. At periastron, the secondary star fills its Roche lobe. The distance to the star has been re-determined as d = 1.3 kpc and a color excess of E(B-V)=0.155 has been derived. This implies an absolute B magnitude of -0.6 for the primary star and -0.7 for the B7 star. We suggest that the discrete absorption components (DACs) observed in the ultraviolet with a periodicity similar to the orbital period may be induced by periastron events.
In weakly magnetized, dilute plasmas in which thermal conduction along magnetic field lines is important, the usual convective stability criterion is modified. Instead of depending on entropy gradients, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) < 0, which we refer to as the Balbus criterion. We refer to the convective instability that results in this regime as the magnetothermal instability (MTI). We use numerical MHD simulations which include anisotropic electron heat conduction to follow the growth and saturation of the MTI in two-dimensional, plane parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured in the simulations agree with the weak field dispersion relation. We investigate the effect of strong fields and isotropic conduction on the linear properties and nonlinear regime of the MTI. In the nonlinear regime, the instability saturates and convection decays away, when the atmosphere becomes isothermal. Sustained convective turbulence can be driven if there is a fixed temperature difference between the top and bottom edges of the simulation domain, and if isotropic conduction is used to create convectively stable layers that prevent the formation of unresolved, thermal boundary layers. The largest component of the time-averaged heat flux is due to advective motions. These results have implications for a variety of astrophysical systems, such as the temperature profile of hot gas in galaxy clusters, and the structure of radiatively inefficient accretion flows.
Airborne and space-based low-resolution spectroscopy in the 1980s discovered tantalizing quantitative relationships between the gas phase C/O abundance ratio in planetary nebulae (PNe) and the fractions of total far-infrared luminosity radiated by the 7.7 and 11.3 micron bands (the C=C stretch and C-H bend, respectively), of polycyclic aromatic hydrocarbons. Only a very small sample of nebulae was studied in this context, limited by airborne observations of the 7.7 micron band, or the existence of adequate IRAS Low Resolution Spectrometer data for the 11.3 micron band. To investigate these trends further, we have expanded the sample of planetaries available for this study using Infrared Space Observatory (ISO) low-resolution spectra secured with the Short Wavelength Spectrometer (SWS) and the Long Wavelength Spectrometer (LWS). The new sample of 43 PNe, of which 17 are detected in PAH emission, addresses the range from C/O=0.2-13 with the objective of trying to delineate the pathways by which carbon dust grains might have formed in planetaries. For the 7.7-micron and 11.3- micron bands, we confirm that the ratio of band strength to total infrared luminosity is correlated with the nebular C/O ratio. Expressed in equivalent width terms, the cut-on C/O ratio for the 7.7-micron band is found to be 0.6(+0.2,-0.4), in good accord with that found from sensitive ground-based measurements of the 3.3-micron band too.
We present in this study flux-calibrated integrated spectra in the range 3600-6800A for 18 concentrated SMC clusters. Cluster reddening values were estimated by interpolation between the extinction maps of Burstein & Heiles (1982, AJ, 87, 1165) and Schlegel et al. (1998, ApJ, 500, 525). The cluster parameters were derived from the template matching procedure by comparing the line strengths and continuum distribution of the cluster spectra with those of template cluster spectra with known parameters and from the equivalent width (EW) method. In this case, new calibrations were used together with diagnostic diagrams involving the sum of EWs of selected spectral lines. A very good agreement between ages derived from both methods was found. The final cluster ages obtained from the weighted average of values taken from the literature and the present measured ones range from 15 Mr (e.g. L51) to 7 Gyr (K3). Metal abundances have been derived for only 5 clusters from the present sample, while metallicity values directly averaged from published values for other 4 clusters have been adopted. Combining the present cluster sample with 19 additional SMC clusters whose ages and metal abundances were put onto a homogeneous scale, we analyse the age and metallicity distributions in order to explore the SMC star formation history and its spatial extent. By considering the distances of the clusters from the SMC centre instead of their projections onto the right ascension and declination axes, the present age-position relation suggests that the SMC inner disk could have been related to a cluster formation episode which reached the peak ~2.5 Gyr ago. Evidence for an age gradient in the inner SMC disk is also presented.
We present an angular power spectrum analysis of the 2MASS full release extended source catalogue. The main sample used includes 518,576 galaxies below an extinction-corrected magnitude of K=13.5 and limited to |b|>20. The power spectrum results provide an estimate of the galaxy density fluctuations at extremely large scales, r<1000 Mpc. We compare this with mock predictions constructed from the LCDM Hubble Volume mock catalogue. We find that over the range 1<l<100 the 2MASS C_l is steeper than that for the Hubble Volume model. However, in the linear regime (l<30) there is good agreement between the two. We investigate in detail the effects of possible sources of systematic error. Converting linear power spectrum predictions for the form of the three-dimensional matter power spectrum, P(k), and assuming a flat CDM cosmology, a primordial n_s=1 spectrum and negligible neutrino mass, we perform fits to the galaxy angular power spectrum at large linear scales (l<30, corresponding to r>50 Mpc). We obtain constraints on the galaxy power spectrum shape of Gamma=0.14+/-0.02, in good agreement with previous estimates inferred at smaller scales. We also constrain the galaxy power spectrum normalisation to sigma_8 b^2=1.36+/-0.10; in combination with previous constraints on sigma_8 we infer a K-band bias of b=1.27+/-0.04. We are also able to provide weak constraints on Omega_m h and Omega_b/Omega_m. These results are based on the usual assumption that the errors derived from the Hubble Volume mocks are applicable to all other models. If we instead assume that the error is proportional to the C_l amplitude then the constraints weaken; for example it becomes more difficult to reject cosmologies with lower Gamma.
An overall analysis of the structure and stellar content of M11 is presented, thanks to the wide-angle 2MASS spatial coverage. We derive photometric and structural parameters and discuss the spatial dependance of the luminosity and mass functions. Photometric parameters basically agree with previous ones mostly based on the optical. We obtained a core radius of 1.23pc and a tidal radius of 29pc. In particular, the cluster is populous enough so that the tidal radius could be obtained by fitting the three-parameter King profile to the radial distribution of stars. We analyzed the spatial distribution of mass functions, finding that the the slope changes from -0.73 in the core to +2.88 in the outer halo. The spatial distribution of mass function slopes derived from 2MASS agrees with that derived from optical CCD data, which further confirms the reliability of 2MASS data for future analyses of this kind at comparable observational limits. We detect mass segregation up to distances from the center of ~20arcmin. We emphasize that the mass function slope in the core is flatter than anywhere else as a consequence of mass segregation. The derived total cluster mass is ~11000 solar masses.
The recent radial velocity Keck data for GJ876 in Laughlin et al. (2005) is shown to be in good agreement, apart from a 6 m/s scatter, with a theoretical calculation (Nauenberg, 2002) based on orbital parameters from a fit to the earlier Keck data. The time variation of the periods of the inner and outer planets, which are locked in a near 2:1 resonance, are evaluated, and their mean values, Pi and Po, are shown to satisfy closely the relation Po/Pi=2+Po/P, where P is the period for the retrograde precesion of the common mean periastron of these planets.
We reexamine the usefulness of fitting blended lightcurve models to microlensing photometric data. We find agreement with previous workers (e.g. Wozniak & Paczynski) that this is a difficult proposition because of the degeneracy of blend fraction with other fit parameters. We show that follow-up observations at specific point along the lightcurve (peak region and wings) of high magnification events are the most helpful in removing degeneracies. We also show that very small errors in the baseline magnitude can result in problems in measuring the blend fraction, and study the importance of non-Gaussian errors in the fit results. The biases and skewness in the distribution of the recovered blend fraction is discussed. We also find a new approximation formula relating the blend fraction and the unblended fit parameters to the underlying event duration needed to estimate microlensing optical depth.
Evolution of stellar bars in disk galaxies is accompanied by dynamical instabilities and secular changes. Following the vertical buckling instability, the bars are known to weaken dramatically and develop a pronounced peanut/boxy shape when observed edge-on. Using high-resolution N-body simulations of stellar disks embedded in live axisymmetric dark matter halos, we have investigated the long-term changes in the bar morphology, specifically the evolution of the bar size, its vertical structure and exchange of angular momentum. We find that following the initial buckling, the bar resumes its growth from deep inside the corotation radius and follows the Ultra-Harmonic resonance thereafter. We also find that this secular bar growth triggers a spectacular secondary vertical buckling instability which leads to the appearance of characteristic peanut/boxy/X-shaped bulges. The secular bar growth is crucial for the recurrent buckling to develop. Furthermore, the secondary buckling is milder, persists over a substantial period of time, ~3 Gyr, and can have observational counterparts. Overall, the stellar bars show recurrent behavior in their properties and evolve by increasing their linear and vertical extents, both dynamically and secularly. We also demonstrate explicitly that the prolonged growth of the bar is mediated by continuous angular momentum transfer from the disk to the surrounding halo, and that this angular momentum redistribution is resonant in nature -- a large number of lower resonances trap disk and halo particles and this trapping is robust.
We used 1627 faint (15.5< R<23) stars in five fields of the Calar Alto Deep Imaging Survey (CADIS) to estimate the structure parameters of the Galaxy. The results were derived by applying two complementary methods: first by fitting the density distribution function to the measured density of stars perpendicular to the Galactic plane, and second by modelling the observed colors and apparent magnitudes of the stars in the field, using Monte Carlo simulations. The best-fitting model of the Galaxy is then determined by minimising the C-statistic, a modified chisquared. Our model includes a double exponential for the stellar disk with scaleheights h_1 and h_2 and a power law halo with exponent alpha. 24480 different parameter combinations have been simulated. Both methods yield consistent results: the best fitting parameter combination is alpha=3.0 (or alpha=2.5, if we allow for a flattening of the halo with an axial ratio of (c/a)=0.6), h_1=300 pc, h_2=900 pc, and the contribution of thick disk stars to the disk stars in the solar neighbourhood is found to be between 4 and 10%.
We report the discovery of a transient radio source 2.7 arcsec (0.1 pc projected distance) South of the Galactic Center massive black hole, Sagittarius A*. The source flared with a peak of at least 80 mJy in March 2004. The source was resolved by the Very Large Array into two components with a separation of ~0.7 arcsec and characteristic sizes of ~0.2 arcsec. The two components of the source faded with a power-law index of 1.1 +/- 0.1. We detect an upper limit to the proper motion of the Eastern component of ~3 x 10^3 km s^-1 relative to Sgr A*. We detect a proper motion of ~10^4 km s^-1 for the Western component relative to Sgr A*. The transient was also detected at X-ray wavelengths with the Chandra X-ray Observatory and the XMM-Newton telescope and given the designation CXOGC J174540.0-290031. The X-ray source falls in between the two radio components. The maximum luminosity of the X-ray source is ~10^36 erg s^-1, significantly sub-Eddington. The radio jet flux density predicted by the X-ray/radio correlation for X-ray binaries is orders of magnitude less than the measured flux density. We conclude that the radio transient is the result of a bipolar jet originating in a single impulsive event from the X-ray source and interacting with the dense interstellar medium of the Galactic Center.
We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field, F_{\mu\nu\lambda}, to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King beta-models, we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fit by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 LSB, HSB and dwarf galaxies including 58 galaxies that were fit photometrically with the single parameter (M/L)_{stars}. The results there were qualitatively similar to those obtained using Milgrom's phenomenological MOND model, although the determined galaxy masses were quantitatively different and MOND does not show a return to Keplerian behavior at extragalactic distances. The results here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.
We report on correlations between radio luminosity and X-ray timing features in X-ray binary systems containing low magnetic field neutron stars and black holes. The sample of neutron star systems consists of 4U 1728-34, 4U 1820-34, Ser X-1, MXB 1730-335, GX 13+1, the millisecond X-ray pulsars SAX J1808.4-3658 and IGR J00291+5934, and these are compared with the black hole system GX 339-4. The analysis has been done using data from pointed observations of the Rossi X-ray Timing Explorer coordinated with radio observations. In the neutron star systems the radio luminosity L_{R} is correlated with the characteristic frequency of the L_{h} Lorentzian component detected contemporaneously in the power spectrum, and anticorrelated with its strength. Similarly, in the black hole system GX 339-4 L_{R} is correlated with the frequency of the L_{\ell} component in the power spectrum, and anti-correlated with its strength. The index of a power-law fit to the correlation is similar in both cases, L_{R} \propto \nu^{~1.4} and L_{R} \propto (rms)^{-2.3}. At lower timing frequencies, the radio luminosity is further found to be correlated with the characteristic (break) frequency of the L_{b} component of the power spectra in the neutron stars and, marginally, with the equivalent break frequency in GX 339-4. We briefly discuss the coupling between the innermost regions of the accretion disc and the production of the jet and, from the behaviour of the ms accreting X-ray pulsars, the possible role of the NS magnetic field.
The power spectrum is traditionally parameterized by a truncated Taylor series: $ln P(k) = ln P_* + (n_*-1) ln(k/k_*) + {1/2} n'_* ln^2(k/k_*)$. It is reasonable to truncate the Taylor series if $|n'_* ln(k/k_*)| << |n_*-1|$, but it is not if $|n'_* ln(k/k_*)| \gtrsim |n_*-1|$. We argue that there is no good theoretical reason to prefer $|n'_*| << |n_*-1|$, and show that current observations are consistent with $|n'_* ln(k/k_*)| ~ |n_*-1|$ even for $|ln(k/k_*)| ~ 1$. Thus, there are regions of parameter space, which are both theoretically and observationally relevant, for which the traditional truncated Taylor series parameterization is inconsistent, and hence it can lead to incorrect parameter estimations. Motivated by this, we propose a simple extension of the traditional parameterization, which uses no extra parameters, but that, unlike the traditional approach, covers well motivated inflationary spectra with $|n'_*| ~ |n_*-1|$. Our parameterization therefore covers not only standard-slow-roll inflation models but also a much wider class of inflation models. We use this parameterization to perform a likelihood analysis for the cosmological parameters.
Stable baryonic Q-balls, which appear in supersymmetric extensions of the Standard Model, could form at the end of cosmological inflation from fragmentation of the Affleck -- Dine condensate. We reconsider astrophysical constraints on such Q-balls as dark matter candidates. Baryonic Q-balls interact with matter by absorbing the baryon number and, effectively, leading to a rapid baryon number non-conservation. We have recently shown that this process can occur at a much faster rate than that used in previous calculations. As a consequence, stability of neutron stars imposes a stringent constraint on the types of Q-balls that can be dark matter. Only the Q-balls that correspond to baryonic flat directions lifted by baryon-number violating operators are allowed as dark-matter candidates.
The young stellar object MWC297 is a B1.5Ve star exhibiting strong hydrogen emission lines. This object has been observed by the AMBER/VLTI instrument in 2-telescope mode in a sub-region of the K spectral band centered around the Br gamma line at 2.1656 microns. The object has not only been resolved in the continuum with a visibility of 0.50+/-0.10, but also in the Br gamma line, where the flux is about twice larger, with a visibility about twice smaller (0.33+/-0.06). The continuum emission is consistent with the expectation of an optically thick thermal emission from dust in a circumstellar disk. The hydrogen emission can be understood by the emission of a halo above the disk surface. It can be modelled as a latitudinal-dependant wind model and it explains the width, the strength and the visibibility through the emission lines. The AMBER data associated with a high resolution ISAAC spectrum constrains the apparent size of the wind but also its kinematics.
We explore observational constraints on a cosmological brane-world scenario in which the bulk is not empty. Rather, exchange of mass-energy between the bulk and the bane is allowed. The evolution of matter fields to an observer on the brane is then modified due to new terms in the energy momentum tensor describing this exchange. We investigate the constraints from various cosmological observations on the flow of matter from the bulk into the brane. Interestingly, we show that it is possible to have a $\Lambda = 0$ cosmology to an observer in the brane which satisfies standard cosmological constraints including the CMB temperature fluctuations, Type Ia supernovae at high redshift, and the matter power spectrum. This model even accounts for the observed suppression of the CMB power spectrum at low multipoles. In this cosmology, the observed cosmic acceleration is attributable to the flow of matter from the bulk to the brane. A peculiar aspect of this cosmology is that the present dark-matter content of the universe may be significantly larger than that of a $\Lambda$CDM cosmology. Its influence, however, is offset by the dark-radiation term. Possible additional observational tests of this new cosmological paradigm are suggested.
Complex period variations of five W UMa type binaries (AB And, OO Aql, DK Cyg, V566 Oph, U Peg) were investigated by analyzing their O--C diagrams, and several common features were found. Four of the five systems show secular period variations at a constant rate on the order of |dP / dt 1/P | ~ 10E-7 1/yr. In the case of AB And, OO Aql, and U Peg a high-amplitude, nearly one-century long quasi-sinusoidal pattern was also found. It might be explained as light-time effect, or by some magnetic phenomena, although the mathematical, and consequently the physical, parameters of these fits are very problematic, as the obtained periods are very close to the length of the total data range. The most interesting feature of the studied O--C diagrams is a low amplitude (2-4 x 10E-3 d) modulation with a period around 18--20 yr in four of the five cases. This phenomenon might be indirect evidence of some magnetic cycle in late-type overcontact binaries as an analog to the observed activity cycles in RS CVn systems.
We study accretion in binaries hosting an intermediate mass black hole (IMBH) of 1000 solar masses, and a donor star more massive than 15 solar masses. These systems experience an active X-ray phase characterized by luminosities varying over a wide interval, from <10^36 erg/s up to a few 10^40 erg/s typical of the ultra luminous X-ray sources (ULXs). Roche lobe overflow on the zero-age main sequence and donor masses above 20 solar masses can maintain a long-lived accretion phase at the level required to feed a ULX source. In wide systems, wind transfer rates are magnified by the focusing action of the IMBH yielding wind luminosities around 10^38 erg/s. These high mass-IMBH binaries can be identified as progenitors of IMBH-radio pulsar (PSR) binaries. We find that the formation of an IMBH-PSR binary does not necessarely require the transit through a ULX phase, but that a ULX can highlight a system that will evolve into an IMBH-PSR, if the mass of the donor star is constrained to lie within 15 to 30 solar masses. We show that binary evolution delivers the pre-exploding helium core in an orbit such that after explosion, the neutron star has a very high probability to remain bound to the IMBH, at distances of 1-10 AU. The detection of an IMBH-PSR binary in the Milky Way has suffered, so far, from the same small number of statistics limit affecting the population of ULXs in our Galaxy. Ongoing deeper surveys or next generation radio telescopes like SKA will have an improved chance to unveil such intriguing systems. Timing analysis of a pulsar orbiting around an IMBH would weigh the black hole in the still uncharted interval of mass around 1000 solar masses
Microturbulence, i.e. enhanced fluctuations of plasma density, electric and magnetic fields, is of great interest in astrophysical plasmas, but occurs on spatial scales far too small to resolve by remote sensing, e.g., at ~ 1-100 cm in the solar corona. This paper reports spatially resolved observations that offer strong support for the presence in solar flares of a suspected radio emission mechanism, resonant transition radiation, which is tightly coupled to the level of microturbulence and provides direct diagnostics of the existence and level of fluctuations on decimeter spatial scales. Although the level of the microturbulence derived from the radio data is not particularly high, <\Delta n^2 >/n^2 ~ 10^{-5}$, it is large enough to affect the charged particle diffusion and give rise to effective stochastic acceleration. This finding has exceptionally broad astrophysical implications since modern sophisticated numerical models predict generation of much stronger turbulence in relativistic objects, e.g., in gamma-ray burst sources.
We report on the observation of rapid particle acceleration in numerical simulations of relativistic jet-plasma interactions and discuss the underlying mechanisms. The dynamics of a charge-neutral, narrow, electron-positron jet propagating through an unmagnetized electron-ion plasma was investigated using a three-dimensional, electromagnetic, particle-in-cell computer code. The interaction excited magnetic filamentation as well as electrostatic plasma instabilities. In some cases, the inductively generated longitudinal electric fields reached the relativistic plasma wave-breaking limit, and the longitudinal momentum of about half the positrons increased by 50% with a maximum gain exceeding a factor of two.
The nature of the viscosity operative in hot, two-temperature accretion disks around AGN has been a long-standing, unsolved problem. It has been previously suggested that protons, in conjunction with the turbulent magnetic field that is likely to exist in the accretion disk, might be crucial in providing this viscosity. Several authors have recently determined diffusion coefficients for charged particles (cosmic rays) propagating in turbulent magnetic fields by means of extensive Monte Carlo simulations. We use the diffusion coefficients for protons determined by these simulations to find the effective mean free path for protons in hot accretion disks. This in turn yields good estimates of the viscosity due to energetic protons embedded in the turbulent magnetic field of a hot, two-temperature accretion disk. We combine this with a simple two-temperature accretion disk model to determine the Shakura-Sunyaev $\alpha$ viscosity parameter arising out of this mechanism. We find that protons diffusing in the turbulent magentic field embedded in a hot accretion d isk provide a physically reasonable source of viscosity in hot accretion disks around AGN.
The VITRUV project has the objective to deliver milli-arcsecond spectro-images of the environment of compact sources like young stars, active galaxies and evolved stars to the community. This instrument of the VLTI second generation based on the integrated optics technology is able to combine from 4 to 8 beams from the VLT telescopes. Working primarily in the near infrared, it will provide intermediate to high spectral resolutions and eventually polarization analysis. This paper summarizes the result from the concept study led within the Joint Research Activity advanced instruments of the OPTICON program.
Resonant oscillations in warped disks are examined in order to explain high-frequency QPOs and horizontal-branch QPOs in low-mass X-ray binaries. Different from our previous work, addressed to the same subjects, we relax in this paper the assumption that the disks are isothermal in the vertical direction. That is, the pressure, density, and temperature are assumed to be distributed in the vertical direction with a polytropic relation, and the polytropic index changes as the disk state changes. By this generalization and by some subsidiaries we can qualitatively explain, within the framework of our resonance model, observed large frequency variations in neutron-star QPOs and little variations in black-hole QPOs. We consider vertical resonances of g-mode oscillations, since they are most appropriate to explain observations.
I present an analysis of the XMM-Newton observations of four millisecond pulsars, J0437-4715, J2124-3358, J1024-0719, and J0034-0534. The new data provide strong evidence of thermal emission in the X-ray flux detected from the first three objects. This thermal component is best interpreted as radiation from pulsar polar caps covered with a nonmagnetic hydrogen atmosphere. A nonthermal power-law component, dominating at energies E>3 keV, can also be present in the detected X-ray emission. For PSR J0437-4715, the timing analysis reveals that the shape and pulsed fraction of the pulsar light curves are energy dependent. This, together with the results obtained from the phase-resolved spectroscopy, supports the two-component (thermal plus nonthermal) interpretation of the pulsar's X-ray radiation. Highly significant pulsations have been found in the X-ray flux of PSRs J2124-3358 and J1024-0719. For PSR J0034-0534, a possible X-ray counterpart of the radio pulsar has been suggested. The inferred properties of the detected thermal emission are compared with predictions of radio pulsar models.
Dusty tori have been suggested to play a crucial role in determining the
physical characteristics of active galactic nuclei (AGN), but investigation of
their properties has stalled for lack of high resolution mid-IR imaging.
Recently, a long-awaited breakthrough in this field was achieved: NGC 1068, a
nearby AGN, was the first extragalactic object to be observed with a mid-IR
interferometer, thereby obtaining the needed angular resolution to study the
alleged torus.
In this proceeding, first the field of AGN research is briefly reviewed, with
an emphasis on models of dusty tori. Second, the general properties of the key
object NGC 1068 are discussed. Third, the MIDI mid-IR interferometric data set
is presented together with a first attempt to interpret this data in the
context of tori models. Fourth, preliminary MIDI interferometric spectra of the
nucleus of the nearby starbursting galaxy Circinus are presented. Finally, we
briefly discuss the prospects of ESA's Darwin mission for observing nearby and
distant AGN. This mission will allow detailed mapping of tori of low luminosity
AGN such as NGC 1068 up to redshifts of 1 - 2 and more luminous AGN up to
redshift of 10 and beyond (abridged).
We propose a strategy for determining the optimal softening length in N-body simulations. This strategy ensures that the ensemble-average acceleration error is minimal at a small radius r_0, and that the error never exceeds the optimal value beyond r_0. Analytic formulae are derived for calculating the acceleration error for any halo profile or softened gravity so that the optimal softening length, epsilon_opt, can be determined without N-body force evaluations. For Navarro-Frenk-White halos with concentration parameter c = 5, we find, by setting r_0 to one percent of the virial radius r_v, that epsilon_opt = 0.11 N^{-0.20} for S2 softening (Hockney & Eastwood) and epsilon_opt = 0.065 N^{-0.26} for Plummer softening, where N is the number of particles within the virial radius, and epsilon_opt is in units of r_v. We also estimate that current state-of-the-art halo simulations suffer >~ 6% acceleration errors at 0.01 r_v, which grow rapidly toward the center of the halo. Given that the mass within 0.01 r_v is only a few thousandths of the halo virial mass, even a few percent acceleration errors may contribute significantly to the uncertainties in the inner slope of the halo.
We propose an alternative formalism to simulate temperature anisotropy maps of CMB in $\Lambda$CDM universes with nontrivial spatial topologies. This formalism avoids the need to explicitly compute the eigenmodes of the Laplacian operator in the spatial sections. Instead, the covariance matrix of the coefficients of the spherical harmonic decomposition of the temperature anisotropies is expressed in terms of the elements of the covering group of the space. We obtain a decomposition of the correlation matrix that isolates the topological contribution to the CMB temperature anisotropies out of the ``noise'' due to the simply connected contribution. An additional decomposition of this topological signature of the correlation matrix for any complicated topology allows us to compute it in terms of correlation matrizes corresponding to simpler topologies, for which closed quadrature formulae might be derived. Another simplification in the calculations is obtained by using invariance properties of the correlation matrix due to some symmetries of the quocient manifolds. We also use the latter decomposition to show that temperature anisotropy maps of the CMB of (not too large) multiply connected universes must show ``patterns of alignement'', and propose a method to look for these patterns, thus opening the door to the development of new methods for detecting topology of our Universe even when the injectivity radius of space is slightly larger than the radius of the last scattering surface. We illustrate all these features with the simplest examples, those of flat homogeneous manifolds, i.e., torii, with special attention given to the cylinder, i.e., $T^1$ topology.
The old open cluster M67 is an ideal testbed for current cluster evolution models because of its dynamically evolved structure and rich stellar populations that show clear signs of interaction between stellar, binary and cluster evolution. Here we present the first truly direct N-body model for M67, evolved from zero age to 4 Gyr taking full account of cluster dynamics as well as stellar and binary evolution. Our preferred model starts with 12000 single stars and 12000 binaries placed in a Galactic tidal field at 8.0 kpc from the Galactic Centre. Our choices for the initial conditions and for the primordial binary population are explained in detail. At 4 Gyr, the age of M67, the total mass has reduced by 90% as a result of mass loss and stellar escapes. The mass and half-mass radius of luminous stars in the cluster are a good match to observations although the model is more centrally concentrated than observations indicate. The stellar mass and luminosity functions are significantly flattened by preferential escape of low-mass stars. We find that M67 is dynamically old enough that information about the initial mass function is lost, both from the current luminosity function and from the current mass fraction in white dwarfs. The model contains 20 blue stragglers at 4 Gyr which is slightly less than the 28 observed in M67. Nine are in binaries. The blue stragglers were formed by a variety of means and we find formation paths for the whole variety observed in M67. Both the primordial binary population and the dynamical cluster environment play an essential role in shaping the population. A substantial population of short-period primordial binaries (with periods less than a few days) is needed to explain the observed number of blue stragglers in M67.
By exploiting the far-infrared(FIR) and radio correlation, we have performed the Likelihood-Ratio analysis to identify optical counterparts to the far-infrared sources in the Lockman Hole. Using the likelihood ratio analysis and the associated reliability, 44 FIR sources have been identified with radio sources. Redshifts have been obtained for 29 out of 44 identified sources. One hyper-luminous infrared galaxy (HyLIRG) with and four ultraluminous infrared galaxies (ULIRGs) are identified in our sample. The space density of the FIR sources at z = 0.3-0.6 is 4.6\times 10^{-5}Mpc^{-3}, implying a rapid evolution of the ULIRG population. Most of \ISO FIR sources have their FIR-radio ratios similar to star-forming galaxies ARP 220 and M82. At least seven of our FIR sources show evidence for the presence of an active galactic nucleus (AGN) in optical emission lines, radio continuum excess, or X-ray activity. Three out of five (60%) of the ULIRG/HyLIRGs are AGN galaxies. Five of the seven AGN galaxies are within the ROSAT X-ray survey field, and two are within the XMM-Newton survey fields. X-ray emission has been detected in only one source, 1EX030, which is optically classified as a quasar. The non-detection in the XMM-Newton 2-10 keV band suggests a very thick absorption obscuring the central source of the two AGN galaxies. Several sources have an extreme FIR luminosity relative to the optical R-band, L(90\mu\mathrm{m})/L(R) > 500, which is rare even among the local ULIRG population. While source confusion or blending might offer an explanation in some cases, they may represent a new population of galaxies with an extreme activity of star formation in an undeveloped stellar system -- i.e., formation of bulges or young ellipticals.
We report time-resolved spectroscopic observations of the SU Ursae Majoris dwarf nova, YZ Cnc, for 2 nights over 11 hrs during its 2002 January superoutburst. The spectra only show absorption-line profiles in the first day. But the lines display blue and red troughs, with ``W'' profiles in the second day. The radial velocity curve of the absorption troughs and emission peaks of H$\beta$ has an amplitude of $49\pm10$ km s$^{-1}$ and a phase offset of $-0.07\pm0.04$, which are very similar to those measured in quiescence; however, the $\gamma$ velocity deviates strongly from the systemic velocity measured in quiescence, showing variation of the order of $\pm$60 km s$^{-1}$. And large shifts of $\sim$70 km s$^{-1}$ and $\sim$0.09, for the orbital-averaged velocity and phase respectively, are also found in our observations. All these phenomena can be well explained with a precession of an eccentric disk and we conclude that these phenomena are the characteristic products of an eccentric accretion disk.
We report the results of two XMM-Newton observations of the ultra-compact low-mass X-ray binary 4U1850-087 located in the galactic globular cluster NGC6712. A broad emission feature at 0.7keV was detected in an earlier ASCA observation and explained as the result of an unusual Ne/O abundance ratio in the absorbing material local to the source. We find no evidence for this feature and derive Ne/O ratios in the range 0.14-0.21, consistent with that of the interstellar medium. During the second observation, when the source was 10% more luminous, there is some evidence for a slightly higher Ne/O ratio and additional absorption. Changes in the Ne/O abundance ratio have been detected from another ultra-compact binary, 4U1543-624. We propose that these changes result from an X-ray induced wind which is evaporated from an O and Ne rich degenerate donor. As the source X-ray intensity increases so does the amount of evaporation and hence the column densities and abundance ratio of Ne and O.
We compute the rates P at which acoustic energy is injected into the solar radial p modes for several solar models. The solar models are computed with two different local treatments of convection: the classical mixing-length theory (MLT hereafter) and Canuto et al (1996)'s formulation (CGM hereafter). Among the models investigated here, our best models reproduce both the solar radius and the solar luminosity at solar age and the observed Balmer line profiles. For the MLT treatment, the rates P do depend significantly on the properties of the atmosphere whereas for the CGM's treatment the dependence of P on the properties of the atmosphere is found smaller than the error bars attached to the seismic measurements. The excitation rates P for modes associated with the MLT models are significantly underestimated compared with the solar seismic constraints. The CGM models yield values for P closer to the seismic data than the MLT models. We conclude that the solar p-mode excitation rates provide valuable constraints and according to the present investigation clearly favor the CGM treatment with respect to the MLT, although neither of them yields values of P as close to the observations as recently found for 3D numerical simulations.
We present the first results of a variable star search in the field and in the globular clusters of the Fornax dwarf spheroidal galaxy. Variable stars were identified using the Image Subtraction Technique (Alard 2000) on time-series data obtained with the ESO 2.2 m and the Magellan 6.5 m telescopes. The variable star sample includes RR Lyrae stars, Dwarf Cepheids and Anomalous Cepheids. The pulsation properties (namely: periods, light curves, period-amplitude relations and classification in Oosterhoff types) of Fornax variables from the present study are discussed in some detail.
We present the results of 11 elliptical galaxies with strong nebular emission lines during our study of star formation history along the Hubble sequence. After removing the dilution from the underlying old stellar populations by use of stellar population synthesis model, we derive the accurate fluxes of all emission lines for these objects, which are later classified with emission line ratios into one Seyfert 2, six LINERs and four HII galaxies. We also identify one HII galaxy (A1216+04) as a hitherto unknown Wolf-Rayet galaxy from the presence of the Wolf-Rayet broad bump at 4650 \AA. We propose that the star-forming activities in elliptical galaxies are triggered by either galaxy-galaxy interaction or the merging of a small satellite/a massive star cluster, as already suggested by recent numerical simulations.
Higher Criticism is a recently developed statistic for non-Gaussian detection, proposed in Donoho & Jin 2004. We find that Higher Criticism is useful for two purposes. First, Higher Criticism has competitive detection power, and non-Gaussianity is detected at the level 99% in the first year WMAP data. We find that the Higher Criticism value of WMAP is outside the 99% confidence region at a wavelet scale of 5 degrees (99.46% of Higher Criticism values based on simulated maps are below the values for WMAP). Second, Higher Criticism offers a way to locate a small portion of data that accounts for the non-Gaussianity. Using Higher Criticism, we have successfully identified a ring of pixels centered at (l\approx 209 deg, b\approx -57 deg), which seems to account for the observed detection of non-Gaussianity at the wavelet scale of 5 degrees. Note that the detection is achieved in wavelet space first. Second, it is always possible that a fraction of pixels within the ring might deviate from Gaussianity even if they do not appear to be above the 99% confidence level in wavelet space. The location of the ring coincides with the cold spot detected in Vielva et al. 2004 and Cruz et al. 2005.
This is the first paper of a series that will present data and scientific results from the WINGS project, a wide-field, multiwavelength imaging and spectroscopic survey of galaxies in 77 nearby clusters. The sample was extracted from the ROSAT catalogs of X-Ray emitting clusters with constraints on the redshift (0.04<z<0.07) and distance from the galactic plane (|b|>20). The global goal of the WINGS project is the systematic study of the local cosmic variance of the cluster population and of the cluster galaxies as a function of cluster properties and local environment. This data collection will allow to define a local, Zero-Point reference to gauge the cosmic evolution when compared to more distant clusters. The core of the project consists of wide-field optical imaging of the selected clusters in the B and V bands. We have also completed a multi-fiber, medium resolution spectroscopic survey for 51 of the clusters in the master sample. The imaging and spectroscopy data were collected using respectively the WFC@INT and WYFFOS@WHT in the northern hemisphere, and the WFI@MPG and 2dF@AAT in the southern one. In addition, a NIR (JK) survey of ~50 clusters and an H_alpha survey of some 10 clusters have been also programmed. In this paper we briefly outline the global objectives and the main characteristics of the WINGS project. Moreover, the observing strategy and the data reduction of the optical imaging survey (WINGS-OPT) are presented. We have achieved a photometric accuracy of ~0.02mag, reaching completeness to V~23.5. Field size and resolution (FWHM) span the absolute intervals (1.6-2.7)Mpc and (0.7-1.7)kpc, respectively, depending on the redshift and on the seeing. This allows the planned studies to get a valuable description of the local properties of clusters and galaxies in clusters.
We present the results of ROSAT and XMM-Newton observations of the recurrent ultraluminous X-ray source (ULX) NGC 253 ULX1. This transient is one of the few ULXs that was detected during several outbursts. The luminosity reached 1.4 10^39 erg s^-1 and 0.5 10^39 erg s^-1 in the detections by ROSAT and XMM-Newton, respectively, indicating a black hole X-ray binary with a mass of the compact object of >11 M_sun. It is highly unlikely that this source is a foreground source or a background AGN. In the ROSAT detection NGC 253 ULX1 showed significant variability, whereas the luminosity was constant in the detection from XMM-Newton. The XMM-Newton EPIC spectra are well-fit by a bremsstrahlung model (kT=2.24 keV, N_H=1.74 10^20 cm^-2), which can be used to describe a comptonized plasma. No counterpart was detected in the optical I, R, B, NUV and FUV bands to limits of 22.9, 24.2, 24.3, 22 and 23 mag, respectively, pointing at a X-ray binary with a low mass companion.
In the first paper of this series we use the publicly available code Gim2D to model the r- and i-band images of all galaxies in a magnitude-limited sample of roughly 1800 morphologically classified galaxies taken from the Sloan Digital Sky Survey. The model is a concentric superposition of two components, each with elliptical isophotes with constant flattening and position angle. The disk luminosity profile is assumed exponential, while the bulge is assumed to have a de Vaucouleurs or a Sersic profile. We find that the parameters returned by Gim2D depend little on the waveband or bulge profile used; their formal uncertainties are usually small. Nevertheless, for bright galaxies the measured distribution of b/a, the apparent disk flattening, deviates strongly from the expected uniform distribution, showing that the `disk' identified by the code frequently corresponds to an intrinsically 3-dimensional structure rather than to a true thin disk. We correct approximately for this systematic problem using the observed statistics of the b/a distribution and estimate, as a function of absolute magnitude, the mean fractions of galaxy light in disks and in `pure bulge' systems (those with no detectable disk). For the brightest galaxies the disk light fraction is about 10% and about 80% are `pure bulge' systems. For faint galaxies most of the light is in disks and we do not detect a `pure bulge' population. Averaging over the galaxy population as a whole, we find that 54 \pm 2% of the local cosmic luminosity density at both r and i comes from disks and 32 \pm 2% from `pure bulge' systems. The remaining 14 \pm 2% comes from bulges in galaxies with detectable disks.
We consider a 5-dimensional Ricci flat bouncing cosmological model in which the 4-dimensional induced matter contains two components at late times - the CDM+baryons and dark energy. We find that the arbitrary function $f(z)$ contained in the solution plays a similar role as the potential $V(\phi)$ in quintessence and phantom dark energy models. To resolve the coincidence problem, it is generally believed that there is a scaling stage in the evolution of the universe. We analyze the condition for this stage and show that a hyperbolic form of the function $f(z)$ can work on it well. We find that during the scaling stage (before $z\approx 2$), the dark energy behaves like (but not identical to) a cold dark matter with an adiabatic sound speed $c_{s}^{2}\approx 0$ and $p_{x}\approx 0$. After $z\approx 2$, the pressure of dark energy becomes negative. The transition from deceleration to acceleration happens at $z_{T}\approx 0.8$ which, as well as other predictions of the $5D$ model, agree with current observations.
The light signal measured by fluorescence telescopes receives - strongly depending on the shower geometry with respect to the detector - a non-negligible contribution from additionally produced Cherenkov light. This Cherenkov contribution has to be accounted for to determine primary parameters properly. In comparison to the previous ansatz used by other experiments, the impact of a new analytical description of Cherenkov light production in EAS on the Auger event reconstruction is investigated.
NGC 4438 is a highly perturbed spiral with a stellar tidal tail and extraplanar molecular gas, now very HI deficient, crossing the center of the Virgo cluster at high speed. Various authors have attributed the perturbed appearance to the ram pressure of the intracluster medium, the tidal interaction with NGC 4435, and an ISM-ISM collision between the ISM of NGC 4438 and NGC 4435. We present new CO observations covering virtually all of NGC 4438 and the center of NGC 4435 and detailed simulations including all of the above effects. For the first time CO is detected in NGC 4435. In NGC 4438 we find double line profiles at distances up to 40" to the west and south-west and redshifted lines with respect to galactic rotation in the south of the center. The lack of gas to the North and East coupled with the large gaseous extent to the West and the redshifted and double line profiles can only be reproduced with a ram pressure wind. NGC 4438 is most probably on its first passage through the cluster center and has been stripped of its HI only over the past 100 Myr. While an ISM-ISM collision between NGC 4435 and NGC 4438 may occur, the effect is not significant compared to ram pressure and tidal forces, not surprising for the passage of an S0 galaxy 5-10 kpc from the center of NGC 4438. We also detect narrow CO lines, in the absence of detected HI, in the northern tidal arm some 15 kpc from the center of NGC 4438. This can be understood from the simulations assuming a few percent of the gas is too dense to experience the ram pressure wind. NGC 4438 has changed greatly over the past 100 Myr due to its plunge through the center of the Virgo cluster and the interaction with the S0 NGC 4435.
We analyze data on pulsation amplitudes and phases for two $\beta$ Cephei stars, $\delta$ Cet and $\nu$ Eri. Str\"omgren photometry and radial velocity measurements are used simultaneously to obtain constraints on mean parameters of the stars and identification of the excited modes. The inference about the radial mode order and mean star parameters is based on comparison of certain complex parameter, $f$, determined from data, with its counterpart derived from linear nonadiabatic modelling of stellar oscillations. The theoretical $f$ values are very sensitive to the adopted opacity data. In our modelling we rely on the data from OPAL and OP projects. Significant differences were found. New seismic models of $\nu$ Eri were constructed with both the OPAL and OP opacities.
The Laser Interferometer Gravitational Wave Observatory (LIGO) has performed a third science run with much improved sensitivities of all three interferometers. We present an analysis of approximately 200 hours of data acquired during this run, used to search for a stochastic background of gravitational radiation. We place upper bounds on the energy density stored as gravitational radiation for three different spectral power laws. For the flat spectrum, our limit of Omega_0<8.4e-4 in the 69-156 Hz band is ~10^5 times lower than the previous result in this frequency range.
The present letter is aimed at exploring the influence of overshooting during the central helium burning in pre-white dwarf progenitors on the pulsational properties of PG1159 stars. To this end we follow the complete evolution an intermediate-mass white dwarf progenitor from the zero age main sequence through the thermally pulsing and born-again phases to the domain of the PG1159 stars. Our results suggest that the presence of mode-trapping features in the period spacings of these hot pulsating stars could result from structure in the carbon-oxygen core. We find in particular that in order to get enough core structure consistent with observational demands, the occurrence of overshoot episodes during the central helium burning is needed. This conclusion is valid for thick helium envelopes like those predicted by our detailed evolutionary calculations. If the envelope thickness were substantially smaller, then the occurrence of core overshooting would be more difficult to disentangle from the effects related to the envelope transition zones.
We have used ESO's new NIR IFS SINFONI during its Science Verification period to observe the central regions of local starburst galaxies. Being Science Verification observations, the aim was two-fold: to demonstrate SINFONI's capabilities while obtaining information on the nature of starclusters in starburst galaxies. The targets chosen include a number of the brighter clusters in NGC1808 and NGC253. Here we present first results.
We argue that a few per cent of "Early Dark Energy" can be detected by the statistics of nonlinear structures. The presence of Dark Energy during linear structure formation is natural in models where the present tiny Dark-Energy density is related to the age of the Universe rather than a new fundamental small parameter. Generalisation of the spherical collapse model shows that the linear collapse parameter delta_c is lowered. The corresponding relative enhancement of weak gravitational lensing on arc-minute scales lowers the value of sigma_8 inferred from a given lensing amplitude as compared to Lambda-CDM. In presence of Early Dark Energy, structures grow slower, such that for given sigma_8 the number of galaxies and galaxy clusters is substantially increased at moderate and high redshift. For realistic models, the number of clusters detectable through their thermal Sunyaev-Zel'dovich effect at redshift unity and above, e.g. with the Planck satellite, can be an order of magnitude larger than for Lambda-CDM.
We present the analysis of the circular polarization, due to Zeeman splitting, of the water masers around a sample of late-type stars to determine the magnetic fields in their circumstellar envelopes. The magnetic field strengths in the water maser regions around the Mira variable stars U Ori and U Her are shown to be several Gauss while those of the supergiants S Per, NML Cyg and VY CMa are several hundred mG. We also show that large scale magnetic fields permeate the CSE of an evolved star; the polarization of the water masers around VX Sgr reveals a dipole field structure. We shortly discuss the coupling of the magnetic field with the stellar outflow, as such fields could possibly be the cause of distinctly aspherical mass-loss and the resulting aspherical planetary nebulae.
Observations made of the X-ray pulsar EXO 1722-363 using the Proportional Counter Array and All Sky Monitor on board the Rossi X-ray Timing Explorer reveal the orbital period of this system to be 9.741 +/- 0.004 d from periodic changes in the source flux. The detection of eclipses, together with the values of the pulse and orbital periods, suggest that this source consists of a neutron star accreting from the stellar wind of an early spectral type supergiant companion. Pulse timing measurements were also obtained but do not strongly constrain the system parameters. The X-ray spectra can be well fitted with a model consisting of a power law with a high energy cutoff and, for some spectra, a blackbody component with a temperature of approximately 0.85 keV.
We use direct-summation N-body integrations to follow the evolution of binary black holes at the centers of galaxy models with large, constant-density cores. Particle numbers as large as 400K are considered. The results are compared with the predictions of loss-cone theory, under the assumption that the supply of stars to the binary is limited by the rate at which they can be scattered into the binary's influence sphere by gravitational encounters. The agreement between theory and simulation is quite good; in particular, we are able to quantitatively explain the observed dependence of binary hardening rate on N. We do not verify the recent claim of Chatterjee, Hernquist & Loeb (2003) that the hardening rate of the binary stabilizes when N exceeds a particular value, or that Brownian wandering of the binary has a significant effect on its evolution. When scaled to real galaxies, our results suggest that massive black hole binaries in gas-poor nuclei would be unlikely to reach gravitational-wave coalescence in a Hubble time.
The vertical as well as horizontal oscillation modes in a thermally excited two-dimensional (2D) dust coulomb cluster were investigated for particle numbers between and using both a box_tree simulation and an analytical method. The horizontal mode spectra is shown to agree with published results while the vertical mode spectra obtained from the box_tree simulation and the analytical method are shown to agree with one another. The maximum frequency of the vertical modes is shown to be the vertical oscillation frequency of the whole system acting as a solid plane and the frequency is shown to decrease as the mode number l decreases. The oscillation patterns of the vertical modes are also investigated. For clusters with large numbers of particles, the high frequency modes are shown to have oscillation patterns similar in shape to Bessel-Fourier functions with various index m and n. For specified values of n, modes with higher m (except for) have lower frequencies and for specified m, modes with higher n have lower frequencies. For low frequency modes, the largest amplitude particle motion is concentrated in a few inner rings with the outer rings remaining almost motionless. This is in contrast to the horizontal modes where the strongest motion of the particles is concentrated in the inner rings at high frequency. For clusters with small numbers of particles (), these modes where the strongest vertical motion is concentrated in the inner rings does not exist resulting in oscillation magnitudes for all the modes appearing in the shape of Bessel-Fourier functions.
The structure and evolution of close binary stars has been studied using the two-dimensional (2D) stellar structure algorithm developed by Deupree (1995). We have calculated a series of solar composition stellar evolution sequences of binary models, where the mass of the 2D model is 8Msun with a point-mass 5Msun companion. We have also studied the structure of the companion in 2D, by considering the zero-age main-sequence (ZAMS) structure of a 5Msun model with an 8Msun point-mass companion. In all cases the binary orbit was assumed to be circular and co-rotating with the rotation rate of the stars. We considered binary models with three different initial separations, a = 10, 14 and 20Rsun. These models were evolved through central hydrogen burning or until the more massive star expanded to fill its critical potential surface or Roche lobe. The calculations show that evolution of the deep interior quantities is only slightly modified from those of single star evolution. Describing the model surface as a Roche equipotential is also satisfactory until very close to the time of Roche lobe overflow, when the self gravity of the model about to lose mass develops a noticeable aspherical component and the surface time scale becomes sufficiently short that it is conceivable that the actual surface is not an equipotential.
The cosmic expansion history tests the dynamics of the global evolution of the universe and its energy density contents, while the cosmic growth history tests the evolution of the inhomogeneous part of the energy density. Precision comparison of the two histories can distinguish the nature of the physics responsible for the accelerating cosmic expansion: an additional smooth component - dark energy - or a modification of the gravitational field equations. With the aid of a new fitting formula for linear perturbation growth accurate to 0.2%, we separate out the growth dependence on the expansion history and introduce a new growth index parameter \gamma that quantifies the gravitational modification.
Using the extensive multi-wavelength data in the GOODS-North field, we construct and draw comparisons between samples of optical and near-IR selected star-forming and passively evolving galaxies at redshifts 1.4<z<2.6. We find overlap at the 70-80% level in samples of z~2 star-forming galaxies selected by their optical (UGR) and near-IR (BzK) colors when subjected to common K-band limits. Deep Chandra data indicate a 25% AGN fraction among near-IR selected objects, much of which occurs among near-IR bright objects (K<20; Vega). Stacking the X-ray emission for the star-forming galaxies, we find the SFR distributions of UGR, BzK, and J-K>2.3 galaxies (Distant Red Galaxies; DRGs) are very similar as a function of K, with K<20 galaxies having <SFR>~120 Msun/yr, a factor of two to three higher than those with K>20.5. The absence of X-ray emission from the reddest DRGs and BzK galaxies with z-K>3 indicates they must have declining star formation histories to explain their red colors and low SFRs. While the M/L ratio of passively-evolving galaxies may be larger on average, the Spitzer/IRAC data indicate that their inferred stellar masses do not exceed the range spanned by optically selected galaxies, suggesting that the disparity in current SFR may not indicate a fundamental difference between optical and near-IR selected massive galaxies (M* > 10^11 Msun). We consider the contribution of UGR, BzK, DRG, and submillimeter galaxies (SMGs) to the SFRD at z~2, taking into account sample overlap. The total SFRD in the interval 1.4<z<2.6 of UGR and BzK galaxies to K=22, and DRGs to K=21 is ~ 0.10+/-0.02 Msun/yr/Mpc^3. Optically-selected galaxies to R=25.5 and K=22 account for ~70% of this total. Submillimeter galaxies absent from the optical and near-IR selected samples contribute <10% to the total SFRD at these redshifts.
We investigate the microlensing effects on a source star surrounded by a circumstellar disk, as a function of wavelength. In the mid- and far-infrared, where the emission of the system is dominated by the cold dusty disk, the microlensing light curve encodes the geometry and surface brightness profile of the disk. For a system located at the Galactic center, we find typical magnifications to be of order 10-20%, and the event lasts over one year. At around 20 microns, where the emission for the star and the disk are comparable, the difference in the emission areas results in a chromatic microlensing event. Finally, in the near-infrared and visible, where the emission of the star dominates, the fraction of star light directly reflected by the disk slightly modifies the light curve of the system which is no longer that of a point source. In each case, the corresponding light curve can be used to probe some of the disk properties. A fraction of 0.1% to 1% optical microlensing events are expected to be associated with circumstellar disk systems. We show that the lensing signal of the disk can be detected with sparse follow-up observations of the next generation space telescopes. While direct imaging studies of circumstellar disks are limited to the solar neighborhood, this microlensing technique can probe very distant disk systems living in various environments and has the potential to reveal a larger diversity of circumstellar disks.
Narrowband searches for Lyman-alpha emission lines are among the most effective ways to identify high-redshift galaxies, especially because they probe not only the galaxies themselves but also the ionization state of the intergalactic medium (IGM). The observed line strengths depend on the amount by which each photon is able to redshift away from line center before encountering neutral gas and hence on the size distribution of HII regions surrounding the sources. Here, we use an analytic model of that size distribution to study the effects of reionization on the luminosity function of Lyman-alpha emitters and their observed spatial distribution. Our model includes the clustering of high-redshift galaxies and thus contains ionized bubbles much larger than those expected around isolated galaxies. As a result, Lyman-alpha emitting galaxies remain visible earlier in reionization: we expect the number counts to decline by only a factor ~2 (or 10) when the mean ionized fraction falls to x_i~0.75 (or 0.5) in the simplest model. Moreover, the absorption is not uniform across the sky: galaxies remain visible only if they sit inside large bubbles, which become increasingly rare as x_i decreases. Thus, the size distribution also affects the apparent clustering of Lyman-alpha selected galaxies. On large scales, it traces that of the large bubbles, which are more biased than the galaxies. On small scales, the clustering increases rapidly as x_i decreases because large HII regions surround strong galaxy overdensities, so a survey automatically selects only those galaxies with neighbours. The transition between these two regimes occurs at the characteristic bubble size. Hence, large Lyman-alpha galaxy surveys have the potential to measure directly the size distribution of HII regions during reionization.
We present the angular power spectrum of the CMB component extracted with FastICA from the Background Emission Anisotropy Scanning Telescope (BEAST) data. BEAST is a 2.2 meter off-axis telescope with a focal plane comprising 8 elements at Q (38-45 GHz) and Ka (26-36 GHz) bands. It operates from the UC White Montain Research Station at an altitude of 3800 meters. The BEAST CMB angular power spectrum has been already calculated by O'Dwyer et.al. using only the Q band data. With two input channels FastICA returns two possible independent component. We found that one of these two is unphysical while the other is a reasonable CMB component. After a detailed calibration procedure based on Monte-Carlo (MC) simulations we extracted the angular power spectrum for the identified CMB component and found a very good agreement with the already published BEAST CMB angular power spectrum and with the WMAP data.
A large body of theoretical and computational work shows that jets - modelled as magnetized disk winds - exert an external torque on their underlying disks that can efficiently remove angular momentum and act as major drivers of disk accretion. These predictions have recently been confirmed in direct HST measurements of the jet rotation and angular momentum transport in low mass protostellar systems. We review the theory of disc winds and show that their physics is universal and scales to jets from both low and high mass star forming regions. This explains the observed properties of outflows in massive star forming regions, before the central massive star generates an ultracompact HII region. We also discuss the recent numerical studies on the formation of massive accretion disks and outflows through gravitational collapse, including our own work on 3D Adaptive Mesh simulations (using the FLASH code) of the hydromagnetic collapse of an initial rotating, and cooling Bonner-Ebert sphere. Magnetized collapse gives rise to outflows. Our own simulations show that both a jet-like disk wind on sub AU scales, and a larger scale molecular outflow occur (Banerjee and Pudritz 2005).
The sum of sunspots number over an odd numbered 11 yr sunspot cycle exceeds that of its preceding even numbered cycle, and it is well known as Gnevyshev and Ohl rule (or G--O rule) after the names of the authors who discovered it in 1948. The G--O rule can be used to predict the sum of sunspot numbers of a forthcoming odd cycle from that of its preceding even cycle. But this is not always possible because occasionally the G--O rule is violated. So far no plausible reason is known either for the G--O rule or the violation of this rule. Here we showed the epochs of the violation of the G--O rule are close to the epochs of the Sun's retrograde orbital motion about the centre of mass of the solar system (i.e., the epochs at which the orbital angular momentum of the Sun is weakly negative). Using this result easy to predict the epochs of violation of the G--O rule well in advance. We also showed that the solar equatorial rotation rate determined from sunspot group data during the period 1879--2004 is correlated/anti-correlated to the Sun's orbital torque during before/after 1945. We have found the existence of a statistically significant $\sim$ 17 yr periodicity in the solar equatorial rotation rate. The implications of these findings for understanding the mechanism behind the solar cycle and the solar-terrestrial relationship are discussed.
In order to search for high-redshift galaxies beyond $z = 6.6$ in the Subaru Deep Field, we have investigated NB921-dropout galaxies where NB921 is the narrowband filter centered at 919.6 nm with FWHM of 13.2 nm for the Suprime-Cam on the Subaru Telescope. There are no secure NB921-dropout candidates brighter than $z^\prime = 25.5$. Based on this result, we discuss the UV luminosity function of star-forming galaxies at $z > 6.6$.
We outline a framework for understanding the X-ray spectra of high mass accretion rate stellar black holes based on X-ray data from RXTE and ASCA. Three spectral regimes can be separated out by the behaviour of the observed disk luminosity and temperature. The well established "standard regime" is seen when the disk dominates the spectrum, where only a small fraction of the luminosity is emitted in the power law tail. These spectra generally satisfy the standard relation expected for thermal emission from a constant area, namely that the disk bolometric luminosity, Ldisk, is proportional to its maximum temperature, Tin^4. However, at higher luminosities this starts to change to Tin^2. This "apparently standard regime" is still dominated by the disk emission, but this difference luminosity-temperature relation and subtle changes in spectral shape may show that another cooling process is required in addition to radiative cooling. At intermediate luminosities there is an anomalous regime (or weak very high state) where the disk temperature and luminosity are less clearly related. These spectra are characterized by the presence of a much stronger comptonized tail indicating high energy electrons. When observed disk emission is corrected for the the effects of comptonisation then these points lie back on the standard relation. The growth of this comptonising corona is also clearly linked to the quasi-periodic oscillations, as these are observed preferentially in the anomalous regime. This presented picture was found to explain the spectral behavior of both black hole binaries in our Galaxy and LMC. Spectral evolution of several bright ULXs observed with ASCA were also successfully explained in the same picture.
We describe a method of disentangling the composite,0.12-5 microns continuum of symbiotic binaries.The observed SED is determined by the IUE/HST archival spectra and flux-points corresponding to the optical UBVRI and infrared JHKLM photometric measurements. The modeled SED is given by superposition of fluxes from the cool giant, hot stellar source and nebula including the effect of the Rayleigh scattering process and considering influence of the iron curtain absorptions.We applied this method to 21 S-type symbiotic stars during quiescence, activity and eclipses.
The magnetic field in our Galaxy is not well known and is difficult to measure. A spiral regular field in the disk between the Galactic arms is favored by observations, however it is still controversial if the field reverses from arm to arm. The parity of the field across the Galactic plane is also not well established. In this letter we demonstrate that cosmic ray protons in the energy range 10^18 to 10^19 eV, if accelerated near the center of the Galaxy, can probe the large scale structure of the Galactic Magnetic Field. In particular if the field is of even parity, and the spiral field reverses direction from arm to arm, i.e. if it is bi-symmetric (BSS), ultra high energy protons will predominantly come from the Southern Galactic hemisphere, and predominantly from the Northern Galactic hemisphere if the field is of even parity and axi-symmetric (ASS). There is no sensitivity to the BSS or ASS configurations if the field is of odd parity.
Based on a new sample of 355 quasars with significant optical polarization and using complementary statistical methods, we confirm that quasar polarization vectors are not randomly oriented over the sky with a probability often in excess of 99.9%. The polarization vectors appear coherently oriented or aligned over huge (~ 1 Gpc) regions of the sky located at both low (z ~ 0.5) and high (z ~ 1.5) redshifts and characterized by different preferred directions of the quasar polarization. In fact, there seems to exist a regular alternance along the line of sight of regions of randomly and aligned polarization vectors with a typical comoving length scale of 1.5 Gpc. Furthermore, the mean polarization angle \bar{\theta} appears to rotate with redshift at the rate of ~ 30\degr per Gpc. The symmetry of the the \bar{\theta} -z relation is mirror-like, the mean polarization angle rotating clockwise with increasing redshift in North Galactic hemisphere and counter-clockwise in the South Galactic one. These characteristics make the alignment effect difficult to explain in terms of local mechanisms, namely a contamination by interstellar polarization in our Galaxy. While interpretations like a global rotation of the Universe can potentially explain the effect, the properties we observe qualitatively correspond to the dichroism and birefringence predicted by photon-pseudoscalar oscillation within a magnetic field. Interestingly, the alignment effect seems to be prominent along an axis not far from preferred directions tentatively identified in the CMB maps. Although many questions and more particularly the interpretation of the effect remain open, alignments of quasar polarization vectors appear as a promising new way to probe the Universe and its dark components at extremely large scales.
The dependence of atmospheric conditions on altitude and time have to be known at the site of an air shower experiment for accurate reconstruction of extensive air showers and their simulations. The height-profile of atmospheric depth is of particular interest as it enters directly into the reconstruction of longitudinal shower development and of the primary energy and mass of cosmic rays. For the southern part of the Auger Observatory, the atmosphere has been investigated in a number of campaigns with meteorological radio soundings and with continuous measurements of ground-based weather stations. Focussing on atmospheric depth and temperature profiles, temporal variations are described and monthly profiles are developed. Uncertainties of the monthly atmospheres that are currently applied in the Auger reconstruction are discussed.
(abridged) The old Star Cluster (SC) systems surrounding any sofar
investigated galaxy represent a powerful tool for the understanding of the
cosmological evolution of their host galaxies. Phases of enhanced cluster
formation can be identified by reliably age dating the clusters. These are
symptoms of violent star formation, results of intense starbursts as commonly
seen in merging nearby galaxies and high-z forming galaxies. However, the young
SC systems of nearby merging galaxies and the old SC systems of nearby passive
galaxies (most likely resulting from high-z starburst events of these galaxies)
appear to be significantly different in terms of their mass functions. Whether
this difference originates in differences in the formation physics/environment
or is only due to dynamical evolution of the cluster systems is currently a hot
issue. A main diagnostic for the survival probability of a newly formed cluster
is most likely its compactness, since less compact, fluffy clusters are more
vulnerable to destructive interactions with the galactic environment of the
cluster. For these reasons, we developed
1) a tool to significantly improve photometry of nearby SCs (which appear
resolved on high-resolution images, provided e.g. by HST), by accurately
measuring and taking the cluster size into account, and
2) a tool to effectively determine ages and masses (plus metallicities and
extinctions) of SCs from multi-wavelength photometry, ranging from the UV to
the NIR.
One of the main problems in the studies of large scale galaxy structures concerns the relation of the correlation properties of a certain population of objects with those of a selected subsample of it, when the selection is performed by considering physical quantities like luminosity or mass. I consider the case where the sampling is defined as in the simplest thresholding selection scheme of the peaks of a Gaussian random field as well as the case of the extraction of point distributions in high density regions from gravitational N-body simulations. I show that an invariant scale under sampling is represented by the zero-crossing scale of \xi(r). By considering recent measurements in the 2dF and SDSS galaxy surveys I note that the zero-point crossing length has not yet been clearly identified, while a dependence on the finite sample size related to the integral constraint is manifest. I show that this implies that other length scales derived from \xi(r) are also affected by finite size effects. I discuss the theoretical implications of these results, when considering the comparison of structures formed in N-body simulations and observed in galaxy samples, and different tests to study this problem.
We study the stellar and wind properties of a sample of Galactic O dwarfs to track the conditions under which weak winds (i.e mass loss rates lower than ~ 1e-8 Msol/yr) appear. The sample is composed of low and high luminosity dwarfs including Vz stars and stars known to display qualitatively weak winds. Atmosphere models including non-LTE treatment, spherical expansion and line blanketing are computed with the code CMFGEN. Both UV and Ha lines are used to derive wind properties while optical H and He lines give the stellar parameters. Mass loss rates of all stars are found to be lower than expected from the hydrodynamical predictions of Vink et al. (2001). For stars with log L/Lsol > 5.2, the reduction is by less than a factor 5 and is mainly due to the inclusion of clumping in the models. For stars with log L/Lsol < 5.2 the reduction can be as high as a factor 100. The inclusion of X-ray emission in models with low density is crucial to derive accurate mass loss rates from UV lines. The modified wind momentum - luminosity relation shows a significant change of slope around this transition luminosity. Terminal velocities of low luminosity stars are also found to be low. The physical reason for such weak winds is still not clear although the finding of weak winds in Galactic stars excludes the role of a reduced metallicity. X-rays, through the change in the ionisation structure they imply, may be at the origin of a reduction of the radiative acceleration, leading to lower mass loss rates. A better understanding of the origin of X-rays is of crucial importance for the study of the physics of weak winds.
Interplanetary scintillation observations of 48 of the 55 Augusto et al. (1998) flat spectrum radio sources were carried out at 111 MHz using the interplanetary scintillation method on the Large Phased Array (LPA) in Russia. Due to the large size of the LPA beam ($1\degr \times 0.5\degr$) a careful inspection of all possible confusion sources was made using extant large radio surveys: 37 of the 48 sources are not confused. We were able to estimate the scintillating flux densities of 13 sources, getting upper limits for the remaining 35. Gathering more or improving extant VLBI data on these sources might significantly improve our results. This proof-of-concept project tells us that compact ($<1''$) flat spectrum radio sources show strong enough scintillations at 111 MHz to establish/constrain their spectra (low-frequency end).
A new view on our Galaxy has recently emerged, with large consequences on its formation scenarios. Not only new dwarf satellites have been detected, still orbiting and tidally disrupting, but also a multitude of stellar streams or tidal debris have been observed suggesting the formation of the halo through successive accretions. The large scatter in the age-metallicity relation in the solar neighborhood points towards several accretion episodes, while the chemical evolution of the disk requires a more or less continuous gas infall. The global kinematics and morphology refined by large surveys such as 2MASS suggest the existence of two embedded bars, as is frequently observed in external galaxies. The mass of the central black hole has been refined through stellar proper motions, and is compatible with the M$_{bh}$-$\sigma$ relation valid for all bulges. The baryonic dark matter is no longer thought to lie in compact objects, and on the contrary, more dark cold gas is revealed by gamma-ray observations. The star forming history can be built, and confronted to numerical models of galaxy evolution both through hierarchical and secular scenarios. Our Galaxy plays thus the role of a prototype to probe galaxy formation theories, and in particular thin and thick disk formation.
The stellar halos of spiral galaxies bear important chemo-dynamical signatures of galaxy formation. We present here the analysis of 89 semi-cosmological spiral galaxy simulations, spanning ~ 4 magnitudes in total galactic luminosity. These simulations sample a wide variety of merging histories and show significant dispersion in halo metallicity at a given total luminosity - more than a factor of ten in metallicity. Our preliminary analysis suggests that galaxies with a more extended merging history possess halos which have younger and more metal rich stellar populations than the stellar halos associated with galaxies with a more abbreviated assembly. A correlation between halo metallicity and its surface brightness has also been found, reflecting the correlation between halo metallicity and its stellar mass. Our simulations are compared with recent Hubble Space Telescope observations of resolved stellar halos in nearby spirals.
We present the results of a systematic, unbiased search for subpulse modulation of 187 pulsars performed with the Westerbork Synthesis Radio Telescope (WSRT) in the Netherlands at an observing wavelength of 21 cm. Using new observations and archival WSRT data we have expanded the list of pulsars which show the drifting subpulse phenomenon by 42, indicating that at least one in three pulsars exhibits this phenomenon. The real fraction of pulsars which show the drifting phenomenon is likely to be larger than some 55%. The majority of the analysed pulsars show subpulse modulation (170), of which the majority were not previously known to show subpulse modulation and 30 show clear systematic drifting. The large number of new drifters we have found allows us, for the first time, to do meaningful statistics on the drifting phenomenon. We find that the drifting phenomenon is correlated with the pulsar age such that drifting is more likely to occur in older pulsars. Pulsars which drift more coherently seem to be older and have a lower modulation index. There is no significant correlation found between P3 and other pulsar parameters (such as the pulsar age), as has been reported in the past. There is no significant preference of drift direction and the drift direction is not found to be correlated with pulsar parameters. None of the four complexity parameters predicted by different emission models (Jenet & Gil 2003) are shown to be inconsistent with the set of modulation indices of our sample of pulsars. Therefore none of the models can be ruled out based on our observations. We also present results on some interesting new individual sources like a pulsar which shows similar subpulse modulation in both the main- and interpulse and six pulsars with opposite drift senses in different components.
We present the XMM-Newton observations obtained during four revolutions in Spring and Summer 2004 of CXOGC J174540.0-290031, a moderately bright transient X-ray source, located at only 2.9" from SgrA*. We report the discovery of sharp and deep X-ray eclipses, with a period of 27,961+/-5 s and a duration of about 1,100+/-100 s, observed during the two consecutive XMM revolutions from August 31 to September 2. No deep eclipses were present during the two consecutive XMM revolutions from March 28 to April 1, 2004. The spectra during all four observations are well described with an absorbed power law continuum. While our fits on the power law index over the four observations yield values that are consistent with Gamma=1.6-2.0, there appears to be a significant increase in the column density during the Summer 2004 observations, i.e. the period during which the eclipses are detected. The intrinsic luminosity in the 2-10 keV energy range is almost constant with 1.8-2.3 x 10^34 (d_8kpc)^2 erg/s over the four observations. In the framework of eclipsing semidetached binary systems, we show that the eclipse period constrains the mass of the assumed main-sequence secondary star to less than 1.0 M_odot. Therefore, we deduce that CXOGC J174540.0-290031 is a low-mass X-ray binary (LMXB). Moreover the eclipse duration constrains the mass of the compact object to less than about 60 M_odot, which is consistent with a stellar mass black hole or a neutron star. The absence of deep X-ray eclipses during the Spring 2004 observations could be explained if the centroid of the X-ray emitting region moves from a position on the orbital plane to a point above the compact object, possibly coincident with the base of the jet which was detected in radio at this epoch. [Abstract truncated].
Distance estimates based on low-resolution spectroscopy and Two Micron All Sky Survey (2MASS) J magnitudes are presented for a large sample of 322 nearby candidates from Luyten's NLTT catalogue. Mainly relatively bright (typically 7 < K_s < 11) and red high proper motion stars have been selected according to their 2MASS magnitudes and optical-to-infrared colours (+1 < R-K_s < +7). Some LHS stars previously lacking spectroscopy have also been included. We have classified the majority of the objects as early-M dwarfs (M2-M5). More than 70% of our targets turned out to lie within the 25 pc horizon of the catalogue of nearby stars, with 50 objects placed within 15 pc and 8 objects being closer than 10 pc. Three objects in the 10 pc sample have no previously published spectral type: LP 876-10 (M4), LP 870-65 (M4.5), and LP 869-26 (M5). A large fraction of the objects in our sample (57%) have independent distance estimates, mainly by the recent efforts of Reid and collaborators. Our distance determinations are generally in good agreement with theirs. 11 rather distant (d > 100 pc) objects have also been identified, including a probable halo, but relatively hot (T_eff = 13000 K) white dwarf (LHS 1200) and 10 red dwarfs with extremely large tangential velocities (250 < v_t < 1150 km/s). Altogether, there are 11 red dwarfs (including one within 70 pc) with tangential velocities larger than about 250 km/s. All these objects are suspected to be in fact subdwarfs, if so, their distances would be only about half of our original estimates. The three most extreme objects in that respect are the K and early M dwarfs LP 323-168, LHS 5343 and LP 552-21 with corrected distances between 180 pc and 400 pc and resulting tangential velocities still larger than about 400 km/s.
We present interferometric imaging at 33 GHz of the Corona Borealis supercluster, using the extended configuration of the Very Small Array. A total area of 24 deg^2 has been imaged, with an angular resolution of 11 arcmin and a sensitivity of 12 mJy/beam. The aim of these observations is to search for Sunyaev-Zel'dovich (SZ) detections from known clusters of galaxies in this supercluster and for a possible extended SZ decrement due to diffuse warm/hot gas in the intercluster medium. We measure negative flux values in the positions of the ten richest clusters in the region. Collectively, this implies a 3.0-sigma detection of the SZ effect. In the clusters A2061 and A2065 we find decrements of approximately 2-sigma. Our main result is the detection of two strong and resolved negative features at -70+-12 mJy/beam (-157+-27 microK) and -103+-10 mJy/beam (-230+-23 microK), respectively, located in a region with no known clusters, near the centre of the supercluster. We discuss their possible origins in terms of primordial CMB anisotropies and/or SZ signals related to either unknown clusters or to a diffuse extended warm/hot gas distribution. Our analyses have revealed that a primordial CMB fluctuation is a plausible explanation for the weaker feature (probability of 37.82%). For the stronger one, neither primordial CMB (probability of 0.33%) nor SZ can account alone for its size and total intensity. The most reasonable explanation, then, is a combination of both primordial CMB and SZ signal. Finally, we explore what characteristics would be required for a filamentary structure consisting of warm/hot diffuse gas in order to produce a significant contribution to such a spot taking into account the constraints set by X-ray data.
We present results of a microlensing survey toward the Andromeda Galaxy (M31) carried out during four observing seasons at the Isaac Newton Telescope (INT). This survey is part of the larger microlensing survey toward M31 performed by the Microlensing Exploration of the Galaxy and Andromeda (MEGA) collaboration. Using a fully automated search algorithm, we indentify 14 candidate microlensing events, three of which are reported here for the first time. Observations obtained at the Mayall telescope are combined with the INT data to produce composite lightcurves for these candidates. The results from the survey are compared with theoretical predictions for the number and distribution of events. These predictions are based on a Monte Carlo calculation of the detection efficiency and disk-bulge-halo models for M31. The models provide the full phase-space distribution functions for the lens and source populations and are motivated by dynamical and observational considerations. They include differential extinction and span a wide range of parameter space characterized primarily by the mass-to-light ratios for the disk and bulge. For most models, the observed event rate is consistent with the rate predicted for self-lensing -- a MACHO halo fraction of 30% or higher can be ruled at the 95% confidence level. The event distribution does show a large near-far asymmetry hinting at a halo contribution to the microlensing signal. Two candidate events are located at particularly large projected radii on the far side of the disk. These events are difficult to explain by self lensing and only somewhat easier to explain by MACHO lensing. A possibility is that one of these is due to a lens in a giant stellar stream.
We investigate HST/ACS and WFPC2 images at the positions of five candidate microlensing events from a large survey of variability in M31 (MEGA). Three closely match unresolved sources, and two produce only flux upper limits. All are confined to regions of the color-magnitude diagram where stellar variability is unlikely to be easily confused with microlensing. Red variable stars cannot explain these events (although background supernova are possible for two). If these lenses arise in M31's halo, they are due to masses 0.08 < m / M_Sun < 0.85 (95% certainty, for a delta-function mass distribution), brown dwarfs for disk/disk, and stellar masses for disk/bulge ``self-lensing''.
We present optical spectra of the peculiar Type Ia supernova (SN Ia) 1999ac. The data extend from -15 to +42 days with respect to B-band maximum and reveal an event that is unusual in several respects. Prior to B-band maximum, the spectra resemble those of SN 1999aa, a slowly declining event, but possess stronger SiII and CaII signatures (more characteristic of a spectroscopically normal SN). Spectra after B-band maximum appear more normal. The expansion velocities inferred from the Iron lines appear to be lower than average; whereas, the expansion velocity inferred from Calcium H and K are higher than average. The expansion velocities inferred from SiII are among the slowest ever observed, though SN 1999ac is not particularly dim. The analysis of the parameters v_10, R(SiII), dv(SiII)/dt, and d_m15 further underlines the unique characteristics of SN 1999ac. We find convincing evidence of CII 6580 in the day -15 spectrum with ejection velocity v > 16,000 km/s, but this signature disappears by day -9. This rapid evolution at early times highlights the importance of extremely early-time spectroscopy.
We investigate the point-parity and mirror-parity handedness of the large angle anisotropy in the cosmic microwave background (CMB). In particular we consider whether the observed low CMB quadrupole could more generally signal odd point-parity, i.e. suppression of even multipoles. Even though this feature is ``visually'' present in most renditions of the WMAP dataset we find that it never supports parity preference beyond the meagre 95% confidence level. This is fortunate as point parity handedness implies almost certainly a high level of galactic contamination. Mirror reflection parity, on the contrary, is related to the emergence of a preferred axis, defining the symmetry plane. We use this technique to make contact with recent claims for an anisotropic Universe, showing that the detected preferred axis is associated with positive (even) mirror parity. This feature may be an important clue in identifying the culprit for this unexpected signal.
The formalism in order to obtain the Dark Energy equation of state is extended to non-flat universes and we consider the inequalities that must be satisfied by Phantom Dark Energy in this case. We show that due to a non-vanishing spatial curvature satisfying the observational bounds, the uncertainty on the determination of the Dark Energy equation of state parameter $w$, when it is taken constant, can be significant and that it is minimal for some redshift $z_{cr}\sim 3$. We consider the potential of future measurements of the gravitational waves emitted by binaries at high redshifts $z>z_{cr}$ to reduce this uncertainty. Results obtained here should also be relevant for a weakly varying equation of state with $w\approx -1$.
The line emission from a growing number of Herbig-Haro jets can be observed and resolved at angular distances smaller than a few arcseconds from the central source. The interpretation of this emission is problematic, since the simplest model of a cooling jet cannot sustain it. It has been suggested that what one actually observes are shocked regions with a filling factor of $\sim 1%$. In this framework, up to now, comparisons with observations have been based on stationary shock models. Here we introduce for the first time the self-consistent dynamics of such shocks and we show that considering their properties at different times, i.e. locations, we can reproduce observational data of the DG Tau microjet. In particular, we can interpret the spatial behavior of the [SII]6716/6731 and [NII]/[OI]6583/6300 line intensity ratios adopting a set of physical parameters that yield values of mass loss rates and magnetic fields consistent with previous estimates. We also obtain the values of the mean ionization fraction and electron density along the jet, compare these values with the ones derived from observations using the sulfur doublet to constrain the electron density (e.g. Bacciotti et al. 1995).
Solid particles with the composition of interstellar dust and power-law size distribution dn/da propto a^{-p} for a < a_max with a_{max} > 3 lambda and 3 < p < 4 will have submm opacity spectral index beta(lambda) = dln(kappa)/dln(nu) approx (p-3) beta_{ism}, where beta_{ism} approx 1.7 is the opacity spectral index of interstellar dust material in the Rayleigh limit. For the power-law index p approx 3.5 that characterizes interstellar dust, and that appears likely for particles growing by agglomeration in protoplanetary disks, grain growth to sizes a > 3 mm will result in beta(1 mm) < ~1. Grain growth can naturally account for beta approx 1 observed for protoplanetary disks, provided that a_{max} > ~ 3 lambda.
Greisen & Calabretta describe a generalized method for specifying the coordinates of FITS data samples. Following that general method, Calabretta & Greisen describe detailed conventions for defining celestial coordinates as they are projected onto a two-dimensional plane. The present paper extends the discussion to the spectral coordinates of wavelength, frequency, and velocity. World coordinate functions are defined for spectral axes sampled linearly in wavelength, frequency, or velocity, linearly in the logarithm of wavelength or frequency, as projected by ideal dispersing elements, and as specified by a lookup table.
Understanding the explosion mechanism of core collapse supernovae is a problem that has plagued nuclear astrophysicists since the first computational models of this phenomenon were carried out in the 1960s. Our current theories of this violent phenomenon center around multi-dimensional effects involving radiation-hydrodynamic flows of hot, dense matter and neutrinos. Modeling these multi-dimensional radiative flows presents a computational challenge that will continue to stress high-performance computing beyond the teraflop to the petaflop level. In this paper we describe a few of the scientific discoveries that we have made via terascale computational simulations of supernovae under the auspices of the SciDAC-funded Terascale Supernova Initiative.
We present results from a survey for gravitationally lensed radio lobes. Lensed lobes are a potentially richer source of information about galaxy mass distributions than lensed point sources, which have been the exclusive focus of other recent surveys. Our approach is to identify radio lobes in the FIRST catalog and then search optical catalogs for coincident foreground galaxies, which are candidate lensing galaxies. We then obtain higher-resolution images of these targets at both optical and radio wavelengths, and obtain optical spectra for the most promising candidates. We present maps of several radio lobes that are nearly coincident with galaxies. We have not found any new and unambiguous cases of gravitational lensing. One radio lobe in particular, FOV J0743+1553, has two hot spots that could be multiple images produced by a z=0.19 spiral galaxy, but the lensing interpretation is problematic.
We study the formation of disks via the cooling flow of gas within galactic haloes using smoothed particle hydrodynamics simulations. These simulations resolve mass scales of a few thousand solar masses in the gas component for the first time. Thermal instabilities result in the formation of numerous warm clouds that are pressure confined by the hot ambient halo gas. The clouds fall slowly onto the disk through non-spherical accretion from material flowing preferentially down the angular momentum axis. The rotational velocity of the infalling cold gas decreases as a function of height above the disk, closely resembling that of the extra-planar gas recently observed around the spiral galaxy NGC 891.
We propose that the anomalously bright white dwarf luminosity function observed in NGC 6791 (Bedin et al 2005) is the consequence of the formation of 0.5 Msun white dwarfs with Helium cores instead of Carbon cores. This may happen if mass loss during the ascent of the Red Giant Branch is strong enough to prevent a star from reaching the Helium flash. Such a model can explain the slower white dwarf cooling (relative to standard models) and fits naturally with scenarios advanced to explain Extreme Horizontal Branch stars, a population of which are also found in this cluster.
The discovery of the Cosmic Infrared Background (CIB) in 1996, together with recent cosmological surveys from the mid-infrared to the millimeter have revolutionized our view of star formation at high redshifts. It has become clear, in the last decade, that a population of galaxies that radiate most of their power in the far-infrared (the so-called ``infrared galaxies'') contributes an important part of the whole galaxy build-up in the Universe. Since 1996, detailed (and often painful) investigations of the high-redshift infrared galaxies have resulted in the spectacular progress covered in this review. We outline the nature of the sources of the CIB including their star-formation rate, stellar and total mass, morphology, metallicity and clustering properties. We discuss their contribution to the stellar content of the Universe and their origin in the framework of the hierarchical growth of structures. We finally discuss open questions for a scenario of their evolution up to the present-day galaxies.
We use a Lucky Imaging system to obtain I-band images with much improved angular resolution on a ground-based 2.5m telescope. We present results from a 10-night assessment campaign on the 2.56m Nordic Optical Telescope and quantify the performance of our system in seeings better than 1.0''. In good seeing we have acquired near diffraction-limited images; in poorer seeing the angular resolution has been routinely improved by factors of 2.5-4. The system can use guide stars as faint as I=16 with full performance and its useful field of view is consistently larger than 40" diameter. The technique shows promise for a number of science programmes, both galactic (eg. binary candidates, brown dwarfs, globular cluster cores) and extragalactic (eg. quasar host galaxies, damped Lyman-alpha absorbers).
Two leading dark matter candidates from supersymmetry and other theories of physics beyond the standard model are WIMPs and weak scale gravitinos. If the lightest stable particle is a gravitino, then a WIMP will decay into it with a natural lifetime of order a month ~ M_{pl}^2/M_{weak}^3. We show that if the bulk of dark matter today came from decays of neutral particles with lifetimes of order a year or smaller, then it could lead to a reduction in the amount of small scale substructure, less concentrated halos and constant density cores in the smallest mass halos. Such beneficial effects may therefore be realized naturally, as discussed by Cembranos, Feng, Rajaraman, and Takayama, in the case of supersymmetry.
Cosmic shear offers a unique opportunity to probe the matter power spectrum directly. Future wide-angle surveys should strongly constrain a number of cosmological parameters. However, conventional wisdom states that the lensing kernel associated with cosmic shear is too broad to detect the baryon oscillations. These features originated prior to recombination, induced by the acoustics of the photon-baryon fluid. Recent galaxy surveys have detected these imprints, and in the future such measurements may even be used to refine our understanding of dark energy. With an angular-dependent selection of photometric redshift bins, we demonstrate that cosmic shear is indeed sensitive to the presence of the baryon oscillations. As an illustration, we show that data from next-generation surveys will be able to rule out a smoothed power spectrum at the two sigma level.
In Sgr A* at the Galactic center, by far the closest and easiest supermassive black hole we can study, the observational evidence is increasingly pointing to the presence of a compact, hot, magnetized disk feeding the accretor. In such low-Mach-number plasmas, forces arising, e.g., from pressure gradients in the plasma, can altogether negate the warping of disks around Kerr black holes caused by the Bardeen-Petterson effect and can lead to coherent precession of the entire disk. In this Letter, we present for the first time highly detailed 3D SPH simulations of the accretion disk evolution in Sgr A*, guided by observational constraints on its physical characteristics, and conclude that indeed the Bardeen-Petterson effect is probably absent in this source. Given what we now understand regarding the emission geometry in this object, we suggest that a ~ 50-500-day modulation in Sgr A*'s spectrum, arising from the disk precession, could be an important observational signature; perhaps the ~ 106-day period seen earlier in its radio flux, if confirmed, could be due to this process. On the other hand, if future observations do not confirm this long modulation in Sgr A*'s spectrum, this would be an indication that either the disk size or orientation is very different from current estimates, or that the black hole is not spinning at all (unlikely), or that our current understanding of how it produces its radiative output is incorrect.
We present results of a study aiming at shedding light on the specific advantages and limitations of methods to derive star formation histories (SFH) in galaxies using resolved stellar populations and using integrated light, respectively. For this purpose, we analyse the integrated-light spectrum of a field in the LMC bar, for which highly resolved HST images are available as well. To be compared with the SFH derived from the color magnitude diagram (CMD) of this field (Smecker-Hane et al. 02), we have performed a set of simulations of galaxies with systematically varying SFHs, but constant metallicity (Z=0.008). We investigate to which extent different SF scenarios can be discriminated on the basis of their photometric and spectral properties, respectively, and determine in how far the detailed SFH obtained by the CMD approach can be reproduced by results based upon integrated properties. Comparing both methods for this nearby field we want to learn about the specific character of both methods and understand to what precision SFHs can be determined for distant galaxies only observable in integrated light. All simulations are performed using our evolutionary synthesis code GALEV.
Cygnus X-3 is a strong X-ray source (L_X about 10^38 erg/s) which is thought to consist of a compact object, accreting matter from a helium star. We find analytically that the estimated ranges of mass-loss rate and orbital-period derivative for Cyg X-3 are consistent with two models: i) the system is detached and the mass loss from the system comes from the stellar wind of a massive helium star, of which only a fraction that allows for the observed X-ray luminosity is accreted, or ii) the system is semidetached and a Roche-lobe-overflowing low- or moderate-mass helium donor transfers mass to the compact object, followed by ejection of its excess over the Eddington rate from the system. These analytical results appear to be consistent with evolutionary calculations. By means of population synthesis we find that currently in the Galaxy there may exist ~1 X-ray binary with a black hole that accretes from a >~ 7 MSun Wolf-Rayet star and ~1 X-ray binary in which a neutron star accretes matter from a Roche-lobe-overflowing helium star with mass <~ 1.5 MSun. Cyg X-3 is probably one of these systems.
We present results of a detailed study aiming at understanding to what precision star formation histories (SFHs) can be determined for distant galaxies observable in integrated light only. Using our evolutionary synthesis code, we have performed a set of simulations of galaxies with a wide range of different SFHs. By analysing the resulting colors, spectra and Lick indices, we investigate to which extent different SF scenarios can be discriminated on the basis of their photometric and spectral properties, respectively. We find the robust result that no later than 4 Gyrs after the latest episode of enhanced star formation all scenarios exhibit very similar colors and indices; in practice, it is not possible to distinguish different scenarios of star formation which have evolved for more than 1, at the utmost 3-4 Gyrs since the last star forming event, even when using spectral indices. For how long different SF scenarios can be disentangled highly depends on the range of colors available and absorption lines considered, as well as on the details of the SFHs to be compared.
We have obtained a high S/N (22.3 hr integration) UV continuum VLT FORS2 spectrum of an extremely bright (z_850 = 24.3) z = 5.515 +/- 0.003 starforming galaxy (BD38) in the field of the z = 1.24 cluster RDCS 1252.9-2927. This object shows substantial continuum (0.41 +/- 0.02 \muJy at \lambda1300) and low-ionization interstellar absorption features typical of LBGs at lower redshift (z ~ 3); this is the highest redshift LBG confirmed via metal absorption spectral features. The equivalent widths of the absorption features are similar to z ~ 3 strong Ly\alpha absorbers. No noticeable Ly\alpha emission was detected (F <= 1.4 * 10^-18 ergs cm^-2 s^-1, 3\sigma). The half-light radius of this object is 1.6 kpc (0\farcs25) and the star formation rate derived from the rest-frame UV luminosity is SFR_UV = 38 h^-2_0.7 M_sun yr^-1 (142 h^-2_0.7 M_sun yr^-1 corrected for dust extinction). In terms of recent determinations of the z ~ 6 UV luminosity function, this object appears to be 6L*. The Spitzer IRAC fluxes for this object are 23.3 and 23.2 AB mag (corrected for 0.3 mag of cluster lensing) in the 3.6\mu and 4.5\mu channels, respectively, implying a mass of 1-6 * 10^10 M_sun from population synthesis models. This galaxy is brighter than any confirmed z ~ 6 i-dropout to date in the z_850 band, and both the 3.6\mu and 4.5\mu channels, and is the most massive starbursting galaxy known at z > 5. -- Abstract Abridged
The coupling of photons and baryons by Thomson scattering in the early universe imprints features in both the Cosmic Microwave Background (CMB) and matter power spectra. The former have been used to constrain a host of cosmological parameters, the latter have the potential to strongly constrain the expansion history of the universe and dark energy. Key to this program is the means to localize the primordial features in observations of galaxy spectra which necessarily involve galaxy bias, non-linear evolution and redshift space distortions. We present calculations, based on mock catalogues produced from particle-mesh simulations, which show the range of behaviors we might expect of galaxies in the real universe. We find that non-linearity, galaxy bias and redshift space distortions all introduce important modifications to the basic picture. Both the galaxy bias and the (isotropic) redshift space distortions lead to relatively smooth modulations in power on the scales of interest to baryon oscillations. The halo and galaxy power spectra exhibit low order structure beyond the acoustic oscillations over the range of scales of relevance. Fitting a cubic polynomial to the ratio of galaxy to dark matter power reduces any remaining structure in the range 0.03<k<0.3 h/Mpc below the 2 per cent level, which is close to the error in our calculations.
Baryon acoustic oscillations, measured through the patterned distribution of galaxies or other baryon tracing objects on very large (100 Mpc) scales, offer a possible geometric probe of cosmological distances. Pluses and minuses in this approach's leverage for understanding dark energy are discussed, as are systematic uncertainties requiring further investigation. Conclusions are that 1) BAO offer promise of a new avenue to distance measurements and further study is warranted, 2) the measurements will need to attain ~1% accuracy (requiring a 10000 square degree spectroscopic survey) for their dark energy leverage to match that from supernovae, but do give complementary information at 2% accuracy. Because of the ties to the matter dominated era, BAO is not a replacement probe of dark energy, but a valuable complement.
In this work we show that Bardeen-Petterson accretion disks can exhibit unique, detectable features in relativistically broadened emission line profiles. Some of the unique characteristics include inverted line profiles with sharper red horns and softer blue horns and even profiles with more than 2 horns from a single rest-frame line. We demonstrate these points by constructing a series of synthetic line profiles using simple two-component disk models. We find that the resultant profiles are very sensitive to the two key parameters one would like to constrain, namely the Bardeen-Petterson transition radius r_{BP} and the relative tilt \beta between the two disk components over a range of likely values [10 < r_{BP}/(GM/c^2) < 40 ; 15deg < \beta < 45deg]. We use our findings to show that some of the ``extra'' line features observed in the spectrum of the Seyfert-I galaxy MCG--6-30-15 may be attributable to a Bardeen-Petterson disk structure. Similarly, we apply our findings to two likely Bardeen-Petterson candidate Galactic black holes - GRO J1655-40 and XTE J1550-564. We provide synthetic line profiles of these systems using observationally constrained sets of parameters. Although we do not formally fit the data for any of these systems, we confirm that our synthetic spectra are consistent with current observations.
This is the fifth paper in a series studying the stellar components, star formation histories, star formation rates and metallicities of a blue compact galaxy (BCG) sample. Based on our high-quality ground-based spectroscopic observations, we have determined the electron temperatures, electron densities, nitrogen abundances and oxygen abundances for 72 star-forming BCGs in our sample, using different oxygen abundance indicators. The oxygen abundance covers the range 7.15 < 12 + log (O/H)< 9.0, and nitrogen is found to be mostly a product of secondary nucleosynthesis for 12 + log (O/H)>8.2 and apparently a product of primary nucleosynthesis for 12 + log (O/H)< 8.2. To assess the possible systematic differences among different oxygen abundance indicators, we have compared oxygen abundances of BCGs obtained with the Te method, R23 method, P method, N2 method and O3N2 method. The oxygen abundances derived from the Te method are systematically lower by 0.1-0.25 dex than those derived from the strong line empirical abundance indicators, consistent with previous studies based on region samples. We confirm the existence of the metallicity-luminosity relation in BCGs over a large range of abundances and luminosities. Our sample of galaxies shows that the slope of the metallicity-luminosity relation for the luminous galaxies (~-0.05) is slightly shallower than that for the dwarf galaxies (~-0.17). An offset was found in the metallicity-luminosity relation of the local galaxies and that of the intermediate redshift galaxies. It shows that the metallicity-luminosity relation for the emission line galaxies at high redshift is displaced to lower abundances, higher luminosities, or both.
The reaction rates of the beta processes for all particles at arbitrary degeneracy are derived, and an {\it analytic} steady state equilibrium condition $\mu_n=\mu_p+2\mu_e$ which results from the equality of electron and positron capture rates in the hot electron-positron plasma with nucleons is also found, if the matter is transparent to neutrinos. This simple analytic formula is valid only if electrons are nondegenerate or mildly degenerate, which is generally satisfied in the hot electron-positron plasma. Therefore, it can be used to efficiently determine the steady state of the hot matter with plenty of positrons. Based on this analytic condition, given the baryon number density and the temperature, if the nucleons are nondegenerate, only one algebraic equation for determining the electron fraction is obtained, which shows the great advantage of the analytic equilibrium condition.
The possibility that neutron stars may contain substantial hyperon populations has important implications for neutron star cooling and, through bulk viscosity, the viability of the r-modes of accreting neutron stars as sources of persistent gravitational waves. In conjunction with laboratory measurements of hypernuclei, astronomical observations were used by Glendenning and Moszkowski [Phys. Rev. Lett. 67, 2414 (1991)] to constrain the properties of hyperonic equations of state within the framework of relativistic mean field theory. We revisit the problem, incorporating recent measurements of a high neutron star mass and a gravitational redshift. We find that only the stiffest of the relativistic hyperonic equations of state commonly used in the literature is compatible with the redshift. However, it is possible to construct a stiffer equation of state within the same framework which produces the observed redshift while satisfying the experimental constraints on hypernuclei. Nonrelativistic potential-based equations of state with hyperons are not constrained by the redshift, primarily due to a smaller stellar radius.
The linear stability theory of Taylor-Couette flows (unbounded in_z_) is described including magnetic fields, Hall effect or a density stratification in order to prepare laboratory experiments to probe the stability of differential rotation in astrophysics. If an axial field is present then the magnetorotational instability (MRI) is investigated also for small magnetic Prandtl numbers. For rotating outer cylinder beyond the Rayleigh line characteristic minima are found for magnetic Reynolds number of the order of 10 and for Lundquist numbers of order 1. The inclusion of extra toroidal current-free fields leads to new oscillating solutions with rather low Reynolds numbers and Hartmann numbers. The Hall effect establishes an unexpected `shear-Hall instability' (SHI) where shear and magnetic field have opposite signs. In this case even rotation laws increasing outwards may become unstable. Recently global solutions have been found for the Taylor-Couette flows with density stratified in axial directions (`stratorotational instability', SRI). They exist beyond the Rayleigh line for Froude numbers of moderate order but also not for too flat rotation laws with \hat\eta^2 < \hat\mu <\hat \eta.
With regard to homeoidally striated Jacobi ellipsoids, a unified theory of systematically rotating and peculiar motions is developed, where both real and imaginary rotation are considered. The effect of positive or negative residual motion excess, is shown to be equivalent to an additional real or imaginary rotation, respectively. Then it is realized that a homeoidally striated Jacobi ellipsoid with assigned velocity distribution, always admits an adjoint configuration i.e. a classical Jacobi ellipsoid of equal mass and axes. In addition, further constraints are established on the amount of residual velocity anisotropy along the principal axes, for triaxial configurations. Special effort is devoted to investigating sequences of virial equilibrium configurations in terms of normalized parameters, including the effects of both density and velocity profile. In particular, it is shown that bifurcation points from axisymmetric to triaxial configurations occur as in classical Jacobi ellipsoids, contrary to earlier results. The reasons of the above mentioned discrepancy are also explained. An interpretation of recent results from numerical simulations on stability (Meza 2002), is provided in the light of the model. A physical interpretation of the early Hubble sequence is shortly reviewed and discussed from the standpoint of the model, according if elliptical galaxies are considered as isolated systems or embedded within dark matter haloes. In any case, a lower limit to the flattening of oblate-like configurations is established. On the other hand, it is found no lower limit to the elongation of prolate-like configurations, and the existence of some sort of instability is predicted, owing to the observed lack of elliptical galaxies more flattened or elongated than E7.
High-velocity galactic outflows, driven by intense bursts of star formation and black hole accretion, are invoked by current theories of galaxy formation to terminate star formation in the most massive galaxies and to deposit heavy elements in the intergalactic medium. From existing observational evidence on high-redshift galaxies, it is unclear whether such outflows are localized to regions of intense star formation just a few kiloparsecs in extent, or whether they instead have a significant impact on the entire galaxy and its surroundings. Here we present two-dimensional spectroscopy of a star-forming galaxy at redshift z=3.09 (seen 11.5 Gyr ago, when the Universe was 20 per cent of its current age): its spatially extended Ly-alpha emission appears to be absorbed by HI in a foreground screen covering the entire galaxy, with a lateral extent of at least 100 kpc and remarkable velocity coherence. It was plausibly ejected from the galaxy during a starburst several 1E8 yr earlier and has subsequently swept up gas from the surrounding intergalactic medium and cooled. This demonstrates the galaxy-wide impact of high-redshift superwinds.
We discuss the contribution of kiloparsec-scale jets in FR I radio galaxies to the diffuse gamma-ray background radiation. The analyzed gamma-ray emission is due inverse-Compton scattering of starlight photon fields by the ultrarelativistic electrons whose synchrotron radiation is detected from such sources at radio, optical and X-ray energies. We find that these objects, under the minimum-power hypothesis (corresponding to a magnetic field of 300 micro G in the brightest knots of these jets), can contribute about one percent to the extragalactic gamma-ray background measured by EGRET. We point out that this result already indicates that the magnetic fields in kpc-scale jets of low-power radio galaxies are not likely to be smaller than 10 micro G on average, as otherwise the extragalactic gamma-ray background would be overproduced.
The biologically damaging solar ultraviolet (UV) radiation (quantified by the DNA-weighted dose) reaches the Martian surface in extremely high levels. Searching for potentially habitable UV-protected environments on Mars, we considered the polar ice caps that consist of a seasonally varying CO2 ice cover and a permanent H2O ice layer. It was found that, though the CO2 ice is insufficient by itself to screen the UV radiation, at 1 m depth within the perennial H2O ice the DNA-weighted dose is reduced to terrestrial levels. This depth depends strongly on the optical properties ofthe H2O ice layers (for instance snow-lile layes). The Earth-like DNA-weighted dose and Photosynthetically Active Radiation (PAR) requirements were used to define the upper and lower limits of the nortern and southern polar radiative habitable zone (RHZ) for which a temporal and spatial mapping was performed. Based on these studies we conclude that phtosynthetic life might be possible within the ice layers of the polar regions. The thickness varies along each Martian polar spring and summer between 1.5 m and 2.4 m for H2= ice-like layers, and a few centimeters for snow-like covers. These Martian Earth-like radiative habitable environments may be primary targets for future Martian astrobiological missions. Special attention should be paid to planetary protection, since the polar RHZ may also be subject to terrestrial contamination by probes.
We apply the thermodynamical model of the cosmological event horizon of the spatially flat FLRW metrics to the study of the recent accelerated expansion phase and to the coincidence problem. This model, called "ehT model" hereafter, led to a dark energy (DE) density $\Lambda $ varying as $r^{-2},$ where $r$ is the proper radius of the event horizon. Recently, another model motivated by the holographic principle gave an independent justification of the same relation between $\Lambda $ and $r$. We probe the theoretical results of the ehT model with respect to the SnIa observations and we compare it to the model deduced from the holographic principle, which we call "LHG model" in the following.Our results are in excellent agreement with the observations for $H\_{0}=64kms^{-1}Mpc^{-1}$, and $\Omega \_{\Lambda }^{0}=0.63\_{-0.01}^{+0.1}$, which leads to $q\_{0}=-0.445$ and $z\_{T}\simeq 0.965$.
Nbody simulations are used to examine the consequences of Neptune's outward migration into the Kuiper Belt, with the simulated endstates being compared rigorously and quantitatively to the observations. These simulations confirm the findings of Chiang et al. (2003), who showed that Neptune's migration into a previously stirred-up Kuiper Belt can account for the Kuiper Belt Objects (KBOs) known to librate at Neptune's 5:2 resonance. We also find that capture is possible at many other weak, high-order mean motion resonances, such as the 11:6, 13:7, 13:6, 9:4, 7:3, 12:5, 8:3, 3:1, 7:2, and the 4:1. The more distant of these resonances, such as the 9:4, 7:3, 5:2, and the 3:1, can also capture particles in stable, eccentric orbits beyond 50 AU, in the region of phase space conventionally known as the Scattered Disk. Indeed, 90% of the simulated particles that persist over the age of the Solar System in the so-called Scattered Disk zone never had a close encounter with Neptune, but instead were promoted into these eccentric orbits by Neptune's resonances during the migration epoch. This indicates that the observed Scattered Disk might not be so scattered. This model also produced only a handful of Centaurs, all of which originated at Neptune's mean motion resonances in the Kuiper Belt. We also report estimates of the abundances and masses of the Belt's various subpopulations (e.g., the resonant KBOs, the Main Belt, and the so-called Scattered Disk), and also provide upper limits on the abundance of Centaurs and Neptune's Trojans, as well as upper limits on the sizes and abundances of hypothetical KBOs that might inhabit the a>50 AU zone.
We combined the ISOCAM Parallel Mode Survey at 6.7 micron (LW2 filter) with the Two Micron All Sky Survey in order to obtain a powerful tool to search for AGN independent of dust extinction. Using moderate colour criteria H-K>0.5 and K-LW2>2.7 we have selected a sample of 77 AGN candidates in an effective area of about 10 square degrees. By means of optical spectroscopy we find 24 (30%) type-1 QSOs at redshifts 0.1<z<2.3; nine of them have z>0.8. About one third of the ISO-2MASS QSOs show so red optical colours, that they are missed in optical and UV AGN surveys like SDSS, 2DF, or HES. With a surface density of about 2 deg^(-2) down to R<18 mag the ISO-2MASS QSOs outnumber the 1.35 deg^(-2) of the SDSS quasar survey by 50%; we find a combined optical-IR QSO surface density of 2.7 deg^(-2). Since only two of the ISO-2MASS QSOs have also J-K>2, the inclusion of the ISO mid-infrared photometry significantly extends the capabilities of the pure 2MASS red AGN survey. We suggest that the newly found red AGN resemble young members of the quasar population, and that quasars spend much of their lifetime in a dust enshrouded phase.
We present a 23 ks XMM-Newton observation of the ULIRG Mrk 273. The hard X-ray spectrum can be modeled by a highly absorbed power law plus an Fe Kalpha emission line. The iron line is broad, suggesting possible superposition of a neutral iron line at 6.4 keV, and a blend of ionized iron lines from Fe XXV and Fe XXVI. Given the relatively short exposure, the three line components can not be singularly resolved with high statistical significance: the neutral component is detected at ~2.5sigma and the Fe XXV line at ~2sigma c.l., while for the Fe XXVI line we can only estimate an upper limit. The broad band spectrum requires, in addition to a highly absorbed power law, at least three collisionally ionized plasma components, which may be associated with star-forming regions. The temperatures of the three plasmas are about 0.3, 0.8 and 6 keV, where the highest of the three is sufficient to produce ionized iron emission lines. An alternative interpretation for the origin of the soft emission might also be given in terms of reflection off some photoionized gas, as has been observed in a number of nearby Compton-thick Seyfert 2 galaxies. A hot gas, photoionized by the primary, continuum can also produce ionized iron lines. Unfortunately, given the limited statistics and the lack of high resolution spectroscopy, it is not possible to distinguish between the two models investigated. We further compare the XMM-Newton findings with the Chandra data obtaining consistent spectral results. The absorption corrected hard X-ray (2-10 keV) luminosity of Mrk 273 corresponds to ~0.2% of the far-IR luminosity, similar to typical values found in pure starbursts. The thermal contribution to the soft X-ray luminosity is ~0.2-0.7x10^42 erg s^-1, comparable to those found in NGC 6240 and other starburst dominated ULIRGs.
In cosmology many dramatically different scenarios with the past (big bang versus bounce) and in the future (de Sitter versus big rip) singularities are compatible with the present day observations. This difficulty is called the degeneracy problem. We use the Akaike and Bayesian information criteria of model selection to overcome this degeneracy and to determine a model with such a set of parameters which gives the most preferred fit to the SNIa data. We consider seven representative scenarios, namely: the CDM models with the cosmological constant, with topological defect, with phantom field, with bounce, with bouncing phantom field, with brane and model with the linear dynamical equation of state parameter. Applying the model selection information criteria we show that AIC indicates the flat phantom model while BIC indicates both flat phantom and flat $\Lambda$CDM models. Finally we conclude that the number of essential parameters chosen by dark energy models which are compared with SNIa data is two.
We observed with XMM-Newton 4 objects selected from the Grossan catalog, with the aim to search for new 'changing-look' AGN. The sample includes all the sources which showed in subsequent observations a flux much lower than the one measured with HEAO A-1: NGC 7674, NGC 4968, IRAS 13218+0552 and NGC 1667. None of the sources was caught in a high flux state during the XMM-Newton observations, whose analysis reveal they are all likely Compton-thick objects. We suggest that, for all the sources, potential problems with the HEAO A-1 source identification and flux measurement prevent us from being certain that the HEAO A-1 data represent a putative 'high' state for these objects. Nonetheless, based on the high flux state and Compton-thin spectrum of its GINGA observation, NGC 7674 represents probably the sixth known case of a 'changing-look' Seyfert 2 galaxy. From the X-ray variability pattern, we can estimate a likely lower limit of a few parsec to the distance of the inner walls of the torus in this object. Remarkably, IRAS 13218+0552 was not detected by XMM-Newton, despite being currently classified as a Seyfert 1 with a large [OIII] flux. However, the original classification was likely to be affected by an extreme velocity outflow component in the emission lines. The object likely harbors an highly obscured AGN and should be re-classified as a Type 2 source.
We present results of numerical simulations of the propagation of ultra high
energy cosmic rays (UHECRs) over cosmological distances, aimed at quantifying
the statistical significance of the highest energy data on the spectrum and
small scale anisotropies as detected by the AGASA experiment. We assess the
significance of the lack of a GZK feature and its compatibility with the
reported small scale anisotropies.
Assuming that UHECRs are protons from extragalactic sources, we find that the
small scale anisotropies are incompatible with the reported spectrum at a
probability level of $2 \times 10^{-5}$. Our analysis of the AGASA results
shows the power of combining spectrum and small scale anisotropy data in future
high statistics experiments, such as Auger.
In this paper we discuss the anti-GZK effect that arises in the framework of the diffusive propagation of Ultra High Energy (UHE) protons. This effect consists in a jump-like increase of the maximum distance from which UHE protons can reach the observer. The position of the jump is independent of the Intergalactic Magnetic Field (IMF) strength and depends only on the energy losses of protons, namely on the transition energy from adiabatic and pair-production energy losses. The Ultra High Energy Cosmic Rays (UHECR) spectrum presents a low-energy steepening approximately at this energy, which is very close to the position of the observed second knee. The dip, seen in the universal spectrum as a signature of the proton interaction with the Cosmic Microwave Background (CMB) radiation, is also present in the case of diffusive propagation in magnetic fields.
Blue Luminescence (BL) was first discovered in a proto-planetary nebula, the Red Rectangle (RR) surrounding the post-AGB star HD 44179. BL has been attributed to fluorescence by small, 3-4 ringed neutral polycyclic aromatic hydrocarbon (PAH) molecules, and was thought to be unique to the RR environment where such small molecules are actively being produced and shielded from the harsh interstellar radiation by a dense circumstellar disk. In this paper we present the BL spectrum detected in several ordinary reflection nebulae illuminated by stars having temperatures between 10,000 -- 23,000 K. All these nebulae are known to also exhibit the infrared emission features called aromatic emission features (AEFs) attributed to large PAHs. We present the spatial distribution of the BL in these nebulae. In the case of Ced~112, the BL is spatially correlated with mid-IR emission structures attributed to AEFs. These observations provide evidence for grain processing and possibly for in-situ formation of small grains and large molecules from larger aggregates. Most importantly, the detection of BL in these ordinary reflection nebulae suggests that the BL carrier is an ubiquitous component of the ISM and is not restricted to the particular environment of the RR.
By using the method presented by Isobe et al. (2002), the non-dimensional reconnection rate Vin/Va has been determined for the impulsive phase of three two-ribbon flares, where Vin is the velocity of the reconnection inflow and Va is the Alfven velocity. The non-dimensional reconnection rate is important to make a constraint on the theoretical models of magnetic reconnection. In order to reduce the uncertainty of the reconnection rate, it is important to determine the energy release rate of the flares from observational data as accurately as possible. To this end, we have carried out one dimensional hydrodynamic simulations of a flare loop and synthesized the count rate detected by the soft X-ray telescope (SXT) aboard Yohkoh satellite. We found that the time derivative of the thermal energy contents in a flare arcade derived from SXT data is smaller than the real energy release rate by a factor of 0.3 - 0.8, depending on the loop length and the energy release rate. The result of simulation is presented in the paper and used to calculate the reconnection rate. We found that reconnection rate is 0.047 for the X2.3 flare on 2000 November 24, 0.015 for the M3.7 flare on 2000 July 14, and 0.071 for the C8.9 flare on 2000 November 16. These values are similar to that derived from the direct observation of the reconnection inflow by Yokoyama et al. (2001), and consistent with the fast reconnection models such as that of Petschek (1964).
Past hydrodynamic simulations have been able to reproduce the high temperatures and densities characteristic of solar flares. These simulations, however, have not been able to account for the slow decay of the observed flare emission or the absence of blueshifts in high spectral resolution line profiles. Recent work has suggested that modeling a flare as an sequence of independently heated threads instead of as a single loop may resolve the discrepancies between the simulations and observations. In this paper we present a method for computing multi-thread, time-dependent hydrodynamic simulations of solar flares and apply it to observations of the Masuda flare of 1992 January 13. We show that it is possible to reproduce the temporal evolution of high temperature thermal flare plasma observed with the instruments on the \textit{GOES} and \textit{Yohkoh} satellites. The results from these simulations suggest that the heating time-scale for a individual thread is on the order of 200 s. Significantly shorter heating time scales (20 s) lead to very high temperatures and are inconsistent with the emission observed by \textit{Yohkoh}.
With the William Herschel Telescope in La Palma we made IJKs observations of an area of about 40'x30' of the Local Group galaxy Draco. This allows us to describe Draco's late-type stellar population across the whole galaxy at a photometric level 2 mag deeper than the 2MASS survey. We detected the red giant branch (RGB) and measured the magnitude of the tip of the RGB in the three bands. From that in the I band we obtain a distance modulus of (m-M)_0=19.49+/-0.06(stat)+/-0.15(sys), in excellent agreement with a measurement from RR Lyrae stars. The peak of the (J-Ks)_0 histogram at different M_Ks suggests that Draco has a mean [Fe/H]=-1.95+/-1.26 while fiducial RGB tracks of Galactic globular clusters indicate a mean [Fe/H]=-1.33+/-0.72 where the error corresponds to the spread around the mean value. There are significant differences between the colour-magnitude diagrams of stars in the inner, medium and outer areas of the galaxy. A metal poor (Z=0.0004) intermediate-age population (about 1.6 Gyr old) is clearly present and emerges in particular between 6' and 12' from the centre of the galaxy. A few additional carbon star candidates have been identified from both their location in the colour-magnitude diagram and from an indication of variability. The large scale distribution of late-type stars is smooth but irregular in shape; this points at a variation of inclination with radius.
We present optical observations and 2D hydrodynamic modeling of an isolated shocked ISM cloud. H$\alpha$ images taken in 1992.6 and 2003.7 of a small optical emission cloud along the southwestern limb of the Cygnus Loop were used to measure positional displacements of $\sim$ $0 \farcs 1$ yr$^{-1}$ for surrounding Balmer dominated emission filaments and $0\farcs025 - \farcs055$ yr$^{-1}$ for internal cloud emission features. These measurements imply transverse velocities of $\simeq$ 250 km s$^{-1}$ and $\simeq$ 80 -- 140 km s$^{-1}$ for ambient ISM and internal cloud shocks respectively. The complex shock structure visible within the cloud indicates that the cloud's internal density distribution is two phased: a smoothly varying background density which is populated by higher density clumps. We present model results for a shock interacting with a non-uniform ISM cloud. We find that this cloud can be well modeled by a smoothly varying power law core surrounded by a low density envelope with a Lorentzian profile. The lack of sharp density gradients in such a model inhibits the growth of Kelvin-Helmholtz instabilities, consistent with the cloud's appearance. Our model results also suggest that cloud clumps have densities $\sim$ 10 times the ambient ISM density and account for $\sim$ 30% of the total cloud volume. Moreover, the observed spacing of internal cloud shocks and model simulations indicate that the distance between clumps is $\sim$ 4 clump radii.
The Pierre Auger Observatory data have been analyzed to search for excesses of events near the direction of the galactic center in several energy ranges around EeV energies. In this region the statistics accumulated by the Observatory are already larger than that of any previous experiment. Using both the data sets from the surface detector and our hybrid data sets (events detected simultaneously by the surface detector and the fluorescence detector) we do not find any significant excess. At our present level of undestanding of the performance and properties of our detector, our results do not support the excesses reported by AGASA and SUGAR experiments. We set an upper bound on the flux of cosmic rays arriving within a few degrees from the galactic center in the energy range from 0.8-3.2 EeV. We also have searched for correlations of cosmic ray arrival directions with the galactic plane and with the super-galactic plane at energies in the range 1-5 EeV and above 5 EeV and have found no significant excess.
Since Bekenstein's (2004) creation of his Tensor-Vector-Scalar theory (TeVeS), the Modified Newtonian dynamics (MOND) paradigm has been redeemed from the embarrassment of lacking a relativistic version. One primary success of TeVeS is that it provides an enhancement of gravitational lensing, which could not be achieved by other MONDian theories. Following Bekenstein's work, we investigate the phenomena of gravitational lensing including deflection angles, lens equations and time delay. We find that the deflection angle would maintain its value while the distance of closest approach vary in the MOND regime, this coincides with the conclusion of Mortlock and Turner's (2001) intuitional approach. Moreover, the scalar field, which is introduced to enhance the deflection angle in TeVeS, contributes a negative effect on the potential time delay. Unfortunately this phenomenon is unmeasurable in lensing systems where we can only observe the time delay between two images for a given source. However, this measurable time delay offers another constraint on the mass ratio of the DM and MOND scenarios, which in general differs from that given by the deflection angle. In other words, for a lensing system, if two masses, m_gN and m_gM, are mutually alternatives for the deflection angles in their own paradigm, regarding the time delay they are in general in an exclusive relation.
We present a calculation of the damping of an ultra-energetic (UHE) cosmic neutrino travelling through the thermal gas of relic neutrinos, using the formalism of finite-temperature field theory. From the self-energy diagram due to Z exchange, we obtain the annihilation cross section for an UHE neutrino interacting with an antineutrino from the background. This method allows us to derive the full expressions for the UHE neutrino transmission probability, taking into account the momentum of relic neutrinos. We compare our results with the approximations in use in the literature. We discuss the effect of thermal motion on the shape of the absorption dips for different UHE neutrino fluxes as well as in the context of relic neutrino clustering. We find that for ratios of the neutrino mass to the relic background temperature $10^2$ or smaller, the thermal broadening of the absorption lines could significantly affect the determination of the neutrino mass and of the characteristics of the population of UHE sources.
The Central Laser Facility is located near the middle of the Pierre Auger Observatory in Argentina. It features a UV laser and optics that direct a beam of calibrated pulsed light into the sky. Light scattered from this beam produces tracks in the Auger optical detectors which normally record nitrogen fluorescence tracks from cosmic ray air showers. The Central Laser Facility provides a "test beam" to investigate properties of the atmosphere and the fluorescence detectors. The laser can send light via optical fiber simultaneously to the nearest surface detector tank for hybrid timing analyses. We describe the facility and show some examples of its many uses.
We report on radiative hydrodynamic simulations of moderate and strong solar flares. The flares were simulated by calculating the atmospheric response to a beam of non-thermal electrons injected at the apex of a one-dimensional closed coronal loop, and include heating from thermal soft X-ray, extreme ultraviolet and ultraviolet (XEUV) emission. The equations of radiative transfer and statistical equilibrium were treated in non-LTE and solved for numerous transitions of hydrogen, helium, and Ca II allowing the calculation of detailed line profiles and continuum emission. This work improves upon previous simulations by incorporating more realistic non-thermal electron beam models and includes a more rigorous model of thermal XEUV heating. We find XEUV backwarming contributes less than 10% of the heating, even in strong flares. The simulations show elevated coronal and transition region densities resulting in dramatic increases in line and continuum emission in both the UV and optical regions. The optical continuum reaches a peak increase of several percent which is consistent with enhancements observed in solar white light flares. For a moderate flare (~M-class), the dynamics are characterized by a long gentle phase of near balance between flare heating and radiative cooling, followed by an explosive phase with beam heating dominating over cooling and characterized by strong hydrodynamic waves. For a strong flare (~X-class), the gentle phase is much shorter, and we speculate that for even stronger flares the gentle phase may be essentially non-existent. During the explosive phase, synthetic profiles for lines formed in the upper chromosphere and transition region show blue shifts corresponding to a plasma velocity of ~120 km/s, and lines formed in the lower chromosphere show red shifts of ~40 km/s.
We have fit synthetic visibilities from 3-D (CO5BOLD + PHOENIX) and 1-D (PHOENIX, ATLAS 12) model stellar atmospheres of Procyon (F5 IV) to high-precision interferometric data from the VLTI Interferometer (K-band) and from the Mark III interferometer (500 nm and 800 nm). These data sets provide a test of theoretical wavelength dependent limb-darkening predictions. The work of Allende Prieto et al. has shown that the temperature structure from a spatially and temporally averaged 3-D hydrodynamical model produces significantly less limb darkening at 500 nm relative to the temperature structure of a 1-D MARCS model atmosphere with a standard mixing-length approximation for convection. Our direct fits to the interferometric data confirm this prediction. A 1-D ATLAS 12 model with ``approximate overshooting'' provides the required temperature gradient. We show, however, that 1-D models cannot reproduce the ultraviolet spectrophotometry below 160 nm with effective temperatures in the range constrained by the measured bolometric flux and angular diameter. We find that a good match to the full spectral energy distribution can be obtained with a composite model consisting of a weighted average of twelve 1-D model atmospheres based on the surface intensity distribution of a 3-D granulation simulation. We emphasize that 1-D models with overshooting may realistically represent the mean temperature structure of F-type stars like Procyon, but the same models will predict redder colors than observed because they lack the multicomponent temperature distribution expected for the surfaces of these stars.
The brightest giant flare from the soft $\gamma$-ray repeater (SGR) 1806-20 was detected on 2004 December 27. The isotropic-equivalent energy release of this burst is at least one order of magnitude more energetic than those of the two other SGR giant flares. Starting from about one week after the burst, a very bright ($\sim 80$ mJy), fading radio afterglow was detected. Follow-up observations revealed the multi-frequency light curves of the afterglow and the temporal evolution of the source size. Here we show that these observations can be understood in a two-component explosion model. In this model, one component is a relativistic collimated outflow responsible for the initial giant flare and the early afterglow, and another component is a subrelativistic wider outflow responsible for the late afterglow. We also discuss triggering mechanisms of these two components within the framework of the magnetar model.
Baryonic acoustic oscillations imprinted in the galaxy power spectrum provide a promising tool for probing the cosmological distance scale and dark energy. We present results from a suite of cosmological N-body simulations aimed at investigating possible systematic errors in the recovery of cosmological distances. We show the robustness of baryonic peaks against nonlinearity, redshift distortions, and mild biases within the linear and quasilinear region at various redshifts. While mildly biased tracers follow the matter power spectrum well, redshift distortions do partially obscure baryonic features in redshift space compared to real space. We calculate the statistical constraints on cosmological distortions from N-body results and compare these to the analytic results from a Fisher matrix formalism. We conclude that the angular diameter distance will be constrained as well as our previous Fisher matrix calculations while the Hubble parameter will be less constrained because of nonlinear redshift distortions.
In this paper, we describe our vectorized and parallelized adaptive mesh refinement (AMR) N-body code with shared time steps, and report its performance on a Fujitsu VPP5000 vector-parallel supercomputer. Our AMR N-body code puts hierarchical meshes recursively where higher resolution is required and the time step of all particles are the same. The parts which are the most difficult to vectorize are loops that access the mesh data and particle data. We vectorized such parts by changing the loop structure, so that the innermost loop steps through the cells instead of the particles in each cell, in other words, by changing the loop order from the depth-first order to the breadth-first order. Mass assignment is also vectorizable using this loop order exchange and splitting the loop into $2^{N_{dim}}$ loops, if the cloud-in-cell scheme is adopted. Here, $N_{dim}$ is the number of dimension. These vectorization schemes which eliminate the unvectorized loops are applicable to parallelization of loops for shared-memory multiprocessors. We also parallelized our code for distributed memory machines. The important part of parallelization is data decomposition. We sorted the hierarchical mesh data by the Morton order, or the recursive N-shaped order, level by level and split and allocated the mesh data to the processors. Particles are allocated to the processor to which the finest refined cells including the particles are also assigned. Our timing analysis using the $\Lambda$-dominated cold dark matter simulations shows that our parallel code speeds up almost ideally up to 32 processors, the largest number of processors in our test.
We present a new analysis of the abundances observed in extremely metal poor stars based on both a new generation of theoretical presupernova models and explosions of zero metallicity massive stars and a new abundance analysis of an homogeneous sample of stars having [Fe/H]<-2.5 (Cayrel et al. 2004).
We present heavy element abundance measurements for 16 DAO white dwarfs, determined from Far-Ultraviolet Spectroscopic Explorer {FUSE} spectra. Evidence of absorption by heavy elements was found in the spectra of all the objects. Measurements were made using models that adopted the temperatures, gravities and helium abundances determined from both optical and FUSE data by Good et al. (2004). Despite the different evolutionary paths that the types of DAO white dwarfs are thought to evolve through, their abundances were not found to vary significantly, apart from for the silicon abundances. Abundances measured when the FUSE derived values of temperature, gravity and helium abundance were adopted were, in general, a factor 1-10 higher than those determined when the optical measure of those parameters was used. Satisfactory fits to the absorption lines were achieved in approximately equal number. The models that used the FUSE determined parameters seemed better at reproducing the strength of the nitrogen and iron lines, while for oxygen, the optical parameters were better.
We present time-resolved optical spectroscopy and photometry of four relatively bright (V~14.0-15.5) long-period cataclysmic variables (CVs) discovered in the Hamburg Quasar Survey: HS0139+0559, HS0229+8016, HS0506+7725 and HS0642+5049. Their respective orbital periods, 243.69+-0.49min, 232.550+-0.049min, 212.7+-0.2min and 225.90+-0.23min are determined from radial velocity and photometric variability studies. HS0506+7725 is characterised by strong Balmer and He emission lines, short-period (~10-20min) flickering and weak X-ray emission in the ROSAT All Sky Survey. The detection of a deep low state (~18.5) identifies HS0506+7725 as a member of the VY Scl stars. HS0139+0559, HS0229+8016 and HS0642+5049 display thick-disc like spectra and no or only weak flickering activity. HS0139+0559 and HS0229+8016 exhibit clean quasi-sinusoidal radial velocity varations of their emission lines but no or very little orbital photometric variability. In contrast, we detect no radial velocity variation in HS0642+5049 but a noticeable orbital brightness variation. We identify all three systems either as UX UMa-type novalike variables or as Z Cam-type dwarf novae. Our identification of these four new systems underlines that the currently known sample of CVs is rather incomplete even for bright objects. The four new systems add to the clustering of orbital periods in the 3-4h range found in the sample of HQS selected CVs, and we discuss the large incidence of magnetic CVs and VY Scl/SW Sex stars found in this period range among the known population of CVs.
We present the evolution of rotation in models of massive single stars
covering a wide range of masses and metallicities. These models reproduce very
well observations during the early stages of the evolution (in particular WR
populations and ratio between type II and type Ib,c at different metallicities,
see Meynet & Maeder 2005).
Our models predict the production of fast rotating black holes. Models with
large initial masses or high metallicity end their life with less angular
momentum in their central remnant with respect to the break-up limit for the
remnant. Many WR star models satisfy the three main criteria (black hole
formation, loss of hydrogen-rich envelope and enough angular momentum to form
an accretion disk around the black hole) for gamma-ray bursts (GRB) production
via the collapsar model (Woosley 1993). Considering all types of WR stars as
GRB progenitors, there would be too many GRBs compared to observations. If we
consider only WO stars (type Ic supernovae as is the case for
SN2003dh/GRB030329, see Matheson et al. 2003) as GRBs progenitors, the GRBs
production rates are in much better agreement with observations. WO stars are
produced only at low metallicities in the present grid of models. This
prediction can be tested by future observations.
The ROSAT all-sky survey discovered several `shrapnels', showing boomerang structures outside the Vela supernova remnant. We observed shrapnel D with the XMM-Newton satellite. There is an X-ray bright ridge structure in our FOV running from north to south. Applying the VNEI model to X-ray spectra of various regions, we find that the plasma in the eastern part from the X-ray ridge is significantly different from that in the western part. The X-ray spectra in the western part can be represented by a single-temperature component. The abundances of heavy elements are almost uniform, whereas they are heavily overabundant, except for Fe; the relative abundances to the solar values are O about 5, Ne about 10, Mg about 10, Fe about 1. This indicates that shrapnel D originated from the ejecta of the supernova. We find that the plasma in the eastern part from the ridge consists of two components with different temperatures; the hot component comes from the ejecta, while the cold component comes from the interstellar matter. These two components are considered to be in contact with each other, forming a contact discontinuity. Around the northern part of the contact discontinuity, we find wave-like structures of which the typical scale are comparable with that of the Rayleigh--Taylor instability.
We investigate which conditions of dust attenuation and stellar populations allow models of dusty, continuously star-forming, bulge-less disk galaxies at 0.8<z<3.2 to meet the different colour selection criteria of high-z ``red'' galaxies (e.g. Rc-K>5.3, Ic-K>4, J-K>2.3). As a main novelty, we use stellar population models that include the thermally pulsating Asymptotic Giant Branch (TP-AGB) phase of stellar evolution. The star formation rate of the models declines exponentially as a function of time, the e-folding time being longer than 3 Gyr. In addition, we use calculations of radiative transfer of the stellar and scattered radiation through different dusty interstellar media in order to explore the wide parameter space of dust attenuation. We find that synthetic disks can exhibit red optical/near-infrared colours because of reddening by dust, but only if they have been forming stars for at least about 1 Gyr. Extremely few models barely exhibit Rc-K>5.3, if the inclination i=90 deg and if the opacity 2*tauV>6. Hence, Rc-K-selected galaxies at 1<z<2 most probably are either systems with an old, passively evolving bulge or starbursts. Synthetic disks at 1<z<2 exhibit 4<Ic-K<4.8, if they are seen edge on (i.e. at i about 90 deg) and if 2*tauV>0.5. This explains the large fraction of observed, edge-on disk-like galaxies with Ks<19.5 and F814W-Ks>4. Finally, models with 2<z<3.2 exhibit 2.3<J-K<3, with no bias towards i about 90 deg and for a large range in opacity (e.g. 2*tauV>1 for i about 70 deg). In conclusion, red disk-like galaxies at 0.8<z<3.2 may not necessarily be dustier than nearby disk galaxies (with 0.5<2*tauV<2) and/or much older than about 1 Gyr. This result is due both to a realistic description of dust attenuation and to the emission contribution by TP-AGB stars... (Abridged)
We review the current cosmological status of neutrinos, with particular emphasis on their effects on Big Bang Nucleosynthesis, Large Scale Structure of the universe and Cosmic Microwave Background Radiation measurements.
Absolute calibration of the Pierre Auger Observatory fluorescence detectors uses a light source at the telescope aperture. The technique accounts for the ombined effects of all detector components in a single measurement. The calibrated 2.5 m diameter light source fills the aperture, providing uniform illumination to each pixel. The known flux from the light source and the response of the acquisition system give the required calibration for each pixel. In the lab, light source uniformity is studied using CCD images and the intensity is measured relative to NIST-calibrated photodiodes. Overall uncertainties are presently 12%, and are dominated by systematics.
We present the results of a systematic study using ASCA of the ionization state for six ``mixed-morphology'' supernova emnants (MMSNRs): IC 443, W49B, W28, W44, 3C391, and Kes 27. MMSNRs show centrally filled thermal X-ray emission, which contrasts to shell-like radio morphology, a set of haracteristics at odds with the standard model of SNR evolution (e.g., the Sedov model). We have therefore studied the evolution of the MMSNRs from the ionization conditions inferred from the X-ray spectra, independent of X-ray morphology. We find highly ionized plasmas approaching ionization equilibrium in all the mmsnrs. The degree of ionization is systematically higher than the plasma usually seen in shell-like SNRs. Radial temperature gradients are also observed in five remnants, with cooler plasma toward the limb. In IC 443 and W49B, we find a plasma structure consistent with shell-like SNRs, suggesting that at least some MMSNRs have experienced similar evolution to shell-like SNRs. In addition to the results above, we have discovered an ``overionized'' ionization state in W49B, in addition to that previously found in IC 443. Thermal conduction can cause the hot interior plasma to become overionized by reducing the temperature and density gradients, leading to an interior density increase and temperature decrease. Therefore, we suggest that the ``center-filled'' X-ray morphology develops as the result of thermal conduction, and should arise in all SNRs. This is consistent with the results that MMSNRs are near collisional ionization equilibrium since the conduction timescale is roughly similar to the ionization timescale. Hence, we conclude that MMSNRs are those that have evolved over$\sim10^4$ yr. We call this phase as the ``conduction phase.''
We extend the rotation measure work of Burn (1966) to the cases of limited sampling of lambda squared space and non-constant emission spectra. We introduce the rotation measure transfer function (RMTF), which is an excellent predictor of n-pi ambiguity problems with the lambda squared coverage. Rotation measure synthesis can be implemented very efficiently on modern computers. Because the analysis is easily applied to wide fields, one can conduct very fast RM surveys of weak spatially extended sources. Difficult situations, for example multiple sources along the line of sight, are easily detected and transparently handled. Under certain conditions, it is even possible to recover the emission as a function of Faraday depth within a single cloud of ionized gas. Rotation measure synthesis has already been successful in discovering widespread, weak, polarized emission associated with the Perseus cluster (De Bruyn and Brentjens, 2005). In simple, high signal to noise situations it is as good as traditional linear fits to polarization angle versus lambda squared plots. However, when the situation is more complex or very weak polarized emission at high rotation measures is expected, it is the only viable option.
Phantom energy can be visualized as a scalar field with a (non-canonical) negative kinetic energy term. In this letter, we use the dynamical system formalism to study the attractor behavior of a cosmological model containing a phantom scalar field $\phi$ endowed with an exponential potential of the form $V(\phi)=V_0 \exp(-\lambda \kappa \phi)$, and a perfect fluid with constant equation of state $\gamma$; the latter can be of the phantom type too. As in the canonical case, three characteristic solutions can be identified. The scaling solution exists but is either unstable or of no physical interest. Thus, there are only two stable critical points which are of physical interest, corresponding to the perfect fluid and scalar field dominated solutions, respectively. The most interesting case arises for $0> \gamma > -\lambda^2/3$, which allows the coexistence of the three solutions. The main features of each solution are discussed in turn.
We report on full-polarization radio observations of the Perseus cluster (Abell 426) using the Westerbork Synthesis Radio Telescope (WSRT) at wavelengths from 81-95 cm. We have employed a novel technique, Rotation Measure synthesis (Brentjens and de Bruyn, 2005) to unravel the polarization properties of the emission across the full field of view and detect polarized emission over a wide range of RM from about 0 to 90 rad m^-2. The low RM emission is associated with our Galaxy, while the high RM emission is associated with the Perseus cluster. The latter reaches typical surface brightness levels of 0.5-1 mJy per beam and must be rather highly polarized. Most of the peripheral polarized emission appears too bright, by about 1-2 orders of magnitude, to be explainable as Thomson scattered emission of the central radio source off the thermal electrons in the cluster. The bulk of the emission associated with the Perseus cluster is probably related to buoyant bubbles of relativistic plasma, probably relics from still active or now dormant AGN within the cluster. A lenticular shaped structure measuring 0.5-1 Mpc is strikingly similar to the structures predicted by Ensslin et al. (1998). At the western edge of the cluster, we detect very long, linear structures that may be related to shocks caused by infall of gas into the Perseus cluster.
We have performed a pilot study of mass loss predictions for late-type Wolf-Rayet (WR) stars as a function of metal abundance, over a range between 10^{-5} < (Z/Zsun) < 10. We find that the winds of nitrogen-rich Wolf-Rayet stars are dominated by iron lines, with a dependence of mass loss on Z similar to that of massive OB stars. For more evolved, carbon-rich, WR stars the wind strength is also found to be dependent on the Fe abundance, so that they depend on the chemical environment of the host galaxy, but with a mass loss metallicity dependence that is less steep than for OB stars. Our finding that WR mass loss is Z-dependent is a new one, with important consequences for black hole formation and X-ray population studies in external galaxies. A further finding of our study is that the Z dependence of C-rich WR stars becomes weaker at metallicities below Z/Zsun < 1/10, and mass loss no longer declines once the metal abundance drops below (Z/Zsun) = 10^{-3}. This is the result of an increased importance of radiative driving by intermediate mass elements, such as carbon. In combination with rapid rotation and/or proximity to the Eddington limit -- likely to be relevant for massive Population III stars -- this effect may indicate a role for mass loss in the appearance and evolution of these objects, as well as a potential role for stellar winds in enriching the intergalactic medium of the early Universe.
Absorbers seen in the spectrum of background quasars are a unique tool to select HI-rich galaxies at all redshifts. In turns, these allow to determine the cosmological evolution of the HI gas, Omega_HI+HeII, a possible indicator of gas consumption as star formation proceeds. The Damped Lyman-alpha systems (DLAs with N(HI) > 10^20.3 cm^-2), in particular, are believed to contain a large fraction of the HI gas but there are also indications that lower column density systems, named ``sub-Damped Lyman-alpha'' systems play a role at high-redshift. Here we present the discovery of high-redshift sub-DLAs based on 17 z>4 quasar spectra observed with the Ultraviolet-Visual Echelle Spectrograph (UVES) on VLT. This sample is composed of 21 new sub-DLAs which, together with another 10 systems from previous ESO archive studies, make up a homogeneous sample. The redshift evolution of the number density of several classes of absorbers is derived and shows that all systems seem to be evolving in the redshift range from z=5 to z~3. This is further used to estimate the redshift evolution of the characteristic radius of these classes of absorbers assuming a Holmberg relation and one unique underlying parent population. DLAs are found to have R_* ~ 20 h_100^-1 kpc, while sub-DLAs have R_* ~ 40 h_100^-1 kpc. The redshift evolution of the column density distribution, f(N,z), down to N(HI) = 10^19 cm^-2 is also presented. A departure from a power law due to a flattening of f(N,z) in the sub-DLA regime is present in the data. f(N,z) is further used to determine the HI gas mass contained in sub-DLAs at z>2. The complete sample shows that sub-DLAs are important at all redshifts from z=5 to z=2. Finally, the possibility that sub-DLAs are less affected by the effects of dust obscuration than classical DLAs are discussed.
We present a detailed analysis of the intrinsic X-ray absorption in the Seyfert 1 galaxy NGC 4151 using Chandra/HETGS data obtained 2002 May, as part of a program which included simultaneous UV spectra using HST/STIS and FUSE. NGC 4151 was in a relatively low flux state during the observations reported here, although roughly 2.5 times as bright in the 2 --10 keV band as during a Chandra observation in 2000. The soft X-ray band was dominated by emission lines, which show no discernible variation in flux between the two observations. The 2002 data show the presence of a very highly ionized absorber, in the form of H-like and He-like Mg, Si, and S lines, as well as lower ionization gas via the presence of inner-shell absorption lines from lower-ionization species of these elements. The former is too highly ionized to be radiatively accelerated in a sub-Eddington source such as NGC 4151. We find that the lower ionization gas had a column density a factor of ~ 3 higher during the 2000 observation. If due to bulk motion, we estimate that this component must have a velocity of more than 1250 km/sec transverse to our line-of-sight. We suggest that these results are consistent with a magneto-hydrodynamic flow.
We report the discovery of a nearly edge-on disk about the A0 star HD 32297 seen in light scattered by the disk grains revealed in NICMOS PSF-subtracted coronagraphic images. The disk extends to a distance of at least 400 AU (3.3") along its major axis with a 1.1 micron flux density of 4.81 +/-0.57 mJy beyond a radius of 0.3" from the coronagraphically occulted star. The fraction of 1.1 micron starlight scattered by the disk, 0.0033 +/- 0.0004, is comparable to its fractional excess emission at 25 + 60 micron of ~ 0.0027 as measured from IRAS data. The disk appears to be inclined 10.5 degrees +/- 2.5 degrees from an edge-on viewing geometry, with its major axis oriented 236.5 degrees +/- 1 degree eastward of north. The disk exhibits unequal brightness in opposing sides and a break in the surface brightness profile along NE-side disk major axis. Such asymmetries might implicate the existence of one or more (unseen) planetary mass companions.