The new CRISP filter at the Swedish 1-m Solar Telescope provides
opportunities for observing the solar atmosphere with unprecedented spatial
resolution and cadence. In order to benefit from the high quality of
observational data from this instrument, we have developed methods for
calibrating and restoring polarized Stokes images, obtained at optical and near
infrared wavelengths, taking into account field-of-view variations of the
filter properties.
In order to facilitate velocity measurements, a time series from a 3D
hydrodynamical granulation simulation is used to compute quiet Sun spectral
line profiles at different heliocentric angles. The synthetic line profiles,
with their convective blueshifts, can be used as absolute references for
line-of-sight velocities.
Observations of the Ca II 8542 {\AA} line are used to study magnetic fields
in chromospheric fibrils. The line wings show the granulation pattern at
mid-photospheric heights whereas the overlying chromosphere is seen in the core
of the line. Using full Stokes data, we have attempted to observationally
verify the alignment of chromospheric fibrils with the magnetic field. Our
results suggest that in most cases fibrils are aligned along the magnetic field
direction, but we also find examples where this is not the case.
Detailed interpretation of Stokes data from spectral lines formed in the
chromospheric data can be made using non-LTE inversion codes. For the first
time, we use a realistic 3D MHD chromospheric simulation of the quiet Sun to
assess how well NLTE inversions recover physical quantities from
spectropolarimetric observations of Ca II 8542 {\AA}. We demonstrate that
inversions provide realistic estimates of depth-averaged quantities in the
chromosphere, although high spectral resolution and high sensitivity are needed
to measure quiet Sun chromospheric magnetic fields.
For more than four decades, scientists have been trying to find an answer to one of the most fundamental questions in paleoclimatology, the `faint young Sun problem'. For the early Earth, models of stellar evolution predict a solar energy input to the climate system which is about 25% lower than today. This would result in a completely frozen world over the first two billion years in the history of our planet, if all other parameters controlling Earth's climate had been the same. Yet there is ample evidence for the presence of liquid surface water and even life in the Archean (3.8 to 2.5 billion years before present), so some effect (or effects) must have been compensating for the faint young Sun. A wide range of possible solutions have been suggested and explored during the last four decades, with most studies focusing on higher concentrations of atmospheric greenhouse gases like carbon dioxide, methane or ammonia. All of these solutions present considerable difficulties, however, so the faint young Sun problem cannot be regarded as solved. Here I review research on the subject, including the latest suggestions for solutions of the faint young Sun problem and recent geochemical constraints on the composition of Earth's early atmosphere. Furthermore, I will outline the most promising directions for future research. In particular I would argue that both improved geochemical constraints on the state of the Archean climate system and numerical experiments with state-of-the-art climate models are required to finally assess what kept the oceans on the Archean Earth from freezing over completely.
We present moderate-resolution spectroscopic observations of red supergiants (RSGs) in the low-metallicity Local Group galaxies NGC 6822 (Z = 0.4Zsun) and Wolf-Lundmark-Melotte (WLM; Z = 0.1Zsun). By combining these observations with reduction techniques for multislit data reduction and flux calibration, we are able to analyze spectroscopic data of 16 RSGs in NGC 6822 and spectrophotometric data of 11 RSGs in WLM. Using these observations we determine spectral types for these massive stars, comparing them to Milky Way and Magellanic Clouds RSGs and thus extending observational evidence of the abundance-dependent shift of RSG spectral types to lower metallicities. In addition, we have uncovered two RSGs with unusually late spectral types (J000158.14-152332.2 in WLM, with a spectral type of M3 I, and J194453.46-144552.6 in NGC 6822, with a spectral type of M4.5 I) and a third RSG (J194449.96-144333.5 in NGC 6822) whose spectral type has varied from a M2.5 in 1997 to a K5 in 2008. All three of these stars could potentially be members of a recently-discovered class of extreme RSG variables.
There is now growing evidence that the soft X-ray excess is almost ubiquitous among unobscured active galaxies. In spite of the various interpretations that have been considered in the last years, the nature of this foremost spectral feature is not firmly established yet. In this context, we review from a reflection perspective the three highest-quality X-ray observations of the narrow-line Seyfert 1 galaxy Tonantzintla S180, obtained by XMM-Newton and Suzaku. The X-ray spectrum of Ton S180 shows only moderate variations over a time span of several years, suggesting that the same physical process accounts for the bulk of the broad-band X-ray emission at the different epochs, and that the properties of the X-ray source are fairly stable. We have successfully applied in our spectral analysis a dual-reflector model, consisting of two separate components: one arises from the surface of the accretion disc, is highly ionized and blurred by relativistic effects; the other is cold, quite faint, and can be associated with a distant reprocessor. Due to the strength and the nearly power-law shape of its soft excess emission, Ton S180 is one of the most challenging sources to test the X-ray reflection scenario. In this work we provide a clear illustration of the great potential and spectral flexibility of blurred reflection models, while discussing some of their current limitations and possible shortcomings.
We present results from a simple general circulation model (GCM) of a GJ 1214b-like super-Earth exoplanet. The dynamical core of our model is a scaled-up version of a shallow atmosphere, terrestrial planet GCM that has previously been used for Mars and therefore employs different boundary conditions and physical processes than downsized gas giant models. We assume the planet is tidally locked and has the observed characteristics of GJ-1214b [Charbonneau et al. 2009] for surface mass, surface radius, orbital period, and surface gravitational acceleration. We assume the atmosphere is composed entirely of water vapor. We assume the planet has a surface (i.e., a density discontinuity at depth), which will provide a mechanical drag and affect the radiative balance at the bottom boundary. We assume a gray atmosphere in the IR. We find that a westerly jet is present aloft at the equator and that the longitude of maximum temperature is shifted eastward of the substellar point. A wavenumber-1 feature is present in the equatorial vertical velocity field, indicative of a standing Kelvin and/or Rossby wave. As such, the circulation does not exhibit a cellular structure as on terrestrial Solar System planets; rather, air parcels move up and down on closed horizontal circuits and always return to the same location. The flow at the midlatitudes and poles has both an easterly component and a component that flows poleward along the evening terminator and equatorward along the morning terminator. Temperature inversions exist in the boundary layer and aloft, which are a result of variations in horizontal heat transport. The flow within the boundary layer is more convergent than in the atmosphere aloft. The surface pressure is higher at the poles than the equator. (abridged)
High-resolution, multi-wavelength, and time-domain observations of the
Galactic centre black hole candidate, Sgr A*, allow for a direct test of
contemporary accretion theory. To date, all models have assumed alignment
between the accretion disc and black hole angular momentum axes, but this is
unjustified for geometrically thick accretion flows like that onto Sgr A*.
Instead, we calculate images and spectra from a set of simulations of accretion
flows misaligned ('tilted') by 15 degrees from the black hole spin axis and
compare them with millimetre (mm) to near-infrared (NIR) observations.
Non-axisymmetric standing shocks from eccentric fluid orbits dominate the
emission, leading to a wide range of possible image morphologies. These effects
invalidate previous parameter estimates from model fitting, including estimates
of the dimensionless black hole spin, except possibly at low values of spin or
tilt. At 1.3mm, the images have crescent morphologies, and the black hole
shadow may still be accessible to future mm-VLBI observations. Shock heating
leads to high energy electrons (T > 10^12 K), which can naturally produce the
observed NIR flux, spectral index, and rapid variability ('flaring'). This NIR
emission is uncorrelated with that in the mm, which also agrees with
observations.
These are the first models to self-consistently explain the time-variable mm
to NIR emission of Sgr A*. Predictions of the model include significant
structural changes observable with mm-VLBI on both the dynamical (hour) and
Lense-Thirring precession (day-year) timescales; and ~30-50 microarcsecond
changes in centroid position from extreme gravitational lensing events during
NIR flares, detectable with the future VLT instrument GRAVITY. If the observed
NIR emission is caused by shock heating in a tilted accretion disc, then the
Galactic centre black hole has a positive, non-zero spin parameter (a > 0).
The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.
We have used the Green Bank Telescope to carry out a deep search for redshifted CO J=2-1 line emission from an extended (>17 kpc) Ly-Alpha blob (LAB), "Himiko", at z~6.595. Our non-detection of CO J=2-1 emission places the strong 3-sigma upper limit of L'_CO < 1.8x10^10 x sqrt(dV/250) K km/s pc^2 on the CO line luminosity. This is comparable to the best current limits on the CO line luminosity in LABs at z~3 and lower-luminosity Lyman-Alpha emitters (LAEs) at z>~6.5. High-z LABs appear to have lower CO line luminosities than the host galaxies of luminous quasars and sub-mm galaxies at similar redshifts, despite their high stellar mass. Although the CO-to-H2 conversion factor is uncertain for galaxies in the early Universe, we assume X_CO = 0.8 Msun (K km/s pc^2)^-1 to obtain the limit M(H_2) < 1.4 x 10^10 Msun on Himiko's molecular gas mass; this is a factor of >2.5 lower than the stellar mass in the z~6.595 LAB.
Intermediate mass black holes (IMBHs) with expected masses M_BH ~ 10^4 M_sun are thought to bridge the gap between stellar mass black holes (M_BH ~ 3 - 100 M_sun) and supermassive black holes found at the centre of galaxies (M_BH > 10^6 M_sun). Until today, no IMBH has been confirmed observationally. The most promising objects to host an IMBH as their central mass are globular clusters. Here, we present high sensitivity multi-epoch 1.6 GHz very long baseline interferometry observations of the globular cluster M15 that has been suggested to host an IMBH. Assuming the IMBH to be accreting matter from its surrounding we expect to detect it as a point source moving with the global motion of the cluster. However, we do not detect any such object within a radius of 6000 AU of the cluster centre in any of the five observations spread over more than one year. This rules out any variability of the putative IMBH on the time scale of one to two months. To get the most stringent upper limit for the flux density of the putative IMBH we concatenate the data of all five epochs. In this data we measure a 3{\sigma} upper flux limit of 10 {\mu}Jy for a central source. We employ the fundamental plane of black hole activity to estimate the mass of the central IMBH candidate. Based on previous X-ray observations of M15 our measurements indicate a 3{\sigma} upper mass limit of ~500 M_sun.
Galactic black hole coronae are composed of a hot, magnetized plasma. The spectral energy distribution produced in this component of X-ray binaries can be strongly affected by different interactions between locally injected relativistic particles and the matter, radiation and magnetic fields in the source. We study the non-thermal processes driven by the injection of relativistic particles into a strongly magnetized corona around an accreting black hole. We compute in a self-consistent way the effects of relativistic bremsstrahlung, inverse Compton scattering, synchrotron radiation, and the pair-production/annihilation of leptons, as well as hadronic interactions. Our goal is to determine the non-thermal broadband radiative output of the corona. The set of coupled kinetic equations for electrons, positrons, protons, and photons are solved and the resulting particle distributions are computed self-consistently. The spectral energy distributions of transient events in X-ray binaries are calculated, as well as the neutrino production. We show that the application to Cygnus X-1 of our model of non-thermal emission from a magnetized corona yields a good fit to the observational data. Finally, we show that the accumulated signal produced by neutrino bursts in black hole coronae might be detectable for sources within a few kpc on timescales of years. Our work leads to predictions for non-thermal photon and neutrino emission generated around accreting black holes, that can be tested by the new generation of very high energy gamma-ray and neutrino instruments.
We present an Ultra-Cool Dwarf (UCD) catalogue compiled from low southern Galactic latitudes and mid-plane, from a cross-correlation of the 2MASS and SuperCOSMOS surveys. The catalogue contains 246 members identified from 5042 sq. deg. within 220 deg. <= l <= 360 deg. and 0 deg. < l <= 30 deg., for |b| <= 15 deg. Sixteen candidates are spectroscopically confirmed in the near-IR as UCDs with spectral types from M7.5V to L9. Our catalogue selection method is presented enabling UCDs from ~M8V to the L-T transition to be selected down to a 2MASS limiting magnitude of Ks ~= 14.5 mag. This method does not require candidates to have optical detections for catalogue inclusion. An optimal set of optical/near-IR and reduced proper-motion selection criteria have been defined that includes: an Rf and Ivn photometric surface gravity test, a dual Rf-band variability check, and an additional photometric classification scheme to selectively limit contaminants. We identify four candidates as possible companions to nearby Hipparcos stars -- observations are needed to identify these as potential benchmark UCD companions. We also identify twelve UCDs within a possible distance 20 pc, three are previously unknown of which two are estimated within 10 pc, complimenting the nearby volume-limited census of UCDs. An analysis of the catalogue spatial completeness provides estimates for distance completeness over three UCD MJ ranges, while Monte-Carlo simulations provide an estimate of catalogue areal completeness at the 75 per cent level. We estimate a UCD space density of Rho (total) = (6.41+-3.01)x10^3/pc^3 over the range of 10.5 <= MJ ~< 14.9, similar to values measured at higher Galactic latitudes (|b| ~> 10 deg.) in the field population and obtained from more robust spectroscopically confirmed UCD samples.
Recent measurements of solar wind turbulence report the presence of intermittent, exponentially distributed angular discontinuities in the magnetic field. In this Letter, we study whether such discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the discontinuities by measuring the fluctuations of the magnetic field direction, Delta theta, across fixed spatial increments Delta x in direct numerical simulations of MHD turbulence with an imposed uniform guide field B_0. A large region of the probability density function (pdf) for Delta theta is found to follow an exponential decay, proportional to exp(-Delta theta/theta_*), with characteristic angle theta_* ~ (14 deg) (b_rms/B_0)^0.65 for a broad range of guide-field strengths. We find that discontinuities observed in the solar wind can be reproduced by MHD turbulence with reasonable ratios of b_rms/B_0. We also observe an excess of small angular discontinuities when Delta x becomes small, possibly indicating an increasing statistical significance of dissipation-scale structures. The structure of the pdf in this case closely resembles the two-population pdf seen in the solar wind. We thus propose that strong discontinuities are associated with inertial-range MHD turbulence, while weak discontinuities emerge from near-dissipation-range turbulence. In addition, we find that the structure functions of the magnetic field direction exhibit anomalous scaling exponents, which indicates the existence of intermittent structures.
We present results from monitoring observations of the gravitationally lensed quasar RX J1131-1231 performed with the Chandra X-ray Observatory. The X-ray observations were planned with relatively long exposures that allowed a search for energy-dependent microlensing in the soft (0.2-2 keV) and hard (2-10 keV) light curves of the images of RX J1131-1231. We detect significant microlensing in the X-ray light-curves of images A and D, and energy-dependent microlensing of image D. The magnification of the soft band appears to be larger than that in the hard band by a factor of ~ 1.3 when image D becomes more magnified. This can be explained by the difference between a compact, softer-spectrum corona that is producing a more extended, harder spectrum reflection component off the disk. This is supported by the evolution of the fluorescent iron line in image D over three consecutive time-averaged phases of the light curve. In the first period, a Fe line at E = 6.36(-0.16,+0.13) keV is detected (at > 99% confidence). In the second period, two Fe lines are detected, one at E = 5.47(-0.08,+0.06) keV (detected at > 99% confidence) and another at E = 6.02(-0.07,+0.09) keV (marginally detected at > 90% confidence), and in the third period, a broadened Fe line at 6.42(-0.15,+0.19) keV is detected (at > 99% confidence). This evolution of the Fe line profile during the microlensing event is consistent with the line distortion expected when a caustic passes over the inner disk where the shape of the fluorescent Fe line is distorted by General Relativistic and Doppler effects.
X-ray surveys contain sizable numbers of star forming galaxies, beyond the AGN which usually make the majority of detections. Many methods to separate the two populations are used in the literature, based on X-ray and multiwavelength properties. We aim at a detailed test of the classification schemes and to study the X-ray properties of the resulting samples. We build on a sample of galaxies selected at 1.4 GHz in the VLA-COSMOS survey, classified by Smolcic et al. (2008) according to their optical colours and observed with Chandra. A similarly selected control sample of AGN is also used for comparison. We review some X-ray based classification criteria and check how they affect the sample composition. The efficiency of the classification scheme devised by Smolcic et al. (2008) is such that ~30% of composite/misclassified objects are expected because of the higher X-ray brightness of AGN with respect to galaxies. The latter fraction is actually 50% in the X-ray detected sources, while it is expected to be much lower among X-ray undetected sources. Indeed, the analysis of the stacked spectrum of undetected sources shows, consistently, strongly different properties between the AGN and galaxy samples. X-ray based selection criteria are then used to refine both samples. The radio/X-ray luminosity correlation for star forming galaxies is found to hold with the same X-ray/radio ratio valid for nearby galaxies. Some evolution of the ratio may be possible for sources at high redshift or high luminosity, tough it is likely explained by a bias arising from the radio selection. Finally, we discuss the X-ray number counts of star forming galaxies from the VLA- and C-COSMOS surveys according to different selection criteria, and compare them to the similar determination from the Chandra Deep Fields. The classification scheme proposed here may find application in future works and surveys.
We present arcsecond resolution 1.4mm observations of the high mass star forming region, Sharpless 2-157, that reveal the cool dust associated with the first stages of star formation. These data are compared with archival images at optical, infrared, and radio wavelengths, and complemented with new arcsecond resolution mid-infrared data. We identify a dusty young HII region, numerous infrared sources within the cluster envelope, and four starless condensations. Three of the cores lie in a line to the south of the cluster peak, but the most massive one is right at the center and associated with a jumble of bright radio and infrared sources. This presents an interesting juxtaposition of high and low mass star formation within the same cluster which we compare with similar observations of other high mass star forming regions and discuss in the context of cluster formation theory.
Recent data from ATIC, CREAM and PAMELA indicate that the cosmic ray energy spectra of protons and nuclei exhibit a remarkable hardening at energies above 100 GeV per nucleon. We propose that the hardening is an interstellar propagation effect that originates from a spatial change of the cosmic ray transport properties in different regions of the Galaxy. The key hypothesis is that the diffusion coefficient is not separable into energy and space variables as usually assumed. Under this scenario, we can reproduce well the observational data. Our model has several implications for the cosmic ray acceleration/propagation physics and can be tested by ongoing experiments such as AMS or Fermi/LAT.
The Pierre Auger Observatory has associated a few ultra high energy cosmic
rays with the direction of Centaurus A. This source has been deeply studied in
radio, infrared, X-ray and $\gamma$-rays (MeV-TeV) because it is the nearest
radio-loud active galactic nuclei.
Its spectral energy distribution or spectrum shows two main peaks, the low
energy peak, at an energy of $10^{-2}$ eV, and the high energy peak, at about
150 keV.
There is also a faint very high energy (E $\geq$ 100 GeV) $\gamma$-ray
emission fully detected by the High Energy Stereoscopic System experiment. In
this work we describe the entire spectrum, the two main peaks with a
Synchrotron/Self-Synchrotron Compton model and, the Very High Energy emission
with a hadronic model. We consider p$\gamma$ and $pp$ interactions. For the
p$\gamma$ interaction, we assume that the target photons are those produced at
150 keV in the leptonic processes. On the other hand, for the pp interaction we
consider as targets the thermal particle densities in the lobes. Requiring a
satisfactory description of the spectra at very high energies with p$\gamma$
interaction we obtain an excessive luminosity in ultra high energy cosmic rays
(even exceeding the Eddington luminosity). However, when considering pp
interaction to describe the $\gamma$-spectrum, the obtained number of ultra
high energy cosmic rays are in agreement with Pierre Auger observations.
Moreover, we calculate the possible neutrino signal from pp interactions on a
Km$^3 $ neutrino telescope using Monte Carlo simulations.
We present observational evidence that leakage of ionising photons from star-forming regions can affect the quantification of the star formation rate (SFR) in galaxies. This effect could partially explain the differences between the SFR estimates using the far ultraviolet (FUV) and the Halpha emission. We find that leakage could decrease the SFR(Ha)/SFR(FUV) ratio by up to a 25 per cent. The evidence is based on the observation that the SFR(Ha)/SFR(FUV) ratio is lower for objects showing a shell Halpha structure than for regions exhibiting a much more compact morphology. The study has been performed on three object samples: low luminosity dwarf galaxies from the Local Volume Legacy survey and star-forming regions in the Large Magellanic Cloud and the nearby Local Group galaxy M33. For the three samples we find differences (1.1-1.4sigma) between the SFR(Ha)/SFR(FUV) for compact and shell objects. Although leakage cannot entirely explain the observed trend of SFR(Ha)/SFR(FUV) ratios for systems with low SFR, we show the mechanism can lead to different SFR estimates when using Halpha and FUV luminosities. Therefore, further study is needed to constrain the contribution of leakage to the low SFR(Ha)/SFR(FUV) ratios observed in dwarf galaxies and its impact on the Halpha flux as a SFR indicator in such objects.
We present long slit spectrophotometry of two HII regions: TOL2146-391 and TOL0357-3915. We performed a detailed analysis that involves abundance determinations relaxing the assumption of homogeneous temperature. The temperature inhomogeneities values, t^2, were obtained through two methods: (i) comparing abundances from oxygen recombination lines to abundances from collisionally excited lines and (ii) by using the line intensity ratios of a set of HeI lines together with the HELIO10 program. We find that the HELIO10 program is a good alternative to obtain a t^2 value in photoionized regions where recombination lines of heavy elements are not available. We have plotted 27 high and low metallicity HII regions in an oxygen degree of ionization versus t^2 diagram; we find areas populated by HII regions and areas void of them; the physical characteristics of each area are discussed. In addition, an average t^2 value can be determined for the objects in each area. We propose to use this <t^2> value for the cases where a direct measurement of t^2 cannot be determined
We present a multi-scale, multi-wavelength source extraction algorithm called getsources. Although it has been designed primarily for use in the far-infrared surveys of Galactic star-forming regions with Herschel, the method can be applied to many other astronomical images. Instead of the traditional approach of extracting sources in the observed images, the new method analyzes fine spatial decompositions of original images across a wide range of scales and across all wavebands. It cleans those single-scale images of noise and background, and constructs wavelength-independent single-scale detection images that preserve information in both spatial and wavelength dimensions. Sources are detected in the combined detection images by following the evolution of their segmentation masks across all spatial scales. Measurements of the source properties are done in the original background-subtracted images at each wavelength; the background is estimated by interpolation under the source footprints and overlapping sources are deblended in an iterative procedure. In addition to the main catalog of sources, various catalogs and images are produced that aid scientific exploitation of the extraction results. We illustrate the performance of getsources on Herschel images by extracting sources in sub-fields of the Aquila and Rosette star-forming regions. The source extraction code and validation images with a reference extraction catalog are freely available.
Observations with interferometric gravitational-wave detectors result in probability sky maps that are multimodal and spread over 10-100 deg^2. We present a scheme for maximizing the probability of imaging optical counterparts to gravitational-wave transients given limited observing resources. Our framework is capable of coordinating many telescopes with different fields of view and limiting magnitudes. We present a case study comparing three different planning algorithms. We find that, with the network of telescopes that was used in the most recent joint LIGO-Virgo science run, a relatively straightforward coordinated approach doubles the detection efficiency relative to each telescope observing independently.
Gamma-ray studies are an essential tool in our search for the origin of cosmic rays. Instruments like the Fermi-LAT, H.E.S.S., MAGIC and VERITAS have revolutionized our understanding of the high energy Universe. This paper describes the status of the very rich field of gamma-ray astrophysics that contains a wealth of data on Galactic and extragalactic particle accelerators. It is the write-up of a rapporteur talk given at the 32nd ICRC in Beijing, China in which new results were presented with an emphasis on the cosmic-ray related studies of the Universe.
We study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity of the local type. We use a large ensemble of simulated WL maps with survey specifications relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary certain parameters around fiducial values: the non-Gaussianity parameter f_NL, the matter density parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the spectral index of the primordial power spectrum n_s, and the dark-energy equation-of-state parameter w_0. We assess the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statistics, and the abundance of shear peaks to these parameters. We find that each of the considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on f_NL. Marginalized constraints from any individual WL probe are much weaker due to strong correlations between parameters. However, the parameter errors can be substantially reduced by combining information from different WL probes. Combining all WL probes yields the following marginal (68% confidence level) uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004, Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by neglecting f_NL on the constraints on the other parameters. We find sigma_8 and w_0 to be the most affected. Moreover, neglecting non-Gaussianity leads to a severe underestimation of the uncertainties in the other cosmological parameters.
A statistical study has been carried out of the relationship between plasma flow velocities and magnetic field parameters during the emergence of active regions at the solar photospheric level with SOHO/MDI data. We have investigated 224 emerging active regions with different spatial scales and position on the solar disk. The following relationships for the first hours of the emergence of active regions have been analysed: i) maximum negative Doppler velocities with the position of the emerging active regions on the solar disk; ii) maximum velocities of upflow and downflow of substance with the flux growth rate and magnetic field strength for the active regions emerging in the central part of the solar disk (the vertical component of plasma flows); iii) maximum positive and negative Doppler velocities with the flux growth rate for the active regions emerging near the limb (the horizontal component of plasma flows); iiii) the flux growth rate with the density of emerging magnetic fields. We have compared magnetic field parameters and the velocities for the first hours of the appearance of active regions with the total magnetic flux at the maximum of their development.
Tidal effects on clumps of material during random non-stationary accretion onto a black hole produce phenomena with distinct temporal characteristics in observed light-curves. During such non-stationary accretion events, the shape of the accreting object evolves in time, and observable quasi-periodic phenomena with variable quasi-periods are produced. A number of characteristic light-curves, obtained with numerical simulations, will be shown. Their relevance to observed phenomena will be briefly discussed.
We report the detection in emission of various highly excited rotational transitions of HNC (J = 6-5 through J =12-11) toward the carbon-star envelope IRC +10 216 using the HIFI instrument on-board the Herschel Space Observatory. Observations of the J = 1-0 and J = 3-2 lines of HNC with the IRAM 30-m telescope are also presented. The lines observed with HIFI have upper level energies corresponding to temperatures between 90 and 340 degrees Kelvin, and trace a warm and smaller circumstellar region than that seen in the interferometric maps of the J = 1-0 transition, whose emission extends up to a radius of 20". After a detailed chemical and radiative transfer modeling, we find that the presence of HNC in the circumstellar envelope of IRC +10 216 is consistent with formation from the precursor ion HCNH+, which in turn is produced through several proton transfer reactions which are triggered by the cosmic-ray ionization. We also find that the radiative pumping through 21 um photons to the first excited state of the bending mode v2 plays a crucial role to populate the high-J HNC levels involved in the transitions observed with HIFI. Emission in these high-J rotational transitions of HNC is expected to be strong in regions which are warm and dense and/or have an intense infrared flux at wavelengths around 21 um.
In this paper we study the dynamics of flavor transformation for neutrinos propagating in the very dense environment of astrophysical compact objects as Type II supernova in post collapse phase and proto-neutron stars. The analysis is based on the generalized Boltzmann equation, which incorporates the neutrinoneutrino interaction. We focus the analysis on the possible collective flavor dynamics, which can displays bipolar and synchronized flavor oscillations, and the associated transition from single-angle to multi-angle regime. To solve the kinetic equations we use an expansion in Legendre polynomials and a two-flavor scheme. On the basis of the numerical simulations we argue that neutrino flavor dynamics is suppressed by the presence of the electron component up to a certain distance from the neutrino sphere, where a transition to multi-angle regime occurs. This distance will move towards the neutrino sphere as the neutrino emission epoch proceeds.
Atmospheric modelling of the components of the visually close binary systems Hip70973 and Hip72479 was used to estimate the individual physical parameters of their components. The model atmospheres were constructed using a grid of Kurucz solar metalicity blanketed models, and used to compute a synthetic spectral energy distribution for each component separately, and hence for the combined system. The total observational spectral energy distributions of the systems were used as a reference for the comparison with the synthetic ones. We used the feedback modified parameters and iteration method to get the best fit between synthetic and observational spectral energy distributions. The physical parameters of the components of the system Hip70973 were derived as: Ta eff = 5700{\pm}75 K, Tb eff = 5400{\pm}75 K, log ga = 4.50{\pm}0.05, log gb = 4.50{\pm}0.05, Ra = 0.98{\pm}0.07R{\odot}, Rb = 0.89{\pm}0.07R{\odot}, and pi = 26.25 {\pm} 1.95 mas, with G4 & G9 spectral types. And those of the system Hip72479 as: Ta eff = 5400 {\pm} 50 K, Tb eff = 5180 {\pm} 50 K, log ga = 4.50 {\pm} 0.05, log gb = 4.60 {\pm} 0.05, Ra = 0.89 {\pm} 0.07R{\odot}, Rb = 0.80 {\pm} 0.07R{\odot}, and pi = 23.59 {\pm} 1.00 mas with G9 & K1 spectral types.
Understanding coronal mass ejection (CME) energetics and dynamics has been a long-standing problem, and although previous observational estimates have been made, such studies have been hindered by large uncertainties in CME mass. Here, the two vantage points of the Solar Terrestrial Relations Observatory (STEREO) COR1 and COR2 coronagraphs were used to accurately estimate the mass of the 2008 December 12 CME. Acceleration estimates derived from the position of the CME front in 3-D were combined with the mass estimates to calculate the magnitude of the kinetic energy and driving force at different stages of the CME evolution. The CME asymptotically approaches a mass of 3.4\pm1.0x10^15 g beyond ~10 R_sun. The kinetic energy shows an initial rise towards 6.3\pm3.7x10^29 erg at ~3 R_sun, beyond which it rises steadily to 4.2\pm2.5x10^30 erg at ~18 R_sun. The dynamics are described by an early phase of strong acceleration, dominated by a force of peak magnitude of 3.4\pm2.2x10^14N at ~3 R_sun, after which a force of 3.8\pm5.4x10^13 N takes affect between ~7-18 R_sun. These results are consistent with magnetic (Lorentz) forces acting at heliocentric distances of <7 R_sun, while solar wind drag forces dominate at larger distances (>7 R_sun).
We have performed accurate iron abundance measurements for 44 red giants
(RGs) in the Carina dwarf spheroidal (dSph) galaxy. We used archival,
high-resolution spectra (R~38,000) collected with UVES at ESO/VLT either in
slit mode (5) or in fiber mode (39, FLAMES/GIRAFFE-UVES). The sample is more
than a factor of four larger than any previous spectroscopic investigation of
stars in dSphs based on high-resolution (R>38,000$) spectra. We did not impose
the ionization equilibrium between neutral and singly-ionized iron lines. The
effective temperatures and the surface gravities were estimated by fitting
stellar isochrones in the V, B-V color-magnitude diagram. To measure the iron
abundance of individual lines we applied the LTE spectrum synthesis fitting
method using MARCS model atmospheres of appropriate metallicity. We found
evidence of NLTE effects between neutral and singly-ionized iron abundances.
Assuming that the FeII abundances are minimally affected by NLTE effects, we
corrected the FeI stellar abundances using a linear fit between FeI and FeII
stellar abundance determinations.
We found that the Carina metallicity distribution based on the corrected FeI
abundances (44 RGs) has a weighted mean metallicity of [Fe/H]=-1.80 and a
weighted standard deviation of sigma=0.24 dex. The Carina metallicity
distribution based on the FeII abundances (27 RGs) gives similar estimates
([Fe/H]=-1.72, sigma=0.24 dex). The current weighted mean metallicities are
slightly more metal poor when compared with similar estimates available in the
literature. Furthermore, if we restrict our analysis to stars with the most
accurate iron abundances, ~20 FeI and at least three FeII measurements (15
stars), we found that the range in iron abundances covered by Carina RGs (~1
dex) agrees quite well with similar estimates based on high-resolution spectra.
Little attention has so far been paid to the division of the observed population between pulsars of the two spin directions that are possible. Almost all pulsars with positive corotational charge density at the polar caps are expected to satisfy space-charge limited flow boundary conditions. Charge separation by blackbody photo-electric transitions in moving ions limits the acceleration potential, analogously with the more usually considered pair creation. But the limitation is more severe so that proton and ion energies can be relativistic but not ultra-relativistic, and these allow the growth of Langmuir-mode induced turbulence that couples directly with the radiation field, as shown by Asseo, Pelletier & Sol. The consequences of this, and of the several possible physical states of the polar cap, are described, qualitatively, as possible explanations for the complex phenomena of nulls, subpulse drift and mode-switching observed in subsets of pulsars.
Using the far-infrared emission, as observed by the Herschel Virgo Cluster Survey (HeViCS), and the integrated HI and CO brightness, we infer the dust and total gas mass for a magnitude limited sample of 35 metal rich spiral galaxies in Virgo. The CO flux correlates tightly and linearly with far-infrared fluxes observed by Herschel. Molecules in these galaxies are more closely related to cold dust rather than to dust heated by star formation or to optical/NIR brightness. We show that dust mass establishes a stronger correlation with the total gas mass than with the atomic or molecular component alone. The dust-to-gas ratio increases as the HI deficiency increases, but in highly HI deficient galaxies it stays constant. Dust is in fact less affected than atomic gas by weak cluster interactions, which remove most of the HI gas from outer and high latitudes regions. Highly disturbed galaxies, in a dense cluster environment, can instead loose a considerable fraction of gas and dust from the inner regions of the disk keeping constant the dust-to-gas ratio. There is evidence that the molecular phase is also quenched. This quencing becomes evident by considering the molecular gas mass per unit stellar mass. Its amplitude, if confirmed by future studies, highlights that molecules are missing in Virgo HI deficient spirals, but to a somewhat lesser extent than dust.
We describe the NASA/Stanford gyroscope relativity mission, Gravity Probe B (GP-B), and provide an overview of the following series of six astrometric and astrophysical papers that report on our radio observations and analyses made in support of this mission. The main goal of this 8.5 year program of differential VLBI astrometry was to determine the proper motion of the guide star of the GP-B mission, the RS CVn binary IM Pegasi (IM Peg; HR 8703). This proper motion is determined with respect to compact, extragalactic reference sources. The results are: -20.833 +- 0.090 mas/yr and -27.267 +- 0.095 mas/yr for, respectively, the right ascension and declination, in local Cartesian coordinates, of IM Peg's proper motion, and 10.370 +- 0.074 mas (i.e., 96.43 +- 0.69 pc) for its parallax (and distance). Each quoted uncertainty is meant to represent an ~70% confidence interval that includes the estimated contribution from systematic error. These results are accurate enough not to discernibly degrade the GP-B estimates of its gyroscopes' relativistic precessions: the frame-dragging and geodetic effects.
We analyse a large number ($> 500$) pointed RXTE observations of Cyg X-1 and model the spectrum of each one. A subset of the observations for which there is simultaneous reliable measure of the hardness ratio by the All Sky Monitor, shows that the sample covers nearly all the spectral shapes of Cyg X-1. The relative strength, width of the Iron line and the reflection parameter are in general correlated with the high energy photon spectral index $\Gamma$. This is broadly consistent with a geometry where for the hard state (low $\Gamma \sim 1.7$) there is a hot inner Comptonizing region surrounded by a truncated cold disk. The inner edge of the disk moves inwards as the source becomes softer till finally in the soft state (high $\Gamma > 2.2$) the disk fills the inner region and active regions above the disk produce the Comptonized component. However, the reflection parameter shows non-monotonic behaviour near the transition region ($\Gamma \sim 2$), suggestive of a more complex geometry or physical state of the reflector. Additionally, the inner disk temperature, during the hard state, is on the average higher than in the soft one, albeit with large scatter. These inconsistencies could be due to limitations in the data and the empirical model used to fit them. The flux of each spectral component is well correlated with $\Gamma$ which shows that unlike some other black hole systems, Cyg X-1 does not show any hysteresis behaviour. In the soft state, the flux of the Comptonized component is always similar to the disk one, which confirms that the ultra-soft state (seen in other brighter black hole systems) is not exhibited by Cyg X-1. The rapid variation of the Compton Amplification factor with $\Gamma$, naturally explains the absence of spectra with $\Gamma < 1.6$, despite a large number having $\Gamma \sim 1.65$.
We made VLBI observations at 8.4 GHz between 1997 and 2005 to estimate the coordinates of the "core" component of the superluminal quasar, 3C 454.3, the ultimate reference point in the distant universe for the NASA/Stanford Gyroscope Relativity Mission, Gravity Probe B. These coordinates are determined relative to those of the brightness peaks of two other compact extragalactic sources, B2250+194 and B2252+172, nearby on the sky, and within a celestial reference frame (CRF), defined by a large suite of compact extragalactic radio sources, and nearly identical to the International Celestial Reference Frame 2 (ICRF2). We find that B2250+194 and B2252+172 are stationary relative to each other, and also in the CRF, to within 1-sigma upper limits of 15 and 30 micro-arcsec/yr in RA and decl., respectively. The core of 3C 454.3 appears to jitter in its position along the jet direction over ~0.2 mas, likely due to activity close to the putative supermassive black hole nearby, but on average is stationary in the CRF within 1-sigma upper limits on its proper motion of 39 micro-arcsec/yr (1.0c) and 30 micro-arcsec/yr (0.8c) in RA and decl., respectively, for the period 2002 - 2005. Our corresponding limit over the longer interval, 1998 - 2005, of more importance to GP-B, is 46 and 56 micro-arcsec/yr in RA and decl., respectively. Some of 3C 454.3's jet components show significantly superluminal motion with speeds of up to ~200 micro-arcsec/yr or 5c in the CRF. The core of 3C 454.3 thus provides for Gravity Probe B a sufficiently stable reference in the distant universe.
We present the principal astrometric results of the very-long-baseline interferometry (VLBI) program undertaken in support of the Gravity Probe B (GP-B) relativity mission. VLBI observations of the GP-B guide star, the RS CVn binary IM Pegasi (HR 8703), yielded positions at 35 epochs between 1997 and 2005. We discuss the statistical assumptions behind these results and our methods for estimating the systematic errors. We find the proper motion of IM Peg in an extragalactic reference frame closely related to the International Celestial Reference Frame 2 (ICRF2) to be -20.83 +- 0.03 +- 0.09 mas/yr in right ascension and -27.27 +- 0.03 +- 0.09 mas/yr in declination. For each component the first uncertainty is the statistical standard error and the second is the total standard error (SE) including plausible systematic errors. We also obtain a parallax of 10.37 +- 0.07 mas (distance: 96.4 +- 0.7 pc), for which there is no evidence of any significant contribution of systematic error. Our parameter estimates for the ~25-day-period orbital motion of the stellar radio emission have SEs corresponding to ~0.10 mas on the sky in each coordinate. The total SE of our estimate of IM Peg's proper motion is ~30% smaller than the accuracy goal set by the GP-B project before launch: 0.14 mas/yr for each coordinate of IM Peg's proper motion. Our results ensure that the uncertainty in IM Peg's proper motion makes only a very small contribution to the uncertainty of the GP-B relativity tests.
When VLBI observations are used to determine the position or motion of a radio source relative to reference sources nearby on the sky, the astrometric information is usually obtained via: (i) phase-referenced maps; or (ii) parametric model fits to measured fringe phases or multiband delays. In this paper we describe a "merged" analysis technique which combines some of the most important advantages of these other two approaches. In particular, our merged technique combines the superior model-correction capabilities of parametric model fits with the ability of phase-referenced maps to yield astrometric measurements of sources that are too weak to be used in parametric model fits. We compare the results from this merged technique with the results from phase-referenced maps and from parametric model fits in the analysis of astrometric VLBI observations of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 nearby on the sky. In these studies we use central-core components of radio sources 3C 454.3 and B2250+194 as our positional references. We obtain astrometric results for IM Peg with our merged technique even when the source is too weak to be used in parametric model fits, and we find that our merged technique yields superior astrometric results to the phase-referenced mapping technique. We used our merged technique to estimate the proper motion and other astrometric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B mission.
We used 8.4 GHz VLBI images obtained at up to 35 epochs between 1997 and 2005
to examine the radio structures of the main reference source, 3C 454.3, and two
secondary reference sources, B2250+194 and B2252+172, for the guide star for
the NASA/Stanford relativity mission Gravity Probe B (GP-B). For one epoch in
2004 May, we also obtained images at 5.0 and 15.4 GHz. The 35 8.4 GHz images
for quasar 3C 454.3 confirm a complex, evolving, core-jet structure. We
identified at each epoch a component, C1, near the easternmost edge of the core
region. Simulations of the core region showed that C1 is located, on average,
0.18 +- 0.06 mas west of the unresolved "core" identified in 43 GHz images. We
also identified in 3C 454.3 at 8.4 GHz several additional components which
moved away from C1 with proper motions ranging in magnitude between 0.9c and
5c. The detailed motions of the components exhibit two distinct bends in the
jet axis located ~3 and ~5.5 mas west of C1. The spectra between 5.0 and 15.4
GHz for the "moving" components are steeper than that for C1. The 8.4 GHz
images of B2250+194 and B2252+172, in contrast to those of 3C 454.3, reveal
compact structures. The spectrum between 5.0 and 15.4 GHz for B2250+194 is
inverted while that for B2252+172 is flat.
Based on its position near the easternmost edge of the 8.4 GHz radio
structure, close spatial association with the 43 GHz core, and relatively flat
spectrum, we believe 3C 454.3 component C1 to be the best choice for the
ultimate reference point for the GP-B guide star. The compact structures and
inverted to flat spectra of B2250+194 and B2252+172 make these objects valuable
secondary reference sources
We present a physical interpretation for the locations of the sources of radio emission in IM Pegasi (IM Peg, HR 8703), the guide star for the NASA/Stanford relativity mission Gravity Probe B. This emission is seen in each of our 35 epochs of 8.4-GHz VLBI observations taken from 1997 to 2005. We found that the mean position of the radio emission is at or near the projected center of the primary to within about 27% of its radius, identifying this active star as the radio emitter. The positions of the radio brightness peaks are scattered across the disk of the primary and slightly beyond, preferentially along an axis with position angle, p.a. = (-38 +- 8) deg, which is closely aligned with the sky projections of the orbit normal (p.a. = -49.5 +- 8.6 deg) and the expected spin axis of the primary. Comparison with simulations suggests that brightness peaks are 3.6 (+0.4,-0.7) times more likely to occur (per unit surface area) near the pole regions of the primary (|latitude| >= 70 deg) than near the equator (|latitude| <= 20 deg), and to also occur close to the surface with ~2/3 of them at altitudes not higher than 25% of the radius of the primary.
We present measurements of the total radio flux density as well as very-long-baseline interferometry (VLBI) images of the star, IM Pegasi, which was used as the guide star for the NASA/Stanford relativity mission Gravity Probe B. We obtained flux densities and images from 35 sessions of observations at 8.4 GHz (wavelength = 3.6 cm) between 1997 January and 2005 July. The observations were accurately phase-referenced to several extragalactic reference sources, and we present the images in a star-centered frame, aligned by the position of the star as derived from our fits to its orbital motion, parallax, and proper motion. Both the flux density and the morphology of IM Peg are variable. For most sessions, the emission region has a single-peaked structure, but 25% of the time, we observed a two-peaked (and on one occasion perhaps a three-peaked) structure. On average, the emission region is elongated by 1.4 +- 0.4 mas (FWHM), with the average direction of elongation being close to that of the sky projection of the orbit normal. The average length of the emission region is approximately equal to the diameter of the primary star. No significant correlation with the orbital phase is found for either the flux density or the direction of elongation, and no preference for any particular longitude on the star is shown by the emission region.
We investigate the possibility that the observed rotation of galaxies can be accounted for by invoking a massive baryonic disc with no need for non-baryonic dark matter or a massive halo. There are 5 primary reasons for suggesting this: 1. there are well known disc surface mass density distributions that naturally produce the observed rotation curves of galaxies. 2. there are a number of rotation curve `puzzles' that cannot be explained by a massive dark matter halo i.e. the success of maximum disc fitting, HI gas scaling to the observed rotation, the disc/halo conspiracy and the interpretation of the Tully-Fisher relation. 3. recent 21cm observations show an almost constant HI surface density and a distinct `cut-off' or edge to galactic discs. We explain this constant surface density in terms of either an optical depth effect or the onset of molecular gas formation and hence the possibility of considerably more gas existing in galaxies. We suggest that the HI cut-off does indeed mark the edge of the galactic disc. 4. there have been an increasing number of recent observations that imply that X_CO may be ten times higher in the outer Galaxy. This `dark' gas may provide adequate mass to account for galaxy rotation. 5. we show that the additional baryonic mass required to account for the rotation of galaxies is just that required to reconcile observed baryons with those predicted by big bang nucleosynthesis. Mestel discs can be used to straight forwardly explain the scaling laws of galaxies, particularly the observed relation between rotation velocity and radius and the oft used Tully-Fisher relation. We discuss observations of the baryonic content of galactic discs and where sufficient `hidden' baryons might be found to account for the rotation.
Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon feVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of feVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 105 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines.
Several models have been proposed to explain the dark energy that is causing universe expansion to accelerate. Here the acceleration predicted by the Holographic Dark Information Energy (HDIE) model is compared to the acceleration that would be produced by a cosmological constant. While identical to a cosmological constant at low redshifts, z<1, the HDIE model results in smaller Hubble parameter values at higher redshifts, z>1, reaching a maximum difference of 2.6\pm0.5% around z \sim 1.7. The next generation of dark energy measurements, both those scheduled to be made in space (ESA's Euclid and NASA's WFIRST missions) and those to be made on the ground (BigBOSS, LSST and Dark Energy Survey), should be capable of determining whether such a difference signature exists or not. The HDIE model is therefore falsifiable.
A number of physical processes show some form of bifurcation or splintering around a given point. Recently, it has been noted that cavity searches for interactions between photons and exotic fields may also show a bifurcation[1]. This paper describes the simulation of bifurcation of an optical beam in the presence of periodic focusing. The work is applicable to searches for exotic particles that employ cavities to extend the optical path of a laser beam through an interacting field. Significantly, simulations of cavity bifurcation reveal both a redistribution of the beam's energy density and a shifting of the beam's profile.
There is a growing trend toward using high-level tools for design and
implementation of radio astronomy digital signal processing (DSP) systems. Such
tools, for example, those from the Collaboration for Astronomy Signal
Processing and Electronics Research (CASPER), are usually platform-specific,
and lack high-level, platform-independent, portable, scalable application
specifications. This limits the designer's ability to experiment with designs
at a high-level of abstraction and early in the development cycle.
We address some of these issues using a model-based design approach employing
dataflow models. We demonstrate this approach by applying it to the design of a
tunable digital downconverter (TDD) used for narrow-bandwidth spectroscopy. Our
design is targeted toward an FPGA platform, called the Interconnect Break-out
Board (IBOB), that is available from the CASPER. We use the term TDD to refer
to a digital downconverter for which the decmation factor and center frequency
can be reconfigured without the need for regenerating the hardware code. Such a
design is currently not available in the CASPER DSP library.
The work presented in this paper focuses on two aspects. Firstly, we
introduce and demonstrate a dataflow-based design approach using the dataflow
interchange format (DIF) tool for high-level application specification, and we
integrate this approach with the CASPER tool flow. Secondly, we explore the
trade-off between the flexibility of TDD designs and the low hardware cost of
fixed-configuration digital downconverter (FDD) designs that use the available
CASPER DSP library. We further explore this trade-off in the context of a
two-stage downconversion scheme employing a combination of TDD or FDD designs.
We study the regular or chaotic nature of motion in a disk galaxy with a dense nucleus and an asymmetric dark halo. Two cases, the 2D model and the 3D model, are investigated. In the 2D model, a considerable fraction of the phase plane is covered by chaotic orbits. Two factors seem to be responsible for the chaotic motion: (i) the dense nucleus and (ii) the asymmetries in the dark halo. Our numerical experiments suggest, that there are several chaotic components on the Poincare phase plane. Different chaotic components are induced by the asymmetries in the halo. Each chaotic component seems to have a different value of the Lyapunov Characteristic Exponent, for small values of the asymmetry parameter and a unique LCE for larger values of the asymmetry parameter. A comparison of the present results with outcomes from previous work is also presented.
We present an analysis of gamma-ray data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope in the region around SNR S147 (G180.0-1.7). A spatially extended gamma-ray source detected in an energy range of 0.2--10 GeV is found to coincide with SNR S147. We confirm its spatial extension at >5sigma confidence level. The gamma-ray flux is (3.8 \pm 0.6) x 10^{-8} photons cm^{-2} s^{-1}, corresponding to a luminosity of 1.3 x 10^{34} (d/1.3 kpc)^2 erg s^{-1} in this energy range. The gamma-ray emission exhibits a possible spatial correlation with prominent Halpha filaments of S147. There is no indication that the gamma-ray emission comes from the associated pulsar PSR J0538+2817. The gamma-ray spectrum integrated over the remnant is likely dominated by the decay of neutral pi mesons produced through the proton--proton collisions in the filaments. Reacceleration of pre-existing CRs and subsequent adiabatic compression in the filaments is sufficient to provide the required energy density of high-energy protons.
We report measurements of the carbon monoxide ground state rotational transition (12C16O J = 1--0) with the Zpectrometer ultra-wideband spectrometer on the 100-m diameter Green Bank Telescope. The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 24 targets identified by Herschel-ATLAS photometric colors from the SPIRE instrument. Nine of the CO measurements are new redshift determinations, substantially adding to the number of detections of galaxies with rest-frame peak submillimeter emission near 100um. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350um-selected galaxies is strikingly similar to the optical redshifts of 850um-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by foreground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding approximate gravitational lens magnifications mu from spectroscopic data alone, yielding mu ~ 3--20. Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 10^13 L_sun.
Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.
We investigate the stability of prograde versus retrograde planets in circular binary systems using numerical simulations. We show that retrograde planets are stable up to distances closer to the perturber than prograde planets. We develop an analytical model to compute the prograde and retrograde mean motion resonances' locations and separatrices. We show that instability is due to single resonance forcing, or caused by nearby resonances' overlap. We validate our results regarding the role of single resonances and resonances' overlap on orbit stability, by computing surfaces of section of the CR3BP. We conclude that the observed enhanced stability of retrograde planets with respect to prograde planets is due to essential differences between the phase-space topology of retrograde versus prograde resonances (at p/q mean motion ratio, prograde resonance is of order p - q while retrograde resonance is of order p + q).
In the light of the recent discovery of a neutron star with a mass accurately determined to be almost two solar masses, it has been suggested that hyperons cannot play a role in the equation of state of dense matter in $\beta$-equilibrium. We re-examine this issue in the most recent development of the quark-meson coupling model. Within a relativistic Hartree-Fock approach and including the full tensor structure at the vector-meson-baryon vertices, we find that not only must hyperons appear in matter at the densities relevant to such a massive star but that the maximum mass predicted is completely consistent with the observation.
It was recently shown that a black hole (BH) is the only compact object that can acquire a scalar charge in scalar Gauss-Bonnet (sGB) theory under the small coupling approximation. This leads to the fact that scalar radiation is emitted from a binary containing at least one BH. In this letter, we find the constraints on this theory from BH low-mass X-ray binaries (BH-LMXBs). First, we calculate the constraint on this theory from the orbital decay rate of A0620-00 and we find it to be more than six orders of magnitude stronger than the solar system bound. Next, we look at XTE J1118+480, whose orbital decay rate has been recently measured with an excess compared to the theoretical prediction in GR due to the radiation reaction. The cause of this excess is currently unknown. Although it is likely that the cause is of astrophysical origin, here we investigate the possibility of explaining this excess with the additional scalar radiation in sGB theory. We find that there still remains a parameter range where the excess can be explained while also satisfying the constraint obtained from A0620-00. The interesting point is that for most of other alternative theories of gravity, it seems difficult to explain this excess with the additional radiation. This is because it would be difficult to evade the constraints from binary pulsars or they have already been constrained rather strongly from other observations such as solar system experiments. We propose several ways to determine whether the excess is caused by the scalar radiation in sGB gravity including future gravitational wave observations with space-borne interferometers, which can give a constraint three orders of magnitude stronger than that from A0620-00.
We present the computation of the loop-induced self-annihilation of dark matter particles into two photons in the framework of the NMSSM. This process is a theoretically clean observable with a "smoking-gun" signature but is experimentally very challenging to detect. The rates were computed with the help of the SloopS program, an automatic code initially designed for the evaluation of processes at the one-loop level in the MSSM. We focused on a light neutralino scenario and discuss how the signal can be enhanced in the NMSSM with respect to the MSSM and then compared with the present limits given by the dedicated search of the FERMI-LAT satellite on the monochromatic gamma lines.
Spherical symmetry for f(R)-gravity is discussed by searching for Noether symmetries. The method consists in selecting conserved quantities in form of currents that reduce dynamics of f(R)-models compatible with symmetries. In this way we get a general method to obtain constants of motion without setting a priori the form of f(R). In this sense, the Noether symmetry results a physical criterium. Relevant cases are discussed.
Links to: arXiv, form interface, find, astro-ph, recent, 1204, contact, help (Access key information)
Observations towards IRC +10216 of CS, SiO, SiS, NaCl, KCl, AlCl, AlF, and NaCN have been carried out with the IRAM 30-m telescope in the 80-357.5 GHz frequency range. A large number of rotational transitions covering a wide range of energy levels, including highly excited vibrational states, are detected in emission and serve to trace different regions of the envelope. Radiative transfer calculations based on the LVG formalism have been performed to derive molecular abundances from the innermost out to the outer layers. The excitation calculations include infrared pumping to excited vibrational states and inelastic collisions, for which up-to-date rate coefficients for rotational and, in some cases, ro-vibrational transitions are used. We find that in the inner layers CS, SiO, and SiS have abundances relative to H$_2$ of 4e-6, 1.8e-7, and 3e-6, respectively, and that CS and SiS have significant lower abundances in the outer envelope, which implies that they actively contribute to the formation of dust. Moreover, in the inner layers, the amount of sulfur and silicon in gas phase molecules is only 27 % for S and 5.6 % for Si, implying that these elements have already condensed onto grains, most likely in the form of MgS and SiC. Metal-bearing molecules lock up a relatively small fraction of metals, although our results indicate that NaCl, KCl, AlCl, AlF, and NaCN, despite their refractory character, are not significantly depleted in the cold outer layers. In these regions a few percent of the metals Na, K, and Al survive in the gas phase, either in atomic or molecular form, and are therefore available to participate in the gas phase chemistry in the outer envelope.
We present the results of a multi-wavelength (FUV to I-band) survey of the stellar populations of the globular cluster NGC 6752, using STIS, ACS and WFC3 on board the Hubble Space Telescope. We have confirmed that two previously identified CV candidates are, in fact, dwarf novae which underwent outbursts during our observations. We have also identified previously unknown optical counterparts to two X-ray sources. We estimate the position of the centre of the cluster, and show that the stellar density profile is not well described by a single King model, indicating that this cluster is in a core-collapsed or post-core collapse phase. The colour-magnitude diagram shows a well-populated horizontal branch, numerous blue stragglers and white dwarfs (WDs), as well as 87 sources in the gap region where we expect to find WD - main sequence binaries, including cataclysmic variables (CVs). The X-ray sources and WD binary systems are the most centrally concentrated populations, with dynamically estimated characteristic masses >1.1Msun and >0.8Msun, respectively.
In this note we revisit the Fisher information content of cosmological power spectra of Gaussian fields, when based on the assumption of a multivariate Gaussian likelihood for estimators. We discuss that while the assumption of a Gaussian likelihood is motivated by the central limit theorem, it leads if used consistently to a Fisher information content that violates the Cramer-Rao inequality, due to the presence of independent information from the parameter dependent covariance matrix. At any fixed multipole, this term is shown to become dominant in the limit of a large number of correlated fields. While the distribution of the estimators does indeed tend to a Gaussian with a large number of modes, it is shown, however, that its Fisher information content does not, in the sense that the covariance matrix never carries independent information content. The reason why the information content of the spectra is correctly described by the usual formula (i.e. without the covariance term) in this estimator perspective is precisely the fact the the estimators have a chi-squared like distribution, and not a Gaussian distribution. The assumption of a Gaussian estimators likelihood is thus from the point of view of the information neither necessary nor really adequate, and we warn against the use of Gaussian likelihoods with parameter dependent covariance matrices for parameter inference from such spectra.
We present a detailed analysis of redshift-space distortions in the two-point correlation function of the 6dF Galaxy Survey (6dFGS). The K-band selected sub-sample which we employ in this study contains 81971 galaxies distributed over 17000deg^2 with an effective redshift z = 0.067. By modelling the 2D galaxy correlation function, xi(r_p,pi), we measure the parameter combination f(z)sigma_8(z) = 0.423 +/- 0.055. Alternatively, by assuming standard gravity we can break the degeneracy between sigma_8 and the galaxy bias parameter, b. Combining our data with the Hubble constant prior from Riess et al (2011), we measure sigma_8 = 0.76 +/- 0.11 and Omega_m = 0.250 +/- 0.022, consistent with constraints from other galaxy surveys and the Cosmic Microwave Background data from WMAP7. Combining our measurement of fsigma_8 with WMAP7 allows us to test the relationship between matter and gravity on cosmic scales by constraining the growth index of density fluctuations, gamma. Using only 6dFGS and WMAP7 data we find gamma = 0.547 +/- 0.088, consistent with the prediction of General Relativity. We note that because of the low effective redshift of 6dFGS our measurement of the growth rate is independent of the fiducial cosmological model (Alcock-Paczynski effect). We also show that our conclusions are not sensitive to the model adopted for non-linear redshift-space distortions. Using a Fisher matrix analysis we report predictions for constraints on fsigma_8 for the WALLABY survey and the proposed TAIPAN survey. The WALLABY survey will be able to measure fsigma_8 with a precision of 4-10%, depending on the modelling of non-linear structure formation. This is comparable to the predicted precision for the best redshift bins of the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that low-redshift surveys have a significant role to play in future tests of dark energy and modified gravity.
The Herschel Reference Survey (HRS) is a guaranteed time Herschel key project aimed at studying the physical properties of the interstellar medium in galaxies of the nearby universe. This volume limited, K-band selected sample is composed of galaxies spanning the whole range of morphological types (from ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster). We present flux density measurements of the whole sample of 323 galaxies of the HRS in the three bands of the Spectral and Photometric Imaging Receiver (SPIRE), at 250, 350 and 500 microns. Aperture photometry is performed on extended galaxies and point spread function (PSF) fitting on timeline data for unresolved objects; we carefully estimate errors and upper limits. The flux densities are found to be in good agreement with those of the HeViCS and KINGFISH key projects in all SPIRE bands, and of the Planck consortium at 350 and 550 microns, for the galaxies in common. This submillimetre catalogue of nearby galaxies is a benchmark for the study of the dust properties in the local universe, giving the zero redshift reference for any cosmological survey.
We have obtained high-resolution data of the z 2 ring-like, clumpy star-forming galaxy (SFG) ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing-limited mode (in H- and J-band). Our data includes all of the main strong optical emission lines: [OII], [OIII], Ha, Hb, [NII] and [SII]. We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's massive, star-forming clumps (sigma \sim 85 km s^-1) and even broader wings (up to 70% of the total Ha flux, with sigma \sim 290 km s^-1) in regions spatially offset from the clumps by \sim 2 kpc. The broad emission likely originates from large-scale outflows with mass outflow rates from individual clumps that are 1-8x the SFR of the clumps. Based on emission line ratio diagnostics ([NII]/Ha and [SII]/Ha) and photoionization and shock models, we find that the emission from the clumps is due to a combination of photoionization from the star-forming regions and shocks generated in the outflowing component, with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6x10^7-10^8 cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm^-3 in and around the clumps, and ionized gas column densities of 1200-8000 Msol/pc^2), and SFR (10-40 Msol/yr), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps are quite different in terms of their internal properties, energetics and relative ages, and thus we are given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant clumps at high-z.
We investigate the pulsational stability of massive (M >~ 120 Msun) main sequence stars of a range of metallicities, including primordial, Population III stars. We include a formulation of convective damping motivated by numerical simulations of the interaction between convection and periodic shear flows. We find that convective viscosity is likely strong enough to stabilize radial pulsations whenever nuclear-burning (the epsilon-mechanism) is the dominant source of driving. This suggests that massive main sequence stars with Z <~ 2 x 10^-3 are pulsationally stable and are unlikely to experience pulsation-driven mass loss on the main sequence. These conclusions are, however, sensitive to the form of the convective viscosity and highlight the need for further high-resolution simulations of the convection-oscillation interaction. For more metal-rich stars (Z >~ 2 x 10^-3), the dominant pulsational driving arises due to the kappa-mechanism arising from the iron-bump in opacity and is strong enough to overcome convective damping. Our results highlight that even for oscillations with periods a few orders of magnitude shorter than the outer convective turnover time, the "frozen-in" approximation for the convection-oscillation interaction is inappropriate, and convective damping should be taken into account when assessing mode stability.
The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3\sigma$ limit). This upper limit is driven by the non-detection at 160 $\mu$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.
A critical challenge in measuring the power spectrum of 21cm emission from cosmic reionization is compensating for the frequency dependence of an interferometer's sampling pattern, which can cause smooth-spectrum foregrounds to appear unsmooth and degrade the separation between foregrounds and the target signal. In this paper, we present an approach to foreground removal that explicitly accounts for this frequency dependence. We apply the delay transformation introduced in Parsons & Backer (2009) to each baseline of an interferometer to concentrate smooth-spectrum foregrounds within the bounds of the maximum geometric delays physically realizable on that baseline. By focusing on delay-modes that correspond to image-domain regions beyond the horizon, we show that it is possible to avoid the bulk of smooth-spectrum foregrounds. We show that delay-modes that are uncorrupted by foregrounds also represent samples of the three-dimensional power spectrum, and can be used to constrain cosmic reionization. Because it uses only spectral smoothness to differentiate foregrounds from the targeted 21cm signature, this per-baseline analysis approach relies on spectrally- and spatially-smooth instrumental responses for foreground removal. For sufficient levels of instrumental smoothness relative to the brightness of interfering foregrounds, this technique substantially reduces the level of calibration previously thought necessary to detect 21cm reionization. As a result, this approach places fewer constraints on antenna configuration within an array, and in particular, facilitates the adoption of configurations that are optimized for power-spectrum sensitivity. Under these assumptions, we demonstrate the potential for the PAPER array to detect 21cm reionization at an amplitude of 10 mK^2 near k~0.2h Mpc^-1 with 128 dipoles in 7 months of observing.
High-velocity halo stars provide important information about the properties of the extreme Galactic halo. The study of unbound and bound Population II stars permits us to better estimate the mass of the halo.} {We carried out a detailed spectroscopic and kinematic study and have significantly refined the distance and the evolutionary state of the star. Its atmospheric parameters, chemical abundances and kinematical properties were determined using high-resolution optical spectroscopy and employing the local-thermodynamic-equilibrium model atmospheres of Kurucz and the spectral analysis code {\sc moog}. We found that CD-62$\degr$1346 is a metal-poor ([Fe/H]=$-$1.7) evolved giant star with $T_{\rm eff}=5300$ K and $\log g=1.7$. The star exhibits high carbon and s-element abundances typical of CH stars. It is also a lead star. Our kinematic analysis of its 3D space motions shows that this star has a highly eccentric ($e=0.91$) retrograde orbit with an apogalactic distance of $\sim 100$ kpc, exceeding by a factor of three the distance to the Magellanic Clouds. The star travels with very high velocity relative to the Galactocentric reference frame ($V_{\rm GRF}=570$ km\,s$^{-1}$). CD-62$\degr$1346 is an evolved giant star and not a subgiant star, as was considered earlier. Whether it is bound or unbound to the Galaxy depends on the assumed mass and on the adopted Galactic potential. We also show that the star HD 5223 is another example of a high-velocity CH star that exceeds the Galactic escape velocity. Possible origins of these two high-velocity stars are briefly discussed. CD-62$\degr$1346 and HD 5223 are the first red giant stars to join the restricted group of hypervelocity stars.
We present theoretical calculations of the quasi-molecular line profiles for the Lyman (Lya, Lyb, Lyg, Lyd) and Balmer (Ha, Hb, Hg, Hd, He, H8, H9, and H10) series, perturbed by collisions with protons. In all calculations we include the dependence of the dipole moments as a function of internuclear distance during the collision. The broadening from ion collisions must be added to the normal electron Stark broadening.
We propose Paraiso, a domain specific language embedded in functional
programming language Haskell, for automated tuning of explicit solvers of
partial differential equations (PDEs) on GPUs as well as multicore CPUs. In
Paraiso, one can describe PDE solving algorithms succinctly using tensor
equations notation. Hydrodynamic properties, interpolation methods and other
building blocks are described in abstract, modular, re-usable and combinable
forms, which lets us generate versatile solvers from little set of Paraiso
source codes.
We demonstrate Paraiso by implementing a compressive hydrodynamics solver. A
single source code less than 500 lines can be used to generate solvers of
arbitrary dimensions, for both multicore CPUs and GPUs. We demonstrate both
manual annotation based tuning and automated tuning of the program.
We report measurement of trigonometric parallax of IRAS 05168+3634 with VERA. The parallax is 0.532 +/- 0.053 mas, corresponding to a distance of 1.88+0.21/-0.17 kpc. This result is significantly smaller than the previous distance estimate of 6 kpc based on kinematic distance. This drastic change in the source distance revises not only physical parameters of IRAS 05168+3634, but also its location of the source, placing it in the Perseus arm rather than the Outer arm. We also measure proper motions of the source. A combination of the distance and the proper motions with systemic velocity yields rotation velocity ({\Theta}) of 227+9/-11 km s-1 at the source, assuming {\Theta}0 = 240 km s-1. Our result combined with previous VLBI results for six sources in the Perseus arm indicates that the sources rotate systematically slower than the Galactic rotation velocity at the LSR. In fact, we show observed disk peculiar motions averaged over the seven sources in the Perseus arm as (Umean, Vmean) = (11 +/- 3, -17 +/- 3) km s-1, indicating that these seven sources are systematically moving toward the Galactic center, and lag behind the Galactic rotation.
High-magnification microlensing events provide an important channel to detect planets. Perturbations near the peak of a high-magnification event can be produced either by a planet or a binary companion. It is known that central perturbations induced by both types of companions can be generally distinguished due to the basically different magnification pattern around caustics. In this paper, we present a case of central perturbations for which it is difficult to distinguish the planetary and binary interpretations. The peak of a lensing light curve affected by this perturbation appears to be blunt and flat. For a planetary case, this perturbation occurs when the source trajectory passes the negative perturbation region behind the back end of an arrowhead-shaped central caustic. For a binary case, a similar perturbation occurs for a source trajectory passing through the negative perturbation region between two cusps of an astroid-shaped caustic. We demonstrate the degeneracy for 2 high-magnification events of OGLE-2011-BLG-0526 and OGLE-2011-BLG-0950/MOA-2011-BLG-336. For OGLE-2011-BLG-0526, the $\chi^2$ difference between the planetary and binary model is $\sim$ 3, implying that the degeneracy is very severe. For OGLE-2011-BLG-0950/MOA-2011-BLG-336, the stellar binary model is formally excluded with $\Delta \chi^2 \sim$ 105 and the planetary model is preferred. However, it is difficult to claim a planet discovery because systematic residuals of data from the planetary model are larger than the difference between the planetary and binary models. Considering that 2 events observed during a single season suffer from such a degeneracy, it is expected that central perturbations experiencing this type of degeneracy is common.
In this paper, we revisit generalized Chaplygin gas (GCG) model as a unified dark matter and dark energy model. The energy density of GCG model is given as $\rho_{GCG}/\rho_{GCG0}=[B_{s}+(1-B_{s})a^{-3(1+\alpha)}]^{1/(1+\alpha)}$, where $\alpha$ and $B_s$ are two model parameters which will be constrained by type Ia supernova as standard candles, baryon acoustic oscillation as standard rulers and the seventh year full WMAP data points. In this paper, we will not separate GCG into dark matter and dark energy parts any more as adopted in the literatures. By using Markov Chain Monte Carlo method, we find the result: $\alpha=0.00126_{- 0.00126- 0.00126}^{+ 0.000970+ 0.00268}$ and $B_s= 0.775_{- 0.0161- 0.0338}^{+ 0.0161+ 0.0307}$.
LWA1 is a new radio telescope operating in the frequency range 10-88 MHz, located in central New Mexico. The telescope consists of 258 pairs of dipole-type antennas whose outputs are individually digitized and formed into beams. Simultaneously, signals from all dipoles can be recorded using one of the instrument's "all dipoles" modes, facilitating all-sky imaging. Notable features of the instrument include high intrinsic sensitivity (about 6 kJy zenith system equivalent flux density), large instantaneous bandwidth (up to 78 MHz), and 4 independently-steerable beams utilizing digital "true time delay" beamforming. This paper summarizes the design of LWA1 and its performance as determined in commissioning experiments. We describe the method currently in use for array calibration, and report on measurements of sensitivity and beamwidth.
In this paper we report on Expanded Very Large Array radio and Chandra and Swift X-ray observations of the outburst decay of the transient black hole candidate MAXI J1659-152 in 2011. We discuss the distance to the source taking the high inclination into account and we conclude that the source distance is probably 6+-2 kpc. The lowest observed flux corresponds to a luminosity of 2x10^31 (d/6 kpc)^2 erg/s This, together with the orbital period of 2.4 hr reported in the literature, suggests that the quiescent X-ray luminosity is higher than predicted on the basis of the orbital period -- quiescent X-ray luminosity relationship. The relation between the accretion and ejection mechanisms can be studied using the observed correlation between the radio and X-ray luminosities as these evolve over an outburst. We determine the behaviour of MAXI J1659-152 in the radio -- X-ray diagram at low X-ray luminosities using the observations reported in this paper and at high X-ray luminosities using values reported in the literature. At high X-ray luminosities the source lies closer to the sources that follow a correlation index steeper than 0.6-0.7. However, when compared to other sources that follow a steeper correlation index, the X-ray luminosity in MAXI J1659-152 is also lower. The latter can potentially be explained by the high inclination of MAXI J1659-152 if the X-ray emission comes from close to the source and the radio emission is originating in a more extended region. However, it is probable that the source was not in the canonical low-hard state during these radio observations and this may affect the behaviour of the source as well. At intermediate X-ray luminosities the source makes the transition from the radio underluminous sources in the direction of the relation traced by the 'standard' correlation similar to what has been reported for H1743-322. (abridged)
Analyzing exoplanets detected by radial velocity or transit observations, we determine the multiplicity of exoplanet host stars in order to study the influence of a stellar companion on the properties of planet candidates. Matching the host stars of exoplanet candidates detected by radial velocity or transit observations with online multiplicity catalogs in addition to a literature search, 57 exoplanet host stars are identified having a stellar companion. The resulting multiplicity rate of at least 12 percent for exoplanet host stars is about four times smaller than the multiplicity of solar like stars in general. The mass and the number of planets in stellar multiple systems depend on the separation between their host star and its nearest stellar companion, e.g. the planetary mass decreases with an increasing stellar separation. We present an updated overview of exoplanet candidates in stellar multiple systems, including 15 new systems (compared to the latest summary from 2009).
Runaway massive stars are O- and B-type stars with high spatial velocities with respect to the interstellar medium. These stars can produce bowshocks in the surrounding gas. Bowshocks develop as arc-shaped structures, with bows pointing to the same direction as the stellar velocity, while the star moves supersonically through the interstellar gas. The piled-up shocked matter emits thermal radiation and a population of locally accelerated relativistic particles is expected to produce non-thermal emission over a wide range of energies. We aim to model the non-thermal radiation produced in these sources. Under some assumptions, we computed the non-thermal emission produced by the relativistic particles and the thermal radiation caused by free-free interactions, for O4I and O9I stars. We applied our model to Zeta Oph (HD 149757), an intensively studied massive star seen from the northern hemisphere. This star has spectral type O9.5V and is a well-known runaway. Spectral energy distributions of massive runaways are predicted for the whole electromagnetic spectrum. We conclude that the non-thermal radiation might be detectable at various energy bands for relatively nearby runaway stars, especially at high-energy gamma rays. Inverse Compton scattering with photons from the heated dust gives the most important contribution to the high-energy spectrum. This emission approaches Fermi sensitivities in the case of Zeta Oph.
We introduce numerical algorithms for initializing multidimensional simulations of stellar explosions with 1D stellar evolution models. The initial mapping from 1D profiles onto multidimensional grids can generate severe numerical artifacts, one of the most severe of which is the violation of conservation laws for physical quantities. We introduce a numerical scheme for mapping 1D spherically-symmetric data onto multidimensional meshes so that these physical quantities are conserved. We validate our scheme by porting a realistic 1D Lagrangian stellar profile to the new multidimensional Eulerian hydro code CASTRO. Our results show that all important features in the profiles are reproduced on the new grid and that conservation laws are enforced at all resolutions after mapping. We also introduce a numerical scheme for initializing multidimensional supernova simulations with realistic perturbations predicted by 1D stellar evolution models. Instead of seeding 3D stellar profiles with random perturbations, we imprint them with velocity perturbations that reproduce the Kolmogorov energy spectrum expected for highly turbulent convective regions in stars. Our models return Kolmogorov energy spectra and vortex structures like those in turbulent flows before the modes become nonlinear. Finally, we describe approaches to determining the resolution for simulations required to capture fluid instabilities and nuclear burning. Our algorithms are applicable to multidimensional simulations besides stellar explosions that range from astrophysics to cosmology.
We examine the duty cycle and the history of star formation (SFH) for high-redshift galaxies at z \geq 6 using cosmological hydrodynamic simulations. We find that, even though individual galaxies have bursty SFH, the averaged SFH between z ~ 15 to z = 6 can be characterized well by either an exponentially increasing functional form with characteristic time-scales of 70 Myr to 200 Myr for galaxies with stellar masses Ms \sim 10^6 M\odot to > 10^10 M\odot respectively, or by a simple power-law form which exhibits a similar mass dependent time-scales. Using the SFH of individual galaxies, we measure the duty cycle of star formation (DC_SFH); i.e., the fraction of time a galaxy of a particular mass spends above a star formation rate (SFR) threshold which would make it observable to the Hubble Space Telescope (HST) during a given epoch. We also examine the fraction of galaxies at a given redshift that are brighter than a rest-frame UV magnitude Muv, which is sufficient enough to make them observable (DC_Muv). We find that both DC_SFH and DC_Muv make a sharp transition from zero (for galaxies with Ms \leq 10^7 M\odot) to unity (for Ms > 10^9 M\odot). The measured duty cycle is also manifested in the intrinsic scatter in the Ms-SFR relationship (\sim 1 dex) and Ms-Muv relationship (\Delta Muv \sim \pm 1 mag). We provide analytic fits to the DC as a function of Ms using a sigmoid function, which can be used in semi-analytic models of galaxy formation. We consider the effects of duty cycle to the observational estimate of galaxy stellar mass functions (GSMF) and the star formation rate density (SFRD), and find that it results in a much shallower low-mass end slopes of the GSMF and a reduction of \geq 70% of our intrinsic SFRD, making our simulation results more compatible with observational estimates.
We report the initial results from an ongoing multi-year spectroscopic survey of novae in M33. The survey resulted in the spectroscopic classification of six novae (M33N 2006-09a, 2007-09a, 2009-01a, 2010-10a, 2010-11a, and 2011-12a) and a determination of rates of decline (t_2 times) for four of them (2006-09a, 2007-09a, 2009-01a, and 2010-10a). When these data are combined with existing spectroscopic data for two additional M33 novae (2003-09a and 2008-02a) we find that 5 of the 8 novae with available spectroscopic class appear to be members of either the He/N or Fe IIb (hybrid) classes, with only two clear members of the Fe II spectroscopic class. This initial finding is very different from what would be expected based on the results for M31 and the Galaxy where Fe II novae dominate, and the He/N and Fe IIb classes together make up only ~20% of the total. It is plausible that the increased fraction of He/N and Fe IIb novae observed in M33 thus far may be the result of the younger stellar population that dominates this galaxy, which is expected to produce novae that harbor generally more massive white dwarfs than those typically associated with novae in M31 or the Milky Way.
With mid-IR and near-IR long-baseline interferometers, we are now mapping the radial distribution of the dusty accreting material in AGNs at sub-pc scales. We currently focus on Type 1 AGNs, where the innermost region is unobscured and its intrinsic structure can be studied directly. As a first systematic study of Type 1s, we obtained mid-/near-IR data for small samples over ~3-4 orders of magnitudes in UV luminosity L of the central engine. Here we effectively trace the structure by observing dust grains that are radiatively heated by the central engine. Consistent with a naive expectation for such dust grains, the dust sublimation radius R_in is in fact empirically known to be scaling with L^1/2 from the near-IR reverberation measurements, and this is also supported by our near-IR interferometry. Utilizing this empirical relationship, we normalize the radial extent by R_in and eliminate the simple L^1/2 scaling for a direct comparison over the samples. We then find that, in the mid-IR, the overall size in units of R_in seems to become more compact in higher luminosity sources. More specifically, the mid-IR brightness distribution is rather well described by a power-law, and this power-law becomes steeper in higher luminosity objects. The near-IR flux does not seem to be a simple inward extrapolation of the mid-IR power-law component toward shorter wavelengths, but it rather comes from a little distinct brightness concentration at the inner rim region of the dust distribution. Its structure is not well constrained yet, but there is tentative evidence that this inner near-IR-emitting structure has a steeper radial distribution in jet-launching objects. All these should be scrutinized with further observations.
We have computed with a fine time grid the evolution of the elemental abundances of He, C, N and O ejected by a massive (M$=10^{6}$\,\Msun) coeval stellar cluster with a Salpeter initial mass function (IMF) over a range of initial abundances. Our computations incorporate the mass loss from massive stars (M >30 Msun) during their wind phase including the Wolf-Rayet phase plus the ejecta from the supernova events. We find that during the Wolf-Rayet phase (t <5 Myr) the cluster ejecta composition suddenly becomes vastly over-solar in He, C, N and O and for all initial abundances with the C and O abundance in the cluster ejecta reaching over 50 times the solar value.
We show that the experiential $E_{\rm p}-L$, $\Gamma-L$, $E_{\rm p}-\Gamma$ and ${\rm \bar{\eta}_\gamma}-E_{\rm p}$ correlations (where $L$ is the luminosity of the prompt emission, $E_{\rm p}$ is the spectral peak energy, $\Gamma$ is the bulk Lorentz factor and $\bar{\eta}_\gamma$ is the emission efficiency of Gamma-ray Burst (GRB)) are well consistent with the relations between the resembling parameters predicted in the photospheric radiation model of the prompt emission of Gamma-ray Bursts. The time-resolved thermal radiation of GRB 090902B does follow the $E_{\rm p}-L$ and $\Gamma-L$ correlations. A reliable interpretation of the four correlations in alternative models is still lack. These facts may point towards a photospheric origin of prompt emission of some Gamma-ray Bursts.
In this paper we examine the effect of the formation and evolution of the disk galaxy on the distribution of dark halo matter. We have made simulations of isolated dark matter (DM) halo and two component (DM + baryons). N-body technique was used for stellar and DM particles and TVD MUSCL scheme for gas-dynamic simulations. The simulations include the processes of star formation, stellar feedback, heating and cooling of the interstellar medium. The results of numerical experiments with high spatial resolution let us to conclude in two main findings. First, accounting of star formation and supernova feedback resolves the so-called problem of cusp in distribution of dark matter predicted by cosmological simulations. Second, the interaction of dark matter with dynamic substructures of stellar and gaseous galactic disk (e.g., spiral waves, bar) has an impact on the shape of the dark halo. In particular, the in-plane distribution of dark matter is more symmetric in runs, where the baryonic component was taken into account.
Phenomenologically, two classes of GRBs (long/soft vs. short/hard) are identified based on their gamma-ray properties. The boundary between the two classes is vague. Multi-wavelength observations lead to identification of two types of GRB progenitor: one related to massive stars (Type II), and another related to compact stars (Type I). Evidence suggests that the majority of long GRBs belong to Type II, while at least the majority of nearby short GRBs belong to Type I. Nonetheless, counter examples do exist. Both long-duration Type I and short-duration Type II GRBs have been observed. In this talk, I review the complications in GRB classification and efforts in diagnosing GRB progenitor based on multiple observational criteria. In particular, I raise the caution to readily accept that all short/hard GRBs detected by BATSE are due to compact star mergers. Finally, I propose to introduce "amplitude" as the third dimension (besides "duration" and "hardness") to quantify burst properties, and point out that the "tip-of-iceberg" effect may introduce confusion in defining the physical category of GRBs, especially for low-amplitude, high-redshift GRBs.
Using stacks of Ly-a images of 2128 Ly-a emitters (LAEs) and 24 protocluster UV-selected galaxies (LBGs) at z=3.1, we examine the surface brightness profiles of Ly-a haloes around high-z galaxies as a function of environment and UV luminosity. We find that the slopes of the Ly-a radial profiles become flatter as the Mpc-scale LAE surface densities increase, but they are almost independent of the central UV luminosities. The characteristic exponential scale lengths of the Ly-a haloes appear to be proportional to the square of the LAE surface densities (r(Lya) \propto Sigma(LAE)^2). Including the diffuse, extended Ly-a haloes, the rest-frame Ly-a equivalent width of the LAEs in the densest regions approaches EW_0(Lya) ~ 200 A, the maximum value expected for young (< 10^7 yr) galaxies. This suggests that Ly-a photons formed via shock compression by gas outflows or cooling radiation by gravitational gas inflows may partly contribute to illuminate the Ly-a haloes; however, most of their Ly-a luminosity can be explained by photo-ionisation by ionising photons or scattering of Ly-a photons produced in HII regions in and around the central galaxies. Regardless of the source of Ly-a photons, if the Ly-a haloes trace the overall gaseous structure following the dark matter distributions, it is not surprising that the Ly-a spatial extents depend more strongly on the surrounding Mpc-scale environment than on the activities of the central galaxies.
We aimed to map the jets and outflows from the Serpens South star forming region and find an empirical relationship between the magnetic field and outflow orientation. Near-infrared H2 v=1-0 S(1) 2.122{\mu}m -line imaging of the \sim 30'-long filamentary shaped Serpens South star forming region was carried out. K s broadband imaging of the same region was used for continuum subraction. Candidate driving sources of the mapped jets/outflows are identified from the list of known protostars and young stars in this region, which was derived from studies using recent Spitzer and Herschel telescope observations. 14 Molecular Hydrogen emission-line objects(MHOs) are identified using our continuum-subtracted images. They are found to constitute ten individual flows. Out of these, nine flows are located in the lower-half(southern) part of the Serpens South filament, and one flow is located at the northern tip of the filament. Four flows are driven by well-identified Class 0 protostars, while the remaining six flows are driven by candidate protostars mostly in the Class I stage, based on the Spitzer and Herschel observations. The orientation of the outflows is systematically perpendicular to the direction of the near-infrared polarization vector, recently published in the literature. No significant correlation was observed between the orientation of the flows and the axis of the filamentary cloud.
The primary focus of this paper is on the particle acceleration mechanism in solar coronal 3D reconnection null-point regions. Starting from a potential field extrapolation of a SOHO magnetogram taken on 2002 November 16, we first performed MHD simulations with horizontal motions observed by SOHO applied to the photospheric boundary of the computational box. After a build-up of electric current in the fan-plane of the null-point, a sub-section of the evolved MHD data was used as initial and boundary conditions for a kinetic particle-in-cell model of the plasma. We find that sub-relativistic electron acceleration is mainly driven by a systematic electric field in the current sheet. A non-thermal population of electrons with a power-law distribution in energy forms, featuring a power-law index of about -1.75. This work provides a first step towards bridging the gap between macroscopic scales on the order of hundreds of Mm and kinetic scales on the order of cm in the solar corona, and explains how to achieve such a cross-scale coupling by utilizing either physical modifications or (equivalent) modifications of the constants of nature. With their exceptionally high resolution --- up to 135 billion particles and 3.5 billion grid cells of size 17.5\,km --- these simulations offer a new opportunity to study particle acceleration in solar-like settings.
We have selelcted 556 Red Horizontal Branch (RHB) stars along the streams of the Sagittarius dwarf galaxy (Sgr) from SDSS DR7 spectroscopic data using a theoretical model. The metallicity and \alpha-elements distributions are investigated for stars in the Sgr streams and for Galactic stars at the same locations. We find that the Sgr stars have two peaks in the metallicity distribution while the Galactic stars have a more prominent metal-poor peak. Meanwhile, [\alpha/Fe] ratios of the Sgr stars are lower than those of the Galactic stars. Among the Sgr stars, we find a difference in the metallicity distribution between the leading and trailing arms of the Sgr tidal tails. The metallicity and [\alpha/Fe] distribution of the leading arm is similar to that of the Galaxy. The trailing arm is composed mainly of a metal rich component and [\alpha/Fe] is obviously lower than that of the Galactic stars. The metallicity gradient is -(1.8 \pm 0.3)\times10^{-3} dex degree^{-1} in the first wrap of the trailing arm and -(1.5 \pm 0.4)\times10^{-3} dex degree^{-1} in the first wrap of the leading arm. No significant gradient exists along the second wraps of the leading or trailing arms. It seems that the Sgr dwarf galaxy initially lost the metal poor component in the second wrap (older) arms due to the tidal force of our Galaxy and then the metal rich component is disrupted in the first wrap (younger) arms. Finally, we found that the velocity dispersion of the trailing arm from 88^\circ<\Lambda_\odot<112^\circ is \sigma = 9.808 \pm 1.0 km s^{-1}, which is consistent with previous work in the literature.
Observations of the cyclotron resonance scattering feature in the X-ray spectrum of GX 301-2 suggest that the surface field of the neutron star is B_CRSF ~ 4 x 10^{12}G. The same value has been derived in modelling the rapid spin-up episodes in terms of the Keplerian disk accretion scenario. However, the spin-down rate observed during the spin-down trends significantly exceeds the value expected in currently used spin-evolution scenarios. This indicates that either the surface field of the star exceeds 50 x B_CRSF, or a currently used accretion scenario is incomplete. We show that the above discrepancy can be avoided if the accreting material is magnetized. The magnetic pressure in the accretion flow increases more rapidly than its ram pressure and, under certain conditions, significantly affects the accretion picture. The spin-down torque applied to the neutron star in this case is larger than that evaluated within a non-magnetized accretion scenario. We find that the observed spin evolution of the pulsar can be explained in terms of the magnetically controlled accretion flow scenario provided the surface field of the neutron star is ~ B_CRSF.
Both observations and numerical simulations show that stellar convective motions are composed of semi-regular flows of convective rolling cells and the fully developed turbulence. Although the convective rolling cells are crucial for the properties of the stellar convection that transports heat and mixes materials in the stellar interior, their contributions have not been included in turbulent convection models proposed up to now. We simplify the structure of the convective rolling cells as a cellular pattern moving circle by circle with different angular velocities around the center, estimating their typical size by solving approximately for the temperature difference over the stationary temperature background and their average shear of velocity by evaluating approximately their kinetic energy transformed by themselves working as thermal engines from the heat involved in the convective rolling cells. We obtain a steady state solution in the fully local equilibrium which is similar to what is obtained in the standard mixing-length theory, by applying such model assumptions to the standard $k$-$\varepsilon$ model and properly choosing the model parameter $c_{\varepsilon 3}$. Accordingly, we propose a $k$-$\omega$ model to include the transport effect of turbulence in stars. Preliminary results of their applications to the sun and other stars with different masses and in different evolutionary stages show good agreements with results of the standard mixing-length theory and results of numerical simulations for the stellar convection.
We study the cosmological information contained in the Minkowski Functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local type primordial non-Gaussianity parameter f_NL. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing, to generate 100 independent convergence maps of 25 deg^2 field-of-view for f_NL = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude sigma_8. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam survey with the proposed survey area of 1500 deg^2, the primordial non-Gaussianity can be constrained with a level of f_NL ~ 80 and w ~ 0.036 by weak lensing MFs. If simply scaled by the effective survey area, a 20000 deg^2 lensing survey using Large Synoptic Survey Telescope will give constraints of f_NL ~ 25 and w ~ 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.
The BLLac object S4 0954+65 is one of the main targets of the Urumqi monitoring program for IntraDay Variable (IDV) sources. Between August 2005 and December 2009, the source was included in 41 observing sessions, carried out at a frequency of 4.8 GHz. The time analysis of the collected light curves, performed through both a structure function analysis and a specifically developed wavelet-based algorithm, disclosed the existence of an annual cycle in the variability timescales, suggesting a fundamental contribution of interstellar scintillation to the IDV pattern of the source. The combined use of the two analysis methods also revealed a dramatic change in the variability characteristics of the source between February and March 2008, at the starting time of a strong outburst phase. The analysis' results suggest that the flaring state of the source coincides with the appearance of multiple timescales in its light curves, indicating that changes in the structure of the relativistically moving emitting region may strongly influence the variability observed on IDV timescales.
Context. SV Psc is an asymptotic giant branch (AGB) star surrounded by an
oxygen-rich dust envelope. The mm-CO line profile of the object's outflow shows
a clear double-component structure. Because of the high angular resolution,
mid-IR interferometry may give strong constraints on the origin of this
composite profile.
Aims. The aim of this work is to investigate the morphology of the
environment around SV Psc using high-angular resolution interferometry
observations in the mid-IR with the Very Large Telescope MID-infrared
Interferometric instrument (VLTI/MIDI).
Methods. Interferometric data in the N-band taken at different baseline
lengths (ranging from 32-64 m) and position angles (73- 142{\deg}) allow a
study of the morphology of the circumstellar environment close to the star. The
data are interpreted on the basis of 2-dimensional, chromatic geometrical
models using the fitting software tool GEM-FIND developed for this purpose.
Results. The results favor two scenarios: (i) the presence of a highly
inclined, optically thin, dusty disk surrounding the central star; (ii) the
presence of an unresolved binary companion at a separation of 13.7 AU and a
position angle of 121.8{\deg} NE. The derived orbital period of the binary is
38.1 yr. This detection is in good agreement with hydrodynamic simulations
showing that a close companion could be responsible for the entrainment of the
gas and dust into a circumbinary structure.
The small-scale CMB temperature we observe on the sky is modulated by perturbations that were super-horizon at recombination, giving differential focussing and lensing that generate a non-zero bispectrum even for single-field inflation where local physics is identical. Understanding this signal is important for primordial non-Gaussianity studies and also parameter constraints from the CMB lensing bispectrum signal. Because of cancellations individual effects can appear larger or smaller than they are in total, so a full analysis may be required to avoid biases. I relate angular scales on the sky to physical scales at recombination using the optical equations, and give full-sky results for the large-scale adiabatic temperature bispectrum from Ricci focussing (expansion of the ray bundle), Weyl lensing (convergence and shear), and temperature redshift modulations of small-scale power. The delta N expansion of the beam is described by the constant temperature 3-curvature, and gives a nearly-observable version of the consistency relation prediction from single-field inflation. I give approximate arguments to quantify the likely importance of dynamical effects, and argue that they can be neglected for modulation scales l <~ 100, which is sufficient for lensing studies and also allows robust tests of local primordial non-Gaussianity using only the large-scale modulation modes. For accurate numerical results early and late-time ISW effects must be accounted for, though I confirm that the late-time non-linear Rees-Sciama contribution is negligible compared to other more important complications. The total corresponds to f_NL ~ 7 for Planck-like temperature constraints and f_NL ~ 11 for cosmic-variance limited data to lmax=2000. Temperature lensing bispectrum estimates are affected at the 0.2 sigma level by Ricci focussing, and up to 0.5 sigma with polarization.
The ANTARES detector is the largest deep sea underwater neutrino telescope in operation. The apparatus comprises a matrix of 885 photomultiplier tubes (PMTs) which detect the Cherenkov light emitted by the charged leptons produced in the charged current interactions of high energy neutrinos with the matter inside or near the detector. Reconstruction of the muon track and energy can be achieved using the time, position and charge information of the hits arriving to the PMTs. A good calibration of the detector is necessary in order to ensure its optimal performance. This contribution reviews the different calibration systems and methods developed by the ANTARES Collaboration.
A variety of periodic phenomena have been observed in conjunction with large solar jets. We aim to find further evidence for {(quasi-)}periodic behaviour in solar jets and determine what the periodic behaviour can tell us about the excitation mechanism and formation process of the large solar jet. Using the 304 {\AA} (He-II), 171 {\AA} (Fe IX), 193 {\AA} (Fe XII/XXIV) and 131 {\AA} (Fe VIII/XXI) filters on-board the Solar Dynamic Observatory (SDO) Atmospheric Imaging Assembly (AIA), we investigate the intensity oscillations associated with a solar jet. Evidence is provided for multiple magnetic reconnection events occurring between a pre-twisted, closed field and open field lines. Components of the jet are seen in multiple SDO/AIA filters covering a wide range of temperatures, suggesting the jet can be classified as a blowout jet. Two bright, elongated features are observed to be co-spatial with the large jet, appearing at the jet's footpoints. Investigation of these features reveal they are defined by multiple plasma ejections. The ejecta display (quasi-)periodic behaviour on timescales of 50 s and have rise velocities of 40-150 km\,s$^{-1}$ along the open field lines. Due to the suggestion that the large jet is reconnection-driven and the observed properties of the ejecta, we further propose that these ejecta events are similar to type-II spicules. The bright features also display (quasi)-periodic intensity perturbations on the timescale of 300 s. Possible explanations for the existence of the (quasi-)periodic perturbations in terms of jet dynamics and the response of the transition region are discussed.
Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles' aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.
Fomalhaut is a young, nearby star that is suspected to harbor an infant
planetary system, interspersed with one or more belts of dusty debris. We
present far-infrared images obtained with the Herschel Space Observatory with
an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and
500 micrometer. The images show the main debris belt in great detail. Even at
high spatial resolution, the belt appears smooth. The region in between the
belt and the central star is not devoid of material; thermal emission is
observed here as well. Also at the location of the star, excess emission is
detected.
We use a dynamical model together with radiative-transfer tools to derive the
parameters of the debris disk. We include detailed models of the interaction of
the dust grains with radiation, for both the radiation pressure and the
temperature determination. Comparing these models to the spatially resolved
temperature information contained in the images allows us to place strong
constraints on the presence of grains that will be blown out of the system by
radiation pressure. We use this to derive the dynamical parameters of the
system.
The appearance of the belt points towards a remarkably active system in which
dust grains are produced at a very high rate by a collisional cascade in a
narrow region filled with dynamically excited planetesimals. Dust particles
with sizes below the blow-out size are abundantly present. The equivalent of
2000 one-km-sized comets are destroyed every day, out of a cometary reservoir
amounting to 110 Earth masses. From comparison of their scattering and thermal
properties, we find evidence that the dust grains are fluffy aggregates, which
indicates a cometary origin. The excess emission at the location of the star
may be produced by hot dust with a range of temperatures, but may also be due
to gaseous free-free emission from a stellar wind.
By actively distorting the Cosmic Microwave Background (CMB) over our past light cone, cosmic strings are unavoidable sources of non-Gaussianity. Developing optimal estimators able to disambiguate a string signal from the primordial type of non-Gaussianity requires calibration over synthetic full sky CMB maps, which till now had been numerically unachievable at the resolution of modern experiments. In this paper, we provide the first high resolution full sky CMB map of the temperature anisotropies induced by a network of cosmic strings since the recombination. The map has about 200 million sub-arcminute pixels in the healpix format which is the standard in use for CMB analyses (Nside=4096). This premiere required about 800,000 cpu hours; it has been generated by using a massively parallel ray tracing method piercing through a thousands of state of art Nambu-Goto cosmic string numerical simulations which pave the comoving volume between the observer and the last scattering surface. We explicitly show how this map corrects previous results derived in the flat sky approximation, while remaining completely compatible at the smallest scales.
We study the angular bispectrum of local type arising from the (possibly correlated) combination of a primordial adiabatic mode with an isocurvature one. Generically, this bispectrum can be decomposed into six elementary bispectra. We estimate how precisely CMB data, including polarization, can enable us to measure or constrain the six corresponding amplitudes, considering separately the four types of isocurvature modes (CDM, baryon, neutrino density, neutrino velocity). Finally, we discuss how the model-independent constraints on the bispectrum can be combined to get constraints on the parameters of multiple-field inflation models.
The Galactic black-hole candidate GX 339-4 exhibited several outbursts at
regular intervals of \sim 2-3 years in the Rossi X-ray Timing Explorer (RXTE)
era. After remaining in an almost quiescent state for 3 years, it again became
X-ray active in 2010 January, continuing to be so over the next \sim 14 months.
We study the timing and spectral properties of the black hole candidate (BHC)
during its recent outburst using RXTE PCA data, starting from 2010 January 12
to 2011 March 6.
Our study provides a comprehensive understanding of the mass accretion
processes and properties of the accretion disk of the black hole candidate. The
PCA spectra of 2.5-25 keV are mainly fitted with a combination of two
components, namely, a disk black body and a power-law. The entire outburst as
observed by RXTE, is divided into 4 spectral states, namely, hard,
hard-intermediate, soft-intermediate, and soft. Quasi-periodic oscillations
(QPOs) were found in 3 out of the 4 states, namely hard, hard-intermediate, and
soft-intermediate. The QPO frequencies increase monotonically from 0.102 Hz to
5.692 Hz in the rising phase of the outburst, while during the declining phase
QPO frequencies decrease monotonically from 6.420 to 1.149 Hz.
The recent outburst of GX 339-4 gives us an opportunity to understand the
evolution of the two-component accretion rates starting from the onset to the
end of the outburst phase. We found that the QPO frequency variation could be
explained by the propagating oscillatory shock model (POS) and the hardness
versus intensity variation can be reproduced if we assume that higher viscosity
causes the conversion of a low angular momentum disk component into a Keplerian
component during the outburst phase. The decline phase starts because of the
reduction in the viscosity.
These lectures attempt to expose the most important ideas, which have been proposed to explain the formation of stars with particular emphasis on the formation of brown dwarfs and low-mass stars. We first describe the important physical processes which trigger the collapse of a self-gravitating piece of fluid and regulate the star formation rate in molecular clouds. Then we review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function paying particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. As it is yet unsettled whether the brown dwarfs form as low-mass stars, we present the theory of brown dwarfs based on disk fragmentation stressing all the uncertainties due to the radiative feedback and magnetic field. Finally, we describe the results of large scale simulations performed to explain the collapse and fragmentation of molecular clouds.
Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.
The ANTARES neutrino telescope was completed in 2008 with the installation of its twelfth line. Its scientific scope is very broad, but the two main goals are the observation of astrophysical sources and the indirect detection of dark matter. The latter is possible through neutrinos produced after the annihilation of WIMPs, which would accumulate in sources like the Sun, the Earth or the Galactic Centre. The neutralino, which arises in Supersymmetry models, is one of the most popular WIMP candidates. KK particles, which appear in Universal Extra Dimension models, are another one. Though in most models these annihilations would not directly produce neutrinos, they are expected from the decay of secondary particles. An important advantage of neutrino telescopes with respect to other indirect searches (like gamma rays or cosmic rays) is that a potential signal (for instance from the Sun) would be very clean, since no other astrophysical explanations could mimic it (like pulsars for the case of the positron excess seen by PAMELA). Moreover, the Galactic Centre is accessible for ANTARES, being in the Northern Hemisphere. In this talk I will present the results of the ANTARES telescope for dark matter searches, which include neutralino and KK particles.
Paper I presents the magnetic structure of a filament that developed in
active region (AR) NOAA 10781. In this paper we complement those results with
the velocities retrieved from Doppler shifts measured at the chromosphere and
the photosphere in the AR filament area. Various inversion methods with
different numbers of atmospheric components and different weighting schemes of
the Stokes profiles were used. The velocities were calibrated on an absolute
scale. A ubiquitous chromospheric downflow is found in the faculae surrounding
the filament, with an average velocity of 1.6 km/s. The filament region,
however, displays upflows in the photosphere on both days, when the linear
polarization (which samples the transverse component of the fields) is given
more weight in the inversions. The upflow speeds of the transverse fields in
the filament region average -0.15 km/s. In the chromosphere, the situation is
different for the two days of observation. On July 3, the chromospheric portion
of the filament is moving upwards as a whole with a mean speed of -0.24 km/s.
However, on July 5 only the section above an orphan penumbra shows localized
upflow patches, while the rest of the filament is dominated by the same
downflows observed elsewhere in the facular region. Photospheric supersonic
downflows that last for tens of minutes are detected below the filament, close
to the PIL.
The observed velocity pattern in this AR filament strongly suggests a
scenario where the transverse fields are mostly dominated by upflows. The
filament flux rope is seen to be emerging at all places and both heights, with
a few exceptions in the chromosphere. This happens within a surrounding facular
region that displays a generalized downflow in the chromosphere and localized
downflows of supersonic character at the photosphere. No large scale downflow
of transverse field lines is observed at the photosphere.
We present seven new transiting hot Jupiters from the WASP-South survey. The planets are all typical hot Jupiters orbiting stars from F4 to K0 with magnitudes of V = 10.3 to 12.5. The orbital periods are all in the range 3.9--4.6 d, the planetary masses range from 0.4--2.3 Mjup and the radii from 1.1--1.4 Mjup. In line with known hot Jupiters, the planetary densities range from Jupiter-like to inflated (rho = 0.13--1.07 rho_jup). We use the increasing numbers of known hot Jupiters to investigate the distribution of their orbital periods and the 3--4-d "pile-up".
We investigate the acceleration of charged particles (both electrons and protons) at collisionless shocks predicted to exist in the vicinity of solar flares. The existence of standing termination shocks has been examined by flare models and numerical simulations e.g., Shibata,Forbes. We study electron energization by numerically integrating the equations of motion of a large number of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large-scale magnetic fluctuations -- known to exist in plasmas and known to have important effects on the nonthermal electron acceleration at shocks -- are also included in our simulations. For the parameters characteristic of the flaring region, our calculations suggest that the termination shock formed in the reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a sufficient turbulence amplitude level ($\delta B^2/B_0^2 \sim 0.3$), about 10% of thermal test-particle electrons are accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy in the available time and the energy spectra are considerably steeper than that of the electrons for the parameters used in our simulations. Our results are qualitatively consistent with the observed hard X-ray emissions in solar flares.
We describe an FPGA based Front-End Data Acquisition Module (FEDAM) for implementing Time-over-Threshold (ToT) Time-to-Digital conversion (TDC) of pulses obtained from the COSMICi astroparticle telescope detector system photomultiplier tubes. The telescope system consists of a minimum of three scintillation detectors configured to detect particle airshowers likely initiated by Ultra High Energy Cosmic Rays (UHECR). The relative time delay of detection events between the detectors is used to estimate the angle of incidence of the shower. The FEDAM provides time-over-threshold measurements with a resolution of 2 ns. This allows determination of shower direction to an error of 0.035 (cos {\theta})-1 radians where {\theta} is the angle between the baseline axis through a pair of detectors and the plane representing the shower front.
We investigate the influence of the initial proto-galaxies over-densities and masses on their evolution, to understand whether the internal properties of the proto-galactic haloes are sufficient to account for the varied properties of the galactic populations. By means of fully hydrodynamical N-body simulations performed with the code EvoL we produce twelve self-similar models of early-type galaxies of different initial masses and over-densities, following their evolution from z \geq 20 down to z \leq 1. The simulations include radiative cooling, star formation, stellar energy feedback, a reionizing photoheating background, and chemical enrichment of the ISM. We find a strong correlation between the initial properties of the proto-haloes and their star formation histories. Massive (10^13M\odot) haloes experience a single, intense burst of star formation (with rates \geq 10^3M\odot/yr) at early epochs, consistently with observations, with a less pronounced dependence on the initial over-density; intermediate mass (10^11M\odot) haloes histories strongly depend on their initial over-density, whereas small (10^9M\odot) haloes always have fragmented histories, resulting in multiple stellar populations, due to the "galactic breathing" phenomenon. The galaxy models have morphological, structural and photometric properties comparable to real galaxies, often closely matching the observed data; even though some disagreement is still there, likely a consequence of some numerical choices. We conclude that internal properties are essentially sufficient to explain many of the observed features of early type galaxies, particularly the complicated and different star formation histories shown by haloes of very different mass. In this picture, nature seems to play the dominant role, whereas nurture has a secondary importance.
In this chapter 2 of the e-book "Self-Organized Criticality Systems" we summarize the classical cellular automaton models, which consist of a statistical aspect that is universal to all SOC systems, and a physical aspect that depends on the physical definition of the observable. Then we derive some general analytical formulations of SOC processes, such as the exponential-growth SOC model and the fractal-diffusive SOC model, which also have universal validity for SOC processes, while specific applications to observations require additional physical scaling laws (e.g., for astrophysical or geophysical observations). Finally we discuss alternative SOC processes, SOC-related, or non-SOC processes, such as: self-organization (without criticality), forced SOC model, Brownian motion or classical diffusion, hyper-diffusion and Levy flight, nonextensive Tsallis entropy, turbulence, percolation, phase transitions, network systems, and chaotic systems. We synthesize a metrics that specifies which observational SOC properties are shared by these processes.
Weak-lensing shear estimates show a troublesome dependence on the apparent
brightness of the galaxies used to measure the ellipticity: In several studies,
the amplitude of the inferred shear falls sharply with decreasing source
significance. This dependence limits the overall ability of upcoming large
weak-lensing surveys to constrain cosmological parameters.
We seek to provide a concise overview of the impact of pixel noise on
weak-lensing measurements, covering the entire path from noisy images to shear
estimates. We show that there are at least three distinct layers, where pixel
noise not only obscures but biases the outcome of the measurements: 1) the
propagation of pixel noise to the non-linear observable ellipticity; 2) the
response of the shape-measurement methods to limited amount of information
extractable from noisy images; and 3) the reaction of shear estimation
statistics to the presence of noise and outliers in the measured ellipticities.
We identify and discuss several fundamental problems and show that each of
them is able to introduce biases in the range of a few tenths to a few percent
for galaxies with typical significance levels. Furthermore, all of these biases
do not only depend on the brightness of galaxies but also on their ellipticity,
with more elliptical galaxies often being harder to measure correctly. We also
discuss existing possibilities to mitigate and novel ideas to avoid the biases
induced by pixel noise. We present a new shear estimator that shows a more
robust performance for noisy ellipticity samples. Finally, we release the
open-source python code to predict and efficiently sample from the noisy
ellipticity distribution and the shear estimators used in this work at
https://github.com/pmelchior/epsnoise.
From optical spectroscopy of X-ray sources observed as part of ChaMP, we present redshifts and classifications for a total of 1569 Chandra sources from our targeted spectroscopic follow up using the FLWO, SAAO, WIYN, CTIO, KPNO, Magellan, MMT and Gemini telescopes, and from archival SDSS spectroscopy. We classify the optical counterparts as 50% BLAGN, 16% NELG, 14% ALG, and 20% stars. We detect QSOs out to z~5.5 and galaxies out to z~3. We have compiled extensive photometry from X-ray to radio bands. Together with our spectroscopic information, this enables us to derive detailed SEDs for our extragalactic sources. We fit a variety of templates to determine bolometric luminosities, and to constrain AGN and starburst components where both are present. While ~58% of X-ray Seyferts require a starburst event to fit observed photometry only 26% of the X-ray QSO population appear to have some kind of star formation contribution. This is significantly lower than for the Seyferts, especially if we take into account torus contamination at z>1 where the majority of our X-ray QSOs lie. In addition, we observe a rapid drop of the percentage of starburst contribution as X-ray luminosity increases. This is consistent with the quenching of star formation by powerful QSOs, as predicted by the merger model, or with a time lag between the peak of star formation and QSO activity. We have tested the hypothesis that there should be a strong connection between X-ray obscuration and star-formation but we do not find any association between X-ray column density and star formation rate both in the general population or the star-forming X-ray Seyferts. Our large compilation also allows us to report here the identification of 81 XBONG, 78 z>3 X-ray sources and 8 Type-2 QSO candidates. Also we have identified the highest redshift (z=5.4135) X-ray selected QSO with optical spectroscopy.
This document describes the exposure time calculator for the Wide-Field Infrared Survey Telescope (WFIRST) high-latitude survey. The calculator works in both imaging and spectroscopic modes. In addition to the standard ETC functions (e.g. background and S/N determination), the calculator integrates over the galaxy population and forecasts the density and redshift distribution of galaxy shapes usable for weak lensing (in imaging mode) and the detected emission lines (in spectroscopic mode). The source code is made available for public use.
We reconsider the strength of the electroweak phase transition (EWPT) in the inert doublet dark matter model, using a quantitatively accurate form for the one-loop finite temperature effective potential, taking into account relevant particle physics and dark matter constraints, focusing on a standard model Higgs mass near 126 GeV, and doing a full scan of the space of otherwise unconstrained couplings. We find that there is a significant (although fine-tuned) space of parameters for achieving an EWPT sufficiently strong for baryogenesis while satisfying the Xenon100 constraints from direct detection and obtaining the correct thermal relic density. We predict that the dark matter mass should be in the range 60-67 GeV, and we discuss possible LHC signatures of the charged and CP-odd Higgs bosons, including small suppression of the h -> 2 photon branching ratio.
We show that in imaginary time quantum metric fluctuations of empty space form a self-consistent De Sitter gravitational instanton that can be thought of as describing the tunneling from "nothing" into De Sitter space of real time (no cosmological constant or scalar fields are needed). The first time, this mechanism is activated to give birth to a flat inflationary Universe. The second time, it is turned on to complete cosmological evolution of the Universe after energy density of matter drops below the threshold (energy density of instantons). An accelerated expansion takes over after the scale factor exceeds this threshold, which marks the birth of dark energy at the redshift 1+z=1.44 and provides a possible solution to the "coincidence problem". The number of gravitons which tunneled into the Universe must be of the order of 10^122 to create the observational value of the Hubble constant. This number has nothing to do with vacuum energy, which is a possible solution to the "old cosmological constant problem". The emptying Universe should possibly complete its evolution by tunneling back to "nothing". After that, the entire scenario is repeated, and it can happen endlessly.
We generalize the previously proposed running vacuum energy model by including a term proportional to \dot{H}, in addition to the existing H^2 term. We show that the added degree of freedom is very constrained if both low redshift and high redshift data are taken into account. Best-fit models are undistinguishable from LCDM at the present time, but could be distinguished in the future with very accurate data at both low and high redshifts. We stress the formal analogy at the phenomenological level of the running vacuum models with recently proposed dark energy models based on the holographic or entropic point of view, where a combination of \dot{H} and H^2 term is also present. However those particular entropic formulations which do not have a constant term in the Friedmann equations are not viable. The presence of this term is necessary in order to allow for a transition from a decelerated to an accelerated expansion. In contrast, the running vacuum models, both the original and the generalized one introduced here contain this constant term in a more natural way. Finally, important conceptual issues common to all these models are emphasized.
We review recent literature on the connection between quantum entanglement and cosmology, with an emphasis on the context of expanding universes. We discuss recent theoretical results reporting on the production of entanglement in quantum fields due to the expansion of the underlying spacetime. We explore how these results are affected by the statistics of the field (bosonic or fermionic), the type of expansion (de Sitter or asymptotically stationary), and the coupling to spacetime curvature (conformal or minimal). We then consider the extraction of entanglement from a quantum field by coupling to local detectors and how this procedure can be used to distinguish curvature from heating by their entanglement signature. We review the role played by quantum fluctuations in the early universe in nucleating the formation of galaxies and other cosmic structures through their conversion into classical density anisotropies during and after inflation. We report on current literature attempting to account for this transition in a rigorous way and discuss the importance of entanglement and decoherence in this process. We conclude with some prospects for further theoretical and experimental research in this area. These include extensions of current theoretical efforts, possible future observational pursuits, and experimental analogues that emulate these cosmic effects in a laboratory setting.
We determine the neutrino spectra arising from low-mass (4-10 GeV) dark matter annihilating in the sun. We also determine the low-mass dark matter capture rates (element by element in the sun), assuming dark matter interacts either through elastic contact interactions, elastic long-range interactions, or inelastic contact interactions. These are the non-detector-specific data needed for determining the sensitivity of a neutrino detector to dark matter annihilating in the sun. As an application, we estimate the sensitivity of a one kiloton liquid scintillation neutrino detector (such as KamLAND) and LBNE (LAr-based) to low-mass dark matter with long-range interactions and compare this to the expected CDMS sensitivity. It is found that KamLAND's sensitivity can exceed that obtainable from the current CDMS data set by two orders of magnitude.
Links to: arXiv, form interface, find, astro-ph, recent, 1204, contact, help (Access key information)
It has been known for a long time that the satellite galaxies of the Milky Way (MW) show a significant amount of phase-space correlation, they are distributed in a highly inclined Disc of Satellites (DoS). We have extended the previous studies on the DoS by analysing for the first time the orientations of streams of stars and gas, and the distributions of globular clusters within the halo of the MW. It is shown that the spatial distribution of MW globular clusters classified as young halo clusters (YH GC) is very similar to the DoS, while 7 of the 14 analysed streams align with the DoS. The probability to find the observed clustering of streams is only 0.3 per cent when assuming isotropy. The MW thus is surrounded by a vast polar structure (VPOS) of subsystems (satellite galaxies, globular clusters and streams), spreading from Galactocentric distances as small as 10 kpc out to 250 kpc. These findings demonstrate that a near-isotropic infall of cosmological sub-structure components onto the MW is essentially ruled out because a large number of infalling objects would have had to be highly correlated, to a degree not natural for dark matter sub-structures. The majority of satellites, streams and YH GCs had to be formed as a correlated population. This is possible in tidal tails consisting of material expelled from interacting galaxies. We discuss the tidal scenario for the formation of the VPOS, including successes and possible challenges. The potential consequences of the MW satellites being tidal dwarf galaxies are severe. If all the satellite galaxies and YH GCs have been formed in an encounter between the young MW and another gas-rich galaxy about 10-11 Gyr ago, then the MW does not have any luminous dark-matter substructures and the missing satellites problem becomes a catastrophic failure of the standard cosmological model.
Using spatially-resolved spectra obtained with the Low Resolution Imaging Spectrometer at the Keck I telescope, we investigate the nature of ionizing sources and kinematic properties of emission-line gas in a LINER galaxy SDSS J091628.05+420818.7, which is a nearby (z = 0.0241) and bright (M_r = -20.2) early-type galaxy. After subtracting stellar absorption features using a combination of simple stellar population models, we measure the flux, line-of-sight velocity, and velocity dispersion of four emission lines, i.e., H{\alpha}, H{\beta}, [O III] {\lambda}5007, and [N II] {\lambda}6584, to study radial change of emission-line fluxes and velocities. Compared to the point-spread-function of the observation, the emission-line region is slightly extended but comparable to the seeing size. The central concentration of emission-line gas suggests that ionization is triggered by a nuclear source, excluding old stellar population as ionizing sources. We find that emission-line gas is counter-rotating with respect to stellar component and that the [O III] {\lambda}5007 line is blueshifted compared to other emission lines, possibly due to an outflow.
We examine the \alpha-element abundance ratio, [\alpha/Fe], of 5620 stars, observed by the Sloan Extension for Galactic Understanding and Exploration survey in the region 6 kpc < R < 16 kpc, 0.15 kpc < |Z| < 1.5 kpc, as a function of Galactocentric radius R and distance from the Galactic plane |Z|. Our results show that the high-\alpha\ thick disk population has a short scale length (L_thick ~ 1.8 kpc) compared to the low-\alpha population, which is typically associated with the thin disk. We find that the fraction of high-\alpha\ stars in the inner disk increases at large |Z|, and that high-\alpha\ stars lag in rotation compared to low-\alpha\ stars. In contrast, the fraction of high-\alpha\ stars in the outer disk is low at all |Z|, and high- and low-\alpha\ stars have similar rotational velocities up to 1.5 kpc from the plane. We interpret these results to indicate that different processes were responsible for the high-\alpha\ populations in the inner and outer disk. The high-\alpha\ population in the inner disk has a short scale length and large scale height, consistent with a scenario in which the thick disk forms during an early gas-rich accretion phase. Stars far from the plane in the outer disk may have reached their current locations through heating by minor mergers. The lack of high-\alpha\ stars at large R and |Z| also places strict constraints on the strength of radial migration via transient spiral structure.
We carry out a time-series analysis of the combined data from three experiments measuring the cosmic muon flux at the Gran Sasso laboratory, at a depth of 3800 m.w.e. These data, taken by the MACRO, LVD and Borexino experiments, span a period of over 20 years, and correspond to muons with a threshold energy, at sea level, of around 1.3 TeV. We compare the best-fit period and phase of the full muon data set with the combined DAMA/NaI and DAMA/LIBRA data, which spans the same time period, as a test of the hypothesis that the cosmic ray muon flux is responsible for the annual modulation detected by DAMA. We find in the muon data a large-amplitude fluctuation with a period of around one year, and a phase that is incompatible with that of the DAMA modulation at 5.2 sigmas. Aside from this annual variation, the muon data also contains a further significant modulation with a period between 10 and 11 years and a power well above the 99.9% C.L threshold for noise, whose phase corresponds well with the solar cycle: a surprising observation for such high energy muons. We see no corresponding long-period oscillation in the stratospheric temperature data.
We develop a method to extract the shape information of line profiles from discrete kinematic data. The Gauss-Hermite expansion, which is widely used to describe the line of sight velocity distributions extracted from absorption spectra of elliptical galaxies, is not readily applicable to samples of discrete stellar velocity measurements, accompanied by individual measurement errors and probabilities of membership. We introduce two parameter families of probability distributions describing symmetric and asymmetric distortions of the line profiles from Gaussianity. These are used as the basis of a maximum likelihood estimator to quantify the shape of the line profiles. Tests show that the method outperforms a Gauss-Hermite expansion for discrete data, with a lower limit for the relative gain of approx 2 for sample sizes N approx 800. To ensure that our methods can give reliable descriptions of the shape, we develop an efficient test to assess the statistical quality of the obtained fit. As an application, we turn our attention to the discrete velocity datasets of the dwarf spheroidals of the Milky Way. In Sculptor, the symmetric deviations are everywhere consistent with velocity distributions more peaked than Gaussian. This may suggest a radially biased orbital structure. In Fornax, there is an evolution in the symmetric deviations of the line profile from a peakier to more flat-topped distribution on moving outwards. This may be an indication of tangential anisotropy in the outer parts. Our methods are sensitive enough to detect evidence for velocity distributions more peaked than Gaussian in Carina and Sextans, which suggest radial anisotropy in the outer parts of these two galaxies. This is all consistent with a picture in which Fornax may have had a different evolutionary history to Sculptor, Carina and Sextans.
We have carefully selected a sample of 60 galaxies that reside in the deepest underdensities of geometrically identified voids within the SDSS. HI imaging of 55 galaxies with the WSRT reveals morphological and kinematic signatures of ongoing interactions and gas accretion. We probe a total volume of 485 Mpc^3 within the voids, with an angular resolution of 8 kpc at an average distance of 85 Mpc. We reach column density sensitivities of 5 x 10^19 cm^-2, corresponding to an HI mass limit of 3 x 10^8 M_sun. We detect HI in 41 galaxies, with total masses ranging from 1.7 x 10^8 to 5.5 x 10^9 M_sun. The upper limits on the 14 non-detections are not inconsistent with their luminosities, given their expected HI mass to light ratios. We find that the void galaxies are generally gas rich, low luminosity, blue disk galaxies, with optical and HI properties that are not unusual for their luminosity and morphology. The sample spans a range of absolute magnitudes (-16.1 > M_r > -20.4) and colors (0.06 < g-r < 0.87), and includes disk and irregular galaxies. We also identify three as early type galaxies, all of which are not detected in HI. All galaxies have stellar masses less than 3 x 10^10 M_sun, and many have kinematic and morphological signs of ongoing gas accretion, suggesting that the void galaxy population is still in the process of assembling. The small scale clustering in the void, within 600 kpc and 200 km/s, is similar to that in higher density regions, and we identify 18 HI rich neighboring galaxies in the voids. Most are within 100 kpc and 100 km/s of the targeted galaxy, and we find no significant population of HI rich low luminosity galaxies filling the voids, contrary to what is predicted by simulations.
The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets.
We use the Marcario Low Resolution Spectrograph (LRS) at the
Hobby-Eberly-Telescope (HET) to study the kinematics of pseudobulges and
classical bulges in the nearby universe. We present major-axis rotational
velocities, velocity dispersions, and h3 and h4 moments derived from
high-resolution (sigma ~ 39 km/s) spectra for 45 S0 to Sc galaxies; for 27 of
the galaxies we also present minor axis data. We combine our kinematics with
bulge-to-disk decompositions. We demonstrate for the first time that purely
kinematic diagnostics of the bulge dichotomy agree systematically with those
based on S\'ersic index. Low S\'ersic index bulges have both increased
rotational support (higher v/sigma values) and on average lower central
velocity dispersions. Furthermore, we confirm that the same correlation also
holds when visual morphologies are used to diagnose bulge type. The previously
noted trend of photometrically flattened bulges to have shallower velocity
dispersion profiles turns to be significant and systematic if the S\'ersic
index is used to distinguish between pseudobulges and classical bulges. The
correlation between h3 and v/sigma observed in elliptical galaxies is also
observed in intermediate type galaxies, irrespective of bulge type. Finally, we
present evidence for formerly undetected counter rotation in the two systems
NGC 3945 and NGC 4736.
Based on observations obtained with the Hobby-Eberly Telescope, which is a
joint project of the University of Texas at Austin, the Pennsylvania State
University, Stanford University, Ludwig-Maximilians-Universit\"at M\"unchen,
and Georg-August-Universit\"at G\"ottingen.
We use distant Blue Horizontal Branch stars with Galactocentric distances 16 < r/kpc < 48 as kinematic tracers of the Milky Way dark halo. We model the tracer density as an oblate, power-law embedded within a spherical power-law potential. Using a distribution function method, we estimate the overall power-law potential and the velocity anisotropy of the halo tracers. We measure the slope of the potential to be gamma ~ 0.4 and the overall mass within 50 kpc is ~ 4 x 10^11 M_sol. The tracer velocity anisotropy is radially biased with beta ~ 0.5, which is in good agreement with local solar neighbourhood studies. Our results provide an accurate outer circular velocity profile for the Milky Way and suggest a relatively high concentration dark matter halo (c_vir ~ 20).
We report hitherto unnoticed patterns in quasar light curves. We characterize segments of quasars' light curves with the slopes of the straight lines fit through them. These slopes appear to be directly related to the quasars' redshifts. Alternatively, using only global shifts in time and flux, we are able to find significant overlaps between the light curves of different pairs of quasars by fitting the ratio of their redshifts. We are then able to reliably determine the redshift of one quasar from another. This implies that one can use quasars as standard clocks, as we explicitly demonstrate by constructing two independent methods of finding the redshift of a quasar from its light curve.
The availability of asteroseismic constraints for a large sample of stars from the missions CoRoT and Kepler paves the way for various statistical studies of the seismic properties of stellar populations. In this paper, we evaluate the impact of rotation-induced mixing and thermohaline instability on the global asteroseismic parameters at different stages of the stellar evolution from the Zero Age Main Sequence to the Thermally Pulsating Asymptotic Giant Branch to distinguish stellar populations. We present a grid of stellar evolutionary models for four metallicities (Z = 0.0001, 0.002, 0.004, and 0.014) in the mass range between 0.85 to 6.0 Msun. The models are computed either with standard prescriptions or including both thermohaline convection and rotation-induced mixing. For the whole grid we provide the usual stellar parameters (luminosity, effective temperature, lifetimes, ...), together with the global seismic parameters, i.e. the large frequency separation and asymptotic relations, the frequency corresponding to the maximum oscillation power {\nu}_{max}, the maximal amplitude A_{max}, the asymptotic period spacing of g-modes, and different acoustic radii. We discuss the signature of rotation-induced mixing on the global asteroseismic quantities, that can be detected observationally. Thermohaline mixing whose effects can be identified by spectroscopic studies cannot be caracterized with the global seismic parameters studied here. But it is not excluded that individual mode frequencies or other well chosen asteroseismic quantities might help constraining this mixing.
We examine the rotation rates, sizes, and star formation (SF) efficiencies of a representative population of simulated disc galaxies extracted from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamic simulations. These simulations include efficient, but energetically feasible supernova feedback, but have not been tuned in any way to produce 'realistic' disc galaxies. Yet, they generate a large number of discs, without requiring extremely high resolution. Over the wide galaxy stellar mass range, 9.0 < log10[Mstar (Msun)] < 10.5, the simulations reproduce the observed Tully-Fisher relation, the rotation curves of disc galaxies in bins of stellar mass, the mass-size relation of disc galaxies, and the SF efficiencies of disc galaxies as inferred from stacked weak lensing and stacked satellite kinematics observations. At higher stellar masses, log10[Mstar (Msun)] > 10.6, the simulated galaxies are too concentrated and have too high SF efficiencies. We conjecture that this shortcoming reflects the neglect of feedback from accreting supermassive black holes in these simulations. We conclude that it is possible to generate a realistic and representative population of disc galaxies using standard numerical hydrodynamic techniques and a plausible implementation of the "subgrid" astrophysical processes thought to be relevant to galaxy formation.
In the present article, we present a galactic gravitational model of three degrees of freedom, in order to investigate and reveal the behavior of orbits in a binary quasar system. The two quasars are hosted in a pair of interacting disk galaxies. We study in detail the regular or chaotic character of motion, in two different cases: the time independent model in both 2D and 3D dynamical systems and the evolving 3D model. Our numerical calculations, indicate that a large fraction of orbits in the 2D system are chaotic in the time independent case. A careful analysis suggest that several Lindblad resonances are also responsible for the chaotic motion of stars in both hosts galaxies. In the time dependent system, we follow the evolution of 3D orbits in our dynamical model, as the two interacting host galaxies develop dense and massive quasars in their cores, by mass transportation from the disks to their nuclei. In this interesting case, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of 3D orbits depends, not only in the galactic interaction but also in the presence of the quasars in the galactic cores of the host galaxies. The outcomes derived from our dynamical model are compared with observational data. Some theoretical arguments to support the numerically derived outcomes are presented, both in 2D and 3D systems, while a comparison with similar earlier work is also made.
Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal' features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.)
Correlations are investigated of the outflow strength of quasars, as measured by the blueshift and asymmetry index (BAI) of the CIV line (Wang et al. 2011), with intensities and ratios of broad emission lines, based on composite quasar spectra built from the Sloan Digital Sky Survey. We find that most of the line ratios of other ions to CIV prominently increases with BAI. These behaviors can be well understood in the context of increasing metallicity with BAI. The strength of dominant coolant, CIV line, decreases and weak collisionally excited lines increase with gas metallicity as a result of the competition between different line coolants. Using SiIV+OIV]/CIV as an indicator of gas metallicity, we present, for the first time, a strong correlation between the metallicitiy and the outflow strength of quasars over a wide range of 1.7 to 6.9 times solar abundance. Our result implies that the metallicity plays an important role in the formation of quasar outflows, likely via affecting outflow acceleration. This effect may have a profound impact on galaxy evolution via momentum feedback and chemical enrichment.
We study the formation of structure in the Universe assuming that dark matter can be described by a scalar field $\tilde{\Phi}$ with a potential $V(\Phi)=-\mathfrak{m}^{2}\tilde{\Phi}^{2}/2+\lambda\tilde{\Phi}^4/4$. We derive the evolution equations of the scalar field in the linear regime of perturbations. We investigate the symmetry breaking and possibly a phase transition of this scalar field in the early Universe. At low temperatures, the scalar perturbations have an oscillating growing mode and therefore, this kind of dark matter could lead to the formation of gravitational structures. In order to study the nonlinear regime, we use the spherical collapse model and show that, in the quadratic potential limit, this kind of dark matter can form virialized structures. The main difference with the traditional Cold Dark Matter paradigm is that the formation of structure in the scalar field model can occur at earlier times. Thus, if the dark matter is of scalar field nature we expect to have large galaxies at high redshifts.
Frequencies of magnetic patch processes on supergranule boundary, namely flux emergence, splitting, merging, and cancellation, are investigated through an automatic detection. We use a set of line of sight magnetograms taken by the Solar Optical Telescope (SOT) on board Hinode satellite. We found 1636 positive patches and 1637 negative patches in the data set, whose time duration is 3.5 hours and field of view is 112" \times 112". Total numbers of magnetic processes are followed: 493 positive and 482 negative splittings, 536 positive and 535 negative mergings, 86 cancellations, and 3 emergences. Total numbers of emergence and cancellation are significantly smaller than those of splitting and merging. Further, frequency dependences of merging and splitting processes on flux content are investigated. Merging has a weak dependence on flux content only with a power- law index of 0.28. Timescale for splitting is found to be independent of parent flux content before splitting, which corresponds to \sim 33 minutes. It is also found that patches split into any flux contents with a same probability. This splitting has a power-law distribution of flux content with an index of -2 as a time independent solution. These results support that the frequency distribution of flux content in the analyzed flux range is rapidly maintained by merging and splitting, namely surface processes. We suggest a model for frequency distributions of cancellation and emergence based on this idea.
We investigate the nature of the inner accretion disk in the neutron star source GX 5-1 by making a detailed study of time lags between X-rays of different energies. Using the cross-correlation analysis, we found anti-correlated hard and soft time lags of the order of a few tens to a few hundred seconds and the corresponding intensity states were mostly the horizontal branch (HB) and upper normal branch (NB). The model independent and dependent spectral analysis showed that during these time lags the structure of accretion disk significantly varied. Both eastern and western approaches were used to unfold the X-ray continuum and systematic changes were observed in soft and hard spectral components. These changes along with a systematic shift in the frequency of quasi-periodic oscillations (QPOs) made it substantially evident that the geometry of the accretion disk is truncated. Simultaneous energy spectral and power density spectral study shows that the production of the horizontal branch oscillations (HBOs) are closely related to the Comptonizing region rather than the disk component in the accretion disk. We found that as the HBO frequency decreases from the hard apex to upper HB, the disk temperature increases along with an increase in the coronal temperature which is in sharp contrast with the changes found in black hole binaries where the decrease in QPO frequency is accompanied by a decrease in the disk temperature and a simultaneous increase in the coronal temperature. We discuss the results in the context of re-condensation of coronal material in the inner region of the disk.
With the Infrared Camera on board AKARI, we carried out near-infrared (2.5-5.0 micron) spectroscopy of the central kiloparsec region of the barred spiral galaxy, NGC1097, categorized as Seyfert 1 with a circumnuclear starburst ring. Our observations mapped the area of ~50"*10" with the resolution of ~5", covering about a half of the ring and the galactic center. As a result, we spatially resolve the starburst ring in the polycyclic aromatic hydrocarbon 3.3 micron, the aliphatic hydrocarbon 3.4-3.6 micron features, and the hydrogen Br alpha 4.05 micron emission. They exhibit spatial distributions significantly different from each other, indicating that the environments vary considerably around the ring. In particular, the aliphatic features are enhanced near the bar connecting the ring with the nucleus, where the structure of hydrocarbon grains seems to be relatively disordered. Near the center, the continuum emission and the CO/SiO absorption features are strong, which indicates that the environments inside the ring are dominated by old stellar populations. The near-infrared spectra do not show any evidence for the presence of nuclear activity.
ANTARES is the first neutrino telescope in the sea. It consists of a three-dimensional array of 885 photomultipliers to collect the Cherenkov light induced by relativistic muons produced in CC interactions of high energy neutrinos. One of the main scientific goals of the experiment is the search for dark matter. We present here the analysis of data taken during 2007 and 2008 to look for a WIMP signal in the Sun. WIMPs are one of the most popular scenarios to explain the dark matter content of the Universe. They would accumulate in massive objects like the Sun or the Galactic Center and their self-annihilation would produce (directly or indirectly) high energy neutrinos detectable by neutrino telescopes. Contrary to other indirect searches (like with gamma rays or positrons), the search for neutrinos in the Sun is free from other astrophysical contributions, so the interpretation of a potential signal in terms of dark matter is much more robust.
The purpose of this study is to analyze the dynamical role of a radiation field on the growth rate of the unstable Kelvin - Helmholtz (KH) perturbations. As a first step toward this purpose, the analyze is done in a general way, irrespective of applying the model to a specific astronomical system. The transition zone between the two layers of the fluid is ignored. Then, we perform a linear analysis and by imposing suitable boundary conditions and considering a radiation field, we obtain appropriate dispersion relation. Unstable modes are studied by solving the dispersion equation numerically, and then growth rates of them are obtained. By analyzing our dispersion relation, we show that for a wide range of the input parameters, the radiation field has a destabilizing effect on KH instability. In eruptions of the galaxies or supermassive stars, the radiation field is dynamically important and because of the enhanced KH growth rates in the presence of the radiation; these eruptions can inject more momentum and energy into their environment and excite more turbulent motions.
Generalised Teleparallel gravity, also referred to as f(T) gravity, has been recently proposed as an extended theory of gravitation able to give rise to an accelerated expansion in a matter only universe. The cosmic speed up is driven by an effective torsion fluid whose equation of state depend on the f(T) function entering the modified gravity Lagrangian. We focus on two particular choices for f(T) which share the nice property to emulate a phantom divide crossing as suggested by some recent data. We check their viability contrasting the predicted background dynamics to the Hubble diagram as traced by both Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs), the measurement of the rate expansion H(z), the Baryon Acoustic Oscillations (BAOs) at different redshifts, and the Cosmic Microwave Background Radiation (CMBR) distance priors. Both f(T) models turn out to be in very good agreement with this large dataset so that we also investigate whether it is possible to discriminate among them relying on the different growth factors.
We consider a coplanar system comprised of a massive central body (a star), a less massive secondary (a planet) on a circular orbit, and a test particle on a bound orbit exterior to that of the secondary. The gravitational pull exerted on the test particle by the secondary acts as a small perturbation, wherefore the trajectory of the particle can be described as an ellipse of a precessing perihelion. While the apsidal motion is defined overwhelmingly by the Newtonian portion of the secondary's gravity, the post-Newtonian portion, too, brings its tiny input. We explore whether this input may be of any astrophysical relevance in the next few decades. We demonstrate that the overall post-Newtonian input of the secondary's gravity can be split into two parts. One can be expressed via the orbital angular momentum of the secondary, another via its orbital radius. Despite some moderately large numerical factors showing up in the expressions for these two parts, the resulting post-Newtonian contributions from the secondary's gravity into the apsidal motion of the test particle turn out to be small enough to be neglected in the near-future measurements.
Motivated by recent discoveries of low-density super-Earths with short orbital periods, we have investigated in-situ accretion of H-He atmospheres on rocky bodies embedded in dissipating warm disks, by simulating quasi-static evolution of atmospheres that connect to the ambient disk. We have found that the atmospheric evolution has two distinctly different outcomes, depending on the rocky body's mass: While the atmospheres on massive rocky bodies undergo runaway disk-gas accretion, those on light rocky bodies undergo significant erosion during disk dispersal. In the atmospheric erosion, the heat content of the rocky body that was previously neglected plays an important role. We have also realized that the atmospheric mass is rather sensitive to disk temperature in the mass range of interest in this study. Our theory is applied to recently-detected super-Earths orbiting Kepler-11 to examine the possibility that the planets are rock-dominated ones with relatively thick H-He atmospheres. The application suggests that the in-situ formation of the relatively thick H-He atmospheres inferred by structure modeling is possible only under restricted conditions; namely, relatively slow disk dissipation and/or cool environments. This study demonstrates that low-density super-Earths provide important clues to understanding of planetary accretion and disk evolution.
We present high-resolution I-band imaging of the core of the globular cluster M15 obtained at the 2.5 m Nordic Optical Telescope with FastCam, a low readout noise L3CCD based instrument. Short exposure times (30 ms) were used to record 200000 images (512 x 512 pixels each) over a period of 2 hours 43 min. The lucky imaging technique was then applied to generate a final image of the cluster centre with FWHM ~ 0".1 and 13" x 13" FoV. We obtained a catalogue of objects in this region with a limiting magnitude of I=19.5. I-band photometry and astrometry are reported for 1181 stars. This is the deepest I-band observation of the M15 core at this spatial resolution. Simulations show that crowding is limiting the completeness of the catalogue. At shorter wavelengths, a similar number of objects has been reported using HST/WFPC observations of the same field. The cross-match with the available HST catalogues allowed us to produce colour-magnitude diagrams where we identify new Blue Straggler star candidates and previously known stars of this class.
We forecast the constraints on modified theories of gravity from the cosmic microwave background (CMB) anisotropies bispectrum that arises from correlations between lensing and the Integrated Sachs-Wolfe effect. In models of modified gravity the evolution of the metric potentials is generally altered and the contribution to the CMB bispectrum signal can differ significantly from the one expected in the standard cosmological model.We adopt a parametrised approach and focus on three different classes of models: Linder's growth index, Chameleon-type models and f(R) theories. We show that the constraints on the parameters of the models will significantly improve with future CMB bispectrum measurements.
Using a simple two-color selection based on $g$-, $z$-, and $K$-band photometry, we pick out 1609 star-forming galaxies (sgzKs) and 422 passively evolving galaxies (pgzKs) at z\sim2$ from a $K$-band-selected sample ($K_{\rm AB} < 22.0$) in an area of $\sim 0.44$ deg$^{2}$ of the All-wavelength Extended Groth Strip International Survey. The number counts of pgzKs\ in our sample turn over at $K_{\rm AB} \sim 21.0$, and both the number of faint and bright objects (including sgzKs\ and pgzKs) exceed the predictions of a recent semi-analytic model of galaxy formation, a more successful model is need to explain this diversity. We also find that the star formation rate (SFR) and specific SFR (sSFR) of sgzKs\ increases with redshift at all masses, implying that star-forming galaxies were much more active on average in the past. Moreover, the sSFR of massive galaxies is lower at all redshifts, suggesting that star formation contributes more to the mass growth of low-mass galaxies than to high-mass galaxies. From {\it Hubble Space Telescope} Wide Field Camera 3 near-infrared imaging data, we find that morphologies of $z\sim2$ galaxies not only have diffuse structures with lower $G$ and higher $M_{20}$ values, but also have single-object morphologies (higher $G$ and lower $M_{20}$), implying that there are morphological variety and different formation process for these galaxies at $z\sim2$. Finally, we also study the fraction of active galactic nuclei (AGNs) in the gzKs, 82 of 828 gzKs\ with four IRAC bands can be classified as AGNs ($\sim$ 10%). Most of these AGN candidates have $L_{\rm 0.5-10\ keV}>10^{41}\,\rm erg\,s^{-1}$.
Blue (metal-poor) globular clusters are observed to have half-light radii that are ~20% larger than their red (metal-rich) counterparts. The origin of this enhancement is not clear and differences in either the luminosity function or in the actual size of the clusters have been proposed. I analyze a set of dynamically self-consistent Monte Carlo globular cluster simulations to determine the origin of this enhancement. I find that my simulated blue clusters have larger half-light radii due to differences in the luminosity functions of metal-poor and metal-rich stars. I find that the blue clusters can also be physically larger, but only if they have a substantial number of black holes heating their central regions. In this case the difference between half-light radii is significantly larger than observed. I conclude that the observed difference in half-light radii between red and blue globular clusters is due to differences in their luminosity functions and that half-light radius is not a reliable proxy for cluster size.
The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.
A way to probe alternative theories of gravitation is to study if they could account for the structures of the universe. We then modified the well-known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. As an application, we simulate the evolution of spiral galaxies to probe alternative theories of gravitation whose weak field limits have a Yukawa-like gravitational potential. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation. It is worth stressing that the recipe given in the present study can be applied to any other alternative theory of gravitation in which the superposition principle is valid.
We study brightness variations in the double lensed quasar UM673 (Q0142-100) with the aim of measuring the time delay between its two images. In the paper we combine our previously published observational data of UM673 obtained during the 2003 - 2005 seasons at the Maidanak Observatory with archival and recently observed Maidanak and CTIO UM673 data. We analyze the V, R and I-band light curves of the A and B images of UM673, which cover ten observational seasons from August 2001 to November 2010. We also analyze the time evolution of the difference in magnitudes between images A and B of UM673 over more than ten years. We find that the quasar exhibits both short-term (with amplitude of \sim 0.1 mag in the R band) and high-amplitude (\sim 0.3 mag) long-term variability on timescales of about several months and several years, respectively. These brightness variations are used to constrain the time delay between the images of UM673. From cross-correlation analysis of the A and B quasar light curves and error analysis we measure the mean time delay and its error of 89 \pm11 days. Given the input time delay of 88 days, the most probable value of the delay that can be recovered from light curves with the same statistical properties as the observed R-band light curves of UM673 is 95{+5/-16}{+14/-29} days (68 and 95 % confidence intervals). Analysis of the V - I color variations and V, R and I-band magnitude differences of the quasar images does not show clear evidence of the microlensing variations between 1998 and 2010.
We report a list of groups consisting of dwarf galaxies only. The sample contains 126 objects, mainly combined in pairs. The most populated group contains six dwarf galaxies. The majority of systems considered reside in the low-density regions and evolve unaffected by massive galaxies. The characteristic sizes and velocity dispersions of groups are 30 kpc and 11 km/s, respectively. They resemble the associations of dwarf galaxies, but are more compact. On the whole, groups and associations form a continuous sequence. Alike the associations, our groups possess high mass-to-luminosity ratios, what is indicative of a large amount of dark matter present in these systems.
The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.
We present a comprehensive study of the periodic variations observed in OGLE I-band light curves of SMC Be/X-ray binaries, discovering new optical periodicities in 9 systems. We find that these periodicities derive from a number of mechanisms, notably disturbance of the decretion disk on the orbital period of the system, and aliased non-radial pulsations. We develop metrics that allow these mechanisms to be distinguished on the basis of the shape of the folded optical light curve, and use these metrics to categorise the periodicities present in \sim 50 SMC binary systems. We conclude that extreme care must be taken in the interpretation of the OGLE light curves since only around 30% of the periodicities present can be unambiguously attributed to orbital periods.
We used the N-body code of Hernquist and Ostriker (1992) to build a dozen cuspy ({\gamma}\approx 1) triaxial models of stellar systems through dissipationless collapses of initially spherical distributions of 10^6 particles. We chose four sets of initial conditions that resulted in models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. Within each set, three different seed numbers were selected for the random number generator used to create the initial conditions, so that the three models of each set are statistically equivalent. We checked the stability of our models using the values of their central densities and of their moments of inertia, which turned out to be very constant indeed. The changes of those values were all less than 3 per cent over one Hubble time and, moreover, we show that the most likely cause of those changes are relaxation effects in the numerical code. We computed the six Lyapunov exponents of nearly 5,000 orbits in each model in order to recognize regular, partially and fully chaotic orbits. All the models turned out to be highly chaotic, with less than 25 per cent of their orbits being regular. We conclude that it is quite possible to obtain cuspy triaxial stellar models that contain large fractions of chaotic orbits and are highly stable. The difficulty to build such models with the method of Schwarzschild (1979) should be attributed to the method itself and not to physical causes.
We perform idealised numerical simulations of magnetic buoyancy instabilities in a model of the solar tachocline. We introduce a simplified model of magnetic flux pumping in an upper layer (the convection zone), and study the effects of its inclusion on the evolution of buoyancy instabilities in a lower layer (the radiative interior). We study its effects on the instability of both a preconceived magnetic slab and of a shear-generated magnetic layer. In the former, we find that in the regime in which the downward pumping velocity is comparable with the Alfven speed of the magnetic layer, flux pumping is able to hold back the bulk of the magnetic field, with only small pockets of strong field able to rise into the upper layer. In simulations in which the magnetic layer is generated by shear, we find that the shear velocity is not necessarily required to exceed that of the pumping (therefore the kinetic energy of the shear is not required to exceed that of the overlying convection), for strong localised pockets of magnetic field to be produced which can rise into the upper layer. This is because magnetic flux pumping acts to store the field below the interface, allowing it to be amplified both by the shear, and by vortical fluid motions, until pockets of field can achieve sufficient strength to rise into the upper layer. In addition, we find that the interface between the two layers is a natural location for the production of strong vertical gradients in the magnetic field. If these gradients are sufficiently strong to allow the development of magnetic buoyancy instabilities, strong shear is not necessarily required to drive them (c.f. previous work by Vasil & Brummell). We find that the addition of magnetic flux pumping appears to be able to assist shear-driven magnetic buoyancy in producing strong flux concentrations that can rise into the convection zone from the radiative interior.
We present the first unambiguous detection of quasi-periodic wave trains within the broad pulse of a global EUV wave (so-called "EIT wave") occurring on the limb. These wave trains, running ahead of the lateral CME front of 2-4 times slower, coherently travel to distances $>R_{sun}/2$ along the solar surface, with initial velocities up to 1400 km/s decelerating to ~650 km/s. The rapid expansion of the CME initiated at an elevated height of 110 Mm produces a strong downward and lateral compression, which may play an important role in driving the primary EUV wave and shaping its front forwardly inclined toward the solar surface. The waves have a dominant 2 min periodicity that matches the X-ray flare pulsations, suggesting a causal connection. The arrival of the leading EUV wave front at increasing distances produces an uninterrupted chain sequence of deflections and/or transverse (likely fast kink mode) oscillations of local structures, including a flux-rope coronal cavity and its embedded filament with delayed onsets consistent with the wave travel time at an elevated (by ~50%) velocity within it. This suggests that the EUV wave penetrates through a topological separatrix surface into the cavity, unexpected from CME caused magnetic reconfiguration. These observations, when taken together, provide compelling evidence of the fast-mode MHD wave nature of the {\it primary (outer) fast component} of a global EUV wave, running ahead of the {\it secondary (inner) slow} component of CME-caused restructuring.
Weakly Interactive Massive Particles are among the most motivated candidates invoked to solve the Dark Matter puzzle. The stability of Dark Matter (DM) typically results from a global protective symmetry. However, the presence of gravitational effects may cause its explicit breaking. We show that this provides a new source of signal for indirect DM searches as well as the embedding of a large class of decaying Dark Matter models into the WIMP paradigm.
We study the first Kaluza-Klein excitation of the Higgs boson in universal extra dimensions as a dark matter candidate. The first-level Higgs boson could be the lightest Kaluza-Klein particle, which is stable due to the conservation of Kaluza-Klein parity, in non-minimal models where boundary localized terms modify the mass spectrum. We calculate the relic abundance and find that it agrees with the observed dark matter density if the mass of the first-level Higgs boson is about 2 TeV. We study also the prospects for detection of this dark matter candidate in direct as well as indirect detection experiments. Although the first-level Higgs boson is a typical weakly interacting massive particle, an observation in any of the conventional experiments is very challenging.
Planck scale physics may influence the evolution of cosmological fluctuations in the early stages of cosmological evolution. Because of the quasi-exponential redshifting, which occurs during an inflationary period, the physical wavelengths of comoving scales that correspond to the present large-scale structure of the Universe were smaller than the Planck length in the early stages of the inflationary period. This trans-Planckian effect was studied before using toy models. The Horava-Lifshitz (HL) theory offers the chance to study this problem in a candidate UV complete theory of gravity. In this paper we study the evolution of cosmological perturbations according to HL gravity assuming that matter gives rise to an inflationary background. As is usually done in inflationary cosmology, we assume that the fluctuations originate in their minimum energy state. In the trans-Planckian region the fluctuations obey a non-linear dispersion relation of Corley-Jacobson type. In the "healthy extension" of HL gravity there is an extra degree of freedom which plays an important role in the UV region but decouples in the IR, and which influences the cosmological perturbations. We find that in spite of these important changes compared to the usual description, the overall scale-invariance of the power spectrum of cosmological perturbations is recovered. However, we obtain oscillations in the spectrum as a function of wavenumber with a relative amplitude of order unity and with an effective frequency which scales nonlinearly with wavenumber. Taking the usual inflationary parameters we find that the frequency of the oscillations is so large as to render the effect difficult to observe.
We propose the cosmological model which allows to describe on equal footing the evolution of matter in the universe on the time interval from the inflation till the domination of dark energy. The matter is considered as a two-component perfect fluid imitated by homogeneous scalar fields between which there is energy exchange. Dark energy is represented by the cosmological constant, which is supposed invariable during the whole evolution of the universe. The matter changes its equation of state with time, so that the era of radiation domination in the early universe smoothly passes into the era of a pressureless gas, which then passes into the late-time epoch, when the matter is represented by a gas of low-velocity cosmic strings. The inflationary phase is described as an analytic continuation of the energy density in the very early universe into the region of small negative values of the parameter which characterizes typical time of energy transfer from one matter component to another. The Hubble expansion rate, energy density of the matter, energy density parameter, and deceleration parameter as functions of time are found.
This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.
In this note I respond to Vilenkin's claim that there must have been a beginning.
While the finite-temperature effective potential in a gauge theory is a gauge-dependent quantity, in several instances a first-order phase transition can be triggered by gauge-independent terms. A particularly interesting case occurs when the potential barrier separating the broken and symmetric vacua of a spontaneously broken symmetry is produced by tree-level terms in the potential. Here, we study this scenario in a simple Abelian Higgs model, for which the gauge-invariant potential is known, augmented with a singlet real scalar. We analyze the possible symmetry breaking patterns in the model, and illustrate in which cases gauge artifacts are expected to manifest themselves most severely. We then show that gauge artifacts can be pronounced even in the presence of a relatively large, tree-level singlet-Higgs cubic interaction. When the transition is strongly first order, these artifacts, while present, are more subtle than in the generic situation.
Links to: arXiv, form interface, find, astro-ph, recent, 1204, contact, help (Access key information)
We perform the largest currently available set of direct N-body calculations of young star cluster models to study the dynamical influence, especially through the ejections of the most massive star in the cluster, on the current relation between the maximum-stellar-mass and the star-cluster-mass. We vary several initial parameters such as the initial half-mass radius of the cluster, the initial binary fraction, and the degree of initial mass segregation. Two different pairing methods are used to construct massive binaries for more realistic initial conditions of massive binaries. We find that lower mass clusters (<= 10^2.5 Msun) do not shoot out their heaviest star. In the case of massive clusters (>= 1000 Msun), no most-massive star escapes the cluster within 3 Myr regardless of the initial conditions if clusters have initial half-mass radii, r_0.5, >= 0.8 pc. However, a few of the initially smaller sized clusters (r_0.5 = 0.3 pc), which have a higher density, eject their most massive star within 3 Myr. If clusters form with a compact size and their massive stars are born in a binary system with a mass-ratio biased towards unity, the probability that the mass of the most massive star in the cluster changes due to the ejection of the initially most massive star can be as large as 20 per cent. Stellar collisions increase the maximum-stellar-mass in a large number of clusters when clusters are relatively dense (M_ecl >= 10^3 Msun and r_0.5 = 0.3 pc) and binary-rich. Overall, we conclude that dynamical effects hardly influence the observational maximum-stellar-mass -- cluster mass relation.
We present broad-band observations of the afterglow and environment of the short GRB 111020A. An extensive X-ray light curve from Swift/XRT, XMM-Newton and Chandra, spanning ~100 seconds to 10 days after the burst, reveals a significant break at t~2 days with pre- and post-break decline rates of alphaX,1 ~ -0.78 and alphaX,2<-1.7, respectively. Interpreted as a jet break, we infer a collimated outflow with an opening angle of ~3-8 degrees. The resulting beaming-corrected gamma-ray (10-1000 keV band) and blastwave kinetic energies are (2-3)e48 erg and (0.3-2)e49 erg, respectively, with the range depending on the unknown redshift of the burst. We report a radio afterglow limit of <39 microJy (3-sigma) from EVLA observations which, along with our finding that vc<vX, constrains the circumburst density to n~0.01-0.1 cm^(-3). Optical observations provide an afterglow limit of i>24.4 mag at 18 hours after the burst, and reveal a potential host galaxy with i~24.3 mag. The sub-arcsecond localization from Chandra provides a precise offset of 0.80"+/-0.11" (1-sigma) from this galaxy corresponding to an offset of 5-7 kpc for z=0.5-1.5. We find a high excess neutral Hydrogen column density of (7.5+/-2.0)e21 cm^(-2) (z=0). Our observations demonstrate that a growing fraction of short GRBs are collimated which may lead to a true event rate of >100-1000 Gpc^(-3) yr^(-1), in good agreement with the NS-NS merger rate of ~200-3000 Gpc^(-3) yr^(-1). This consistency is promising for coincident short GRB-gravitational wave searches in the forthcoming era of Advanced LIGO/VIRGO.
The observed distribution of neutron star (NS) masses reflects the physics of core-collapse supernova explosions and the structure of the massive stars that produce them at the end of their evolution. We present a Bayesian analysis that directly compares the NS mass distribution observed in double NS systems to theoretical models of NS formation. We find that models with standard binary mass ratio distributions are strongly preferred over independently picking the masses from the initial mass function, although the strength of the inference depends on whether current assumptions for identifying the remnants of the primary and secondary stars are correct. Second, NS formation models with no mass fallback are favored because they reduce the dispersion in NS masses. The double NS system masses thus directly point to the mass coordinate where the supernova explosion was initiated, making them an excellent probe of the supernova explosion mechanism. If we assume no fallback and simply vary the mass coordinate separating the remnant and the supernova ejecta, we find that for solar metallicity stars the explosion most likely develops at the edge of the iron core at a specific entropy of about 2.8 k_B. The primary limitations of our study are the poor knowledge of the supernova explosion mechanism and the lack of broad range of SN model explosions of LMC to solar metallicity.
We present new detections of the CO(5-4), CO(7-6), [CI](1-0) and [CI](2-1) molecular and atomic line transitions towards the unlensed, obscured quasar AMS12 (z=2.7672), observed with the IRAM PdBI. This is the first unlensed, high redshift source to have both [CI] transitions detected. Continuum measurements between 70 $\mu$m and 3 mm are used to constrain the FIR SED, and we find a best fit FIR luminosity of log[Lfir/Lsol] = 13.5+/-0.1, dust temperature T_d = 88+/-8 K and emissivity index {\beta} = 0.6+/-0.1. The highly-excited molecular gas probed by CO(3-2), (5-4) and (7-6), is modelled with large velocity gradient (LVG) models. The gas kinetic temperature T_g, density n(H2), and the characteristic size r0, are determined using the dust temperature from the FIR SED as a prior for the gas temperature. The best fitting parameters are T_g = 90+/-8 K, n(H2) = 10^(3.9+/-0.1) cm^(-3) and r0 = 0.8+/-0.04 kpc. The ratio of the [CI] lines gives a [CI] excitation temperature of 43+/-10 K, indicating the [CI] and the high-excitation CO are not in thermal equilibrium. The [CI] excitation temperature is below that of T_d and T_g of the high-excitation CO, perhaps because [CI] lies at a larger radius where there may also be a large reservoir of CO at a cooler temperature, perhaps detectable through the CO(1-0). Using the [CI](1-0) line we can estimate the strength of the CO(1-0) line and hence the gas mass. This suggests that a significant fraction (~30%) of the molecular gas is missed from the high-excitation line analysis. The Eddington limited black hole mass is found from the bolometric luminosity to be Mbh >~ 1.5x10^9 Msol. Along with the stellar mass of 3x10^11 Msol, these give a black hole - bulge mass ratio of Mbh/Mbulge >~ 0.005. This is in agreement with studies on the evolution of the Mbh/Mbulge relationship at high redshifts, which find a departure from the local value ~0.002.
We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the ~5 Myr old starburst cluster Westerlund 1 and ~2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.
In this mini-review we discuss first why we should investigate cosmological models beyond LCDM. We then show how to describe dark energy or modified gravity models in a fluid language with the help of one background and two perturbation quantities. We review a range of dark energy models and study how they fit into the phenomenological framework, including generalizations like phantom crossing, sound speeds different from c and non-zero anisotropic stress, and how these effective quantities are linked to the underlying physical models. We also discuss the limits of what can be measured with cosmological data, and some challenges for the framework.
We test two different methods of using near-infrared extinction to estimate distances to dark clouds in the first quadrant of the Galaxy using large near infrared (2MASS and UKIDSS) surveys. VLBI parallax measurements of masers around massive young stars provide the most direct and bias-free measurement of the distance to these dark clouds. We compare the extinction distance estimates to these maser parallax distances. We also compare these distances to kinematic distances, including recent re-calibrations of the Galactic rotation curve. The extinction distance methods agree with the maser parallax distances (within the errors) between 66% and 100% of the time (depending on method and input survey) and between 85% and 100% of the time outside of the crowded Galactic center. Although the sample size is small, extinction distance methods reproduce maser parallax distances better than kinematic distances; furthermore, extinction distance methods do not suffer from the kinematic distance ambiguity. This validation gives us confidence that these extinction methods may be extended to additional dark clouds where maser parallaxes are not available.
We discuss the relevance of directional detection experiments in the post-discovery era and propose a method to extract the local dark matter phase space distribution from directional data. The first feature of this method is a parameterization of the dark matter distribution function in terms of integrals of motion, which can be analytically extended to infer properties of the global distribution if certain equilibrium conditions hold. The second feature of our method is a decomposition of the distribution function in moments of a model independent basis, with minimal reliance on the ansatz for its functional form. We illustrate our method using the Via Lactea II N-body simulation as well as an analytical model for the dark matter halo. We conclude that O(1000) events are necessary to measure deviations from the Standard Halo Model and constrain or measure the presence of anisotropies.
We review the different frameworks in which Galileon scalar fields have been seen to emerge such an in DGP, New Massive Gravity and Ghost-free massive Gravity and emphasize their relation with the Lovelock invariant in braneworld models. The existence of a non-renormalization theorem for Galileon scalar fields makes them especially attractive candidates for inflation as well as for late-time acceleration. In particular we review the self-accelerating and degravitating branches of solutions present in Galileon models when arising from Massive Gravity and discuss their phenomenology.
The discovery of cosmic acceleration is one of the most important developments in modern cosmology. The observation, thirteen years ago, that type Ia supernovae appear dimmer that they would have been in a decelerating universe followed by a series of independent observations involving galaxies and cluster of galaxies as well as the cosmic microwave background, all point in the same direction: we seem to be living in a flat universe whose expansion is currently undergoing an acceleration phase. In this paper, we review the various observational evidences, most of them gathered in the last decade, and the improvements expected from projects currently collecting data or in preparation.
We present the 1.9--4.2um spectra of the five bright (L<11.2) young stars associated with silhouette disks with moderate to high inclination angle of 39--80deg in the M42 and M43 regions. The water ice absorption is seen toward d121-1925 and d216-0939, while the spectra of d182-316, d183-405, and d218-354 show no water ice feature around 3.1um within the detection limits. By comparing the water ice features toward nearby stars, we find that the water ice absorption toward d121-1925 and d216-0939 most likely originates from the foreground material and the surrounding disk, respectively. The angle of the disk inclination is found to be mainly responsible for the difference of the optical depth of the water ice among the five young stars. Our results suggest that there is a critical inclination angle between 65deg and 75deg for the circumstellar disk where the water ice absorption becomes strong. The average density at the disk surface of d216-0939 was found to be 6.38x10^(-18) g cm^(-3). The water ice absorption band in the d216-0939 disk is remarkable in that the maximum optical depth of the water ice band is at a longer wavelength than detected before. It indicates that the primary carrier of the feature is purely crystallized water ice at the surface of the d216-0939 disk with characteristic size of ~0.8um, which suggests grain growth. This is the first direct detection of purely crystallized water ice in a silhouette disk.
Assuming the universe is spatially homogeneous on the largest scales lays the foundation for almost all cosmology. This idea is based on the Copernican principle, that we are not at a particularly special place in the universe. Surprisingly, this philosophical assumption has yet to be rigorously demonstrated independently of the standard paradigm. This issue has been brought to light by cosmological models which can potentially explain apparent acceleration by spatial inhomogeneity rather than dark energy. These models replace the temporal fine tuning associated with Lambda with a spatial fine tuning, and so violate the Copernican assumption. While is seems unlikely that such models can really give a realistic solution to the dark energy problem, they do reveal how poorly constrained radial inhomogeneity actually is. So the bigger issue remains: How do we robustly test the Copernican principle independently of dark energy or theory of gravity?
We present for the first time metallicity maps generated using data from the Wide Field Spectrograph (WiFeS) on the ANU 2.3m of 9 Luminous Infrared Galaxies (LIRGs) and discuss the abundance gradients and distribution of metals in these systems. We have carried out optical integral field spectroscopy (IFS) of several several LIRGs in various merger phases to investigate the merger process. In a major merger of two spiral galaxies with preexisting disk abundance gradients, the changing distribution of metals can be used as a tracer of gas flows in the merging system as low metallicity gas is transported from the outskirts of each galaxy to their nuclei. We employ this fact to probe merger properties by using the emission lines in our IFS data to calculate the gas-phase metallicity in each system. We create abundance maps and subsequently derive a metallicity gradient from each map. We compare our measured gradients to merger stage as well as several possible tracers of merger progress and observed nuclear abundances. We discuss our work in the context of previous abundance gradient observations and compare our results to new galaxy merger models which trace metallicity gradient. Our results agree with the observed flattening of metallicity gradients as a merger progresses. We compare our results with new theoretical predictions that include chemical enrichment. Our data show remarkable agreement with these simulations.
We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery from errors. The framework can be efficiently implemented and incurs an acceptable overhead. In practice, we measure the time spent in the framework to be less than 1% of the wall-clock time. Due to the unified structure of the interface, incorporating multiple modules addressing the same physics in different ways is relatively straightforward. Different modules may be advanced serially or in parallel. Despite initial concerns, we have encountered relatively few problems with this strict separation between modules, and the results of our simulations are consistent with earlier results using more traditional monolithic approaches. This framework provides a platform to combine existing simulation codes or develop new physical solver codes within a rich "ecosystem" of interchangeable modules.
In this paper, we investigate the structure of condensation fronts from warm diffuse gas to cold neutral medium (CNM) under the plane parallel geometry. The solutions have two parameters, the pressure of the cold neutral medium (CNM) and the mass flux across the transition front and their ranges are much wider than previous thought. First, we consider the pressure range where the three phases, the CNM, the unstable phase, and the warm neutral medium, can coexist in the pressure equilibrium. In a wide range of the mass flux, we find solutions connecting the CNM and unstable phase. Moreover, we find solutions in larger pressure range where there is only one thermal equilibrium state, or the CNM. These solutions can be realized in shock-compressed regions that is promising sites of the molecular cloud formation. We also find remarkable properties in our solutions. The heat conduction becomes less important with increasing mass flux, and the thickness of the transition layer is characterized by the cooling length instead of the Field length.
We present the statistical properties of molecular clumps in the Galactic center 50 km s$^{-1}$ molecular cloud (GCM-0.02-0.07) based on observations of the CS $J=1-0$ emission line with the Nobeyama Millimeter Array. In the cloud, 37 molecular clumps with local thermal equilibrium (LTE) masses of $2\times10^2-6\times10^3 M_\odot$ were identified by using the {\it clumpfind} algorithm. The velocity widths of the molecular clumps are about five-fold those of Galactic disk molecular clouds with the same radius. The virial-theorem masses are three-fold the LTE masses. The mass and size spectra can be described by power laws of $dN/dM\propto M^{-2.6\pm0.1}$ ($M\gtrsim 900M_\odot$) and $dN/dR\propto R^{-5.9\pm0.3}$ ($R\gtrsim 0.35$ pc), respectively. The statistical properties of the region interacting with the Sgr A East shell and those of the non-interacting part of the cloud are significantly different. The interaction probably makes the mass function steeper, from $dN/dM\propto M^{-2.0\pm0.1}$ in the non-interacting part to $dN/dM\propto M^{-4.0\pm0.2}$ in the interacting region. On the other hand, the interaction presumably truncates the size spectrum on the larger side of $R\sim 0.4$ pc.
Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.
We compile from literature the broadband SEDs of twelve TeV blazars observed simultaneously or quasi-simultaneously with Fermi/LAT and other instruments. Two SEDs are available for each of the objects and the state is identified as a low or high state according to its flux density at GeV/TeV band. The observed SEDs of BL Lac objects (BL Lacs) are fitted well with the synchrotron + synchrotron-self-Compton (syn+SSC) model, whereas the SEDs of the two flat spectrum radio quasars (FSRQs) need to include the contributions of external Compton scattering. In this scenario, it is found that the Doppler factor delta of FSRQs is smaller than that of BL Lacs, but the magnetic field strength B of FSRQs is larger than that of BL Lacs. The increase of the peak frequency of the SEDs is accompanied with the increase of the flux for the individual sources, which seems opposite to the observational phenomena of the blazar sequence. We refer this phenomenon to blazar anti-sequence of spectral variability for individual TeV blazars. However, both the blazar sequence from FSRQs to BL Lacs and blazar anti-sequence of the spectral variability from low state to high state are accompanied by an increase of the break Lorentz factor of the electron's spectrum gamma_b and a decrease of B. We propose a model in which the mass accretion rate is the driving force behind both the blazar sequence for ensembles of blazars and the blazar anti-sequence for individual blazars. Specifically we suggest that the differences in mass accretion rate of different blazars produce the observed blazar sequence, but \Delta M in each blazar results in the observed blazar anti-sequence.
It has been suggested that type II radio bursts are due to energetic electrons accelerated at coronal shocks. Radio observations, however, have poor or no spatial resolutions to pinpoint the exact acceleration locations of these electrons. In this paper, we discuss a promising approach to infer the electron acceleration location by combining radio and white light observations. The key assumption is to relate specific morphological features (e.g. spectral bumps) of the dynamic spectra of type II radio bursts, to imaging features (e.g. CME going into a streamer) along the CME (and its driven shock) propagation. In this study, we examine the CME-streamer interaction for the solar eruption dated on 2003 November 1. The presence of spectral bump in the relevant type II radio burst is identified, which is interpreted as a natural result of the shock-radio emitting region entering the dense streamer structure. The study is useful for further determinations of the location of type II radio burst and the associated electron acceleration by CME-driven shock.
We present the first numerical simulations that self-consistently follow the formation of dense molecular clouds in colliding flows. Our calculations include a time-dependent model for the H2 and CO chemistry that runs alongside a detailed treatment of the dominant heating and cooling processes in the ISM. We adopt initial conditions characteristic of the warm neutral medium and study two different flow velocities - a slow flow with v = 6.8 km/s and a fast flow with v = 13.6 km/s. The clouds formed by the collision of these flows form stars, with star formation beginning after 16 Myr in the case of the slower flow, but after only 4.4 Myr in the case of the faster flow. In both flows, the formation of CO-dominated regions occurs only around 2 Myr before the onset of star formation. Prior to this, the clouds produce very little emission in the J = 1 -> 0 transition line of CO, and would probably not be identified as molecular clouds in observational surveys. In contrast, our models show that H2-dominated regions can form much earlier, with the timing depending on the details of the flow. In the case of the slow flow, small pockets of gas become fully molecular around 10 Myr before star formation begins, while in the fast flow, the first H2-dominated regions occur around 3 Myr before the first prestellar cores form. Our results are consistent with models of molecular cloud formation in which the clouds are dominated by "dark" molecular gas for a considerable proportion of their assembly history.
A modified Chaplygin gas model (MCG), $\rho_{MCG}/\rho_{MCG0}=[B_{s}+(1-B_{s})a^{-3(1+B)(1+\alpha)}]^{1/(1+\alpha)}$, as a unified dark matter model and dark energy model is constrained by using current available cosmic observational data points which include type Ia supernovae, baryon acoustic oscillation and the seventh year full WMAP data points. As a contrast to the consideration in the literatures, we {\it do not} separate the MCG into two components, i.e. dark mater and dark energy component, but we take it as a whole energy component-a unified dark sector. By using Markov Chain Monte Carlo method, a tight constraint is obtained: $\alpha= 0.000727_{- 0.00140- 0.00234}^{+ 0.00142+ 0.00391}$, $B=0.000777_{- 0.000302- 0.000697}^{+ 0.000201+ 0.000915}$ and $B_s= 0.782_{- 0.0162- 0.0329}^{+ 0.0163+ 0.0307}$ .}
We present an analytic model for the local bias of dark matter halos in a LCDM universe. The model uses the halo mass density instead of the halo number density and is searched for various halo mass cuts, smoothing lengths, and redshift epoches. We find that, when the logarithmic density is used, the second-order polynomial can fit the numerical relation between the halo mass distribution and the underlying matter distribution extremely well. In this model the logarithm of the dark matter density is expanded in terms of log halo mass density to the second order. The model remains excellent for all halo mass cuts (from M_{cut}=3\times10^{11}$ to $3\times10^{12}h^{-1}M_{\odot}$), smoothing scales (from $R=5h^{-1}$Mpc to $50h^{-1}$Mpc), and redshift ranges (from z=0 to 1.0) considered in this study. The stochastic term in the relation is found not entirely random, but a part of the term can be determined by the magnitude of the shear tensor.
Using Hinode EUV Imaging Spectrometer (EIS) spectra recorded daily at Sun center from the end of 2006 to early 2011, we studied the long-term evolution of the quiet corona. The light curves of the higher temperature emission lines exhibit larger variations in sync with the solar activity cycle while the cooler lines show reduced modulation. Our study shows that the high temperature component of the corona changes in quiet regions, even though the coronal electron density remains almost constant there. The results suggest that heat input to the quiet corona varies with the solar activity cycle.
Main aim of this paper is the first light curve and apsidal motion analysis of thirteen eccentric eclipsing binaries and to reveal their basic physical properties. All of the systems were studied by method of period analysis of times of minima and the light curve analysis. Many new times of minima for all of the systems were derived and collected from the data obtained by the automatic, robotic or satellite telescopes. This allow us to study the apsidal motion in these systems in detail for the first time. From the light curve analysis the first rough estimations of the physical properties of these systems were done. The analyzed systems undergo an apsidal motion with the following periods in years: AR CMa (44 +/- 10), OZ Hya (117 +/- 53), V498 Mon (62 +/- 4), V521 Mon (217 +/- 37), V684 Mon (74.5 +/- 20), V730 Mon (39 +/- 12), GV Nor (197 +/- 67), NS Nor (516 +/- 230), TZ Pyx (157 +/- 37), V385 Sco (1926 +/- 980), V629 Sco (56 +/- 17), V881 Sco (131 +/- 48), V1082 Sco (186 +/- 280). Also 190 new minima times were derived and included. The period and light curve analyses were done, however a more detailed spectroscopic analysis is needed to confirm the physical parameters of the components with higher conclusiveness.
B-type stars are promising targets for asteroseismic modelling, since their
frequency spectrum is relatively simple.
We deduce and summarise observational constraints for the hybrid pulsator,
HD50230, earlier reported to have deviations from a uniform period spacing of
its gravity modes. The combination of spectra and a high-quality light curve
measured by the CoRoT satellite allow a combined approach to fix the position
of HD50230 in the HR diagram.
To describe the observed pulsations, classical Fourier analysis was combined
with short-time Fourier transformations and frequency spacing analysis
techniques. Visual spectra were used to constrain the projected rotation rate
of the star and the fundamental parameters of the target. In a first
approximation, the combined information was used to interpret multiplets and
spacings to infer the true surface rotation rate and a rough estimate of the
inclination angle.
We identify HD50230 as a spectroscopic binary and characterise the two
components. We detect the simultaneous presence of high-order g modes and
low-order p and g-modes in the CoRoT light curve, but were unable to link them
to line profile variations in the spectroscopic time series. We extract the
relevant information from the frequency spectrum, which can be used for seismic
modelling, and explore possible interpretations of the pressure mode spectrum.
Deep, near-infrared JHK-maps were observed for 10 nearby, grand-design,
spiral galaxies using HAWK-I/VLT to study the distribution of young stellar
clusters in them and thereby determine whether strong spiral perturbations can
influence star formation. Complete, magnitude-limited candidate lists of
star-forming complexes were obtained by searching within the K-band maps. The
properties of the complexes were derived from (H-K)-(J-H) diagrams including
the identification of the youngest complexes (i.e. <7 Myr) and the estimation
of their extinction.
Young stellar clusters with ages <7 Myr have significant internal extinction
in the range of Av=3-7m, while older ones typically have Av<1m. The cluster
luminosity function (CLF) is well-fitted by a power law with an exponent of
around -2 and displays no evidence of a high luminosity cut-off. The brightest
cluster complexes in the disk reach luminosities of Mk = -15.5m or estimated
masses of 10^6 Mo. At radii with a strong, two-armed spiral pattern, the star
formation rate in the arms is higher by a factor of 2-5 than in the inter-arm
regions. The CLF in the arms is also shifted towards brighter Mk by at least
0.4m. We also detect clusters with colors compatible with Large Magellanic
Cloud intermediate age clusters and Milky Way globular clusters. The (J-K)-Mk
diagram of several galaxies shows, for the brightest clusters, a clear
separation between young clusters that are highly attenuated by dust and older
ones with low extinction.
The gap in the (J-K)-Mk diagrams implies that there has been a rapid
expulsion of dust at an age around 7 Myr, possibly triggered by supernovae.
Strong spiral perturbations concentrate the formation of clusters in the arm
regions and shifts their CLF towards brighter magnitudes.
We have previously reported on chemical abundance trends with evolutionary
state in the globular cluster NGC 6397 discovered in analyses of spectra taken
with FLAMES at the VLT. Here, we reinvestigate the FLAMES-UVES sample of 18
stars, ranging from just above the turnoff point (TOP) to the red giant branch
below the bump. Inspired by new calibrations of the infrared flux method, we
adopt a set of hotter temperature scales. Chemical abundances are determined
for six elements (Li, Mg, Ca, Ti, Cr, and Fe). Signatures of cluster-internal
pollution are identified and corrected for in the analysis of Mg.
On the modified temperature scales, evolutionary trends in the abundances of
Mg and Fe are found to be significant at the 2{\sigma} and 3{\sigma} levels,
respectively. The detailed evolution of abundances for all six elements agrees
with theoretical isochrones, calculated with effects of atomic diffusion and a
weak to moderately strong efficiency of turbulent mixing. The age of these
models is compatible with the external determination from the white dwarf
cooling sequence. We find that the abundance analysis cannot be reconciled with
the strong turbulent-mixing efficiency inferred elsewhere for halo field stars.
A weak mixing efficiency reproduces observations best, indicating a
diffusion-corrected primordial lithium abundance of log {\epsilon}(Li) = 2.57
+- 0.10. At 1.2{\sigma}, this value agrees well with WMAP-calibrated Big-Bang
nucleosynthesis predictions.
We report on sensitive observations of the CO(7-6) and CI(2-1) transitions in the z=2.79 QSO host galaxy RXJ0911.4+0551 using the IRAM Plateau de Bure interferometer (PdBI). Our extremely high signal to noise spectra combined with the narrow CO line width of this source (FWHM = 120 km/s) allows us to estimate sensitive limits on the space-time variations of the fundamental constants using two emission lines. Our observations show that the CI and CO line shapes are in good agreement with each other but that the CI line profile is of order 10% narrower, presumably due to the lower opacity in the latter line. Both lines show faint wings with velocities up to +/-250 km/s, indicative of a molecular outflow. As such the data provide direct evidence for negative feedback in the molecular gas phase at high redshift. Our observations allow us to determine the observed frequencies of both transitions with so far unmatched accuracy at high redshift. The redshift difference between the CO and CI lines is sensitive to variations of dF/F with F=alpha^2/mu where alpha is the fine structure constant and mu the proton-to-electron mass ratio. We find dF/F=6.9 +/-3.7 x 10^-6 at a lookback time of 11.3 Gyr, which within the uncertainties, is consistent with no variations of the fundamental constants.
Observations of neutral hydrogen in spiral galaxies reveal a sharp cutoff in the radial density profile at some distance from the center. Using 22 galaxies with known HI distributions as an example, we discuss the question of whether this effect can be associated exclusively with external ionizing radiation, as is commonly assumed. We show that before the surface density reaches $\sigma_{\textrm{HI}}\le 0.5 {\cal M}_\odot/{\textrm {pc}}^2$(the same for galaxies of different types), it is hard to expect the gas to be fully ionized by background radiation. For two of 13 galaxies with a sharp drop in the HI profile, the "steepening" can actually be caused by ionization. At the same time, for the remaining galaxies, the observed cutoff in the radial HI profile is closer to the center than if it was a consequence of ionization by background radiation and, therefore, it should be caused by other factors.
We survey the job situation of women in astronomy in Germany and of German women abroad and review indicators for their career development. Our sample includes women astronomers from all academic levels from doctoral students to professors, as well as female astronomers who have left the field. We find that networking and human support are among the most important factors for success. Experience shows that students should carefully choose their supervisor and collect practical knowledge abroad. We reflect the private situation of female German astronomers and find that prejudices are abundant, and are perceived as discriminating.We identify reasons why women are more likely than men to quit astronomy after they obtain their PhD degree. We give recommendations to young students on what to pay attention to in order to be on the successful path in astronomy.
Recent observational results indicate that the functional shape of the spatially-resolved star formation-molecular gas density relation depends on the spatial scale considered. These results may indicate a fundamental role of sampling effects on scales that are typically only a few times larger than those of the largest molecular clouds. To investigate the impact of this effect, we construct simple models for the distribution of molecular clouds in a typical star-forming spiral galaxy, and, assuming a power-law relation between SFR and cloud mass, explore a range of input parameters. We confirm that the slope and the scatter of the simulated SFR-molecular gas surface density relation depend on the size of the sub-galactic region considered, due to stochastic sampling of the molecular cloud mass function, and the effect is larger for steeper relations between SFR and molecular gas. There is a general trend for all slope values to tend to ~unity for region sizes larger than 1-2 kpc, irrespective of the input SFR-cloud relation. The region size of 1-2 kpc corresponds to the area where the cloud mass function becomes fully sampled. We quantify the effects of selection biases in data tracing the SFR, either as thresholds (i.e., clouds smaller than a given mass value do not form stars) or backgrounds (e.g., diffuse emission unrelated to current star formation is counted towards the SFR). Apparently discordant observational results are brought into agreement via this simple model, and the comparison of our simulations with data for a few galaxies supports a steep (>1) power law index between SFR and molecular gas.
We have carried out survey observations toward the Galactic plane at l~38 deg in the 12CO and 13CO J=1-0 lines using the Nobeyama Radio Observatory 45-m telescope. A wide area (0.8 x 0.8 deg) was mapped with high spatial resolution (17"). The line of sight samples the gas in both the Sagittarius arm and the inter-arm regions. The present observations reveal how the structure and physical conditions vary across a spiral arm. We classify the molecular gas in the line of sight into two distinct components based on its appearance: the bright and compact B component and the fainter and diffuse (i.e., more extended) D component. The B component is predominantly seen at the spiral arm velocities, while the D component dominates at the inter-arm velocities and is also found at the spiral arm velocities. We introduce the brightness distribution function and the brightness distribution index (BDI, which indicates the dominance of the B component) in order to quantify the map's appearance. The radial velocities of BDI peaks coincide with those of high 12CO J=3-2/12CO J=1-0 intensity ratio (i.e., warm gas) and H II regions, and tend to be offset from the line brightness peaks at lower velocities (i.e., presumably downstream side of the arm). Our observations reveal that the gas structure at small scales changes across a spiral arm: bright and spatially confined structures develop in a spiral arm, leading to star formation at downstream side, while extended emission dominates in the inter-arm region.
The massive binary system Eta Carinae and the surrounding HII complex, the Carina Nebula, are potential particle acceleration sites from which very-high-energy (VHE; E > 100 GeV) \gamma-ray emission could be expected. This paper presents data collected during VHE \gamma-ray observations with the H.E.S.S. telescope array from 2004 to 2010, which cover a full orbit of Eta Carinae. In the 33.1-hour data set no hint of significant \gamma-ray emission from Eta Carinae has been found and an upper limit on the \gamma-ray flux of 7.7 x 10-13 ph cm-2 s-1 (99% confidence level) is derived above the energy threshold of 470 GeV. Together with the detection of high-energy (HE; 0.1 GeV > E > 100 GeV) \gamma-ray emission by the Fermi-LAT up to 100 GeV, and assuming a continuation of the average HE spectral index into the VHE domain, these results imply a cut-off in the \gamma-ray spectrum between the HE and VHE \gamma-ray range. This could be caused either by a cut-off in the accelerated particle distribution or by severe \gamma-\gamma\ absorption losses in the wind collision region. Furthermore, the search for extended \gamma-ray emission from the Carina Nebula resulted in an upper limit on the \gamma-ray flux of 4.2 x 10-12 ph cm-2 s-1 (99% confidence level). The derived upper limit of ~23 on the cosmic-ray enhancement factor is compared with results found for the old-age mixed-morphology supernova remnant W 28.
We investigate the possibilities of wavelet-based RM Synthesis for the recognition of structures of regular and turbulent magnetic fields in extended magnetized objects, like galaxies and galaxy clusters. Wavelets allow to reformulate the RM Synthesis method in a scale-dependent way and to visualize the data as a function of Faraday depth and scale. We present observational tests to recognize regular magnetic fields without and with one or two reversals along the line of sight and imprints of turbulent magnetic fields. A region with a regular magnetic field generates a broad "disk" in Faraday space ("Faraday spectrum"), with two "horns" if the distribution of cosmic-ray electrons is broader than that of the thermal electrons. Each field reversal generates one asymmetric "horn" on top of the "disk". A region with a turbulent field can be recognized by a "Faraday forest" of many components. We argue that the ratio of maximum to minimum wavelengths is an important parameter because it determines the range of scales which can be identified in Faraday space. Full recognition of magnetic field structures in spiral galaxies or galaxy clusters requires analysis of data cubes in position-position-Faraday depth ("PPF cubes"), observed over a wide and continuous wavelength range, from about 100 MHz to several GHz, providing good resolution as well as recognition of a wide range of scales in Faraday space. The planned SKA fulfils this condition and will be close to a perfect "Faraday telescope". The combination of data from the present telescopes LOFAR (low frequencies) and EVLA (high frequencies) looks promising for the recognition of magnetic structures at all scales. The addition of WSRT or GMRT data at intermediate frequencies to those from LOFAR and EVLA data fills the gap between the LOFAR and EVLA wavelength ranges.
The ionizing star of the planetary nebula NGC 2392 is too cool to explain the high excitation of the nebular shell, and an additional ionizing source is necessary. We use photoionization modeling to estimate the temperature and luminosity of the putative companion. Our results show it is likely to be a very hot (Teff ~ 250kK), dense white dwarf. If the stars form a close binary, they may merge within a Hubble time, possibly producing a Type Ia supernova.
The effect of rarefaction acceleration on the propagation dynamics and structure of relativistically hot jets is studied through relativistic hydrodynamic simulations. We emphasize the nonlinear interaction of rarefaction waves excited at the interface between a cylindrical jet and the surrounding medium. From simplified one-dimensional models with radial jet structure, we find that a decrease in the relativistic pressure due to the interacting rarefaction waves in the central zone of the jet transiently yields a more powerful boost of the bulk jet than that expected from single rarefaction acceleration. This leads to a cyclic in-situ energy conversion between thermal and bulk kinetic energies which induces radial oscillating motion of the jet. The oscillation timescale is characterized by the initial pressure ratio of the jet to the ambient medium, and follows a simple scaling relation $\tau_{\rm oscillation} \propto (P_{\rm jet,0}/P_{\rm amb,0})^{1/2}$. It is confirmed from extended two-dimensional simulations that this radial oscillating motion in the one-dimensional system manifests as modulation of the structure of the jet in a more realistic situation where a relativistically hot jet propagates through an ambient medium. It is found that when the ambient medium has a power law pressure distribution, the size of the reconfinement region along the propagation direction of the jet in the modulation structure $\lambda$ evolves according to a self-similar relation $\lambda \propto t^{\alpha/2}$ where $\alpha$ is the power-law index of the pressure distribution.
The activity of the Sun as a result of cyclic changes of the global magnetic field is studied. As a consequence of the analysis of magnetic activity of solar-type stars the following power dependencies were found: the dependence between the rotation periods and the effective temperatures, the dependence between the duration of the "11-year" cycles of activity and the effective temperatures, and the dependence between the duration of quasi-biennial cycles and the effective temperatures. It is shown that the physical nature of these dependencies associated with the observed properties of solar-type stars and can be explained by the existence of internal Rossby waves around the base of convective shells of these stars.
We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.
It has been suggested that the transition between magnetorotationally active and dead zones in protoplanetary disks should be prone to the excitation of vortices via Rossby wave instability (RWI). However, the only numerical evidence for this has come from alpha disk models, where the magnetic field evolution is not followed, and the effect of turbulence is parametrized by Laplacian viscosity. We aim to establish the phenomenology of the flow in the transition in 3D resistive-magnetohydrodynamical models. We model the transition by a sharp jump in resistivity, as expected in the inner dead zone boundary, using the Pencil Code to simulate the flow. We find that vortices are readily excited in the dead side of the transition. We measure the mass accretion rate finding similar levels of Reynolds stress at the dead and active zones, at the $\alpha\approx 10^{-2}$ level. The vortex sits in a pressure maximum and does not migrate, surviving until the end of the simulation. A pressure maximum in the active zone also triggers the RWI. The magnetized vortex that results should be disrupted by parasitical magneto-elliptic instabilities, yet it subsists in high resolution. This suggests that either the parasitic modes are still numerically damped, or that the RWI supplies vorticity faster than they can destroy it. We conclude that the resistive transition between the active and dead zones in the inner regions of protoplanetary disks, if sharp enough, can indeed excite vortices via RWI. Our results lend credence to previous works that relied on the alpha-disk approximation, and caution against the use of overly reduced azimuthal coverage on modeling this transition.
We have spectroscopically identified ?about 100 G-, K- and M-type members of the Scorpius Centaurus complex. To deduce the age of these young stars we compare their Li ?6708 absorption line strengths against those of stars in the TW Hydrae association and ?beta Pictoris moving group. These line strengths indicate that Sco-Cen stars are younger than ?beta Pic stars whose ages of ?~12 Myr have previously been derived from a kinematic traceback analysis. Our derived age, ?~10Myr, for stars in the LCC and UCL subgroups of ScoCen is younger than previously published ages based on the moving cluster method and upper main sequence fitting. The discrepant ages are likely due to an incorrect (or lack of) cross-calibration between model-dependent and model-independent age-dating methods.
We have used a combination of broad-band near-infrared and optical Johnson-Cousins photometry to study the dust properties in the line of sight to the Galactic globular cluster M4. We have investigated the reddening effects in terms of absolute strength and variation across the cluster field, as well as the form of the reddening law defined by the type of dust. Here, we determine the ratio of absolute to selective extinction (R_V) in the line of sight towards M4, which is known to be a useful indicator for the type of dust and therefore characterizes the applicable reddening law. Our method is independent of age assumptions and appears to be significantly more precise and accurate than previous approaches. We obtain A_V/E(B-V)=3.76\pm0.07 (random error) for the dust in the line of sight to M4 for our set of filters. With this value, the distance to M4 is found to be 1.80\pm0.05 kpc (random error). A reddening map for M4 has been created, which reveals a spatial differential reddening of delta E(B-V)>0.2 mag across a field within 10' around the cluster centre and a total mean reddening of E(B-V)=0.37\pm0.01. In order to provide accurate zero points for the extinction coefficients of our photometric filters, we investigated the impact of stellar parameters such as temperature, surface gravity and metallicity on the extinction properties in different bandpasses. Using both synthetic ATLAS9 spectra and observed spectral energy distributions, we found similar sized effects for the range of temperature, surface gravity, and metallicity typical of globular cluster stars: each causes a change of about 3% in the necessary correction factor for each filter combination. From our calculations, we provide extinction zero points for Johnson-Cousins and 2MASS filters, spanning a wide range of stellar parameters and dust types, suited for accurate, object-specific extinction corrections.
We present continuum high resolution Submillimeter Array (SMA) observations of the transition disk object RX J1633.9-2442, which is located in the Ophiuchus molecular cloud and has recently been identified as a likely site of ongoing giant planet formation. The observations were taken at 340 GHz (880 micron) with the SMA in its most extended configuration, resulting in an angular resolution of 0.3" (35 AU at the distance of the target). We find that the disk is highly inclined (i ~50 deg) and has an inner cavity ~25 AU in radius, which is clearly resolved by our observations. We simultaneously model the entire optical to millimeter wavelength spectral energy distribution (SED) and SMA visibilities of RX J1633.9-2442 in order to constrain the structure of its disk. We find that an empty cavity ~25 AU in radius is inconsistent with the excess emission observed at 12, 22, and 24 micron. Instead, the mid-IR excess can be modeled by either a narrow, optically thick ring at ~10 AU or an optically thin region extending from ~7 AU to ~25 AU. The inner disk (r < 5 AU) is mostly depleted of small dust grains as attested by the lack of detectable near-IR excess. We also present deep Keck aperture masking observations in the near-IR, which rule out the presence of a companion up to 500 times fainter than the primary star (in K-band) for projected separations in the 5-20 AU range. We argue that the complex structure of the RX J1633.9-2442 disk is best explained by multiple planets embedded within the disk. We also suggest that the properties and incidence of objects such as RX J1633.9-2442, T Cha, and LkCa 15 (and those of the companions recently identified to these two latter objects) are most consistent with the runaway gas accretion phase of the core accretion model, when giant planets gain their envelopes and suddenly become massive enough to open wide gaps in the disk.
This is a very informal report that gives further details on the evidence for a bubble universe based on an anomaly in the angular distribution of quasar magnitudes that was presented in a short paper in arXiv:1202.4433. This report addresses some concerns of two reviewers. It is meant to be read in conjunction with 1202.4433. There is very little overlap between the two articles. This extended discussion is, by necessity, somewhat more technical in nature. I am grateful for the reviewers' comments that forced me to understand these issues more thoroughly.
We present a novel approach to modified theories of gravity that consists of adding to the Einstein-Hilbert Lagrangian an f(R) term constructed a la Palatini. Using the respective dynamically equivalent scalar-tensor representation, we show that the theory can pass the Solar System observational constraints even if the scalar field is very light. This implies the existence of a long-range scalar field, which is able to modify the cosmological and galactic dynamics, but leaves the Solar System unaffected. We also verify the absence of instabilities in perturbations and provide explicit models which are consistent with local tests and lead to the late-time cosmic acceleration.
Brown has recently introduced a covariant formulation of the BSSN equations which is well suited for curvilinear coordinate systems. This is particularly desirable as many astrophysical phenomena are symmetric with respect to the rotation axis or are such that curvilinear coordinates adapt better to their geometry. However, the singularities associated with such coordinate systems are known to lead to numerical instabilities unless special care is taken (e.g., regularization at the origin). Cordero-Carrion will present a rigorous derivation of partially implicit Runge-Kutta methods in forthcoming papers, with the aim of treating numerically the stiff source terms in wave-like equations that may appear as a result of the choice of the coordinate system. We have developed a numerical code solving the BSSN equations in spherical symmetry and the general relativistic hydrodynamic equations written in flux-conservative form. A key feature of the code is that it uses a second-order partially implicit Runge-Kutta method to integrate the evolution equations. We perform and discuss a number of tests to assess the accuracy and expected convergence of the code, namely a pure gauge wave, the evolution of a single black hole, the evolution of a spherical relativistic star in equilibrium, and the gravitational collapse of a spherical relativistic star leading to the formation of a black hole. We obtain stable evolutions of regular spacetimes without the need for any regularization algorithm at the origin.
Chameleon fields, which are scalar field dark energy candidates, can evade fifth force constraints by becoming massive in high-density regions. However, this property allows chameleon particles to be trapped inside a vacuum chamber with dense walls. Afterglow experiments constrain photon-coupled chameleon fields by attempting to produce and trap chameleon particles inside such a vacuum chamber, from which they will emit an afterglow as they regenerate photons. Here we discuss several theoretical and systematic effects underlying the design and analysis of the GammeV and CHASE afterglow experiments. We consider chameleon particle interactions with photons, Fermions, and other chameleon particles, as well as with macroscopic magnetic fields and matter. The afterglow signal in each experiment is predicted, and its sensitivity to various properties of the experimental apparatus is studied. Finally, we use CHASE data to exclude a wide range of photon-coupled chameleon dark energy models.
Axino, related to the SUSY transformation of axion, can mix with goldstino in principle. This case is realized when some superfields carrying nonvanishing Peccei-Quinn charges develop both scalar VEVs and F-terms. In this case, we present a proper definition of axion and axino. With this definition, we present the QCD axino mass in the most general framework, including non-minimal K\"ahler potential. The axino mass is known to have a hierarchical mass structure depending on accidental symmetries. If $G_A=0$ where $G=K+M_P^2\ln|W|^2$, we obtain $m_{\tilde a}=(1+ \epsilon^2)m_{3/2}$ with the axino-gravitino mixing parameter $\epsilon$ in the K\"ahler potential. For $G_A\ne 0$, the axino mass depends on the details of the K\"ahler potential. In the gauge mediation scenario, the gaugino mass is the dominant axino mass parameter. Therefore, we can take the theoretical QCD axino mass as a free parameter in the study of its cosmological effects, ranging from eV to multi-TeV scales, without a present knowledge on its ultraviolet completion.
Quantum and thermal fluctuations of an irrotational fluid are studied across the transition regime connecting a protoinflationary phase of decelerated expansion to an accelerated epoch driven by a single inflaton field. The protoinflationary inhomogeneities are suppressed when the transition to the slow roll phase occurs sharply over space-like hypersurfaces of constant energy density. If the transition is delayed, the interaction of the quasi-normal modes related, asymptotically, to fluid phonons and inflaton quanta leads to an enhancement of curvature perturbations. It is shown that the dynamics of the fluctuations across the protoinflationary boundaries is determined by the monotonicity properties of the pump fields controlling the energy transfer between the background geometry and the quasi-normal modes of the fluctuations. After corroborating the analytical arguments with explicit numerical examples, general lessons are drawn on the classification of the protoinflationary transition.
Links to: arXiv, form interface, find, astro-ph, recent, 1204, contact, help (Access key information)
This work aims to search for the observational footprints of the interactions between the interstellar medium (ISM) and stellar winds in the Wolf-Rayet (WR) nebula NGC6888 in order to understand its ionization structure, chemical composition, and kinematics. We have collected a set of integral field spectroscopy observations across NGC6888, obtained with PPAK in the optical range performing both 2D and 1D analyses. Attending to the 2D analysis in the northeast part of NGC6888, we have generated maps of the extinction structure and electron density. We produced statistical frequency distributions of the radial velocity and diagnostic diagrams. We have found that the spectra of a localized region to the southwest of this pointing can be represented well by shock models. Furthermore, we performed a thorough study of integrated spectra in nine regions over the whole nebula. We derived electron densities ranging from <100 to 360 cm^(-3). The electron temperature varies from ~7700 K to ~10200 K. A strong variation of up to a factor 10 between different regions in the nitrogen abundance has been found: N/H appears lower than the solar abundance in those positions observed at the edges and very enhanced in the observed inner parts. Oxygen appears slightly underabundant with respect to solar value, whereas the helium abundance is found to be above it. Finally, we provide a scenario for the evolution of NGC6888 to explain the features observed. This scheme consists of a structure of multiple shells.
We present the discovery of a source with broadband infrared photometric characteristics similar to Sakurai's Object. WISE J180956.27-330500.2 (hereafter, J1810-3305) shows very red WISE colors, but a very blue 2MASS [K] vs. WISE [W1 (3.4 micron)] color. It was not visible during the IRAS era, but now has a 12 micron flux well above the IRAS point source catalog detection limit. There are also indications of variability in historical optical photographic plates, as well as in multi-epoch AKARI mid-infrared measurements. The broadband infrared spectral energy distribution, post-IRAS brightening and multiwavelength variability are all characteristics also shared by Sakurai's Object - a post asymptotic giant branch (post-AGB) star which underwent a late thermal pulse and recently ejected massive envelopes of dust that are currently expanding and cooling. Optical progenitor colors suggest that J1810-3305 may have been of late spectral class. Its dramatic infrared brightening, and the detection of a late-type optical counterpart are consistent with a scenario in which we have caught an extremely massive dust ejection event (in 1998 or shortly before) during the thermal pulse of an AGB star, thus providing a unique opportunity to observe stellar evolution in this phase. J1810-3305 is the only source in the entire WISE preliminary data release with similar infrared SED and variability, emphasizing the rarity of such sources. Confirmation of its nature is of great importance.
The analysis of relativistically broadened X-ray spectral features from the inner accretion disk provides a powerful tool for measuring the spin of supermassive black holes (SMBH) in active galactic nuclei (AGN). However, AGN spectra are often complex and careful analysis employing appropriate and self-consistent models are required if one is to obtain robust results. In this paper, we revisit the deep July-2009 Suzaku observation of the Seyfert galaxy NGC3783 in order to study in a rigorous manner the robustness of the inferred black hole spin parameter. Using Monte Carlo Markov Chain (MCMC) techniques, we identify a (partial) modeling degeneracy between the iron abundance of the disk and the black hole spin parameter. We show that the data for NGC3783 strongly require both supersolar iron abundance (Z_Fe=2-4Zsun) and a rapidly spinning black hole (a>0.88). We discuss various astrophysical considerations that can affect the measured abundance. We note that, while the abundance enhancement inferred in NGC3783 is modest, the X-ray analysis of some other objects has found extreme iron abundances. We introduce the hypothesis that the radiative levitation of iron ions in the innermost regions of radiation-dominated AGN disks can enhance the photospheric abundance of iron. We show that radiative levitation is a plausible mechanism in the very inner regions of high accretion rate AGN disks.
Theoretical studies of structure formation find an inverse proportionality between the concentration of dark matter haloes and virial mass. This trend has been recently confirmed for virial masses Mvir > ~6e12 Msun by the observation of the X-ray emission from the hot halo gas. We present an alternative approach to this problem, exploring the concentration of dark matter haloes over galaxy scales on a sample of 18 early-type systems. Our c-Mvir relation is consistent with the X-ray analysis, extending towards lower virial masses, covering the range from ~4e11 Msun up to 5e12 Msun. A combination of the lensing analysis along with photometric data allows us to constrain the baryon fraction within a few effective radii, which is compared with prescriptions for adiabatic contraction (AC) of the dark matter haloes. We find that the standard methods for AC are strongly disfavored, requiring additional mechanisms -- such as mass loss during the contraction process -- to play a role during the phases following the collapse of the haloes.
We report the discovery of a new ultra-faint globular cluster in the constellation of Ursa Minor, based on stellar photometry from the MegaCam imager at the Canada-France-Hawaii Telescope (CFHT). We find that this cluster, Munoz 1, is located at a distance of 45 +/- 5 kpc and at a projected distance of only 45 arcmin from the center of the Ursa Minor dSph galaxy. Using a Maximum Likelihood technique we measure a half-light radius of 0.5 arcmin, or equivalently 7 pc and an ellipticity consistent with being zero. We estimate its absolute magnitude to be M_V=-0.4 +/- 0.9, which corresponds to L_V=120 (+160, -65) L_sun and we measure a heliocentric radial velocity of -137 +/- 4 km/s based on Keck/DEIMOS spectroscopy. This new satellite is separate from Ursa Minor by ~30 kpc and 110 km/s suggesting the cluster is not obviously associated with the dSph, despite the very close angular separation. Based on its photometric properties and structural parameters we conclude that Munoz 1 is a new ultra-faint stellar cluster. Along with Segue 3 this is one of the faintest stellar clusters known to date.
A deep narrow-band survey for Ly-alpha emission carried out on the VLT-FORS2 has revealed 98 Ly-alpha candidates down to a flux limit of 4.e-18 erg/s/cm^2 in a volume of 5500 comoving Mpc^3 at z=2.4 centered on the hyperluminous quasar HE0109-3518. The properties of the detected sources in terms of their i) equivalent width distribution, ii) luminosity function, and iii) the average luminosity versus projected distance from the quasar, all suggest that a large fraction of these objects have been fluorescently "illuminated" by HE0109-3518. This conclusion is supported by comparison with detailed radiative transfer simulations of the effects of the quasar illumination. 18 objects have a rest-frame Equivalent Width (EW0) larger than 240A, the expected limit for Ly-alpha emission powered by Population II star formation and 12 sources among these do not have any continuum counterpart in a deep V-band imaging of the same field. For these, a stacking analysis indicates EW0>800A, effectively ruling out Ly-alpha powered by internal star formation. These sources are thus the best candidates so far for proto-galactic clouds or "dark" galaxies at high-redshift, whose existence has recently been suggested by several theoretical studies. Assuming they are mostly ionized by the quasar radiation, we estimate that their gas masses would be about 10^9 Msun implying that their star formation efficiencies (SFE) are less than 4.e-12 yr^-1 one order of magnitude below the SFE of the most gas-rich dwarf galaxies locally, and five hundred times lower than typical massive star-forming galaxies at z~2. We have also discovered extended, filamentary gas, also likely illuminated by the quasar, around some of the brightest continuum-detected sources with EW0>240A. This emission is compatible with the expectations for circum-galactic cold streams but other origins, including tidal stripping, are also possible.
We have observed a sample of 19 carbon stars in the Sculptor, Carina, Fornax, and Leo I dwarf spheroidal galaxies with the Infrared Spectrograph on the Spitzer Space Telescope. The spectra show significant quantities of dust around the carbon stars in Sculptor, Fornax, and Leo I, but little in Carina. Previous comparisons of carbon stars with similar pulsation properties in the Galaxy and the Magellanic Clouds revealed no evidence that metallicity affected the production of dust by carbon stars. However, the more metal-poor stars in the current sample appear to be generating less dust. These data extend two known trends to lower metallicities. In more metal-poor samples, the SiC dust emission weakens, while the acetylene absorption strengthens. The bolometric magnitudes and infrared spectral properties of the carbon stars in Fornax are consistent with metallicities more similar to carbon stars in the Magellanic Clouds than in the other dwarf spheroidals in our sample. A study of the carbon budget in these stars reinforces previous considerations that the dredge-up of sufficient quantities of carbon from the stellar cores may trigger the final superwind phase, ending a star's lifetime on the asymptotic giant branch.
We investigate where brightest cluster galaxies (BCGs) sit on the fundamental plane of black hole (BH) activity, an established relation between the X-ray luminosity, the radio luminosity and the mass of a BH. Our sample mostly consists of BCGs that lie at the centres of massive, strong cooling flow clusters, therefore requiring extreme mechanical feedback from their central active galactic nucleus (AGN) to offset cooling of the intracluster plasma (L_mech>10^44-45 erg/s). Based on the BH masses derived from the M_BH-sigma and M_BH-M_K correlations, we find that all of our objects are offset from the plane such that they appear to be less massive than predicted from their X-ray and radio luminosities (to more than a 99 per cent confidence level). For these objects to be consistent with the fundamental plane, the M_BH-sigma and M_BH-M_K correlations therefore seem to underestimate the BH masses of BCGs, on average by a factor of 10. Our results suggest that the standard relationships between BH mass and host galaxy properties no longer hold for these extreme galaxies. Furthermore, our results imply that if these BHs follow the fundamental plane, then many of those that lie in massive, strong cool core clusters must be ultramassive with M_BH>10^10M_sun. This exceeds the largest BH masses known and has important ramifications for our understanding of the formation and evolution of BHs.
The Alcock-Paczynski (AP) effect uses the fact that, when analyzed with the correct geometry, we should observe structure that is statistically isotropic in the Universe. For structure undergoing cosmological expansion with the background, this constrains the product of the Hubble parameter and the angular diameter distance. However, the expansion of the Universe is inhomogeneous and local curvature depends on density. We argue that this distorts the AP effect on small scales. After analyzing the dynamics of galaxy pairs in the Millennium simulation, we find an interplay between peculiar velocities, galaxy properties and local density that affects how pairs trace cosmological expansion. We find that only low mass, isolated galaxy pairs trace the average expansion with a minimum "correction" for peculiar velocities. Other pairs require larger, more cosmology and redshift dependent peculiar velocity corrections and, in the small-separation limit of being bound in a collapsed system, do not carry cosmological information.
We develop a new test of local bias, by constructing a locally biased halo density field from sampling the dark matter-halo distribution. Our test differs from conventional tests in that it preserves the full scatter in the bias relation and it does not rely on perturbation theory. We put forward that bias parameters obtained using a smoothing scale R can only be applied to computing the halo power spectrum at scales k ~ 1/R. Our calculations can automatically include the running of bias parameters and give vanishingly small loop corrections at low-k. Our proposal results in much better agreement of the sampling and perturbation theory results with the original simulations. In particular, unlike the standard interpretation of local bias in the literature, our treatment of local bias does not generate a constant power at very large scales, in agreement with our measurements in the simulations of the halo power spectrum at wavenumbers below its turn-over. Using perturbation theory and our non-perturbative sampling technique we also demonstrate that nonlocal bias effects discovered recently in simulations impact the power spectrum only at the few percent level in the weakly nonlinear regime.
Stellar streams provide unique probes of galactic potentials, with the longer streams normally providing the cleaner measurements. In this paper, we show an example of a short tidal stream that is particularly sensitive to the shape of the Milky Way's dark matter halo: the globular cluster tidal stream NGC 5466. This stream has an interesting deviation from a smooth orbit at its western edge. We show that such a deviation favours an underlying oblate or triaxial halo (irrespective of plausible variations in the Milky Way disc properties and the specific halo parametrisation chosen); spherical or prolate halo shapes can be excluded at a high confidence level. Therefore, more extensive data sets along the NGC 5466 tidal stream promise strong constraints on the Milky Way halo shape.
A search for emission lines in foreground galaxies in quasar spectra (z(gal) < z(QSO)) of the Sloan Digital Sky Survey (SDSS) data release 5 (DR5) reveals 23 examples of quasars shining through low redshift, foreground galaxies at small impact parameters (< 10 kpc). About 74,000 quasar spectra were examined by searching for narrow H{\alpha} emission lines at z < 0.38, at a flux level greater than 5 \times 10^-17 ergs cm^-2 s^-1, then confirming that other expected emission lines of the H II regions in the galaxy are detected. The galaxies were deblended from the quasar images to get colors and morphologies. For cases that allow the galaxy and the quasar to be deblended, the galaxies are blue (0.95 <(u-r)< 1.95). Extinction and reddening through the galaxies is determined from the (g-i) color excesses of the quasars. These reddening values are compared with the flux ratio of H{\alpha} to H{\beta}, which reflect the extinction for an undetermined fraction of the sightline through each galaxy. No trends were found relating E(B-V)_(g-i), impact parameter (b), and (u-r) for the galaxies or between E(B-V) derived from (g-i) and that derived from H{\alpha}/H{\beta}. Comparison with previous studies of quasar absorption systems indicate our sample is more reddened, suggesting disk-dominated absorber galaxies. Measurement or limits on galactic, interstellar Ca II and Na I absorption lines are given from the quasar spectrum. No trends were found relating Ca II equivalent width (W (Ca II)) or Na I equivalent width (W (Na I)) to b, but a correlation of r_s = -0.77 ({\alpha} = 0.05) was found relating W (Ca II) and E(B-V)(g-i) .
We have searched for [OIII] 5007 emission in high resolution spectroscopic data from Flames/Giraffe VLT observations of 174 massive globular clusters (GCs) in NGC4472. No planetary nebulae (PNe) are observed in these clusters, constraining the number of PNe per bolometric luminosity, \alpha<0.8*10^{-7}PN/L_{\odot}. This is significantly lower than the rate predicted from stellar evolution, if all stars produce PNe. Comparing our results to populations of PNe in galaxies, we find most galaxies have a higher \alpha than these GCs (more PNe per bolometric luminosity - though some massive early-type galaxies do have similarly low \alpha). The low \alpha required in these GCs suggests that the number of PNe per bolometric luminosity does not increase strongly with decreasing mass or metallicity of the stellar population. We find no evidence for correlations between the presence of known GC PNe and either the presence of low mass X-ray binaries (LMXBs) or the stellar interaction rates in the GCs. This, and the low \alpha observed, suggests that the formation of PNe may not be enhanced in tight binary systems. These data do identify one [OIII] emission feature, this is the (previously published) broad [OIII] emission from the cluster RZ 2109. This emission is thought to originate from the LMXB in this cluster, which is accreting at super-Eddington rates. The absence of any similar [OIII] emission from the other clusters favors the hypothesis that this source is a black hole LMXB, rather than a neutron star LMXB with significant geometric beaming of its X-ray emission.
A model for strong kinetic Alfv\'en plasma turbulence at scales smaller than the ion gyroscale is proposed. It is argued that magnetic and density fluctuations are concentrated mostly at two-dimensional structures, which leads to their Fourier energy spectra $E(k_\perp)\propto k_{\perp}^{-8/3}$, where $k_\perp$ is the wave-vector component normal to the strong background magnetic field. The results are shown to be in good agreement with numerical simulations, and they can explain recent observations of magnetic and density fluctuations in the solar wind at sub-proton scales.
We have analysed optical, near-, and mid-IR images of the bright planetary nebula (PN) IC 418. These images probe unambiguously for the first time a number of low surface brightness structures and shells around the bright main nebula, including radial filaments or rays, a system of three concentric rings, and two detached haloes with sizes ~ 150" and 220"\times250", respectively. The main nebula is slightly off-centered with respect to the elliptical outer halo. The time-lapse between the two haloes is 10,000-50,000 yr, whereas the time-lapse between the three concentric rings is ~ 630 yr. We emphasize the advantages of near- and mid-IR imaging for the detection of faint structures and shells around optically bright nebulae.
We have succeeded in obtaining magnetized star models that have extremely strong magnetic fields in the interior of the star. In our formulation, arbitrary functions of the magnetic flux function appear in the expression for the current density. By appropriately choosing the functional form for one of the arbitrary functions that corresponds to the distribution of the toroidal current density, we have obtained configurations with magnetic field distributions that are highly localized within the central part and near the magnetic axis region. The absolute values of the central magnetic fields are stronger than those of the surface region by two orders of magnitude. By applying our results to magnetars, the internal magnetic poloidal fields could be 10^17 G, although the surface magnetic fields are about 10^15 G in the case of magnetars. For white dwarfs, the internal magnetic poloidal fields could be 10^12 G, when the surface magnetic fields are 10^9-10^10 G .
We report multiwavelength observations of the black hole transient GX 339-4 during its outburst decay in 2011 using the data from RXTE, Swift and SMARTS. Based on the X-ray spectral, temporal, and the optical/infrared (OIR) properties, the source evolved from the soft-intermediate to the hard state. Twelve days after the start of the transition towards the hard state, a rebrightening was observed simultaneously in the optical and the infrared bands. Spectral energy distributions (SED) were created from observations at the start, and close to the peak of the rebrightening. The excess OIR emission above the smooth exponential decay yields flat spectral slopes for these SEDs. Assuming that the excess is from a compact jet, we discuss the possible locations of the spectral break that mark the transition from optically thick to optically thin synchrotron components. Only during the rising part of the rebrightening, we detected fluctuations with the binary period of the system. We discuss a scenario that includes irradiation of the disk in the intermediate state, irradiation of the secondary star during OIR rise and jet emission dominating during the peak to explain the entire evolution of the OIR light curve.
MAXI J1659-152 is a bright X-ray transient black-hole candidate binary system discovered in September 2010. We report here on MAXI, RXTE, Swift, and XMM-Newton observations during its 2010/2011 outburst. We find that during the first one and a half week of the outburst the X-ray light curves display drops in intensity at regular intervals, which we interpret as absorption dips. About three weeks into the outbursts, again drops in intensity are seen. These dips have, however, a spectral behaviour opposite to that of the absorption dips, and are related to fast spectral state changes (hence referred to as transition dips). The absorption dips recur with a period of 2.414+/-0.005 hrs, which we interpret as the orbital period of the system. This implies that MAXI J1659-152 is the shortest period black-hole candidate binary known to date. The inclination of the accretion disk with respect to the line of sight is estimated to be 65-80 degrees. We suggest the companion to the black-hole candidate to be an M5 dwarf star, with a mass and radius of about 0.2-0.3 M_sun and 0.2-0.25 R_sun, respectively. We find that the companion had an initial mass of about 1.5 M_sun, which evolved to its current mass in about 5-6 billion years. The system is rather compact (orbital separation of about 1.33 R_sun), and is located at a distance of about 8.6 kpc, with a height above the Galactic plane of about 2.4 kpc. The characteristics of short orbital period and high Galactic scale height are shared with two other transient black-hole candidate X-ray binaries, i.e., XTE J1118+480 and Swift J1735.5-0127. All three are kicked out of the Galactic plane into the halo, rather than being formed in a globular cluster.
Microflares are small activities in solar low atmosphere, some are in the low corona, and others in the chromosphere. Observations show that some of the microflares are triggered by magnetic reconnection between emerging flux and a pre-existing background magnetic field. We perform 2.5D compressible resistive MHD simulations of magnetic reconnection with gravity considered. The background magnetic field is a canopy-type configuration which is rooted at the boundary of the solar supergranule. By changing the bottom boundary conditions in the simulation, new magnetic flux emerges up at the center of the supergranule and reconnects with the canopy-type magnetic field. We successfully simulate the coronal and chromospheric microflares, whose current sheets are located at the corona and the chromosphere, respectively. The microflare of coronal origin has a bigger size and a higher temperature enhancement than that of chromospheric origin. In the microflares of coronal origin, we also found a hot jet ($\sim$$1.8 \times 10^6$ K), which is probably related to the observational EUV/SXR jets, and a cold jet ($\sim$$10^4$ K), which is similar to the observational H$\alpha$/Ca surges, whereas there is only an H$\alpha$/Ca bright point in the microflares of chromospheric origin. The study of parameter dependence shows that the size and strength of the emerging magnetic flux are the key parameters which determine the height of the reconnection location, and further determine the different observational features of the microflares.
A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for computing the MHD part, namely, modified Mac Cormack Scheme (MMC), Lax-Friedrichs scheme (LF) and weighted essentially non-oscillatory (WENO) scheme. All of them are second order, two-step, component-wise schemes for hyperbolic conservative equations. The total variation diminishing (TVD) limiters and approximate Riemann solvers are also equipped. A high resolution can be achieved by the hierarchical block-structured AMR mesh. We use the extended generalized Lagrange multiplier (EGLM) MHD equations to reduce the non-divergence free error produced by the scheme in the magnetic induction equation. The numerical algorithms for the non-ideal terms, e.g., the resistivity and the thermal conduction, are also equipped in the MAP code. The details of the AMR and MPI algorithms are described in the paper.
The spin of Cygnus X-1 is measured by fitting reflection models to Suzaku data covering the energy band 0.9-400 keV. The inner radius of the accretion disc is found to lie within 2 gravitational radii (r_g=GM/c^2) and a value for the dimensionless black hole spin is obtained of 0.97^{+0.014}_{-0.02}. This agrees with recent measurements using the continuum fitting method by Gou et al. and of the broad iron line by Duro et al. The disc inclination is measured at 23.7^{+6.7}_{-5.4} deg, which is consistent with the recent optical measurement of the binary system inclination by Orosz et al of 27+/-0.8 deg. We pay special attention to the emissivity profile caused by irradiation of the inner disc by the hard power-law source. The X-ray observations and simulations show that the index q of that profile deviates from the commonly used, Newtonian, value of 3 within 3r_g, steepening considerably within 2r_g, as expected in the strong gravity regime.
The High Energy Stereoscopic System (H.E.S.S.) is an array of four imaging
atmospheric-Cherenkov telescopes located in Namibia and designed to detect
extensive air showers initiated by gamma-rays in the very-high-energy domain.
It is an ideal instrument for surveying the Galactic plane in search of new
sources, thanks to its location in the Southern Hemisphere, its excellent
sensitivity, and its large field-of-view. The efforts of the H.E.S.S. Galactic
Plane Survey, the first comprehensive survey of the inner Galaxy at TeV
energies, have contributed to the discovery of an unexpectedly large and
diverse population of over 60 sources of VHE gamma rays within its current
range of l=250 to 65 degrees in longitude and |b|<=3.5 degrees in latitude. The
population of VHE gamma-ray emitters is dominated by the pulsar wind nebula and
supernova remnant source classes, although nearly a third remain unidentified
or confused.
The sensitivity of H.E.S.S. to sources in the inner Galaxy has improved
significantly over the past two years, from continued survey observations,
dedicated follow-up observations of interesting source candidates, and from the
development of advanced methods for discrimination of gamma-ray-induced showers
from the dominant background of hadron-induced showers. The latest maps of the
Galaxy at TeV energies will be presented, and a few remarkable new sources will
be highlighted.
Recent cosmological data favour additional relativistic degrees of freedom beyond the three active neutrinos and photons, often referred to as 'dark' radiation. Light sterile neutrinos is one of the prime candidates for such additional radiation. However, constraints on sterile neutrinos based on the current cosmological data have been derived using simplified assumptions about thermalisation of the sterile neutrino at the Big Bang Nucleosynthesis (BBN) epoch. These assumptions are not necessarily justified and here we solve the full quantum kinetic equations in the (1 active + 1 sterile) scenario and derive the number of thermalised species just before BBN begins (T~1MeV) for null (L=0) and large (L=0.01) initial lepton asymmetry and for a range of possible mass-mixing parameters. We find that the full thermalisation assumption during the BBN epoch is justified for initial small lepton asymmetry only. Partial or null thermalisation occurs when the initial lepton asymmetry is large.
I'm trying to summarize the science communicated via oral presentations and by posters at the IAU Symposium 285 "Cosmic Masers - from OH to H_0", which took place from January 29 to February 3, 2012 in Stellenbosch, South Africa.
The ANTARES telescope is the largest underwater neutrino telescope existing at present. It is based on the detection of Cherenkov light produced in sea water by neutrino-induced muons. The detector, consisting of a tri-dimensional array of 885 photomultipliers arranged on twelve vertical lines, is located at a depth of 2475 m in the Mediterranean Sea, 40 km off the French coast. The main goal of the experiment is to probe the Universe by means of neutrino events in an attempt to investigate the nature of high energy astrophysical sources, to contribute to the identification of cosmic ray sources, and to explore the nature of dark matter. In this contribution we will review the status of the detector, illustrate its operation and performance, and present the first results from the analysis carried out on atmospheric muons and neutrinos, as well as from the search for astrophysical neutrino sources.
The ANTARES neutrino telescope is installed at a depth of 2.5 km of the Mediterranean Sea and consists of a three-dimensional array of 885 photomultipliers arranged on twelve detector lines. The prime objective is to detect high-energy neutrinos from extraterrestrial origin. Relativistic muons emerging from charged-current muon neutrino interactions in the detector surroundings produce a cone of Cerenkov light which allows the reconstruction of the original neutrino direction. The collaboration has implemented different methods to search for neutrino point sources in the data collected since 2007. Results obtained with these methods as well as the sensitivity of the telescope are presented.
Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.
The observation that old, massive galaxies have a larger fraction of low-mass stars than their younger, lower-mass counterparts adds to mounting evidence that star formation may have been different in the early Universe.
Any non-gravitational coupling between dark matter and dark energy modifies the cosmological dynamics. Interactions in the dark sector are considered to be relevant to address the coincidence problem. Moreover, in various models the observed accelerated expansion of the Universe is a pure interaction phenomenon. Here we review recent approaches in which a coupling between both dark components is crucial for the evolution of the Universe.
The possibility of the Ly{\beta} fluorescence mechanism being operational in classical Be stars and thereby contributing to the strength of the OI 8446 {\AA} line has been recognized for long. However this supposition needs to be quantified by comparing observed and predicted OI line ratios. In the present work, optical and near-infrared spectra of classical Be stars are presented. We analyse the observed strengths of the OI 7774, 8446, 11287 and 13165 {\AA} lines which have been theoretically proposed as diagnostics for identifying the excitation mechanism. We have considered and examined the effects of Ly{\beta} fluorescence, collisional excitation, recombination and continuum fluorescence on these OI line strengths. From our analysis it appears that the Ly{\beta} fluorescence process is indeed operative in Be stars.
In the course of the XMM-Newton survey of the Small Magellanic Cloud (SMC), two new bright X-ray sources were discovered exhibiting the spectral characteris- tics of High Mass X-ray Binaries - but revealing only weak evidence for pulsations in just one of the objects(at 153s in XMMUJ010743.1-715953). The accurate X- ray source locations permit the identi?cation of these X-ray source with Be stars, thereby strongly suggesting these systems are new Be/X-ray binaries. From blue spectra the proposed classi?cation for XMMUJ010633.1-731543 is B0.5-1Ve and for XMMUJ010743.1-715953 it is B2IV-Ve.
We analysed data from five XMM-Newton observations of GX 13+1 to investigate the variability of the photo-ionised absorber present in this source. We fitted EPIC and RGS spectra obtained from the "least-variable" intervals with a model consisting of disc-blackbody and blackbody components together with a Gaussian emission feature at ~6.55-6.7 keV modified by absorption due to cold and photo-ionised material. We found a significant correlation between the hard, ~6-10 keV, flux, the ionisation and column density of the absorber and the equivalent width of the broad iron line. We interpret the correlation in a scenario in which a disc wind is thermally driven at large, ~10^{10} cm, radii and the broad line results from reprocessed emission in the wind and/or hot atmosphere. The breadth of the emission line is naturally explained by a combination of scattering, recombination and fluorescence processes. We attribute the variations in the absorption and emission along the orbital period to the view of different parts of the wind, possibly located at slightly different inclination angles. We constrain the inclination of GX 13+1 to be between 60 and 80 degrees from the presence of strong absorption in the line of sight, that obscures up to 80% of the total emission in one observation, and the absence of eclipses. We conclude that the presence of a disc wind and/or a hot atmosphere can explain the current observations of narrow absorption and broad iron emission features in neutron star low mass X-ray binaries as a class.
On various occasions, several authors suggested to model bodily tides with
superposition of two symmetrical bulges. One bulge is always aimed at the
secondary, and thus implements the instantaneous reaction of the primary's
shape and potential to the tide-rising gravitational pull exerted on it by the
secondary. This portion of the tide is called "adiabatic tide" (Zahn 1966a,b)
or "elastic tide" (Ferraz Mello 2012; Krasinsky 2006). The second bulge is
assumed to align orthogonally to the direction to the tide-raising secondary.
So this bulge is set to implement the entire nonelastic portion of the
primary's deformation. This, second bulge is called "dissipative tide" (Zahn
1966a,b; Krasinsky 2006) or "creep tide" (Ferraz Mello 2012).
We demonstrate that the two-bulge method is not a separate approximation, but
ensues directly from the Fourier expansion of a linear tidal theory equipped
with an arbitrary rheological model involving a departure from elasticity.
While less economical mathematically, the two-bulge approach has a good
illustrative power, and may be employed (like, e.g., in Remus et al. 2012a,b)
on a par with a more concise method of complex amplitudes.
We show how the correlation between the thermal Sunyaev Zel'dovich effect (tSZ) from galaxy clusters and dust emission from cosmic infrared background (CIB) sources can be calculated in a halo model framework. Using recent tSZ and CIB models, we find that the size of the tSZ x CIB cross-correlation is approximately 10 per cent at 150 GHz. The contribution to the total angular power spectrum is of order -1 \mu K^2 at ell=3000, however, this value is uncertain by a factor of two to three, primarily because of CIB source modelling uncertainties. We expect the large uncertainty in this component to degrade upper limits on the kinematic Sunyaev Zel'dovich effect (kSZ), due to similarity in the frequency dependence of the tSZ x CIB and kSZ across the frequency range probed by current Cosmic Microwave Background missions. We also find that the degree of tSZ x CIB correlation is higher for mm x sub-mm spectra than mm x mm, because more of the sub-mm CIB originates at lower redshifts (z<2), where most tSZ clusters are found.
ANTARES is the first undersea neutrino detector ever built and presently the neutrino telescope with the largest effective area operating in the Northern Hemisphere. A three- dimensional array of photomultiplier tubes detects the Cherenkov light induced by the muons produced in the interaction of high energy neutrinos with the matter surrounding the detector. The detection of astronomical neutrino sources is one of the main goals of ANTARES. The search for point-like neutrino sources with the ANTARES telescope is described and the preliminary results obtained with data collected from 2007 to 2010 are shown. No cosmic neutrino source has been observed and neutrino flux upper limits have been calculated for the most promising source candidates.
GAIA will provide observations of some multiple asteroid and dwarf systems. These observations are a way to determine and improve the quantification of dynamical parameters, such as the masses and the gravity fields, in these multiple systems. Here we investigate this problem in the cases of Pluto's and Eugenia's system. We simulate observations reproducing an approximate planning of the GAIA observations for both systems, as well as the New Horizons observations of Pluto. We have developed a numerical model reproducing the specific behavior of multiple asteroid system around the Sun and fit it to the simulated observations using least-square method, giving the uncertainties on the fitted parameters. We found that GAIA will improve significantly the precision of Pluto's and Charon's mass, as well as Petit Prince's orbital elements and Eugenia's polar oblateness.
This paper develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating, especially impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. The new features in EBTEL are important for accurate tracking of, in particular, the density. The 0D results are compared to a 1D hydro code (Hydrad) with generally good agreement. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function, (b) extensive parameter surveys and (c) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.
Using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope, we have obtained high-resolution ultraviolet spectra of the C I absorption toward two stars behind the Local Leo Cold Cloud (LLCC). At a distance (~20 pc) that places it well inside the Local Bubble, the LLCC is the nearest example of the coldest known (T~20 K) diffuse interstellar clouds. The STIS measurements of the C I fine-structure excitation toward HD 85259 and HD 83023 indicate that the thermal gas pressure of the LLCC is much greater than that of the warm clouds in the Local Bubble. The mean LLCC pressure measured toward these two stars (60,000 cm$^{-3}$ K) implies an H I density of ~3000 cm$^{-3}$ and a cloud thickness of ~200 AU at the 20 K cloud temperature. Such a thin, cold, dense structure could arise at the collision interface between converging flows of warm gas. However, the measured LLCC pressure is appreciably higher than that expected in the colliding cloud interpretation given the velocity and column density constraints on warm clouds in the HD 85259 and HD 83023 sightlines. Additional STIS measurements of the Zn II, Ni II, and Cr II column densities toward HD 85259 indicate that the LLCC has a modest "warm cloud" dust depletion pattern consistent with its low dust-to-gas ratio determined from H I 21 cm and 100 micron observations. In support of the inferred sheet-like geometry for the LLCC, a multi-epoch comparison of the Na I absorption toward a high-proper-motion background star reveals a 40% column density variation indicative of LLCC Na I structure on a scale of ~50 AU.
We have identified and quantified semi-major axis drifts in Near-Earth Asteroids (NEAs) by performing orbital fits to optical and radar astrometry of all numbered NEAs. We focus on a subset of 54 NEAs that exhibit some of the most reliable and strongest drift rates. Our selection criteria include a Yarkovsky sensitivity metric that quantifies the detectability of semi-major axis drift in any given data set, a signal-to-noise metric, and orbital coverage requirements. In 42 cases, the observed drifts (~10^-3 AU/Myr) agree well with numerical estimates of Yarkovsky drifts. This agreement suggests that the Yarkovsky effect is the dominant non-gravitational process affecting these orbits, and allows us to derive constraints on asteroid physical properties. In 12 cases, the drifts exceed nominal Yarkovsky predictions, which could be due to inaccuracies in our knowledge of physical properties, faulty astrometry, or modeling errors. If these high rates cannot be ruled out by further observations or improvements in modeling, they would be indicative of the presence of an additional non-gravitational force, such as that resulting from a loss of mass of order a kilogram per second. We define the Yarkovsky efficiency f_Y as the ratio of the change in orbital energy to incident solar radiation energy, and we find that typical Yarkovsky efficiencies are ~10^-5.
The solar chromosphere is weakly ionized and interactions between ionized particles and neutral particles likely have significant consequences for the thermodynamics of the plasma. We investigate the importance of introducing neutral particles using numerical 2.5D radiative MHD simulations obtained with the Bifrost code. The models span from the upper layers of the convection zone to the low corona, and solve the full MHD equations with non-grey and NLTE radiative transfer, and thermal conduction. The effects of partial ionization are implemented using the generalized Ohm's law. The approximations required in going from three fluids to the generalized Ohm's law are tested in our simulations. The Ohmic diffusion, the Hall term, and ambipolar diffusion show strong variations in the chromosphere. These strong variations of the various magnetic diffusivities are absent or significantly underestimated when, as has been common for these types of studies, using the VAL-C model as a basis for estimates. In addition, we find that differences in estimating the magnitude of ambipolar diffusion arise depending on which method is used to calculate the ion-neutral collision frequency. These differences cause uncertainties in the different magnetic diffusivity terms. In the chromosphere, we find that the ambipolar diffusion is of the same order of magnitude or even larger than the numerical diffusion used to stabilize our code. As a consequence, ambipolar diffusion produces a strong impact on the modeled atmosphere. Perhaps more importantly, it suggests that at least in the chromospheric domain, self-consistent simulations of the solar atmosphere driven by magneto-convection can accurately describe the impact of the dominant form of resistivity (ambipolar diffusion). This suggests that such simulations may be more realistic in their approach to the lower solar atmosphere than previously assumed.
Cold fronts in cluster cool cores should be erased on short timescales by thermal conduction, unless protected by magnetic fields that are "draped" parallel to the front surfaces, suppressing conduction perpendicular to the fronts. We present MHD simulations of cold front formation in the core of a galaxy cluster with anisotropic thermal conduction, exploring a parameter space of conduction strengths parallel and perpendicular to the field lines. Including conduction has a strong effect on the temperature of the core and the cold fronts. Though magnetic field lines are draping parallel to the front surfaces, the temperature jumps across the fronts are nevertheless reduced. The field geometry is such that the cold gas below the front surfaces can be connected to hotter regions outside via field lines along directions perpendicular to the plane of the sloshing motions and along sections of the front which are not perfectly draped. This results in the heating of this gas below the front on a timescale of a Gyr, but the sharpness of the density and temperature jumps may still be preserved. By modifying the density distribution below the front, conduction may indirectly aid in suppressing Kelvin-Helmholtz instabilities. If conduction along the field lines is unsuppressed, we find that the characteristic sharp jumps in X-ray emission seen in observations of clusters do not form. This suggests that the presence of sharp cold fronts in hot clusters could be used to place upper limits on conduction in the {\it bulk} of the ICM. Finally, the combination of sloshing and anisotropic thermal conduction can result in a larger flux of heat to the core than either process in isolation. While still not sufficient to prevent a cooling catastrophe in the very central ($r \sim$ 5 kpc) regions of the cool core, it reduces significantly the mass of cool gas that accumulates outside those radii.
We report on the first large-scale production of low radioactivity argon from underground gas wells. Low radioactivity argon is of general interest, in particular for the construction of large scale WIMP dark matter searches and detectors of reactor neutrinos for non-proliferation efforts. Atmospheric argon has an activity of about 1 Bq/kg from the decays of 39Ar; the concentration of 39Ar in the underground argon we are collecting is at least a factor of 100 lower than this value. The argon is collected from a stream of gas from a CO2 well in southwestern Colorado with a Vacuum Pressure Swing Adsorption (VPSA) plant. The gas from the well contains argon at a concentration of 400-600 ppm, and the VPSA plant produces an output stream with an argon concentration at the level of 30,000-50,000 ppm (3-5%) in a single pass. This gas is sent for further processing to Fermilab where it is purified by cryogenic distillation. The argon production rate is presently 0.5 kg/day.
Based on the probability distributions of products of random variables, we propose a simple stringy mechanism that prefers the meta-stable vacua with a small cosmological constant. We state some relevant properties of the probability distributions of functions of random variables. We then illustrate the mechanism within the flux compactification models in Type IIB string theory. As a result of the stringy dynamics, we argue that the generic probability distribution for the meta-stable vacua typically peaks with a divergent behavior at the zero value of the cosmological constant. However, its suppression in the single modulus model studied here is modest.
We discuss the possibility of cogenesis generating the ratio of baryon asymmetry to dark matter in a Stueckelberg U(1) extension of the standard model and of the minimal supersymmetric standard model. For the U(1) we choose $L_{\mu}-L_{\tau}$ which is anomaly free and can be gauged. The dark matter candidate arising from this extension is a singlet of the standard model gauge group but is charged under $L_{\mu}-L_{\tau}$. Solutions to the Boltzmann equations for relics in the presence of asymmetric dark matter are discussed. It is shown that the ratio of the baryon asymmetry to dark matter consistent with the current WMAP data, i.e., the cosmic coincidence, can be successfully explained in this model with the depletion of the symmetric component of dark matter from resonant annihilation via the Stueckelberg gauge boson. For the extended MSSM model it is shown that one has a two component dark matter picture with asymmetric dark matter being the dominant component and the neutralino being the subdominant component (i.e., with relic density a small fraction of the WMAP cold dark matter value). Remarkably, the subdominant component can be detected in direct detection experiments such as SuperCDMS and XENON-100. Further, it is shown that the class of Stueckelberg models with a gauged $L_{\mu}-L_{\tau}$ will produce a dramatic signature at a muon collider with the $\sigma(\mu^+\mu^-\to \mu^+\mu^-,\tau^+\tau^-)$ showing a detectable $Z'$ resonance while $\sigma(\mu^+\mu^-\to e^+e^-)$ is devoid of this resonance. Asymmetric dark matter arising from a $U(1)_{B-L}$ Stueckelberg extension is also briefly discussed. Finally, in the models we propose the asymmetric dark matter does not oscillate and there is no danger of it being washed out from oscillations.
We develop a systematic approach to the calculation of scattering cross sections in theories with violation of the Lorentz invariance taking into account the whole information about the theory Lagrangian. As an illustration we derive the Feynman rules and formulas for sums over polarizations in spinor electrodynamics with Lorentz violating operators of dimensions four and six. These rules are applied to compute the probabilities of several astrophysically relevant processes. We calculate the rates of photon decay and vacuum Cherenkov radiation along with the cross sections of electron-positron pair production on background radiation and in the Coulomb field. The latter process is essential for detection of photon-induced air showers in the atmosphere.
We introduce an effective one-fluid description of the interacting dark sector in a spatially flat Friedmann-Robertson-Walker space-time and investigate the stability of the power-law solutions. We find the "source equation" for the total energy density and determine the energy density of each dark component. We study linear and nonlinear interactions which depend on the dark matter and dark energy densities, their first derivatives, the total energy density with its derivatives up to second order and the scale factor. We solve the evolution equations of the dark components for both interactions, examine exhaustively several examples and show cases where the problem of the coincidence is alleviated. We show that a generic nonlinear interaction gives rise to the "relaxed Chaplygin gas model" whose effective equation of state includes the variable modified Chaplygin gas model while some others nonlinear interactions yield de Sitter and power-law scenarios.
We expose a new symmetry for linear perturbations around a solution of non-linear Fierz-Pauli massive gravity plus a bare cosmological constant. The cosmological constant is chosen such that the background metric is flat while the Stuckelberg fields have a non-trivial profile. Around this background, at linear order the new symmetry reduces the propagating degrees of freedom to those of General Relativity, namely the massless helicity 2 modes only. We discuss the physical consequences and possible applications of these findings.
The nucleation process of quark matter in cold (T = 0) stellar matter is investigated using the microscopic Brueckner-Hartree-Fock approach to describe the hadronic phase, and the MIT bag model, the Nambu-Jona-Lasinio, and the Chromo Dielectric models to describe the deconfined phase of quark matter. The consequences of the nucleation process for neutron star physics are outlined. Hyperonic stars are metastable only for some of the quark matter equations of state considered. The effect of an hyperonic three body force on the metastability of compact stars is estimated, and it is shown that, except for the Nambu-Jona-Lasinio model and the MIT bag model with a large bag pressure, the other models predict the formation of hybrid stars with a maximum mass not larger than \sim 1.62 M\odot .
The formation of non-relativistic collisionless shocks in laboratory with ultrahigh intensity lasers is studied via \emph{ab initio} multi-dimensional particle-in-cell simulations. The microphysics behind shock formation and dissipation, and the detailed shock structure are analyzed, illustrating that the Weibel instability plays a crucial role in the generation of strong subequipartition magnetic fields that isotropize the incoming flow and lead to the formation of a collisionless shock, similarly to what occurs in astrophysical scenarios. The possibility of generating such collisionless shocks in laboratory opens the way to the direct study of the physics associated with astrophysical shocks.
Links to: arXiv, form interface, find, astro-ph, recent, 1204, contact, help (Access key information)