We measure the clustering of DEEP2 galaxies at z=1 as a function of luminosity on scales 0.05 Mpc/h to 20 Mpc/h. Drawing from a parent catalog of 25,000 galaxies at 0.7<z<1.3 in the full DEEP2 survey, we create volume-limited samples having upper luminosity limits between M_B=-19 and M_B=-20.5, roughly 0.2-1 L^* at z=1. We find that brighter galaxies are more strongly clustered than fainter galaxies and that the slope of the correlation function does not depend on luminosity for L<L^*. The brightest galaxies, with L>L^*, show a steeper slope and a strong rise in clustering amplitude on small scales. The relative bias of galaxies as a function of L/L^* is steeper than the relation found locally for SDSS galaxies (Zehavi et al. 2005) over the luminosity range that we sample. The absolute bias of galaxies at z=1 is scale-dependent on scales r_p<1 Mpc/h, and rises most significantly on small scales for the brightest samples. For a concordance cosmology, the large-scale bias varies from 1.26 +/-0.04 to 1.54 +/-0.05 as a function of luminosity and implies that DEEP2 galaxies reside in dark matter halos with a minimum mass of ~1-3 10^12 h^-1 M_sun.
We employ high-resolution dissipationless simulations of the concordance LambdaCDM cosmology to model the observed luminosity dependence and evolution of galaxy clustering through most of the age of the universe, from z~5 to z~0. We use a simple, non-parametric model which monotonically relates galaxy luminosities to the maximum circular velocity of dark matter halos (V_max) by preserving the observed galaxy luminosity function in order to match the halos in simulations with observed galaxies. The novel feature of the model is the use of the maximum circular velocity at the time of accretion, V_max,acc, for subhalos, the halos located within virial regions of larger halos. We argue that for subhalos in dissipationless simulations, V_max,acc reflects the luminosity and stellar mass of the associated galaxies better than the circular velocity at the epoch of observation, V_max,now. The simulations and our model L-V_max relation predict the shape, amplitude, and luminosity dependence of the two-point correlation function in excellent agreement with the observed galaxy clustering in the SDSS data at z~0 and in the DEEP2 samples at z~1 over the entire probed range of projected separations, 0.1<r_p/(h^-1 Mpc)<10.0. In particular, the small-scale upturn of the correlation function from the power-law form in the SDSS and DEEP2 luminosity-selected samples is reproduced very well. At z~3-5, our predictions also match the observed shape and amplitude of the angular two-point correlation function of Lyman-break galaxies (LBGs) on both large and small scales, including the small-scale upturn. (ABRIDGED)
We describe a high resolution Smoothed Particle Hydrodynamics (SPH) simulation of a spheroidal galaxy starting from \LCDM initial conditions which is intentionally simple, and includes photoionization and cooling of the intergalactic medium but not feedback from AGN or supernovae. The galaxy undergoes an initial burst of star formation at $z \approx 5$, accompanied by the formation of a bubble of heated gas and does not experience a major merger after $z=3$. Heating from shocks and -PdV work dominates over cooling so that for most of the gas the temperature is an increasing function of time. By $z \approx 1$, 80% of the final stellar spheroid is in place and the spectral energy distribution resembles those of observed extremely red objects (EROs). By the present day, the simulated galaxy is an old ($\approx 10 {\rm Gyrs}$), kinematically hot stellar system with a stellar mass of $\approx 1.2 \times 10^{11} M_{\odot}$, surrounded by a hot gaseous halo containing 40% of the baryonic matter. Stars dominate the mass of the galaxy up to $\approx 4$ effective radii ($\approx 10$ kpc). Projected photometric and kinematic properties are in good agreement with observed field elliptical galaxies. In particular, unlike some recent simulations, the final stellar system has a concentration that is quite typical of real elliptical galaxies. Our simulation shows that a realistic intermediate mass giant elliptical galaxy with a plausible formation history can be formed from \LCDM initial conditions without requiring recent major mergers or feedback from supernovae or AGN.
The brightnesses of supernovae are commonly understood to indicate that
cosmological expansion is accelerating due to dark energy. However the entire
discussion presumes a perfectly transparent universe because no effects of
reddening associated with the interstellar extinction law are seen. We note
that with two kinds of dark matter (baryonic and non-baryonic) strongly
dominating the known mass of the universe, it is seriously premature to assume
that these dark matter components have not reduced the transmission of the
universe for cosmological sources.
We show that the long-known $Lyman-\alpha$ clouds, if nucleated by the
population of baryonic dark matter primordial planetoids indicated by quasar
microlensing, would act as spherical lenses and achromatically fade
cosmologically distant sources. We attempt to estimate the amount of this
cosmological fading, but ultimately the calculation is limited by lack of a
satisfactory model for the tenuous outer parts of a primordial planetoid. We
also consider the effects of such cosmological fading on the light of quasars.
Protoplanetary nebulae typically present non-spherical envelopes. The origin of such geometry is still controversial. There are indications that it may be carried over from an earlier phase of stellar evolution, such as the AGB phase. But how early in the star's evolution does the non-spherical envelope appear? Li-rich giants show dusty circumstellar envelopes that can help answer that question. We study a sample of fourteen Li-rich giants using optical polarimetry in order to detect non-spherical envelopes around them. We used the IAGPOL imaging polarimeter to obtain optical linear polarization measurements in V band. Foreground polarization was estimated using the field stars in each CCD frame. After foreground polarization was removed, seven objects presented low intrinsic polarization (0.19 - 0.34)% and two (V859 Aql and GCSS 557) showed high intrinsic polarization values (0.87 - 1.16)%. This intrinsic polarization suggests that Li-rich giants present a non-spherical distribution of circumstellar dust. The intrinsic polarization level is probably related to the viewing angle of the envelope, with higher levels indicating objects viewed closer to edge-on. The correlation of the observed polarization with optical color excess gives additional support to the circumstellar origin of the intrinsic polarization in Li-rich giants. The intrinsic polarization correlates even better with the IRAS 25 microns far infrared emission. Analysis of spectral energy distributions for the sample show dust temperatures for the envelopes tend to be between 190 and 260 K. We suggest that dust scattering is indeed responsible for the optical intrinsic polarization in Li-rich giants. Our findings indicate that non-spherical envelopes may appear as early as the red giant phase of stellar evolution.
We consider geometries and possible physical models for weak low ionization absorbers based on the relative incidence of low and high ionization absorption systems. We found a total of 16 metal-line systems, with low and/or high ionization absorption detected in our survey of weak low ionization absorption systems from the archive of HST/STIS data. The weak low ionization absorbers trace an abundant population of metal-enriched regions (close to solar metallicity). Generally, models show that these systems have a ~10pc region of higher density gas and a ~1kpc region of lower density phase of higher ionization absorption. We compare absorption systems detected in low and/or high ionization gas and find the following: 1) All but 1 of the 10 weak low ionization systems have a related high ionization phase. In 3 cases the high ionization gas has only a single component, kinematically centered on the low ionization absorption, and in the other 6 cases there are additional high ionization components offset in velocity. The system, toward quasar 3C 273, do not have a high ionization cloud; 2) There are just 6 systems with only a high ionization phase as compared to the 9 systems with both low and high ionization phases; 3) The high ionization absorption in weak low ionization systems is, on average, stronger than in systems with only high ionization absorption; 4) The kinematic structure of the high ionization in weak low ionization systems is similar to that in high ionization only systems. We find that filamentary and sheetlike geometries are favored, due to the relatively small observed cross-section of high ionization only systems. Although low ionization absorbers are not closely associated with luminous galaxies, they arise in their immediate environments within the cosmic web.
We report the discoveries of two, two-image gravitationally lensed quasars selected from the Sloan Digital Sky Survey: SDSS J0806+2006 at z_s=1.540 and SDSS J1353+1138 at z_s=1.629 with image separations of 1.40" and 1.41" respectively. Spectroscopic and optical/near-infrared imaging follow-up observations show that the quasar images have identical redshifts and possess extended objects between the images that are likely to be lens galaxies at z_l~0.6 in SDSS J0806+2006 and z_l~0.3 in SDSS J1353+1138. The field of SDSS J0806+2006 contains several nearby galaxies that may significantly perturb the system, and SDSS J1353+1138 has an extra component near its Einstein ring that is probably a foreground star. Simple mass models with reasonable parameters reproduce the quasar positions and fluxes of both systems.
We report on new X-ray observations of the large-scale jets recently discovered in X-rays from the black hole candidate 4U 1755-33. Our observations in 2004 show that the jets found in 2001 are still present in X-rays. However, sensitive radio observations in 2004 failed to detect the jets. We suggest that synchrotron radiation is a viable emission mechanism for the jets and that thermal bremsstrahlung and inverse-Compton emission are unlikely on energetic grounds. In the synchrotron interpretation, the production of X-rays requires acceleration of electrons up to ~ 60 TeV, the jet power is ~ 4x10^35 erg s^-1, and the radio non-detection requires a spectral index alpha > -0.65 (S_\nu \propto \nu^\alpha) which is similar to the indexes found in lobes surrounding some other compact objects. We find an upper limit on the flux of 4U 1755-33 in quiescence of 5x10^-16 erg cm^-2 s^-1 in the 0.3-8 keV band.
Non-thermal TeV $\gamma$-ray emission within multi-pc scale has been observed from the center region of our galaxy. We argue that these $\gamma$-rays are the result of a transient activity of the massive black hole Sgr A$^*$ which resides at the Galactic center. About thousands of years ago, the black hole may have experienced an active phase by capturing a red giant star and forming an accretion disk, temporarily behaving like an active galactic nuclear. A powerful jet, which contains plenty of high speed protons, was launched during the process. These runaway protons interact with the dense ambient medium, producing TeV $\gamma$-ray emission through the $\pi^\circ$-decay process. We show that the total energy deposited in this way is large enough to account for observations. The diffusion length of protons is also consistent with the observed size of the TeV source.
We discuss a new class of Micro Pattern Gas Detectors, the Gas Pixel Detector (GPD), in which a complete integration between the gas amplification structure and the read-out electronics has been reached. An Application-Specific Integrated Circuit (ASIC) built in deep sub-micron technology has been developed to realize a monolithic device that is, at the same time, the pixelized charge collecting electrode and the amplifying, shaping and charge measuring front-end electronics. The CMOS chip has the top metal layer patterned in a matrix of 80 micron pitch hexagonal pixels, each of them directly connected to the underneath electronics chain which has been realized in the remaining five layers of the 0.35 micron VLSI technology. Results from tests of a first prototype of such detector with 2k pixels and a full scale version with 22k pixels are presented. The application of this device for Astronomical X-Ray Polarimetry is discussed. The experimental detector response to polarized and unpolarized X-ray radiation is shown. Results from a full MonteCarlo simulation for two astronomical sources, the Crab Nebula and the Hercules X1, are also reported.
We study the physical properties derived from interstellar cloud complexes having a fractal structure. We first generate fractal clouds with a given fractal dimension and associate each clump with a maximum in the resulting density field. Then, we discuss the effect that different criteria for clump selection has on the derived global properties. We calculate the masses, sizes and average densities of the clumps as a function of the fractal dimension (D_f) and the fraction of the total mass in the form of clumps (epsilon). In general, clump mass does not fulfill a simple power law with size of the type M_cl ~ (R_cl)**(gamma), instead the power changes, from gamma ~ 3 at small sizes to gamma<3 at larger sizes. The number of clumps per logarithmic mass interval can be fitted to a power law N_cl ~ (M_cl)**(-alpha_M) in the range of relatively large masses, and the corresponding size distribution is N_cl ~ (R_cl)**(-alpha_R) at large sizes. When all the mass is forming clumps (epsilon=1) we obtain that as D_f increases from 2 to 3 alpha_M increases from ~0.3 to ~0.6 and alpha_R increases from ~1.0 to ~2.1. Comparison with observations suggests that D_f ~ 2.6 is roughly consistent with the average properties of the ISM. On the other hand, as the fraction of mass in clumps decreases (epsilon<1) alpha_M increases and alpha_R decreases. When only ~10% of the complex mass is in the form of dense clumps we obtain alpha_M ~ 1.2 for D_f=2.6 (not very different from the Salpeter value 1.35), suggesting this a likely link between the stellar initial mass function and the internal structure of molecular cloud complexes.
The vertical profiles of disc galaxies are built by the material trapped around stable periodic orbits, which form their "skeletons". According to this, the knowledge of the stability of the main families of periodic orbits in appropriate 3D models, can predict possible morphologies for edge-on disc galaxies. In a pilot survey we compare the orbital structures which lead to the appearance of "peanuts" and "X"-like features with the edge-on profiles of three disc galaxies (IC 2531, NGC 4013 and UGC 2048). The subtraction from the images of a model representing the axisymmetric component of the galaxies reveals the contribution of the non-axisymmetric terms. We find a direct correspondence between the orbital profiles of 3D bars in models and the observed main morphological features of the residuals. We also apply a simple unsharp masking technique in order to study the sharpest features of the images. Our basic conclusion is that the morphology of the boxy "bulges" of these galaxies can be explained by considering disc material trapped around stable 3D periodic orbits. In most models these building-blocks periodic orbits are bifurcated from the planar central family of a non-axisymmetric component, usually a bar, at low order vertical resonances. In such a case the boxy "bulges" are parts of bars seen edge-on. For the three galaxies we study the families associated with the "peanut" or "X"-shape morphology are most probably bifurcations at the vertical 2/1 or 4/1 resonance.
We present the current performances of the AMBER / VLTI instrument in terms of differential observables (differential phase and differential visibility) and show that we are already able to reach a sufficient precision for very low mass companions spectroscopy and mass characterization. We perform some extrapolations with the knowledge of the current limitations of the instrument facility. We show that with the current setup of the AMBER instrument, we can already reach $3\sigma = 10^{-3}$ radians and have the potential to some low mass companions characterization (Brown dwarves or hypothetical very hot Extra Solar Giant Planets). With some upgrades of the VLTI infrastructure, improvements of the instrument calibration and improvements of the observing strategy, we will be able to reach $3\sigma = 10^{-4}$ radians and will have the potential to perform Extra Solar Giant Planets spectroscopy and mass characterization.
We have proposed a new High Resolution Shock Capturing (HRSC) scheme for Special Relativistic Hydrodynamics (SRHD) based on the semidiscrete central Godunov-type schemes and a modified Weighted Essentially Non-oscillatory (WENO) data reconstruction algorithm. This is the first application of the semidiscrete central schemes with high order WENO data reconstruction to the SRHD equations. This method does not use a Riemann solver for flux computations and a number of one and two dimensional benchmark tests show that the algorithm is robust and comparable in accuracy to other SRHD codes.
Jeans showed analytically that, in an infinite uniform-density isothermal gas, plane-wave perturbations collapse to dense sheets if their wavelength, $\lambda$, satisfies $\lambda > \lambda_{_{\rm JEANS}} = (\pi a^2 / G \rho_{_0})^{1/2}$ (where $a$ is the isothermal sound speed and $\rho_{_0}$ is the unperturbed density); in contrast, perturbations with smaller $\lambda$ oscillate about the uniform density state. Here we show that Smoothed Particle Hydrodynamics reproduces these results well, even when the diameters of the SPH particles are twice the wavelength of the perturbation. Our simulations are performed in 3-D with initially settled (i.e. non-crystalline) distributions of particles. Therefore there exists the seed noise for artificial fragmentation, but it does not occur. We conclude that, although there may be -- as with any numerical scheme -- `skeletons in the SPH cupboard', a propensity to fragment artificially is evidently not one of them.
The development of fast numerical methods for multilevel radiative transfer (RT) applications often leads to important breakthroughs in astrophysics, because they allow the investigation of problems that could not be properly tackled using the methods previously available. Probably, the most familiar example is the so-called Multilevel Accelerated $\Lambda$-Iteration (MALI) technique of Rybicki & Hummer for the case of a local approximate operator, which is based on Jacobi iteration. However, there are superior operator-splitting methods, based on Gauss-Seidel (GS) and Successive Overrelaxation (SOR) iteration, which provide a dramatic increase in the speed with which non-LTE multilevel transfer problems can be solved in one, two and three-dimensional geometries. Such RT methods, which were introduced by Trujillo Bueno & Fabiani Bendicho ten years ago, are the main subject of the first part of this paper. We show in some detail how they can be applied for solving multilevel RT problems in spherical geometry, for both atomic and molecular line transitions. The second part of the article addresses the issue of the calculation of the molecular number densities when the approximation of instantaneous chemical equilibrium turns out to be inadequate, which happens to be the case whenever the dynamical time scales of the astrophysical plasma under consideration are much shorter than the time needed by the molecules to form.
We present the results of two XMM-Newton observations of Jupiter carried out in 2003 for 100 and 250 ks (or 3 and 7 planet rotations) respectively. X-ray images from the EPIC CCD cameras show prominent emission from the auroral regions in the 0.2 - 2.0 keV band: the spectra are well modelled by a combination of emission lines, including most prominently those of highly ionised oxygen (OVII and OVIII). In addition, and for the first time, XMM-Newton reveals the presence in both aurorae of a higher energy component (3 - 7 keV) which is well described by an electron bremsstrahlung spectrum. This component is found to be variable in flux and spectral shape during the Nov. 2003 observation, which corresponded to an extended period of intense solar activity. Emission from the equatorial regions of Jupiter's disk is also observed, with a spectrum consistent with that of solar X-rays scattered in the planet's upper atmosphere. Jupiter's X-rays are spectrally resolved with the RGS which clearly separates the prominent OVII contribution of the aurorae from the OVIII, FeXVII and MgXI lines, originating in the low-latitude disk regions of the planet.
We present multi-frequency, multi-epoch radio imaging of the complex radio
source B2151+174 in the core of the cluster, Abell 2390 (z~0.23). From new and
literature data we conclude that the FRII-powerful radio source is the
combination of a compact, core-dominated `medium-symmetric object' (MSO) with a
more extended, steeper spectrum mini-halo. B2151+174 is unusual in a number of
important aspects: i) it is one of the most compact and flat spectrum sources
in a cluster core known; ii) it shows a complex, compact twin-jet structure in
a north-south orientation; iii) the orientation of the jets is 45deg misaligned
with apparent structure (ionization cones and dust disk) of the host galaxy on
larger scales. Since the twin-jet of the MSO has its northern half with an
apparent `twist', it might be that precession of the central supermassive black
hole explains this misalignement.
B2151+174 may be an example of the early stage (10^3-10^4 yrs duration) of a
`bubble' being blown into the ICM where the plasma has yet to expand.
We used the 8 June 2004 transit of Venus (ToV) as a surrogate to test observing methods, strategies and techniques that are being contemplated for future space missions to detect and characterize extrasolar terrestrial planets (ETPs) as they transit their host stars, notably NASA's Kepler mission planned for 2008. As an analog to "Kepler-like" photometric transit observations, we obtained (spatially unresolved) radiometric observations with the ACRIM 3 instrument on ACRIMSAT to follow the effect of the ToV on the total solar irradiance (TSI). Contemporaneous high-resolution broadband imagery with NASA's TRACE spacecraft provided, directly, measures of the stellar (solar) astrophysical noise that can intrinsically limit such transit observations. During the ~ 5.5 h transit, the planet's angular diameter was approximately 1/32 the solar diameter, thus covering ~ 0.1 of the stellar surface. With our ACRIM 3 data, we measure temporal changes in TSI with a 1 sigma per sample (unbinned) uncertainty of approximately 100 mW m^-2 (0.007%). A diminution in TSI of ~ 1.4 W m^-2 (~ 0.1%, closely corresponding to the geometrically occulted area of the photosphere) was measured at mid-transit compared with a mean pre/post transit TSI of ~ 1365.9 W m^-2. These observations serve as a surrogate to future photometric observations of ETPs such as Kepler will deliver. Detailed analysis of the ToV, a rare event within our own solar system, with time-resolved radiometry augmented with high-resolution imagery provides a useful analogue for investigating the detectability and characterization of ETPs from observations that are anticipated in the near future.
We have mapped the SiO J=5-4 line at 217GHz from the HH211 molecular outflow with the Submillimeter Array (SMA). The high resolution map (1.6''x0.9'') shows that the SiO J=5-4 emission comes from the central narrow jet along the outflow axis with a width of ~0.8'' (~250 AU) FWHM. The SiO jet consists of a chain of knots separated by 3-4'' (~1000 AU) and most of the SiO knots have counterparts in shocked H_2 emission seen in a new, deep VLT near-infrared image of the outflow. A new, innermost pair of knots are discovered at just +/-2'' from the central star. The line ratio between the SiO J=5-4 data and upper limits from the SiO J=1-0 data of Chandler & Richer (2001) suggests that these knots have a temperature in excess of 300-500 K and a density of (0.5-1) x10^7 cm^{-3}. The radial velocity measured for these knots is ~30 km/s, comparable to the maximum velocity seen in the entire jet. The high temperature, high density, and velocity structure observed in this pair of SiO knots suggest that they are closely related to the primary jet launched close to the protostar.
We have computed, based on 17 infrared radial velocities, the first set of orbital elements for the M giant in the symbiotic binary V2116 Ophiuchi. The giant's companion is a neutron star, the bright X-ray source GX 1+4. We find an orbital period of 1161 days by far the longest of any known X-ray binary. The orbit has a modest eccentricity of 0.10 with an orbital circularization time of less than 10^6 years. The large mass function of the orbit significantly restricts the mass of the M giant. Adopting a neutron-star mass of 1.35M(Sun), the maximum mass of the M giant is 1.22M(Sun), making it the less massive star. Derived abundances indicate a slightly subsolar metallicity. Carbon and nitrogen are in the expected ratio resulting from the red-giant first dredge-up phase. The lack of O-17 suggests that the M-giant has a mass less than 1.3M(Sun), consistent with our maximum mass. The red giant radius is 103R(Sun), much smaller than the estimated Roche lobe radius. Thus, the mass loss of the red giant is via a stellar wind. Although the M giant companion to the neutron star has a mass similar to the late-type star in low-mass X-ray binaries, its near-solar abundances and apparent runaway velocity are not fully consistent with the properties of this class of stars.
We report spectroscopic observations of HD 27638B, the secondary in a visual binary in which the physically associated primary (separation approximately 19 arcsec) is a B9V star. The secondary shows strong Li 6708 absorption suggesting youth, and has attracted attention in the past as a candidate post-T Tauri star although this has subsequently been ruled out. It was previously known to be a double-lined spectroscopic binary (F8+G6) with a period of 17.6 days, and to show velocity residuals indicating a more distant massive third companion with a period of at least 8 years. Based on our radial velocity measurements covering more than two cycles of the outer orbit, along with other measurements, we derive an accurate triple orbital solution giving an outer period of 9.447 +/- 0.017 yr. The third object is more massive than either of the other two components of HD 27638B, but is not apparent in the spectra. We derive absolute visual magnitudes and effective temperatures for the three visible stars in HD 27638. Isochrone fitting based on those properties gives an age of 200 +/- 50 Myr for the system. We infer also an inclination angle of about 53.3 degrees for the inner orbit of HD 27638B. We detect a near-infrared excess in HD 27638B which we attribute to the third star being a close binary composed of late-type stars. This explains its large mass and lack of a visible signature. Modeling of this excess allows us to infer not only the masses of the components of the unseen companion, but also the inclination angle of the outer orbit (approximately 73 deg). The HD 27638 system is thus at least quintuple.
We consider the structure of self-gravitating marginally stable accretion disks in galactic centers in which a small fraction of the disk mass has been converted into proto-stars. We find that proto-stars accrete gaseous disk matter at prodigious rates. Mainly due to the stellar accretion luminosity, the disk heats up and geometrically thickens, shutting off further disk fragmentation. The existing proto-stars however continue to gain mass by gas accretion. As a results, the initial mass function for disk-born stars at distances R ~ 0.03-3 parsec from the super-massive black hole should be top-heavy. The effect is most pronounced at around R ~ 0.1 parsec. We suggest that this result explains observations of rings of young massive stars in our Galaxy and in M31, and predict that more of such rings will be discovered.
This text was originally written to accompany a series of lectures that I gave at the `35th Saas-Fee advanced course' in Switzerland and at the Institute for Astronomy of the University of Hawaii. It reviews my current understanding of the dynamics of comets and of the origin and primordial sculpting of their reservoirs. It starts discussing the structure of the Kuiper belt and the current dynamics of Kuiper belt objects, including scattered disk objects. Then it discusses the dynamical evolution of Jupiter family comets from the trans-Neptunian region, and of long period comets from the Oort cloud. The formation of the Oort cloud is then reviewed, as well as the primordial sculpting of the Kuiper belt. Finally, these issues are revisited in the light of a new model of giant planets evolution that has been developed to explain the origin of the late heavy bombardment of the terrestrial planets.
In an investigative 16 hour L band observation using the MERLIN radio interferometric array, we have resolved both the pulsar PSR B1951+32 and structure within the flat spectral radio continuum region, believed to be the synchrotron nebula associated with the interaction of the pulsar and its `host' supernova remnant CTB 80. The extended structure we see, significant at $\sim$ 4.5 $\sigma$, is of dimensions 2.5" $\times$ 0.75", and suggests a sharp bow shaped arc of shocked emission, which is correlated with similar structure observed in lower resolution radio maps and X-ray images. Using this MERLIN data as a new astrometric reference for other multiwavelength data we can place the pulsar at one edge of the HST reported optical synchrotron knot, ruling out previous suggested optical counterparts, and allowing an elementary analysis of the optical synchrotron emission which appears to trail the pulsar. The latter is possibly a consequence of pulsar wind replenishment, and we suggest that the knot is a result of magnetohydrodynamic (MHD) instabilities. These being so, it suggests a dynamical nature to the optical knot, which will require high resolution optical observations to confirm.
We report the initial results of a VLT/NACO high spatial resolution imaging survey for multiple systems among 58 M-type members of the nearby Upper Scorpius OB association. Nine pairs with separations below 100 have been resolved. Their small angular separations and the similarity in the brightness of the components (DMagK <1 for all of them), indicate there is a reasonable likelihood several of them are true binaries rather than chance projections. Follow-up imaging observations with WHT/LIRIS of the two widest binaries confirm that their near-infrared colours are consistent with physical very low mass binaries. For one of these two binaries, WHT/LIRIS spectra of each component were obtained. We find that the two components have similar M6-M7 spectral types and signatures of low-gravity, as expected for a young brown dwarf binary in this association. Our preliminary results indicate a possible population of very low-mass binaries with semimajor axis in the range 100 AU 150 AU, which has not been seen in the Pleiades open cluster. If these candidates are confirmed (one is confirmed by this work), these results would indicate that the binary properties of very low-mass stars and brown dwarfs may depend on the environment where they form.
We investigate the dissociation equilibrium of $\rm H_2$ in very cool, helium-rich white dwarf atmospheres. We present the solution of the non-ideal chemical equilibrium for the dissociation of molecular hydrogen in a medium of dense helium. We find that at the photosphere of cool white dwarfs of $T_{\rm eff}\rm=4000 K$, the non-ideality results in an increase of the mole fraction of molecular hydrogen by up to a factor of $\sim 10$, compared to the equilibrium value for the ideal gas. This increases the $\rm H_{2}-He$ CIA opacity by an order of magnitude and will affect the determination of the abundance of hydrogen in very cool, helium-rich white dwarfs.
We present the results of X-ray observations of four bright transients sources detected in the July 2004 XMM-Newton observations of the central bulge of M31. Two X-ray sources, XMMU J004315.5+412440 and XMMU J004144.7+411110, were discovered for the first time. Two other sources, CXOM31 J004309.9+412332 and CXOM31 J004241.8+411635, were previously detected by Chandra. The properties of the sources suggest their identification with accreting binary systems in M31. The X-ray spectra and variability of two sources, XMMU J004144.7+411110 and CXOM31 J004241.8+411635, are similar to that of the Galactic black hole transients, making them a good black hole candidates. The X-ray source XMMU J004315.5+412440 demonstrates a dramatic decline of the X-ray flux on a time scale of three days, and a remarkable flaring behavior on a short time scales. The X-ray data on XMMU J004315.5+412440 and CXOM31 J004309.9+412332 suggest that they can be either black hole or neutron star systems. Combining the results of 2000-2004 XMM observations of M31, we estimate a total rate of the bright transient outbursts in the central region of M31 to be 6-12 per year, in agreement with previous studies.
Curves in a family derived from powers of the polar coordinate formula for ellipses are found to provide good fits to bound orbits in a range of power-law potentials. This range includes the well-known $1/r$ (Keplerian) and logarithmic potentials. These approximate orbits, called p-ellipses, retain some of the basic geometric properties of ellipses. They satisfy and generalize Newton's apsidal precession formula, which is one of the reasons for their surprising accuracy. Because of their simplicity the p-ellipses make very useful tools for studying trends among power-law potentials, and especially the occurence of closed orbits. The occurence of closed or nearly closed orbits in different potentials highlights the possibility of period resonances between precessing, eccentric orbits and circular orbits, or between the precession period of multi-lobed closed orbits and satellite periods. These orbits and their resonances promise to help illuminate a number of problems in galaxy and accretion disk dynamics.
The first generation of epoch of reionization experiments have the potential to characterize the processes and history of reionization for redshifts 6<z<12 by observing the power spectrum of the redshifted 21 cm radiation. These power spectrum observations could also contribute to our knowledge of the underlying cosmology, particularly at redshifts before the effects of reionization strongly dominate the cosmic radio background. In this paper we explore the cosmological constraints of 21 cm power spectrum observations with the Mileura Widefield Array (MWA) and its potential successor, the MWA5000, assuming a flat LCDM cosmology and that reionization occurs below a redshift of z=8. In addition to a standard set of cosmological parameters, we include two nuisance parameters related to the systematics and foregrounds faced by radio observations. While the MWA observations cannot constrain the underlying cosmology, MWA5000 could provide useful constraints, in particular on the slope of the inflationary power spectrum, n_s, and the running of the spectral index, alpha_s.
A number of radio interferometers are currently being planned or constructed to observe 21 cm emission from reionization. Not only will such measurements provide a detailed view of the Epoch of Reionization, but, since the 21 cm emission also traces the distribution of matter in the Universe, this signal can be used to constrain cosmological parameters at redshifts 6 ~< z ~< 20. The sensitivity of an interferometer to the cosmological information in the signal may depend on how precisely the angular dependence of the 21 cm 3-D power spectrum can be measured. Here, we quantify all the effects that break the spherical symmetry of the 3-D 21 cm power spectrum and produce physically motivated predictions for this power spectrum, utilizing an analytic model for reionization. We find that upcoming observatories will be sensitive to the 21 cm signal over a wide range of scales, from larger than 100 comoving Mpc to as small as 1 comoving Mpc. We discuss how cosmological parameters can be measured, and show that the first generation of 21 cm observations should moderately improve existing constraints on cosmological parameters for certain low-redshift reionization scenarios. A two year observation with the second generation interferometer MWA5000 in combination with the CMB telescope Planck can improve constraints on Omega_w (to +- 0.019, a 30% improvement over Planck alone), Omega_m h^2 (+- 0.0011, 50%), Omega_b h^2 (+- 0.00013, 30%), Omega_nu (+-0.003, 300%), n_s (+- 0.0034, 30%), and alpha_s (+- 0.004, 100%). Larger interferometers, such as SKA, have the potential to do even better. If the Universe is substantially ionized by z ~ 12 or if spin temperature fluctuations are important, we show that it will be difficult to place competitive constraints on cosmological parameters with the 21 cm signal.
We present our recent results from numerical simulations of a magnetized flow in the vicinity of a black hole in the context of the collapsar model for GRBs. The simulations show that after an initial transient, the flow settles into a complex convolution of several distinct, time-dependent flow components including an accretion torus, its corona and outflow, an inflow and an outflow in the polar funnel. We focus on studying the nature and connection between these components, in particular between the inflows and related outflows. We find that rotational and MHD effects launch, accelerate, and sustain the outflows. We also find that an outflow can be formed even when the collapsing envelope has initially a very weak magnetic field and a very small angular momentum. Our main conclusion is that even for a relatively weak initial magnetic field and a slow rotation, a gravitational collapse of a stellar envelope can lead to formation of a very strong and very fast jet.
We report new limits on cosmic neutrino fluxes from the test flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment, which completed an 18.4 day flight of a prototype long-duration balloon (LDB) payload in early 2004. We search for impulsive events that could be associated with ultra-high energy neutrino interactions in the ice, and derive limits that constrain several models for ultra-high energy neutrino fluxes, including the long-standing Z-burst model.
We report on the identification of 54 embedded clusters around 217 massive protostellar candidates of which 34 clusters are new detections. The embedded clusters are identified as stellar surface density enhancements in the 2 $\mu$m All Sky Survey (2MASS) data. Because the clusters are all associated with massive stars in their earliest evolutionary stage, the clusters should also be in an early stage of evolution. Thus the properties of these clusters should reflect properties associated with their formation rather than their evolution. For each cluster, we estimate the mass, the morphological type, the photometry and extinction. The clusters in our study, by their association with massive protostars and massive outflows, reinstate the notion that massive stars begin to form after the first generation of low mass stars have completed their accretion phase. Further, the observed high gas densities and accretion rates at the centers of these clusters is consistent with the hypothesis that high mass stars form by continuing accretion onto low mass stars.
Temperature maps of the Cosmic Microwave Background (CMB) radiation, as those obtained by the Wilkinson Microwave Anisotropy Probe (WMAP), provide one of the most precise data sets to test fundamental hypotheses of modern cosmology. One of these issues is related to the statistical properties of the CMB temperature fluctuations, which would have been produced by Gaussian random density fluctuations when matter and radiation were in thermal equilibrium in the early Universe. We analysed here the WMAP data and found that the distribution of the CMB temperature fluctuations P^{CMB}(Delta T) can be quite well fitted by the anomalous temperature distribution emerging within nonextensive statistical mechanics. This theory is based on the nonextensive entropy S_q = k (1 - \int dx [P_q(x)]^q) /(q-1), with the Boltzmann-Gibbs expression as the limit case q -> 1. For the frequencies investigated (\nu= 40.7, 60.8, and 93.5 GHz), we found that P^{CMB}(Delta T) is well described by P_q(Delta T) \propto 1/[1 + (q-1) B(\nu) (Delta T)^2]^{1/(q-1)}, with q = 1.055 \pm 0.002, which exclude, at the 99% confidence level, exact Gaussian temperature distributions P^{Gauss}(Delta T) \propto e^{- B(\nu) Delta T^2}, corresponding to the q -> 1 limit, to properly represent the CMB temperature fluctuations measured by WMAP.
We present infrared echelle spectroscopy of three Herbig-Haro (HH) driving sources (SVS 13, B5-IRS 1 and HH 34 IRS) using Subaru-IRCS. The large diameter of the telescope and wide spectral coverage of the spectrograph allowed us to detect several H_2 and [Fe II] lines in the H- and K-bands. These include H_2 lines arising from v=1-3 and J=1-11, and [Fe II] lines with upper level energies of E/k=1.1-2.7x10^4 K. For all objects the outflow is found to have two velocity components: (1) a high-velocity (-70 to -130 km s^-1) component (HVC), seen in [Fe II] or H_2 emission and associated with a collimated jet; and (2) a low-velocity (-10 to -30 km s^-1) component (LVC), which is seen in H_2 emission only and is spatially more compact. Such a kinematic structure resembles optical forbidden emission line outflows associated with classical T Tauri stars, whereas the presence of H_2 emission reflects the low-excitation nature of the outflowing gas close to these protostars. The observed H_2 flux ratios indicate a temperature of 2-3x10^3 K, and a gas density of 10^5 cm^-3 or more, supporting shocks as the heating mechanism. B5-IRS 1 exhibits faint extended emission associated with the H_2-LVC, in which the radial velocity slowly increases with distance from the protostar (by ~20 km s^-1 at ~500 AU). This is explained as warm molecular gas entrained by an unseen wide-angled wind. The [Fe II] flux ratios indicate electron densities to be ~10^4 cm^-3 or greater, similar to forbidden line outflows associated with classical T Tauri stars. Finally the kinematic structure of the [Fe II] emission associated with the base of the B5-IRS 1 and HH 34 IRS outflows is shown to support disk-wind models.
Evidence of small-scale condensations in the magnetic molecular clouds has been accumulating over the past decades through radio and optical/ultraviolet observations. The origin and shape of these small-scale condensations is a disputable issue. Nejad-Asghar & Ghanbari (2004 hereafter NG) have recently studied the effect of the linear thermal instability on the formation of fluctuations in molecular clouds. The authors inferred that under certain conditions (e.g., depending on expansion or contraction of the background) thermal instability and ambipolar diffusion can produce spherical, oblate, or prolate condensations.
We present a detail analysis on the spectral lags of the short gamma-ray bursts (GRBs) and compare them with that of the long GRBs by using the CGRO/BATSE GRB Catalog. Our results are as follows. (1)The spectral lag distribution of the short GRBs is significantly different from that of the long GRBs. Excluding the statistical fluctuation effect, a proportion of ~17% of the short GRBs have a negative spectral lag, i.e., the hard photons being lag behind the soft photons. We do not find any peculiar features from their light curves to distinguish these bursts from those with a positive spectral lag. We argue that a more physical mechanism dominated the hard lag may be hid behind the morphological features of the light curves. This should be a great challenge to the current GRB models. We notice that this proportion is consistent with the proportion of short GRBs correlated with nearby galaxies newly discovered by Tanvir et al., although it is unclear if these short GRBs are indeed associated with the sources originated at low redshift. (2)While the spectral lags of the long GRBs are strongly correlated with the pulse durations, they are not for the short GRBs. However, the ratios of the spectral lag to the pulse duration for the short and long GRBs normally distributed at 0.023 and 0.046, respectively, with the sample width, indicating that the curvature effect alone could not explain the difference of the spectral lags between the two types of GRBs. The hydrodynamic timescales of the outflows and the radiative processes at work in GRBs might also play an important role as suggested by Daigne and Mochkovitch.
We explore the consequences of new observational and theoretical evidence that long gamma-ray bursts prefer low metallicity environments. Using recently derived mass-metallicity correlations and the mass function from SDSS studies, and adopting an average cosmic metallicity evolution from \citet{kewley2005} and \citet{savaglio2005} we derive expressions for the the relative number of massive stars formed below a given fraction of solar metallicity, $\epsilon$, as function of redshift. We demonstrate that about 1/10th of all stars form with $\epsilon < 0.1$. Therefore, a picture where the majority of GRBs form with $\epsilon < 0.1$ is not inconsistent with an empirical global SN/GRB ratio of 1/1000. It implies that (1) GRB's peak at a significantly higher redshift than supernovae; (2) massive star evolution at low metallicity may be qualitatively different and; (3) the larger the low-metallicity bias of GRBs the less likely binary evolution channels can be significant GRB producers.
We perform numerical simulations of the Kelvin-Helmholtz instability in the mid-plane of a protoplanetary disc. A two-dimensional corotating slice in the azimuthal--vertical plane of the disc is considered where we include the Coriolis force and the radial advection of the Keplerian rotation flow. Dust grains, treated as individual particles, move under the influence of friction with the gas, while the gas is treated as a compressible fluid. The friction force from the dust grains on the gas leads to a vertical shear in the gas rotation velocity. As the particles settle around the mid-plane due to gravity, the shear increases, and eventually the flow becomes unstable to the Kelvin-Helmholtz instability. The Kelvin-Helmholtz turbulence saturates when the vertical settling of the dust is balanced by the turbulent diffusion away from the mid-plane. The azimuthally averaged state of the self-sustained Kelvin-Helmholtz turbulence is found to have a constant Richardson number in the region around the mid-plane where the dust-to-gas ratio is significant. Nevertheless the dust density has a strong non-axisymmetric component. We identify a powerful clumping mechanism, caused by the dependence of the rotation velocity of the dust grains on the dust-to-gas ratio, as the source of the non-axisymmetry. Our simulations confirm recent findings that the critical Richardson number for Kelvin-Helmholtz instability is around unity or larger, rather than the classical value of 1/4.
We present a short atlas illustrating the unusual Ca {\sc ii} K-line profiles in upper main sequence stars with anomalous abundances. Slopes of the profiles for 10 cool, magnetic chemically peculiar (CP) stars change abruptly at the very core, forming a deep "nib." The nibs show the same or nearly the same radial velocity as the other atomic lines. The near wings are generally more shallow than in normal stars. In three magnetic CP stars, the K-lines are too weak to show this shape, though the nibs themselves are arguably present. The Ca {\sc ii} H-lines also show deep nibs, but the profiles are complicated by the nearby, strong H$\epsilon$ absorption. The K-line structure is nearly unchanged with phase in $\beta$ CrB and $\alpha$ Cir. Calculations, including NLTE, show that other possibilities in addition to chemical stratification may yield nib-like cores.
Results from a study of Fe K-alpha emission lines for a sample of six non-magnetic Cataclysmic Variables (CVs) using the high resolution X-ray data from the Chandra High Energy Transmission Grating (HETG) are presented. Two of the sources, SS Cyg and U Gem are observed in both quiescent and outburst states whereas V603 Aql, V426 Oph, WX Hyi and SU UMa are observed only in quiescence. The fluorescent Fe line is prominent in V603 Aql, V426 Oph and SS Cyg during quiescence indicating the presence of a conspicuous reflection component in these sources. The observed equivalent width of the fluorescent Fe line is consistent with reflection from a white dwarf surface that subtends 2pi solid angle at the X-ray source. During the outburst in SS Cyg, the fluorescent line is red-shifted by about 2300 km/s. The Fe XXV triplet at 6.7 keV is found to be dominant in all sources. The value of the G-ratio derived from the Fe XXV triplet indicates that the plasma is in collisional ionization equilibrium during the quiescent state. The Fe XXV line is significantly broadened in U Gem and SS Cyg during the outbursts compared to quiescence, indicating the presence of a high velocity material near the white dwarf during the outburst. The ratio of Fe XXVI/XXV indicates a higher ionization temperature during quiescence than in outburst in U Gem and SS Cyg.
High energy neutrino emission from GRBs is discussed. In this paper, by using the simulation kit GEANT4, we calculate proton cooling efficiency including pion-multiplicity and proton-inelasticity in photomeson production. First, we estimate the maximum energy of accelerated protons in GRBs. Using the obtained results, neutrino flux from one burst and a diffuse neutrino background are evaluated quantitatively. We also take account of cooling processes of pion and muon, which are crucial for resulting neutrino spectra. We confirm the validity of analytic approximate treatments on GRB fiducial parameter sets, but also find that the effects of multiplicity and high-inelasticity can be important on both proton cooling and resulting spectra in some cases. Finally, assuming that the GRB rate traces the star formation rate, we obtain a diffuse neutrino background spectrum from GRBs for specific parameter sets. We introduce the nonthermal baryon-loading factor, rather than assume that GRBs are sources of UHECRs. We find that the obtained neutrino background is comparable with the prediction of Waxman & Bahcall, although our ground in estimation is different from theirs. In this paper, we study on various parameters since there are many parameters in the model. The observations of a neutrino background has a possibility not only to test the internal shock model of GRBs but also to give us information about parameters in the model and whether GRBs are sources of UHECRs or not.
The three-dimensional space distribution of 48921 quasars and 16113 Seyfert galaxies with redshifts z <= 6 is investigated. The global anisotropy caused by the shift of the observer place by \simeq 50 Mpc from a center of their symmetry (supposed center of the Metagalaxy) to the side of the vector with equatorial coordinates \alpha \simeq 13 degrees and \delta \simeq 70 degrees has been found in the placement of these objects. In the opposite direction there exists the extensive region where the progressive decrease of redshifts up to a minimum (near \alpha \simeq 193 degrees, \delta \simeq -70 degrees) is observed. The influence of gravitational potential and possible rotation of the Metagalaxy on the anisotropy of redshifts of the cosmological objects has been considered.
The relation between X-ray luminosity and near-infrared luminosity for
early-type galaxies has been examined. Near-infrared (NIR) luminosities should
provide a superior measure of stellar mass compared to optical luminosities
used in previous studies, especially if there is significant star-formation or
dust present in the galaxies. However, we show that the X-ray-NIR relations are
remarkably consistent with the X-ray-optical relations. This indicates that the
large scatter of the relations is dominated by scatter in the X-ray properties
of early-type galaxies, and is consistent with early-types consisting of old,
quiescent stellar populations.
We have investigated scatter in terms of environment, surface brightness
profile, Mg2, H_beta, H_gamma line strength indices, spectroscopic age, and
nuclear H_alpha emission. We found that galaxies with high Mg2 index, low
H_beta and H_gamma indices or a `core' profile have a large scatter in Lx,
whereas galaxies with low Mg2, high H_beta and H_gamma indices or `power-law'
profiles, generally have Lx<10^41 erg/s. There is no clear trend in the scatter
with environment or nuclear H_alpha emission.
Most, perhaps all, stars go through a phase of vigorous outflow during formation. We examine, through 3D MHD simulation, the effects of protostellar outflows on cluster formation. We find that the initial turbulence in the cluster-forming region is quickly replaced by motions generated by outflows. The protostellar outflow-driven turbulence (``protostellar turbulence'' for short) can keep the region close to a virial equilibrium long after the initial turbulence has decayed away. We argue that there exist two types of turbulence in star-forming clouds: a primordial (or ``interstellar'') turbulence and a protostellar turbulence, with the former transformed into the latter mostly in embedded clusters such as NGC 1333. Since the majority of stars are thought to form in clusters, an implication is that the stellar initial mass function is determined to a large extent by the stars themselves, through outflows which individually limit the mass accretion onto forming stars and collectively shape the environments (density structure and velocity field) in which most cluster members form. We speculate that massive cluster-forming clumps supported by protostellar turbulence gradually evolve towards a highly centrally condensed ``pivotal'' state, culminating in rapid formation of massive stars in the densest part through accretion.
A new method of measuring cosmology with gamma-ray bursts(GRBs) has been proposed by Liang and Zhang recently. In this method, only observable quantities including the rest frame peak energy of the nu F_nu spectrum E'_p, the isotropic energy of GRB (E_{gamma,iso}), and the break time of the optical afterglow light curves in the rest frame (t'_b) are used. By considering this method we constrain the cosmological parameters and the redshift at which the universe expanded from the deceleration to acceleration phase. We add five recently-detected GRBs to the sample and derive E_{gamma, iso}/10^{52} ergs=(0.93\pm0.25)\times (E'_p/100 keV)^{1.91\pm0.32}\times (t'_b/1 day)^{-0.93\pm 0.38} for a flat cosmology with Omega_M=0.28 and H_0=71.0 km s^{-1} Mpc^{-1}. This relation is independent of the medium density around bursts and the efficiency of conversion of the explosion energy to gamma-ray energy that are difficult to measure. We regard the E_{gamma,iso}(E'_p, t'_b) relationship as a standard candle and find 0.05<Omega_M<0.48 and Omega_Lambda<1.15 (at the 1 sigma confidence level). In a flat universe with the cosmological constant we obtain 0.25<Omega_M<0.46 and 0.54<Omega_Lambda<0.78 at the 1 sigma confidence level. The transition redshift is z_T=0.69_{-0.06}^{+0.05}. Combining 20 GRBs and 157 type Ia supernovae, we find Omega_M=0.29\pm 0.03 for a flat universe and the transition redshift is z_T=0.61\pm 0.05, which is slightly larger than the value found by considering SNe Ia alone. In particular, We also discuss several dark-energy models in which the equation of state w(z) includes one-parameter and two-parameter parameterization and investigate the constraints on the model parameters in details.
We present an analysis of fine-structure transitions of FeII and SiII detected in a high-resolution optical spectrum of the afterglow of GRB 051111 (z=1.54948). The fine-structure absorption features arising from FeII* to FeII****, as well as SiII*, are confined to a narrow velocity structure extending over +/-30 km/s, which we interpret as the burst local environment, most likely a star forming region. We investigate two scenarios for the excitation of the fine-structure levels by collisions with electrons and radiative pumping by an infra-red or ultra-violet radiation field produced by intense star formation in the GRB environment, or by the GRB afterglow itself. We find that the conditions required for collisional excitation of FeII fine-structure states cannot be easily reconciled with the relatively weak SiII* absorption. Radiative pumping by either IR or UV emission requires >10^3 massive hot OB stars within a compact star-forming region a few pc in size, and in the case of IR pumping a large dust content. On the other hand, it is possible that the GRB itself provides the source of IR and/or UV radiation, in which case we estimate that the excitation takes place at a distance of ~10-20 pc from the burst. Detailed radiative transfer calculations are required in order to verify that excitation of the low-ionization fine-structure states is possible given the intense UV flux from the burst. Still, it is clear that GRB absorption spectroscopy can provide direct information on the mode and conditions of star formation at high redshift.
We investigate the impact finite simulation box size has on the structural and kinematic properties of Cold Dark Matter haloes forming in cosmological simulations. Our approach involves generating a single realisation of the initial power spectrum of density perturbations and studying how truncation of this power spectrum on scales larger than L_cut affects the structure of dark matter haloes at z=0. In particular, we have examined the cases of L_cut = f_cut L_box with f_cut=1 (i.e. no truncation), 1/2, 1/3 and 1/4. In common with previous studies, we find that the suppression of long wavelength perturbations reduces the strength of clustering, as measured by a suppression of the 2-point correlation function xi(r), and reduces the numbers of the most massive haloes, as reflected in the depletion of the high mass end of the mass function n(M). Interestingly, we find that truncation has little impact on the internal properties of haloes. The masses of high mass haloes decrease in a systematic manner as L_cut is reduced, but the distribution of concentrations is unaffected. On the other hand, the median spin parameter is ~50% lower in runs with f_cut<1. We argue that this is an imprint of the linear growth phase of the halo's angular momentum by tidal torquing, and that the absence of any measurable trend in concentration and the weak trend observed in halo shape reflect the importance of virialisation and complex mass accretion histories for these quantities. These results are of interest for studies that require high mass resolution and statistical samples of simulated haloes, such as simulations of the population of first stars. Our analysis shows that large-scale tidal fields have relatively little effect on the internal properties of Cold Dark Matter haloes and hence may be ignored in such studies.
Various results are obtained for a Friedmann-Robertson-Walker cosmology. We derive an exact equation that determines Hubble's law, clarify issues concerning the speeds of faraway objects and uncover a "tail-light angle effect" for distant luminous sources. The latter leads to a small, previously unnoticed correction to the parallax distance formula.
Recently, the HESS collaboration has reported the detection of gamma-ray emission above a few hundred GeV from eight new sources located close to the Galactic Plane. The source HESS J1813-178 has sparked particular interest, as subsequent radio observations imply an association with SNR G12.82-0.02. Triggered by the detection in VHE gamma-rays, a positionally coincident source has also been found in INTEGRAL and ASCA data. In this Letter we present MAGIC observations of HESS J1813-178, resulting in the detection of a differential gamma-ray flux consistent with a hard-slope power law, described as dN/(dA dt dE) = (3.3+/-0.5)*10^{-12} (E/TeV)^{-2.1+/-0.2} cm^(-2)s^(-1)TeV^(-1). We briefly discuss the observational technique used, the procedure implemented for the data analysis, and put this detection in the perspective of multifrequency observations.
We analyse a light curve of the symbiotic star BF Cyg, covering 114 years of its photometric history. The star had a major outburst around the year 1894. Since then the mean optical brightness of the system is in steady decline, reaching only in the last few years its pre-outburst value. Superposed on this general decline are some 6 less intense outbursts of 1-2 magnitude and duration of 2000-5000 days. We find a cycle of ~6376 days, or possibly twice this period, in the occurrence of these outbursts. We suggest that the origin of the system outbursts is in some magnetic cycle in the outer layers of the giant star of the system, akin to the less intense ~8000 days magnetic cycle of our Sun. We further find, that in addition to its well known binary period of 757.3 days, BF Cyg possesses also another photometric period of 798.8 days. This could be the rotation period of the giant star of the system. If it is, the beat period of these two periodicities, 14580 days, is the rotation period of a tidal wave on the surface of the giant. A 4th period of 4436 days, the beat period of the 14580 and the 6376 cycles is possibly also present in the LC. We predict that BF Cyg will be at the peak of its next outburst around the month of May in the year 2007. The newly discovered 798.8 days period explains the disappearance of the orbital modulation at some epochs in the light curve. The 757.3 oscillations will be damped again around the year 2013.
We have used INTEGRAL & RXTE data to investigate the timing properties of the source in correlation with its spectral states as defined by different positions in the colour-colour diagram. The source shows two distinct branches in the colour-colour diagram that resemble those of the Z sources. The hard branch (similar to the horizontal branch of Z sources) is characterised by a low-amplitude change of the hard colour compared to the change in the soft colour. In the soft branch (analogue to the normal branch) the amplitude of variability of the hard colour is about three times larger than that of the soft colour. As the count rate decreases the source moves up gradually through the soft to the hard branch. The aperiodic variability (excluding the pulse noise) consists of band-limited noise (represented by three broad Lorentzian components) and two QPOs at 0.05 Hz and 0.22 Hz. The strength of the lower frequency QPO increases as the source approaches the hard branch (similar to HBOs in Z sources). The higher frequency QPO reaches maximum significance when the source is in the middle of the branch (like NBOs). We have performed the first measurements of phase lags in the band limited noise below 8 Hz in an accreting X-ray pulsar and found that soft lags dominate at high frequencies. Above the pulse frequency (0.23 Hz), the amplitude of the lag increases as the X-ray flux increases. The Z topology appears to be a signature of the neutron star binaries as it is present in all types of neutron-star binaries (Z, atoll and, as we show here, in accreting pulsars as well). However, the motion along this track, the time scales through the different branches of the diagram and the aperiodic variability associated with portions of the Z track differ for each subclass of neutron-star binary.
Polarized diffuse emission observations at 1.4-GHz in a high Galactic latitude area of the northern Celestial hemisphere are presented. The 3.2 X 3.2 deg^2 field, centred at RA = 10h 58m, Dec = +42deg 18' (B1950), has Galactic coordinates l~172deg, b~+63deg and is located in the region selected as northern target of the BaR-SPOrt experiment. Observations have been performed with the Effelsberg 100-m telescope. We find that the angular power spectra of the E- and B-modes have slopes of beta_E = -1.79 +/- 0.13 and beta_B = -1.74 +/- 0.12, respectively. Because of the very high Galactic latitude and the smooth emission, a weak Faraday rotation action is expected, which allows both a fair extrapolation to Cosmic Microwave Background Polarization (CMBP) frequencies and an estimate of the contamination by Galactic synchrotron emission. We extrapolate the E-mode spectrum up to 32-GHz and confirm the possibility to safely detect the CMBP E-mode signal in the Ka band found in another low emission region (Carretti et al. 2005b). Extrapolated up to 90-GHz, the Galactic synchrotron B-mode looks to compete with the cosmic signal only for models with a tensor-to-scalar perturbation power ratio T/S < 0.001, which is even lower than the T/S value of 0.01 found to be accessible in the only other high Galactic latitude area investigated to date. This suggests that values as low as T/S = 0.01 might be accessed at high Galactic latitudes. Such low emission values can allow a significant red-shift of the best frequency to detect the CMBP B-mode, also reducing the contamination by Galactic dust, and opening interesting perspectives to investigate Inflation models.
Verification of theoretical predictions of an oscillating behavior of the fine-structure constant, alpha, with cosmic time requires high precision measurements at individual redshifts, while in earlier studies the mean Delta alpha/alpha values averaged over wide redshift intervals were usually reported. This requirement can be met via the Single Ion Differential alpha Measurement (SIDAM) procedure. We apply SIDAM to the FeII lines associated with the damped Ly-alpha system observed at z=1.15 in the spectrum of HE0515-4414. The weighted mean calculated on base of carefully selected 34 FeII pairs is <Delta alpha/alpha>=(-0.07+/-0.84)10^{-6}. The precision of this estimate represents the absolute improvement with respect to what has been done in the measurements of Delta alpha/alpha.
We determine the characteristics of the 7mm to 20cm wavelength radio variability in Sgr A* on time scales from days to three decades. The amplitude of the intensity modulation is between 30 and 39% at all wavelengths. Analysis of uniformly sampled data with proper accounting of the sampling errors associated with the lightcurves shows that Sgr A* exhibits no 57- or 106-day quasi-periodic oscillations, contrary to previous claims. The cause of the variability is investigated by examining a number of plausible scintillation models, enabling those variations which could be attributed to interstellar scintillation to be isolated from those that must be intrinsic to the source. Thin-screen scattering models do not account for the variability amplitude on most time scales. However, models in which the scattering region is extended out to a radius of 50-500pc from the Galactic Center account well for the broad characteristics of the variability on >4-day time scales. The ~ 10% variability on <4-day time scales at 0.7-3cm appears to be intrinsic to the source. The degree of scintillation variability expected at millimeter wavelengths depends sensitively on the intrinsic source size; the variations, if due to scintillation, would require an intrinsic source size smaller than that expected.
We present the results of a photometric multisite campaign on the $\delta$ Scuti Pre-Main-Sequence star IP Per. Nine telescopes have been involved in the observations, with a total of about 190 hours of observations over 38 nights. Present data confirms the multiperiodic nature of this star and leads to the identification of at least nine pulsational frequencies. Comparison with the predictions of linear non-adiabatic radial pulsation models allowed us to identify only five of the nine observed frequencies, and to constrain the position of IP Per in the HR diagram. The latter is in good agreement with the empirical determination of the stellar parameters obtained by Miroshnichenko et al. (2001). An initial interpretation of the observed frequencies using the Aarhus non-radial pulsation code suggests that three frequencies could be associated with non-radial ($l$=2) modes. Finally, we present new evolutionary and pulsation models at lower metallicity (Z=0.008) to take into account the possibility that IP Per is metal deficient, as indicated by Miroshnichenko et al. (2001).
The debate on the oxygen abundances of metal-poor stars has its origin in contradictory results obtained using different abundance indicators. To achieve a better understanding of the problem we have acquired high quality spectra with the Ultraviolet and Visual Echelle Spectrograph at VLT, with a signal-to-noise of the order of 100 in the near ultraviolet and 500 in the optical and near infrared wavelength range. Three different oxygen abundance indicators, OH ultraviolet lines around 310.0 nm, the [OI] line at 630.03 nm and the OI lines at 777.1-5 nm were observed in the spectra of 13 metal-poor subgiants with -3.0<=[Fe/H]<=-1.5. Oxygen abundances were obtained from the analysis of these indicators which was carried out assuming local thermodynamic equilibrium and plane-parallel model atmospheres. Abundances derived from OI were corrected for departures from local thermodynamic equilibrium. Stellar parameters were computed using Teff-vs-color calibrations based on the infrared flux method and Balmer line profiles, Hipparcos parallaxes and FeII lines. [O/Fe] values derived from the forbidden line at 630.03 nm are consistent with an oxygen/iron ratio that varies linearly with [Fe/H] as [O/Fe]}=-0.09(+/-0.08)[Fe/H]+0.36(+/-0.15). Values based on the OI triplet are on average 0.19+/-0.22 dex(s.d.) higher than the values based on the forbidden line while the agreement between OH ultraviolet lines and the forbidden line is much better with a mean difference of the order of -0.09+/-0.25 dex(s.d.). In general, our results follow the same trend as previously published results with the exception of the ones based on OH ultraviolet lines. In that case our results lie below the values which gave rise to the oxygen abundance debate for metal-poor stars.
Ultraviolet radiation is known to inhibit photosynthesis, induce DNA destruction and cause damage to a wide variety of proteins and lipids. In particular, UV radiation between 200-300 nm becomes energetically very damaging to most of the terrestrial biological systems. On the other hand, UV radiation is usually considered one of the most important energy source on the primitive Earth for the synthesis of many biochemical compounds and, therefore, essential for several biogenesis processes. In this work, we use these properties of the UV radiation to define the bounderies of an ultraviolet habitable zone. We also analyze the evolution of the UV habitable zone during the main sequence stage of the star. We apply these criteria to study the UV habitable zone for those extrasolar planetary systems that were observed by the International Ultraviolet Explorer (IUE). We analyze the possibility that extrasolar planets and moons could be suitable for life, according to the UV constrains presented in this work and other accepted criteria of habitability (liquid water, orbital stability, etc.).
We analyze the gas accretion flow through a planet-produced gap in a protoplanetary disk. We adopt the alpha disk model and ignore effects of planetary migration. We develop a semi-analytic, one-dimensional model that accounts for the effects of the planet as a mass sink and also carry out two-dimensional hydrodynamical simulations of a planet embedded in a disk. The predictions of the mass flow rate through the gap based on the semi-analytic model generally agree with the hydrodynamical simulations at the 25% level. Through these models, we are able to explore steady state disk structures and over large spatial ranges. The presence of an accreting Jupiter-mass planet significantly lowers the density of the disk within a region of several times the planet's orbital radius. The mass flow rate across the gap (and onto the central star) is typically 10% to 25% of the mass accretion rate outside the orbit of the planet, for planet-to-star mass ratios that range from 5e-5 to 1e-3.
A sample of metal-poor subgiants has been observed with the UVES spectrograph at the Very Large Telescope and abundances of Li and Be have been determined. Typical signal-to-noise per spectral bin values for the co-added spectra are of the order of 500 for the Li I line (670.78 nm) and 100 for the Be II doublet lines (313.04 nm). The spectral analysis of the observations was carried out using the Uppsala suite of codes and MARCS (1D-LTE) model atmospheres with stellar parameters from photometry, parallaxes, isochrones and Fe II lines. Abundance estimates of the light elements were corrected for departures from local thermodynamic equilibrium in the line formation. Effective temperatures and Li abundances seem to be correlated and Be abundances correlate with [O/H]. Standard models predict Li and Be abundances approximately one order of magnitude lower than main-sequence values which is in general agreement with the observations. On average, our observed depletions seem to be 0.1 dex smaller and between 0.2 and 0.4 dex larger (depending on which reference is taken) than those predicted for Li and Be, respectively. This is not surprising since the initial Li abundance, as derived from main-sequence stars on the Spite plateau, may be systematically in error by 0.1 dex or more, and uncertainties in the spectrum normalisation and continuum drawing may affect our Be abundances systematically.
We report abundance excesses of 1.2 and 2.6 dex, respectively, for vanadium and bromine in the hot, peculiar star 3 Cen A. Abundances for these two odd-Z elements have not been previously reported for this star. Taken with previous work, they strengthen the case of the origin of the abundance peculiarities by diffusion.
(abridged) We introduce and discuss the properties of a theoretical (B-K)-(J-K) integrated colour diagram for single-age, single-metallicity stellar populations. This combination of integrated colours is able to largely disentangle the well known age-metallicity degeneracy when the age of the population is greater than ~300 Myr. We discuss in detail the effect on this colour-colour diagram of alpha-enhanced metal abundance ratios, the presence of blue horizontal branch stars unaccounted for in the theoretical calibration, and of statistical colour fluctuations in low mass stellar systems. In the case of populations with multiple stellar generations, the luminosity-weighted mean age obtained from this diagram is shown to be heavily biased towards the youngest stellar components. We apply this method to several datasets for which optical and near-IR photometry are available in the literature. For the two Local Group dwarf galaxies NGC185 and NGC6822, the mean ages derived from the integrated colours are consistent with the star formation histories inferred independently from photometric observations of their resolved stellar populations. A sample of bright field and Virgo cluster elliptical galaxies is found to exhibit a range of luminosity-weighted mean ages from 3 to 14 Gyr, with a mean of 8 Gyr, independent of environment, and mean metallicities at or just above the solar value. Colour gradients are found in all of the galaxies studied, in the sense that central regions are redder. Aperture data for five Virgo early-type dwarf galaxies show that these galaxies appear to be shifted to lower mean metallicities and lower mean ages (range 1 to 6 Gyr) than their higher luminosity counterparts.
We analysed the fine structure of the phase space distribution function of nearby subdwarfs using data extracted from various catalogues. Applying a new search strategy based on Dekker's theory of galactic orbits, we found four overdensely populated regions in phase space. Three of them were correlated with previously known star streams: the Hyades-Pleiades and Hercules streams in the thin disk of the Milky Way and the Arcturus stream in the thick disk. In addition we find evidence for another stream in the thick disk, which resembles closely the Arcturus stream and probably has the same extragalactic origin.
We introduce a method to relate a possible truncation of the star cluster mass function at the high mass end to the shape of the cluster luminosity function (LF). We compare the observed LFs of five galaxies containing young star clusters with synthetic cluster population models with varying initial conditions. The LF of the SMC, the LMC and NGC 5236 are characterized by a power-law behavior NdL~L^-a dL, with a mean exponent of <a> = 2.0 +/- 0.2. This can be explained by a cluster population formeda with a constant cluster formation rate, in which the maximum cluster mass per logarithmic age bin is determined by the size-of-sample effect and therefore increases with log(age/yr). The LFs of NGC 6946 and M51 are better described by a double power-law distribution or a Schechter function. When a cluster population has a mass function that is truncated below the limit given by the size-of-sample effect, the total LF shows a bend at the magnitude of the maximum mass, with the age of the oldest cluster in the population, typically a few Gyr due to disruption. For NGC 6946 and M51 this implies a maximum mass of M_max = 5*10^5 M_sun. Faint-ward of the bend the LF has the same slope as the underlying initial cluster mass function and bright-ward of the bend it is steeper. This behavior can be well explained by our population model. We compare our results with the only other galaxy for which a bend in the LF has been observed, the ``Antennae'' galaxies (NGC 4038/4039). There the bend occurs brighter than in NGC 6946 and M51, corresponding to a maximum cluster mass of M_max = 2*10^6 M_sun (abridged).
We present the luminosity function (LF) of star clusters in M51 based on HST/ACS observations taken as part of the Hubble Heritage project. The clusters are selected based on their size and with the resulting 5990 clusters we present one of the largest cluster samples of a single galaxy. We find that the LF can be approximated with a double power-law distribution with a break around M_V = -8.9. On the bright side the index of the power-law distribution is steeper (a = 2.75) than on the faint-side (a = 1.93), similar to what was found earlier for the ``Antennae'' galaxies. The location of the bend, however, occurs about 1.6 mag fainter in M51. We confront the observed LF with the model for the evolution of integrated properties of cluster populations of Gieles et al., which predicts that a truncated cluster initial mass function would result in a bend in, and a double power-law behaviour of, the integrated LF. The combination of the large field-of view and the high star cluster formation rate of M51 make it possible to detect such a bend in the LF. Hence, we conclude that there exists a fundamental upper limit to the mass of star clusters in M51. Assuming a power-law cluster initial mass function with exponentional cut-off of the form NdM ~ M^-b * exp(-M/M_C)dM, we find that M_C = 10^5 M_sun. A direct comparison with the LF of the ``Antennae'' suggests that there M_C = 4*10^5 M_sun.
We report a detailed attempt to model the ortho-to-para abundance ratio of c-C3H2 so as to reproduce observed values in the cores of the well-known source TMC-1. According to observations, the ortho-to-para ratios vary, within large uncertainties, from a low of near unity to a high of approximately three depending on the core. We used the osu.2003 network of gas-phase chemical reactions augmented by reactions that specifically consider the formation, depletion, and interconversion of the ortho and para forms of the c-C3H2 and its precursor ion c-C3H3+. We investigated the sensitivity of the calculated ortho-to-para ratio for c-C3H2 to a large number of factors. For the less evolved cores C, CP, and D, we had no difficulty reproducing the observed ortho-to-para ratios of 1-2.
In the framework of the Chandrasekhar-mass deflagration model for Type Ia supernovae (SNe Ia), a persisting free parameter is the initial morphology of the flame front, which is linked to the ignition process in the progenitor white dwarf. Previous analytical models indicate that the thermal runaway is driven by temperature perturbations (''bubbles'') that develop in the white dwarf's convective core. In order to probe the conditions at ignition (diameters, temperatures and evolutionary timescales), we have performed hydrodynamical 2D simulations of buoyant bubbles in white dwarf interiors. Our results show that fragmentation occurring during the bubble rise affects the outcome of the bubble evolution. Possible implications for the ignition process of SNe Ia are discussed.
Most star complexes are in fact complexes of stars, clusters and gas clouds; term "star complexes" was introduced as general one disregarding the preferential content of a complex. Generally the high rate of star formation in a complex is accompanied by the high number of bound clusters, including massive ones, what was explained by the high gas pressure in such regions. However, there are also complexes, where clusters seems to be more numerous in relation to stars than in a common complex. The high rate of clusters - but not isolated stars - formation seems to be typical for many isolated bursts of star formation, but deficit of stars might be still explained by the observational selection. The latter cannot, however, explain the complexes or the dwarf galaxies, where the high formation rate of only stars is observed. The possibility of the very fast dissolution of parental clusters just in such regions should itself be explained. Some difference in the physical conditions (turbulence parameters ?) within the initial gas supercloud might be a reason for the high or low stars/clusters number ratio in a complex.
Most star complexes are in fact complexes of stars, clusters and gas clouds; term "star complexes" was introduced as general one disregarding the preferential content of a complex. Generally the high rate of star formation in a complex is accompanied by the high number of bound clusters, including massive ones, what was explained by the high gas pressure in such regions. However, there are also complexes, where clusters seems to be more numerous in relation to stars than in a common complex. The high rate of clusters - but not isolated stars - formation seems to be typical for many isolated bursts of star formation, but deficit of stars might be still explained by the observational selection. The latter cannot, however, explain the complexes or the dwarf galaxies, where the high formation rate of only stars is observed. The possibility of the very fast dissolution of parental clusters just in such regions should itself be explained. Some difference in the physical conditions (turbulence parameters ?) within the initial gas supercloud might be a reason for the high or low stars/clusters number ratio in a complex.
We calculate the spectra of ultra-high-energy cosmic rays (UHECRs) in an explicit top-down model based on the decays of metastable neutral `crypton' states in a flipped SU(5) string model. For each of the eight specific 10th-order superpotential operators that might dominate crypton decays, we calculate the spectra of both protons and photons, using a code incorporating supersymmetric evolution of the injected spectra. For all the decay operators, the total UHECR spectra are compatible with the available data. Also, the fractions of photons are compatible with all the published upper limits, but may be detectable in future experiments.
We present the stellar and gas kinematics of a sample of 18 nearby late-type spiral galaxies (Hubble types ranging from Sb to Sd), observed with the integral-field spectrograph SAURON at the 4.2-m William Herschel Telescope. SAURON covers the spectral range 4800-5380 A, allowing us to measure the Hbeta, Fe, Mgb absorption features and the emission in the Hbeta line and the [OIII], and [NI] doublets over a 33x41 arcsec field of view. The maps cover the nuclear region of these late-type galaxies and in all cases include the entire bulge. In many cases the stellar kinematics suggests the presence of a cold inner region, as visible from a central drop in the stellar velocity dispersion. The ionised gas is almost ubiquitous and behaves in a complicated fashion: the gas velocity fields often display more features than the stellar ones, including wiggles in the zero-velocity lines, irregular distributions, ring-like structures. The line ratio [OIII]/Hbeta often takes on low values over most of the field, probably indicating a wide-spread star formation.
We have obtained 2D spectral data for a sample of 58 nearby S0 galaxies with the Multi-Pupil Spectrograph of the 6m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. The Lick indices H-beta, Mgb, and <Fe> are calculated separately for the nuclei and for the bulges taken as the rings between R=4" and 7"; and the luminosity-weighted ages, metallicities, and Mg/Fe ratios of the stellar populations are estimated by confronting the data to SSP models. Four types of galaxy environments are considered: clusters, centers of groups, other places in groups, and field. The nuclei are found to be on average slightly younger than the bulges in any types of environments, and the bulges of S0s in the sparse environments are younger than those in the dense environments. The effect can be partly attributed to the well-known age correlation with the stellar velocity dispersion in early-type galaxies (in our sample the galaxies in sparse environements are on average less massive than those in dense environments). However for the most massive S0s, with the stellar velocity dispersion of 170-220 km/s, the age dependence on the environment is still significant at the confidence level of 1.5 sigma.
We observed the HH 211 jet in the submillimeter continuum and the CO(3-2) and SiO(8-7) transitions with the Submillimeter Array. The continuum source detected at the center of the outflow shows an elongated morphology, perpendicular to the direction of the outflow axis. The high-velocity emission of both molecules shows a knotty and highly collimated structure. The SiO(8-7) emission at the base of the outflow, close to the driving source, spans a wide range of velocities, from -20 up to 40 km s^{-1}. This suggests that a wide-angle wind may be the driving mechanism of the HH 211 outflow. For distances greater than 5" (1500 AU) from the driving source, emission from both transitions follows a Hubble-law behavior, with SiO(8-7) reaching higher velocities than CO(3-2), and being located upstream of the CO(3-2) knots. This indicates that the SiO(8-7) emission is likely tracing entrained gas very close to the primary jet, while the CO(3-2) is tracing less dense entrained gas. From the SiO(5-4) data of Hirano et al. we find that the SiO(8-7)/SiO(5-4) brightness temperature ratio along the jet decreases for knots far from the driving source. This is consistent with the density decreasing along the jet, from (3-10)x10^6 cm^{-3} at 500 AU to (0.8-4)x10^6 cm^{-3} at 5000 AU from the driving source.
We study the colors, structural properties, and star formation histories of a sample of ~1600 dwarfs over look-back times of ~3 Gyr (z=0.002-0.25). The sample consists of 401 distant dwarfs drawn from the Galaxy Evolution from Morphologies and SEDs (GEMS) survey, which provides high resolution Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) images and accurate redshifts, and of 1291 dwarfs at 10-90 Mpc compiled from the Sloan Digitized Sky Survey (SDSS). We find that the GEMS dwarfs are bluer than the SDSS dwarfs, which is consistent with star formation histories involving starbursts and periods of continuous star formation. The full range of colors cannot be reproduced by single starbursts or constant star formation alone. We derive the star formation rates of the GEMS dwarfs and estimate the mechanical luminosities needed for a complete removal of their gas. We find that a large fraction of luminous dwarfs are likely to retain their gas, whereas fainter dwarfs are susceptible to a significant gas loss, if they would experience a starburst.
We present an analysis of Chandra observations of the galaxy clusters A2670 and A2107. Their cD galaxies have large peculiar velocities (>200km/s) and thus the clusters appear to be undergoing mergers. In A2670, we find a comet-like structure around one of the brightest galaxies. At the leading edge of the structure, there is a cold front. The mass of the X-ray gas in the comet-like structure suggests that the galaxy was in a small cluster or group, and its intracluster medium (ICM) is being stripped by ram-pressure. The regions of cool interstellar medium (ISM) of the cD galaxies in A2670 and A2107 are very compact. This is similar to the brightest galaxies in the Coma cluster, which is also a merging cluster. In each galaxy, the short cooling time of the ISM requires a heating source; the compact nature of the ISM makes it unlikely that the heating source is a central active galactic nucleus (AGN).
The collapse of non-collisional dark matter and the formation of pancake structures in the Universe are investigated approximately. Collapse is described by a system of ordinary differential equations, in the model of a uniformly rotating, three-axis, uniform density ellipsoid. Violent relaxation, mass and angular momentum losses are taken into account phenomenologically. The formation of the equilibrium configuration, secular instability and the transition from a spheroid to a three-axis ellipsoid are investigated numerically and analytically in this dynamical model.
We present a wide grid of models for the structure and transmission properties of warm absorbers in active galactic nuclei (AGN). Contrary to commonly used constant density models, our absorbing cloud is assumed to be under constant total (gas plus radiation) pressure. This assumption implies the coexistence of material at different temperatures and ionization states, which is a natural consequence of pressure and thermal equilibrium. Our photoionization code allows us to compute the profiles of the density, the temperature, the gas pressure, the radiation pressure and the ionization state across the cloud, and to calculate the radiative transfer of continuum and lines including Compton scattering. Therefore, equivalent widths of both saturated and unsaturated lines are properly modeled. For each pair of the incident spectrum slope and the ionization parameter at the cloud surface there is a natural upper limit to the total column densities of the cloud due to thermal instabilities. These maximum values are comparable to the observational constraints on the column density of warm absorbers which may give support to constant total pressure models. In all models we note considerable absorption around 6.4 keV which modifies the intrinsic relativistically broadened iron line profile originating in an accretion disk illuminated atmosphere. Our models can be applied to fitting the spectroscopic data from the {\it XMM-Newton} and {\it Chandra} satellites.
We investigate the stability of differentially rotating proto-neutron stars (PNSs) with a toroidal magnetic field. Stability criteria for nonaxisymmetric MHD instabilities are derived using a local linear analysis. PNSs are expected to have much stronger radial shear in the rotation velocity compared to normal stars. We find that nonaxisymmetric magnetorotational instability (NMRI) with a large azimuthal wavenumber $m$ is dominant over the kink mode ($m=1$) in differentially rotating PNSs. The growth rate of the NMRI is of the order of the angular velocity $\Omega$ which is faster than that of the kink-type instability by several orders of magnitude. The stability criteria are analogous to those of the axisymmetric magnetorotational instability with a poloidal field, although the effects of leptonic gradients are considered in our analysis. The NMRI can grow even in convectively stable layers if the wavevectors of unstable modes are parallel to the restoring force by the Brunt-V\"ais\"al\"a oscillation. The nonlinear evolution of NMRI could amplify the magnetic fields and drive MHD turbulence in PNSs, which may lead to enhancement of the neutrino luminosity.
We present interstellar elemental abundance measurements derived from Space Telescope Imaging Spectrograph echelle observations of 47 sight lines extending up to 6.5 kpc through the Galactic disk. These paths probe a variety of interstellar environments, covering ranges of nearly four orders of magnitude in molecular hydrogen fraction f(H_2) and more than two in mean hydrogen sight line density <n_H>. Coupling the current data with Goddard High Resolution Spectrograph data from 17 additional sight lines and the corresponding Far Ultraviolet Spectroscopic Explorer and Copernicus observations of H_2 absorption features, we explore magnesium, phosphorus, manganese, nickel, copper, and germanium gas-phase abundance variations as a function of <n_H>: density-dependent depletion is noted for each element, consistent with a smooth transition between two abundance plateaus identified with warm and cold neutral interstellar medium depletion levels. The observed scatter with respect to an analytic description of these transitions implies that total elemental abundances are homogeneous on length scales of hundreds of parsecs, to the limits of abundance measurement uncertainty. The probable upper limit we determine for intrinsic variability at any <n_H> is 0.04 dex, aside from an apparent 0.10 dex deficit in copper (and oxygen) abundances within 800 pc of the Sun. Magnesium dust abundances are shown to scale with the amount of silicon in dust and, in combination with a similar relationship between iron and silicon, these data appear to favor the young F and G star values of Sofia & Meyer (2001; ApJ 554, L221) as an elemental abundance standard for the Galaxy.
We examine the variability of the high-ionizaton C IV line in a sample of 105 quasars observed at multiple epochs by the Sloan Digital Sky Survey. We find a strong correlation between the change in the C IV line flux and the change in the line width, but no correlations between the change in flux and changes in line center and skewness. The relation between line flux change and line width change is consistent with a model in which a broad line base varies with greater amplitude than the line core. The objects studied here are more luminous and at higher redshift than those normally studied for variability, ranging in redshift from 1.65 to 4.00 and in absolute r-band magnitude from roughly -24 to -28. Using moment analysis line-fitting techniques, we measure line fluxes, centers, widths and skewnesses for the C IV line at two epochs for each object. The well-known Baldwin Effect is seen for these objects, with a slope beta = -0.22. The sample has a median intrinsic Baldwin Effect slope of beta = -0.85; the C IV lines in these high-luminosity quasars appear to be less responsive to continuum variations than those in lower luminosity AGN. Additionally, we find no evidence for variability of the well known blueshift of the C IV line with respect to the low-ionization Mg II line in the highest flux objects, indicating that this blueshift might be useful as a measure of orientation.
We discuss the effects of an enhanced interstellar radiation field (ISRF) on the observables of protostellar cores in the Orion cloud region. Dust radiative transfer is used to constrain the envelope physical structure by reproducing SCUBA 850 micron emission. Previously reported 13CO, C17O and H2CO line observations are reproduced through detailed Monte Carlo line radiative transfer models. It is found that the 13CO line emission is marginally optically thick and sensitive to the physical conditions in the outer envelope. An increased temperature in this region is needed in order to reproduce the 13CO line strengths and it is suggested to be caused by a strong heating from the exterior, corresponding to an ISRF in Orion 10^3 times stronger than the "standard" ISRF. The typical temperatures in the outer envelope are higher than the desorption temperature for CO. The C17O emission is less sensitive to this increased temperature but rather traces the bulk envelope material. The data are only fit by a model where CO is depleted, except in the inner and outermost regions where the temperature increases above 30-40 K. The fact that the temperatures do not drop below approximately 25 K in any of the envelopes whereas a significant fraction of CO is frozen-out suggest that the interstellar radiation field has changed through the evolution of the cores. The H2CO lines are successfully reproduced in the model of an increased ISRF with constant abundances of 3-5x10^{-10}.
Goldreich-Sridhar model of incompressible turbulence provides an elegant approach to describing strong MHD turbulence. It relies on the fact that interacting Alfvenic waves are independent and have random polarization. However, in case of strong interaction a spontaneous local assymetry can arise. We used direct numerical simulations to show that polarization alignment occurs and it grows larger at smaller scales. This effect would lead to a shallower spectrum and stronger anisotropy. Even small changes in these two properties will have important astrophysical consequencies, e.g. for the cosmic ray physics.
We use the formalism of finite-temperature field theory to study the interactions of ultra-high energy (UHE) cosmic neutrinos with the thermal background of relic neutrinos. From the imaginary part of the neutrino self-energy, calculated in terms of the Z boson propagator near the resonance, we derive general expressions for the UHE neutrino transmission probability. This allows us to take into account the thermal effects introduced by the momentum distribution of the relic neutrinos. We compare our results with the approximate expressions existing in the literature and discuss the influence of thermal effects on the absorption dips in the context of realistic UHE neutrino fluxes and favoured neutrino mass schemes.
The study of Be in stars of differing metal content can elucidate the formation mechanisms and the Galactic chemical evolution of Be. We have obtained high-resolution, high signal-to-noise spectra of the resonance lines of Be II in eight stars with the high-dispersion spectrograph (HDS) on the Subaru 8.2 m telescope on Mauna Kea. Abundances of Be have been determined through spectrum synthesis. The stars with [Fe/H] values > -1.1 conform to the published general trend of Be vs Fe. We have confirmed the high Be abundance in HD 94028 and have found a similarly high Be abundance in another star, HD 132475, at the same metallicity: [Fe/H] = -1.5. These two stars are 0.5 - 0.6 dex higher in Be than the Be-Fe trend. While that general trend contains the evidence for a Galaxy-wide enrichment in Be and Fe, the higher Be abundances in those two stars indicates local Be enrichments. Possible enrichment mechanisms include hypernovae and multiple supernova explosions contained in a superbubble. The star G 64-37 has [Fe/] = -3.2; we have determined its Be abundance to look for evidence of a Be plateau. It's Be abundance appears to extend the Be-Fe trend to lower Fe abundances without any evidence for a plateau as had been indicated by a high Be abundance in another very metal-poor star, G 64-12. Although these two stars have similar Be abundances within the errors, it could be that their different Be values may be indicating that a Be dispersion exists even at the lowest metallicities.
The effect of the length of inflation on the power spectra of scalar and tensor perturbations is estimated using the power-law inflation model with a scale factor of a(t) = t^q. Considering various pre-inflation models with radiation-dominated or scalar matter-dominated periods before inflation in combination with two matching conditions, the temperature angular power spectrum (TT) and temperature-polarization cross-power spectrum (TE) are calculated and a likelihood analysis is performed. It is shown that the discrepancies between the Wilkinson Microwave Anisotropy Probe (WMAP) data and the LCDM model, such as suppression of the spectrum at l = 2,3 and oscillatory behavior, may be explained by the finite length of inflation model if the length of inflation is near 60 e-folds and q > 300. The proposed models retain similar values of chi^2 to that achieved by the LCDM model with respect to fit to the WMAP data, but display different characteristics of the angular TE power spectra at l < 20.
We report a discovery of the Zeeman resolved spectral lines, corresponding to the extremely large magnetic field modulus <Bs>=17.5 kG, in the cool Ap star HD 178892. The mean longitudinal field of this star reaches 7.5 kG, and its rotational modulation implies the strength of the dipolar magnetic component Bp>=23 kG. We have revised rotation period of the star using the All Sky Automated Survey photometry and determined P=8.2478 d. Rotation phases of the magnetic and photometric maxima of the star coincide with each other. We obtained Geneva photometric observation of HD 178892 and estimated Teff=7700+/-250 K using photometry and the hydrogen Balmer lines. Preliminary abundance analysis reveals abundance pattern typical of rapidly oscillating Ap stars.
We review the observed properties of extremely hot hydrogen-deficient post-AGB stars of spectral type [WC] and PG1159. Their H-deficiency is probably caused by a (very) late helium-shell flash or a AGB final thermal pulse, laying bare interior stellar regions which are usually kept hidden below the hydrogen envelope. Thus, the photospheric element abundances of these stars allow to draw conclusions about details of nuclear burning and mixing processes in the precursor AGB stars. We summarize the state-of-the-art of stellar evolution models which simulate AGB evolution and the occurrence of a late He-shell flash. We compare predicted element abundances to those determined by quantitative spectral analyses performed with advanced non-LTE model atmospheres. A good qualitative and quantitative agreement is found. Future work can contribute to an even more complete picture of the nuclear processes in AGB stars.
We present the latest results of the Clusters in the Zone of Avoidance (CIZA) survey, which is mapping the large-scale matter distribution behind the Milky Way by performing the first systematic search for X-ray luminous galaxy clusters at low Galactic latitudes. The survey's approach, which uses X-ray emission to locate cluster candidates, overcomes the problems faced by optically selected cluster surveys which have traditionally avoided this region of the sky due to the severe extinction present along the Galactic plane. We here present the second flux-limited CIZA cluster catalog containing 60 X-ray luminous galaxy clusters, 88% of which are new discoveries. We also examine the degree to which known superclusters extend into the Zone of Avoidance and highlight newly discovered structures which have previously gone unnoticed. We show that the survey has found far fewer rich clusters in the Great Attractor region than would be expected given the region's proposed mass. Instead, we find a significant increase in the number of clusters behind the Great Attractor, with the most notable being an association of clusters near the Shapley supercluster. We propose these clusters trace an extension of the large-scale filament network in which the Shapley concentration is embedded. We also highlight an association of clusters near the Galactic anticenter, which is the first supercluster found to be completely hidden by the Milky Way. Our finding of a less massive Great Attractor and the detection of significant structures behind the complex supports studies which suggest the motion of nearby galaxies, including that of the Local Group, is due, in part, to a large-scale bulk flow which is induced by overdensities beyond the Great Attractor region.
I review here the results of the first RV survey for spectroscopic companions to very young brown dwarfs (BDs) and (very) low-mass stars in the ChaI star-forming cloud with UVES at the VLT. This survey studies the binary fraction in an as yet unexplored domain not only in terms of primary masses (substellar regime) and ages (a few Myr) but also in terms of companion masses (sensitive down to planetary masses) and separations (< 1 AU). The UVES spectra obtained so far hint at spectroscopic companions of a few Jupiter masses around one BD and around one low-mass star (M4.5) with orbital periods of at least several months. Furthermore, the data indicate a multiplicity fraction consistent with field BDs and stellar binaries for periods < 100 days.
A beginning end-point for galaxy motions within the 10-Mpc Local Volume is constructed by requiring a smooth distribution of (luminous) mass at the time of recombination, which is shown to be equivalent to a smooth Hubble flow at early times. It is found, by this purely kinematical method, that present peculiar motions are too small by a factor of at least several (and largely in the wrong direction) to have produced the observed structures within the age of the universe. Known dynamical effects are inadequate to remove the discrepancy. This result is different in origin from previously known ``cold flow'' problems. The simple dynamical picture often used within the Local Volume (for instance, in deriving masses through calculation of a zero-velocity surface) is thus called into question. The most straightforward explanation (though not the only possible) is that there exists a large quantity of baryonic matter in this region so far undetected, and unassociated with galaxies or groups.
It has been pointed out in Ref. [1] that in the presence of an external magnetic field the axion mass receives an electromagnetic contribution. We show that if a magnetic field with energy density larger than ~10^{-8} times the energy density of the Universe existed at temperatures of a few GeV, that contribution would be dominant and consequently the cosmic evolution of the axion field would change substantially. In particular, the expected axion relic abundance would be lowered, allowing a small relaxation of the present cosmological bound on the Peccei-Quinn constant.
A number of deep, wide-field, near-infrared surveys employing new infrared cameras on 4m-class telescopes are about to commence. These surveys have the potential to determine the fraction of luminous dust-obscured quasars that may have eluded surveys undertaken at optical wavelengths. In order to understand the new observations it is essential to make accurate predictions of surface densities and number-redshift relations for unobscured quasars in the near-infrared based on information from surveys at shorter wavelengths. The accuracy of the predictions depends critically on a number of key components. The commonly used single power-law representation for quasar SEDs is inadequate and the use of an SED incorporating the upturn in continuum flux at lambda~12000 A is essential. The presence of quasar host galaxies is particularly important over the restframe wavelength interval 8000 < lambda < 16000 A and we provide an empirical determination of the magnitude distribution of host galaxies using a low redshift sample of quasars from the SDSS DR3 quasar catalogue. A range of models for the dependence of host galaxy luminosity on quasar luminosity is investigated, along with the implications for the near-infrared surveys. Even adopting a conservative model for the behaviour of host galaxy luminosity the number counts for shallow surveys in the K-band increase by a factor of two. The degree of morphological selection applied to define candidate quasar samples in the near-infrared is found to be an important factor in determining the fraction of the quasar population included in such samples.
Over much of the initial mass function, stars are destined to become luminous and cool red giants. They may then be able to produce dust in an atmosphere which has been elevated by strong radial pulsations, and hence drive a wind. The amount of mass that is lost in this way can be a very significant fraction of the stellar mass, and especially in the case of intermediate-mass stars it is highly enriched. The delay between a star's birth and its feedback into the environment varies from several million years for massive stars to almost the age of the Universe for the least massive red giants we see today. I here present a review on the metallicity dependence of red giant winds. I show that recent measurements not only confirm theoretical expectations, but also admonish of common misconceptions with implications for feedback at low initial metallicity.
The physical nature of the currently observed dark energy in the universe is completely unclear, and many different theoretical models co-exist. Nevertheless, if dark energy is produced by vacuum fluctuations then there is a chance to probe some of its properties by simple laboratory tests based on Josephson junctions. These electronic devices can be used to perform `vacuum fluctuation spectroscopy', by directly measuring a noise spectrum induced by vacuum fluctuations. One would expect to see a cutoff near 1.7 THz in the measured power spectrum, provided the new physics underlying dark energy couples to electric charge. The effect exploited by the Josephson junction is a subtile nonlinear mixing effect and has nothing to do with the Casimir effect or other effects based on van der Waals forces. A Josephson experiment of the suggested type will now be built, and we should know the result within the next 3 years.
We report on ESO--VLT adaptive optics imaging of one radio-loud quasar at z $\sim$ 3. In spite of the large distance of the object we are able to detect its surrounding extended nebulosity the properties of which are consistent with an underlying massive galaxy of M$_K$ $\sim$ --27 and effective radius R$_e$ = 7 kpc. As far as we know this is the clearest detection of a radio loud quasar host at high redshift. The host luminosity is indicative of the existence of massive spheroids even at this early cosmic epoch. The host luminosity is about 1 magnitude fainter than the expected value based on the average trend of the host galaxies of RLQ at lower redshift. The result, which however is based on a single object, suggests that at z $\sim$ 3 there is a deviation from a luminosity--redshift dependence regulated only by passive evolution.
Below the knee in the cosmic ray spectrum, balloon and spacecraft experiments offer the capability of direct composition and energy measurements on the primary particles. A major difficulty is obtaining enough exposure to extend the range of direct measurements sufficiently high in energy to permit overlap with ground-based observations. Presently, balloon and space measurements extend only up to ~100 TeV, well below the range of ground-based experiments. The prospect of Ultra-Long Duration Balloon missions offers the promise of multiple long flights that can build up exposure. The status of balloon measurements to measure the high energy proton and nuclear composition and spectrum is reviewed, and the statistical considerations involved in searching for a steepening in the spectrum are discussed. Given the very steeply falling spectrum, it appears unlikely that balloon experiments will be able to extend the range of direct measurements beyond 1000 TeV any time in the near future. Especially given the recent suggestions from KASCADE that the proton spectrum steepens only at 4000-5000 TeV, the chance of detecting the knee with direct measurements of protons to iron on balloons is not likely to occur without significant increases in the payload and flight duration capabilities of high altitude balloons.
This paper studies the dynamical evolution of young stellar clusters with $N$ = 100 - 1000 members. We use N-body simulations to explore how evolution depends on system size $N$ and the initial conditions. Motivated by recent observations of extremely young systems, this study compares subvirial and virial starting states. Multiple realizations of equivalent cases (100 simulations per case) are used to build up a robust statistical description of these systems, e.g., distributions of closest approaches, mass profiles, and distributions of radial locations. These results provide a framework from which to assess the effects of these clusters on star and planet formation. The distributions of radial positions are used in conjunction with distributions of FUV luminosities (also calculated here) to determine the radiation exposure of circumstellar disks. The distributions of closest approaches are used in conjunction with scattering cross sections (calculated here from $10^5$ scattering experiments) to determine the probability of solar system disruption. We also use the nearby cluster NGC 1333 as a test case. Our main conclusion is that clusters in this size range have only a modest effect on forming planetary systems: Interaction rates are low so that the typical solar system experiences a single encounter within 1000 AU. Radiation exposure is low, with median FUV flux $G_0$ = 900, so that photoevaporation of disks is only important beyond 30 AU. Given the low interaction rates and modest radiation levels, we suggest that solar system disruption is a rare event in these clusters.
We use line-of-sight velocity information on the filamentary emission-line nebula of NGC 1275 to infer a dynamical model of the nebula's flow through the surrounding intracluster gas. We detect outflowing gas and flow patterns that match simulations of buoyantly rising bubbles from which we deduce that some of the nebula filaments have been drawn out of NGC 1275. We find a radial gradient of the ratio [NII]6584/H-alpha which may be due to a variation in metallicity, interactions with the surrounding intracluster medium or a hardening of the excitation mechanism. We find no preferred spatial correlation of stellar clusters within the filaments and there is a notable lack of [OIII]5007 emission, therefore it is unlikely that the filaments are ionized by stellar UV.
We study the evolution of the fine-structure constant, $\alpha$, induced by non-linear density perturbations in the context of the simplest class of quintessence models with a non-minimal coupling to the electromagnetic field, in which the two available free functions (potential and gauge kinetic function) are Taylor-expanded up to linear order. We show that the results obtained using the spherical infall model for an infinite wavelength inhomogeneity are inconsistent with the results of a local linearized gravity study and we argue in favour of the second approach. We also discuss recent claims that the value of $\alpha$ inside virialised regions could be significantly different from the background one on the basis of these findings.
We investigate models that suggest that the vacuum energy decays into cold dark matter (CDM) and into a homogeneous distribution of a thermalized cosmic microwave background (CMB), which is characteristic of an adiabatic vacuum energy decay into photons. We show that the density fluctuation spectrum obtained from the CMB and galaxy distribution data agreement put strong limits on the rate of decay of the vacuum energy. A vacuum energy decaying into CDM increases its total density \rho, diluting (\delta \rho/\rho)^2. The observed temperature fluctuations of the CMB photons (\delta T/T)^2 are approximately proportional to CDM density fluctuations (\delta \rho/\rho)^2. In both case, when evaluating (\delta \rho/\rho)^2 at the recombination era, its present measured value must be increased by a factor F. Since the (\delta \rho/\rho)^2 derived from the CMB and galaxy distribution data agree to \sim 10%, the maximum value for F is F_{\rm max}\aprox 1.1. Our results indicate that the rate of the decay of the vacuum energy into CDM and CMB photons are extremely small.
We performed a spectroscopic study of the 1809 keV gamma-ray line from 26Al decay in the Galaxy using the SPI imaging spectrometer with its high-resolution Ge detector camera on the INTEGRAL observatory. We analyzed observations of the first two mission years, fitting spectra from all 7130 telescope pointings in narrow energy bins to models of instrumental background and the 26Al sky. Instrumental background is estimated from independent tracers of cosmic-ray activation. The shape of the 26Al signal is compared to the instrumental response to extract the width of the celestial line. We detect the 26Al line at \~16sigma significance. The line is broadened only slightly, if at all; we constrain the width to be below 2.8 keV (FWHM, 2 sigma). The average Doppler velocities of 26Al at the time of its decay in the interstellar medium (decay time~1.04 My) therefore are probably around 100 km/s, in agreement with expectations from Galactic rotation and interstellar turbulence. The flux and spatial distribution of the emission are found consistent with previous observations. The derived amount of 26Al in the Galaxy is 2.8 (+/-0.8) M_solar.
Angular differential imaging is a high-contrast imaging technique that reduces quasi-static speckle noise and facilitates the detection of nearby companions. A sequence of images is acquired with an altitude/azimuth telescope while the instrument field derotator is switched off. This keeps the instrument and telescope optics aligned and allows the field of view to rotate with respect to the instrument. For each image, a reference PSF is constructed from other appropriately-selected images of the same sequence and subtracted to remove quasi-static PSF structure. All residual images are then rotated to align the field and are combined. Observed performances are reported for Gemini North data. It is shown that quasi-static PSF noise can be reduced by a factor \~5 for each image subtraction. The combination of all residuals then provides an additional gain of the order of the square root of the total number of acquired images. A total speckle noise attenuation of 20-50 is obtained for one-hour long observing sequences compared to a single 30s exposure. A PSF noise attenuation of 100 was achieved for two-hour long sequences of images of Vega, reaching a 5-sigma contrast of 20 magnitudes for separations greater than 8". For a 30-minute long sequence, ADI achieves 30 times better signal-to-noise than a classical observation technique. The ADI technique can be used with currently available instruments to search for ~1MJup exoplanets with orbits of radii between 50 and 300 AU around nearby young stars. The possibility of combining the technique with other high-contrast imaging methods is briefly discussed.
(abridged)We present a model which predicts inflation without the presence of inflaton fields, based on the \epsilon R^2 and Starobinsky models. It links the above models to the observable universe, in particular, to the ratio r of tensor to scalar fluctuations. In our model, we assume the existence of particles with the mass M that have a long decay time. These particles which were gravitationally produced \sim 60e-folds before the end of inflation produced the nearly scale invariant scalar density fluctuations which are observed. Gravitational waves (tensor fluctuations) were also produced at this epoch. The ratio of tensor to scalar fluctuations r (which are to be measured in the near future to good accuracy) determines M, which together with H_0, determine the time at the end of inflation, t_end. At t_end, the Hubble parameter begins to oscillate rapidly, gravitationally producing the bulk of the M particles, which we identify with the matter content of the universe today. The time required for the universe to dissipate its vacuum energy into M particles is found to be t_dis \simeq 6M_Pl^2/M^3. We assume that the time t_RH, (called the reheating time) needed for the M particles to decay into relativistic particles, is very much greater than that necessary to create the M particles, t_dis. From the ratio f\equiv t_dis/t_RH and g_\ast (the total number of degrees of freedom of the relativistic particles) we can, then, evaluate the maximum temperature of the universe, T_max, and the reheat temperature, T_RH, at t_RH. Our model, thus, predicts M, t_dis, t_end, T_max, T_RH, t_max, and t_RH as a function of r, f, and g_\ast (and to a weaker extent the particle content of the vacuum near the Planck epoch).
We study the frequency and angular dependences of Cherenkov radio pulses originated by the excess of electrons in electromagnetic showers in different dense media. We develop a simple model to relate the main characteristics of the electric field spectrum to the properties of the shower such as longitudinal and lateral development. This model allows us to establish the scaling of the electric field spectrum with the properties of the medium such as density, radiation length, Moliere radius, critical energy and refraction index. We normalize the predictions of the scaling relations to the numerical results obtained in our own developed GEANT4-based Monte Carlo simulation, and we give a unified parameterization of the frequency spectrum and angular distribution of the electric field in ice, salt, and the lunar regolith, in terms of the relevant properties of the media. Our parameterizations are valid for electromagnetic showers below the energy at which the Landau-Pomeranchuk-Migdal effect starts to be relevant in these media. They also provide an approximate estimate of radio emission in hadronic showers induced by high energy cosmic rays or neutrinos.
(abridged) We have investigated 136 Chandra extragalactic sources without broad optical emission lines, including 93 galaxies with narrow emission lines (NELG) and 43 with only absorption lines (ALG). Based on fx/fo, Lx, X-ray spectral hardness and optical emission line diagnostics, we have conservatively classified 36 normal galaxies (20 spirals and 16 ellipticals) and 71 AGNs. We found no statistically significant evolution in Lx/LB, within the limited z range. We have built log(N)-log(S), after correcting for completeness based on a series of simulations. The best-fit slope is -1.5 for both S and B energy bands, which is considerably steeper than that of the AGN-dominated cosmic background sources, but slightly flatter than the previous estimate, indicating normal galaxies will not exceed the AGN population until fx ~ 2 x 10-18 erg s-1 cm-2 (a factor of ~5 lower than the previous estimate). A group of NELGs appear to be heavily obscured in X-rays, i.e., a typical type 2 AGN. After correcting for intrinsic absorption, their X-ray luminosities could be Lx > 10^44 erg s-1, making them type 2 quasar candidates. While most X-ray luminous ALGs (XBONG - X-ray bright, optically normal galaxy candidates) do not appear to be significantly absorbed, we found two heavily obscured objects, which could be as luminous as an unobscured broad-line quasar. Among 43 ALGs, we found two E+A galaxy candidates with strong Balmer absorption lines, but no [OII] line. The X-ray spectra of both galaxies are soft and one of them has a nearby close companion galaxy, supporting the merger/interaction scenario rather than the dusty starburst hypothesis.
We present phase-resolved low resolution $JHK$ and higher resolution $K$-band spectroscopy of the polar VV Pup. All observations were obtained when VV Pup was in a low accretion state having a K magnitude near 15. The low resolution observations reveal cyclotron emission in the $J$ band during some phases, consistent with an origin near the active 30.5 MG pole on the white dwarf. The secondary in VV Pup appears to be a normal M7V star and we find that the $H$ and $K$ band fluxes are entirely due to this star at all orbital phases during the low accretion state. We use our higher resolution Keck spectroscopy to produce the first $K$-band radial velocity curve for VV Pup. Our orbital solution yields $K_2$=414$\pm27$ km sec$^{-1}$ and leads to mass estimates of M$_1$=0.73$\pm$0.05 M$_{\odot}$ and M$_2$=0.10$\pm$0.02 M$_{\odot}$. We find that the mass accretion rates during the normal low states of the polars VV Pup, EF Eri, and EQ Cet are near 10$^{-13}$ M$_{\odot}$ yr$^{-1}$. The fact that \.M is not zero in low state polars indicates active secondary stars in these binary systems, including the sub-stellar donor star present in EF Eri.
We present a high-resolution, high signal-to-noise optical spectrum of the afterglow of GRB 051111 obtained with the HIRES spectrograph on the Keck I 10-m telescope. The spectrum exhibits three redshifted absorption systems with the highest, at z=1.54948, arising in the GRB host galaxy. While the Ly-alpha feature is outside the range of our spectrum, the high column density of weakly-depleted Zn suggests that the host is a damped Lyman-alpha system with N(HI)>10^21(Z/Z_sun)^-1. The bulk of the gas (>80%) is confined to a narrow velocity range of |v|<30 km/s exhibiting strong dust depletion of refractory elements such as Fe and Cr. The depletion pattern is similar to that observed in warm disk clouds of the Milky Way. We also detect absorption from all ground-level fine-structure states of FeII, the first such example in a QSO-DLA or GRB-absorption spectrum, which indicate conditions that are consistent with the "warm disk" depletion pattern. The absorption profiles of FeII and MgII extend over several hundred km/s, with a depletion pattern that more closely resembles that of QSO-DLAs, suggesting that the sight line to GRB 051111 probes the halo of the host galaxy in addition to the dense disk. Thus, detailed diagnostics of the interstellar medium of GRB host galaxies continue to provide insight into regions which are generally missed in quasar surveys.
We have conducted a search for giant pulses from four millisecond pulsars using the 100m Green Bank Telescope. Coherently dedispersed time-series from PSR J0218+4232 were found to contain giant pulses of very short intrinsic duration whose energies follow power-law statistics. The giant pulses are in phase with the two minima of the radio integrated pulse profile but are phase aligned with the peaks of the X-ray profile. Historically, individual pulses more than 10-20 times the mean pulse energy have been deemed to be ``giant pulses''. As only 4 of the 155 pulses had energies greater than 10 times the mean pulse-energy, we argue the emission mechanism responsible for giant pulses should instead be defined through: (a) intrinsic timescales of microsecond or nanosecond duration; (b) power-law energy statistics; and (c) emission occurring in narrow phase-windows coincident with the phase windows of non-thermal X-ray emission. Four short-duration pulses with giant-pulse characteristics were also observed from PSR B1957+20. As the inferred magnetic fields at the light cylinders of the millisecond pulsars that emit giant pulses are all very high, this parameter has previously been considered to be an indicator of giant pulse emissivity. However, the frequency of giant pulse emission from PSR~B1957+20 is significantly lower than for other millisecond pulsars that have similar magnetic fields at their light cylinders. This suggests that the inferred magnetic field at the light cylinder is a poor indicator of the rate of emission of giant pulses.
I describe electromagnetic model of gamma ray bursts and contrast its main properties and predictions with hydrodynamic fireball model and its magnetohydrodynamical extension. The electromagnetic model assumes that rotational energy of a relativistic, stellar-mass central source (black-hole--accretion disk system or fast rotating neutron star) is converted into magnetic energy through unipolar dynamo mechanism, propagated to large distances in a form of relativistic, subsonic, Poynting flux-dominated wind and is dissipated directly into emitting particles through current-driven instabilities. Thus, there is no conversion back and forth between internal and bulk energies as in the case of fireball model. Collimating effects of magnetic hoop stresses lead to strongly non-spherical expansion and formation of jets. Long and short GRBs may develop in a qualitatively similar way, except that in case of long bursts ejecta expansion has a relatively short, non-relativistic, strongly dissipative stage inside the star. Electromagnetic and fireball models (as well as strongly and weakly magnetized fireballs) lead to different early afterglow dynamics, before deceleration time. Finally, I discuss the models in view of latest observational data in the Swift era.
The XBootes Survey is a 5-ks Chandra survey of the Bootes Field of the NOAO Deep Wide-Field Survey (NDWFS). This survey is unique in that it is the largest (9.3 deg^2), contiguous region imaged in X-ray with complementary deep optical and near-IR observations. We present a catalog of the optical counterparts to the 3,213 X-ray point sources detected in the XBootes survey. Using a Bayesian identification scheme, we successfully identified optical counterparts for 98% of the X-ray point sources. The optical colors suggest that the optically detected galaxies are a combination of z<1 massive early-type galaxies and bluer star-forming galaxies whose optical AGN emission is faint or obscured, whereas the majority of the optically detected point sources are likely quasars over a large redshift range. Our large area, X-ray bright, optically deep survey enables us to select a large sub-sample of sources (773) with high X-ray to optical flux ratios (f_x/f_o>10). These objects are likely high redshift and/or dust obscured AGN. These sources have generally harder X-ray spectra than sources with 0.1<f_x/f_o<10. Of the 73 X-ray sources with no optical counterpart in the NDWFS catalog, 47 are truly optically blank down to R~25.5 (the average 50% completeness limit of the NDWFS R-band catalogs). These sources are also likely to be high redshift and/or dust obscured AGN.
Hypervelocity stars have been recently discovered in the outskirts of galaxies, such as the unbound star in the Milky Way halo, or the three anomalously fast intracluster planetary nebulae (ICPNe) in the Virgo Cluster. These may have been ejected by close 3-body interactions with a binary supermassive black hole (SMBBH), where a star which passes within the semimajor axis of the SMBBH can receive enough energy to eject it from the system. Stars ejected by SMBBHs may form a significant sub-population with very different kinematics and mean metallicity than the bulk of the intracluster stars. The number, kinematics, and orientation of the ejected stars may constrain the mass ratio, semimajor axis, and even the orbital plane of the SMBBH. We investigate the evolution of the ejected debris from a SMBBH within a clumpy and time-dependent cluster potential using a high resolution, self-consistent cosmological N-body simulation of a galaxy cluster. We show that the predicted number and kinematic signature of the fast Virgo ICPNe is consistent with 3-body scattering by a SMBBH with a mass ratio $10:1$ at the center of M87.
We study the phenomenology of cosmic-rays (CRs) at the galactic/extragalactic transition, focusing on two opposite models for the composition of the extragalactic (EG) component. Model A assumes a mixed source composition, with nuclear abundances similar to that of the low-energy CRs, while model B assumes that the EG sources accelerate only protons. We study the limits within which both scenarios can reproduce the high-energy CR data, from the point of view of the energy spectrum, each with its own interpretation of the ankle (as the GCR/EGCR transition for model A; as a pair-production dip for model B) and its own source spectrum (in E?2.3 model A; in E?2.6 or E?2.7 for model B). We then focus on the energy evolution of the high-energy CR composition and compare our predictions for two main composition-related observables (Xmax and hlnAi) with the available data, to conclude that model A is currently favoured. Uncertainties are discussed and distinctive features of the two models are identified, which should allow one to distinguish between the models in the near future when more precise measurements are available with higher-statistics experiments.
We present line strengths in the bulges and inner disks of 38 galaxies in the local universe, including several galaxies whose bulges were previously identified as being disk-like in their colors or kinematics, to see if their spectral properties reveal evidence for secular evolution. We find that red bulges of all Hubble types are similar to luminous ellipticals in their central stellar populations. They have large luminosity-weighted ages, metallicities, and alpha/Fe ratios. Blue bulges can be separated into a metal-poor class that is restricted to late-types with small velocity dispersion and a young, metal-rich class that includes all Hubble types and velocity dispersions. Luminosity-weighted metallicities and alpha/Fe ratios are sensitive to central velocity dispersion and maximum disk rotational velocity. Red bulges and ellipticals follow the same scaling relations. We see differences in some scaling relations between blue and red bulges and between bulges of barred and unbarred galaxies. Most bulges have decreasing metallicity with increasing radius; galaxies with larger central metallicities have steeper gradients. Where positive age gradients (with the central regions being younger) are present, they are invariably in barred galaxies. The metallicities of bulges are correlated with those of their disks. While this and the differences between barred and unbarred galaxies suggest that secular evolution cannot be ignored, our results are generally consistent with the hypothesis that mergers have been the dominant mechanism responsible for bulge formation.
We use new deep NIR and MIR observations to analyze the 850 um image of the GOODS-N region around the HDF-N. We show that much of the 850 um background is picked out by sources with H(AB) or 3.6 um(AB)<23.25 (1.8 uJy). These sources contribute an 850 um background of 24+/-2 Jy/deg^2. This is a much higher fraction of the measured background (31-45 Jy/deg^2) than is found with current 20 cm or 24 um samples. Roughly one-half of these NIR-selected sources have spectroscopic identifications, and we can assign robust photometric redshifts to nearly all of the remaining sources using their UV to MIR SEDs. We use the redshift and spectral type information to show that a large fraction of the background light comes from sources with z=0-1.5 and that the sources responsible have intermediate spectral types. The redshift distribution of the NIR-selected 850 um light lies well below that of the much smaller amount of light traced by the more luminous, radio-selected submillimeter sources. We therefore require a revised star-formation history with a lower star-formation rate at high redshifts. We use a stacking analysis of the 20 cm light in the NIR sample to show that the star-formation history of the total 850 um sample peaks at z~1. This is considerably lower than the currently accepted value of z~2-3.
We present metallicities, [Fe/H], and elemental abundance ratios, [X/Fe], for
a sample of 24 Cepheids in the outer Galactic disk. The sample have
Galactocentric distances covering 12 < R_GC (kpc) < 17.2 making them the most
distant Galactic Cepheids upon which detailed abundance analyses have been
performed. We find sub-solar ratios of [Fe/H] and overabundances of [alpha/Fe],
[La/Fe], and [Eu/Fe]. All abundance ratios exhibit a dispersion that exceeds
the measurement uncertainties. As seen in our previous studies of old open
clusters and field giants, enhanced ratios of [alpha/Fe] and [Eu/Fe] reveal
that recent star formation has taken place in the outer disk.
The outer disk Cepheids appear to exhibit a bimodal distribution for [Fe/H]
and [alpha/Fe]. Most of the Cepheids continue the trends with Galactocentric
distance exhibited by Andrievsky's larger Cepheid sample and we refer to these
stars as the "Galactic Cepheids". A minority of the Cepheids show considerably
lower [Fe/H] and higher [alpha/Fe] and we refer to these stars as the "Merger
Cepheids". One signature of a merger event would be composition differences
between the "Galactic" and "Merger" Cepheids. The Cepheids satisfy this
requirement and we speculate that the distinct compositions suggest that the
"Merger Cepheids" may have formed under the influence of significant merger or
accretion events. The short lifetimes of the Cepheids reveal that the merger
event may be on-going with the Monoceros ring and Canis Major galaxy being
possible merger candidates. [Abridged]
Galaxy clusters are the most massive gravitational structures in the universe, and also provide unique probes of its expansion. Their strong X-ray emission and unique cosmic microwave background scattering signature afford an independent method to obtain distances based on the physics of ionized plasmas. In this article we determine the distance to 39 clusters of galaxies in the redshift range 0.14<z<0.89 using Chandra X-ray data and radio observations from the OVRO and BIMA interferometric arrays. We analyze the plasma and dark matter distribution in clusters using a hydrostatic equilibrium model that accounts for radial variations in density, temperature and abundance, and quantify the statistical and systematic errors of this method. The analysis is performed via a Markov chain Monte Carlo technique that provides simultaneous estimation of all model parameters. We measure a Hubble constant of H_0= 77.1 +3.8-3.4 +10.0-8.0 Km/s/Mpc (statistical followed by systematic uncertainty at 68% confidence) for an Omega_M=0.3, Omega_Lambda=0.7 cosmology. We also analyze the data using an isothermal beta model that is free of the hydrostatic equilibrium assumption, and find H_0=73.4 +4.5-3.8 +9.5-7.6 Km/s/Mpc; to avoid effects from cool cores in clusters, we repeated this analysis excluding the central 100 kpc from the X-ray data, and find H_0= 77.2 +4.8-4.2 +10.1-8.2 km/s/Mpc (statistical followed by systematic uncertainty at 68% confidence). The consistency between the models illustrates the relative insensitivity of SZE/X-ray determinations of H_0 to the details of the cluster model. Our determination of the Hubble parameter in the distant universe agrees with the recent measurement from the Hubble Space Telescope key project that probes the nearby universe.
This was my final paper for the AST 308 Galaxies class at Michigan State University. Using many sources I was able to compile a moderate amount of information concerning the evidence for, and the formation of Supermassive Black Holes.
The L1551 molecular cloud contains two small clusters of Class 0 and I protostars, as well as a halo of more evolved Class II and III YSOs, indicating a current and at least one past burst of star formation. We present here new, sensitive maps of 850 and 450 um dust emission covering most of the L1551 cloud, new CO J=2-1 data of the molecular cloud, and a new, deep, optical image of [SII] emission. No new Class 0/I YSOs were detected. Compact sub-millimetre emitters are concentrated in two sub-clusters: IRS5 and L1551NE, and the HL~Tauri group. Both stellar groups show significant extended emission and outflow/jet activity. A jet, terminating at HH 265 and with a very weak associated molecular outflow, may originate from LkHa 358, or from a binary companion to another member of the HL Tauri group. Several Herbig Haro objects associated with IRS5/NE were clearly detected in the sub-mm, as were faint ridges of emission tracing outflow cavity walls. We confirm a large-scale molecular outflow originating from NE parallel to that from IRS5, and suggest that the "hollow shell" morphology is more likely due to two interacting outflows. We confirm the presence of a prestellar core (L1551-MC) of mass 2-3 Mo north-west of IRS5. The next generation cluster may be forming in this core. The L1551 cloud appears cometary in morphology, and appears to be illuminated and eroded from the direction of Orion, perhaps explaining the multiple episodes of star formation in this cloud. The full paper (including figures) can be downloaded at this http URL, or viewed at this http URL
We describe a method of estimating the abundance of short-period extrasolar planets based on the results of a photometric survey for planetary transits. We apply the method to a 21-night survey with the 2.5m Isaac Newton Telescope of \~32000 stars in a ~0.5 deg by 0.5 deg square field including the open cluster NGC 7789. From the colour-magnitude diagram we estimate the mass and radius of each star by comparison with the cluster main sequence. We search for injected synthetic transits throughout the lightcurve of each star in order to determine their recovery rate, and thus calculate the expected number of transit detections and false alarms in the survey. We take proper account of the photometric accuracy, time sampling of the observations and criteria (signal-to-noise and number of transits) adopted for transit detection. Assuming that none of the transit candidates found in the survey will be confirmed as real planets, we place conservative upper limits on the abundance of planets as a function of planet radius, orbital period and spectral type.
We show that contrary to what is expected from 1D stationary model atmospheres, 3D hydrodynamical modeling predicts a considerable influence of convection on the spectral properties of late-type giants. This is due to the fact that convection overshoots into the formally stable outer atmospheric layers producing a notable granulation pattern in the 3D hydrodynamical models, which has a direct influence on the observable spectra and colors. Within the framework of standard 1D model atmospheres the average thermal stratification of the 3D hydro model can not be reproduced with any reasonable choice of the mixing length parameter and formulation of the turbulent pressure. The differences in individual photometric colors -- in terms of 3D versus 1D -- reach up to ~0.2 mag, or \Delta Teff~70K. We discuss the impact of full 3D hydrodynamical models on the interpretation of observable properties of late-type giants, briefly mentioning problems and challenges which need to be solved for bringing these models to a routine use within the astronomical community in 5-10 years from now.
To assess the current status in the theoretical modeling of the spectral properties of late-type giants, we provide a comparison of synthetic photometric colors of late-type giants (calculated with PHOENIX, MARCS and ATLAS model atmospheres) with observations, at [M/H]=0.0 and -2.0. Overall, there is a good agreement between observed and synthetic colors, and synthetic colors and published Teff-color relations, both at [M/H]=0.0 and -2.0. Deviations from the observed trends in Teff-color planes are generally within \pm 150K (or less) in the effective temperature range Teff=3500-4800K. Synthetic colors calculated with different stellar atmosphere models typically agree to ~100K, within a large range of effective temperatures and gravities. Some discrepancies are seen in the Teff-(B-V) plane below Teff~3800K at [M/H]=0.0, due to difficulties in reproducing the 'turn-off' to the bluer colors which is seen in the observed data at Teff~3600K. Note that at [M/H]=-2.0 effective temperatures given by the scale of Alonso et al. (1999) are generally lower than those resulting from other Teff-color relations based both on observed and synthetic colors.
The expectation of explaining cosmological observations without requiring new energy sources is forsooth worthy of investigation. In this letter, a new kind of Cardassian models, called exponential Cardassian models, for the late-time universe are investigated in the context of the spatially flat FRW universe scenario. We fit the exponential Cardassian models to current type Ia supernovae data and find they are consistent with the observations. Furthermore, we point out that the equation-of-state parameter for the effective dark fluid component in exponential Cardassian models can naturally cross the cosmological constant divide $w=-1$ that observations favor mildly without introducing exotic material that destroy the weak energy condition.
We present VLT/ISAAC near-infrared imaging of the host galaxies of 15 low luminosity quasars at 1 < z < 2. This work complements our studies to trace the cosmological evolution of the host galaxies of high luminosity quasars. The radio-loud (RLQ) and radio-quiet (RQQ) quasars have similar distribution of redshift and luminosity, and together the high and low luminosity quasars cover a large range of the quasar luminosity function. Both RLQ and RQQ hosts resemble massive inactive ellipticals undergoing passive evolution. However, RLQ hosts are systematically more luminous than RQQ hosts, as also found for the high luminosity quasars. The difference in the host luminosity remains the same from z = 2 to z = 0. For the entire set of quasars, we find a correlation between the nuclear and the host luminosities, albeit with a large scatter. The correlation is less apparent for the RQQs than for the RLQs.
We investigate the redshift dependence of the contribution to the extragalactic far-infrared/sub-millimeter background from galaxies detected by the Spitzer Space Telescope at 8um and 24um. Using seven-band optical to mid-infrared photometry, we estimate photometric redshifts for the Spitzer sources which appear to be mostly L* galaxies at a median redshift of z=1.0. These sources, extracted from deep 8um and 24um mosaics of the CUDSS 14-hour field with 5sigma limits of 5.8uJy and 70uJy respectively, exhibit significant 850um and 450um emission as observed by SCUBA. At 850um, after removing >=4sigma sources and those securely identified in our previous cross-matching paper, we measure stacked flux at the significance level of 4.4sigma and 2.9sigma from the full 8um and 24um galaxy catalogue respectively. At 450um, flux is detected from all 8um galaxies at the level of 3.5sigma, while there is no significant emission from the 24um galaxies. We find that the 850um flux is emitted almost exclusively at z>~1.3 from the Spitzer sources with 0.44mJy (4.7sigma) per 8um source and 0.51mJy (2.8sigma) per 24um source. This corresponds to a contribution of (16+/-3)% toward the 850um extra-galactic background from the 8um sources and (5.0+/-1.8)% from the 24um sources. At 450um, only the 8um sources within the redshift interval 1<z<2 exhibit significant emission with an average flux per source of 3.35mJy (3.0sigma). This is a contribution of (37+/-12)% to the 450um background.
GRB 050911, discovered by the Swift Burst Alert Telescope, was not seen 4.6 hr later by the Swift X-ray Telescope, making it one of the very few X-ray non-detections of a Gamma-Ray Burst (GRB) afterglow at early times. The gamma-ray light-curve shows at least three peaks, the first two of which (~T_0 - 0.8 and T_0 + 0.2 s, where T_0 is the trigger time) were short, each lasting 0.5 s. This was followed by later emission 10-20 s post-burst. The upper limit on the unabsorbed X-ray flux was 1.7 x 10^-14 erg cm^-2 s^-1 (integrating 46 ks of data taken between 11 and 18 September), indicating that the decay must have been rapid. All but one of the long bursts detected by Swift were above this limit at ~4.6 hr, whereas the afterglows of short bursts became undetectable more rapidly. Deep observations with Gemini also revealed no optical afterglow 12 hr after the burst, down to r=24.0 (5-sigma limit). We speculate that GRB 050911 may have been formed through a compact object (black hole-neutron star) merger, with the later outbursts due to a longer disc lifetime linked to a large mass ratio between the merging objects. Alternatively, the burst may have occured in a low density environment, leading to a weak, or non-existent, forward shock - the so-called 'naked GRB' model.
The detectors designed for Gaia, the next ESA space astrometry mission to be launched in 2011, will allow to observe repeatedly stars very close to Jupiter's limb. This will open a unique opportunity to test General Relativity by performing many Eddington-like experiments through the comparison between the pattern of a starfield observed with or without Jupiter. We have derived the main formulas relevant for the monopole and quadrupole light deflection by an oblate planet and developed a simulator to investigate the processing of the Gaia astrometric observation in the vicinity of the planet. The results show that such an experiment carried out with the Gaia data will provide a new fully independent determination of the PPN parameter gamma by means of differential astrometric measurements and, more importantly, for the first time will evidence the bending effect due to the quadrupole moment with a 3-sigma confidence level. Given the accuracy of the experiment for the monopole deflection, this will permit to test alternative modelling of the light bending by moving masses.
At which masses does the regime of globular clusters end and the one of dwarf
galaxies begin? And what separates these two classes of hot stellar systems?
We examine to what extend very massive (>10^7 Mo) young star clusters are
similar to their lower mass counter parts and to which degree they resemble
other objects in their mass regime (dwarf--globular transition objects (DGTOs),
ultra compact dwarf galaxies (UCDs), galaxy nuclei).
The comparison is performed by placing the recently observed very massive
young clusters onto known scaling relation defined by globular clusters (with
typical masses <10^6 Mo) and/or by hot stellar systems with sizes up to those
of giant galaxies.
The very massive (>10^{6.5-7} Mo) young clusters seem to show a mass--radius
relation compatible with the one defined by hot stellar systems of galaxy mass.
This, in turn, can explain their location on the other scaling relations
investigated. It contrasts with the behaviour of the less massive young
clusters and of globular clusters, which do not exhibit any mass-radius
relation. However, the behaviour of the most massive clusters is similar to
that of most other objects in that mass regime (10^6-10^8 Mo).
We show that the properties of young massive clusters are compatible with
other objects in the same mass regime such as DGTOs/UCDs. They present a
possible direct avenue of formation for those objects, which does not require
the transformation of a previously existing stellar system. Simulations and
observations support the possibility of the formation of such very massive
young clusters by early mergers of lower mass stellar clusters, which could
explain the emergence of a mass--radius relation.
We show that aperiodic and quasiperiodic variability of bright LMXBs - atoll
and Z- sources, on ~sec - msec time scales is caused primarily by variations of
the boundary layer luminosity. The accretion disk emission is less variable on
these time scales and its power density follows 1/f law, contributing to
observed flux variation at low frequencies and low energies only. The kHz QPOs
have the same origin as variability at lower frequencies - independent of the
nature of the "clock", the actual luminosity modulation takes place on the NS
surface. The boundary layer spectrum remains nearly constant during luminosity
variations and can be represented by the Fourier frequency resolved spectrum.
In the range of Mdot~(0.1-1)*Mdot_Edd it depends weakly on the global mass
accretion rate and in the limit Mdot~Mdot_Edd is close to Wien spectrum with
kT~2.4 keV. Its independence on the Mdot lends support to the suggestion by
Inogamov & Sunyaev (1999) that the boundary layer is radiation pressure
supported.
Based on the knowledge of the boundary layer spectrum we attempt to relate
the motion along the Z-track to changes of physically meaningful parameters.
Our results suggest that the contribution of the boundary layer to the observed
emission decreases along the Z-track from conventional ~50% on the horizontal
branch to a rather small number on the normal branch. This decrease can be
caused, for example, by obscuration of the boundary layer by the geometrically
thickened accretion disk at Mdot~Mdot_Edd. Alternatively, this can indicate
significant change of the structure of the accretion flow at Mdot~Mdot_Edd and
disappearance of the boundary layer as a distinct region of the significant
energy release associated with the NS surface.
We present a JHKL survey of the massive star forming region RCW 57 (NGC 3576) based on L-band data at 3.5 micron taken with SPIREX (South Pole Infrared Explorer), and 2MASS JHK data at 1.25-2.2 micron. This is the second of two papers, the first one concerning a similar JHKL survey of 30 Doradus. Colour-colour and colour-magnitude diagrams are used to detect sources with infrared excess. This excess emission is interpreted as coming from circumstellar disks, and hence gives the cluster disk fraction (CDF). Based on the CDF and the age of RCW 57, it is possible to draw conclusions on the formation and early evolution of massive stars. The infrared excess is detected by comparing the locations of sources in JHKL colour-colour and L vs. (K-L) colour-magnitude diagrams to the reddening band due to interstellar extinction. A total of 251 sources were detected. More than 50% of the 209 sources included in the diagrams have an infrared excess. Comparison with other JHKL surveys, including the results on 30 Doradus from the first paper, support a very high initial disk fraction (>80%) even for massive stars, although there is an indication of a possible faster evolution of circumstellar disks around high mass stars. 33 sources only found in the L-band indicate the presence of heavily embedded, massive Class I protostars. We also report the detection of diffuse PAHs emission throughout the RCW 57 region.
We build an accurate database of 5200 HCN and HNC rotation-vibration energy
levels, determined from existing laboratory data. 20~000 energy levels in the
Harris et al. (2002) linelist are assigned approximate quantum numbers. These
assignments, lab determined energy levels and Harris et al (2002) energy levels
are incorporated in to a new energy level list. A new linelist is presented, in
which frequencies are computed using the lab determined energy levels where
available, and the ab initio energy levels otherwise.
The new linelist is then used to compute new model atmospheres and synthetic
spectra for the carbon star WZ Cas. This results in better fit to the spectrum
of WZ Cas in which the absorption feature at 3.56 micron is reproduced to a
higher degree of accuracy than has previously been possible. We improve the
reproduction of HCN absorption features by reducing the abundance of Si to
[Si/H] = --0.5 dex, however, the strengths of the $\Delta v=2$ CS band heads
are over-predicted.
We have modelled the process of reionization in the high redshift Universe to determine the epoch of reionization. Reionization is driven by star formation in high redshift galaxies. We employ a semi-analytic model of galaxy formation to track the formation of these galaxies, their influence on the intergalactic medium (IGM) and the back-reaction of the IGM on further galaxy formation. This represents a more complete and physical modelling of reionization than has been conducted in the past. In particular, compared to our previous work our new calculations contain significant improvements in the modelling of the effects of reionization on the galaxy population and in the cooling model used to compute the rate at which galaxies form (our new model includes photoheating from a self-consistently computed ionizing background and also includes cooling due to molecular hydrogen). We find that reionization can be achieved by z~14-15 in a cosmological model motivated by the results from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite, consistent with the optical depth to reionization measured by WMAP. However, a cosmological model with a running spectral index is able to achieve reionization before z~9 only with very extreme assumptions about the physics of feedback at high redshifts. These results assume that all ionizing photons in galaxies are able to escape and ionize the IGM. If this is not the case then the redshift of reionization could be substantially reduced. Under the assumption that early star formation leads to the formation of very massive stars we find that extended periods of partial reionization and double reionizations can occur. Such models do not fully reionize until z~6-7 but obtain an optical depth which is consistent with the results from the WMAP satellite. (abridged)
We present 3D hydrodynamical simulations of ram pressure stripping of massive
disc galaxies in clusters. Studies of galaxies that move face-on have predicted
that in such a geometry the galaxy can lose a substantial amount of its
interstellar medium. But only a small fraction of galaxies is moving face-on.
Therefore, in this work we focus on a systematic study of the effect of the
inclination angle between the direction of motion and the galaxy's rotation
axis.
In agreement with some previous works, we find that the inclination angle
does not play a major role for the mass loss as long as the galaxy is not
moving close to edge-on. We can predict this behaviour by extending Gunn &
Gott's estimate of the stripping radius, which is valid for face-on geometries,
to moderate inclinations.
The inclination plays a role as long as the ram pressure is comparable to
pressures in the galactic plane, which can span two orders of magnitude. For
very strong ram pressures, the disc will be stripped completely, and for very
weak ram pressures, mass loss is negligible independent of inclination. We show
that in non-edge-on geometries the stripping proceeds remarkably similar. A
major difference between different inclinations is the degree of asymmetry
introduced in the remaining gas disc.
We demonstrate that the tail of gas stripped from the galaxy does not
necessarily point in a direction opposite to the galaxy's direction of motion.
Therefore, the observation of a galaxy's gas tail may be misleading about the
galaxy's direction of motion.
We study the red giant populations of two dE galaxies, AM 1339-445 and AM 1343-452, with the aim of investigating the number and luminosity of any upper asymptotic giant branch (AGB) stars present. The galaxies are members of the Centaurus A group (D~3.8 Mpc) and are classified as outlying (R~350 kpc) satellites of Cen A. The analysis is based on near-IR photometry for individual red giant stars, derived from images obtained with ISAAC on the VLT. The photometry, along with optical data derived from WFPC2 images retrieved from the HST science archive, enable us to investigate the stellar populations of the dEs in the vicinity of the red giant branch (RGB) tip. In both systems we find stars above the RGB tip, which we interpret as intermediate-age upper-AGB stars. The presence of such stars is indicative of extended star formation in these dEs similar to that seen in many, but not all, dEs in the Local Group. For AM 1339-445, the brightest of the upper-AGB stars have Mbol~-4.5 while those in AM 1343-452 have Mbol~-4.8 mag. These luminosities suggest ages of approximately 6.5+/-1 and 4+/-1 Gyr as estimates for the epoch of the last episode of significant star formation in these systems. In both cases the number of upper-AGB stars suggests that ~15% of the total stellar population is in the form of intermediate-age stars, considerably less than is the case for outlying dE satellites of the Milky Way such as Fornax and LeoI.
The effects of mass-varying neutrinos on cosmic microwave background (CMB) anisotropies and large scale structures (LSS) are studied. In these models, dark energy and neutrinos are coupled such that the neutrino masses are functions of the scalar field playing the role of dark energy. We begin by describing the cosmological background evolution of such a system. It is pointed out that, similar to models with a dark matter/dark energy interaction, the apparent equation of state measured with SNIa can be smaller than -1. We then discuss the effect of mass-varying neutrinos on the CMB anisotropies and the matter power spectrum. A suppression of power in the CMB power spectrum at large angular scales is usually observed. We give an explanation for this behaviour and discuss different couplings and quintessence potentials to show the generality of the results obtained. We perform a likelihood analysis using wide-ranging SNIa, CMB and LSS observations to assess whether such theories are viable. Treating the neutrino mass as a free parameter we find that the constraints on the coupling are weak, since CMB and LSS surveys give only upper bounds on the neutrino mass. However, fixing a priori the neutrino masses, we find that there is some evidence that the existence of such a coupling is actually preferred by current cosmological data over the standard LambdaCDM cosmology.
In this paper we derive a generic expression, which is valid for scales larger than Hubble radius and contains only the local terms, for the second order curvature perturbations for more than one field, provided the expansion is sourced by the energy density of a single field. As an application, motivated by our previous paper [1], we apply our formalism to two fields during preheating, where the inflaton oscillations are sourced by $\lambda\phi^4$ potential which is governing the expansion of the Universe. A second field $\sigma$, coupled to the inflaton through $g^2\phi^2\sigma^2$, is excited from the vacuum fluctuations. The excited modes of $\sigma$ amplify the super-Hubble isocurvature perturbations, which seed the second order curvature perturbations to give rise to a significantly large non-Gaussianity. Our results show that within 3 inflaton oscillations for a range of parameters, $1< g^2/\lambda < 3$, the non-Gaussianity parameter becomes: $f_{NL}\geq {\cal O}(1000)$, which is already ruled out by the current WMAP observation.
We study the possibility of using the entire probability distribution function (PDF) of the aperture mass Map and its related cumulative probability distribution function (CPDF) to obtain meaningful constraints on cosmological parameters. Deriving completely analytic expressions for the associated covariance matrices, we construct the Fisher matrix and use it to estimate the accuracy with which various cosmological parameters can be recovered from future surveys using such statistics. This formalism also includes the effect of various noises such as intrinsic ellipticity distribution of galaxies and finite survey size. The estimation errors are then compared with the ones derived from low order moments of the PDF (variance and skewness) to check how efficiently the high Map tail can be used to constrain cosmological parameters such as Omega_m, sigma_8 and dark energy equation of state w_de. We find that for future surveys such as JDEM the full PDF does not bring significant tightening of constraints on cosmology beyond what is already achievable by the joint use of second and third order moments.
Published B and V fluxes from nearby Type Ia supernovae are fitted to light-curve templates with 4-6 adjustable parameters. Separately, B magnitudes from the same sample are fitted to a linear dependence on B-V color within a post-maximum time window prescribed by the CMAGIC method. These fits yield two independent SN magnitude estimates B_max and B_BV. Their difference varies systematically with decline rate Delta m_15 in a form that is compatible with a bilinear but not a linear dependence; a nonlinear form likely describes the decline-rate dependence of B_max itself. A Hubble fit to the average of B_max and B_BV requires a systematic correction for observed B-V color that can be described by a linear coefficient R = 2.59 +- 0.24, well below the coefficient R_B ~ 4.1 commonly used to characterize the effects of Milky Way dust. At 99.9% confidence the data reject a simple model in which no color correction is required for SNe that are clustered at the blue end of their observed color distribution. After systematic corrections are performed, B_max and B_BV exhibit mutual rms intrinsic variation equal to 0.074 +- 0.019 mag, of which at least an equal share likely belongs to B_BV. SN magnitudes measured using maximum-luminosity or CMAGIC methods show comparable rms deviations of order ~ 0.14 mag from the Hubble line. The same fit also establishes a 95% confidence upper limit of 486 km/s on the rms peculiar velocity of nearby SNe relative to the Hubble flow.
We report spectra obtained with the Spitzer Space Telescope in the 5 to 35 micron range of HD 233517, an evolved K2 III giant with circumstellar dust. At wavelengths longer than 13 microns, the flux is a smooth continuum that varies approximately as frequency to the -5/3 power. For wavelengths shorter than 13 microns, although the star is oxygen-rich, PAH features produced by carbon-rich species at 6.3 microns, 8.2 microns, 11.3 microns and 12.7 microns are detected along with likely broad silicate emission near 20 microns. These results can be explained if there is a passive, flared disk orbiting HD 233517. Our data support the hypothesis that organic molecules in orbiting disks may be synthesized in situ as well as being incorporated from the interstellar medium.
We report the detection of broad absorptions due to Si IV 4088-4116 A in the Luminous Blue Variable (LBV) AG Carinae during its last hot phase (2001-2003). Our NLTE spectral analysis, with the radiative transfer code CMFGEN, revealed the photospheric nature of these lines predicting, however, much narrower and deeper absorption profiles than observed. Using a recently-developed code to compute synthetic spectra in 2D geometry allowing for the effects of rotation, we could match the broad absorptions with a high projected rotational velocity of 190 +/- 30 km/s on 2001 April. Analysis of spectra obtained on 2002 March and 2003 January, when the star was cooling, yielded to a projected rotational velocity of 110 +/- 10 km/s and 85 +/- 10 km/s, respectively. The derived rotational velocities are proportional to R^-1, as expected from angular momentum conservation. We discuss the effects of such high rotation on the spectral analysis of AG Car, and on the wind terminal velocity. Our results show direct spectroscopic evidence, for the first time, that a LBV may rotate at a significant fraction of its break-up velocity. Thus, AG Car (and possibly other LBVs) is indeed close to the Gamma-Omega limit, as predicted by theoretical studies of LBVs.
A great deal of study has been carried out over the last twenty years on the origin of the magnetic activity in the Galactic center. One of the most popular hypotheses assumes milli-Gauss magnetic field with poloidal geometry, pervading the inner few hundred parsecs of the Galactic-center region. However, there is a growing observational evidence for the large-scale distribution of a much weaker field of B \lesssim 10 micro G in this region. Here, we propose that the Galactic-center magnetic field originates from turbulent activity that is known to be extreme in the central hundred parsecs. In this picture the spatial distribution of the magnetic field energy is highly intermittent, and the regions of strong field have filamentary structures. We propose that the observed nonthermal radio filaments appear in (or, possibly, may be identified with) such strongly magnetized regions. At the same time, the large-scale diffuse magnetic field is weak. Both results of our model can explain the magnetic field measurements of the the Galactic-center region. In addition, we discuss the role of ionized outflow from stellar clusters in producing the long magnetized filaments perpendicular to the Galactic plane.
(abridged) High-z galaxy redshift surveys open up exciting possibilities for precision determinations of neutrino masses and inflationary models. The high-z surveys are more useful for cosmology than low-z ones owing to much weaker non-linearities in matter clustering, redshift-space distortion and galaxy bias. We can then utilize the two-dimensional information of the linear power spectrum in angular and redshift space to measure the scale-dependent suppression of matter clustering due to neutrino free-streaming as well as the shape of the primordial power spectrum. To illustrate capabilities of high-z surveys for constraining neutrino masses and the primordial power spectrum, we compare three future redshift surveys covering 300 square degrees at 0.5<z<2, 2<z<4, and 3.5<z<6.5. We find that, combined with the CMB data expected from the Planck satellite, these surveys allow precision determination of the total neutrino mass with the projected errors of sigma(m_nu) =0.059, 0.043, and 0.025 eV, respectively, which are actually smaller than the lower limits to the neutrino masses implied from the neutrino oscillation experiments, by up to a factor of 4 for the highest redshift survey. The accuracies of constraining the tilt and running index of the primordial power spectrum, sigma(n_s)=(3.8, 3.7, 3.0)x10^-3, and sigma(alpha_s)=(5.9, 5.7, 2.4)x10^-3, respectively, are smaller than the current uncertainties by more than an order of magnitude, which will allow us to discriminate between candidate inflationary models. In particular, the error on alpha_s from the highest redshift survey is not very far away from the prediction of a class of simple inflationary models driven by a massive scalar field with self-coupling, alpha_s=-(0.8-1.2)x10^-3.
We explore the growth of super-massive black holes and host galaxy bulges in the galaxy population using the Millennium Run LCDM simulation coupled with a model of galaxy formation. We find that, if galaxy mergers are the primary drivers for both bulge and black hole growth, then in the simplest picture one should expect the mBH-mBulge relation to evolve with redshift, with a larger black hole mass associated with a given bulge mass at earlier times relative to the present day. This result is independent of an evolving cold gas fraction in the galaxy population. The evolution arises from the disruption of galactic disks during mergers that make a larger fractional mass contribution to bulges at low redshift than at earlier epochs. There is no comparable growth mode for the black hole population. Thus, this effect produces evolution in the mBH-mBulge relation that is driven by bulge mass growth and not by black holes.
We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from high-resolution cosmological simulations to investigate the effects of interloper contamination and show that interlopers significantly bias the estimated slope of the projected radial distribution of satellites. We also show that the distribution of interlopers around galaxies is expected to be non-uniform in velocity space because galaxies are clustered and reside in crowded environments. Successful methods of interloper contamination correction should therefore take into account environments of the host galaxies. Two such new methods are presented and the most reliable of them is used to correct for interloper contamination in analyses of the SDSS galaxy sample. The best fit power-law slope of the interloper-corrected surface density distribution of satellites, Sigma(R) ~ R^alpha, in the volume-limited SDSS sample is alpha ~ -1.7 +/- 0.1, independent of the galaxy and satellite luminosities. Comparison with Lambda-CDM simulations shows that the radial distribution of the SDSS satellites is more concentrated than that of subhalos around galaxy-sized halos, especially at R < 100 kpc/h. The predicted dark matter radial distribution is somewhat more concentrated than the profile of the SDSS satellites, but the difference is not statistically significant for our sample.
The apparent shapes of spiral galaxies in the 2-Micron All Sky Survey Large Galaxy Atlas are used to constrain the intrinsic shape of their disks. When the distribution of apparent axis ratios is estimated using a nonparametric kernel method, the shape distribution is inconsistent with axisymmetry at the 90% confidence level in the B band and at the 99% confidence level in the K band. If spirals are subdivided by Hubble type, the late-type spirals (Sc and later) are consistent with axisymmetry, while the earlier spirals are strongly inconsistent with axisymmetry. The distribution of disk ellipticity can be fitted adequately with either a Gaussian or a lognormal distribution. The best fits for the late spirals imply a median ellipticity of epsilon = 0.07 in the B band and epsilon = 0.02 in the K band. For the earlier spirals, the best fits imply a median ellipticity of epsilon = 0.18 in the B band and epsilon = 0.30 in the K band. The observed scatter in the Tully-Fisher relation, for both late and early spirals, is consistent with the disk ellipticity measured in the B band. This indicates that excluding spirals of Hubble type earlier than Sc will reduce the intrinsic scatter in the Tully-Fisher relation used as a distance indicator.
The cosmological potential of large-scale structure observations for cosmology have been extensively discussed in the litterature. In particular, it has recently been shown how Sunyaev-Zel'dovich (SZ) cluster surveys can be used to constrain dark energy parameters. In this paper, we study whether selection and systematics effects will limit future wide-field SZ surveys from achieving their cosmological potential. For this purpose, we use a sky simulation and an SZ-cluster detection software presented in Pires et al. (2005), using the future Olimpo, APEX and Planck surveys as a concrete examples. We show that the SZ-cluster selection function and contamination of SZ-cluster catalogues are more complex than is usually assumed. In particular, the simulated field-to-field detected cluster counts is a factor 3 larger than the expected Poisson fluctuations. We also study the impact of missing redshift information and of the uncertainty of the scaling relations for low mass clusters. We quantify, through hypothesis tests, how near-future SZ experiments can be used to discriminate between different structure formation models. Using a maximum likelihood approach, we then study the impact of these systematics on the joint measurement of cosmological models and of cluster scaling relations.
We present the discovery of the first X-ray counterpart to a Rotating RAdio Transient (RRAT) source. RRAT J1819--1458 is a relatively highly magnetized (B $\sim 5\times10^{13}$ G) member of a new class of unusual pulsar-like objects discovered by their bursting activity at radio wavelengths. The position of RRAT J1819--1458 was serendipitously observed by the {\sl Chandra} ACIS-I camera in 2005 May. At that position we have discovered a pointlike source, CXOU J181934.1--145804, with a soft spectrum well fit by an absorbed blackbody with $N_H = 7^{+7}_{-4} \times 10^{21}$ cm$^{-2}$ and temperature $kT=0.12 \pm 0.04$ keV, having an unabsorbed flux of $\sim2 \times 10^{-12}$ ergs cm$^{-2}$ s$^{-1}$ between 0.5 and 8 keV. No optical or infrared (IR) counterparts are visible within $1''$ of our X-ray position. The positional coincidence, spectral properties, and lack of an optical/IR counterpart make it highly likely that CXOU J181934.1--145804 is a neutron star and is the same object as RRAT J1819--1458. The source showed no variability on any timescale from the pulse period of 4.26~s up to the five-day window covered by the observations, although our limits (especially for pulsations) are not particularly constraining. The X-ray properties of CXOU J181934.1--145804, while not yet measured to high precision, are similar to those of comparably-aged radio pulsars and are consistent with thermal emission from a cooling neutron star.
The intergalactic neutral hydrogen which is responsible for the Lyman alpha forest of quasar absorption is a tracer of much larger amounts of ionised hydrogen. The ionised component has yet to be detected directly, but is expected to scatter CMB photons via the Sunyaev-Zel'dovich (SZ) effect. We use hydrodynamic simulations of a LambdaCDM universe to create mock quasar spectra and CMB sky maps. We find that the high-z Lya forest gas causes temperature fluctuations of the order of 1 muK rms in the CMB on arcmin scales. The kinetic and thermal SZ effects have a similar magnitude at z=3, with the thermal effect becoming relatively weaker as expected at higher z. The CMB signal associated with lines of sight having HI column densities > 10^18 cm^-2 is only marginally stronger than that for lower column densities. The strong dependence of rms temperature fluctuation on mean Lya absorbed flux, however, suggests that the CMB signal effectively arises in lower density material. We investigate the use of the cross-correlation of the Lya forest and the microwave background to detect the SZ effect at redshifts 2-4. In so doing we are able to set direct limits on the density of diffuse ionised intergalactic baryons. We carry out a preliminary comparison at a mean redshift z=3 of 3488 quasar spectra from SDSS Data Release 3 and the WMAP first year data. Assuming that the baryons are clustered as in a LambdaCDM cosmology, and have the same mean temperature, the cross-correlation yields a weak limit on the cosmic density of ionised baryons Omega_(b,I), which is Omega_(b,I) < 0.8 at 95% confidence. With data from upcoming CMB telescopes, we anticipate that a direct detection of the high redshift ionised IGM will soon be possible, providing an important consistency check on cosmological models.
We investigate statistically the dynamical consequences of cosmological fluxes of matter and related moments on progenitors of today's dark matter haloes. Their dynamics is described via canonical perturbation theory which accounts for two types of perturbations: the tidal field corresponding to fly-bys and accretion of dark matter through the halo's outer boundary. he dynamical equations are solved linearly, order by order, projecting on a biorthogonal basis to consistently satisfy the field equation. Since our solution of the Boltzmann Poisson equations is explicit, it allows statistical predictions for the ensemble distribution of the inner dynamical features of haloes. The secular evolution of open galactic haloes is investigated: we derive the kinetic equation which governs the quasi-linear evolution of dark matter profile induced by infall and its corresponding gravitational correlations. This yields a Fokker Planck-like equation for the angle-averaged underlying distribution function. We show how these extensions to the classical theory could be used to (i) observationally constrain the statistical nature of the infall (ii) predict the observed distribution and correlations of substructures in upcoming surveys, (iii) predict the past evolution of the observed distribution of clumps, and finally (iv) weight the relative importance of the intrinsic (via the unperturbed distribution function) and external (tidal and/ or infall) influence of the environment in determining the fate of galaxies.
We present results from a survey of a 1300 arcmin^2 region of the Orion B South molecular cloud, including NGC 2024, NGC 2023, and the Horsehead Nebula (B33), obtained using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. Submillimeter continuum observations at 450 microns and 850 microns are discussed. Using an automated algorithm, 57 discrete emission features (``clumps'') are identified in the 850 micron map. The physical conditions within these clumps are investigated under the assumption that the objects are in quasi-hydrostatic equilibrium. The best fit dust temperature for the clumps is found to be T_d = 18 +/- 4 K, with the exception of those associated with the few known far infrared sources residing in NGC 2024. The latter internally heated sources are found to be much warmer. In the region surrounding NGC 2023, the clump dust temperatures agree with clump gas temperatures determined from molecular line excitation measurements of the CO molecule. The bounding pressure on the clumps lies in the range log(k^-1 P cm^3 K^-1) = 6.1 +/- 0.3. The cumulative mass distribution is steep at the high mass end, as is the stellar Initial Mass Function. The distribution flattens significantly at lower masses, with a turn-over around 3 -- 10 M_sun.
We present new optical observations of young massive star clusters in Arp 220, the nearest ultraluminous infrared galaxy, taken in UBVI with the Hubble Space Telescope ACS/HRC camera. We find a total of 206 probable clusters whose spatial distribution is centrally concentrated toward the nucleus of Arp 220. We use model star cluster tracks to determine ages, luminosities, and masses for 14 clusters with complete UBVI indices or previously published near-infrared data. We estimate rough masses for 24 additional clusters with I < 24 mag from BVI indices alone. The clusters with useful ages fall into two distinct groups: a ``young'' population (< 10 Myr) and an intermediate-age population (~300 Myr). There are many clusters with masses clearly above 10^6 Msun and possibly even above 10^7 Msun in the most extreme instances. These masses are high enough that the clusters being formed in the Arp 220 starburst can be considered as genuine young globular clusters. In addition, this study allows us to extend the observed correlation between global star formation rate and maximum cluster luminosity by more than an order of magnitude in star formation rate.
I summarize current knowledge of galaxy formation with emphasis on the initial conditions provided by the Lambda CDM cosmology, integral constraints from cosmological quantities, and the demographics of high-redshift protogalaxies. Tables are provided summarizing the number density, star formation rate and stellar mass per object, cosmic star formation rate and stellar mass densities, clustering length and typical dark matter halo masses for Lyman break galaxies, Lyman alpha emitting galaxies, Distant red galaxies, Sub-millimeter galaxies, and Damped Lyman alpha absorption systems. I also discuss five key unsolved problems in galaxy formation and prognosticate advances that the near future will bring.
Recent observational studies of $\omega$ Centauri by {\it Hubble Space Telescope} have discovered a double main sequence in the color magnitude diagrams (CMDs) of its stellar populations. The stellar population with the blue main sequence (bMS) is observationally suggested to have a helium abundance much larger, by $\Delta Y\sim 0.12$, than that of the red main sequence (rMS). By using somewhat idealized models in which stars of the bMS are formed from gas ejected from those of the rMS, we quantitatively investigate whether the helium overabundance of the bMS can result from self-enrichment from massive AGB stars, from mass loss of very massive young stars, or from type II supernovae within $\omega$ Cen. We show that as long as the helium enrichment is due to ejecta from the rMS formed earlier than the bMS, none of the above three enrichment scenarios can explain the observed properties of the bMS self-consistently for reasonable IMFs. The common, serious problem in all cases is that the observed number fraction of the bMS can not be explained without assuming unusually top-heavy IMFs. This failure of the self-enrichment scenarios implies that most of the helium-enriched gas necessary for the formation of the bMS originated from other external sources. We thus suggest a new scenario that most of the second generation of stars (i.e., the bMS) in $\omega$ Cen could be formed from gas ejected from field stellar populations that surrounded $\omega$ Cen when it was a nucleus of an ancient dwarf galaxy.
Recent observations have discovered star formation activities in the extreme outer regions of disk galaxies. However it remains unclear what physical mechanisms are responsible for triggering star formation in such low-density gaseous environments of galaxies. In order to understand the origin of these outer star-forming regions, we numerically investigate how the impact of dark matter subhalos orbiting a gas-rich disk galaxy embedded in a massive dark matter halo influences the dynamical evolution of outer HI gas disk of the galaxy. We find that if the masses of the subhalos ($M_{\rm sb}$) in a galaxy with an extended HI gas disk are as large as $10^{-3} \times M_{\rm h}$, where $M_{\rm h}$ is the total mass of the galaxy's dark halo, local fine structures can be formed in the extended HI disk. We also find that the gas densities of some apparently filamentary structures can exceed a threshold gas density for star formation and thus be likely to be converted into new stars in the outer part of the HI disk in some models with larger $M_{\rm sb}$. These results thus imply that the impact of dark matter subhalos (``dark impact'') can be important for better understanding the origin of recent star formation discovered in the extreme outer regions of disk galaxies. We also suggest that characteristic morphologies of local gaseous structures formed by the dark impact can indirectly prove the existence of dark matter subhalos in galaxies. We discuss the origin of giant HI holes observed in some gas-rich galaxies (e.g., NGC 6822) in the context of the dark impact.
Multi-wavelength observations of Galactic black hole transients during outburst decay are instrumental for our understanding of the accretion geometry and the formation of outflows around black hole systems. H1743-322, a black hole transient observed intensely in X-rays and also covered in the radio band during its 2003 decay, provides clues about the changes in accretion geometry during state transitions and also the general properties of X-ray emission during the intermediate and the low-hard states. In this work, we report on the evolution of spectral and temporal properties in X-rays and the flux in the radio band with the goal of understanding the nature of state transitions observed in this source. We concentrate on the transition from the thermal dominant state to the intermediate state that occurs on a timescale of one day. We show that the state transition is associated with a sudden increase in power-law flux. We determine that the ratio of the power-law flux to the overall flux in the 3--25 keV band must exceed 0.6 to observe strong timing noise. Even after the state transition, once this ratio was below 0.6, the system transited back to the thermal dominant state for a day. We show that the emission from the compact radio core does not turn on during the transition from the thermal dominant state to the intermediate state but does turn on when the source reaches the low-hard state, as seen in 4U 1543-47 and GX 339-4. We find that the photon index correlates strongly with the QPO frequency and anti-correlates with the rms amplitude of variability. We also show that the variability is more likely to be associated with the power-law emission than the disk emission.
In the course of our search for double degenerate binaries as potential progenitors of type Ia supernovae with the ESO VLT (ESO SN Ia Progenitor SurveY SPY) several new subdwarf B (sdB) binaries with sdB primary components were discovered. In this paper, we present detailed analyses of six radial velocity variable sdB stars. Radial velocity curves have been measured. From the mass functions we derive lower limits to the masses of the unseen companions and we discuss their nature. In addition, stellar parameters like effective temperatures, surface gravities and helium abundances were determined as well as metal abundances.
Modeling the multiwavelength emission of successive regions in the jet of the quasar PKS 1136-135 we find indication that the jet suffers deceleration near its end on a (deprojected) scale of ~400 kpc. Adopting a continuous flow approximation we discuss the possibility that the inferred deceleration from a Lorentz factor Gamma=6.5 to Gamma=2.5 is induced by entrainment of external gas. Some consequences of this scenario are discussed.
In this paper we report a new method for determining the thickness of the non-edge-on disk galaxies. Our method is based on the observation comparing with the theoretical researches of the distribution of the vertical velocity dispersion, which is obtained from the solutions of three dimensional Poisson's equation and the galactic dynamical equation. This method also allows us to investigate the mass-to-light ratio of the disk. As examples, the thickness and mass-to-light ratio of two disk galaxies, NGC 1566 and NGC 5247, which have been extensively studied in spectroscopy, have been calculated and the results are presented.
We study how the proportion of star-forming galaxies evolves between z=0.8 and z=0 as a function of galaxy environment, using the [OII] line in emission as a signature of ongoing star formation. Our high-z dataset comprises 16 clusters, 10 groups and another 250 galaxies in poorer groups and the field at z=0.4-0.8 from the ESO Distant Cluster Survey, plus another 9 massive clusters at similar redshifts. As a local comparison, we use samples of galaxy systems selected from the Sloan Digital Sky Survey at 0.04< z < 0.08. At high-z most systems follow a broad anticorrelation between the fraction of star-forming galaxies and the system velocity dispersion. At face value, this suggests that at z=0.4-0.8 the mass of the system largely determines the proportion of galaxies with ongoing star formation. At these redshifts the strength of star formation (as measured by the [OII] equivalent width) in star-forming galaxies is also found to vary systematically with environment. Sloan clusters have much lower fractions of star-forming galaxies than clusters at z=0.4-0.8 and, in contrast with the distant clusters, show a plateau for velocity dispersions $ \ge 550 km s^-1$, where the fraction of galaxies with [OII] emission does not vary systematically with velocity dispersion. We quantify the evolution of the proportion of star-forming galaxies as a function of the system velocity dispersion and find it is strongest in intermediate-mass systems (sigma ~ 500-600 km s^-1 at z=0). To understand the origin of the observed trends, we use the Press-Schechter formalism and the Millennium Simulation and show that galaxy star formation histories may be closely related to the growth history of clusters and groups. We propose a scheme that is able to account for the observed relations between the star-forming fraction and \sigma [abridged].
We present Swift and XMM observations of GRB 050326, detected by Swift-BAT. The fluence was 7.7x10^-6 erg cm^-2 (20-150 keV), and its spectrum was hard, with a power law photon index 1.25. The afterglow light curve did not show any break nor flares between ~1 hr and ~6 d after the burst, and decayed with a slope 1.70. The afterglow spectrum is well fitted by a power-law model, suffering absorption both in the Milky Way and in the host galaxy. The rest-frame Hydrogen column density is significant, N_H_z > 4x10^21 cm^-2, and the redshift of the absorber is z > 1.5. There was good agreement between the Swift-XRT and XMM results. By comparing the prompt and afterglow fluxes, we found that an early break occurred before the XRT observation. The properties of the GRB 050326 afterglow are well described by a spherical fireball expanding in a uniform external medium, so a further steepening is expected at later times. The lack of such a break constrains the jet angle to be >7 deg. Using the redshift constraints provided by the X-ray analysis, we also estimated that the beaming-corrected gamma-ray energy was >3x10^51 erg, at the high end of GRB energies. Despite the brightness in X rays, only deep limits could be placed by Swift-UVOT at optical/UV wavelengths. Thus, this GRB was "truly dark", with the optical-to-X-ray spectrum violating the synchrotron limit. The optical and X-ray observations are consistent either with an absorbed event or with a high-redshift one. To obey the Ghirlanda relation, a moderate/large redshift z>4.5 is required. (abridged)
The ``central engine'' of AGN is thought to be powered by accretion on a central nucleus believed to be a super-massive black hole. The localization and exact mechanism of the energy release in AGN are still not well understood. We present observational evidence for the link between variability of the radio emission of the compact jet, optical and X-ray continua emission and ejections of new jet components in the radio galaxy 3C 390.3. The time delays between the light curves of the individual jet components and the light curve of the optical continuum are estimated by using minimization methods and the discret correlation function. We find that the variations of the optical continuum are correlated with radio emission from a stationary feature in the jet. This correlation indicates that the source of variable non-thermal continuum radiation is located in the innermost part of the relativistic jet. We suggest that the continuum emission from the jet and counterjet ionizes material in a subrelativistic outflow surrounding the jet, which results in a formation of two conical regions with broad emission lines (in addition to the conventional broad line region around the central nucleus) at a distance more or equal to 0.4 parsecs from the central engine. Implications for modeling of the broad-line regions are discussed.
The fluctuating accretion model of Lyubarskii (1997) and its extension by Kotov et al. (2001), seeks to explain the spectral-timing properties of the X-ray variability of accreting black holes in terms of inward-propagating mass accretion fluctuations produced at a broad range of radii. The fluctuations modulate the X-ray emitting region as they move inwards and can produce temporal-frequency-dependent lags between energy bands, and energy-dependent power spectral densities (PSDs) as a result of the different emissivity profiles, which may be expected at different X-ray energies. Here we use a simple numerical implementation to investigate in detail the X-ray spectral-timing properties of the model and their relation to several physically interesting parameters, namely the emissivity profile in different energy bands, the geometrical thickness and viscosity parameter of the accretion flow, the strength of damping on the fluctuations and the temporal coherence (measured by the `quality-factor', Q) of the fluctuations introduced at each radius. We find that a geometrically thick flow with large viscosity parameter is favoured, and confirm that the predicted lags are quite robust to changes in the emissivity profile, and physical parameters of the accretion flow, which may help to explain the similarity of the lag spectra in the low/hard and high/soft states of Cyg X-1. We also demonstrate the model regime where the light curves in different energy bands are highly spectrally coherent. We compare model predictions directly to X-ray data from the Narrow Line Seyfert~1 galaxy NGC 4051 and the BHXRB Cyg X-1 in its high/soft state and show that this general scheme can reproduce simultaneously the time lags and energy-dependence of the PSD.
Khomenko et al. estimate the mean magnetic field strength of the quiet Sun to be 20 G. The figure is smaller than several existing estimates, and it comes from the comparison between observed Zeeman polarization signals and synthetic signals from numerical simulations of magneto-convection. The numerical simulations require an artificially large magnetic diffusivity, which smears out magnetic structures smaller than the grid scale. Assuming a turbulent cascade for the unresolved artificially smeared magnetic fields, we find that their unsigned magnetic flux is at least as important as that explicitly shown in the simulation. The unresolved fields do not produce Zeeman polarization but contribute to the unsigned flux.Since they are not considered by Khomenko et al., their mean magnetic field strength has to be regarded as a lower limit. This kind of bias is not specific of a particular numerical simulation or a spectral line. It is to be expected when observed quiet Sun Zeeman signals are compared with synthetic signals from simulations.
We report on the discovery of a z = 3.16 Lyman-alpha emitting blob in the GOODS South field. The blob has a total Ly-alpha luminosity of ~ 10^(43) erg s^(-1) and a diameter larger than 60 kpc. The available multi-wavelength data in the GOODS field consists of 13 bands from X-rays (Chandra) to infrared (Spitzer). Unlike other discovered Ly-alpha blobs, this blob shows no obvious continuum counter-part in any of the broad-bands. In particular, no optical counter-parts are found in the deep HST/ACS imaging available. For previously published blobs, AGN (Active Galactic Nuclei) or 'superwind' models have been found to provide the best match with the data. We here argue that the most probable origin of the extended Ly-alpha emission from the blob in the GOODS South field is cold accretion onto a dark matter halo.
Several items on the diagnostics and interpretation of coronal loop observations are under debate. In this work, we analyze a well-defined loop system detected in a time-resolved observation in several spectral bands. The dataset includes simultaneous images in the TRACE 171 A, 195 A and 284 A bands, and Yohkoh/SXT, and two rasters taken with SoHO/CDS in twelve relevant lines. The loop is initially best visible in the TRACE 195 A filter band, and later in the 171 A filter band, with correspondence with the CDS raster images at log T \~ 6.0-6.1. We have taken as pixel-by-pixel background the latest TRACE, Yohkoh and CDS images where the loop has faded out. We examine the loop morphology evolution, the light curves, the TRACE filter ratio distribution and evolution, the images and emission measure from the CDS spectral lines. Our analysis detects that, after background subtraction, the emission along the loop and its evolution are non-uniform, especially in the 171 A filter band, and that the TRACE 195/171 filter ratio has a moderately non-uniform distribution along the loop and evolves in time. Both the light curves and the filter ratio evolution indicate a globally cooling loop. Relatively hot plasma may be present at the beginning while, during the first CDS raster, the data indicate a rather moderate thermal structuring of the loop. Our data analysis supports a coherent scenario across the different bands and instruments, points out difficulties in diagnostic methods and puts quantitative basis for detailed forward modeling.
A new type of cosmic-ray experiment is under construction in the Pyh\"asalmi mine in the underground laboratory of the University of Oulu, Finland. It aims to study the composition of cosmic rays at and above the knee region. The experiment, called EMMA, will cover approximately 150 square-metres of detector area. The array is capable of measuring the multiplicity and the lateral distribution of underground muons, and the arrival direction of the air shower. The full-size detector is expected to run by the end of 2007.
In our continuing optical spectroscopic campaign to identify the longer-wavelength counterparts of newly-discovered hard X-ray sources detected by INTEGRAL, we observed the putative optical counterparts of seven southern sources at the South African Astronomical Observatory and at the European Southern Observatory. For two of these objects, optical photometry was also acquired. These observations firmly established the nature of four of them: we found that IGR J10404-4625 (=LEDA 93974), 4U 1344-60 and IGR J16482-3036 are Active Galactic Nuclei (AGNs) at redshifts z = 0.0237, 0.013 and 0.0313, respectively, and that 2RXP J130159.6-635806 is a Galactic High-Mass X-ray Binary (HMXB). We also give possible optical identifications for three further objects, namely IGR J11215-5952, IGR J11305-6256 and IGR J16207-5129, which are consistent with being Galactic HMXBs. Physical parameters for these objects are also evaluated by collecting and discussing the available multiwavelength information. The detection of four definite or likely HMXBs out of seven objects in our sample further stresses INTEGRAL's crucial contribution in hunting this class of object. Also, the determination of the extragalactic nature of a substantial fraction of the INTEGRAL survey sources underlines the importance of hard X-ray observations for the study of background AGNs located beyond the `Zone of Avoidance' of the Galactic Plane.
A prolonged timing of millisecond pulsars has revealed low-frequency
uncorrelated noise, presumably of astrophysical origin, in the pulse arrival
time (PAT) residuals for some of them. In most cases, pulsars in globular
clusters show a low-frequency modulation of their rotational phase and spin
rate. The relativistic time delay of the pulsar signal in the curved space time
of randomly distributed and moving globular cluster stars (the Shapiro effect)
is suggested as a possible cause of this modulation.
Given the smallness of the aberration corrections that arise from the
nonstationarity of the gravitational field of the randomly distributed ensemble
of stars under consideration, a formula is derived for the Shapiro effect for a
pulsar in a globular cluster. The derived formula is used to calculate the
autocorrelation function of the low-frequency pulsar noise, the slope of its
power spectrum, and the behavior of the $\sigma_z$ statistic that characterizes
the spectral properties of this noise in the form of a time function. The
Shapiro effect under discussion is shown to manifest itself for large impact
parameters as a low-frequency noise of the pulsar spin rate with a spectral
index of n=-1.8 that depends weakly on the specific model distribution of stars
in the globular cluster. For small impact parameters, the spectral index of the
noise is n=-1.5.
We present here the results of a deep (130 ks) XMM-Newton observation of the cluster of galaxies 2A 0335+096. The deep exposure allows us to study in detail its temperature structure and its elemental abundances. We fit three different thermal models and find that the multi-temperature wdem model fits our data best. We find that the abundance structure of the cluster is consistent with a scenario where the relative number of Type Ia supernovae contributing to the enrichment of the intra-cluster medium is ~25%, while the relative number of core collapse supernovae is ~75%. Comparison of the observed abundances to the supernova yields does not allow us to put any constrains on the contribution of Pop III stars to the enrichment of the ICM. Radial abundance profiles show a strong central peak of both Type Ia and core collapse supernova products. Both the temperature and iron abundance maps show an asymmetry in the direction of the elongated morphology of the surface brightness. In particular the temperature map shows a sharp change over a brightness edge on the southern side of the core, which was identified as a cold front in the Chandra data. This suggests that the cluster is in the process of a merger with a subcluster. Moreover, we find that the blobs or filaments discovered in the core of the cluster by Chandra are, contrary to the previous results, colder than the ambient gas and they appear to be in pressure equilibrium with their environment.
We study axisymmetric models of layered protoplanetary discs taking radiative transfer effects into account, and allowing for a residual viscosity in the dead zone. We also explore the effect of different viscosity prescriptions. In addition to the ring instability reported in the first paper of the series we find an oscillatory instability of the dead zone, accompanied by variations of the accretion rate onto the central star. We provide a simplified analytical description explaining the mechanism of the oscillations. Finally, we find that the residual viscosity enables stationary accretion in large regions of layered discs. Based on results obtained with the help of a simple 1-D hydrocode we identify these regions, and discuss conditions in which layered discs can give rise to FU~Orionis phenomena.
Kashlinsky et al. (2005) find a significant cosmic infrared background fluctuation excess on angular scales >50 arcsec that cannot be explained by instrumental noise or local foregrounds. The excess has been tentatively attributed to emission from primordial very massive (PopIII) stars formed <200 Myr after the Big Bang. Using an evolutionary model motivated by independent observations and including various feedback processes, we find that PopIII stars can contribute <40% of the total background intensity, \nu J_\nu ~ 1-2 nW m^-2 sr^-1 in the 0.8-8 \mum range, produced by all galaxies at z>5. The clustering of these high redshift sources accounts very precisely for the infrared fluctuation excess. We conclude that the origin of such fluctuations can be attributed to galaxies at redshift z>5 predominantly hosting stars with masses and properties similar to the present ones.
Archival observations of 18 starburst galaxies that span a wide range in
metallicity reveal for the first time a correlation between the ratio of
emission line fluxes of [FeII] at 26 microns and [NeII] at 12.8 microns and the
7.7 micron PAH strength, with the [FeII]/[NeII] flux ratio decreasing with
increasing PAH strength. We also find a strong correlation between the
[FeII]/[NeII] flux ratio and the host galaxy metallicity, with the flux ratio
decreasing with increasing metallicity.
Since [FeII] emission has been linked primarily to supernova shocks, we
attribute the high [FeII]/[NeII] ratios in low-metallicity galaxies to enhanced
supernova activity. We consider this to be a dominant mechanism for PAH
destruction, rather than grain destruction in photoionized regions surrounding
young massive stars. We also consider whether the extreme youth of the
low-metallicity galaxies is responsible for the lack of PAH emission.
We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe in a region of radius 26Mpc/h, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter haloes identified in a high resolution N-body simulation of the standard LCDM cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of dark matter particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover present-day halo distances with typical errors of less than 3 per cent, compared to 5 per cent errors assuming Hubble flow distances. The total halo mass and the linear bias were both found to be constained at the 50 per cent level. The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions, and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by SIM and GAIA.
We present a model independent method of reconstructing the Lagrangian for the k-essence field driving the present acceleration of the universe. We consider the simplest k-essence model for which the potential is constant. Later we use three parametrizations for the Hubble parameter $H(z)$, consistent with the recent SN1a data, to yield the Lagrangian $F$. Our reconstruction program does not generate any physically realistic Lagrangian for models that allow phantom crossing, whereas models without phantom crossing, yield well behaved Lagrangian.
We present the results of time-resolved spectroscopy of 13 O-type stars in the Cas OB6 stellar association. We conducted a survey for radial velocity variability in search of binary systems, which are expected to be plentiful in young OB associations. Here we report the discovery of two new single-lined binaries, and we present new orbital elements for three double-lined binaries (including one in the multiple star system HD 17505). One of the double-lined systems is the eclipsing binary system DN Cas, and we present a preliminary light curve analysis that yields the system inclination, masses, and radii. We compare the spectra of the single stars and the individual components of the binary stars with model synthetic spectra to estimate the stellar effective temperatures, gravities, and projected rotational velocities. We also make fits of the spectral energy distributions to derive E(B-V), R=A_V/E(B-V), and angular diameter. A distance of 1.9 kpc yields radii that are consistent with evolutionary models. We find that 7 of 14 systems with spectroscopic data are probable binaries, consistent with the high binary frequency found for other massive stars in clusters and associations.
Interpretation of the energy spectrum and arrival distribution of cosmic rays is complicated by lack of knowledge of the nature of the primaries. We review claims for the mass composition above 1017 eV where it can be determined only indirectly from air-shower observables. Difficulties in comparisons between data arise because, inevitably, a set of measurements is interpreted using the best model of hadronic interactions available at the time of analysis. We discuss the situation and conclude that the evidence for a proton-dominated mass composition, even at the highest energies, is unconvincing. However, it may be that there are consistent differences between mass measurements from optical techniques and those based upon other shower observables. We also find that iron nuclei of ultra high energy can probably escape from the galaxies that host GRBs, possible cosmic ray accelerators. The accelerators must lie nearby.
Aims: To investigate the stability and angular momentum transport by the strato-rotational instability in the nonlinear regime. Methods: The hydrodynamic compressible equations are solved in a cartesian box in which the outer cylinder is embedded. Gravity along the rotation axis leads to density stratification. No-slip boundary conditions are used in the radial direction, while free-slip conditions are used on the two ends of the cylinders. Results: The strato-rotational instability is confirmed and the Reynolds stress is shown to transport angular momentum away from the axis. However, the growth rate decreases with increasing Reynolds number. This, as well as the presence of boundaries renders this instability less relevant for astrophysical applications.
We provide improved atomic calculation of wavelengths, oscillator strengths, and autoionization rates relevant to the $2\to 3$ inner-shell transitions of Fe~VI--XVI, the so-called Fe~M-shell unresolved transition array (UTA). A second order many-body perturbation theory is employed to obtain accurate transition wavelengths, which are systematically larger than previous theoretical results by 15--45~m{\AA}. For a few transitions of Fe~XVI and Fe~XV where laboratory measurements exist, our new wavelengths are accurate to within a few m{\AA}. Using these new calculations, the apparent discrepancy in the velocities between the Fe~M-shell UTA and other highly ionized absorption lines in the outflow of NGC 3783 disappears. The oscillator strengths in our new calculation agree well with the previous theoretical data, while the new autoionization rates are significantly larger, especially for lower charge states. We attribute this discrepancy to the missing autoionization channels in the previous calculation. The increased autoionization rates may slightly affect the column density analysis of the Fe~M-shell UTA for sources with high column density and very low turbulent broadening. The complete set of atomic data is provided as an electronic table.
The intense 0.511 MeV gamma-ray line emission from the Galactic Center observed by INTEGRAL requires a large annihilation rate of nonrelativistic positrons. If these positrons are injected at even mildly relativistic energies, higher-energy gamma rays will also be produced. We calculate the gamma-ray spectrum due to inflight annihilation and compare to the observed diffuse Galactic gamma-ray data. Even in a simplified but conservative treatment, we find that the positron injection energies must be $\lesssim 3$ MeV, which strongly constrains models for Galactic positron production.