Showing papers in "arXiv: Astrophysics in 2009"

RXTE Monitoring of the Anomalous X-ray Pulsar 1E 1048.1-5937: Long-Term Variability and the 2007 March Event

Abstract: After three years of no unusual activity, Anomalous X-ray Pulsar 1E 1048.1-5937 reactivated in 2007 March. We report on the detection of a large glitch (deltav/v = 1.63(2) x 10(exp -5)) on 2007 March 26 (MJD 54185.9), contemporaneous with the onset of a pulsed-flux flare, the third flare observed from this source in 10 years of monitoring with the Rossi X-ray Timing Explorer. Additionally, we report on a detailed study of the evolution of the timing properties, the pulsed flux, and the pulse profile of this source as measured by RXTE from 1996 July to 2008 January. In our timing study, we attempted phase coherent timing of all available observations. We show that in 2001, a timing anomaly of uncertain nature occurred near the rise of the first pulsed flux flare; we show that a likely glitch (deltav/v = 2.91(9) x 10(exp -6)) occurred in 2002, near the rise of the second flare, and we present a detailed description of the variations in the spin-down. In our pulsed flux study, we compare the decays of the three flares and discuss changes in the hardness ratio. In our pulse profile study, we show that the profile exhibited large variations near the peak of the first two flares, and several small short-term profile variations during the most recent flare. Finally, we report on the discovery of a small burst 27 days after the peak of the last flare, the fourth burst discovered from this source. We discuss the relationships between the observed properties in the framework of the magnetar model.

Lyman alpha line formation in starbursting galaxies I. Moderately thick, dustless, and static HI media

Abstract: The Lya line transfer in an extremely thick medium of neutral hydrogen is investigated by adopting an accelerating scheme in our Monte Carlo code to skip a large number of core or resonant scatterings. This scheme reduces computing time significantly with no sacrifice in the accuracy of the results. We applied this numerical method to the Lya transfer in a static, uniform, dustless, and plane-parallel medium. Two types of photon sources have been considered, the midplane source and the uniformly distributed sources. The emergent profiles show double peaks and absorption trough at the line-center. We compared our results with the analytic solutions derived by previous researchers, and confirmed that both solutions are in good agreement with each other. We investigated the directionality of the emergent Lya photons and found that limb brightening is observed in slightly thick media while limb darkening appears in extremely thick media. The behavior of the directionality is noted to follow that of the Thomson scattered radiation in electron clouds, because both Lya wing scattering and Thomson scattering share the same Rayleigh scattering phase function. The mean number of wing scatterings just before escape is in exact agreement with the prediction of the diffusion approximation. The Lya photons constituting the inner part of the emergent profiles follow the relationship derived from the diffusion approximation. We present a brief discussion on the application of our results to the formation of Lya broad absorption troughs and P-Cygni type Lya profiles seen in the UV spectra of starburst galaxies.

Are ring galaxies the ancestors of giant low surface brightness galaxies

Abstract: Giant low surface brightness galaxies (GLSBs), such as Malin 1, have unusually large and flat discs. Their formation is a puzzle for cosmological simulations in the cold dark matter scenario. We suggest that GLSBs might be the final product of the dynamical evolution of collisional ring galaxies. In fact, our simulations show that, approximately 0.5-1.5 Gyr after the collision which lead to the formation of a ring galaxy, the ring keeps expanding and fades, while the disc becomes very large (~100 kpc) and flat. At this stage, our simulated galaxies match many properties of GLSBs (surface brightness profile, morphology, HI spectrum and rotation curve).

A spatially resolved study of photoelectric heating and [C II] cooling in the LMC

Abstract: Context. Photoelectric heating is a dominant heating mechanism for many phases of the interstellar medium. We study this mechanism throughout the Large Magellanic Cloud (LMC). Aims. We aim to quantify the importance of the [C II] cooling line and the photoelectric heating process of various environments in the LMC and to investigate which parameters control the extent of photoelectric heating. Methods. We use the BICE [C II] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of photoelectric heating: photoelectric heating rate over power absorbed by grains, i.e. the observed [C II] line strength over the integrated infrared emission. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study the variations as a function of H emission, dust temperatures, and the total infrared luminosity. The observed variations are interpreted in a theoretical framework. From this we estimate radiation field, gas temperature, and electron density. Results. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while the lowest coincide with bright star-forming regions (~1.4 times lower). The [C II] line emission constitutes 1.32% of the far infrared luminosity across the whole of the LMC. We find correlations between the [C II] emission and ratios of the mid infrared and far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total infrared surface brightness, S_(TIR), the [C II] surface brightness can be described as: S_([C II]) = 1.25 S^(0.69)_(TIR)[10^(-3) erg s^(-1) cm^(-2) sr^(-1)], for S_(TIR) ≳ 3.2 x 10^(-4) erg s^(-1) cm^(-2) sr^(-1). We provide a simple model of the photoelectric efficiency as a function of the total infrared luminosity. We find a power-law relation between radiation field and electron density, consistent with other studies. The [C II] emission is well-correlated with the 8 µm emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric heating process.

A Dynamical Analysis of the Poor Galaxy Clusters A2626 and A2440

Abstract: We use 189 new radial velocities, an $R$ band CCD mosiac, archival X-ray images, and gas temperatures to study substructure, masses, luminosity functions, and mass-to-light ratios in A2626 and A2440. The symmetric X-ray emission from A2626 suggests a relaxed, single-component system; however, 159 redshifts reveals a complex, three component cluster at cz~17,500km/s. One components is a typical X-ray bright cluster, a second has a cluster-like galaxy population with lower central galaxy and gas density, and the third is a background structure with field galaxy composition. A comparison of the magnitude distributions within the two subclusters suggests that A2626 is a merger in progress. Virial masses and a composite luminosity function for a region with projected radius r=1.5h^{-1}Mpc yield a mass to light ratio M/L_R~610h. Analysis of the X-ray emission from the primary component yields a gas mass fraction of 2.2h^{-3/2}% and a baryon fraction of 3.4%. A radial infall model implies that the virial mass may be an underestimate. In A2440, deep CCD $R$ band photometry indicates a striking correspondence between the galaxy and gas distributions. The galaxy distribution has three main components, each associated with a giant elliptical galaxy. The two larger peaks in the galaxy distribution coincide with the primary peaks in the X-ray emission, and the third is associated with a significant X-ray surface brightness enhancement. We use this galaxy-gas correspondence and 48 redshifts to argue that the subclusters are bound within a single system, and that the two primary components are beginning to merge. The composite luminosity function and estimates of the subcluster virial masses indicate a mass-to-light ratio from 660h to 880h.

The ART of Cosmological Simulations

Abstract: We describe the basic ideas of MPI parallelization of the N-body Adaptive Refinement Tree (ART) code. The code uses self-adaptive domain decomposition where boundaries of the domains (parallelepipeds) constantly move—with many degrees of freedom—in the search of the minimum of CPU time. The actual CPU time spent by each MPI task on previous time-step is used to adjust boundaries for the next time-step. For a typical decomposition of 53 domains, the number of possible changes in boundaries is 384≈1040. We describe two algorithms of finding minimum of CPU time for configurations with a large number of domains. Each MPI task in our code solves the N-body problem where the large-scale distribution of matter outside of the boundaries of a domain is represented by relatively few temporary large particles created by other domains. At the beginning of a zero-level time-step, domains create and exchange large particles. Then each domain advances all its particles for many small time-steps. At the end of the large step, the domains decide where to place new boundaries and re-distribute particles. The scheme requires little communications between processors and is very efficient for large cosmological simulations.

Cosmological Blastwaves and the Intergalactic Medium

Abstract: Winds from protogalactic starbursts and quasars can drive shocks that heat, ionize, and enrich the intergalactic medium The Sedov-Taylor solution for point-like explosions adequately describes these blastwaves early in their development, but as the time since the explosion ($t - t_1$) approaches the age of the universe ($t$), cosmological effects begin to alter the blastwave's structure and growth rate This paper presents an analytical solution for adiabatic blastwaves in an expanding universe, valid when the IGM is homogeneous and contains only a small fraction of the total mass density ($\Omega_{\rm IGM} << \Omega_0$) Using this analytical solution, we examine the role protogalactic explosions might play in determining the state of intergalactic gas at $z \sim 2 - 4$

Period-Luminosity-Colour distribution and classification of Galactic oxygen-rich LPVs II

Abstract: The kinematic and Period-Luminosity-Colour distribution of O-rich Long-Period Variable stars of the solar neighbourhood is interpreted in terms of pulsation modes, masses and metallicities. It is first shown that, because of input physics imperfections, the periods and mean colours derived from the existing linear and nonlinear nonadiabatic models must significantly depart from the actual behaviour of the stars. As a consequence systematic corrections have to be applied, as a first approximation, to our linear model grid. These free parameters, as well as the mixing length, are calibrated on the LPVs of the LMC and of some globular clusters, assuming a mean mass of one solar mass for the LMC Mira-like stars. Then, the masses and metallicities corresponding to the four kinematic/photometric populations of local LPVs are evaluated. The possibility of a varying mixing-length parameter is discussed and taken into account.

Morphokinematic Properties of the 21 Micron Source IRAS 07134+1005

Abstract: We report the results of a Submillimeter Array (SMA) interferometric observation of 21-micron source IRAS 07134+1005 in the CO J=3-2 line. In order to determine the morpho-kinematic properties of the molecular envelope of the object, we constructed a model using the Shape software to model the observed CO map. We find that the molecular gas component of the envelopes can be interpreted as a geometrically thick expanding torus with an expanding velocity of 8 km/s. The inner and outer radii of the torus determined by fitting Shape models are 1.2" and 3.0", respectively. The inner radius is consistent with the previous values determined by radiative transfer modeling of the spectral energy distribution and mid-infrared imaging of the dust component. The radii and expansion velocity of the torus suggest that the central star has left the asymptotic giant branch about 1140-1710 years ago, and that the duration of the equatorial enhanced mass loss is about 2560-3130 years. From the absence of an observed jet, we suggest that the formation of a bipolar outflow may lack behind in time from the creation of the equatorial torus.

New CMB Power Spectrum Constraints from MSAMI

Abstract: We present new cosmic microwave background (CMB) anisotropy results from the combined analysis of the three flights of the first Medium Scale Anisotropy Measurement (MSAM1). This balloon-borne bolometric instrument measured about 10 square degrees of sky at half-degree resolution in 4 frequency bands from 5.2 icm to 20 icm with a high signal-to-noise ratio. Here we present an overview of our analysis methods, compare the results from the three flights, derive new constraints on the CMB power spectrum from the combined data and reduce the data to total-power Wiener-filtered maps of the CMB. A key feature of this new analysis is a determination of the amplitude of CMB fluctuations at $\ell \sim 400$. The analysis technique is described in a companion paper by Knox.

Analysis of extremely low signal-to-noise ratio data from INTEGRAL/PICsIT

Abstract: Context. The PICsIT detector onboard the INTEGRAL satellite was designed to provide information about emission in the soft γ-ray band for many bright sources. Due to strong and variable instrumental background, only 4 objects have been detected so far using standard software. Aims. The moderate sensitivity of PICsIT can be compensated for in the case of many objects by adopting a long exposure time, thanks to INTEGRAL’s large field of view and an observing strategy focused on the Galactic plane. With angular resolution far higher than that of all other instruments operating in a similar energy band, PICsIT is suitable for fields too crowded or too significantly affected by Galactic diffuse emission. Therefore, it is desirable to improve the spectral extraction software to both obtain more reliable results and enlarge the number of objects that can be studied. Methods. The new PICsIT spectral extraction method is based on three elements: careful modelling of the background, an energydependent pixel-illumination function, and the computation of the probability density of the source count rate. Background maps of the detector plane are prepared for short periods and relatively narrow energy bands to insure that the background dependence on time and energy is modelled well. The most important element of the new spectral extraction method is the proper treatment of the Poisson-distributed data, developed within a Bayesian framework. Results. The new method was tested extensively on both a large true data set and simulated data. Results assumed in simulations were reproduced perfectly, without any bias and with high precision. Count rates measured for Crab were far more stable than those obtained with the standard software. For weaker sources, the new method produced spectra of far higher quality and allows us to detect at least 8 additional objects. Comparison with other INTEGRAL instruments demonstrated that PICsIT is well calibrated and provides valuable information about the continuum emission in the 250 keV – 1 MeV band, detectable currently only by INTEGRAL.

Section on prospects for dark matter detection of the white paper on the status and future of ground-based TeV gamma-ray astronomy.

Abstract: This is a report on the findings of the dark matter science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses the prospects for dark matter detection with future gamma-ray experiments, and the complementarity of gamma-ray measurements with other indirect, direct or accelerator-based searches. We conclude that any comprehensive search for dark matter should include gamma-ray observations, both to identify the dark matter particle (through the characteristics of the gamma-ray spectrum) and to measure the distribution of dark matter in galactic halos.

Kinematic sub-populations in dwarf spheroidal galaxies

Abstract: We present new spectroscopic data for twenty six stars in the recently-discovered Canes Venatici I (CVnI) dwarf spheroidal galaxy. We use these data to investigate the recent claim of the presence of two dynamically inconsistent stellar populations in this system (Ibata et al. 2006). We do not find evidence for kinematically distinct pop ulations in our sample and we are able to obtain a mass estimate for CVnI that is consistent with all available data, including previously published data. We discuss possible differences between our sample and the earlier data set and study the general detectability of sub-populat ions in small kinematic samples. We conclude that in the absence of supporting observational evidence (for example, metallicity gradients), sub-populations in small kinematic samples (typically fewer than 100 stars) should be treated with extreme caution, as their detection depends on multiple parameters and rarely produces a signal at the 3σ confidence level. It is therefore essential to determine exp licitly the statistical significance of any suggested sub-populati on.

Hydro-Gravitational Dynamics of Planets and Dark Energy

Abstract: Self gravitational fluid mechanical methods termed hydro-gravitational-dynamics (HGD) predict plasma fragmentation 0.03 Myr after the turbulent big bang to form protosuperclustervoids, turbulent protosuperclusters, and protogalaxies at the 0.3 Myr transition from plasma to gas. Linear protogalaxyclusters fragment at 0.003 Mpc viscous-inertial scales along turbulent vortex lines or in spirals, as observed. The plasma protogalaxies fragment on transition into white-hot planet-mass gas clouds (PFPs) in million-solar-mass clumps (PGCs) that become globular-star-clusters (GCs) from tidal forces or dark matter (PGCs) by freezing and diffusion into 0.3 Mpc halos with 97% of the galaxy mass. The weakly collisional non-baryonic dark matter diffuses to > Mpc scales and fragments to form galaxy cluster halos. Stars and larger planets form by binary mergers of the trillion PFPs per PGC, mostly on 0.03 Mpc galaxy accretion disks. Stars deaths depend on rates of planet accretion and internal star mixing. Moderate accretion rates pro-duce white dwarfs that evaporate surrounding gas planets by spin-radiation to form planetary nebulae before Supernova Ia events, dimming some events to give systematic distance errors, the dark energy hypothesis, and overestimates of the universe age.

Cosmology from Antarctica

Abstract: Observation of the CMB is central to observational cosmology, and the Antarctic Plateau is an exceptionally good site for this work. The first attempt at CMB observations from the Plateau was an expedition to the South Pole in December 1986 by the Radio Physics Research group at Bell Laboratories. Sky noise and opacity were measured. The results were sufficiently encouraging that in the Austral summer of 1988-1989, three CMB groups participated in the "Cucumber" campaign, where a temporary site dedicated to CMB anisotropy measurements was set up 2 km from South Pole Station. Winter-time observations became possible with the establishment in 1990 of the Center for Astrophysical Research in Antarctica (CARA), a National Science Foundation Science and Technology Center. CARA developed year-round observing facilities in the "Dark Sector", a section of Amundsen-Scott South Pole Station dedicated to astronomical observations. CARA scientists fielded several astronomical instruments: AST/RO, SPIREX, White Dish, Python, Viper, ACBAR, and DASI. By 2001, data from CARA, together with BOOMERANG, a CMB experiment on a long-duration balloon launched from McMurdo Station on the coast of Antarctica, showed clear evidence that the overall geometry of the Universe is flat, as opposed to being positively or negatively curved. In 2002, the DASI group reported the detection of polarization in the CMB. These observations strongly support the concordance model of cosmology, where the dynamics of a flat Universe are dominated by forces exerted by the Dark Energy and Dark Matter. The South Pole Telescope (SPT) is a newly-operational 10 m diameter offset telescope designed to rapidly measure anisotropies on scales much smaller than 1 degree.

Gravitational radiation driven supermassive black hole binary inspirals as periodically variable electromagnetic sources

Abstract: Supermassive black hole binaries (SMBHBs) produced in galaxy mergers are thought to complete their coalescence, below separations of r_GW=10^{-3} (M_BH/10^8 M_sun)^{3/4} pc, as their orbit decays due to the emission of gravitational waves (GWs). It may be possible to identify such GW-driven inspirals statistically in an electromagnetic (EM) survey for variable sources. A GW-driven binary spends a characteristic time T_GW at each orbital separation r_orb < r_GW that scales with the corresponding orbital time t_orb as T_GW = (const) t_orb^{8/3}. If the coalescing binary produces variations in the EM emission on this timescale, then it could be identified as a variable source with a characteristic period t_var = t_orb. The incidence rate of sources with similar inferred BH masses, showing near-periodic variability on the time-scale t_var, would then be proportional to t_var^{8/3}. Luminosity variations corresponding to a fraction f_Edd<0.01 of the Eddington luminosity would have been missed in current surveys. However, if the binary inspirals are associated with quasars, we show that a dedicated survey could detect the population of SMBHBs with a range of periods around tens of weeks. The discovery of a population of periodic sources whose abundance obeys N_var = (const) t_var^{8/3} would confirm (i) that the orbital decay is indeed driven by GWs, and (ii) that circumbinary gas is present at small orbital radii and is being perturbed by the BHs. Deviations from the t_var^{8/3} power-law could constrain the structure of the circumbinary gas disk and viscosity-driven orbital decay. We discuss constraints from existing data, and quantify the sensitivity and sky coverage that could yield a detection in future surveys.

Age-Related Observations of Low Mass Pre-Main and Young Main Sequence Stars (Invited Review)

Abstract: This overview summarizes the age dating methods available for young sub-solar mass stars. Pre-main sequence age diagnostics include the Hertzsprung-Russell (HR) diagram, spectroscopic surface gravity indicators, and lithium depletion; asteroseismology is also showing recent promise. Near and beyond the zero-age main sequence, rotation period or vsini and activ- ity (coronal and chromospheric) diagnostics along with lithium depletion serve as age proxies. Other authors in this volume present more detail in each of the aforementioned areas. Herein, I focus on pre-main sequence HR diagrams and address the questions: Do empirical young clus- ter isochrones match theoretical isochrones? Do isochrones predict stellar ages consistent with those derived via other independent techniques? Do the observed apparent luminosity spreads at constant effective temperature correspond to true age spreads? While definitive answers to these questions are not provided, some methods of progression are outlined.

Symbiotic Self-Organized Criticality and Intermittent Turbulence in an MHD Current Sheet with a Threshold Instability

Abstract: Department of Physics and Astronomy, University of Calgary, Calgary, Alberta, Canada T2N 1N4(Dated: January 8, 2009)We report numerical evidence of a self-organized criticality (SOC) and intermittent turbulence(IT) symbiosis in a resistive magnetohydrodynamic (MHD) current sheet model that includes alocal hysteretic switch to capture plasma physical processes outside of MHD that are describedin the model as current-dependent resistivity. Results from numerical simulations show scale-freeavalanches of magnetic energy dissipation characteristic of SOC, as well as multi-scaling in thevelocity field numerically indistinguishable from certain hierarchical turbulence theories. We suggestthat SOC and IT may be complementary descriptions of dynamical states realized by driven currentsheets – which occur ubiquitously in astrophysical and space plasmas.

On a systematic bias in SBF-based distances due to gravitational microlensing

Abstract: The effect of gravitational microlensing on the determination of extragalactic distances using surface brightness fluctuations (SBF) technique is considered and a method to calculate SBF amplitudes in presence of microlensing is presented. With a simple approximation for the magnification power spectrum at low optical depth the correction to the SBF-based luminosity distance is calculated. The results suggest the effect can be safely neglected at present but may become important for SBF-based Hubble diagrams at luminosity distances of about 1Gpc and beyond.

The NuMoon experiment: first results

Abstract: The NuMoon project uses the Westerbork Synthesis Radio Telescope to search for short radio pulses from the Moon. These pulses are created when an ultra high energy cosmic ray or neutrino initiates a particle cascade inside the Moon's regolith. The cascade has a negative charge excess and moves faster than the local speed of light, which causes coherent Cherenkov radiation to be emitted. With 100 hours of data, a limit on the neutrino flux can be set that is an order of magnitude better than the current one (based on FORTE). We present an analysis of the first 10 hours of data.

A Search for Deep Impact’s Large Particle Ejecta

Abstract: The Deep Impact encounter with the nucleus of 9P/Tempel ejected small grains (a ~ 100 micron) ejected by the impact will likely have the lowest ejection velocities and will weakly interact with solar radiation pressure. Therefore, large particles may remain near the nucleus for weeks or months after ejection by Deep Impact. We present initial highlights of our Spitzer Space Telescope/MIPS 24 micron camera program to image comet 9P/Tempel at 30, 80, 420, and 560 days after the Deep Impact encounter. The MIPS data, combined with our dynamics model, enable detection of large dust grains potentially ejected by Deep Impact.

Empirical Earth Rotation Model: a Consistent Way to Evaluate Earth Orientation Parameters

Abstract: It is customary to perform analysis of the Earth’s rotation in two steps: first, to present results of estimation of the Earth orientation parameters in the form of time series based on a simplified model of variations of the Earth’s rotation for a short period of time, and then to process this time series of adjustmentsby applying smoothing, re-sampling and other numerical algorithms. This approach suffers from self-inconsistency: the total set of Earth orientation parameters depends on a subjective choice of the apriori Earth orientation model, crosscorrelations between points of the time series are lost, and results of an operational analysis per se have a limited use for end users. Here, an alternative approach of direct estimation of the coefficients of expansion of Euler angle perturbations into basis functions is developed. These coefficients describe the Earth’s rotation over the entire period of observations and are evaluated simultaneously with stationpositions, source coordinates and other parameters in a single LSQ solution. In the framework of this approach considerably larger errors in apriori EOP model are tolerated. This approachgives a significant conceptual simplification of representation and estimation of the Earth’s rotation.