The WorldWide Telescope computer program, released to researchers and the public as a free resource in 2008 by Microsoft Research, has changed the way the ever-growing Universe of online astronomical data is viewed and understood. The WWT program can be thought of as a scriptable, interactive, richly visual browser of the multi-wavelength Sky as we see it from Earth, and of the Universe as we would travel within it. In its web API format, WWT is being used as a service to display professional research data. In its desktop format, WWT works in concert (thanks to SAMP and other IVOA standards) with more traditional research applications such as ds9, Aladin and TOPCAT. The WWT Ambassadors Program (founded in 2009) recruits and trains astrophysically-literate volunteers (including retirees) who use WWT as a teaching tool in online, classroom, and informal educational settings. Early quantitative studies of WWTA indicate that student experiences with WWT enhance science learning dramatically. Thanks to the wealth of data it can access, and the growing number of services to which it connects, WWT is now a key linking technology in the Seamless Astronomy environment we seek to offer researchers, teachers, and students alike.
Photoionization modeling of certain low-ionization broad absorption lines in quasars implies very compact (Delta R~0.01 pc), galaxy-scale (R kpc) absorbers blueshifted by several 1000 km s^-1. While these are likely signatures of quasar outflows, the lifetimes of such compact absorbers are too short for them to be direct ejecta from a nuclear wind. Instead, I argue that the absorbing clouds must be transient and created in situ. Following arguments detailed by Faucher-Giguere, Quataert, & Murray (2011), I show that a model in which the cool absorbers form in radiative shocks arising when a quasar blast wave impacts an interstellar cloud along the line of sight successfully explains the key observed properties. Using this radiative shock model, the outflow kinetic luminosities for three luminous quasars are estimated to be Edot,k~2-5% L_AGN (with corresponding momentum fluxes Pdot~2-15 L_AGN/c), consistent with feedback models of the M-sigma relation. These energetics are similar to those recently inferred of molecular outflows in local ultra-luminous infrared galaxies and in post-starburt winds, suggesting that active galactic nuclei (AGN) are capable of driving such outflows. Radiative shocks probably affect the multiphase structure of outflows in a range of other systems, including narrower and higher-ionization quasar absorption lines, and compact intergalactic absorbers ejected by star formation and AGN activity.
(Abridged) Compact groups, with their high number densities, small velocity dispersions, and an interstellar medium that has not been fully processed, provide a local analog to conditions of galaxy interactions in the earlier universe. The frequent and prolonged gravitational encounters that occur in compact groups affect the evolution of the constituent galaxies in a myriad of ways, for example gas processing and star formation. Recently, a statistically significant "gap" has been discovered mid-infrared IRAC colorspace of compact group galaxies. This gap is not seen in field samples and is a new example of how the compact group environment may affect the evolution of member galaxies. In order to investigate the origin and nature of this gap, we have compiled a sample of 49 compact groups. We find that a statistically significant deficit of galaxies in this gap region of IRAC colorspace is persistant in this sample, lending support to the hypothesis that the compact group environment inhibits moderate SSFRs. We note a curvature in the colorspace distribution, which is fully consistent with increasing dust temperature as the activity in a galaxy increases. This full sample of 49 compact groups allows us to subdivide the data according to physical properties of the groups. An analysis of these subsamples indicates that neither projected physical diameter nor density show a trend in colorspace within the values represented by this sample. We hypothesize that the apparent lack of a trend is due to the relatively small range of properties in this sample. Thus, the relative influence of stochastic effects becomes dominant. We analyze spectral energy distributions of member galaxies as a function of their location in colorspace and find that galaxies in different regions of MIR colorspace contain dust with varying temperatures and/or PAH emission.
We present a sample of 2865 emission line galaxies with strong nebular He II {\lambda}4686 emissions in Sloan Digital Sky Survey Data Release 7 and use this sample to investigate the origin of this line in star-forming galaxies. We show that star-forming galaxies and galaxies dominated by an active galactic nucleus form clearly separated branches in the He II {\lambda}4686/H{\beta} versus [N II] {\lambda}6584/H{\alpha} diagnostic diagram and derive an empirical classification scheme which separates the two classes. We also present an analysis of the physical properties of 189 star forming galaxies with strong He II {\lambda}4686 emissions. These star-forming galaxies provide constraints on the hard ionizing continuum of massive stars. To make a quantitative comparison with observation we use photoionization models and examine how different stellar population models affect the predicted He II {\lambda}4686 emission. We confirm previous findings that the models can predict He II {\lambda}4686 emission only for instantaneous bursts of 20% solar metallicity or higher, and only for ages of ~ 4 - 5 Myr, the period when the extreme-ultraviolet continuum is dominated by emission from Wolf-Rayet stars. We find however that 83 of the star-forming galaxies (40%) in our sample do not have Wolf-Rayet features in their spectra despite showing strong nebular He II {\lambda}4686 emission. We discuss possible reasons for this and possible mechanisms for the He II {\lambda}4686 emission in these galaxies.
The role of environmentally induced gas stripping in driving galaxy evolution in groups remains poorly understood. Here we present extensive Chandra and Very Large Array mosaic observations of the hot and cold interstellar medium within the members of the nearby, X-ray bright NGC 2563 group, a prime target for studies of the role of gas stripping and interactions in relatively small host halos. Our observations cover nearly all group members within a projected radius of 1.15 Mpc (~1.4 R_vir) of the group center, down to a limiting X-ray luminosity and HI mass of 3e39 erg/s and 2e8 M_sun, respectively. The X-ray data are consistent with efficient ram pressure stripping of the hot gas halos of early-type galaxies near the group core, but no X-ray tails are seen and the limited statistics preclude strong conclusions. The HI results suggest moderate HI mass loss from the group members when compared to similar field galaxies. Six of the 20 HI-detected group members show HI evidence of ongoing interactions with other galaxies or with the intragroup medium. Suggestive evidence is further seen for galaxies with close neighbors in position-velocity space to show relatively low HI content, consistent with tidal removal of HI. The results thus indicate removal of both hot and cold gas from the group members via a combination of ram pressure stripping and tidal interactions. We also find that 16 of the 20 HI detections occur on one side of the group, reflecting an unusual morphological segregation whose origin remains unclear.
We present broadband multi-wavelength observations of GRB 080310 at redshift z = 2.43. This burst was bright and long-lived, and unusual in having extensive optical and near IR follow-up during the prompt phase. Using these data we attempt to simultaneously model the gamma-ray, X-ray, optical and IR emission using a series of prompt pulses and an afterglow component. Initial attempts to extrapolate the high energy model directly to lower energies for each pulse reveal that a spectral break is required between the optical regime and 0.3 keV to avoid over predicting the optical flux. We demonstrate that afterglow emission alone is insufficient to describe all morphology seen in the optical and IR data. Allowing the prompt component to dominate the early-time optical and IR and permitting each pulse to have an independent low energy spectral indices we produce an alternative scenario which better describes the optical light curve. This, however, does not describe the spectral shape of GRB 080310 at early times. The fit statistics for the prompt and afterglow dominated models are nearly identical making it difficult to favour either. However one enduring result is that both models require a low energy spectral index consistent with self absorption for at least some of the pulses identified in the high energy emission model.
We present far-infrared (FIR) analysis of 68 Brightest Cluster Galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500um), we calculate the obscured star formation rate (SFR). 22(+6.2,-5.3)% of the BCGs are detected in the far-infrared, with SFR= 1-150 M_sun/yr. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Ha emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2x10^11 L_sun), only a small (<0.4 mag) reddening correction is required for SFR(Ha) to agree with SFR_FIR. The relatively low Ha extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate form normal stellar mass loss.
We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ~ 50, 000 galaxy systems. We derive accurate projected number density profiles of satellites down to 4 magnitudes fainter than their primaries. We find the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and the colour and concentration of their primaries.
We simulate the viscous evolution of an accretion disc around a spinning black hole. In general any such disc is misaligned, and warped by the Lense-Thirring effect. Unlike previous studies we use effective viscosities constrained to be consistent with the internal fluid dynamics of the disc. We find that nonlinear fluid effects, which reduce the effective viscosities in warped regions, can promote the breaking of the disc into two distinct planes. This occurs when the Shakura & Sunyaev dimensionless viscosity parameter alpha is <~ 0.3 and the initial angle of misalignment between the disc and hole is >~ 45 degrees. The break can be a long-lived feature, propagating outwards in the disc on the usual alignment timescale, after which the disc is fully co- or counter-aligned with the hole. Such a break in the disc may be significant in systems where we know the inclination of the outer accretion disc to the line of sight, such as some X-ray binaries: the inner disc, and so any jets, may be noticeably misaligned with respect to the orbital plane.
We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high precision proper motions (Kallivayalil et al. 2006a,b). Given the new velocities, cosmological simulations of structure formation favor a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy (Boylan-Kolchin et al. 2011, Busha et al. 2011). We illustrate here that the observed irregular morphology and internal kinematics of the MCs (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. This picture further supports a first infall scenario (Besla et a. 2007). In particular, we demonstrate that the Magellanic Stream, a band of HI gas trailing behind the MCs 150 degrees across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally explained by a recent direct collision with the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars (de Vaucouleurs & Freeman 1972), the majority of which are not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are important not only for the evolution of dwarf galaxies but also have direct consequences for the buildup of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type halos.
We develop a radiation pressure-balanced model for the interstellar medium of high-redshift galaxies that describes many facets of galaxy formation at z>~6, including star formation rates and distributions and gas accretion onto central black holes. We first show that the vertical gravitational force in the disk of such a model is dominated by the disk self-gravity but that both radiation pressure on dust grains and turbulent pressure from dense clumps and disk instabilities are negligible compared with the radiation pressure of starlight on gas. Constraining our model to reproduce the UV luminosity function of Lyman-break galaxies (LBGs), we limit the available parameter-space to wind mass-loading factors 1--4 times the canonical value for momentum-driven winds. We then focus our study by exploring the effects of different angular momentum transport mechanisms in the galactic disk and find that viscosity driven by gravitational torques, such as from linear spiral waves or non-linear orbit crossings, can build up black hole masses by z=6 consistent with canonical M-sigma relations with a duty cycle of unity, while infall mediated by a local viscosity such as in an alpha-disk results in negligible BH accretion. Both gravitational torque models produce X-ray emission from active galactic nuclei in high redshift LBGs in excess of the estimated contribution from high-mass X-ray binaries and consistent with a recent analysis of deep Chandra observations by Cowie et al. We find that future observations with larger sample sizes may be able to distinguish between these different angular momentum transport mechanisms.
In this article, I summarize and discuss the body of evidence which has accumulated in favor of dark matter in the form of approximately 10 GeV particles. This evidence includes the spectrum and angular distribution of gamma rays from the Galactic Center, the synchrotron emission from the Milky Way's radio filaments, the diffuse synchrotron emission from the Inner Galaxy (the "WMAP Haze") and low-energy signals from the direct detection experiments DAMA/LIBRA, CoGeNT and CRESST-II. This collection of observations can be explained by a relatively light dark matter particle with an annihilation cross section consistent with that predicted for a simple thermal relic (sigma v ~ 10^-26 cm^3/s) and with a distribution in the halo of the Milky Way consistent with that predicted from simulations. Astrophysical explanations for the gamma ray and synchrotron signals, in contrast, have not been successful in accommodating these observations. Similarly, the phase of the annual modulation observed by DAMA/LIBRA (and now supported by CoGeNT) is inconsistent with all known or postulated modulating backgrounds, but are in good agreement with expectations for dark matter scattering. This scenario is consistent with all existing indirect and collider constraints, as well as the constraints placed by CDMS. Consistency with xenon-based experiments can be achieved if the response of liquid xenon to very low-energy nuclear recoils is somewhat suppressed relative to previous evaluations, or if the dark matter possesses different couplings to protons and neutrons.
We have recently developed an extended merger-tree model that efficiently follows hierarchical evolution of galaxy clusters and provides a quantitative description of both their dark matter and gas properties. We employed this diagnostic tool to calculate the thermal SZ power spectrum and cluster number counts, accounting explicitly for uncertainties in the relevant statistical and intrinsic cluster properties, such as the halo mass function and the gas equation of state. Results of these calculations are compared with those obtained from a direct analytic treatment and from hydrodynamical simulations. We show that under certain assumptions on the gas mass fraction our results are consistent with the latest SPT measurement. Our approach can be particularly useful in predicting cluster number counts and their dependence on cluster and cosmological parameters.
The existence of predominantly cold non-baryonic dark matter is unambiguously demonstrated by several observations (e.g., structure formation, big bang nucleosynthesis, gravitational lensing, and rotational curves of spiral galaxies). A candidate well motivated by particle physics is a weakly interacting massive particle (WIMP). Self-annihilating WIMPs would affect the stellar evolution especially in the early universe. Stars powered by self-annihilating WIMP dark matter should possess different properties compared with standard stars. While a direct detection of such dark matter powered stars seems very challenging, their cumulative emission might leave an imprint in the diffuse metagalactic radiation fields, in particular in the mid-infrared part of the electromagnetic spectrum. In this work the possible contributions of dark matter powered stars (dark stars; DSs) to the extragalactic background light (EBL) are calculated. It is shown that existing data and limits of the EBL intensity can already be used to rule out some DS parameter sets.
We have collated multiplicity data for five clusters (Taurus, Chamaeleon I,
Ophiuchus, IC348, and the Orion Nebula Cluster). We have applied the same mass
ratio (flux ratios of delta K <= 2.5) and primary mass cuts (~0.1-3.0 Msun) to
each cluster and therefore have directly comparable binary statistics for all
five clusters in the separation range 62-620 au, and for Taurus, Chamaeleon I,
and Ophiuchus in the range 18-830 au. We find that the trend of decreasing
binary fraction with cluster density is solely due to the high binary fraction
of Taurus, the other clusters show no obvious trend over a factor of nearly 20
in density.
With N-body simulations we attempt to find a set of initial conditions that
are able to reproduce the density, morphology and binary fractions of all five
clusters. Only an initially clumpy (fractal) distribution with an initial total
binary fraction of 73 per cent (17 per cent in the range 62-620 au) is able to
reproduce all of the observations (albeit not very satisfactorily). Therefore,
if star formation is universal the initial conditions must be clumpy and with a
high (but not 100 per cent) binary fraction. This could suggest that most
stars, including M-dwarfs, form in binaries.
The issue of which stars may reach the conditions of electron/positron pair formation instability is of importance to understand the final evolution both of the first stars and of contemporary stars. The criterion to enter the pair instability regime in density and temperature is basically controlled by the mass of the oxygen core. The main sequence masses that produce a given oxygen core mass are, in turn, dependent on metallicity, mass loss, and convective and rotationally-induced mixing. We examine the evolution of massive stars to determine the minimum main sequence mass that can encounter pair-instability effects, either a pulsational pair instability (PPISN) or a full-fledged pair-instability supernova (PISN). We concentrate on zero-metallicity stars with no mass loss subject to the Schwarzschild criterion for convective instability, but also explore solar metallicity and mass loss and the Ledoux criterion. As expected, for sufficiently strong rotationally-induced mixing, the minimum main sequence mass is encountered for conditions that induce effectively homogeneous evolution such that the original mass is converted almost entirely to helium and then to oxygen. For this case, we find that the minimum main sequence mass is about 40Msun to encounter PPISN and about 65Msun to encounter a PISN. The implications of these results for the first stars and for contemporary supernovae is discussed.
We present a catalogue of 12.2-GHz methanol masers detected towards 6.7-GHz methanol masers observed in the unbiased Methanol Multibeam (MMB) survey in the longitude range 330\circ (through 360\circ) to 10\circ. This is the first portion of the catalogue which, when complete, will encompass all of the MMB detections. We report the detection of 184 12.2-GHz sources towards 400 6.7-GHz methanol maser targets, equating to a detection rate of 46 per cent. Of the 184 12.2-GHz detections, 117 are reported here for the first time. We draw attention to a number of 'special' sources, particularly those with emission at 12.2-GHz stronger than their 6.7-GHz counterpart and conclude that these unusual sources are not associated with a specific evolutionary stage.
This document describes a convention for compressing n-dimensional images and storing the resulting byte stream in a variable-length column in a FITS binary table. The FITS file structure outlined here is independent of the specific data compression algorithm that is used. The implementation details for 4 widely used compression algorithms are described here, but any other compression technique could also be supported by this convention. The general principle used in this convention is to first divide the n-dimensional image into a rectangular grid of subimages or 'tiles'. Each tile is then compressed as a block of data, and the resulting compressed byte stream is stored in a row of a variable length column in a FITS binary table. By dividing the image into tiles it is generally possible to extract and uncompress subsections of the image without having to uncompress the whole image.
This document describes a convention for compressing FITS binary tables that is modeled after the FITS tiled-image compression method (White et al. 2009) that has been in use for about a decade. The input table is first optionally subdivided into tiles, each containing an equal number of rows, then every column of data within each tile is compressed and stored as a variable-length array of bytes in the output FITS binary table. All the header keywords from the input table are copied to the header of the output table and remain uncompressed for efficient access. The output compressed table contains the same number and order of columns as in the input uncompressed binary table. There is one row in the output table corresponding to each tile of rows in the input table. In principle, each column of data can be compressed using a different algorithm that is optimized for the type of data within that column, however in the prototype implementation described here, the gzip algorithm is used to compress every column.
The checksum keywords described here provide an integrity check on the information contained in FITS HDUs. (Header and Data Units are the basic components of FITS files, consisting of header keyword records followed by optional associated data records). The CHECKSUM keyword is defined to have a value that forces the 32-bit 1's complement checksum accumulated over all the 2880-byte FITS logical records in the HDU to equal negative 0. (Note that 1's complement arithmetic has both positive and negative zero elements). Verifying that the accumulated checksum is still equal to -0 provides a fast and fairly reliable way to determine that the HDU has not been modified by subsequent data processing operations or corrupted while copying or storing the file on physical media.
The dynamics of magnetic field decay with Hall drift is investigated. Assuming that axisymmetric magnetic fields are located in a spherical crust with uniform conductivity and electron number density, long-term evolution is calculated up to Ohmic dissipation. The nonlinear coupling between poloidal and toroidal components is explored in terms of their energies and helicity. Nonlinear oscillation by the drift in strongly magnetized regimes is clear only around the equipartition between two components. Significant energy is transferred to the poloidal component when the toroidal component initially dominates. However, the reverse is not true. Once the toroidal field is less dominant, it quickly decouples due to a larger damping rate. The polar field at the surface is highly distorted from the initial dipole during the Hall drift timescale, but returns to the initial dipole in a longer dissipation timescale, since it is the least damped one.
The Cold Dark Matter paradigm predicts vast numbers of dark matter sub-halos to be orbiting in galactic halos. The sub-halos are detectable through the gaps they create gaps in stellar streams. The gap-rate is an integral over the density of sub-halos, their mass function, velocity distribution and the dynamical age of the stream. The rate of visible gap creation is a function of the width of the stream. The available data for four streams: the NW stream of M31, the Pal~5 stream, the Orphan Stream and the Eastern Banded Structure, are compared to the LCDM predicted relation. We find a remarkably good agreement, although there remains much to be done to improve the quality of the result. The narrower streams require that there is a total population of order 10^5 sub-halos above 10^5 M_sun to create the gaps.
Observed spectra of R Coronae Borealis (RCB) and hydrogen-deficient carbon (HdC) stars are analyzed by synthesizing the C2 Swan bands (1,0), (0,0), and (0,1) using our detailed line list and the Uppsala model atmospheres. The (0,1) and (0,0) C2 bands are used to derive the 12C abundance, and the (1,0) 12C13C band to determine the 12C/13C ratios. The carbon abundance derived from the C2 Swan bands is about the same for the adopted models constructed with different carbon abundances over the range: 8.5 (C/He = 0.1%), to 10.5 (C/He = 10%). Carbon abundances derived from C I lines are about a factor of 4 lower than the carbon abundance of the adopted model atmosphere over the same C/He interval, as reported by Asplund et al. (2000), who dubbed the mismatch between adopted and derived C abundance the 'carbon problem'. In principle, the carbon abundances obtained from C2 Swan bands and that assumed for the model atmosphere can be equated for a particular choice of C/He that varies from star to star. Then, the carbon problem for C2 bands is eliminated. However, such C/He ratios are in general less than those of the EHe stars, the seemingly natural relatives to the RCB and HdC stars. A more likely solution to the C2 carbon problem may lie in a modification of the model atmosphere's temperature structure. The derived carbon abundances and the 12C/13C ratios are discussed in light of the double degenerate (DD) and the final flash (FF) scenarios.
X-ray loud M dwarfs are a major source of by-product (contamination) in the X-ray band of the multiwavelength quasar survey (MWQS). As a by-product, the low dispersion spectra of 22 M dwarfs are obtained in which the spectra of 16 sources are taken for the first time. The spectral types and distance of the sample are given based on spectral indices CaH2, CaH3, and TiO5. The parameter {\zeta}TiO/CaH is calculated to make the metallicity class separation among dwarfs, subdwarfs and extreme subdwarfs. We also discuss the distributions in the diagrams of Log(Lx/Lbol) versus spectral type and infrared colors.
We study relations between the X-ray luminosity, orbital period and absolute near-infrared magnitude of persistent low-mass X-ray binaries (LMXBs). We show that often optical and near-infrared spectral energy distribution of LMXBs can be adequately described by a simple model of an accretion disc and a secondary star reprocessing X-ray emission of a central compact object. This gives us an evidence that using an X-ray luminosity and an absolute infrared magnitude of a persistent LMXB one can make reliable estimate of its orbital period. Using a sample of well-known LMXBs, we have constructed a correlation of L_X, P_orb and M_K values which can be approximated by a straight line with the RMS scatter at the level of ~0.3 mag. Such a correlation, being to some extent an analogous to the correlation, found by van Paradijs & McClintock 1994, might be helpful for future population studies especially in the light of forthcoming surveys of the Galaxy in X-ray and infrared spectral domains.
This chapter reviews the recent progress made mainly during the last two decades on wave turbulence in magnetized plasmas (MHD, Hall MHD and electron MHD) in the incompressible and compressible cases. The emphasis is made on homogeneous and anisotropic turbulence which usually provides the best theoretical framework to investigate space and laboratory plasmas. The interplanetary medium and the solar atmosphere are presented as two examples of media where anisotropic wave turbulence is relevant. The most important results of wave turbulence are reported and discussed in the context of space and simulated magnetized plasmas. Important issues and possible spurious interpretations are eventually discussed.
The mass function of superclusters is derived fully analytically with the help of the extended excursion set theory and shown to be in excellent agreement with the numerical results from various publicly available N-body simulation database. We introduce a new multi-dimensional barrier model in which the formation of superclusters occurs when the initial shear eigenvalues that perform non-Markovian random walks enter a clustering zone surrounded by one re ecting and two absorbing barriers. The multi-dimensional barrier heights are determined from the first order Lagrangian perturbation theory and found to be independent of redshift and background cosmology. With the help of our analytic model for the supercluster mass function, the relative abundance of the rich superclusters is analytically evaluated at a given epoch and found to be sensitive to the growth rate of the cosmic web. Our result implies that the relative abundance of the rich superclusters at a given epoch may be useful as a cosmological test of gravity.
The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely the existence of exoplanets not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing planetary system, created in a tight cluster, and is captured as a new planet. In the present work we study this interaction process by integrating trajectories of planet-sized bodies, which encounter a binary system consisting of a Jupiter-sized planet revolving around a Sun-like star. To simplify the problem we assume coplanar orbits for the bound and the free-floating planet and an initially parabolic orbit for the free-floating planet. By calculating the uncertainty exponent, a quantity that measures the dependence of the final state of the system on small changes of the initial conditions, we show that the interaction process is a fractal classical scattering. In this way we see that the statistical approach we follow to tackle the problem is justified. The possible final outcomes of this interaction are only four, namely flyby, planet exchange, capture or disruption. We give the probability of each outcome as a function of the incoming planet's mass. %We also give the final %elements of the orbits of both planets as functions of the masses %and the elements of the orbit of the initially bound planet. We find that the probability of exchange or capture (in prograde as well as retrograde orbits and for very long times) is non-negligible, a fact that might explain the possible future observations of planetary systems with orbits either retrograde or tight and highly eccentric.
Russell (1948) famously described eclipses as the "royal road" to stellar astrophysics. From photometric and spectroscopic observations it is possible to measure the masses and radii (to 1% or better!), and thus surface gravities and mean densities, of stars in eclipsing binary systems using nothing more than geometry. Adding an effective temperature subsequently yields luminosity and then distance (or vice versa) to high precision. This wealth of directly measurable quantities makes eclipsing binaries the primary source of empirical information on the properties of stars, and therefore a cornerstone of stellar astrophysics. In this review paper I summarise the current standing of eclipsing binary research, present an overview of useful analysis techniques, and conclude with a glance to the future.
Neutrino emission rate from baryonic matter in neutron stars via weak neutral vector interaction is computed up to order O(v_F^6), where v_F is the Fermi velocity in units of speed of light. The vector current polarization tensors are evaluated with full vertices which include re-summed series in the particle-hole channel. The neutrino emissivity is enhanced compared to the O(v_F^4) order up to 10% for values v_F < 0.4 characteristic to baryons in compact stars.
Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and here, we investigate the kinematic structure and energy content. We find that IRDCs overall display organized velocity fields, with only localized disruptions due to embedded star formation. The local effects seen in NH3 emission are not high velocity outflows but rather moderate (few km/s) increases in the line width that exhibit maxima near or coincident with the mid-infrared emission tracing protostars. These line width enhancements could be the result of infall or (hidden in NH3 emission) outflow. Not only is the kinetic energy content insufficient to support the IRDCs against collapse, but also the spatial energy distribution is inconsistent with a scenario of turbulent cloud support. We conclude that the velocity signatures of the IRDCs in our sample are due to active collapse and fragmentation, in some cases augmented by local feedback from stars.
The origin, structure and evolution of the small gas cloud, G2, is investigated, that is on an orbit almost straight into the Galactic central supermassive black hole (SMBH). G2 is a sensitive probe of the hot accretion zone of Sgr A*, requiring gas temperatures and densities that agree well with models of captured shock-heated stellar winds. Its mass is equal to the critical mass below which cold clumps would be destroyed quickly by evaporation. Its mass is also constrained by the fact that at apocenter its sound crossing timescale was equal to its orbital timescale. Our numerical simulations show that the observed structure and evolution of G2 can be well reproduced if it formed in pressure equilibrium with the surrounding in 1995 at a distance from the SMBH of 7.6e16 cm. If the cloud would have formed at apocenter in the 'clockwise' stellar disk as expected from its orbit, it would be torn into a very elongated spaghetti-like filament by 2011 which is not observed. This problem can be solved if G2 is the head of a larger, shell-like structure that formed at apocenter. Our numerical simulations show that this scenario explains not only G2's observed kinematical and geometrical properties but also the Br_gamma observations of a low surface brightness gas tail that trails the cloud. In 2013, while passing the SMBH G2 will break up into a string of droplets that within the next 30 years mix with the surrounding hot gas and trigger cycles of AGN activity.
Motivated by recent results on lognormal statistics showing that the moment hierarchy of a lognormal variable completely fails at capturing its information content in the large variance regime, we discuss in this work the inadequacy of the hierarchy of correlation functions to describe a correlated lognormal field, which provides a roughly accurate description of the non-linear cosmological matter density field. We present families of fields having the same hierarchy of correlation functions than the lognormal field at all orders. This explicitly demonstrates the little studied though known fact that the correlation function hierarchy never provides a complete description of a lognormal field, and that it fails to capture information in the non-linear regime, where other simple observables are left totally unconstrained. We discuss why perturbative, Edgeworth-like approaches to statistics in the non-linear regime, common in cosmology, can never reproduce or predict that effect, and why it is however generic for tailed fields, hinting at a breakdown of the perturbation theory based on the field fluctuations. We make a rough but successful quantitative connection to N-body simulations results, that showed that the spectrum of the log-density field carries more information than the spectrum of the field entering the non-linear regime.
The disruption of a binary star system by the massive black hole at the Galactic Centre, SgrA*, can lead to the capture of one star around SgrA* and the ejection of its companion as a hypervelocity star (HVS). We consider the possibility that these stars may have planets and study the dynamics of these planets. Using a direct $N$-body integration code, we simulated a large number of different binary orbits around SgrA*. For some orbital parameters, a planet is ejected at a high speed. In other instances, a HVS is ejected with one or more planets orbiting around it. In these cases, it may be possible to observe the planet as it transits the face of the star. A planet may also collide with its host star. In such cases the atmosphere of the star will be enriched with metals. In other cases, a planet is tidally disrupted by SgrA*, leading to a bright flare.
Dust settling and grain growth are the first steps in the planet-formation
process in protoplanetary disks. These disks are observed around stars with
different spectral types, and there are indications that the disks around lower
mass stars are significantly flatter, which could indicate that they settle and
evolve faster, or in a different way.
We aim to test this assumption by modeling the median spectral energy
distributions (SEDs) of three samples of protoplanetary disks: around Herbig
stars, T Tauri stars and brown dwarfs. We focus on the turbulent mixing
strength to avoid a strong observational bias from disk and stellar properties
that depend on stellar mass.
We generated SEDs with the radiative transfer code MCMax, using a hydrostatic
disk structure and settling the dust in a self-consistent way with the
alpha-prescription to probe the turbulent mixing strength.
We are able to fit all three samples with a disk with the same input
parameters, scaling the inner edge to the dust evaporation radius and disk mass
to millimeter photometry. The Herbig stars require a special treatment for the
inner rim regions, while the T-Tauri stars require viscous heating, and the
brown dwarfs lack a good estimate of the disk mass because only few millimeter
detections exist.
We find that the turbulent mixing strength does not vary across the stellar
mass range for a fixed grain size distribution and gas-to-dust ratio. Regions
with the same temperature have a self-similar vertical structure independent of
stellar mass, but regions at the same distance from the central star appear
more settled in disks around lower mass stars. We find a relatively low
turbulent mixing strength of alpha = 10^(-4) for a standard grain size
distribution, but our results are also consistent with alpha = 0.01 for a grain
size distribution with fewer small grains or a lower gas-to-dust ratio.
We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^(\nu) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^(\nu) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field.
[Abridged]We present analysis of a pair of unusually energetic coronal hard X-ray (HXR) sources detected by RHESSI during the impulsive phase of an X3.9 class solar flare on 2003 November 3, which simultaneously shows two intense footpoint (FP) sources. A distinct loop top (LT) coronal source is detected up to ~150 keV and a second (upper) coronal source up to ~80 keV. These photon energies are much higher than commonly observed in coronal sources and pose grave modeling challenges. The LT source in general appears higher in altitude with increasing energy and exhibits a more limited motion compared to the expansion of the thermal loop. The high energy LT source shows an impulsive time profile and its nonthermal power law spectrum exhibits soft-hard-soft evolution during the impulsive phase, similar to the FP sources. The upper coronal source exhibits an opposite spatial gradient and a similar spectral slope compared to the LT source. These properties are consistent with the model of stochastic acceleration of electrons by plasma waves or turbulence. However, the LT and FP spectral index difference (varying from ~0-1) is much smaller than commonly measured and than that expected from a simple stochastic acceleration model. Additional confinement or trapping mechanisms of high energy electrons in the corona are required. Comprehensive modeling including both kinetic effects and the macroscopic flare structure may shed light on this behavior. These results highlight the importance of imaging spectroscopic observations of the LT and FP sources up to high energies in understanding electron acceleration in solar flares. Finally, we show that the electrons producing the upper coronal HXR source may very likely be responsible for the type III radio bursts at the decimetric/metric wavelength observed during the impulsive phase of this flare.
FS Tau B is one of the few T Tauri stars that possess a jet and a counterjet as well as an optically-visible cavity wall. We obtained images and spectra of its jet-cavity system in the near-infrared H and K bands using Subaru/IRCS and detected the jet and the counterjet in the [Fe II] 1.644 \mu m line for the first time. Within the inner 2" the blueshifted jet is brighter, whereas beyond ~ 5" the redshifted counterjet dominates the [Fe II] emission. The innermost blueshifted knot is spectrally resolved to have a large line width of ~ 110 km/s, while the innermost redshifted knot appears spectrally unresolved. The velocity ratio of the jet to the counterjet is ~ 1.34, which suggests that FS Tau B is driving an asymmetric jet, similar to those found in several T Tauri Stars. Combining with optical observations in the literature, we showed that the blueshifted jet has lower density and higher excitation than the redshifted counterjet. We suggest that the asymmetry in brightness and velocity is the manifestation of a bipolar outflow driving at different mass-loss rates, while maintaining balance of linear momentum. A full explanation to the asymmetry in the FS Tau B system awaits detail modeling and further investigation of the kinematic structure of the wind-associated cavity walls.
An extension of the generalized relativistic mean-field (gRMF) model with density dependent couplings is introduced in order to describe thermodynamical properties and the composition of dense nuclear matter for astrophysical applications. Bound states of light nuclei and two-nucleon scattering correlations are considered as explicit degrees of freedom in the thermodynamical potential. They are represented by quasiparticles with medium dependent properties. The model interpolates between the correct low-density limit, the model independent virial equation of state (VEoS), and the RMF description around nuclear saturation density where clusters are dissolved. A comparison between the fugacity expansions of the VEoS and the gRMF model provides consistency relations between the quasiparticles properties, the nucleon-nucleon scattering phase shifts and the meson-nucleon couplings of the gRMF model at zero density. Relativistic effects are found to be important at temperatures that are typical in astrophysical applications. The case of neutron matter is studied in detail.
In this paper we focus on the gravitational thermodynamics of the far future. Cosmological observations suggest that most matter will be diluted away by the cosmological expansion, with the rest collapsing into supermassive black holes. The likely future state of our local universe is a supermassive black hole slowly evaporating in an empty universe dominated by a positive cosmological constant. We describe some overlooked features of how the cosmological horizon responds to the black hole evaporation. The presence of a black hole depresses the entropy of the cosmological horizon by an amount proportional to the geometric mean of the entropies of the black hole and cosmological horizons. As the black hole evaporates and loses its mass in the process, the total entropy increases obeying the second law of thermodynamics. The entropy is produced by the heat from the black hole flowing across the extremely cold cosmological horizon. Once the evaporation is complete, the universe becomes empty de Sitter space that (in the presence of a true cosmological constant) is the maximum entropy thermodynamic equilibrium state. We propose that flat Minkowski space is an improper limit of this process which obscures the thermodynamics. The cosmological constant should be regarded not only as an energy scale, but also as a scale for the maximum entropy of a universe. In this context, flat Minkowski space is indistinguishable from de Sitter with extremely small cosmological constant, yielding a divergent entropy. This introduces an unregulated infinity in black hole thermodynamics calculations, giving possibly misleading results.
We compare the performance of a recently proposed empirical climate model based on astronomical harmonics against all available general circulation climate models (GCM) used by the IPCC (2007) to interpret the 20th century global surface temperature. The proposed model assumes that the climate is resonating with, or synchronized to a set of natural harmonics that have been associated to the solar system planetary motion, mostly determined by Jupiter and Saturn. We show that the GCMs fail to reproduce the major decadal and multidecadal oscillations found in the global surface temperature record from 1850 to 2011. On the contrary, the proposed harmonic model is found to well reconstruct the observed climate oscillations from 1850 to 2011, and it is able to forecast the climate oscillations from 1950 to 2011 using the data covering the period 1850-1950, and vice versa. The 9.1-year cycle is shown to be likely related to a decadal Soli/Lunar tidal oscillation, while the 10-10.5, 20-21 and 60-62 year cycles are synchronous to solar and heliospheric planetary oscillations. Finally, we show how the presence of these large natural cycles can be used to correct the IPCC projected anthropogenic warming trend for the 21st century. By combining this corrected trend with the natural cycles, we show that the temperature may not significantly increase during the next 30 years mostly because of the negative phase of the 60-year cycle. The same IPCC projected anthropogenic emissions would imply a global warming by about 0.3-1.2 K by 2100, contrary to the IPCC 1.0-3.6 K projected warming. The results of this paper reinforce previous claims that the relevant physical mechanisms that explain the detected climatic cycles are still missing in the current GCMs and that climate variations at the multidecadal scales are astronomically induced and, in first approximation, can be forecast.
Links to: arXiv, form interface, find, astro-ph, recent, 1201, contact, help (Access key information)
We present a study of the ICM of the galaxy cluster CIZA J2242.8+5301 using deep XMM-Newton observations. The cluster hosts an extremely elongated (2 Mpc), narrow (~50 kpc) radio relic that has been nicknamed the "Sausage". Additionally, a counter-relic is also present, along with a faint, extended radio halo. We have searched for evidence of shock fronts in the surface brightness, temperature, density, and pressure, and we studied the cluster morphology using power ratios. The surface brightness profiles to the north and south of the centre are almost identical in shape. This symmetry supports the hypothesis that the two merging clusters have almost equal masses and a small impact parameter. The ICM on the inner side of the relics (the side towards the cluster centre), has a relatively low temperature of ~5 keV and only jumps to temperatures >10 keV after about 500 kpc. The jumps in temperature and pressure coincide with two symmetric "bumps" in the X-ray surface brightness profiles. We discuss possible causes for this. The inner shocks do not cause detectable radio relics. Surprisingly, the temperatures and pressures determined from the spectral fits depend strongly on the assumed abundance table and the absorption model. We speculate that this is characteristic of all clusters at low Galactic latitudes.
Observations of ionized and neutral gas outflows in radio-galaxies (RGs) suggest that AGN radio jet feedback has a galaxy-scale impact on the host ISM, but it is still unclear how the molecular gas is affected. We present deep Spitzer IRS spectroscopy of 8 RGs that show fast HI outflows. All of these HI-outflow RGs have bright H2 mid-IR lines that cannot be accounted for by UV or X-ray heating. This suggests that the radio jet, which drives the HI outflow, is also responsible for the shock-excitation of the warm H2 gas. In addition, the warm H2 gas does not share the kinematics of the ionized/neutral gas. The mid-IR ionized gas lines are systematically broader than the H2 lines, which are resolved by the IRS (with FWHM up to 900km/s) in 60% of the detected H2 lines. In 5 sources, the NeII line, and to a lesser extent the NeIII and NeV lines, exhibit blue-shifted wings (up to -900km/s with respect to the systemic velocity) that match the kinematics of the outflowing HI or ionized gas. The H2 lines do not show broad wings, except tentative detections in 3 sources. This shows that, contrary to the HI gas, the H2 gas is inefficiently coupled to the AGN jet-driven outflow of ionized gas. While the dissipation of a small fraction (<10%) of the jet kinetic power can explain the dynamical heating of the molecular gas, our data show that the bulk of the warm molecular gas is not expelled from these galaxies.
Recent high resolution spectroscopic analysis of nearby FGK stars suggests that a high C/O ratio of greater than 0.8, or even 1.0, is relatively common. Two published catalogs of measurements find C/O>0.8 in 25-30% of systems, and C/O>1.0 in ~6-10% of systems. It has been suggested that in protoplanetary disks with C/O>0.8 that the condensation pathways to refractory planet-making material will differ from what occurred in our solar system, where C/O=0.55. The carbon-rich disks are calculated to make carbon-dominated rocky planets, rather than oxygen-dominated ones, which would be very unlike the Earth. Here we suggest that the derived stellar C/O ratios are overestimated, given the extreme paucity of carbon dwarfs stars (<0.1%) found in large samples of low mass stars, such as from the Sloan Digital Sky Survey. The reason for this overestimation is not immediately apparent, but could be due to the choices of lines used, or limitations of abundance analysis from one-dimensional LTE stellar atmosphere models. Furthermore, from the estimated errors on the measured stellar C/O ratios, we find that the significance of the high C/O tail is weakened, with a true measured fraction of C/O>0.8 in 10-15% of stars, and C/O>1.0 in 1-5%, athough these are still likely overestimates. We suggest that infrared T-dwarf spectra could show how common high C/O is in the stellar neighborhood, as the chemistry and spectra of such objects would differ compared to those with solar-like abundances. We expect that carbon-dominated rocky planets are rarer than others have suggested.
H and He features in photospheric spectra have seldom been used to infer quantitatively the properties of Type IIb, Ib and Ic supernovae (SNe IIb, Ib and Ic) and their progenitor stars. Most radiative transfer models ignored NLTE effects, which are extremely strong especially in the He-dominated zones. In this paper, a comprehensive set of model atmospheres for low-mass SNe IIb/Ib/Ic is presented. Long-standing questions such as how much He can be contained in SNe Ic, where He lines are not seen, can thus be addressed. The state of H and He is computed in full NLTE, including the effect of heating by fast electrons. The models are constructed to represent iso-energetic explosions of the same stellar core with differently massive H/He envelopes on top. The synthetic spectra suggest that 0.06 - 0.14 M_sun of He and even smaller amounts of H suffice for optical lines to be present, unless ejecta asymmetries play a major role. This strongly supports the conjecture that low-mass SNe Ic originate from binaries where progenitor mass loss can be extremely efficient.
Modified gravity theories capable of genuine self-acceleration typically invoke a galileon scalar which mediates a long range force, but is screened by the Vainshtein mechanism on small scales. In such theories, non-relativistic stars carry the full scalar charge (proportional to their mass), while black holes carry none. Thus, for a galaxy free-falling in some external gravitational field, its central massive black hole is expected to lag behind the stars. To look for this effect, and to distinguish it from other astrophysical effects, one can correlate the gravitational pull from the surrounding structure with the offset between the stellar center and the black hole. The expected offset depends on the central density of the galaxy, and ranges up to ~0.1 kpc for small galaxies. The observed offset in M87 cannot be explained by this effect unless the scalar force is significantly stronger than gravity. We also discuss the systematic offset of compact objects from the galactic plane as another possible signature.
Two decades ago "transitional disks" described spectral energy distributions (SEDs) of T Tauri stars with small near-IR excesses, but significant mid- and far-IR excesses. Many inferred this indicated dust-free holes in disks, possibly cleared by planets. Recently, this term has been applied disparately to objects whose Spitzer SEDs diverge from the expectations for a typical full disk. Here we use irradiated accretion disk models to fit the SEDs of 15 such disks in NGC 2068 and IC 348. One group has a "dip" in infrared emission while the others' continuum emission decreases steadily at all wavelengths. We find that the former have an inner disk hole or gap at intermediate radii in the disk and we call these objects "transitional" and pre-transitional" disks, respectively. For the latter group, we can fit these SEDs with full disk models and find that millimeter data are necessary to break the degeneracy between dust settling and disk mass. We suggest the term "transitional" only be applied to objects that display evidence for a radical change in the disk's radial structure. Using this definition, we find that transitional and pre-transitional disks tend to have lower mass accretion rates than full disks and that transitional disks have lower accretion rates than pre-transitional disks. These reduced accretion rates onto the star could be linked to forming planets. Future observations of transitional and pre-transitional disks will allow us to better quantify the signatures of planet formation in young disks.
Impulsive radio-frequency signals from astronomical sources are dispersed by the frequency dependent index of refraction of the interstellar media and so appear as chirped signals when they reach earth. Searches for dispersed impulses have been limited by false detections due to radio frequency interference (RFI) and, in some cases, artifacts of the instrumentation. Many authors have discussed techniques to excise or mitigate RFI in searches for fast transients, but comparisons between different approaches are lacking. This work develops RFI mitigation techniques for use in searches for dispersed pulses, employing data recorded in a "Fly's Eye" mode of the Allen Telescope Array as a test case. We gauge the performance of several RFI mitigation techniques by adding dispersed signals to data containing RFI and comparing false alarm rates at the observed signal-to-noise ratios of the added signals. We find that Huber filtering is most effective at removing broadband interferers, while frequency centering is most effective at removing narrow frequency interferers. Neither of these methods is effective over a broad range of interferers. A method that combines Huber filtering and adaptive interference cancellation provides the lowest number of false positives over the interferers considered here. The methods developed here have application to other searches for dispersed pulses in incoherent spectra, especially those involving multiple beam systems.
Galaxy clusters, the most massive collapsed structures, have been routinely used to determine cosmological parameters. When using clusters for cosmology, the crucial assumption is that they are relaxed. However, subarcminute resolution Sunyaev-Zel'dovich (SZ) effect images compared with high resolution X-ray images of some clusters show significant offsets between the two peaks. We have carried out self-consistent N-body/hydrodynamical simulations of merging galaxy clusters using FLASH to study these offsets quantitatively. We have found that significant displacements result between the SZ and X-ray peaks for large relative velocities for all masses used in our simulations as long as the impact parameters were about 100-250 kpc. Our results suggest that the SZ peak coincides with the peak in the pressure times the line-of-sight characteristic length and not the pressure maximum (as it would for clusters in equilibrium). The peak in the X-ray emission, as expected, coincides with the density maximum of the main cluster. As a consequence, the morphology of the SZ signal and therefore the offset between the SZ and X-ray peaks change with viewing angle. As an application, we compare the morphologies of our simulated images to observed SZ and X-ray images and mass surface densities derived from weak lensing observations of the merging galaxy cluster CL0152-1357. We find that a large relative velocity of 4800 km/s is necessary to explain these observations. We conclude that an analysis of the morphologies of multi-frequency observations of merging clusters can be used to put meaningful constraints on the initial parameters of the progenitors.
We reconstruct the shape of the primordial power spectrum of curvature perturbations in extended cosmological models, including addition of massive neutrinos, extra relativistic species or varying primordial helium abundance, from the latest cosmic microwave background data from the Wilkinson Microwave Anisotropy Probe, the Atacama Cosmology Telescope and the South Pole Telescope. We find that a scale-invariant primordial spectrum is disfavored by the data at 95% confidence level even in the presence of massive neutrinos, however it can lie within the 95% confidence region if the effective number of relativistic species or the primordial helium abundance is allowed to vary freely. The constraints on the extension parameters from WMAP7+ACT+H0+BAO, are the total mass of neutrinos sum(m_nu) < 0.48 eV (95% CL), the effective number of relativistic species N_eff = 4.50 +/- 0.81 and the primordial helium abundance Y_p = 0.303 +/- 0.075. The constraints from WMAP7+SPT+H0+BAO, are sum(m_nu) < 0.45 eV (95% CL), N_eff = 3.86 +/- 0.63 and Y_p = 0.277 +/- 0.050.
We have performed our spectroscopic and multi-color photometric observations of the recurrent nova U Scorpii at the earliest stage of the 2010 outburst. 0.37 days after the discovery of the outburst, we can see broad and prominent emission lines of Balmer series, He I, N II, N III, O I and Mg II on the spectra. The FWHM of H$\alpha$ line yields an expansion velocity of approximately 6200 km/s. This line also accompanies a blue shifted absorption line (so called P Cygni profile). 1.37 days after the discovery, H$\alpha$ line shows a nearly flat-topped profile in contrast to the previous day. From our multi-color photometry, we can see its rapid decline (one magnitude per day) of the brightness in each color band.
We perform a statistical analysis of the temporal and spectral properties of the latest Fermi gamma-ray bursts (GRBs) to revisit the classification of GRBs. We find that the bimodalities of duration and the energy ratio (Epeak/Fluence) and the anti-correlation between spectral hardness (hardness ratio (HR), peak energy and spectral index) and duration (\Delta T) support the long/soft-short/hard classification scheme for Fermi GRBs. The HR - \Delta T anti-correlation strongly depends upon the spectral shape of GRBs and energy bands, and the bursts with the curved spectra in the typical BATSE energy bands show a tighter anti-correlation than those with the power-law spectra in the typical BAT energy bands. This might explain why the HR - \Delta T correlation is not evident for those GRB samples detected by instruments with a narrower/softer energy bandpass. We also analyze the intrinsic energy correlation for the GRBs with measured redshift and well defined peak energies. The current sample suggests Ep,rest=2455 \times (Eiso/10^{52})^{0.59} for short GRBs, significantly different from that for long GRBs. However, both the long and short GRBs comply with the same Ep,rest-Liso correlation.
GRB 090618 was simultaneously detected by Swift-BAT and Fermi-GBM. Its light curve shows two emission episodes consisting of four prominent pulses. The pulse in the first episode (episode A) has a smoother morphology than the three pulses in the second episode (episode B). Using the pulse peak-fit method, we have performed a detailed analysis of the temporal and spectral characteristics of these four pulses and found out that the first pulse (pulse A) exhibits distinctly different properties than the others in episode B (pulses B1, B2 and B3) in the following aspects. (i) Both the pulse width ($w$) and the rise-to-decay ratio of pulse ($r/d$, pulse asymmetry) in GRB 090618 are found to be energy-dependent. The indices of the power-law correlation between $w$ and $E$ for the pulses in episode B however are larger than that in episode A. Moreover the pulses B1, B2 and B3 tend to be more symmetric at the higher energy bands while the pulse A displays a reverse trend. (ii) Pulse A shows a hard-to-soft spectral evolution pattern, while the three pulses in the episode B follow the light curve trend. (iii) Pulse A has a longer lag than the pulses B1, B2 and B3. The mechanism which causes the different pulse characteristics within one single GRB is unclear.
We have compared the oxygen and nitrogen abundances derived from global emission-line SDSS spectra of galaxies using (1) the Te method and (2) two recent strong line calibrations: the ON and NS calibrations. Using the Te method, anomously high N/O abundances ratios have been found in some SDSS galaxies. To investigate this, we have Monte Carlo simulated the global spectra of composite nebulae by a mix of spectra of individual components, based on spectra of well-studied HII regions in nearby galaxies. We found that the Te method results in an underestimated oxygen abundance (and hence in an overestimated nitrogen-to-oxygen ratio) if HII regions with different physical properties contribute to the global spectrum of composite nebulae. This effect is somewhat similar to the small-scale temperature fluctuations in HII regions discussed by Peimbert. Our work thus suggests that the high Te-based N/O abundances ratios found in SDSS galaxies may not be real. However, such an effect is not expected to be present in dwarf galaxies since they have generally an uniform chemical composition. The ON and NS calibrations give O and N abundances in composite nebulae which agree with the mean luminosity-weighted abundances of their components to within 0.2 dex.
I review work on modelling the infrared and submillimetre SEDs of galaxies. The underlying physical assumptions are discussed and spherically symmetric, axisymmetric, and 3-dimensional radiative transfer codes are reviewed. Models for galaxies with Spitzer IRS data and for galaxies in the Herschel-Hermes survey are discussed. Searches for high redshift infrared and submillimetre galaxies, the star formation history, the evolution of dust extinction, and constraints from source-counts, are briefly discussed.
[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.
Planets with several Earth masses and a few day orbital periods have been discovered through radial velocity and transit surveys. Regardless of their formation mechanism, a key evolution issue is the efficiency of their retention near their host stars. If these planets attained their present-day orbits during or shortly after the T Tauri phase of their host stars, a large fraction would have encountered intense stellar magnetic field. Since these planets have a higher conductivity than the atmosphere of their stars, the magnetic flux tube connecting the planet and host star would slip though the envelope of the star faster than across the planet. The induced electro-motive force across the planet's diameter leads to a potential drop which propagates along a flux tube away from the planet with an Alfven speed. The foot of the flux tube sweeps across the stellar surface and the potential drop drives a DC current analogous to that proposed for the Io-Jupiter electrodynamic interaction. The ohmic dissipation of this current produces potentially observable hot spots in the star envelope. The current heats the planet and leads to a Lorrentz torque which drives the planet's orbit to evolve toward circularization and synchronization with the star's spin. The net effect is the damping of the planet's orbital eccentricity. Around slowly (rapidly) spinning stars, this process also causes rocky planets with periods less than a few days to undergo orbital decay (expansion/stagnation) within a few Myr. In principle, this effect can determine the retention efficiency of short-period hot Earths. We also estimate the ohmic dissipation in these planets and show that it can lead to severe structure evolution and potential loss of volatile material. However, these effects may be significantly weakened by the reconnection of the induced field [Slightly shortened abstract].
We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.
We investigate for the first time the effects of a Warm Dark Matter (WDM) power spectrum on the statistical properties of galaxies using a semi-analytic model of galaxy formation. The WDM spectrum we adopt as a reference case is suppressed - compared to the standard Cold Dark Matter (CDM) case - below a cut-off scale ~ 1 Mpc corresponding (for thermal relic WDM particles) to a mass m_X=0.75 keV. This ensures consistency with present bounds provided by the microwave background WMAP data and by the comparison of hydrodynamical N-body simulations with observed Lyman-{\alpha} forest. We run our fiducial semi-analytic model with such a WDM spectrum to derive galaxy luminosity functions (in B, UV, and K bands) and the stellar mass distributions over a wide range of cosmic epochs, to compare with recent observations and with the results in the CDM case. The predicted color distribution of galaxies in the WDM model is also checked against the data. When compared with the standard CDM case, the luminosity and stellar mass distributions we obtain assuming a WDM spectrum are characterized by: i) a flattening of the faint end slope and ii) a sharpening of the cutoff at the bright end for z \lesssim 0.8. We discuss how the former result is directly related to the smaller number of low-mass haloes collapsing in the WDM scenario, while the latter is related to the smaller number of satellite galaxies accumulating in massive haloes at low redshift, thus suppressing the accretion of small lumps on the central, massive galaxies. These results shows how a adopting a WDM power spectrum may contribute to solve two major problems of CDM galaxy formation scenarios, namely, the excess of predicted faint (low mass) galaxies at low and - most of all - high redshifts, and the excess of bright (massive) galaxies at low redshifts.
We combine the estimated metallicities [Fe/H], abundances [\alpha/Fe], positions and motions of a sample of 27,500 local (7<R/kpc<9, 0.5<|z|/kpc<2.5) SDSS/SEGUE G-type dwarf stars to investigate the chemo-orbital properties of the Milky Way's disk around the Sun. When we derive the orbital properties reflecting angular momentum, circularity, and thickness as function of [\alpha/Fe] vs. [Fe/H], we find that there is a smooth variation with [\alpha/Fe], a proxy for age. At the same time, the orbital properties of the old stars with [\alpha/Fe]$\gtrsim$0.25 do show a transition with [Fe/H]: below [Fe/H]$\simeq$-0.6 the orbital angular momentum decreases, and the orbits become significantly non-circular and thicker. Radial migration of stars into the Solar neighborhood would naturally result in a smooth variation in the orbital properties, but the latter old metal-poor stars form a clear challenge, in particular because a basic feature of radial migration is that stars remain on near-circular orbits. When we next select stars on near-circular orbits, we indeed find besides the \alpha-young 'thin-disk' stars a significant contribution to the \alpha-old 'thick-disk' metal-rich stars. However, the remaining \alpha-old 'thick-disk' stars on eccentric orbits, including nearly all old metal-poor stars, are difficult to explain with radial migration alone, but might have formed through early-on gas-rich mergers. We thus find chemo-orbital evidence that the thicker component of the Milky Way disk is not distinct from the thin component as expected from smooth internal evolution through radial migration, except for the old metal-poor stars with different orbital properties which could be part of a distinct thick-disk component formed through an external mechanism.
The role of neutrinos in stars is introduced for students with little prior astrophysical exposure. We begin with neutrinos as an energy-loss channel in ordinary stars and conversely, how stars provide information on neutrinos and possible other low-mass particles. Next we turn to the Sun as a measurable source of neutrinos and other particles. Finally we discuss supernova (SN) neutrinos, the SN 1987A measurements, and the quest for a high-statistics neutrino measurement from the next nearby SN. We also touch on the subject of neutrino oscillations in the high-density SN context.
With the advent of large-collecting-area instruments, the number of objects that can be reached by optical long-baseline interferometry is steadily increasing. We present here a few results on massive binary stars, showing the interest of using this technique for studying the insight of interactions in these systems. Indeed, many massive stars with extended environments host, or are suspected to host, companion stars. These companions could have an important role in shaping the circumstellar environment of the system. These examples provide a view in which binarity could be an ingredient, among many others, for the activity of these stars.
We consider a possibility of identification of sources of cosmic rays (CR) of the energy above 1 TeV via observation of degree-scale extended gamma-ray emission which traces the locations of recent sources in the Galaxy. Such emission in the energy band above 100 GeV is produced by CR nuclei and electrons released by the sources and spreading into the interstellar medium. We use the data from the Fermi gamma-ray telescope to locate the degree-scale 100 GeV gamma-ray sources in the Galactic Plane. We find that a large fraction of the sources is associated to pulsars with spin down age less than ~30 kyr and hence to the recent supernova explosions. This supports the hypothesis of supernova origin of Galactic CRs. We notice that the degree-scale extended emission does not surround shell-like supernova remnants without pulsars. Based on this observation, we argue that the presence of the pulsar is essential for the CR acceleration process.
Using the group catalog obtained from zCOSMOS spectroscopic data and the complementary photometric data from the COSMOS survey, we explore segregation effects occurring in groups of galaxies at intermediate/high redshifts. We built two composite groups at intermediate (0.2 <= z <= 0.45) and high (0.45 < z <= 0.8) redshifts, and we divided the corresponding composite group galaxies into three samples according to their distance from the group center. We explored how galaxy stellar masses and colors - working in narrow bins of stellar masses - vary as a function of the galaxy distance from the group center. We found that the most massive galaxies in our sample (Log(M_gal/M_sun) >= 10.6) do not display any strong group-centric dependence of the fractions of red/blue objects. For galaxies of lower masses (9.8 <= Log(M_gal/M_sun) <= 10.6) there is a radial dependence in the changing mix of red and blue galaxies. This dependence is most evident in poor groups, whereas richer groups do not display any obvious trend of the blue fraction. Interestingly, mass segregation shows the opposite behavior: it is visible only in rich groups, while poorer groups have a a constant mix of galaxy stellar masses as a function of radius. We suggest a simple scenario where color- and mass-segregation originate from different physical processes. While dynamical friction is the obvious cause for establishing mass segregation, both starvation and galaxy-galaxy collisions are plausible mechanisms to quench star formation in groups at a faster rate than in the field. In poorer groups the environmental effects are caught in action superimposed to secular galaxy evolution. Their member galaxies display increasing blue fractions when moving from the group center to more external regions, presumably reflecting the recent accretion history of these groups.
We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction within the critical curve. We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of the stellar mass fraction that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample (lensing velocity dispersion less than 230 km/s), we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with the heavier IMF inferred from the lensing and dynamical analysis of more massive early-type lens galaxies from the SLACS sample, this result provides strong evidence against a universal stellar IMF.
We present an intensive radio and X-ray monitoring campaign on the 2009 outburst of the Galactic black hole candidate X-ray binary H1743-322. With the high angular resolution of the Very Long Baseline Array, we resolve the jet ejection event and measure the proper motions of the jet ejecta relative to the position of the compact core jets detected at the beginning of the outburst. This allows us to accurately couple the moment when the jet ejection event occurred with X-ray spectral and timing signatures. We find that X-ray timing signatures are the best diagnostic of the jet ejection event in this outburst, which occurred as the X-ray variability began to decrease and the Type C quasi-periodic oscillations disappeared from the X-ray power density spectrum. However, this sequence of events does not appear to be replicated in all black hole X-ray binary outbursts, even within an individual source. In our observations of H1743-322, the ejection was contemporaneous with a quenching of the radio emission, prior to the start of the major radio flare. This contradicts previous assumptions that the onset of the radio flare marks the moment of ejection. The jet speed appears to vary between outbursts, with a possible positive correlation with outburst luminosity. The compact core radio jet reactivated on transition to the hard intermediate state at the end of the outburst, and not when the source reached the low hard spectral state. Comparison with the known near-infrared behaviour of the compact jets suggests a gradual evolution of the compact jet power over a few days near the beginning and end of an outburst.
The radius of neutron stars can in principle be measured via the normalisation of a blackbody fitted to the X-ray spectrum during thermonuclear (type-I) X-ray bursts, although few previous studies have addressed the reliability of such measurements. Here we examine the apparent radius in a homogeneous sample of long, mixed H/He bursts from the low-mass X-ray binaries GS 1826-24 and KS 1731-26. The measured blackbody normalisation (proportional to the emitting area) in these bursts is constant over a period of up to 60s in the burst tail, even though the flux (blackbody temperature) decreased by a factor of 60-75% (30-40%). The typical rms variation in the mean normalisation from burst to burst was 3-5%, although a variation of 17% was found between bursts observed from GS 1826-24 in two epochs. A comparison of the time-resolved spectroscopic measurements during bursts from the two epochs shows that the normalisation evolves consistently through the burst rise and peak, but subsequently increases further in the earlier epoch bursts. The elevated normalisation values may arise from a change in the anisotropy of the burst emission, or alternatively variations in the spectral correction factor, f_c, of order 10%. Since burst samples observed from systems other than GS 1826-24 are more heterogeneous, we expect that systematic uncertainties of at least 10% are likely to apply generally to measurements of neutron-star radii, unless the effects described here can be corrected for.
Based on the analysis of the microwave observations at frequency of 2.60 -- 3.80 GHz in a solar X1.3 flare event observed at Solar Broadband RadioSpectrometer in Huairou (SBRS/Huairou) on 2005 July 30, an interesting reversed drifting quasi-periodic pulsating structure (R-DPS) is confirmed. The R-DPS is mainly composed of two drifting pulsating components: one is a relatively slow very short-period pulsation (VSP) with period of about 130 -- 170 ms, the other is a relatively fast VSP with period of about 70 -- 80 ms. The R-DPS has a weak left-handed circular polarization. Based on the synthetic investigations of Reuven Ramaty High Energy Solar Spectroscopic Imaging (RHESSI) hard X-ray, Geostationary Operational Environmental Satellite (GOES) soft X-ray observation, and magnetic field extrapolation, we suggest the R-DPS possibly reflects flaring dynamic processes of the emission source regions.
Methanol masers can provide valuable insight into the processes involved in high-mass star formation, however the local environment in which they form is still unclear. Four primary, yet conflicting, models have emerged to explain the commonly observed methanol maser structures at 6.67 GHz. These suggest that masers trace accretion disks, outflows, shock fronts or disks dominated by infall/outflows. One proposed means of testing these models is through mapping the local magnetic field structures around maser sources, which were predicted to lie parallel to shock and outflows and perpendicular to accretion disks. To follow up this suggestion we have determined magnetic field directions from full polarisation observations of 10 6.67-GHz sources. We find morphology that is parallel to the source structure, indicative of shocks or outflows, in five sources and perpendicular morphology indicative of disks in three. These results do not support any of the expected models and the diverse morphologies observed indicate that the masers could be emitting from different evolutionary stages or environments, or from a common local environment with complex associated magnetic fields. To resolve this conflict we suggest a new approach that will search the simulations of massive star formation, which are just becoming available, for suitable sites for maser emission.
We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]
We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same approach as XML but with a simpler syntax. Together with a specification of PSDF, we provide background and motivation, as well as specific examples in a variety of computer languages. We also offer a web site from which these examples can be retrieved, in order to make it easy to augment existing codes in order to give them the option to produce PSDF output.
There are few data available with which to constrain the thermal history of the intergalactic medium (IGM) following global recombination. Thus far, most constraints flow from analyses of the Cosmic Microwave Background and optical spectroscopy along a few lines of sight. However, direct study of the IGM in emission or absorption against the CMB via the 1S hyperfine transition of Hydrogen would enable broad characterization thermal history and source populations. New generations of radio arrays are in development to measure this line signature. Bright foreground emission and the complexity of instrument calibration models are significant hurdles. How to optimize these is uncertain, resulting in a diversity in approaches. We discuss recent limits on line brightness, array efforts including the new Large Aperture Experiment to Detect the Dark Ages (LEDA), and the next generation Hydrogen Reionization Array (HERA) concept.
The Infrared Array Camera (IRAC) on the Spitzer Space Telescope has revealed that a number of high-mass protostars are associated with extended mid-infrared emission, particularly prominent at 4.5-micron. These are called "Green Fuzzy" emission or "Extended Green Objects". We present color analysis of this emission toward six nearby (d=2-3 kpc) well-studied high-mass protostars and three candidate high-mass protostars identified with the Spitzer GLIMPSE survey. In our color-color diagrams most of the sources show a positive correlation between the [3.6]-[4.5] and [3.5]-[5.8] colors along the extinction vector in all or part of the region. We compare the colors with those of scattered continuum associated with the low-mass protostar L 1527, modeled scattered continuum in cavities, shocked emission associated with low-mass protostars, modeled H2 emission for thermal and fluorescent cases, and modeled PAH emission. Of the emission mechanisms discussed above, scattered continuum provides the simplest explanation for the observed linear correlation. In this case, the color variation within each object is attributed to different foreground extinctions at different positions. Alternative possible emission mechanisms to explain this correlation may be a combination of thermal and fluorescent H2 emission in shocks, and a combination of scattered continuum and thermal H2 emission, but detailed models or spectroscopic follow-up are required to further investigate this possibility. Our color-color diagrams also show possible contributions from PAHs in two objects. However, none of our sample show clear evidence for PAH emission directly associated with the high-mass protostars, several of which should be associated with ionizing radiation. This suggests that those protostars are heavily embedded even at mid-infrared wavelengths.
A study of stellar ages of a sample from the Sloan Digital Sky Survey (SDSS) is presented. The results are consolidated with a set of globular clusters (GCs) and show that this stellar sample is composed by one dominant population of 10-12 Gyr old. This supports the Eggen's scenario claiming that the inner halo of the Milky Way formed rapidly, probably during the collapse of the proto-Galactic cloud.
Energy exchange processes play a crucial role in the early Universe, affecting the thermal balance and the dynamical evolution of the primordial gas. In the present work we focus on the consequences of a non-thermal distribution of the level populations of H$_2$: first, we determine the excitation temperatures of vibrational transitions and the non-equilibrium heat transfer; second, we compare the modifications to chemical reaction rate coefficients with respect to the values obtained assuming local thermodynamic equilibrium; third, we compute the spectral distortions to the cosmic background radiation generated by the formation of H$_2$ in vibrationally excited levels. We conclude that non-equilibrium processes cannot be ignored in cosmological simulations of the evolution of baryons, although their observational signatures remain below current limits of detection. New fits to the equilibrium and non-equilibrium heat transfer functions are provided.
The gamma-ray binary system PSR B1259-63 has recently passed through periastron and has been of particular interest as it was observed by Fermi near the December 2010 periastron passage. The system has been detected at very high energies with H.E.S.S. The most probable production mechanism is inverse Compton scattering between target photons from the optical companion and disc, and relativistic electrons in the pulsar wind. We present results of a full anisotropic inverse Compton scattering model of the system, taking into account the IR excess from the extended circumstellar disc around the optical companion.
Aims. With our time-dependent model atmosphere code PHOENIX, our goal is to simulate light curves and spectra of hydrodynamical models of all types of supernovae. In this work, we simulate near-infrared light curves of SNe Ia and confirm the cause of the secondary maximum. Methods. We apply a simple energy solver to compute the evolution of an SN Ia envelope during the free expansion phase. Included in the solver are energy changes due to expansion, the energy deposition of {\gamma}-rays and interaction of radiation with the material. Results. We computed theoretical light curves of several SN Ia hydrodynamical models in the I, J, H, and K bands and compared them to the observed SN Ia light curves of SN 1999ee and SN 2002bo. By changing a line scattering parameter in time, we obtained quite reasonable fits to the observed near-infrared light curves. This is a strong hint that detailed NLTE effects in IR lines have to be modeled, which will be a future focus of our work. Conclusions. We found that IR line scattering is very important for the near-infrared SN Ia light curve modeling. In addition, the recombination of Fe III to Fe II and of Co III to Co II is responsible for the secondary maximum in the near-infrared bands. For future work the consideration of NLTE for all lines (including the IR subordinate lines) will be crucial.
We present optical spectroscopy of MWC 656 and MWC 148, the proposed optical counterparts of the gamma-ray sources AGL J2241+4454 and HESS J0632+0 57, respectively. The main parameters of the Halpha emission line (EW, FWHM and centroid velocity) in these stars are modulated on the proposed orbital periods of 60.37 and 321 days, respectively. These modulations are likely produced by the resonant interaction of the Be discs with compact stars in eccentric orbits. We also present radial velocity curves of the optical stars folded on the above periods and obtain the first orbital elements of the two gamma-ray sources thus confirming their binary nature. Our orbital solution support eccentricities e~0.4 and 0.83+-0.08 for MWC 656 and MWC 148, respectively. Further, our orbital elements imply that the X-ray outbursts in HESS J0632+057/MWC 148 are delayed ~0.3 orbital phases after periastron passage, similarly to the case of LS I +61 303. In addition, the optical photometric light curve maxima in AGL J2241+4454/MWC 656 occur ~0.25 phases passed periastron, similar to what is seen in LS I +61 303. We also find that the orbital eccentricity is correlated with orbital period for the known gamma-ray binaries. This is explained by the fact that small stellar separations are required for the efficient triggering of VHE radiation. Another correlation between the EW of Halpha and orbital period is also observed, similarly to the case of Be/X-ray binaries. These correlations are useful to provide estimates of the key orbital parameters Porb and e from the Halpha line in future Be gamma-ray binary candidates.
The importance of nuclear reactions in low-density astrophysical plasmas with ion temperatures $T \geq 10^{10}$ K has been recognized for more than thirty years. However, the lack of comprehensive data banks of relevant nuclear reactions and the limited computational power have not previously allowed detailed theoretical studies. Recent developments in these areas make it timely to conduct comprehensive studies on the nuclear properties of very hot plasmas formed around compact relativistic objects such as black holes and neutron stars. Such studies are of great interest in the context of scientific programs of future low-energy cosmic $\gamma$-ray spectrometry. In this work, using the publicly available code TALYS, we have built a large nuclear network relevant for temperatures exceeding $10^{10}$ K. We have studied the evolution of the chemical composition and accompanying prompt gamma-ray emission of such high temperature plasmas. We present the results on the abundances of light elements D, T, $^3$He, $^4$He, $^{6}$Li, $^{7}$Li $^{9}$Be, $^{10}$B, $^{11}$B, and briefly discuss their implications on the astrophysical abundances of these elements.
We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar GX 1+4 and a possible interpretation within a model of quasi-spherical accretion. New period measurements from BeppoSAX, INTEGRAL and Fermi observations have been combined with previously published data from four decades of observations. During the 1970s GX 1+4 was spinning up with the fastest rate among the known X-ray pulsars at the time. In the mid 1980s it underwent a change during a period of low X-ray ux and started to spin down with a rate similar in magnitude to the previous spin up rate. The spin period has changed from ~110 s to ~160 s within the last three decades. Our results demonstrate that the overall spin down trend continues and is stronger than ever. We compare the observations with predictions from a model assuming quasi-spherical accretion from the slow wind of the M giant companion.
The cosmic star formation rate (CSFR), is an important clue to investigate the history of the assembly and evolution of galaxies. Here, we develop a method to study the CSFR from a purely theoretical point of view. Starting from detailed models of chemical evolution, we obtain the histories of star formation of galaxies of different morphological types. These histories are then used to determine the luminosity functions of the same galaxies by means of a spectro-photometric code. We obtain the CSFR under different hypothesis. First, we study the hypothesis of a pure luminosity evolution scenario, in which all galaxies are supposed to form at the same redshift and then evolve only in luminosity. Then we consider scenarios in which the number density or the slope of the LFs are assumed to vary with redshift. After comparison with available data we conclude that a pure luminosity evolution does not provide a good fit to the data, especially at very high redshift, although many uncertainties are still present in the data. On the other hand, a variation in the number density of ellipticals and spirals as a function of redshift can provide a better fit to the observed CSFR. We also explore cases of variable slope of the LFs with redshift and variations of number density and slope at the same time. We cannot find any of those cases which can improve the fit to the data respect to the solely number density variation. Finally, we compute the evolution of the average cosmic metallicity in galaxies with redshift.
We provide measurements of the Ba isotopic fractions for five metal-poor stars derived with an LTE analysis using 1D model stellar atmospheres. We use high resolution (R\equiv{\lambda}/\Delta{\lambda}=90000-95000), very high signal-to-noise (S/N>500) spectra to determine the fraction of odd Ba isotopes (fodd) by measuring subtle asymmetries in the profile of the Ba ii line at 4554 {\AA}. We also use two different macroturbulent broadening techniques, Gaussian and radial-tangential, to model the Fe lines of each star, and propagate each technique to model macroturbulent broadening in the Ba 4554 {\AA} line. We conduct a 1D non-LTE (NLTE) treatment of the Fe lines in the red giant HD122563 and the subgiant HD140283 in an attempt to improve the fitting. We determine [Ba/Eu] ratios for the two giants in our study, HD122563 and HD88609, which can also be used to determine the relative contribution of the s- and r-processes to heavy-element nucleosynthesis, for comparison with fodd. We find fodd for HD122563, HD88609 and HD84937, BD+26\circ3578 and BD-04\circ3208 to be -0.12\pm0.07, -0.02\pm0.09, and -0.05\pm0.11, 0.08\pm0.08 and 0.18\pm0.08 respectively. This means that all stars examined here show isotopic fractions more compatible with an s-process dominated composition. The [Ba/Eu] ratios in HD122563 and HD88609 are found to be -0.20\pm0.15 and -0.47\pm0.15 respectively, which indicate instead an r-process signature. We report a better statistical fit to the majority of Fe profiles in each star when employing a radial-tangential broadening technique during our 1D LTE investigation. We have shown that, from a statistical point of view, one must consider using a radial-tangential broadening technique rather than a Gaussian one to model Fe line macroturbulences when working in 1D. No improvement to Fe line fitting is seen when employing a NLTE treatment.
We studied the Na-O anti-correlation from moderately high resolution spectra for 35 stars on the blue HB (BHB), one RR Lyrae, and 55 stars are on the red HB (RHB) of NGC1851. We also derived abundances for He and N in BHB stars, and Ba and upper limits for N in RHB stars. The RHB stars clearly separate into two groups: the vast majority are O-rich and Na-poor, while about 10-15% are Na-rich and moderately O-poor. Most Na-rich RHB stars are also Ba-rich and there is an overall correlation between Na and Ba abundances within the RHB. The group of Ba-rich RHB stars resides on the warmer edge and includes ~10% of the RHB stars. We propose that they are the descendant of the stars on the RGB sequence with very red v-y colour. This sequence is known also to consist of Ba and perhaps CNO-rich stars. However, the upper limit we obtain for N ([N/Fe]<1.55) for one of the Ba-rich stars coupled with the low C-abundances for RGB Ba-rich stars from the literature suggests that the total CNO might not be particularly high ([(C+N+O)/Fe]<=0.15). The other Na-rich RHB stars are also at the warm edge of the RHB and the only RR Lyrae is Na-rich and moderately O-poor. We also find a Na-O anticorrelation among BHB stars, partially overlapping that found among RHB stars, though generally BHB stars are more Na-rich and O-poor. However, there is no clear correlation between temperature and Na and O abundances within the BHB. The average He abundance in BHB stars is Y=0.29+/-0.05. N abundances are quite uniform at [N/Fe]=1.16+/-0.14 among BHB stars, with a small trend with temperature. This value is consistent with normal CNO abundance and excludes that BHB stars are very CNO-rich: this leaves an age spread of ~1.5 Gyr as the only viable explanation for the split of the SGB. [Abridged]
We investigate the structure of the Milky Way by determining how features in a spatial map correspond to CO features in a velocity map. We examine structures including logarithmic spiral arms, a ring and a bar. We explore the available parameter space, including the pitch angle of the spiral arms, radius of a ring, and rotation curve. We show that surprisingly, a spiral arm provides a better fit to the observed molecular ring than a true ring feature. This is because both a spiral arm, and the observed feature known as the molecular ring, are curved in velocity longitude space. We find that much of the CO emission in the velocity longitude map can be fitted by a nearly symmetric 2 armed spiral pattern. One of the arms corresponds to the molecular ring, whilst the opposite arm naturally reproduces the Perseus arm. Multiple arms also contribute to further emission in the vicinity of the molecular ring and match other observed spiral arms. Whether the Galactic structure consists primarily of two, or several spiral arms, the presence of 2 symmetric logarithmic spirals, which begin in the vicinity of the ends of the bar, suggest a spiral density wave associated with the bar.
Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS16293-2422 and learn more about the water formation mechanisms by the determination of the HDO/H2O abundance ratio. A spectral survey of the source IRAS16293-2422 has been carried out in the framework of the CHESS Herschel Key program with the HIFI instrument, allowing the detection of numerous HDO lines. Other transitions were previously observed with ground-based telescopes in the framework of TIMASSS. The spherical Monte Carlo radiative transfer code RATRAN has been used to reproduce the observed line profiles of HDO assuming an abundance jump, corresponding to the sublimation of the molecules trapped on the icy grain mantles in the hot corino. To determine the H2O abundance throughout the envelope, a similar study has been applied to the H2-18O observed lines, as the H2O main isotope lines are contaminated by the outflows. We derive an inner HDO abundance of 1.7e-7 and an outer HDO abundance of 8e-11. To reproduce the HDO absorption lines, it is necessary to add an absorbing layer in front of the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. The HDO/H2O ratio is ~1.4-5.8% in the hot corino whereas it is ~0.2-2.2% in the outer envelope. It is estimated at ~4.8% in the added absorbing layer. Although it is clearly higher than the cosmic D/H abundance, the HDO/H2O ratio remains lower than the D/H ratio derived for other deuterated molecules observed in the same source. The similar ratios derived in the hot corino and in the added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and methanol which formed later once the CO molecules have depleted on the grains.
SHARDS, an ESO/GTC Large Program, is an ultra-deep (26.5 mag) spectro-photometric survey with GTC/OSIRIS designed to select and study massive passively evolving galaxies at z=1.0-2.3 in the GOODS-N field using a set of 24 medium-band filters (FWHM~17 nm) covering the 500-950 nm spectral range. Our observing strategy has been planned to detect, for z>1 sources, the prominent Mg absorption feature (at rest-frame ~280 nm), a distinctive, necessary, and sufficient feature of evolved stellar populations (older than 0.5 Gyr). These observations are being used to: (1) derive for the first time an unbiased sample of high-z quiescent galaxies, which extends to fainter magnitudes the samples selected with color techniques and spectroscopic surveys; (2) derive accurate ages and stellar masses based on robust measurements of spectral features such as the Mg(UV) or D(4000) indices; (3) measure their redshift with an accuracy Delta(z)/(1+z)<0.02; and (4) study emission-line galaxies (starbursts and AGN) up to very high redshifts. The well-sampled optical SEDs provided by SHARDS for all sources in the GOODS-N field are a valuable complement for current and future surveys carried out with other telescopes (e.g., Spitzer, HST, and Herschel).
The dust in protoplanetary disks is influenced by a lot of different processes. Besides others, heating processes are the most important ones: they change not only the physical and chemical properties of dust particles, but also their emission spectra. In order to compare observed infrared spectra of young stellar systems with laboratory data of hot (up to 700{\deg}C) circumstellar dust analogues, we investigate materials, which are important constituents of dust in protoplanetary disks. We calculated the optical constants by means of a simple Lorentzian oscillator fit and apply them to simulations of small-particle emission spectra in order to compare our results with real astronomical spectra of AGB-stars and protoplanetary disks.
Statistical measures of galaxy clusters are sensitive to neutrino masses in the sub-eV range. We explore the possibility of using cluster number counts from the ongoing PLANCK/SZ and future cosmic-variance-limited surveys to constrain neutrino masses from CMB data alone. The precision with which the total neutrino mass can be determined from SZ number counts is limited mostly by uncertainties in the clusters mass function; these are explicitly accounted for in our forecast. We find that projected results from the PLANCK/SZ survey can be used to determine the total neutrino mass with a ($1\sigma$) uncertainty of 0.08-0.11 eV, assuming fiducial neutrino mass in the range 0.1-0.3 eV, if the survey detection limit is set at the $5\sigma$ significance level. This improves on the limits expected from PLANCK/CMB lensing measurements by a factor 1.5-2. With a cosmic-variance-limited SZ survey we obtain ($1\sigma$) uncertainty of 0.05-0.07 eV. Combined PLANCK and the X-ray RASS cluster catalogs could constrain $M_{\nu}$ at the 0.04 and 0.07, assuming fiducial neutrino masses 0.1 and 0.3 eV, respectively. A few percent uncertainty in the mass function parameters could result in a factor of up to $\sim 50%$ degradation of our PLANCK, CVL and PLANCK+RASS forecasts. The latter two, which significantly benefit from cluster surveys, are more prone to mass function uncertainty. This degradation becomes less significant for neutrino masses 0.3 eV or higher. This highlights the relevance of mass function uncertainties for cosmological parameter estimation. Our analysis shows that if the (total) neutrino mass is close to the lower limits deduced from neutrino oscillation experiments, cluster number counts provide a viable complimentary cosmological probe to CMB lensing constraints on $M_{\nu}$.
Planck's law describes the radiation of black bodies. The study of its properties is of special interest, as black bodies are a good description for the behavior of many phenomena. In this work a new mathematical study of Planck's law is performed and new properties of this old acquaintance are obtained. As a result, the exact form for the locus in a color-color diagrams has been deduced, and an analytical formula to determine with precision the black body temperature of an object from any pair of measurements has been developed. Thus, using two images of the same field obtained with different filters, one can compute a fast estimation of black body temperatures for every pixel in the image, that is, a new image of the black body temperatures for all the objects in the field. Once these temperatures are obtained, the method allows, as a consequence, a quick estimation of their emission in other frequencies, assuming a black body behavior. These results provide new tools for data analysis.
This document describes a FITS convention developed by the IRAF Group (D. Tody, R. Seaman, and N. Zarate) at the National Optical Astronomical Observatory (NOAO). This convention is implemented by the fgread/fgwrite tasks in the IRAF fitsutil package. It was first used in May 1999 to encapsulate preview PNG-format graphics files into FITS files in the NOAO High Performance Pipeline System. A FITS extension of type 'FOREIGN' provides a mechanism for storing an arbitrary file or tree of files in FITS, allowing it to be restored to disk at a later time.
Astrophysical plasmas in accretion discs are usually treated in the framework of fluid or MHD approaches but there are some situations where these treatments become inadequate and one needs to revert to the more fundamental underlying kinetic theory. This occurs when the plasma becomes effectively collisionless or weakly-collisional such as, for example, in radiatively inefficient accretion flows onto black holes. In this paper, we lay down the basics of kinetic theory in these contexts. In particular, we formulate the kinetic theory for quasi-stationary collisionless accretion disc plasmas in the framework of a Vlasov-Maxwell description, taking the plasma to be non-relativistic, axisymmetric, gravitationally-bound and subject to electromagnetic fields. Quasi-stationary solutions for the kinetic distribution functions are constructed which are shown to admit temperature anisotropies. The physical implications of the theory are then investigated and the equations of state and angular momentum conservation law are discussed. Analysis of the Ampere equation reveals the existence of a quasi-stationary kinetic dynamo which gives rise to self-generation of poloidal and azimuthal magnetic fields and operates even in the absence of turbulence and/or instability phenomena.
The physical mechanism responsible for driving accretion flows in astrophysical accretion disks is commonly thought to be related to the development of plasma instabilities and turbulence. A key question is therefore the determination of consistent equilibrium configurations for accretion-disk plasmas and investigation of their stability properties. In the case of collisionless plasmas kinetic theory provides the appropriate theoretical framework. This paper presents a kinetic description of low-frequency and long-wavelength axisymmetric electromagnetic perturbations in non-relativistic, strongly-magnetized and gravitationally-bound axisymmetric accretion-disk plasmas in the collisionless regime. The analysis, carried out within the framework of the Vlasov-Maxwell description, relies on stationary kinetic solutions of the Vlasov equation which allow for the simultaneous treatment of non-uniform fluid fields, stationary accretion flows and temperature anisotropies. It is demonstrated that, for such solutions, no axisymmetric unstable perturbations can exist occurring on characteristic time and space scales which are long compared with the Larmor gyration time and radius. Hence, these stationary configurations are actually stable against axisymmetric kinetic instabilities of this type. As a fundamental consequence, this rules out the possibility of having the axisymmetric magneto-rotational or thermal instabilities to arise in these systems.
Determining the mass of stars is crucial both to improving stellar evolution
theory and to characterising exoplanetary systems. Asteroseismology offers a
promising way to estimate stellar mean density. When combined with accurate
radii determinations, such as is expected from GAIA, this yields accurate
stellar masses. The main difficulty is finding the best way to extract the mean
density from a set of observed frequencies.
We seek to establish a new method for estimating stellar mean density, which
combines the simplicity of a scaling law while providing the accuracy of an
inversion technique.
We provide a framework in which to construct and evaluate kernel-based linear
inversions which yield directly the mean density of a star. We then describe
three different inversion techniques (SOLA and two scaling laws) and apply them
to the sun, several test cases and three stars.
The SOLA approach and the scaling law based on the surface correcting
technique described by Kjeldsen et al. (2008) yield comparable results which
can reach an accuracy of 0.5 % and are better than scaling the large frequency
separation. The reason for this is that the averaging kernels from the two
first methods are comparable in quality and are better than what is obtained
with the large frequency separation. It is also shown that scaling the large
frequency separation is more sensitive to near-surface effects, but is much
less affected by an incorrect mode identification. As a result, one can
identify pulsation modes by looking for an l and n assignment which provides
the best agreement between the results from the large frequency separation and
those from one of the two other methods. Non-linear effects are also discussed
as is the effects of mixed modes. In particular, it is shown that mixed modes
bring little improvement as a result of their poorly adapted kernels.
We present a new model for the behavior of scattered time-dependent, asymmetric near-UV emission from the nearby ejecta of {\eta} Car. Using a 3-D hydrodynamical simulation of {\eta} Car's binary colliding winds, we show that the 3-D binary orientation derived by Madura et al. (2012) is capable of explaining the asymmetric near-UV variability observed in the Hubble Space Telescope Advanced Camera for Surveys/High Resolution Camera (HST ACS/HRC) F220W images of Smith et al. (2004b). Models assuming a binary orientation with i ~ 130 to 145 degrees, {\omega} ~ 230 to 315 degrees, PAz ~ 302 to 327 degrees are consistent with the observed F220W near-UV images. We find that the hot binary companion does not significantly contribute to the near-UV excess observed in the F220W images. Rather, we suggest that a bore-hole effect and the reduction of Fe II optical depths inside the wind-wind collision cavity carved in the extended photosphere of the primary star lead to the time-dependent directional illumination of circum-binary material as the companion moves about in its highly elliptical orbit.
The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.
The bright Seyfert 1 galaxy Mrk 509 was monitored by XMM-Newton and other satellites in 2009 to constrain the location of the outflow. We have studied the response of the photoionised gas to changes in the ionising flux produced by the central regions. We used the 5 discrete ionisation components A-E detected in the time-averaged spectrum taken with the RGS. Using the ratio of fluxed EPIC and RGS spectra, we put tight constraints on the variability of the absorbers. Monitoring with the Swift satellite started 6 weeks before the XMM-Newton observations, allowing to use the ionising flux history and to develop a model for the time-dependent photoionisation. Components A and B are too weak for variability studies, but the distance for component A is known from optical imaging of the [O III] line to be ~3 kpc. During the 5 weeks of the XMM-Newton observations we found no evidence of changes in the 3 X-ray dominant ionisation components C-E, despite a huge soft X-ray intensity increase of 60% in the middle of our campaign. This excludes high-density gas close to the black hole. Instead, using our time-dependent modelling, we find low density and derive firm lower limits to the distance of these components. Component D shows evidence for variability on longer time scales, yielding an upper limit to the distance. For component E we derive an upper limit to the distance based on the argument that the thickness of the absorbing layer must be less than its distance to the black hole. Combining these results, at the 90% confidence level, component C has a distance of >70 pc, component D between 5-33 pc, and component E >5 pc but smaller than 21-400 pc, depending upon modelling details. These results are consistent with the upper limits from the HST/COS observations of our campaign and point to an origin of the dominant, slow (v<1000 km/s) outflow components in the NLR or torus-region of Mrk 509.
Astronomy, as many other scientific disciplines, is facing a true data deluge which is bound to change both the praxis and the methodology of every day research work. The emerging field of astroinformatics, while on the one end appears crucial to face the technological challenges, on the other is opening new exciting perspectives for new astronomical discoveries through the implementation of advanced data mining procedures. The complexity of astronomical data and the variety of scientific problems, however, call for innovative algorithms and methods as well as for an extreme usage of ICT technologies.
We analyze the deviations of transit times from a linear ephemeris for the Kepler Objects of Interest (KOI) through Quarter six (Q6) of science data. We conduct two statistical tests for all KOIs and a related statistical test for all pairs of KOIs in multi-transiting systems. These tests identify several systems which show potentially interesting transit timing variations (TTVs). Strong TTV systems have been valuable for the confirmation of planets and their mass measurements. Many of the systems identified in this study should prove fruitful for detailed TTV studies.
Harmonic analysis is a tool to infer cosmic topology from the measured astrophysical cosmic microwave background CMB radiation. For overall positive curvature, Platonic spherical manifolds are candidates for this analysis. We combine the specific point symmetry of the Platonic manifolds with their deck transformations. This analysis in topology leads from manifolds to orbifolds. We discuss the deck transformations of the orbifolds and give basis functions for the harmonic analysis as linear combinations of Wigner polynomials on the 3-sphere. They provide new tools for detecting cosmic topology from the CMB radiation.
Transit timing variations provide a powerful tool for confirming and characterizing transiting planets, as well as detecting non-transiting planets. We report the results an updated TTV analysis for 822 planet candidates (Borucki et al. 2011; Batalha et al. 2012) based on transit times measured during the first seventeen months of Kepler observations (Rowe et al 2012). We present 35 TTV candidates (4.1% of suitable data sets) based on long-term trends and 153 mostly weaker TTV candidates (18% of suitable data sets) based on excess scatter of TTV measurements about a linear ephemeris. We anticipate that several of these planet candidates could be confirmed and perhaps characterized with more detailed TTV analyses using publicly available Kepler observations. For many others, Kepler has observed a long-term TTV trend, but an extended Kepler mission will be required to characterize the system via TTVs. We find that the occurence rate of planet candidates that show TTVs is significantly increased (~60%-76%) for planet candidates transiting stars with multiple transiting planet candidate when compared to planet candidates transiting stars with a single transiting planet candidate.
We present X-ray observations of the new transient magnetar Swift J1834.9-0846, discovered with Swift BAT on 2011 August 7. The data were obtained with Swift, RXTE, CXO, and XMM-Newton both before and after the outburst. Timing analysis reveals singe peak pulsations with a period of 2.4823 s and an unusually high pulsed fraction, 85+/-10%. Using the RXTE and CXO data, we estimated the period derivative, dot{P}=8\times 10^{-12} s/s, and confirmed the high magnetic field of the source, B=1.4\times 10^{14} G. The decay of the persistent X-ray flux, spanning 48 days, is consistent with a power law, t^{-0.5}. In the CXO/ACIS image, we find that the highly absorbed point source is surrounded by extended emission, which most likely is a dust scattering halo. Swift J1834.9-0846 is located near the center of the radio supernova remnant W41 and TeV source HESS J1834-087. An association with W41 would imply a source distance of about 4 kpc; however, any relation to the HESS source remains unclear, given the presence of several other candidate counterparts for the latter source in the field. Our search for an IR counterpart of Swift J1834.9-0846 revealed no source down to K_s=19.5 within the 0.6' CXO error circle.
Weakly Interacting Massive Particles (WIMPs) can be captured by the Sun and the Earth, sink to their cores, annihilate and produce neutrinos that can be searched for with neutrino telescopes. The calculation of the capture rates of WIMPs in the Sun and especially the Earth are affected by large uncertainties coming mainly from effects of the planets in the solar system, reducing the capture rates by up to an order of magnitude (or even more in some cases). We show that the WIMPs captured by weak scatterings in the Sun also constitute an important bound WIMP population in the solar system. Taking this population and its interplay with the population bound through gravitational diffusion into account cancel the planetary effects on the capture rates, and the capture essentially proceeds as if the Sun and the Earth were free in the galactic halo. The neutrino signals from the Sun and the Earth are thus significantly higher than claimed in the scenarios with reduced capture rates.
Using the results of a previous X-ray photo-ionization modelling of blue-shifted Fe K absorption lines on a sample of 42 local radio-quiet AGNs observed with XMM-Newton, in this letter we estimate the location and energetics of the associated ultra-fast outflows (UFOs). Due to significant uncertainties, we are essentially able to place only lower/upper limits. On average, their location is in the interval ~0.0003-0.03pc (~10^2-10^4 r_s) from the central black hole, consistent with what is expected for accretion disk winds/outflows. The mass outflow rates are constrained between ~0.01-1 M_{\odot} yr^{-1}, corresponding to >5-10% of the accretion rates. The average lower-upper limits on the mechanical power are log\dot{E}_K~42.6-44.6 erg s^{-1}. However, the minimum possible value of the ratio between the mechanical power and bolometric luminosity is constrained to be comparable or higher than the minimum required by simulations of feedback induced by winds/outflows. Therefore, this work demonstrates that UFOs are indeed capable to provide a significant contribution to the AGN cosmological feedback, in agreement with theoretical expectations and the recent observation of interactions between AGN outflows and the interstellar medium in several Seyferts galaxies.
Two new families of self-consistent axisymmetric truncated equilibrium models for the description of quasi-relaxed rotating stellar systems are presented. The first extends the spherical King models to the case of solid-body rotation. The second is characterized by differential rotation, designed to be rigid in the central regions and to vanish in the outer parts, where the energy truncation becomes effective. The models are constructed by solving the nonlinear Poisson equation for the self-consistent mean-field potential. For rigidly rotating configurations, the solutions are obtained by an asymptotic expansion on the rotation strength parameter. The differentially rotating models are constructed by means of an iterative approach based on a Legendre series expansion of the density and the potential. The two classes of models exhibit complementary properties. The rigidly rotating configurations are flattened toward the equatorial plane, with deviations from spherical symmetry that increase with the distance from the center. For models of the second family, the deviations from spherical symmetry are strongest in the central region, whereas the outer parts tend to be quasi-spherical. The relevant parameter spaces are explored and the intrinsic and projected structural properties are described. Special attention is given to the effect of different options for the truncation of the distribution function in phase space. Models in the moderate rotation regime are best suited to applications to globular clusters. For general interest in stellar dynamics, at high values of the rotation strength the differentially rotating models exhibit a toroidal core embedded in a quasi-spherical configuration. Physically simple analytical models of the kind presented here provide insights into dynamical mechanisms and may be a basis for more realistic investigations with the help of N-body simulations.
We study magnetic dipole moments of right-handed neutrinos in a keV neutrino dark matter model. This model is a simple extension of the standard model with only right-handed neutrinos and a pair of charged particles added. One of the right-handed neutrinos is the candidate of dark matter with a keV mass. Some bounds on the dark matter magnetic dipole moment and model parameters are obtained from cosmological observations.
The ultimate climate emergency is a "runaway greenhouse": a hot and water vapour rich atmosphere limits the emission of thermal radiation to space, causing runaway warming. Warming ceases only once the surface reaches ~1400K and emits radiation in the near-infrared, where water is not a good greenhouse gas. This would evaporate the entire ocean and exterminate all planetary life. Venus experienced a runaway greenhouse in the past, and we expect that Earth will in around 2 billion years as solar luminosity increases. But could we bring on such a catastrophe prematurely, by our current climate-altering activities? Here we review what is known about the runaway greenhouse to answer this question, describing the various limits on outgoing radiation and how climate will evolve between these. The good news is that almost all lines of evidence lead us to believe that is unlikely to be possible, even in principle, to trigger full a runaway greenhouse by addition of non-condensible greenhouse gases such as carbon dioxide to the atmosphere. However, our understanding of the dynamics, thermodynamics, radiative transfer and cloud physics of hot and steamy atmospheres is weak. We cannot therefore completely rule out the possibility that human actions might cause a transition, if not to full runaway, then at least to a much warmer climate state than the present one. High climate sensitivity might provide a warning. If we, or more likely our remote descendants, are threatened with a runaway greenhouse then geoengineering to reflect sunlight might be life's only hope. ...[2 sentences cut to meet arXiv char limit]... The runaway greenhouse also remains relevant in planetary sciences and astrobiology: as extrasolar planets smaller and nearer to their stars are detected, some will be in a runaway greenhouse state.
We have calculated the ground state electronic structure of He under pressure from 0 to 1500 GPa using both all-electron full-potential and pseudopotential methods based on the density functional theory (DFT). We find that throughout this pressure range, pseudopotentials yield essentially the same energy-volume curve for all of bcc, fcc, and hcp configurations as does the full-potential method, a strong indication that pseudopotential approximation works well for He both as the common element in some giant planets and as detrimental impurities in fusion reactor materials. The hcp lattice is always the most stable structure and bcc the least stable one. Since the energy preference of hcp over fcc and bcc is within 0.01 eV below 100 GPa and about 0.1 eV at 1500 GPa, on the same order of the error bar in local or semi-local density approximations in DFT, phase transitions can only be discussed with more precise description of electron correlation in Quantum Monte Carlo or DFT-based GW methods.
Astrophysical black hole candidates are thought to be the Kerr black holes predicted by General Relativity, as these objects cannot be explained otherwise without introducing new physics. However, there is no observational evidence that the space-time around them is really described by the Kerr solution. The Kerr black hole hypothesis can be tested with the already available X-ray data by extending the continuum-fitting method, a technique currently used by astronomers to estimate the spins of stellar-mass black hole candidates. In general, we cannot put a constraint on possible deviations from the Kerr geometry, but only on some combination between these deviations and the spin. The measurement of the radio power of transient jets in black hole binaries can potentially break this degeneracy, thus allowing for testing the Kerr-nature of these objects.
The combination of GR and the Standard Model disagrees with numerous observations on scales from our Solar System up. In the concordance model of cosmology, these contradictions are removed or alleviated by the introduction of three completely independent new components of stress-energy -- the inflaton, dark matter, and dark energy. Each of these in its turn is meant to have (or to currently) dominate the dynamics of the universe. There is still no non-gravitational evidence for any of these dark sectors; nor for the required extensions of the standard model. An alternative is to imagine that GR itself must be modified. Certain coincidences of scale even suggest that one might expect not to have to make three independent. Because they must address the most different types of data, attempts to replace dark matter with modified gravity are the most controversial. A phenomenological model (or family of models), Modified Newtonian Dynamics, has, over the last few years seen several covariant realizations. We discuss a number of challenges that any model that seeks to replace dark matter with modified gravity must face: the loss of Birkhoff's Theorem, and the calculational simplifications it implies; the failure to explain clusters, whether static or interacting, and the consequent need to introduce dark matter of some form, whether hot dark matter neutrinos, or dark fields that arise in new sectors of the modified gravity theory; the intrusion of cosmological expansion into the modified force law, that arises precisely because of the coincidence in scale between the centripetal acceleration at which Newtonian gravity fails in galaxies, and the cosmic acceleration. We conclude with the observation that, although modified gravity may indeed manage to replace dark matter, it is likely to do so by becoming or incorporating, a dark matter theory itself.
In these comments, we clarify the problematic aspects of gravitational interaction in a weak-field limit of Kaluza-Klein models. We explain why some models meet the classical gravitational tests, while others do not. We show that variation of the total volume of the internal spaces generates the fifth force. This is the main reason of the problem. It happens for all considered models (linear with respect to the scalar curvature and non-linear f(R), with toroidal and spherical compactifications). In the case of models with toroidal compactification, we demonstrate how tension (with and without effects of non-linearity) of the gravitating source can eliminate the fifth force, resulting in agreement with the observations. It takes place for latent solitons, black strings and black branes. In the case of spherical compactification, the fifth force is replaced by the Yukawa interaction for models with the stabilized internal space. For large Yukawa masses, the effect of this interaction is negligibly small, and considered models satisfy the gravitational tests at the same level of accuracy as General Relativity.
Links to: arXiv, form interface, find, astro-ph, recent, 1201, contact, help (Access key information)
Magnetic Doppler imaging is currently the most powerful method of interpreting high-resolution spectropolarimetric observations of stars. This technique has revealed the presence of unexpected small-scale magnetic fields on the surfaces of Ap stars. These studies were recently criticisied by Stift et al. (2012), who claimed that magnetic inversions are not robust and are undermined by neglecting a feedback on the Stokes line profiles from the local atmospheric structure in the regions of enhanced metal abundance. We show that Stift et al. misinterpreted published magnetic Doppler imaging results and neglected some of the most fundamental principles behind magnetic mapping. We demonstrate that the variation of atmospheric structure across the surface of a star with chemical spots affects the local continuum intensity but is negligible for the normalised local Stokes profiles. For the disk-integrated spectra of an Ap star with extreme abundance variations, we find that the assumption of a mean model atmosphere leads to moderate errors in Stokes I but is negligible for polarisation spectra. Employing a new magnetic inversion code, which incorporates the horizontal variation of atmospheric structure, we reconstructed new maps of magnetic field and Fe abundance for the Ap star alpha^2 CVn. The resulting distribution of chemical spots changes insignificantly compared to the previous modelling based on a single model atmosphere, while the magnetic field geometry does not change at all. This shows that the assertions by Stift et al. are exaggerated as a consequence of unreasonable assumptions and extrapolations, as well as methodological flaws and inconsistencies of their analysis. Our discussion proves that published magnetic inversions based on a mean stellar atmosphere are robust and reliable, and that the presence of small-scale magnetic field structures on the surfaces of Ap stars is real.
We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL (dNL) AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that dNL AGN induced by the relative motion of SMBH pairs are a generic but short-lived feature of gaseous major mergers. dNL AGN are most likely to be observed in late-stage mergers, during the kpc-scale phase of SMBH inspiral or soon after the SMBH merger. However, even within the kpc-scale phase, only a minority of dNL AGN are directly induced by SMBH motion; their lifetimes are typically a few Myr. Most double peaks arise from gas kinematics near the SMBH, although prior to the SMBH merger up to 80% of all dNL profiles may be influenced by SMBH motion via altered peak ratios or velocity offsets. The total lifetimes of dNL AGN depend strongly on viewing angle and on properties of the merging galaxies. Also, in a typical merger, at least 10-40% of the double peaks induced by SMBH motion have small projected separations, 0.1-1 kpc, such that dual peaks of stellar surface brightness are not easily resolved. Diffuse tidal features can indicate late-stage galaxy mergers, although they do not distinguish SMBH pairs from merged SMBHs. We show that dNL profiles with peak velocity splittings > 500 km s^-1 or with measurable overall velocity shifts are often associated with SMBH pairs. Our results support the notion that selection of dNL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multi-wavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs. (Abridged)
(Abridged) Since its launch in October 2002, the INTEGRAL satellite has revolutionized our knowledge of the hard X-ray sky thanks to its unprecedented imaging capabilities and source detection positional accuracy above 20 keV. Nevertheless, many of the newly-detected sources in the INTEGRAL sky surveys are of unknown nature. The combined use of available information at longer wavelengths (mainly soft X-rays and radio) and of optical spectroscopy on the putative counterparts of these new hard X-ray objects allows us to pinpoint their exact nature. Continuing our long-standing program that has been running since 2004, and using 6 different telescopes of various sizes, we report the classification through optical spectroscopy of 22 more unidentified or poorly studied high-energy sources detected with the IBIS instrument onboard INTEGRAL. We found that 16 of them are active galactic nuclei (AGNs), while the remaining 6 objects are within our Galaxy. Among the identified extragalactic sources, 14 are Type 1 AGNs; of these, 6 lie at redshift larger than 0.5 and one has z = 3.12, which makes it the second farthest object detected in the INTEGRAL surveys up to now. The remaining AGNs are of type 2, and one of them is a pair of interacting Seyfert 2 galaxies. The Galactic objects are identified as two cataclysmic variables, one high-mass X-ray binary, one symbiotic binary and two chromospherically active stars. We thus still find that AGNs are the most abundant population among hard X-ray objects identified through optical spectroscopy. Moreover, we note that the higher sensitivity of the more recent INTEGRAL surveys is now enabling the detection of high-redshift AGNs, thus allowing the exploration of the most distant hard X-ray emitting sources and possibly of the most extreme blazars.
We present and analyse optical and ultra-violet colours for passive and optically-red Coma cluster galaxies for which we have spectroscopic age and element abundance estimates. Our sample of 150 objects covers a wide range in mass, from giant ellipticals to the bright end of the dwarf-galaxy regime. We focus on the colours FUV-i, NUV-i, FUV-NUV, u*-g and g-i. We find that all of these colours are correlated with both luminosity and velocity dispersion at the >5 sigma level, with FUV-i and FUV-NUV becoming bluer with increasing `mass' while the other colours become redder. We perform an empirical analysis to assess what fraction of the variation in each colour can be accounted for by variations in the average stellar populations, as traced by the optical spectra. For u*-g and g-i, most of the observed scatter (~80% after allowing for measurement errors and for systematic errors in u*-g) is attributable to stellar population variations, with colours becoming redder with increasing age and metallicity (Mg/H). The FUV-i colour becomes bluer with increasing age and with increasing Mg/H, favouring a `metal-rich single-star' origin for the UV upturn. However, correlations with the optically-dominant stellar populations account for only about half of the large observed scatter. We propose that the excess scatter in FUV-i may be due to a varying proportion of ancient stars in galaxies with younger average ages. The NUV-i colour is sensitive to age and Mg/H, but exhibits excess scatter that can be attributed to `leakage' of the FUV upturn. Applying a correction based on the FUV-i colour we account for ~80% of the variance in NUV-i, as for the optical colours. The FUV-NUV colour shows strong correlations with age and Mg/H, and little residual scatter. Interpreting this colour is complicated however, since it mixes the effects from the main-sequence turn-off with those from the hot evolved stars.
We present neutrino mass bounds using 900,000 luminous galaxies with photometric redshifts measured from Sloan Digital Sky Survey III Data Release Eight (SDSS DR8). The galaxies have photometric redshifts between $z = 0.45$ and $z = 0.65$, and cover 10,000 square degrees and thus probe a volume of 3$h^{-3}$Gpc$^3$, enabling tight constraints to be derived on the amount of dark matter in the form of massive neutrinos. A new bound on the sum of neutrino masses $\sum m_\nu < 0.26$ eV, at 95% confidence level (CL), is obtained after combining our sample of galaxies, which we call "CMASS", with WMAP 7 year Cosmic Microwave Background (CMB) data and the most recent measurement of the Hubble parameter from the Hubble Space Telescope (HST). This constraint is obtained with a conservative multipole range choice of $30 < \ell < 200$ in order to minimize non-linearities, and a free bias parameter in each of the four redshift bins. We study the impact of assuming this linear galaxy bias model using mock catalogs, and find that this model causes a small ($\sim 1-1.5 \sigma$) bias in $\Omega_{\rm DM} h^2$. For this reason, we also quote neutrino bounds based on a conservative galaxy bias model containing additional, shot noise-like free parameters. In this conservative case, the bounds are significantly weakened, e.g. $\sum m_\nu < 0.36$ eV (95% confidence level) for WMAP+HST+CMASS ($\ell_{\rm max}=200$). We also study the dependence of the neutrino bound on multipole range ($\ell_{\rm max}=150$ vs $\ell_{\rm max}=200$) and on which combination of data sets is included as a prior. The addition of supernova and/or Baryon Acoustic Oscillation data does not significantly improve the neutrino mass bound once the HST prior is included. [abridged]
The interstellar and intra-cluster medium in giant elliptical galaxies and clusters of galaxies is often assumed to be in hydrostatic equilibrium. Numerical simulations, however, show that about 5-30% of the pressure in a cluster is provided by turbulence induced by, for example, the central AGN and merger activity. We aim to put constraints on the turbulent velocities and turbulent pressure in the ICM of the giant elliptical galaxies NGC 5044 and NGC 5813 using XMM-Newton RGS observations. The magnitude of the turbulence is estimated using the Fe XVII lines at 15.01 A, 17.05 A, and 17.10 A in the RGS spectra. At low turbulent velocities, the gas becomes optically thick in the 15.01 A line due to resonant scattering, while the 17 A lines remain optically thin. By comparing the (I(17.05)+I(17.10))/I(15.01) line ratio from RGS with simulated line ratios for different Mach numbers, the level of turbulence is constrained. The measurement is limited by systematic uncertainties in the atomic data, which are at the 20-30% level. We find that the line ratio in NGC 5813 is significantly higher than in NGC 5044. This difference can be explained by a higher level of turbulence in NGC 5044. The high turbulent velocities and the fraction of the turbulent pressure support of >40% in NGC 5044, assuming isotropic turbulence, confirm that it is a highly disturbed system, probably due to an off-axis merger. The turbulent pressure support in NGC 5813 is more modest at 15-45%. The (I(17.05)+I(17.10))/I(15.01) line ratio in an optically thin plasma, calculated using AtomDB v2.0.1, is 2 sigma above the ratio measured in NGC 5044, which cannot be explained by resonant scattering. This shows that the discrepancies between theoretical, laboratory, and astrophysical data on Fe XVII lines need to be reduced to improve the accuracy of the determination of turbulent velocities using resonant scattering.
We present resolved Herschel images of a circumbinary debris disk in the 99 Herculis system. The primary is a late F-type star. The binary orbit is well characterised and we conclude that the disk is misaligned with the binary plane. Two different models can explain the observed structure. The first model is a ring of polar orbits that move in a plane perpendicular to the binary pericenter direction. We favour this interpretation because it includes the effect of secular perturbations and the disk can survive for Gyr timescales. The second model is a misaligned ring. Because there is an ambiguity in the orientation of the ring, which could be reflected in the sky plane, this ring either has near-polar orbits similar to the first model, or has a 30 degree misalignment. The misaligned ring, interpreted as the result of a recent collision, is shown to be implausible from constraints on the collisional and dynamical evolution. Because disk+star systems with separations similar to 99 Herculis should form coplanar, possible formation scenarios involve either a close stellar encounter or binary exchange in the presence of circumstellar and/or circumbinary disks. Discovery and characterisation of systems like 99 Herculis will help understand processes that result in planetary system misalignment around both single and multiple stars.
To investigate the universality hypothesis of the initial mass function in the substellar regime, the population of the rho Ophiuchi molecular cloud is analysed by including a new sample of low-mass spectroscopically confirmed members. We have conducted a large spectroscopic follow-up of young substellar candidates uncovered in our previous photometric survey. It is shown that the study of the substellar population of the rho Ophiuchi molecular cloud is hampered only by the high extinction in the cluster ruling out an apparent paucity of brown dwarfs. The spectral types and extinction are derived for a newly found population of substellar objects, and its masses estimated by comparison to evolutionary models. A thoroughly literature search is conducted to provide an up-to-date census of the cluster, which is then used to derive the luminosity and mass functions, as well as the ratio of brown dwarfs to stars in the cluster. These results are compared to other young clusters. The discovery of 16 new members of rho Ophiuchi, 13 of them in the substellar regime, reveals the low mass end of its population and shows the success of our photometric candidate selection with the WIRCam survey. The study of the brown dwarf population of the cluster reveals a high disk fraction of 76 (+5-8)%. Taking into account the characteristic peak mass of the derived mass function and the ratio of brown dwarfs to stars, we conclude that the mass function of rho Ophiuchi is similar to other nearby young clusters.
In this work we investigate the link between galaxy velocity dispersion, mass and other properties (color, morphology) with the properties of dark matter halos by comparing the clustering of galaxies at both fixed mass and velocity dispersion. We use the Sloan Digital Sky Survey to define a volume limited sample of massive galaxies complete in both stellar mass (>6e10 Msun) and velocity dispersion (>75 km/s). Using this sample we show that at fixed velocity dispersion there is no dependence of the clustering amplitude on stellar or dynamical mass. Conversely when stellar or dynamical mass are fixed there is a clear dependence of the clustering amplitude on velocity dispersion with higher dispersion galaxies showing a higher clustering amplitude. We also show that whilst when stellar or dynamical mass are fixed there remains a dependence of clustering amplitude on morphology, there is no such dependency when dispersion is fixed. However, we do see a dependence of the clustering amplitude on color when both mass and dispersion are fixed. Despite this, even when we restrict our samples to only elliptical or red galaxies the relationship between dispersion and clustering amplitude at fixed mass remains. It seems likely that the residual correlation with color is driven by satellite galaxies in massive halos being redder at fixed dispersion. The lack of a similar morphology dependence implies that the mechanism turning satellites red is not changing their morphology. Our central result is that velocity dispersion is more closely related to the clustering amplitude of galaxies than either stellar or dynamical mass. This implies that velocity dispersion is more tightly correlated with the halo properties that determine clustering, either halo mass or age, and supports the notion that the star formation history of a galaxy is more closely related to its halo properties than its overall mass.
The clumping of massive star winds is an established paradigm confirmed by multiple lines of evidence and supported by stellar wind theory. The purpose of this paper is to bridge the gap between detailed models of inhomogeneous stellar winds in single stars and the phenomenological description of donor winds in supergiant high-mass X-ray binaries (HMXBs). We use results from time-dependent hydrodynamical models of the instability in the line-driven wind of a massive supergiant star to derive the time-dependent accretion rate onto a compact object in the Bondi-Hoyle-Lyttleton approximation. The strong density and velocity fluctuations in the wind result in strong variability of the synthetic X-ray light curves. The model predicts a large scale X-ray variability, up to eight orders of magnitude, on relatively short timescales. The apparent lack of evidence for such strong variability in the observed HMXBs indicates that the details of accretion process act to reduce the variability due to the stellar wind velocity and density jumps. We, also, study the absorption of X-rays in the clumped stellar wind by means of a 2-D stochastic wind model and find that absorption of X-rays changes strongly at different orbital phases. Furthermore, we address the photoionization in the clumped wind, and show that the degree of ionization is affected by the wind clumping. A correction factor for the photoionization parameter is derived. It is shown that the photoionization parameter is reduced by a factor Xi compared to the smooth wind models with the same mass-loss rate, where Xi is the wind inhomogeneity parameter. We conclude that wind clumping must also be taken into account when comparing the observed and model spectra of the photoionized stellar wind.
We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the sub-mm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into sub-mm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Sub-mm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.
We report transient quasi-periodic oscillations (QPOs) on minute timescales in relativistic, radiative models of the galactic center source Sgr A*. The QPOs result from nonaxisymmetric $m=1$ structure in the accretion flow excited by MHD turbulence. Near-infrared (NIR) and X-ray power spectra show significant peaks at frequencies comparable to the orbital frequency at the innermost stable circular orbit (ISCO) $f_o$. The excess power is associated with inward propagating magnetic filaments inside the ISCO. The amplitudes of the QPOs are sensitive to the electron distribution function. We argue that transient QPOs appear at a range of frequencies in the neighborhood of $f_o$ and that the power spectra, averaged over long times, likely show a broad bump near $f_o$ rather than distinct, narrow QPO features.
Giant molecular clouds contain supersonic turbulence and simulations of MHD turbulence show that these supersonic motions decay in roughly a crossing time, which is less than the estimated lifetimes of molecular clouds. Such a situation requires a significant release of energy. We run models of C-type shocks propagating into gas with densities around 10^3 cm^(-3) at velocities of a few km / s, appropriate for the ambient conditions inside of a molecular cloud, to determine which species and transitions dominate the cooling and radiative energy release associated with shock cooling of turbulent molecular clouds. We find that these shocks dissipate their energy primarily through CO rotational transitions and by compressing pre-existing magnetic fields. We present model spectra for these shocks and by combining these models with estimates for the rate of turbulent energy dissipation, we show that shock emission should dominate over emission from unshocked gas for mid to high rotational transitions (J >5) of CO. We also find that the turbulent energy dissipation rate is roughly equivalent to the cosmic ray heating rate and that the ambipolar diffusion heating rate may be significant, especially in shocked gas.
The XMM-Newton Distant Cluster Project is a serendipitous survey for clusters of galaxies at redshifts z>=0.8 based on deep archival XMM-Newton observations. ... Low-significance candidate high-z clusters are followed up with the seven-channel imager GROND (Gamma-Ray Burst Optical and Near-Infrared Detector) that is mounted at a 2m-class telescope. ... The test case is XMMU J0338.7+0030, suggested to be at z~1.45+/-0.15 from the analysis of the z-H vs H colour-magnitude diagram obtained from the follow-up imaging. Later VLT-FORS2 spectroscopy enabled us to identify four members, which set this cluster at z=1.097+/-0.002. To reach a better knowledge of its galaxy population, we observed XMMU J0338.7+0030 with GROND for about 6 hr. The publicly available photo-z code le Phare was used. The Ks-band number counts of the non-stellar sources out of the 832 detected down to z'~26 AB-mag in the 3.9x4.3 square arcmin region of XMMU J0338.7+0030 imaged at all GROND bands clearly exceed those computed in deep fields/survey areas at ~20.5 - 22.5 AB-mag. The photo-z's of the three imaged spectroscopic members yield z=1.12+/-0.09. The spatial distribution and the properties of the GROND sources with a photo-z in the range 1.01 - 1.23 confirm the correspondence of the X-ray source with a galaxy over-density at a significance of at least 4.3 sigma. Candidate members that are spectro-photometrically classified as elliptical galaxies define a red locus in the i'-z' vs z' colour-magnitude diagram that is consistent with the red sequence of the cluster RDCS J0910+5422 at z=1.106. XMMU J0338.7+0030 hosts also a population of bluer late-type spirals and irregulars. The starbursts among the photometric members populate both loci, consistently with previous results. The analysis of the available data set indicates that XMMU J0338.7+0030 is a low-mass cluster (M_200 ~ 1E14 M_sun) at z=1.1. (Abridged)
The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) were previously unknown. We present high quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have lightcurves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10,000 asteroids mainly in the magnitude range ~14 to ~20.
Simulations have shown that a diverse range of extrasolar terrestrial planet bulk compositions are likely to exist, based on the observed variations in host star elemental abundances. Based on recent studies, it is expected that a significant proportion of host stars may have Mg/Si ratios below 1. Here we examine this previously neglected group of systems. Planets simulated as forming within these systems are found to be Mg-depleted (compared to the Earth), consisting of silicate species such as pyroxene and various feldspars. Planetary carbon abundances also vary in accordance with the host stars C/O ratio. The predicted abundances are in keeping with observations of polluted white dwarfs, lending validity to this approach. Further studies are required to determine the full planetary impacts of the bulk compositions predicted here.
[abridged] We introduce the Phoenix Project, a set of LCDM simulations of the dark matter component of nine rich galaxy clusters. Each cluster is simulated at least at two different numerical resolutions. For eight of them, the highest resolution corresponds to ~1.3e8 particles within the virial radius, while for one this number is over one billion. Because of their recent assembly, these cluster haloes are significantly less relaxed than galaxy haloes, leading to decreased regularity, increased halo-to-halo variations, and systematic differences in concentration and substructure fraction. All density profiles steepen gradually from the centre outwards, but there is considerable scatter in the dependence of logarithmic slope, gamma on radius. At the innermost convergence radius, r_conv ~3 kpc/h (~ 0.2% of the virial radius) the mean and rms scatter is gamma=1.05+-0.19 for the nine haloes. As for galaxy haloes, there is little indication of an approach to an asymptotic inner power law. For individual clusters, strongly aspherical mass distributions can produce projected surface density variations at given radius spanning up to a factor of three, depending on projection direction. This may in part explain the high apparent concentration of some observed strong-lensing clusters. The shape of the surface density profile, gamma_p(R) depends only weakly on projection direction, however, and is quite well approximated in the inner regions by the NFW formula. Substructure in the Phoenix haloes is slightly more abundant, especially in the inner regions, than in the galaxy haloes of the Aquarius Project. The subhalo mass function is also steeper: dN/dM \propto M^{-1.98} in the range 1e-6<M_sub/M200<1e-3, compared to M^{-1.94} for Aquarius haloes. Resolved subhaloes nevertheless contribute only 11 +-3 of the virial mass in the Phoenix clusters.
The detection of clusters of galaxies in large surveys plays an important part in extragalactic astronomy, and particularly in cosmology, since cluster counts can give strong constraints on cosmological parameters. X-ray imaging is in particular a reliable means to discover new clusters, and large X-ray surveys are now available. Considering XMM-Newton data for a sample of 40 Abell clusters, we show that their analysis with a Kolmogorov distribution can provide a distinctive signature for galaxy clusters. The Kolmogorov method is sensitive to the correlations in the cluster X-ray properties and can therefore be used for their identification, thus allowing to search reliably for clusters in a simple way.
We study moderate gravitational lensing where a background galaxy is magnified substantially, but not multiply imaged, by an intervening galaxy. We focus on the case where both the lens and source are elliptical galaxies. The signatures of moderate lensing include isophotal distortions and systematic shifts in the fundamental plane and Kormendy relation, which can potentially be used to statistically determine the galaxy mass profiles. These effects are illustrated using Monte Carlo simulations of galaxy pairs where the foreground galaxy is modelled as a singular isothermal sphere model and observational parameters appropriate for the Large Synoptic Survey Telescope (LSST). The range in radius probed by moderate lensing will be larger than that by strong lensing, and is in the interesting regime where the density slope may be changing.
We investigate the effect of photodesorption on the snow line position at the surface of a protoplanetary disk around a Herbig Ae/Be star, motivated by the detection of water ice particles at the surface of the disk around HD142527 by Honda et al. For this aim, we obtain the density and temperature structure in the disk with a 1+1D radiative transfer and determine the distribution of water ice particles in the disk by the balance between condensation, sublimation, and photodesorption. We find that photodesorption induced by the far-ultraviolet radiation from the central star depresses the ice-condensation front toward the mid-plane and pushes the surface snow line outward significantly when the stellar effective temperature exceeds a certain critical value. This critical effective temperature depends on the stellar luminosity and mass, the water abundance in the disk, and the yield of photodesorption. We present an approximate analytic formula for the critical temperature. We separate Herbig Ae/Be stars into two groups on the HR diagram according to the critical temperature; one is the disks where photodesorption is effective and from which we may not find ice particles at the surface, and the other is the disks where photodesorption is not effective. We estimate the snow line position at the surface of the disk around HD142527 to be 100--300 AU which is consistent with the water ice detection at >140 AU in the disk. All results depend on the dust grain size by a complex way and this point requires more work in future.
We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.
We present high-cadence observations and simulations of the solar photosphere, obtained using the Rapid Oscillations in the Solar Atmosphere imaging system and the MuRAM magneto-hydrodynamic code, respectively. Each dataset demonstrates a wealth of magneto-acoustic oscillatory behaviour, visible as periodic intensity fluctuations with periods in the range 110-600 s. Almost no propagating waves with periods less than 140s and 110s are detected in the observational and simulated datasets, respectively. High concentrations of power are found in highly magnetised regions, such as magnetic bright points and intergranular lanes. Radiative diagnostics of the photospheric simulations replicate our observational results, confirming that the current breed of magneto-hydrodynamic simulations are able to accurately represent the lower solar atmosphere. All observed oscillations are generated as a result of naturally occurring magnetoconvective processes, with no specific input driver present. Using contribution functions extracted from our numerical simulations, we estimate minimum G-band and 4170 Angstrom continuum formation heights of 100 km and 25 km, respectively. Detected magneto-acoustic oscillations exhibit a dominant phase delay of -8 degrees between the G-band and 4170 Angstrom continuum observations, suggesting the presence of upwardly propagating waves. More than 73% of MBPs (73% from observations, 96% from simulations) display upwardly propagating wave phenomena, suggesting the abundant nature of oscillatory behaviour detected higher in the solar atmosphere may be traced back to magnetoconvective processes occurring in the upper layers of the Sun's convection zone.
Recent estimates of the Cepheid distance modulus of NGC 6822 differ by 0.18 mag. To investigate this we present new multi-epoch JHKs photometry of classical Cepheids in the central region of NGC 6822 and show that there is a zero-point difference from earlier work. These data together with optical and mid-infrared observations from the literature are used to derive estimates of the distance modulus of NGC 6822. A best value of 23.40 mag is adopted, based on an LMC distance modulus of 18.50 mag. The standard error of this quantity is ~0.05 mag. We show that to derive consistent moduli from Cepheid observations at different wavelengths, it is necessary that the fiducial LMC period-luminosity relations at these wavelengths should refer to the same subsample of stars. Such a set is provided. A distance modulus based on RR Lyrae variables agrees with the Cepheid result.
Using spectroscopic data taken with Keck II DEIMOS by the Z-PAndAS team in the Andromeda-Triangulum region, I present a comparison of the disc and satellite systems of Andromeda with those of our own Galaxy. I discuss the observed discrepancies between the masses and scale radii of Andromeda dwarf spheroidal galaxies of a given luminosity with those of the Milky Way. I also also present an analysis of the newly discovered M31 thick disc, which is measured to be hotter, more extended and thicker than that seen in the Milky Way.
Factors contributing to the scatter around the ridge-line period-luminosity relationship are listed, followed by a discussion how to eliminate the adverse effects of these factors (mode of pulsation, crossing number, temperature range, reddening, binarity, metallicity, non-linearity of the relationship, blending), in order to reduce the dispersion of the P-L relationship.
Based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan telescopes we present detailed elemental abundances for 64 red giant stars in the inner and outer Galactic disk. For the inner disk sample (4-7 kpc from the Galactic centre) we find that stars with both thin and thick disk abundance patterns are present while for Galactocentric distances beyond 10 kpc, we only find chemical patterns associated with the local thin disk, even for stars far above the Galactic plane. Our results show that the relative densities of the thick and thin disks are dramatically different from the solar neighbourhood, and we therefore suggest that the radial scale length of the thick disk is much shorter than that of the thin disk. A thick disk scale-length of L_{thick}=2.0 kpc, and L_{thin}=3.8 kpc for the thin disk, better match the data.
We highlight some results from our high-resolution spectroscopic elemental abundance survey of F and G dwarf stars in the solar neighbourhood.
We have determined detailed elemental abundances and stellar ages for a sample of now 38 microlensed dwarf and subgiant stars in the Galactic bulge. Stars with sub-solar metallicities are all old and have enhanced alpha-element abundances -- very similar to what is seen for local thick disk stars. The metal-rich stars on the other hand show a wide variety of stellar ages, ranging from 3-4 Gyr to 12 Gyr, and an average around 7-8 Gyr. The existence of young and metal-rich stars are in conflict with recent photometric studies of the bulge which claim that the bulge only contains old stars.
In order to complete Stark broadening data for Ne II, and O III lines, needed for analysis of stellar atmospheres, we determined, within the semiclassical perturbation method, the missing Stark broadening parameters for the broadening by collisions with protons and ionized helium, for 15 Ne II and 5 O III multiplets. Also, electron, proton, and ionized helium impact broadening parameters for an important Ne II multiplet in the visible part of the spectrum, and for three Ne III multiplets, were calculated. The obtained data will be included in the STARK-B database, which is a part of Virtual Atomic and Molecular Data Center.
In order to complete data on Stark broadening parameters for Ar I line in the visible spectrum, we determined Stark widths and shifts due to electron, proton, and ionized helium impacts, for nine lines (4191.0, 4259.4, 5912.1, 6043.2, 6045.0, 6752.9, 7503.9, 7514.6, 7724.2 {\AA}), using jK coupling and semiclassical-perturbation theory. The obtained results will enter the STARK-B database, which is a part of Virtual Atomic and Molecular Data Center.
White dwarf and pre-white dwarf atmospheres are one of the best examples for the application of Stark broadening research results in astrophysics, due to plasma conditions very favorable for this line broadening mechanism. For example in hot hydrogen-deficient (pre-) white dwarf stars Teff = 75 000 K - 180 000 K and log g = 5.5-8 [cgs]. Even for much cooler DA and DB white dwarfs with typical effective temperatures of 10 000 K - 20 000 K, Stark broadening is usually the dominant broadening mechanism. In this review, Stark broadening in white dwarf spectra is considered and the attention is drawn to the STARK-B database (this http URL), containing Stark broadening parameters needed for white dwarf spectra analysis and synthesis, as well as to the new search facilities which will provide the collective effort to develop Virtual Atomic and Molecular Data Center (VAMDC - this http URL).
Stark broadening theory is currently operated for calculating widths and shifts of spectral lines that are needed for spectroscopic diagnostics and modelling in astrophysics, laboratory and technological plasmas. We have calculated a great number of data, obtained through the impact semi- classical perturbation theory: tables have been published for neutral atom and ion emitters, and typical temperatures, electron and ion densities. They are currently implemented in the STARK-B database which participates to the European effort within the VAMDC (Virtual Atomic and Molecular data Centre). Despite of that, a great number of data are still missing and their orders of magnitude would at least be welcome. In the present paper, we will revisit and compare the orders of magnitudes and trends of the impact Stark widths and shifts, by considering their semiclassical perturbation expressions. We will also provide fitting formulae which are essential for the modelling codes of stellar atmospheres and envelope
Besides the need of Stark broadening parameters for a number of problems in physics, and plasma technology, in hot star atmospheres the conditions exist where Stark widths are comparable and even larger than the thermal Doppler widths. Using the semiclassical perturbation method we investigated here the influence of collisions with charged particles for In III spectral lines. We determined a number of Stark broadening parameters important for the investigation of plasmas in the atmospheres of A-type stars and white dwarfs. Also, we have compared the obtained results with existing experimental data. The results will be included in the STARK-B database, the Virtual Atomic and Molecular Data Center and the Serbian Virtual Observatory.
Based on two models, we investigate the molecular-to-atomic gas ratio in Virgo cluster galaxies in comparison with field galaxies. We show that the enhanced metallicity for cluster members and the ram pressure stripping of atomic gas from the disk periphery cannot fully explain the observed gas component ratios. The additional environmental factors affecting the interstellar medium and leading to an increase in the molecular gas fraction should be taken into account for cluster galaxies.
The aim of Dawn mission is the acquisition of data from orbits around two
bodies, (4)Vesta and (1)Ceres, the two most massive asteroids. Due to the low
thrust propulsion, Dawn will slowly cross and transit through ground-track
resonances, where the perturbations on Dawn orbit may be significant. In this
context, to safety go the Dawn mission from the approach orbit to the lowest
science orbit, it is essential to know the properties of the crossed
resonances. This paper analytically investigates the properties of the major
ground-track resonances (1:1, 1:2, 2:3 and 3:2) appearing for Vesta orbiters:
location of the equilibria, aperture of the resonances and period at the stable
equilibria. We develop a general method using an averaged Hamiltonian
formulation with a spherical harmonic approximation of the gravity field. If
the values of the gravity field coefficient change, our method stays correct
and applicable. We also discuss the effect of one uncertainty on the C20 and
C22 coefficients on the properties of the 1:1 resonance. These results are
checked by numerical tests. We determine that the increase of the eccentricity
appearing in the 2:3 resonance is due to the C22 and S22 coefficients.
Our method can be easily adapted to missions similar to Dawn because,
contrarily to the numerical results, the analytical formalism stays the same
and is valid for a wide range of physical parameters of the asteroid (namely
the shape and the mass) as well as for different spacecraft orbits.
Finally we numerically study the probability of the capture in resonance 1:1.
Our paper reproduces, explains and supplements the results of Tricarico and
Sykes (2010).
The planetary nebula (PN) NGC 6778 harbors a binary central star with a short orbital period and displays two systems of fast collimated outflows. In order to assess the influence of the evolution through a common-envelope phase of the binary system of NGC 6778 on its formation and shaping, we have used narrow-band images and high-dispersion long-slit spectra of the nebula to investigate its detailed morphology and kinematics. We find that the overall structure of NGC 6778 can be described as a bipolar PN. The equatorial ring is highly disrupted and many radial features (filamentary wisps and cometary knots) also evidence strong dynamical effects. There are clear connections between the bipolar lobes and the fast collimated outflows: the collimated outflows seem to arise from bright knots at the tips of the bipolar lobes, whereas the kinematics of the bipolar lobes is distorted. We suggest that the interaction of the fast collimated outflows of NGC 6778 with its nebular envelope has resulted in the disruption of the nebular shell and equatorial ring.
1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high energy (VHE, E>100 GeV) gamma-rays due to its high X-ray and radio flux. Observations in the VHE gamma-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic particles present in the source and the radiation processes at work. Because of the distance of the source, the gamma-ray spectrum might provide further limits on the level of the Extragalactic Background Light (EBL). We report observations made between October 2005 and December 2009 with H.E.S.S., an array of four imaging atmospheric Cherenkov telescopes. Observations at high energies (HE, 100 MeV - 100 GeV) with the Fermi-LAT instrument in the first 20 months of its operation are also reported. To complete the multi-wavelength picture, archival UV and X-ray observations with the Swift satellite and optical observations with the ATOM telescope are also used. Based on the observations with H.E.S.S., 1ES 0414+009 is detected for the first time in the VHE band. An excess of 224 events is measured, corresponding to a significance of 7.8 sigma. The photon spectrum of the source is well described by a power law, with photon index of 3.45 \pm 0.25stat \pm 0.20syst. The integral flux above 200 GeV is (1.88 \pm 0.20stat \pm 0.38syst) \times10-12 cm-2 s-1. Observations with the Fermi-LAT in the first 20 months of operation show a flux between 200 MeV and 100 GeV of (2.3 \pm 0.2stat) \times 10-9 erg cm-2 s-1, and a spectrum well described by a power-law function with a photon index 1.85 \pm 0.18. Swift/XRT observations show an X-ray flux between 2 and 10 keV of (0.8 - 1) \times 10-11 erg cm-2 s-1, and a steep spectrum (2.2 - 2.3). Combining X-ray with optical-UV data, a fit with a log-parabolic function locates the synchrotron peak around 0.1 keV. ...
In previous works of this series, we have shown that late B- and early A-type stars have genuine bimodal distributions of rotational velocities and that late A-type stars lack slow rotators. The distributions of the surface angular velocity ratio \Omega/\Omega_crit (\Omega_crit is the critical angular velocity) have peculiar shapes according to spectral type groups, which can be caused by evolutionary properties. We aim to review the properties of these rotational velocity distributions in some detail as a function of stellar mass and age. We have gathered v sin i for a sample of 2014 B6- to F2-type stars. We have determined the masses and ages for these objects with stellar evolution models. The (Teff, log L/Lsun)-parameters were determined from the uvby-\beta photometry and the HIPPARCOS parallaxes. The velocity distributions show two regimes that depend on the stellar mass. Stars less massive than 2.5 Msun have a unimodal equatorial velocity distribution and show a monotonical acceleration with age on the main sequence (MS). Stars more massive have a bimodal equatorial velocity distribution. Contrarily to theoretical predictions, the equatorial velocities of stars from about 1.7 Msun to 3.2 Msun undergo a strong acceleration in the first third of the MS evolutionary phase, while in the last third of the MS they evolve roughly as if there were no angular momentum redistribution in the external stellar layers. The studied stars might start in the ZAMS not necessarily as rigid rotators, but with a total angular momentum lower than the critical one of rigid rotators. The stars seem to evolve as differential rotators all the way of their MS life span and the variation of the observed rotational velocities proceeds with characteristic time scales \delta(t)\sim 0.2 t_MS, where t_MS is the time spent by a star in the MS.
I briefly review the progress of accretion disc theory since the seminal paper of Shakura and Sunyaev.
The 6 cm and 20 cm radio continuum properties of all 85 galaxies with reported 22 GHz H2O maser emission and luminosity distance D > 0.5 Mpc are studied. For the total of 55 targets for which both 6 cm and 20 cm measurements exist and for the subsample of 42 sources with masers related to active galactic nuclei (AGN), a spectral index could be determined from an assumed power-law dependence. The mean value of the resulting spectral index is in both cases 0.66+-0.07. Comparing radio properties of the maser galaxies with a sample of Seyferts without detected H2O maser, we find that (1) the spectral indices agree within the error limits, and (2) maser host galaxies have higher nuclear radio continuum luminosities, exceeding those of the comparison sample by factors of order 5. Only considering the subsample of galaxies with masers associated with AGN, there seems to be a trend toward rising maser luminosity with nuclear radio luminosity (both at 6 cm and 20 cm). However, when accounting for the Malmquist effect, the correlation weakens to a level, which is barely significant. Overall, the study indicates that nuclear radio luminosity is a suitable indicator to guide future AGN maser searches and to enhance detection rates, which are otherwise quite low (<10%).
The transit timing variation (TTV) method allows the detection of non-transiting planets through their gravitational perturbations. Since TTVs are strongly enhanced in systems close to mean-motion resonances (MMR), even a low mass planet can produce an observable signal. This technique has thus been proposed to detect terrestrial planets. In this letter, we analyse TTV signals for systems in or close to MMR in order to illustrate the difficulties arising in the determination of planetary parameters. TTVs are computed numerically with an n-body integrator for a variety of systems close to MMR. The main features of these TTVs are also derived analytically. Systems deeply inside MMR do not produce particularly strong TTVs, while those close to MMR generate quasiperiodic TTVs characterised by a dominant long period term and a low amplitude remainder. If the remainder is too weak to be detected, then the signal is strongly degenerate and this prevents the determination of the planetary parameters. Even though an Earth mass planet can be detected by the TTV method if it is close to a MMR, it may not be possible to assert that this planet is actually an Earth mass planet. On the other hand, if the system is right in the center of a MMR, the high amplitude oscillation of the TTV signal vanishes and the detection of the perturber becomes as difficult as it is far from MMR.
In setting up initial conditions for cosmological N-body simulations there are, fundamentally, two choices: either maximizing the correspondence of the initial density field to the assumed fourier-space clustering or, instead, matching to the real-space clustering. As a stringent test of both approaches, I perform ensembles of simulations using power law models and exploit the self-similarity of these initial conditions to quantify the accuracy of the results. Originally proposed by Pen 1997 and implemented by Sirko 2005, I show that the real-space motivated approach, which allows the DC mode to vary, performs well in exhibiting the expected self-similar behavior in the mean xi(r) and P(k) and in both methods this behavior extends below the scale of the initial mean interparticle spacing. I also test the real-space method with simulations of a simplified, powerlaw model for baryon acoustic oscillations, again with success, and mindful of the need to generate mock catalogs using simulations I show extensive powerlaw tests for the halo mass function and halo bias in our simulations. Although requiring a few to many times more simulations than the standard, fourier-space method to reach a given certainty on the correlation function or power spectrum, I find that the real-space method is more reliable for modeling P(k) when the clustering level becomes significant on the scale of the simulation box. As such a carefully-constructed real-space approach could be optimal for simulating extremely red power spectra (n_eff < -2), as in excessively small box simulations to model the "end" of the CDM hierarchy. I conclude by discussing some possibilities for optimizing the real-space method for more general use and an appendix demonstrates the potential for using perturbation theory to model the effect of the box scale on the simulated growth of structure.
The hadronic model of Active Galactic Nuclei and other compact high energy astrophysical sources assumes that ultra-relativistic protons, electron-positron pairs and photons interact via various hadronic and electromagnetic processes inside a magnetized volume, producing the multiwavelength spectra observed from these sources. A less studied property of such systems is that they can exhibit a variety of temporal behaviours due to the operation of different feedback mechanisms. We investigate the effects of one possible feedback loop, where \gamma-rays produced by photopion processes are being quenched whenever their compactness increases above a critical level. This causes a spontaneous creation of soft photons in the system that result in further proton cooling and more production of \gamma-rays, thus making the loop operate. We perform an analytical study of a simplified set of equations describing the system, in order to investigate the connection of its temporal behaviour with key physical parameters. We also perform numerical integration of the full set of kinetic equations verifying not only our analytical results but also those of previous numerical studies. We find that once the system becomes `supercritical', it can exhibit either a periodic behaviour or a damped oscillatory one leading to a steady state. We briefly point out possible implications of such a supercriticality on the parameter values used in Active Galactic Nuclei spectral modelling, through an indicative fitting of the VHE emission of blazar 3C 279.
As part of an ongoing collaboration between student groups at high schools and professional astronomers, we have searched for the presence of circum-binary planets in a bona-fide unbiased sample of twelve post-common envelope binaries (PCEBs) from the Catalina Sky Survey (CSS) and the Sloan Digital Sky Survey (SDSS). Although the present ephemerides are significantly more accurate than previous ones, we find no clear evidence for orbital period variations between 2005 and 2011 or during the 2011 observing season. The sparse long-term coverage still permits O-C variations with a period of years and an amplitude of tens of seconds, as found in other systems. Our observations provide the basis for future inferences about the frequency with which planet-sized or brown-dwarf companions have either formed in these evolved systems or survived the common envelope (CE) phase.
Rapidly rotating Neutron Stars in Low Mass X-ray Binaries (LMXBs) may be an interesting source of Gravitational Waves (GWs). In particular, several modes of stellar oscillation may be driven unstable by GW emission, and this can lead to a detectable signal. Here we illustrate how current X-ray and ultra-violet (UV) observations can constrain the physics of the r-mode instability. We show that the core temperatures inferred from the data would place many systems well inside the unstable region predicted by standard physical models. However, this is at odds with theoretical expectations. We discuss different mechanisms that could be at work in the stellar interior, and we show how they can modify the instability window and make it consistent with the inferred temperatures.
We report the discovery of 8.5 sigma high-frequency quasi-periodic oscillations (HFQPOs) at 66 Hz in the RXTE data of the black hole candidate IGR J17091-3624, a system whose X-ray properties are very similar to those of microquasar GRS 1915+105. The centroid frequency of the strongest peak is ~66 Hz, its quality factor above 5 and its rms is between 4 and 10%. We found a possible additional peak at 164 Hz when selecting a subset of data; however, at 4.5 sigma level we consider this detection marginal. These QPOs have hard spectrum and are stronger in observations performed between September and October 2011, during which IGR J17091-3624 displayed for the first time light curves which resemble those of the gamma variability class in GRS 1915+105. We find that the 66 Hz QPO is also present in previous observations (4.5 sigma), but only when averaging ~235 ksec of relatively high count rate data. The fact that the HFQPOs frequency in IGR J17091-3624 matches surprisingly well that seen in GRS 1915+105 raises questions on the mass scaling of QPOs frequency in these two systems. We discuss some possible interpretations, however, they all strongly depend on the distance and mass of IGR J17091-3624, both completely unconstrained today.
The cause for the observed acceleration in the expansion of the universe is unknown, and referred to as "dark energy" for convenience. Dark energy could be an unknown energy component, or a modification of Einstein's general relativity. This dictates the measurements that are optimal in unveiling the nature of dark energy: the cosmic expansion history, and the growth history of cosmic large scale structure. I will examine Type Ia supernovae and galaxy clustering as dark energy probes, and discuss recent results and future prospects.
Unification Models of Active Galactic Nuclei postulate that all the observed differences between Type 1 and Type 2 objects are due to orientation effects with respect to the line-of-sight to the observer. The key ingredient of these models is the obscuring medium, historically envisaged as a toroidal structure on a parsec scale. However, many results obtained in the last few years are clearly showing the need for a more complex geometrical distribution of the absorbing media. In this paper we review the various pieces of evidence for obscuring media on different scales, from the vicinity of the black hole to the host galaxy, in order to picture an updated unification scenario explaining the complex observed phenomenology. We conclude by mentioning some of the open issues.
A new EAS Cherenkov light array, Tunka-133, with ~1 km^2 geometrical area has been installed at the Tunka Valley (50 km from Lake Baikal) in 2009. The array permits a detailed study of cosmic ray energy spectrum and mass composition in the energy range 10^16 - 10^18 eV with a uniform method. We describe the array construction, DAQ and methods of the array calibration.The method of energy reconstruction and absolute calibration of measurements are discussed. The analysis of spatial and time structure of EAS Cherenkov light allows to estimate the depth of the EAS maximum X_max. The results on the all particles energy spectrum and the mean depth of the EAS maximum X_max vs. primary energy derived from the data of two winter seasons (2009 -- 2011), are presented. Preliminary results of joint operation of the Cherenkov array with antennas for detection of EAS radio signals are shown. Plans for future upgrades -- deployment of remote clusters, radioantennas and a scintillator detector network and a prototype of the HiSCORE gamma-telescope -- are discussed.
The BATSE catalogue is searched for evidence of spindown of black holes or proto-neutron stars (PNS) by extracting normalized light curves (nLC). The nLC are obtained by matched filtering, to suppress intermediate time scales such as due to the shock break-out of GRB jets through a remnant stellar envelope. We find consistency within a few percent of the nLC and the model template for spindown of an initially extremal black hole against high-density matter at the ISCO. The large BATSE size enables a study of the nLC as a function of durations $T_{90}$. The resulting $\chi^2_{red}$ is within a $2.35\sigma$ confidence interval for durations $T_{90}>20$ s, which compares favorably with the alternative of spindown against matter further out and spindown of a PNS, whose $\chi^2$ fits are, respectively, outside the 4$\sigma$ and 12$\sigma$ confidence intervals. We attribute spindown against matter at the ISCO to cooling by gravitational-wave emission from non-axisymmetric instabilities in the inner disk or torus as the result of a Hopf bifurcation in response to energetic input from the central black hole. This identification gives an attractive outlook for chirps in quasi-periodic gravitational waves lasting tens of seconds of interest to LIGO, Virgo and the LCGT.
Astronomical instruments generally possess spatially variant point-spread functions, which determine the amount by which an image pixel is blurred as a function of position. Several techniques have been devised to handle this variability in the context of the standard image deconvolution problem. We have developed an iterative gravitational lens modeling code called Mirage that determines the parameters of pixelated source intensity distributions for a given lens model. We are able to include the effects of spatially variant point-spread functions using the iterative procedures in this lensing code. In this paper, we discuss the methods to include spatially variant blurring effects and test the results of the algorithm in the context of gravitational lens modeling problems.
The dense molecular clump at the center of the Barnard 59 (B59) complex is the only region in the Pipe Nebula that has formed a small,stellar cluster. The previous analysis of a high resolution near-IR dust extinction map revealed that the nuclear region in B59 is a massive, mostly quiescent clump of 18.9 solar masses. The clump shows a monolithic profile, possibly indicating that the clump is on the way to collapse, with no evident fragmentation that could lead to another group of star systems. In this paper we present new analysis that compares the dust extinction map with a new dust emission radio-continuum map of higher spatial resolution. We confirm that the clump does not show any significant evidence for prestellar fragmentation at scales smaller than those probed previously.
We study the plasma correlation effects on nonresonant thermonuclear reactions of carbon and oxygen in the interiors of white dwarfs and liquid envelopes of neutron stars. We examine the effects of electron screening on thermodynamic enhancement of thermonuclear reactions in dense plasmas beyond the linear mixing rule. Using these improved enhancement factors, we calculate carbon and oxygen ignition curves in white dwarfs and neutron stars. The energy balance and ignition conditions in neutron star envelopes are evaluated, taking their detailed thermal structure into account. The result is compared to the simplified "one-zone model," which is routinely used in the literature. We also consider the effect of strong magnetic fields on the ignition curves in the ocean of magnetars.
The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)
Context. A number of millimeter and submillimeter facilities with linear polarization observing capabilities have started operating during last years. These facilities, as well as other previous millimeter telescopes and interferometers, require bright and stable linear polarization calibrators to calibrate new instruments and to monitor their instrumental polarization. The current limited number of adequate calibrators implies difficulties in the acquisition of these calibration observations. Aims. Looking for additional linear polarization calibrators in the millimeter spectral range, in mid-2006 we started monitoring 3C 286, a standard and highly stable polarization calibrator for radio observations. Methods. Here we present the 3 and 1mm monitoring observations obtained between September 2006 and October 2011 with the XPOL polarimeter on the IRAM 30m Millimeter Telescope. Results. Our observations show that 3C 286 is a bright source of constant total flux with 3mm flux density S_3mm = (0.90 \pm 0.02) Jy. The 3mm linear polarization degree (p_3mm = [13.6\pm 0.3]%) and polarization angle (chi_3mm = [37.5\pm 0.8] deg, expressed in the equatorial coordinate system) are also constant during the time span of our observations. Although with poorer time sampling and data quality, our 1mm observations of 3C 286 are also consistent with a constant 1mm source with flux density S_1mm = (0.29 \pm 0.03) Jy. Indeed, our statistical analysis demonstrate that the linear polarization degree and polarization angle of 3C 286 at 1mm (p_1mm = [16.1\pm1.6]%, and chi_1mm = [33.0 \pm 6.8] deg., respectively) are also constant. Conclusions. This, together with the previously known compact structure of 3C 286 -extended by \sim 3.5" in the sky- allow us to propose 3C 286 as a new calibrator for both single dish and interferometric polarization observations at 3mm, and possibly at shorter wavelengths.
Analytic estimates of the viscous time-scale due to cloud-cloud collisions have been as high as thousands of Gyr. Consequently, cloud collisions are widely ignored as a source of viscosity in galactic disks. However, capturing the hydrodynamics of discs in simple analytic models is a challenge, both because of the wide dynamic range and importance of 2D and 3D effects. To test the validity of analytic models we present estimates for the viscous time-scale that are measured from three dimensional SPH simulations of disc formation and evolution. We have deliberately removed uncertainties associated with star-formation and feedback thereby enabling us to place lower bounds on the time-scale for this process. We also contrast collapse simulations with results from simulations of initially stable discs and examine the impact of numerical parameters and assumptions on our work, to constrain possible systematics in our estimates. We find that cloud-collision viscous time-scales are in the range of 0.6-16 Gyr, considerably shorter than previously estimated. This large discrepency can be understood in terms of how the efficiency of collisions is included in the analytical estimates. We find that the viscous time-scale only depends weakly on the number of clouds formed, and so while the viscous time-scale will increase with increasing resolution, this effect is too weak to alter our conclusions.
Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2-3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The non-radiative simulations overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation is responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high-redshift z > 2.
We present a detailed analysis of the X-ray spectrum of the Seyfert 2 galaxy NGC454E, belonging to the interacting system NGC454. Observations performed with Suzaku, XMM-Newton and Swift allowed us to detect a dramatic change in the curvature of the 2-10 keV spectrum, revealing a significant variation of the absorbing column density along the line of sight (from ~ 1 x10^{24}cm^{-2} to ~ 1x10^{23}cm^{-2}). Consequently, we propose this source as a new member of the class of "changing look" AGN, i.e. AGN that have been observed both in Compton-thin (NH =10^{23 cm^{-2}) and reflection dominated states (Compton-thick, NH >10^{24} cm^{-2}). Due to the quite long time lag (6 months) between the Suzaku and XMM-Newton observations we cannot infer the possible location of the obscuring material causing the observed variability. In the 6-7 keV range the XMM-Newton observation also shows a clear signature of the presence of an ionized absorber. Since this feature is not detected during the Suzaku observation (despite its detectability), the simplest interpretation is that the ionized absorber is also variable; its location is estimated to be within ~10^{-3} pc from the central black hole, probably much closer in than the rather neutral absorber.
We explore the newly discovered "hangup-kick" effect, which greatly amplifies the recoil for configuration with partial spin- orbital-angular momentum alignment, by studying a set of 48 new simulations of equal-mass, spinning black-hole binaries. We propose a phenomenological model for the recoil that takes this new effect into account and then use this model, in conjunction with statistical distributions for the spin magnitude and orientations, based on accretion simulations, to find the probabilities for observing recoils of several thousand km/s. In addition, we provide initial parameters, eccentricities, radiated linear and angular momentum, precession rates and remnant mass, spin, and recoils for all 48 configurations. Our results indicate that surveys exploring peculiar (redshifted or blueshifted) differential line-of-sight velocities should observe at least one case above 2000 km/s out of four thousand merged galaxies. The probability that a remnant BH receives a total recoil exceeding the ~2000 km/s escape velocity of large elliptical galaxies is ten times larger. Probabilities of recoils exceeding the escape velocity quickly rise to 5% for galaxies with escape velocities of 1000 km/s and nearly 20% for galaxies with escape velocities of 500 km/s. In addition the direction of these large recoils is strongly peaked toward the angular momentum axis, with very low probabilities of recoils exceeding 350 km/s for angles larger than 45 deg. with respect to the orbital angular momentum axis.
We study oscillons, extremely long-lived localized oscillations of a scalar field, with three different potentials: quartic, sine-Gordon model and in a new class of convex potentials. We use an absorbing boundary at the end of the lattice to remove emitted radiation. The energy and the frequency of an oscillon evolve in time and are well fitted by a constant component and a decaying, radiative part obeying a power law as a function time. The power spectra of the emitted radiation show several distinct frequency peaks where oscillons release energy. In two dimensions, and with suitable initial conditions, oscillons do not decay within the range of the simulations, which in quartic theory reach 10^8 time units. While it is known that oscillons in three-dimensional quartic theory and sine-Gordon model decay relatively quickly, we observe a surprising persistence of the oscillons in the convex potential with no sign of demise up to 10^7 time units. This leads us to speculate that an oscillon in such a potential could actually live infinitely long both in two and three dimensions.
The aim of this paper is to report on the existence of a wide variety of exact solutions, ranging from black holes to wormholes, when a conformally coupled scalar field with a self interacting potential containing a linear, a cubic and a quartic self interaction is taken as a source of the energy-momentum tensor, in the Einstein theory with a cosmological constant. Among all the solutions there are two particularly interesting. On the one hand, the spherically symmetric black holes when the cosmological constant is positive; they are shown to be everywhere regular, namely there is no singularity neither inside nor outside the event horizon. On the other hand, there are spherically symmetric and topological wormholes that connect two asymptotically (anti) de Sitter regions with a different value for the cosmological constant. The regular black holes and the wormholes are supported by everywhere regular scalar field configurations.
We report on the multipolar equations of motion for spinning test bodies in the de Sitter spacetime of constant positive curvature. The dynamics of spinning particles is discussed for the two supplementary conditions of Frenkel and Tulczyjew. Furthermore, the 4-momentum and the spin are explicitly expressed in terms of the spacetime coordinates with the help of the 10 Killing vectors available in de Sitter spacetime.
The cosmological consequences of the f(R) gravity are reviewed in the framework of recent data obtained by PAMELA (Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics) experiment. This collaboration has reported an excess of positron events that cannot be explained by conventional cosmology and particle physics, and are usually ascribed to the dark matter presence (in particular, weak interacting massive particles). The dark matter interpretation of PAMELA data has motivated the study of alternative cosmological models (with respect to the standard cosmology) owing to the fact that they predict an enhancement of the Hubble expansion rate, giving rise in such a way to thermal relics with a larger relic abundance. In this paper we face this problem for f(R) gravity models assuming a power law correction to the standard action of General Relativity, i.e. f(R) = R + \alpha R^n. In the regime in which the energy density induced by (higher order) curvature terms are greater than (or of the order of) the energy density of standard radiation, we find that the exponent n must assume values n > 1. In the opposite regime, taking into account the constraints provided by big bang nucleosynthesis, we get n < 1 + 2 \times 10^-3. The latter bound excludes the model f(R) = R + \alpha R^2. Our analysis shows that the considered model allows to explain the PAMELA puzzle for dark matter relic particles with masses greater or lesser than 10^2 GeV.
We determine the sensitivity to a possible variation of the proton-to-electron mass ratio \mu for torsion-wagging-rotation transitions in the ground state of methylamine (CH3NH2). Our calculation uses an effective Hamiltonian based on a high-barrier tunneling formalism combined with extended-group ideas. The \mu-dependence of the molecular parameters that are used in this model are derived and the most important ones of these are validated using the spectroscopic data of different isotopologues of methylamine. We find a significant enhancement of the sensitivity coefficients due to energy cancellations between internal rotational, overall rotational and inversion energy splittings. The sensitivity coefficients of the different transitions range from -19 to +24. The sensitivity coefficients of the 78.135, 79.008, and 89.956 GHz transitions that were recently observed in the disk of a z = 0.89 spiral galaxy located in front of the quasar PKS 1830-211 [S. Muller et al. Astron. Astrophys. 535, A103 (2011)] were calculated to be -0.87 for the first two and -1.4 for the third transition, respectively. From these transitions a preliminary upper limit for a variation of the proton to electron mass ratio of \Delta \mu/\mu< 9 x 10^{-6} is deduced.
We study the dynamics in the neighborhood of simple and double unstable periodic orbits in a rotating 3D autonomous Hamiltonian system of galactic type. In order to visualize the four dimensional spaces of section we use the method of color and rotation. We investigate the structure of the invariant manifolds that we found in the neighborhood of simple and double unstable periodic orbits in the 4D spaces of section. We consider orbits in the neighborhood of the families x1v2, belonging to the x1 tree, and the z-axis (the rotational axis of our system). Close to the transition points from stability to simple instability, in the neighborhood of the bifurcated simple unstable x1v2 periodic orbits we encounter the phenomenon of stickiness as the asymptotic curves of the unstable manifold surround regions of the phase space occupied by rotational tori existing in the region. For larger energies, away from the bifurcating point, the consequents of the chaotic orbits form clouds of points with mixing of color in their 4D representations. In the case of double instability, close to x1v2 orbits, we find clouds of points in the four dimensional spaces of section. However, in some cases of double unstable periodic orbits belonging to the z-axis family we can visualize the associated unstable eigensurface. Chaotic orbits close to the periodic orbit remain sticky to this surface for long times (of the order of a Hubble time or more). Among the orbits we studied we found those close to the double unstable orbits of the x1v2 family having the largest diffusion speed.
We present an updated analysis of the constrained minimal supersymmetric standard model with mu>0 supplemented by an `asymptotic' Yukawa coupling quasi-unification condition, which allows an acceptable b-quark mass. Imposing constraints from the cold dark matter abundance in the universe, B physics, the muon anomalous magnetic moment, and the mass m_h of the lightest neutral CP-even Higgs boson, we find that the lightest neutralino cannot act as a cold dark matter candidate. This is mainly because the upper bound on the lightest neutralino relic abundance from cold dark matter considerations, despite the fact that this abundance is drastically reduced by neutralino-stau coannihilations, is incompatible with the recent data on the branching ratio of B_s --> mu^+ mu^-. Allowing for a different particle, such as the axino or the gravitino, to be the lightest supersymmetric particle and, thus, constitute the cold dark matter in the universe, we find that the predicted m_h's in our model favor the range (119-126) GeV.
The gravitational wave signature emitted from a merging binary depends on the orientation of an observer relative to the binary. Previous studies suggest that emission along the total initial or total final angular momenta leads to both the strongest and simplest signal from a precessing compact binary. In this paper we describe a concrete counterexample: a binary with $m_1/m_2=4$, $a_1=0.6 \hat{x} = -a_2$, placed in orbit in the x,y plane. We extract the gravitational wave emission along several proposed emission directions, including the initial (Newtonian) orbital angular momentum; the final (~ initial) total angular momentum; and the dominant principal axis of $<L_a L_b>_M$. Using several diagnostics, we show that the suggested preferred directions are not representative. For example, only for a handful of other directions (< 15%) will the gravitational wave signal have comparable shape to the one extracted along each of these fiducial directions, as measured by a generalized overlap (>0.95). We conclude that the information available in just one direction (or mode) does not adequately encode the complexity of orientation-dependent emission for even short signals from merging black hole binaries. Future investigations of precessing, unequal-mass binaries should carefully explore and model their orientation-dependent emission.
Links to: arXiv, form interface, find, astro-ph, recent, 1201, contact, help (Access key information)
We measure the acoustic scale from the angular power spectra of the Sloan Digital Sky Survey III (SDSS-III) Data Release 8 imaging catalog that includes 872,921 galaxies over ~ 10,000 deg^2 between 0.45<z<0.65. The extensive spectroscopic training set of the Baryon Oscillation Spectroscopic Survey (BOSS) luminous galaxies allows precise estimates of the true redshift distributions of galaxies in our imaging catalog. Utilizing the redshift distribution information, we build templates and fit to the power spectra of the data, which are measured in our companion paper, Ho et al. 2011, to derive the location of Baryon acoustic oscillations (BAO) while marginalizing over many free parameters to exclude nearly all of the non-BAO signal. We derive the ratio of the angular diameter distance to the sound horizon scale D_A/r_s= 9.212 + 0.416 -0.404 at z=0.54, and therefore, D_A= 1411+- 65 Mpc at z=0.54; the result is fairly independent of assumptions on the underlying cosmology. Our measurement of angular diameter distance D_A is 1.4 \sigma higher than what is expected for the concordance LCDM (Komatsu et al. 2011), in accordance to the trend of other spectroscopic BAO measurements for z >~ 0.35. We report constraints on cosmological parameters from our measurement in combination with the WMAP7 data and the previous spectroscopic BAO measurements of SDSS (Percival et al. 2010) and WiggleZ (Blake et al. 2011). We refer to our companion papers (Ho et al. 2011; de Putter et al. 2011) for investigations on information of the full power spectrum.
We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.
Recent gravitational microlensing observations predict a vast population of free-floating giant planets that outnumbers main sequence stars almost twofold. A frequently-invoked mechanism for generating this population is a dynamical instability that incites planet-planet scattering and the ejection of one or more planets in isolated main sequence planetary systems. Here, we demonstrate that this process alone probably cannot represent the sole source of these galactic wanderers. By using straightforward quantitative arguments and N-body simulations, we argue that the observed number of exoplanets exceeds the plausible number of ejected planets per system from scattering. Thus, other potential sources of free-floaters, such as planetary stripping in stellar clusters and post-main-sequence ejection, must be considered.
SHARDS (Survey for High-z Absorption Red & Dead Sources) is an unbiased ultra-deep spectro-photometric survey with GTC@OSIRIS aimed at selecting and studying massive passively evolving galaxies at z=1.0-2.3 using a set of 24 medium-band filters (FWHM\sim17 nm) at 500-950 nm in GOODS-N. Our observing strategy is optimized to detect at z>1 the prominent Mg absorption feature at rest-frame ~280 nm, a distinctive, necessary, and sufficient feature of evolved stellar populations. Nonetheless, the data quality allow a plethora of studies on galaxy populations, including Emission Lines Galaxies (ELGs) about which we have started our first science verification project presented in this contribution.
We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of Herschel, we examine the dust temperature and dust mass surface density distribution. We find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses compared to interarm regions. However, the distribution of dust temperature and mass are not spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the spiral arms lie ahead of the dust surface density peaks. The dust mass surface density correlates well with the distribution of molecular gas as traced by CO (J=3-2) images (JCMT) and the star formation rate as traced by H?2 with a correction for obscured star formation using 24 micron emission. Using HI images from THINGS to trace the atomic gas component, we make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean gas-to-dust ratio of 84 \pm 4 with higher values in the inner region assuming a constant CO-to-H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor that varies with metallicity.
Magnetospheres of pulsars are thought to be filled with plasma, and variations in plasma supply can affect both pulsar emission properties and spin-down rates. A number of recently discovered "intermittent" pulsars switch between two distinct states: an "on", radio-loud state, and an "off", radio-quiet state. Spin-down rates in the two states differ by a large factor, $\sim 1.5-2.5$, which is not easily understood in the context of current models. In this Letter we present self-consistent numerical solutions of "on" and "off" states of intermittent pulsar magnetospheres. We model the "on" state as a nearly ideal force-free magnetosphere with abundant magnetospheric plasma supply. The lack of radio emission in the "off" state is associated with plasma supply disruption that results in lower plasma density on the open field lines. We model the "off" state using nearly vacuum conditions on the open field lines and nearly ideal force-free conditions on the closed field lines, where plasma can remain trapped even in the absence of pair production. The toroidal advection of plasma in the closed zone in the "off" state causes spin-downs that are a factor of $\sim 2$ higher than vacuum values, and we naturally obtain a range of spin-down ratios between the "on" and "off" states, $\sim 1.2-2.9$, which corresponds to a likely range of pulsar inclination angles of $30{-}90^\circ$. We consider the implications of our model to a number of poorly understood but possibly related pulsar phenomena, including nulling, timing noise, and rotating radio transients.
The deflection of ultra-high energy cosmic rays depends on the shape of the injection spectrum of the source and the pervasive cosmic magnetic fields. In this work it is applied the wavelet transform on the sphere to search for energy ordered filamentary structures arisen from magnetic bending. These structures, the so-called multiplets, can bring relevant information concerning the intervening magnetic field
We present the characterization of the star KOI 961, an M dwarf with transit signals indicative of three short-period exoplanets, originally discovered by the Kepler Mission. We proceed by comparing KOI 961 to Barnard's Star, a nearby, well-characterized mid-M dwarf. By comparing colors, optical and near-infrared spectra, we find remarkable agreement between the two, implying similar effective temperatures and metallicities. Both are metal-poor compared to the Solar neighborhood, have low projected rotational velocity, high absolute radial velocity, large proper motion and no quiescent H-alpha emission--all of which is consistent with being old M dwarfs. We combine empirical measurements of Barnard's Star and expectations from evolutionary isochrones to estimate KOI 961's mass (0.13 +/- 0.05 Msun), radius (0.17 +/- 0.04 Rsun) and luminosity (2.40 x 10^(-3.0 +/- 0.3) Lsun). We calculate KOI 961's distance (38.7 +/- 6.3 pc) and space motions, which, like Barnard's Star, are consistent with a high scale-height population in the Milky Way. We perform an independent multi-transit fit to the public Kepler light curve and significantly revise the transit parameters for the three planets. We calculate the false-positive probability for each planet-candidate, and find a less than 1% chance that any one of the transiting signals is due to a background or hierarchical eclipsing binary, validating the planetary nature of the transits. The best-fitting radii for all three planets are less than 1 Rearth, with KOI 961.03 being Mars-sized (Rp = 0.57 +/- 0.18 Rearth), and they represent some of the smallest exoplanets detected to date.
We introduce a new implementation of the FastICA algorithm on simulated LOFAR EoR data with the aim of accurately removing the foregrounds and extracting the 21-cm reionization signal. We find that the method successfully removes the foregrounds with an average fitting error of 0.5 per cent and that the 2D and 3D power spectra are recovered across the frequency range. We find that for scales above several PSF scales the 21-cm variance is successfully recovered though there is evidence of noise leakage into the reconstructed foreground components. We find that this blind independent component analysis technique provides encouraging results without the danger of prior foreground assumptions.
Type Ia supernova (SNe Ia) are thought to originate in the explosion of a white dwarf. The explosion could be triggered by the merger of two white dwarfs ('double-degenerate' origin), or by mass transfer from a companion star (the 'single-degenerate' path). The identity of the progenitor is still controversial; for example, a recent argument against the single-degenerate origin has been widely rejected. One way to distinguish between the double- and single-degenerate progenitors is to look at the center of a known SN Ia remnant to see whether any former companion star is present. A likely ex-companion star for the progenitor of Tycho's supernova has been identified, but that claim is still controversial. Here we report that the central region of the supernova remnant SNR 0509-67.5 (the site of a Type Ia supernova 400+-50 years ago, based on its light echo) in the Large Magellanic Cloud contains no ex-companion star to a limit of V=26.9 magnitude (M_V=+8.4) within the extreme 99.73% region with radius 1.43". The lack of any ex-companion star to deep limits rules out all published single-degenerate models. The only remaining possibility is that the progenitor for this particular SN Ia was a double-degenerate system.
We analyze several data sets obtained by Hinode/EIS and find various types of flows during CMEs and EUV jet eruptions. CME-induced dimming regions are found to be characterized by significant blueshift and enhanced line width by using a single Gaussian fit. While a red-blue (RB) asymmetry analysis and a RB-guided double Gaussian fit of the coronal line profiles indicate that these are likely caused by the superposition of a strong background emission component and a relatively weak (~10%) high-speed (~100 km s-1) upflow component. This finding suggests that the outflow velocity in the dimming region is probably of the order of 100 km s-1, not ~20 km s-1 as reported previously. Density and temperature diagnostics suggest that dimming is primarily an effect of density decrease rather than temperature change. The mass losses in dimming regions as estimated from different methods are roughly consistent with each other and they are 20%-60% of the masses of the associated CMEs. With the guide of RB asymmetry analysis, we also find several temperature-dependent outflows (speed increases with temperature) immediately outside the (deepest) dimming region. In an erupted CME loop and an EUV jet, profiles of emission lines formed at coronal and transition region temperatures are found to exhibit two well-separated components, an almost stationary component accounting for the background emission and a highly blueshifted (~200 km s-1) component representing emission from the erupting material. The two components can easily be decomposed through a double Gaussian fit and we can diagnose the electron density, temperature and mass of the ejecta. Combining the speed of the blueshifted component and the projected speed of the erupting material derived from simultaneous imaging observations, we can calculate the real speed of the ejecta.
We present a precise photometric calibration of the first 1.5 years of science imaging from the Pan-STARRS1 survey (PS1), an ongoing optical survey of the entire sky north of declination -30 degrees in five bands. Building on the techniques employed by Padmanabhan et al. (2008) in the Sloan Digital Sky Survey (SDSS), we use repeat PS1 observations of stars to perform the relative calibration of PS1 in each of its five bands, solving simultaneously for the system throughput, the atmospheric transparency, and the large-scale detector flat field. Both internal consistency tests and comparison against the SDSS indicate that we achieve relative precision of <10 mmag in g, r, and i_P1, and ~10 mmag in z and y_P1. The spatial structure of the differences with the SDSS indicates that errors in both the PS1 and SDSS photometric calibration contribute similarly to the differences. The analysis suggests that both the PS1 system and the Haleakala site will enable <1% photometry over much of the sky.
We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)
We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.
We demonstrate a unified modeling method using the spectral synthesis code CLOUDY to self-consistently compute the structure of an H II and PDR that are in pressure equilibrium. With the unified model, emission contributions from PDRs and H II regions are computed simultaneously. We explore the effect of varying the spectral energy distribution of incident continuum, and apply the results to extra-galaxies. We focus on PDR diagnostics that also emerge from H II regions, and far-infrared (FIR) thermal emission from dust grains. We first present plots of FIR continuum ratios and the contribution from H II regions for dust FIR emission. With these plots, we successfully explain the FIR continuum ratios of M 82 observed by Herschel.
We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S_n) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of 15Msun star and compare them with the results in the latter computation. We also demonstrate multi-D capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.
The study of the morphology of galaxies is important in order to understand the formation and evolution of galaxies and their sub-components as a function of luminosity, environment, and star-formation and galaxy assembly over cosmic time. Disentangling the many variables that affect galaxy evolution and morphology, requires large galaxy samples and automated ways to measure morphology. The advent of large digital sky surveys, with unprecedented depth and resolution, coupled with sophisticated quantitative methods for morphology measurement are providing new insights in this fast evolving field of astronomical research.
Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.
One of the methods of detection and analysis of solar flares is observing the time variations of certain solar spectral lines. During solar flares, a raise of electron concentration occurs in Earth's ionosphere which results in amplitude and phase variations of the recorded very low frequency (VLF) waves. We compared the data obtained by the analysis of recorded VLF signals and line spectra for different solar flares. In this paper we treated the DHO VLF signal transmitted from Germany at the frequency of 23.4 kHz recorded by the AWESOME system in Belgrade (Serbia) during solar flares in the period between 10:40 UT and 13:00 UT on 2011 April 22.
Phase-targeted RXTE observations have allowed us to detect a transient 71.49 \pm 0.02 s signal that is most likely to be originating from the supergiant fast X-ray transient IGR J17544-2619. The phase-folded light curve shows a possible double-peaked structure with a pulsed flux of ~4.8*10^-12 erg cm^-2 s^-1 (3-10 keV). Assuming the signal to indicate the spin period of the neutron star in the system, the provisional location of IGR J17544-2619 on the Corbet diagram places the system within the classical wind-fed supergiant XRB region. Such a result illustrates the growing trend of supergiant fast X-ray transients to span across both of the original classes of HMXB in Porb - Pspin space.
Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the disc's radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the disc's longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($\pi$). The disc's radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.
MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov Telescope) is a system of two 17 meters Cherenkov telescopes, sensitive to very high energy (VHE; $> 10^{11}$ eV) gamma radiation above an energy threshold of 50 GeV. The first telescope was built in 2004 and operated for five years in stand-alone mode. A second MAGIC telescope (MAGIC-II), at a distance of 85 meters from the first one, started taking data in July 2009. Together they integrate the MAGIC stereoscopic system. Stereoscopic observations have improved the MAGIC sensitivity and its performance in terms of spectral and angular resolution, especially at low energies. We report on the status of the telescope system and highlight selected recent results from observations of galactic and extragalactic gamma-ray sources. The variety of sources discussed includes pulsars, galactic binary systems, clusters of galaxies, radio galaxies, quasars, BL Lacertae objects and more.
We present preliminary results of the long term spectral monitoring of two active galactic nuclei with different broad line shapes: Ark 564 and Arp 102B. Ark 564 is a bright nearby narrow line Syfert 1 (NLS1) galaxy with relatively narrow permitted optical emission lines and a high FeII/H${\beta}$ ratio, while Arp 102B is a nearby broad-line radio galaxy with broad double-peaked Balmer emission lines. The spectra of Ark 564 were observed during 11-year period (1999-2009) and the spectra of Arp 102B in the 12-year period (1998-2009), with SAO 6-m and 1-m telescopes (Russia) and the GHAO 2.1-m telescope (Cananea, Mexico).
We demonstrate that habitable Earth-type planets and moons can exist in the Kepler-16 system by investigating their orbital stability in the standard and extended habitable zone (HZ). We find that Earth-type planets in S-type orbits are possible within the standard HZ in direct vicinity of Kepler-16b, thus constituting habitable exomoons. However, Earth-mass planets cannot exist in P-type orbits around the two stellar components within the standard HZ. Yet, P-type Earth-mass planets can exist superior to the giant planet in the extended HZ pertaining to considerably enhanced back-warming in the planetary atmosphere if facilitated. We briefly discuss the potential detectability of such habitable Earth-type moons and planets positioned in S-type and P-type orbits, respectively.
We study the far infrared (60-500 $\mu$m) colours of late-type galaxies in the $Herschel$ Reference Survey, a K-band selected, volume limited sample of nearby galaxies. The far infrared colours are correlated with each other, with tighter correlations for the indices that are closer in wavelength. We also compare the different colour indices to various tracers of the physical properties of the target galaxies, such as the surface brightness of the ionising and non-ionising stellar radiation, the dust attenuation and the metallicity. The emission properties of the cold dust dominating the far infrared spectral domain are regulated by the properties of the interstellar radiation field. Consistent with that observed in nearby, resolved galaxies, our analysis shows that the ionising and the non-ionising stellar radiation, including that emitted by the most evolved, cold stars, both contribute to the heating of the cold dust component. This work also shows that metallicity is another key parameter characterising the cold dust emission of normal, late-type galaxies. A single modified black body with a grain emissivity index $\beta$=1.5 better fits the observed SPIRE flux density ratios $S250/S350$ vs. $S350/S500$ than $\beta$=2, although values of $\beta$ $\simeq$ 2 are possible in metal rich, high surface brightness galaxies. Values of $\beta$ $\lesssim$ 1.5 better represent metal poor, low surface brightness objects. This observational evidence provides strong constraints for dust emission models of normal, late type galaxies.
We provide a detailed theoretical study aimed at the observational finding about the nu Octantis binary system that indicates the possible existence of a Jupiter-type planet in this system. If a prograde planetary orbit is assumed, it has earlier been argued that the planet, if existing, should be located outside the zone of orbital stability. However, a previous study by Eberle & Cuntz (2010) [ApJ 721, L168] concludes that the planet is most likely stable if assumed to be in a retrograde orbit with respect to the secondary system component. In the present work, we significantly augment this study by taking into account the observationally deduced uncertainty ranges of the orbital parameters for the stellar components and the suggested planet. Furthermore, our study employs additional mathematical methods, which include monitoring the Jacobi constant, the zero velocity function, and the maximum Lyapunov exponent. We again find that the suggested planet is indeed possible if assumed to be in a retrograde orbit, but it is virtually impossible if assumed in a prograde orbit. Its existence is found to be consistent with the deduced system parameters of the binary components and of the suggested planet, including the associated uncertainty bars given by observations.
The Tayler instability (TI) is experimentally realized in a liquid-metal flow confined in a columnar container with an insulating outer cylinder. To predict the critical electrical current and the expected growth rates, simulations with MHD codes are used for a container with an inner cylinder whose radius is small or even zero. The very small magnetic Prandtl number of the gallium alloy (Pm\simeq 10^{-6}) only influences the growth rates rather than the critical field amplitudes. It is thus allowed to calculate the critical Hartmann numbers for marginal instability also with direct numerical simulations. The theoretical value of the critical electric current of 2.8 kA, resulting from both linear theory and simulations, is well confirmed by the experiment. Also the predicted (small) growth rates of the nonaxisymmetric kink-type perturbations are certified by the observed data. Due to the rather long growth times of order of minutes, the resulting Joule heating might excite convective modes forming a saturation mechanism of the observed TI.
Real-time seeing estimation at the focus of a telescope is nowadays strongly emphasized as this knowledge virtually drives the dimensioning of adaptive optics systems and instrument operational aspects. In this context we study the interest of using active optics Shack-Hartmann (AOSH) sensor images to provide accurate estimate of the seeing. The AOSH practically delivers long exposure spot PSFs -- at the critical location of the telescope focus -- being directly related to the atmospheric seeing in the line of sight. Although AOSH sensors are not specified to measure spot sizes but slopes, we show that accurate seeing estimation from AOSH images can be obtained with a dedicated algorithm. The sensitivity and comparison of two algorithms to various parameters is analyzed in a systematic way, demonstrating that efficient estimation of the seeing can be obtained by adequate means.
Fringing in CCD images is troublesome from the aspect of photometric quality and image flatness in the final reduced product. Additionally, defringing during calibration requires the inefficient use of time during the night to collect and produce a "supersky" fringe frame. The fringe pattern observed in a CCD image for a given near-IR filter is dominated by small thickness variations across the detector with a second order effect caused by the wavelength extent of the emission lines within the bandpass which produce the interference pattern. We show that essentially any set of emission lines which generally match the wavelength coverage of the night sky emission lines within a bandpass will produce an identical fringe pattern. We present an easy, inexpensive, and efficient method which uses a neon lamp as a flat field source and produces high S/N fringe frames to use for defringing an image during the calibration process.
We present the ORIGAMI method of identifying structures, particularly halos, in cosmological N-body simulations. Structure formation can be thought of as the folding of an initially flat three-dimensional manifold in six-dimensional phase space. ORIGAMI finds the outer folds that delineate these structures. Halo particles are identified as those that have undergone shell-crossing along 3 orthogonal axes, providing a dynamical definition of halo regions that is independent of density. ORIGAMI also identifies other morphological structures: particles that have undergone shell-crossing along 2, 1, or 0 orthogonal axes correspond to filaments, walls, and voids respectively. We compare this method to a standard Friends-of-Friends halo-finding algorithm and find that ORIGAMI halos are somewhat larger, more diffuse, and less spherical, though the global properties of ORIGAMI halos are in good agreement with other modern halo-finding algorithms.
We review our knowledge of the mixing properties of magnetic OB stars and discuss whether the observational data presently available support, as predicted by some theoretical models, the idea that magnetic phenomena favour the transport of the chemical elements. A (likely statistical) relationship between enhanced mixing and the existence of a field has been emerging over the last few years. As discussed in this contribution, however, a clear answer to this question is presently hampered by the lack of large and well-defined samples of magnetic and non-magnetic stars.
The mini-spiral is a feature of the interstellar medium in the central ~2 pc of the Galactic center. It is composed of several streamers of dust and ionised and atomic gas with temperatures between a few 100 K to 10^4 K. There is evidence that these streamers are related to the so-called circumnuclear disk of molecular gas and are ionized by photons from massive, hot stars in the central parsec. We attempt to constrain the emission mechanisms and physical properties of the ionized gas and dust of the mini-spiral region with the help of our multiwavelength data sets. Our observations were carried out at 1.3 mm and 3 mm with the mm interferometric array CARMA in California in March and April 2009, with the MIR instrument VISIR at ESO's VLT in June 2006, and the NIR Br-gamma with VLT NACO in August 2009. We present high resolution maps of the mini-spiral, and obtain a spectral index of 0.5 for Sgr A*, indicating an inverted synchrotron spectrum. We find electron densities within the range 0.8-1.5x10^4 cm-3 for the mini-spiral from the radio continuum maps, along with a dust mass contribution of ~0.25 solar masses from the MIR dust continuum, and extinctions ranging from 1.8-3 at 2.16 micron in the Br-gamma line. We observe a mixture of negative and positive spectral indices in our 1.3 mm and 3 mm observations of the extended emission of the mini-spiral, which we interpret as evidence that there are a range of contributions to the thermal free-free emission by the ionized gas emission and by dust at 1.3 mm.
[Abridged] We present a new grid of massive population III star models including the effects of rotation on the stellar structure and chemical mixing, and magnetic torques for the transport of angular momentum. Based on the grid, we also present a phase diagram for the expected final fates of rotating massive Pop III stars. Our non-rotating models become redder than the previous models in the literature, given the larger overshooting parameter adopted in this study. In particular, convective dredge-up of the helium core material into the hydrogen envelope is observed in our non-rotating very massive star models (>~200 Msun), which is potentially important for the chemical yields. On the other hand, the stars become bluer and more luminous with a higher rotational velocity. With the Spruit-Tayler dynamo, our models with a sufficiently high initial rotational velocity can reach the critical rotation earlier and lose more mass as a result, compared to the previous models without magnetic fields. The most dramatic effect of rotation is found with the so-called chemically homogeneous evolution (CHE), which is observed for a limited mass and rotational velocity range. CHE has several important consequences: 1) Both primary nitrogen and ionizing photons are abundantly produced. 2) Conditions for gamma-ray burst progenitors are fulfilled for an initial mass range of 13 - 84 Msun. 3) Pair instability supernovae of type Ibc are expected for 84 -190 Msun and 4) Both a pulsational pair instability supernova and a GRB may occur from the same progenitor of about 56 - 84 Msun, which might significantly influence the consequent GRB afterglow. We find that CHE does not occur for very massive stars (> 190 Msun), in which case the hydrogen envelope expands to the red-supergiant phase and the final angular momentum is too low to make any explosive event powered by rotation.
In the quest for the formation and evolution of galaxy clusters, Rakos and co-workers introduced a spectrophotometric method using the modified Str\"omgren photometry. But with the considerable debate toward the project's abilities, we re-introduce the system after a thorough testing of repeatability of colors and reproducibility of the ages and metallicities for six common galaxies in the three A779 data sets. A fair agreement has been found between the modified Str\"omgren and Str\"omgren filter systems to produce similar colors (with the precision of 0.09 mag in (uz-vz), 0.02 mag in (bz-yz), and 0.03 mag in (vz-vz)), ages and metallicities (with the uncertainty of 0.36 Gyr and 0.04 dex from the PCA and 0.44 Gyr and 0.2 dex using the GALEV models). We infer that the technique is able to relieve the age-metallicity degeneracy by separating the age effects from the metallicity effects, but still unable to completely break. We further extend this paper to re-study the evolution of galaxies in the low mass, dynamically poor A779 cluster by correlating the luminosity (mass), density, radial distance with the estimated age, metallicity, and the star formation history. Our results distinctly show the bimodality of the young, low-mass, metal-poor population with the mean age of 6.7 Gyr (\pm 0.5 Gyr) and the old, high-mass, metal-rich galaxies with the mean age of 9 Gyr (\pm 0.5 Gyr). The method also observes the color evolution of the blue cluster galaxies to red, and the downsizing phenomenon. Our analysis shows that the modified Str\"omgren photometry is very well suited for studying low- and intermediate-z clusters, as it is capable of observing deeper with better spatial resolution at spectroscopic redshift limits, and the narrowband filters estimate the age and metallicity with lesser uncertainties compared to other methods that study stellar population scenarios.
Testing the cosmic distance duality relation (CDDR) constitutes an important task for cosmology and fundamental physics since any violation of it would be a clear evidence of new physics. In this {\it Letter}, we propose a new test for the CDDR using only measurements of the gas mass fraction of galaxy clusters from Sunyaev-Zeldovich ($f_{SZE}$) and X-ray surface brightness ($f_{X-ray}$) observations. We show that the relation between current $f_{X-ray}$ and $f_{SZE}$ observations is given by $f_{SZE}=\eta f_{X-ray}$, where $\eta$ quantifies deviations from the CDDR. Since this latter expression is valid for the same object in a given galaxy cluster sample, the method proposed removes possible contaminations from different systematics error sources and redshift differences involved in luminosity and angular diameter distance measurements. We apply this cosmological model-independent methodology to the most recent $f_{X-ray}$ and $f_{SZE}$ data and show that no significant violation of the CDDR is found.
That are many factors that contribute to breaking the spherical symmetry of a neutron star. Most notably amongst those is rotation, magnetic field and/or accretion. Such phenomena are known to influence not only the macroscopic structure of the star, but also its microscopic composition. The breaking of symmetry in the microscopic and macroscopic realms might lead to anisotropic heat transport, which in turn might alter the thermal evolution of the object. The purpose of this paper is to investigate the cooling of rapidly rotating neutron stars, with the ultimate goal of understanding the effects of a 2D structure on thermal evolution. The equations that govern the thermal evolution of rotating neutron stars are presented, and the cooling of objects with different frequencies is computed numerically. We show that rotation can significantly influence the thermal evolution of neutron stars. Amongst the major modifications introduced by the 2D structure is the appearance of a hot spot on the poles, and an increase of the thermal coupling time of the core and the crust of the star. We also show that this increase is independent of the microscopic properties of the core, and depends only on the frequency of the object.
We present optical light curves of twenty-nine novae in M31 during the 2009 and 2010 observing seasons of the Palomar Transient Factory (PTF). The dynamic and rapid cadences in PTF monitoring of M31, from one day to even ten minutes, provide excellent temporal coverage of nova light curves, enabling us to record the photometric evolution of M31 novae in unprecedented detail. We also detect eight of these novae in the near ultraviolet (UV) band with the Galaxy Evolution Explorer (GALEX) satellite. Novae M31N2009-10b and 2010-11a show prominent UV emission peaking a few days prior to their optical maxima, possibly implying aspherical outbursts. Additionally, our blue-shifted spectrum of the recent outburst of PT And (M31N2010-12a) indicates that it is a recurrent nova in M31 and not a dwarf nova in the Milky Way as was previously assumed. Finally, we systematically searched for novae in all confirmed globular clusters of M31 and found only M31N 2010-10f associated with Bol 126. The specific nova rate in the M31 globular cluster system is thus about one per year which is not enhanced relative to the rate outside the globular cluster system.
A sample of 18286 radio-loud AGN is presented, constructed by combining the SDSS DR7 with the NVSS and FIRST radio surveys. Using this sample, the differences between `high-excitation' (or `quasar-mode'; HERG) and `low-excitation' (`radio-mode'; LERG) radio galaxies are investigated. A primary difference is the distinct nature of the Eddington-scaled accretion rate onto their central black holes: HERGs typically have accretion rates between 1 and 10% of Eddington, whereas LERGs predominatly accrete at a rate below 1% Eddington. This is consistent with models where the population dichotomy is caused by a switch between radiatively efficient and inefficient accretion modes at low accretion rates. Local radio luminosity functions are derived separately for the two populations, showing that although LERGs dominate at low luminosity and HERGs above 1e26 W/Hz, examples of both classes are found at all radio luminosities. Using the V/Vmax test it is shown that the populations show differential cosmic evolution at fixed radio luminosity: HERGs evolve strongly at all luminosities, while LERGs show weak or no evolution. This suggests that the luminosity-dependent evolution of the radio luminosity function is driven, at least in part, by the changing relative contributions of these two populations with luminosity. The host galaxies of the sources are also distinct: HERGs are typically of lower stellar mass, with lower black hole masses, bluer colours and weaker 4000-Ang breaks indicating younger stellar populations. These results offer strong support to the picture in which HERGs are fuelled at high rates through radiative accretion disks by cold gas, perhaps from mergers and interactions, while LERGs are fuelled via radiatively inefficient flows at low accretion rates, often by gas associated with the hot X-ray haloes of their host galaxy/cluster, as part of a radio-AGN feedback loop (abridged).
DARWIN (dark matter wimp search with noble liquids) is a design study for a next-generation, multi-ton dark matter detector in Europe. Liquid argon and/or liquid xenon are the target media for the direct detection of dark matter candidates in the form of weakly interacting massive particles (WIMPs). Light and charge signals created by particle interactions in the active detector volume are observed via the time projection chamber technique. DARWIN is to probe the spin-independent, WIMP-nucleon cross section down 1e-48 cm2 and to measure WIMP-induced nuclear recoil spectra with high-statistics, should they be discovered by an existing or near-future experiment. After a brief introduction, I will describe the project, selected R&D topics, expected backgrounds and the physics reach.
The graphical representation of the equation of state of a relativistic fermion system with a four-fermion interaction in the strong coupling regime is obtained as a function of the four-fermion coupling constant. It is shown, by increasing the coupling constant strength, how the crossover from a superconducting BCS regime to a Bose-Einstein-condensate (BEC) regime is manifested in the nature of the quasiparticles' energy spectrum. It is found that when the system has a distinguishable BEC nature its pressure becomes negative. We discuss the implications for astrophysics of the unstable BEC regime with associated negative pressure.
Following on after three previous papers discussing the formation of primordial black holes during the radiation-dominated era of the early universe, we present here a further investigation of the critical nature of the collapse. In particular, we focus on the long-lived intermediate state, which appears in collapses of perturbations close to the critical limit, and examine the extent to which this follows a similarity solution, as seen for critical collapse under more idealized circumstances (rather than within the context of an expanding universe, as studied here). We find that a similarity solution is indeed realised, to good approximation, for a region contained within the past light-cone of the forming black hole (and eventual singularity). The self-similarity is not exact, however, and this is explained by the presence within the light-cone of some outer matter still coupled to the expanding universe, which does not participate in the self-similarity. Our main interest, from a cosmological point of view, is in a radiative fluid with equation of state parameter $w=1/3$ (when the pressure $p$ and energy density $e$ are taken to be related by $p = we$). Other values of $w$, in the range $0 - 1$, have also been considered in the literature on critical collapse and we have looked at some of these too, within the context of our calculations, with the aim of gaining further insight into our main case of interest. As expected, we find that the features of scaling-law behaviour, intermediate state and similarity solution are preserved in all of the cases studied but with some interesting variations in the details. As in our previous work, we have started our simulations with initial supra-horizon scale perturbations of a type which could have come from inflation.
Quantum fields written on noncommutative spacetime (Groenewold - Moyal plane) obey twisted commutation relations. In this paper we show that these twisted commutation relations result in Hanbury-Brown Twiss (HBT) correlations that are distinct from that for ordinary bosonic or fermionic fields, and hence can provide us useful information about underlying noncommutative nature of spacetime. The deviation from usual bosonic/fermionic statistics becomes pronounced at high energies, suggesting that a natural place is to look at Ultra High Energy Cosmic Rays (UHECRs). Since the HBT correlations are sensitive only to the statistics of the particles, observations done with UHECRs are capable of providing unambiguous signatures of noncommutativity, without any detailed knowledge of the mechanism and source of origin of UHECRs.
Links to: arXiv, form interface, find, astro-ph, recent, 1201, contact, help (Access key information)
We report the discovery of ten proplyd-like objects in the vicinity of the massive OB association Cygnus OB2. They were discovered in IPHAS H-Alpha images and are clearly resolved in broad-band HST/ACS, near-IR and Spitzer mid-IR images. All exhibit the familiar tadpole shape seen in photoevaporating objects such as the Orion proplyds, with a bright ionization front at the head facing the central cluster of massive stars, and a tail stretching in the opposite direction. Many also show secondary ionization fronts, complex tail morphologies or multiple heads. We consider the evidence that these are either proplyds or `evaporating gaseous globules' (EGGs) left over from a fragmenting molecular cloud, but find that neither scenario fully explains the observations. Typical sizes are 50,000--100,000 AU, larger than the Orion proplyds, but in agreement with the theoretical scaling of proplyd size with distance from the ionizing source. These objects are located at projected separations of 6-14pc from the OB association, compared to 0.1pc for the Orion proplyds, but are clearly being photoionized by the 65 O-type stars in Cyg OB2. Central star candidates are identified in near- and mid-IR images, supporting the proplyd scenario, though their large sizes and notable asymmetries is more consistent with the EGG scenario. A third possibility is therefore considered, that these are a unique class of photoevaporating partially-embedded young stellar objects that have survived the destruction of their natal molecular cloud. This has implications for the properties of stars that form in the vicinity of massive stars.
UV and optical surveys are essential to gain insight into the processes driving galaxy formation and evolution. The rest-frame UV emission is key to measure the cosmic SFR. However, UV light is strongly reddened by dust. In starburst galaxies, the UV colour and the attenuation are linked, allowing to correct for dust extinction. Unfortunately, evidence has been accumulating that the relation between UV colour and attenuation is different for normal star-forming galaxies when compared to starburst galaxies. It is still not understood why star-forming galaxies deviate from the UV colour-attenuation relation of starburst galaxies. Previous work and models hint that the role of the shape of the attenuation curve and the age of stellar populations have an important role. In this paper we aim at understanding the fundamental reasons to explain this deviation. We have used the CIGALE SED fitting code to model the far UV to the far IR emission of a set of 7 reasonably face-on spiral galaxies from the HRS. We have explored the influence of a wide range of physical parameters to quantify their influence and impact on the accurate determination of the attenuation from the UV colour, and why normal galaxies do not follow the same relation as starburst galaxies. We have found that the deviation can be best explained by intrinsic UV colour differences between different regions in galaxies. Variations in the shape of the attenuation curve can also play a secondary role. Standard age estimators of the stellar populations prove to be poor predictors of the intrinsic UV colour. These results are also retrieved on a sample of 58 galaxies when considering their integrated fluxes. When correcting the emission of normal star-forming galaxies for the attenuation, it is crucial to take into account possible variations in the intrinsic UV colour as well as variations of the shape of the attenuation curve.
Despite their unique astrophysical relevance, the outcome of white dwarf binary mergers has so far only been studied for a very restricted number of systems. Here we present the results of a survey with more than two hundred simulations systematically scanning the white dwarf binary parameter space. We consider white dwarf masses ranging from 0.2 to 1.2 $M_\odot$ and account for their different chemical compositions. We find excellent agreement with the orbital evolution predicted by mass transfer stability analysis. Much of our effort in this paper is dedicated to determining which binary systems are prone to a thermonuclear explosion just prior to merger or at surface contact. We find that a large fraction of He-accreting binary systems explode: all dynamically unstable systems with accretor masses below 1.1 $M_\odot$ and donor masses above $\sim$ 0.4 $M_\odot$ are found to trigger a helium detonation at surface contact. A substantial fraction of these systems could explode at earlier times via detonations induced by instabilities in the accretion stream, as we have demonstrated in our previous work. We do not find definitive evidence for an explosion prior to merger or at surface contact in any of the studied double carbon-oxygen systems. Although we cannot exclude their occurrence if some helium is present, the available parameter space for a successful detonation in a white dwarf binary of pure carbon-oxygen composition is small. We demonstrate that a wide variety of dynamically unstable systems are viable type Ia candidates. The next decade thus holds enormous promise for the study of these events, in particular with the advent of wide-field synoptic surveys allowing a detailed characterization of their explosive properties.
We introduce a new technique that uses galaxy clustering to constrain how satellite galaxies lose stellar mass and contribute to the diffuse "intrahalo light" (IHL). We implement two models that relate satellite galaxy stellar mass loss to the detailed knowledge of subhalo dark matter mass loss. Model 1 assumes that the fractional stellar mass loss of a galaxy is proportional to the fractional amount of dark matter mass loss of its subhalo. Model 2 accounts for a delay in the time that stellar mass is lost since the galaxy resides deep in the potential well of the subhalo which may experience dark matter mass loss for some time before the galaxy is affected. We use these models to predict the stellar masses of a population of galaxies and use abundance matching to predict the clustering of several r-band luminosity threshold samples from the Sloan Digital Sky Survey. Abundance matching assuming no stellar mass loss (akin to abundance matching at the time of subhalo infall) over-estimates the correlation function on small scales (<~ 1 Mpc), while allowing too much stellar mass loss leads to an under-estimate. For each sample, we are thus able to constrain the amount of stellar mass loss required to match the observed clustering. We find that less luminous satellite galaxies experience more efficient stellar mass loss than luminous satellites. From these models, we can infer the amount of stellar mass that is deposited into the IHL. We find that both of our model predictions for the mean amount of IHL as a function of halo mass are consistent with current observational measurements. However, our two models predict a different amount of scatter in the IHL from halo to halo, with Model 2 being favored by observations. This demonstrates that a comparison to IHL measurements provides independent verification of our stellar mass loss models. (Abridged)
We have combined deep photometry in the B,V and I bands from CTIO/MOSAIC of the Sculptor dwarf spheroidal galaxy, going down to the oldest Main Sequence Turn-Offs, with spectroscopic metallicity distributions of Red Giant Branch stars. This allows us to obtain the most detailed and complete Star Formation History to date, as well as an accurate timescale for chemical enrichment. The Star Formation History shows that Sculptor is dominated by old ($>$10 Gyr), metal-poor stars, but that younger, more metal-rich populations are also present. Using Star Formation Histories determined at different radii from the centre we show that Sculptor formed stars with an increasing central concentration with time. The old, metal-poor populations are present at all radii, while more metal-rich, younger stars are more centrally concentrated. We find that within an elliptical radius of 1 degree, or 1.5 kpc from the centre, a total mass in stars of 7.8$\times10^{6}$ M$_{\odot}$ was formed, between 14 and 7 Gyr ago, with a peak at 13$-$14 Gyr ago. We use the detailed Star Formation History to determine age estimates for individual Red Giant Branch stars with high resolution spectroscopic abundances. Thus, for the first time, we can directly determine detailed timescales for the evolution of individual chemical elements. We find that the trends in alpha-elements match what is expected from an extended, relatively uninterrupted period of star formation continuing for 6$-$7 Gyr. The knee in the alpha-element distribution occurs at an age of 10.9$\pm$1Gyr, suggesting that SNe Ia enrichment began $\approx2\pm$1Gyr after the start of star formation in Sculptor.
We revisit the scaling relations and star-forming histories of local elliptical galaxies using a novel selection method applied to the Sloan Digital Sky Survey DR7. We combine two probability-based automated spectroscopic and morphological classifications of about 600000 galaxies with z<0.25 to isolate true elliptical galaxies. Our sample selection method does not introduce artificial cuts in the parameters describing the galaxy but instead it associates to every object a weight measuring the probability of being in a given spectro-morphological class. Thus the sample minimizes the selection biases. We show that morphologically defined ellipticals are basically distributed in 3 spectral classes, which dominate at different stellar masses. The bulk of the population (about 50%) is formed by a well defined class of galaxies with old stellar populations that formed their stars at very early epochs in a short episode of star formation. They dominate the scaling relations of elliptical galaxies known from previous works and represent the canonical elliptical class. At the low mass end, we find a population of slightly larger ellipticals, with smaller velocity dispersions at fixed stellar mass, which seem to have experienced a more recent episode of star formation probably triggered by gas-rich minor mergers. The high mass end tends to be dominated by a third spectral class, slightly more metal rich and with more efficient stellar formation than the reference class. This third class contributes to the curvature of the mass-size relation at high masses reported in previous works. Our method is therefore able to isolate typical spectra of elliptical galaxies following different evolutive pathways.
[Abridged] We make use of a semi-analytic cosmological model that includes simple prescriptions for dust attenuation and emission to make predictions for the observable and physical properties of galaxies that may be detected by the recently launched Herschel Space Observatory in deep fields such as GOODS-Herschel. We compare our predictions for differential galaxy number counts in the PACS (100 & 160) and SPIRE (250, 350, and 500 micron) bands with available observations. We find very good agreement with the counts in the PACS bands, for the overall counts and for galaxies binned by redshift at z< 2. At z > 2 our model underpredicts the number of bright galaxies by a factor of ten. The agreement is much worse for all three SPIRE bands, and becomes progressively worse with increasing wavelength. We discuss a number of possible reasons for these discrepancies, and hypothesize that the effect of blending on the observational flux estimates is likely to be the dominant issue. We note that the PACS number counts are relatively robust to changes in the dust emission templates, while the predicted SPIRE number counts are more template dependent. We present quantitative predictions for the relationship between the observed PACS 160 and SPIRE 250 micron fluxes and physical quantities such as halo mass, stellar mass, cold gas mass, star formation rate, and total infrared (IR) luminosity, at different redshifts. We also present quantitative predictions for the correlation between PACS 160 micron flux and the probability that a galaxy has experienced a recent major or minor merger. Although our models predict a strong correlation between these quantities, such that more IR-luminous galaxies are more likely to be merger-driven, we find that more than half of all high redshift IR-luminous galaxies detected by Herschel are able to attain their high star formation rates without enhancement by a merger.
We investigate the role of the environment on the colour and stellar population gradients in a local sample of ~3500 central and ~1150 satellite SDSS early-type galaxies (ETGs). The environment is parameterized in terms of the number of satellite galaxies, N_gal in each group. For central galaxies, we find that both optical colour and mass-to-light (M/L) ratio gradients are shallower in central galaxies residing in denser environments (higher N_gal). This trend is driven by metallicity gradients, while age gradients appear to be less dependent on the environment and to have a larger scatter. On the other hand, satellites do not show any differences in terms of the environment. The same results are found if galaxies are classified by central age, and both central and satellite galaxies have shallower gradients if they are older and steeper gradients if younger, satellites being independent of ages. In central galaxies, we show that the observed trends can be explained with the occurrence of dry mergings, which are more numerous in denser environments and producing shallower colour gradients because of more uniform metallicity distributions due to the mixing of stellar populations, while no final clues about merging occurrence can be obtained for satellites. Finally we discuss all systematics on stellar population fitting and their impact on the final results.
The Sun will eventually lose about half of its current mass nonlinearly over several phases of post-main sequence evolution. This mass loss will cause any surviving orbiting body to increase its semimajor axis and perhaps vary its eccentricity. Here, we use a range of Solar models spanning plausible evolutionary sequences and assume isotropic mass loss to assess the possibility of escape from the Solar System. We find that the critical semimajor axis in the Solar System within which an orbiting body is guaranteed to remain bound to the dying Sun due to perturbations from stellar mass loss alone is approximately 1,000 AU - 10,000 AU. The fate of objects near or beyond this critical semimajor axis, such as the Oort Cloud, outer scattered disc and specific bodies such as Sedna, will significantly depend on their locations along their orbits when the Sun turns off of the main sequence. These results are applicable to any exoplanetary system containing a single star with a mass, metallicity and age which are approximately equal to the Sun's, and suggest that few extrasolar Oort Clouds could survive post-main sequence evolution intact.
We performed an integrated optical polarization survey of 70 nearby galaxies to study the relationship between linear polarization and galaxy properties. To date this is the largest survey of its kind. The data were collected at McDonald Observatory using the Imaging Grism Polarimeter on the Otto Struve 2.1m telescope. Most of the galaxies did not have significant level of linear polarization, where the bulk is <1%. A fraction of the galaxies showed a loose correlation between the polarization and position angle of the galaxy, indicating that dust scattering is the main source of optical polarization. The unbarred spiral galaxies are consistent with the predicted relationship with inclination from scattering models of ~sin^2i.
Accretion of minor satellites has been postulated as the most likely mechanism to explain the significant size evolution of the massive galaxies over cosmic time. Using a sample of 629 massive (Mstar~10^11 Msun) galaxies from the near-infrared Palomar/DEEP-2 survey, we explore which fraction of these objects has satellites with 0.01 Msat < Mcentral < 1 (1:100) up to z=1 and which fraction has satellites with 0.1 Msat < Mcentral < 1 (1:10) up to z=2 within a projected radial distance of 100 kpc. We find that the fraction of massive galaxies with satellites, after the background correction, remains basically constant and close to ~30% for satellites with a mass ratio down to 1:100 up to z=1, and ~15% for satellites with a 1:10 mass ratio up to z=2. The family of spheroid-like massive galaxies presents a 2-3 times larger fraction of objects with satellites than the group of disk-like massive galaxies. A crude estimation of the number of 1:3 mergers a massive spheroid-like galaxy experiences since z~2 is around 2. For a disk-like galaxy this number decreases to ~1.
The acceleration of electrons at 1-10 keV energies is the cause of the polar aurora displays, and an important factor of magnetic energy transfer from the solar wind to the Earth. Two main families of acceleration processes are observed: those based on coherent quasi-static structures called double layers, and those based of the propagation of Alfv\'en Waves (AW). This paper is a review of the Alfv\'enic acceleration processes, and of their role in the global dynamics of the auroral zone.
We investigate the constraints that can be placed on the cosmic string tension by using the current Pulsar Timing Array limits on the stochastic gravitational wave background (SGWB). We have developed a code to compute the spectrum of gravitational waves (GWs) based on the widely accepted one-scale model. In its simplest form the one-scale model allows one to vary: (i) the string tension, G\mu/c^2; (ii) the size of cosmic string loops relative to the horizon at birth, \alpha; (iii) the spectral index of the emission spectrum, q; (iii) the cut-off in the emission spectrum, n_*; and (v) the intercommutation probability, p. The amplitude and slope of the spectrum in the nHz frequency range is very sensitive to these unknown parameters. We have also investigated the impact of more complicated scenarios with multiple initial loop sizes, in particular the 2-\alpha models proposed in the literature and a log-normal distribution for \alpha. We have computed the constraint on G\mu/c^2 due to the limit on a SGWB imposed by data from the European Pulsar Timing Array. Taking into account all the possible uncertainties in the parameters we find a conservative upper limit of G\mu/c^2<5.3x 10^{-7} which typically occurs when the loop production scale is close to the gravitational backreaction scale, \alpha\approx\Gamma G\mu/c^2. Stronger limits are possible for specific values of the parameters which typically correspond to the extremal cases \alpha\ll \Gamma G\mu/c^2 and \alpha\gg \Gamma G\mu/c^2. This limit is less stringent than the previously published limits which are based on cusp emission, an approach which does not necessarily model all the possible uncertainties. We discuss the prospects for lowering this limit by two orders of magnitude, or even a detection of the SGWB, in the very near future in the context of the Large European Array for Pulsars and the Square Kilometre Array.
Belonging to the group of B[e] stars, V921 Scorpii is associated with a strong infrared excess and permitted and forbidden line emission, indicating the presence of low- and high-density circumstellar gas and dust. Many aspects of V921 Sco and other B[e] stars still remain mysterious, including their evolutionary state and the physical conditions resulting in the class-defining characteristics. In this paper, we employ VLTI/AMBER spectro-interferometry in order to reconstruct high-resolution (lambda/2B=0.0013") model-independent interferometric images for three wavelength bands around 1.65, 2.0, and 2.3 micrometer. In our images, we discover a close (25.0+/-0.8 milliarcsecond, corresponding to 29+/-0.9 AU at 1.15 kpc) companion around V921 Sco. Between two epochs in 2008 and 2009, we measure orbital motion of 7 degrees, implying an orbital period of about 35 years (for a circular orbit). Around the primary star, we detect a disk-like structure with indications for a radial temperature gradient. The polar axis of this AU-scale disk is aligned with the arcminute-scale bipolar nebula in which V921 Sco is embedded. Using Magellan/IMACS imaging, we detect multi-layered arc-shaped sub-structure in the nebula, suggesting episodic outflow activity from the system with a period of about 25 years, roughly matching the estimated orbital period of the companion. Our study supports the hypothesis that the B[e] phenomenon is related to dynamical interaction in a close binary system.
We present a study of protoplanetary disks in spatially resolved low-mass binary stars in the well-known Orion Nebula Cluster (ONC) in order to assess the impact of binarity on the properties of circumstellar disks and its relation to the cluster environment. This is the currently largest such study in a clustered high stellar density star forming environment. We particularly aim at determining the presence of magnetospheric accretion and dust disks for each binary component, and at measuring the overall disk frequency. We carried out spatially resolved Adaptive Optics assisted near-IR photometry and spectroscopy of 26 binaries in the ONC, and determine stellar parameters such as effective temperatures and spectral types, luminosities, masses, as well as accretion properties and near-infrared excess for individual binary components. A fraction of 40(+10/-9)% of the binary components in the sample can be inferred to be T Tauri stars possesing an accretion disk. This is marginally lower than the disk fraction of single stars of ~50% in the ONC. We find that disks in wide binaries of >200AU separation are consistent with random pairing, while the evolution of circumprimary and circumsecondary disks is observed to be synchronized in closer binaries. Circumbinary disks appear to be not suited to explain this difference. Further, we identify several mixed pairs of accreting and non-accreting components, suggesting that these systems are common, and without preference for the more or less massive component to evolve faster. The derived mass accretion rates of the ONC binary components are of similar magnitude as those for ONC single stars and for binaries in the Taurus star forming region. The paper concludes with a discussion of the (presumably weak) connection between the presence of inner accretion disks in young binary systems and the existence of planets in stellar multiples.(abridged)
In this paper we explore the effect of decaying dark matter (DDM) on large-scale structure and possible constraints from galaxy imaging surveys. DDM models have been studied, in part, as a way to address apparent discrepancies between the predictions of standard cold dark matter models and observations of galactic structure. Our study is aimed at developing independent constraints on these models. In such models, DDM decays into a less massive, stable dark matter (SDM) particle and a significantly lighter particle. The small mass splitting between the parent DDM and the daughter SDM provides the SDM with a recoil or "kick" velocity vk, inducing a free-streaming suppression of matter fluctuations. This suppression may be probed via weak lensing power spectra measured by a number of forthcoming imaging surveys that aim primarily to constrain dark energy. Using scales on which linear perturbation theory alone is valid (multipoles < 300), surveys like Euclid or LSST can be sensitive to vk > 90 km/s for lifetimes ~ 1-5 Gyr. To estimate more aggressive constraints, we model nonlinear corrections to lensing power using a simple halo evolution model that is in good agreement with numerical simulations. In our most ambitious forecasts, using multipoles < 3000, we find that imaging surveys can be sensitive to vk ~ 10 km/s for lifetimes < 10 Gyr. Lensing will provide a particularly interesting complement to existing constraints in that they will probe the long lifetime regime far better than contemporary techniques. A caveat to these ambitious forecasts is that the evolution of perturbations on nonlinear scales will need to be well calibrated by numerical simulations before they can be realized. This work motivates the pursuit of such a numerical simulation campaign to constrain dark matter with cosmological weak lensing.
We discuss the theory and implementation of statistically rigorous fits to synchrotron self Compton models for datasets obtained from multi-wavelength observations of active galactic nuclei spectral energy distributions. The methods and techniques that we present are, then, exemplified reporting on a recent study of a nearby and well observed extragalactic source, Markarian 421.
The accelerating expansion of the universe is the most surprising cosmological discovery in many decades, implying that the universe is dominated by some form of "dark energy" with exotic physical properties, or that Einstein's theory of gravity breaks down on cosmological scales. The profound implications of cosmic acceleration have inspired ambitious experimental efforts to measure the history of expansion and growth of structure with percent-level precision or higher. We review in detail the four most well established methods for making such measurements: Type Ia supernovae, baryon acoustic oscillations (BAO), weak gravitational lensing, and galaxy clusters. We pay particular attention to the systematic uncertainties in these techniques and to strategies for controlling them at the level needed to exploit "Stage IV" dark energy facilities such as BigBOSS, LSST, Euclid, and WFIRST. We briefly review a number of other approaches including redshift-space distortions, the Alcock-Paczynski test, and direct measurements of H_0. We present extensive forecasts for constraints on the dark energy equation of state and parameterized deviations from GR, achievable with Stage III and Stage IV experimental programs that incorporate supernovae, BAO, weak lensing, and CMB data. We also show the level of precision required for other methods to provide constraints competitive with those of these fiducial programs. We emphasize the value of a balanced program that employs several of the most powerful methods in combination, both to cross-check systematic uncertainties and to take advantage of complementary information. Surveys to probe cosmic acceleration produce data sets with broad applications, and they continue the longstanding astronomical tradition of mapping the universe in ever greater detail over ever larger scales.
We present a survey of atomic hydrogen HI) emission in the direction of the Galactic Center conducted with the CSIRO Australia Telescope Compact Array (ATCA). The survey covers the area -5 deg < l < +5, -5 deg < b <+5 deg over the velocity range -309 < v_{LSR} < 349 km/s with a velocity resolution of 1 km/s. The ATCA data are supplemented with data from the Parkes Radio Telescope for sensitivity to all angular scales larger than the 145 arcsec angular resolution of the survey. The mean rms brightness temperature across the field is 0.7 K, except near (l,b)=(0 deg, 0 deg) where it increases to ~2 K. This survey complements the Southern Galactic Plane Survey to complete the continuous coverage of the inner Galactic plane in HI at ~2 arcmin resolution. Here we describe the observations and analysis of this Galactic Center survey and present the final data product. Features such as Bania's Clump 2, the far 3 kiloparsec arm and small high velocity clumps are briefly described.
In solar active regions (ARs), Doppler shifts measured along polarity inversion lines (PILs) of the line-of-sight (LOS) magnetic field determine one component of the velocity perpendicular to the magnetic field. Along PILs, these velocities can be used to: (i) improve estimates of photospheric electric fields, which can, in turn, be used to derive the Poynting flux of magnetic energy across the photosphere; and (ii) constrain the physical processes underlying flux cancellation, the mutual apparent loss of magnetic flux as closely spaced, opposite-polarity magnetogram features approach each other. Unfortunately, at least two factors introduce uncertainties into the zero point of measured Doppler velocities. First, instrumental effects (e.g., thermal variations in instrument components) can cause drifts in calibrations. Second, the convective blueshift, a well-known correlation between intensity and upflows, can bias estimates of the plasma's rest wavelength. Here, we present a method to absolutely calibrate LOS velocities using three successive vector magnetograms and one Dopplergram coincident with the central magnetogram. The method assumes ideal electric fields govern magnetic field evolution along PILs, and enforces consistency between changes in LOS flux near PILs and the transport of transverse magnetic flux by LOS velocities. For a subset of the initial vector magnetograms released by the SDO/HMI Team, we find clear evidence of offsets in the Doppler zero point, of 150 -- 300 m s$^{-1}$, which exhibit variations on the timescale of the SDO orbit. Estimation of the Doppler zero point from modeling the center-to-limb variation of convective blueshifts cannot account for such instrumental biases. Our method could also be used to identify episodes of flux cancellation or flux emergence in which non-ideal electric fields are present, and can characterize those fields.
In this paper we point out that redshift surveys can break the degeneracy between the galaxy bias, the power spectrum normalization, \sigma_{8,0} and the growth factor, without the need for external information by using a simple and rather general parameterization for the growth rate, the well known \gamma-parameterization. We find that in next-generation surveys like Euclid \sigma_{8,0} and \gamma can be measured to within 1%, 5%, respectively, while the bias b(z) can be measured to within 1-2% in every redshift bin.
The 3D structure of an active region (AR) filament is studied using nonlinear force-free field (NLFFF) extrapolations based on simultaneous observations at a photospheric and a chromospheric height. To that end, we used the Si I 10827 \AA\ line and the He I 10830 \AA\ triplet obtained with the Tenerife Infrared Polarimeter (TIP) at the VTT (Tenerife). The two extrapolations have been carried out independently from each other and their respective spatial domains overlap in a considerable height range. This opens up new possibilities for diagnostics in addition to the usual ones obtained through a single extrapolation from, typically, a photospheric layer. Among those possibilities, this method allows the determination of an average formation height of the He I 10830 \AA\ signal of \approx 2 Mm above the surface of the sun. It allows, as well, to cross-check the obtained 3D magnetic structures in view of verifying a possible deviation from the force- free condition especially at the photosphere. The extrapolations yield a filament formed by a twisted flux rope whose axis is located at about 1.4 Mm above the solar surface. The twisted field lines make slightly more than one turn along the filament within our box, which results in 0.055 turns/Mm. The convex part of the field lines (as seen from the solar surface) constitute dips where the plasma can naturally be supported. The obtained 3D magnetic structure of the filament depends on the choice of the observed horizontal magnetic field as determined from the 180\circ solution of the azimuth. We derive a method to check for the correctness of the selected 180\circ ambiguity solution.
We discuss selected large-scale anomalies in the maps of temperature anisotropies in the cosmic microwave background. Specfically, these include alignments of the largest modes of CMB anisotropy with one another and with the geometry and direction of motion of the Solar System, and the unexpected absence of two-point angular corellations especially outside the region of the sky most contaminated by the Galaxy. We discuss these findings in relation to expectations from standard inflationary cosmology. This paper is adapted from a talk given by one of us (GDS) at the SEENET-2011 meeting in August 2011 on the Serbian bank of the Danube River.
We present limits for the compactification scale in the theory of Large Extra Dimensions (LED) proposed by Arkani-Hamed, Dimopoulos, and Dvali. We use 11 months of data from the Fermi Large Area Telescope (Fermi-LAT) to set gamma ray flux limits for 6 gamma-ray faint neutron stars (NS). To set limits on LED we use the model of Hannestad and Raffelt (HR) that calculates the Kaluza-Klein (KK) graviton production in supernova cores and the large fraction subsequently gravitationally bound around the resulting NS. The predicted decay of the bound KK gravitons to {\gamma}{\gamma} should contribute to the flux from NSs. Considering 2 to 7 extra dimensions of the same size in the context of the HR model, we use Monte Carlo techniques to calculate the expected differential flux of gamma-rays arising from these KK gravitons, including the effects of the age of the NS, graviton orbit, and absorption of gamma-rays in the magnetosphere of the NS. We compare our Monte Carlo-based differential flux to the experimental differential flux using maximum likelihood techniques to obtain our limits on LED. Our limits are more restrictive than past EGRET-based optimistic limits that do not include these important corrections. Additionally, our limits are more stringent than LHC based limits for 3 or fewer LED, and comparable for 4 LED. We conclude that if the effective Planck scale is around a TeV, then for 2 or 3 LED the compactification topology must be more complicated than a torus.
We present the first results from our high-precision infrared (IR) astrometry program at the Canada-France-Hawaii Telescope. We measure parallaxes for 79 ultracool dwarfs (spectral types M6--T9) in 45 systems, with a median uncertainty of 1.3 mas (2.5%) while the best are 0.7 mas (1.0%). We provide the first parallaxes for 46 objects in 27 systems, and for another 25 objects in 15 systems, we significantly improve upon published results, with a median (best) improvement of 2.0x (5x). Three systems show astrometric perturbations indicative of orbital motion; two are known binaries (2MASS J0518-2828AB and 2MASS J1404-3159AB) and one is spectrally peculiar (SDSS J0805+4812). In addition, we present here a large set of Keck adaptive optics imaging that more than triples the number of binaries with L6--T5 components that have both multi-band photometry and distances. Our data enable an unprecedented look at the photometric properties of brown dwarfs as they cool through the L/T transition. Going from ~L8 to ~T4.5, flux in the Y and J bands increases by ~0.9 mag and ~0.7 mag, respectively (the Y- and J-band "bumps"), while the H-band flux increases by ~0.0--0.3 mag, and the K- and L'-band fluxes decline monotonically. This wavelength dependence is consistent with cloud clearing over a narrow range of temperature, since condensate opacity is expected to dominate at 1.0--1.3 micron. Interestingly, despite more than doubling the near-IR census of L/T transition objects, we find a conspicuous paucity of objects on the color--magnitude diagram just blueward of the late-L/early-T sequence. This "L/T gap" occurs at J-H = 0.1--0.3 mag, J-K = 0.0--0.4 mag, and implies that the last phases of cloud evolution occur rapidly. Finally, we provide a comprehensive update to the absolute magnitudes of ultracool dwarfs as a function of spectral type.
We report the discovery of a Type Ia supernova (SNIa) at redshift z=1.55 with the infrared detector of the Wide Field Camera 3 (WFC3-IR) on the Hubble Space Telescope (HST). This object was discovered in CANDELS imaging data of the Hubble Ultra Deep Field, and followed as part of the CANDELS+CLASH Supernova project, comprising the SN search components from those two HST multi-cycle treasury programs. This is the highest redshift SNIa with direct spectroscopic evidence for classification. It is also the first SN Ia at z>1 found and followed in the infrared, providing a full light curve in rest-frame optical bands. The classification and redshift are securely defined from a combination of multi-band and multi-epoch photometry of the SN, ground-based spectroscopy of the host galaxy, and WFC3-IR grism spectroscopy of both the SN and host. This object is the first of a projected sample at z>1.5 that will be discovered by the CANDELS and CLASH programs. The full CANDELS+CLASH SN Ia sample will enable unique tests for evolutionary effects that could arise due to differences in SN Ia progenitor systems as a function of redshift. This high-z sample will also allow measurement of the SN Ia rate out to z~2, providing a complementary constraint on SN Ia progenitor models.
Small levels of turbulence can be present in stellar radiative interiors due to, e.g., instability of rotational shear. In this paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification induces strong anisotropy with a very small ratio of radial-to-horizontal turbulence intensities. Angular momentum is transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. This non-diffusive transport known as the Lambda-effect has outward direction in radius and is much more efficient compared to the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed using direct numerical simulations combined with the test-scalar method. As a consequence of the non-diffusive transport of angular momentum, the estimated characteristic time of rotational coupling (< 100 Myr) between radiative core and convective envelope in young solar-type stars is much shorter compared to the time-scale of Lithium depletion (~1 Gyr).
The preliminary design of the new space gamma-ray telescope GAMMA-400 for the energy range 100 MeV - 3 TeV is presented. The angular resolution of the instrument, 1-2{\deg} at E{\gamma} ~100 MeV and ~0.01^{\circ} at E{\gamma} > 100 GeV, its energy resolution ~1% at E{\gamma} > 100 GeV, and the proton rejection factor ~10E6 are optimized to address a broad range of science topics, such as search for signatures of dark matter, studies of Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission, gamma-ray bursts, as well as high-precision measurements of spectra of cosmic-ray electrons, positrons, and nuclei.
The difference between vacuum energy of quantum fields in Minkowski space and in Friedmann-Robterson-Walker universe might be related to the observed dark energy. The vacuum energy of the Veneziano ghost field introduced to solve the $U(1)_A$ problem in QCD is of the form, $ H+ {\cal O}(H^2)$. Based on this, we study the dynamical evolution of a phenomenological dark energy model whose energy density is of the form $\alpha H+\beta H^2$. In this model, the universe approaches to a de Sitter phase at late times. We fit the model with current observational data including SnIa, BAO, CMB, BBN, Hubble parameter and growth rate of matter perturbation. It shows that the universe begins to accelerate at redshift $z\sim 0.75$ and this model is consistent with current data. In particular, this model fits the data of growth factor well as the $\Lambda CDM$ model.
We solve for the velocity fields of momentum-conserving supershells driven by feedback from supermassive black holes or nuclear star clusters (central massive objects: CMOs). We treat, for the first time, the case of CMOs embedded in gaseous protogalaxies with non-isothermal dark matter haloes having peaked circular-speed profiles. We find the CMO mass that is sufficient to drive any shell to escape any such halo. In the limit of large halo mass, relevant to real galaxies, this critical CMO mass depends only on the peak circular speed in the halo, scaling as M_crit \propto V_c,pk^4.
UGC 12281 has been classified as having a pure disk and being a low surface brightness galaxy (LSBG), thus being an obvious member of the so-called superthin galaxies. At the same time it represents an extremely untypical type of LSBG due to its remarkable amount of current star formation and evidence for extraplanar ionized gas. This makes it become a perfect tool to investigate the triggering of star formation in LSB galaxies, located in an alleged isolated area. By means of deep photometry and long-slit spectroscopy we analyse the H$\alpha$ halo and verify the existence of a potential dwarf companion which we found on processed SDSS images.
The integrated brightness of the Sun shows variability on time-scales from minutes to decades. This variability is mainly caused by pressure mode oscillations, by granulation and by dark spots and bright faculae on the surface of the Sun. By analyzing the frequency spectrum of the integrated brightness we can obtain greater knowledge about these phenomena. It is shown how the frequency spectrum of the integrated brightness of the Sun in the frequency range from 0.1 to 3.2 mHz shows clear signs of both granulation, faculae and p-mode oscillations and that the measured characteristic time-scales and amplitudes of the acoustic signals from granulation and faculae are consistent with high-resolution observations of the solar surface. Using 13 years of observations of the Sun's integrated brightness from the VIRGO instrument on the SOHO satellite it is shown that the significance of the facular component varies with time and that it has a significance above 0.99 around half the time. Furthermore, an analysis of the temporal variability in the measured amplitudes of both the granulation, faculae and p-mode oscillation components in the frequency spectrum reveals that the amplitude of the p-mode oscillation component shows variability that follows the solar cycles, while the amplitudes of the granulation and facular components show signs of quasi-annual and quasi-biennial variability, respectively.
In this letter we compute stringent astrophysical and cosmological constraints on a recently proposed Eddington-inspired Born-Infeld theory of gravity. We find, using a generalized version of the Zel'dovich approximation, that in this theory a pressureless cold dark matter fluid has a non-zero effective sound speed. We compute the corresponding effective Jeans length and show that it is approximately equal to the fundamental length of the theory $R_*=\kappa^{1/2} G^{-1/2}$, where $\kappa$ is the only additional parameter of theory with respect to general relativity and $G$ is the gravitational constant. This scale determines the minimum size of compact objects which are held together by gravity. We also estimate the critical mass above which pressureless compact objects are unstable to colapse into a black hole, showing that it is approximately equal to the fundamental mass $M_* = \kappa^{1/2} c^2 G^{-3/2}$, and we show that the maximum density attainable inside stable compact stars is roughly equal to the fundamental density $\rho_*=\kappa^{-1} c^2$, where $c$ is the speed of light in vacuum. We find that the mere existence of astrophysical objects of size $R$ which are held together by their own gravity leads to the constraint $\kappa < G R^2$. In the case of neutron stars this implies that $\kappa < 10^{-2} \, {\rm m^5 \, kg^{-1} \, s^{-2}}$, a limit which is stronger by about 10 orders of magnitude than big bang nucleosynthesis constraints and by more than 7 orders of magnitude than solar constraints.
Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.
The discovery of new objects in modern wide-field asteroid and comet surveys can be enhanced by first identifying observations belonging to known solar system objects. The assignation of new observations to a known object is an attribution problem that occurs when a least squares orbit already exists for the object but a separate fit is not possible to just the set of new observations. In this work we explore the strongly asymmetric attribution problem in which the existing least squares orbit is very well constrained and the new data are sparse. We describe an attribution algorithm that introduces new quality control metrics in the presence of strong biases in the astrometric residuals. The main biases arise from the stellar catalogs used in the reduction of asteroid observations and we show that a simple debiasing with measured regional catalog biases significantly improves the results. We tested the attribution algorithm using data from the PS1 survey that used the 2MASS star catalog for the astrometric reduction. We found small but statistically significant biases in the data of up to 0.1 arcsec that are relevant only when the observations reach the level of accuracy made possible by instruments like PS1. The false attribution rate was measured to be < 1/1000 with a simple additional condition that can reduce it to zero while the attribution efficiency is consistent with 100%.
We present a brief description of a model for the broad emission line region (BELR) in quasars, which is supported by analysis of CIV and other emission lines in the spectra of high-z SDSS quasars. Specifically we consider a two-component BELR with a disk and wind where the relative strength of each component is a function of luminosity. The implications of such a model for our understanding of quasar outflows and estimates of their black hole masses and accretion rates are discussed.
Studies are under way to propose a new generation of post-VLTI interferometers. The Carlina concept studied at the Haute- Provence Observatory is one of the proposed solutions. It consists in an optical interferometer configured like a diluted version of the Arecibo radio telescope: above the diluted primary mirror made of fixed cospherical segments, a helium balloon (or cables suspended between two mountains), carries a gondola containing the focal optics. Since 2003, we have been building a technical demonstrator of this diluted telescope. First fringes were obtained in May 2004 with two closely-spaced primary segments and a CCD on the focal gondola. We have been testing the whole optical train with three primary mirrors. The main aim of this article is to describe the metrology that we have conceived, and tested under the helium balloon to align the primary mirrors separate by 5-10 m on the ground with an accuracy of a few microns. The servo loop stabilizes the mirror of metrology under the helium balloon with an accuracy better than 5 mm while it moves horizontally by 30 cm in open loop by 10-20 km/h of wind. We have obtained the white fringes of metrology; i.e., the three mirrors are aligned (cospherized) with an accuracy of {\approx} 1 micron. We show data proving the stability of fringes over 15 minutes, therefore providing evidence that the mechanical parts are stabilized within a few microns. This is an important step that demonstrates the feasibility of building a diluted telescope using cables strained between cliffs or under a balloon. Carlina, like the MMT or LBT, could be one of the first members of a new class of telescopes named diluted telescopes.
The AGILE space mission (whose instrument is sensitive in the energy ranges 18-60 keV, and 30 MeV - 50 GeV) has been operating since 2007. Assessing the statistical significance of time variability of gamma-ray sources above 100 MeV is a primary task of the AGILE data analysis. In particular, it is important to check the instrument sensitivity in terms of Poisson modeling of the data background, and to determine the post-trial confidence of detections. The goals of this work are: (i) evaluating the distributions of the likelihood ratio test for "empty" fields, and for regions of the Galactic plane; (ii) calculating the probability of false detection over multiple time intervals. In this paper we describe in detail the techniques used to search for short-term variability in the AGILE gamma-ray source database. We describe the binned maximum likelihood method used for the analysis of AGILE data, and the numerical simulations that support the characterization of the statistical analysis. We apply our method to both Galactic and extra-galactic transients, and provide a few examples. After having checked the reliability of the statistical description tested with the real AGILE data, we obtain the distribution of p-values for blind and specific source searches. We apply our results to the determination of the post-trial statistical significance of detections of transient gamma-ray sources in terms of pre-trial values. The results of our analysis allow a precise determination of the post-trial significance of {\gamma}-ray sources detected by AGILE.
We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2 and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+ indicating it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine there is a probable 'dead zone' where accretion is suppressed, present in a layer, Z/R <~ 0.1 - 0.2, in the disk midplane, within R \approx 200 AU.
We performed a spatially resolved spectral X-ray study of the pulsar wind nebula (PWN) in the supernova remnant G0.9+0.1. Furthermore we modeled its nonthermal emission in the X-ray and very high energy (VHE, E > 100 GeV) gamma-ray regime. Using Chandra ACIS-S3 data, we investigated the east-west dependence of the spectral properties of G0.9+0.1 by calculating hardness ratios. We analyzed the EPIC-MOS and EPIC-pn data of two on-axis observations of the XMM-Newton telescope and extracted spectra of four annulus-shaped regions, centered on the region of brightest emission of the source. A radially symmetric leptonic model was applied in order to reproduce the observed X-ray emission of the inner part of the PWN. Using the optimized model parameter values obtained from the X-ray analysis, we then compared the modeled inverse Compton (IC) radiation with the published H.E.S.S. gamma-ray data. The spectral index within the four annuli increases with growing distance to the pulsar, whereas the surface brightness drops. With the adopted model we are able to reproduce the characteristics of the X-ray spectra. The model results for the VHE gamma radiation, however, strongly deviate from the H.E.S.S. data.
Recently there has been tremendous increase in the number of identified extra-solar planetary systems. Our understanding of their formation is tied to exoplanet internal structure models, which rely upon equations of state of light elements and compounds like water. Here we present shock compression data for water with unprecedented accuracy that shows water equations of state commonly used in planetary modeling significantly overestimate the compressibility at conditions relevant to planetary interiors. Furthermore, we show its behavior at these conditions, including reflectivity and isentropic response, is well described by a recent first-principles based equation of state. These findings advocate this water model be used as the standard for modeling Neptune, Uranus, and "hot Neptune" exoplanets, and should improve our understanding of these types of planets.
We study the non-thermal jet emission of the BL Lac object B3 2247+381 during a high optical state. The MAGIC telescopes observed the source during 13 nights between September 30th and October 30th 2010, collecting a total of 14.2 hours of good quality very high energy (VHE) $\gamma$-ray data. Simultaneous multiwavelength data was obtained with X-ray observations by the Swift satellite and optical R-band observations at the KVA-telescope. We also use high energy $\gamma$-ray (HE, 0.1 GeV-100 GeV) data from the Fermi satellite. The BL Lac object B3 2247+381 (z=0.119) was detected, for the first time, at VHE $\gamma$-rays at a statistical significance of 5.6 $\sigma$. A soft VHE spectrum with a photon index of -3.2 $\pm$ 0.6 was determined. No significant short term flux variations were found. We model the spectral energy distribution using a one-zone SSC-model, which can successfully describe our data.
The double pulsar PSR J0737-3039A/B displays short, 30 s eclipses that arise around conjunction when the radio waves emitted by pulsar A are absorbed as they propagate through the magnetosphere of its companion pulsar B. These eclipses offer a unique opportunity to probe directly the magnetospheric structure and the plasma properties of pulsar B. We have performed a comprehensive analysis of the eclipse phenomenology using multi-frequency radio observations obtained with the Green Bank Telescope. We have characterized the periodic flux modulations previously discovered at 820 MHz by McLaughlin et al., and investigated the radio frequency dependence of the duration and depth of the eclipses. Based on their weak radio frequency evolution, we conclude that the plasma in pulsar B's magnetosphere requires a large multiplicity factor (~ 10^5). We also found that, as expected, flux modulations are present at all radio frequencies in which eclipses can be detected. Their complex behavior is consistent with the confinement of the absorbing plasma in the dipolar magnetic field of pulsar B as suggested by Lyutikov & Thompson and such a geometric connection explains that the observed periodicity is harmonically related to pulsar B's spin frequency. We observe that the eclipses require a sharp transition region beyond which the plasma density drops off abruptly. Such a region defines a plasmasphere which would be well inside the magnetospheric boundary of an undisturbed pulsar. It is also two times smaller than the expected standoff radius calculated using the balance of the wind pressure from pulsar A and the nominally estimated magnetic pressure of pulsar B.
We present a detailed analysis of the global and fine structure of four middle-mass disc galaxies obtained from simulations in a $\Lambda$CDM scenario. These objects have photometric D/T ratios in good agreement with those observed for late-type spirals, as well as kinematic properties in agreement with the observational Tully-Fisher relation. We identify the different dynamical components at z=0 on the basis of both orbital parameters and the binding energy of stars in the galaxy. In this way, we recognize a slowly rotating centrally concentrated spheroid, and two disc components supported by rotation: a thin disc with stars in nearly circular orbits, and a thick disc with orbital parameters transitional between the thin disc and the spheroid. The spheroidal component is composed mainly by old, metal-poor and {\alpha}-enhanced stars. The distribution of metals in this component shows, however, a clear bimodality with a low-metallicity peak, which could be related to a classical bulge, and a high-metallicity peak, which could be related to a pseudo-bulge. The thin disc appears in our simulations as the youngest and most metal-rich component. The radial distribution of ages and colours in this component are U-shaped: the new stars are forming in the inner regions, and then migrate through secular processes. Finally, we also find a thick disc containing about 16% of the total stellar mass and with properties that are intermediate between those of the thin disc and the spheroid. Its low-metallicity stars are {\alpha}-enhanced when compared to thin disc stars of the same metallicity. The structural parameters (e.g., the scale height) of the simulated thick discs suggest that such a component could result from the combination of different thickening mechanisms that include merger-driven processes, but also long-lived internal perturbations of the thin disc. [Abridged]
We develop and apply an enhanced regularization algorithm, used in RHESSI X-ray spectral analysis, to constrain the ill-posed inverse problem that is determining the DEM from solar observations. We demonstrate this computationally fast technique applied to a range of DEM models simulating broadband imaging data from SDO/AIA and high resolution line spectra from Hinode/EIS, as well as actual active region observations with Hinode/EIS and XRT. As this regularization method naturally provides both vertical and horizontal (temperature resolution) error bars we are able to test the role of uncertainties in the data and response functions. The regularization method is able to successfully recover the DEM from simulated data of a variety of model DEMs (single Gaussian, multiple Gaussians and CHIANTI DEM models). It is able to do this, at best, to over four orders of magnitude in DEM space but typically over two orders of magnitude from peak emission. The combination of horizontal and vertical error bars and the regularized solution matrix allows us to easily determine the accuracy and robustness of the regularized DEM. We find that the typical range for the horizontal errors is $\Delta$log$T\approx 0.1 -0.5$ and this is dependent on the observed signal to noise, uncertainty in the response functions as well as the source model and temperature. With Hinode/EIS an uncertainty of 20% greatly broadens the regularized DEMs for both Gaussian and CHIANTI models although information about the underlying DEMs is still recoverable. When applied to real active region observations with Hinode/EIS and XRT the regularization method is able to recover a DEM similar to that found via a MCMC method but in considerably less computational time.
Supersonic turbulence is an essential element in understanding how structure within interstellar gas is created and shaped. In the context of star formation, many computational studies show that the mass spectrum of density and velocity fluctuations within dense clouds, as well as the distribution of their angular momenta, trace their origin to the statistical and physical properties of gas that is lashed with shock waves. In this article, we review the observations, simulations, and theories of how turbulent-like processes can account for structures we see in molecular clouds. We then compare traditional ideas of supersonic turbulence with a simpler physical model involving the effects of multiple shock waves and their interaction in the interstellar medium. Planar intersecting shock waves produce dense filaments, and generate vortex sheets that are essential to create the broad range of density and velocity structure in clouds. As an example, the lower mass behaviour of the stellar initial mass function can be traced to the tendency of a collection of shock waves to build-up a log-normal density distribution (or column density). Vorticity - which is essential to produce velocity structure over a very broad range of length scales in shocked clouds - can also be generated by the passage of curved shocks or intersecting planar shocks through such media. Two major additional physical forces affect the structure of star forming gas - gravity and feedback processes from young stars. Both of these can produce power-law tails at the high mass end of the IMF.
We report here on the discovery of stellar occultations, observed with Kepler, that recur periodically at 15.685 hour intervals, but which vary in depth from a maximum of 1.2% to a minimum that can be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a V = 16 magnitude K dwarf with T_eff = 4400 K. The out-of-occultation behavior shows no evidence for ellipsoidal light variations, indicating that the mass of the orbiting object is less than ~3 M_J. Because the eclipse depths are highly variable, they cannot be due solely to transits of a single planet with a fixed size. We discuss but dismiss a scenario involving a binary giant planet whose mutual orbit plane precesses, bringing one of the planets into and out of a grazing transit. This scenario seems ruled out by the dynamical instability that would result from such a configuration. The much more likely explanation involves macroscopic particles - e.g., dust, possibly in the form of micron-sized pyroxene grains - escaping the atmosphere of a slowly disintegrating planet not much larger than Mercury in size. The planetary surface is hot enough to sublimate; the resultant silicate vapor accelerates off the planet via a Parker-type thermal wind, dragging dust grains with it. We infer a mass loss rate from the observations of order ~1 M_earth/Gyr, with a dust-to-gas ratio possibly of order unity. For our fiducial 0.1 M_earth planet (twice the mass of Mercury), the evaporation timescale may be ~0.2 Gyr. Smaller mass planets are disfavored because they evaporate still more quickly, as are larger mass planets because they have surface gravities too strong to sustain outflows with the requisite mass-loss rates. The occultation profile evinces an ingress-egress asymmetry that could reflect a comet-like dust tail trailing the planet; we present simulations of such a tail.
We present the results of Rossi X-ray Timing Explorer (RXTE) and Swift monitoring observations of the magnetar 1E 1547.0-5408 following the pulsar's radiative outbursts in 2008 October and 2009 January. We report on a study of the evolution of the timing properties and the pulsed flux from 2008 October 4 through 2009 December 26. We show that the pulsed flux decrease which followed an initial rise in the 2008 outburst was interrupted by a spike ~9 days after the initial outburst. In our timing study, a phase-coherent analysis shows that for the first 29 days following the 2008 outburst, there was a very fast increase in the magnitude of the rotational frequency derivative nudot, such that the second derivative was a factor of ~60 larger than that reported in data from 2007. This nudot magnitude increase occurred in concert with the decay of the pulsed flux following the start of the 2008 event. Following the 2009 outburst, for the first 23 days, the second derivative was consistent with zero, and nudot had returned to close to its 2007 value. In contrast to the 2008 event, the 2009 outburst showed a major increase in persistent flux, relatively little change in the pulsed flux, and sudden significant spectral hardening ~15 days after the outburst. We show that, excluding the month following each of the outbursts, and because of the noise and the sparsity in the data, multiple plausible timing solutions fit the pulsar's frequency behavior. We note similarities in the behavior of 1E 1547.0-5408 following the 2008 outburst to that seen in the AXP 1E 1048.1-5937 following its 2001-2002 outburst and discuss this in terms of the magnetar model.
Context. Gravitational lensing is one of the leading tools in understanding the dark side of the Universe. The need for accurate, efficient and effective methods which are able to extract this information along with other cosmological parameters from cosmic shear data is ever growing. COSEBIs, Complete Orthogonal Sets of E-/B-Integrals, is a recently developed statistical measure that encompasses the complete E-/B-mode separable information contained in the shear correlation functions measured on a finite angular range. Aims. The aim of the present work is to test the properties of this newly developed statistics for a higher-dimensional parameter space and to generalize and test it for shear tomography. Methods. We use Fisher analysis to study the effectiveness of COSEBIs. We show our results in terms of figure-of-merit quantities, based on Fisher matrices. Results. We find that a relatively small number of COSEBIs modes is always enough to saturate to the maximum information level. This number is always smaller for 'logarithmic COSEBIs' than for 'linear COSEBIs', and also depends on the number of redshift bins, the number and choice of cosmological parameters, as well as the survey characteristics. Conclusions. COSEBIs provide a very compact way of analyzing cosmic shear data, i.e., all the E-/B-mode separable second-order statistical information in the data is reduced to a small number of COSEBIs modes. Furthermore, with this method the arbitrariness in data binning is no longer an issue since the COSEBIs modes are discrete. Finally, the small number of modes also implies that covariances, and their inverse, are much more conveniently obtainable, e.g., from numerical simulations, than for the shear correlation functions themselves.
We present preliminary results from an investigation into broad absorption line (BAL) variability within a sample of 41 radio-loud quasars (RLQs). Using 28 new Hobby-Eberly Telescope (HET) spectra along with earlier Sloan Digital Sky Survey (SDSS) or other archival data, we generate a total set of 50 pairs of BAL equivalent width measurements. Absorption variability in BAL RLQs typically consists of modest changes in the depth of trough segments, and variability is more common on longer rest-frame timescales; these tendencies are similar to previous findings for BAL radio-quiet quasars (RQQs). BAL variability in RLQs does not show any obvious dependence upon radio luminosity or loudness, but there is suggestive support for greater fractional variability within lobe-dominated RLQs.
We analyze the general relativistic oscillations of thin accretion disks around compact astrophysical objects interacting with the surrounding medium through non-gravitational forces. The interaction with the external medium (a thermal bath) is modeled via a friction force, and a random force, respectively. The general equations describing the stochastically perturbed disks are derived by considering the perturbations of trajectories of the test particles in equatorial orbits, assumed to move along the geodesic lines. By taking into account the presence of a viscous dissipation and of a stochastic force we show that the dynamics of the stochastically perturbed disks can be formulated in terms of a general relativistic Langevin equation. The stochastic energy transport equation is also obtained. The vertical oscillations of the disks in the Schwarzschild and Kerr geometries are considered in detail, and they are analyzed by numerically integrating the corresponding Langevin equations. The vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases.
From the viewpoint of no-cloning theorem we postulate a relation between the current accelerated expansion of our universe and the inflationary expansion in the very early universe. It implies that the fate of our universe should be in a state with accelerated expansion. Quantitatively we find that the no-cloning theorem leads to a lower bound on the cosmological constant which is compatible with observations.
Self-consistent account of the most simple non-gauge vector fields leads to a broad spectrum of regular scenarios of temporal evolution of the Universe completely within the frames of the Einstein's General relativity. The longitudinal non-gauge vector field is "the missing link in the chain", displaying the repulsive elasticity and allowing the macroscopic description of the main features of the Universe evolution. The singular Big Bang turns into a regular inflation-like state of maximum compression with the further accelerated expansion at late times. The parametric freedom of the theory allows to forget the troubles of fine tuning. In the most interesting cases the analytical solutions of the Einstein's equations are found.
It has been pointed out recently that the presence of dilaton field in the early Universe can dilute the neutralino dark matter (DM) abundance, if Universe is not radiation dominated at DM decoupling, due to its dissipative-like coupling to DM. In this scenario two basic mechanisms compete, the modified Hubble expansion rate tending to increase the relic density and a dissipative force that tends to decrease it. The net effect can lead to an overall dramatic decrease of the predicted relic abundance, sometimes by amounts of the order of O(10^2) or so. This feature is rather generic, independent of any particular assumption on the underlying string dynamics, provided dilaton dominates at early eras after the end of inflation but before Big Bang Nucleosynthesis (BBN). The latter ensures that BBN is not upset by the presence of the dilaton. In this paper, within the context of such a scenario, we study the phenomenology of the constrained minimal supersymmetric model (CMSSM) by taking into account all recent experimental constraints, including those from the LHC searches. We find that the allowed parameter space is greatly enlarged and includes regions that are beyond the reach of LHC. The allowed regions are compatible with Direct Dark Matter searches since the small neutralino annihilation rates, that are now in accord with the cosmological data on the relic density, imply small neutralino-nucleon cross sections below the sensitivities of the Direct Dark Matter experiments. It is also important that the new cosmologically accepted regions are compatible with Higgs boson masses larger than 120 GeV, as it is indicated from the LHC experimental data.
The observational relation between the density of baryon and dark matter in the Universe, $\Omega_{\rm DM}/\Omega_B\simeq 5$, is one of the most difficult problems to solve in modern cosmology. We discuss a scenario that explains this relation by combining the asymmetric dark matter scenario and the spontaneous baryogenesis associated with the flat direction in the supersymmetric standard model. A part of baryon asymmetry is transferred to charge asymmetry $D$ that dark matter carries, if a symmetry violating interaction that works at high temperature breaks not only $B-L$ but also $D$ symmetries simultaneously. In this case, the present number density of baryon and dark matter can be same order if the symmetric part of dark matter annihilates sufficiently. Moreover, the baryon number density can be enhanced as compared to that of dark matter if another $B-L$ violating interaction is still in thermal equilibrium after the spontaneous genesis of dark matter, which accommodates a TeV scale asymmetric dark matter model.
We argue that the local violation of P and CP invariance in heavy ion collisions and the universal thermal aspects observed in high energy collisions are in fact two sides of the same coin, and both are related to quantum anomalies of QCD. We argue that the low energy relations representing the quantum anomalies of QCD are saturated by coherent low dimensional vacuum configurations as observed in Monte Carlo lattice studies. The thermal spectrum and approximate universality of the temperature with no dependence on energy of colliding particles in this framework is due to the fact that the emission results from the distortion of these low dimensional vacuum sheets rather than from the colliding particles themselves. The emergence of the long- range correlations of P odd domains (a feature which is apparently required for explanation of the asymmetry observed at RHIC and LHC) is also a result of the same distortion of the QCD vacuum configurations. We formulate the corresponding physics using the effective low energy effective Lagrangian. We also formulate the same physics in terms of the dual holographic picture when low-dimensional sheets of topological charge embedded in 4d space, as observed in Monte Carlo simulations, are identified with D2 branes. Finally, we argue that study of these long range correlations in heavy ion collisions could serve as a perfect test of a proposal that the observed dark energy in present epoch is a result of a tiny deviation of the QCD vacuum energy in expanding universe from its conventional value in Minkowski spacetime.
Links to: arXiv, form interface, find, astro-ph, recent, 1201, contact, help (Access key information)