The parametric resonance responsible for preheating after inflation will end when self-interactions of the resonating field and interactions of this field with secondary degrees of freedom become important. In many cases, the effect may be quantified in terms of an effective mass and the resulting shifting out of the spectrum of the strongest resonance band. In certain curvaton models, such thermal blocking can even occur before preheating has begun, delaying or even preventing the decay of the curvaton. We investigate numerically to what extent this thermal blocking is realised in a specific scenario, and whether the effective mass is well approximated by the perturbative leading order thermal mass. We find that the qualitative behaviour is well reproduced in this approximation, and that the end of preheating can be confidently estimated.
Poorly understood feedback processes associated with highly-luminous black hole accretion in quasars may dramatically affect the properties of their host galaxies. We search for the effect of quasar feedback on surrounding gas using Planck maps of the thermal Sunyaev-Zel'dovich effect (tSZ). By stacking tSZ Compton-y maps centered on the locations of 26,686 spectroscopic quasars from the Sloan Digital Sky Survey, we detect a strong but unresolved tSZ Compton-y signal at >5 sigma significance that likely originates from a combination of virialized halo atmosphere gas and quasar feedback effects. We show that the feedback contribution to our detected quasar tSZ signal is likely to dominate over virialized halo gas by isolating the feedback tSZ component for high- and low-redshift quasars. We find that this quasar tSZ signal also scales with black hole mass and bolometric luminosity, all consistent with general expectations of quasar feedback. We estimate the mean angularly-integrated Compton-y of quasars at z~1.5 to be 3.5x10^-6 Mpc^2, corresponding to mean total thermal energies in feedback and virialized halo gas of 1.1(+/- 0.2) x 10^62 erg, and discuss the implications for quasar feedback. If confirmed, the large total thermal feedback energetics we estimate of 5% (+/-1% statistical uncertainty) of the black hole mass will have important implications for the effects of quasar feedback on the host galaxy, as well as the surrounding intergalactic medium.
The solar corona is orders of magnitude hotter than the underlying photosphere, but how the corona attains such high temperatures is still not understood. Soft X-ray (SXR) emission provides important diagnostics for thermal processes in the high-temperature corona, and is also an important driver of ionospheric dynamics at Earth. There is a crucial observational gap between ~0.2 and ~4 keV, outside the ranges of existing spectrometers. We present observations from a new SXR spectrometer, the Amptek X123-SDD, which measured the spatially-integrated solar spectral irradiance from ~0.5 to ~5 keV, with ~0.15 keV FWHM resolution, during sounding rocket flights on 2012 June 23 and 2013 October 21. These measurements show that the highly variable SXR emission is orders of magnitude greater than that during the deep minimum of 2009, even with only weak activity. The observed spectra show significant high-temperature (5-10 MK) emission and are well fit by simple power-law temperature distributions with indices of ~6, close to the predictions of nanoflare models of coronal heating. Observations during the more active 2013 flight indicate an enrichment of low first-ionization potential (FIP) elements of only ~1.6, below the usually-observed value of ~4, suggesting that abundance variations may be related to coronal heating processes. The XUV Photometer System Level 4 data product, a spectral irradiance model derived from integrated broadband measurements, significantly overestimates the spectra from both flights, suggesting a need for revision of its non-flare reference spectra, with important implications for studies of Earth ionospheric dynamics driven by solar SXRs.
We measure the duty cycles for an existing sample of well observed, nearby dwarf novae using data from AAVSO, and present a quantitative empirical relation between the duty cycle of dwarf novae outbursts and the X-ray luminosity of the system in quiescence. We have found that $\log DC=0.63(\pm0.21)\times(\log L_{X}({\rm erg\,s^{-1}})-31.3)-0.95(\pm0.1)$, where DC stands for duty cycle. We note that there is intrinsic scatter in this relation greater than what is expected from purely statistical errors. Using the dwarf nova X-ray luminosity functions from \citet{Pretorius12} and \citet{Byckling10}, we compare this relation to the number of dwarf novae in the Galactic Bulge Survey which were identified through optical outbursts during an 8-day long monitoring campaign. We find a specific frequency of X-ray bright ($L_{X}>10^{31}\,{\rm erg\,s^{-1}}$) Cataclysmic Variables undergoing Dwarf Novae outbursts in the direction of the Galactic Bulge of $6.6\pm4.7\times10^{-5}\,M_{\odot}^{-1}$. Such a specific frequency would give a Solar neighborhood space density of long period CVs of $\rho=5.6\pm3.9\times10^{-6}\,$pc$^{-3}$. We advocate the use of specific frequency in future work, given that projects like LSST will detect DNe well outside the distance range over which $\rho\approx{\textrm const}$.
Multicolor photometric data are presented for the CALSPEC stars P177D and P330E. Together with previously published photometry for nine other CALSPEC standards, the photometric observations and synthetic photometry from HST/STIS spectrophotometry agree in the B, V, R, and I bands to better than $\sim$1\% (10 mmag).
Recent calculations using non-linear relativistic cosmological perturbation theory show biases in the mean luminosity distance and distance modulus at low redshift. We show that these effects may be understood very simply as a non-relativistic, and purely kinematic, Malmquist-like bias, and we describe how the effect changes if one averages over sources that are limited by apparent magnitude. This effect is essentially identical to the distance bias from small-scale random velocities that has previously been considered by astronomers, though we find that the standard formula overestimates the homogeneous bias by a factor 2.
We present the results from coordinated X-ray observations of the ultraluminous X-ray source NGC 5204 X-1 performed by NuSTAR and XMM-Newton in early 2013. These observations provide the first detection of NGC 5204 X-1 above 10 keV, extending the broadband coverage to 0.3-20 keV. The observations were carried out in two epochs separated by approximately 10 days, and showed little spectral variation, with an observed luminosity of Lx = (4.95+/-0.11)e39 erg/s. The broadband spectrum confirms the presence of a clear spectral downturn above 10 keV, only hinted at by previous observations. This cutoff is inconsistent with the standard low/hard state seen in Galactic black hole binaries, as would be expected from an intermediate mass black hole accreting at significantly sub-Eddington rates given the observed luminosity. The continuum is apparently dominated by two optically thick thermal-like components, potentially accompanied by a faint non-thermal tail at high energies. The broadband spectrum is likely associated with an accretion disk that differs from a standard Shakura & Sunyaev thin disk.
We have designed a realistic simulation of astronomical observing using a relatively low-cost commercial CCD camera and a microcontroller-based circuit that drives LEDs inside a light-tight box with time-varying intensities. As part of a laboratory experiment, students can acquire sequences of images using the camera, and then perform data analysis using a language such as MATLAB or Python to: (a) extract the intensity of the imaged LEDs, (b) perform basic calibrations on the time-series data, and (c) convert their data into the frequency domain where they can then identify the frequency structure. The primary focus is on studying light curves produced by the pulsating white dwarf stars. The exercise provides an introduction to CCD observing, a framework for teaching concepts in numerical data analysis and Fourier techniques, and connections with the physics of white dwarf stars.
The Fermi bubbles are gigantic gamma-ray structure in our Galaxy. The physical origin of the bubbles is still under debate. The leading scenarios can be divided into two categories. One is the nuclear star forming activity like extragalactic starburst galaxies and the other is the past active galactic nucleus (AGN) like activity of the Galactic center supermassive black hole. In this letter, we propose that metal abundance measurements will provide an important clue to probe their origin. Based on a simple spherically symmetric bubble model, we find that the generated metallicity and abundance pattern of the bubbles gas strongly depend on assumed star formation or AGN activities. Star formation scenarios predict higher metallicities and abundance ratios of [O/Fe] and [Ne/Fe] than AGN scenarios do because of supernovae ejecta. Furthermore, the resultant abundance depends on the gamma-ray emission process because different mass injection histories are required for different the gamma-ray emission processes due to the acceleration and cooling time scale of non-thermal particles. Future X-ray missions such as Astro-H and Athena will give a clue to probe the origin of the bubbles through the abundance measurements with their high energy resolution instruments.
NaSt1 (aka Wolf-Rayet 122) is a peculiar emission-line star embedded in an extended nebula of [N II] emission with a compact dusty core. This object was characterized by Crowther & Smith (1999) as a Wolf-Rayet (WR) star cloaked in an opaque nebula of CNO-processed material, perhaps analogous to Eta Car and its Homunculus nebula, albeit with a hotter central source. To discern the morphology of the [N II] nebula we performed narrowband imaging using the Hubble Space Telescope and Wide-field Camera 3. The images reveal that the nebula has a disk-like geometry tilted 12 degrees from edge-on, composed of a bright central ellipsoid surrounded by a larger clumpy ring. Ground-based spectroscopy reveals radial velocity structure (~10 km/s) near the outer portions of the nebula's major axis, which is likely to be the imprint of outflowing gas. Near-infrared adaptive-optics imaging with Magellan AO has resolved a compact ellipsoid of Ks-band emission aligned with the larger [N II] nebula, which we suspect is the result of scattered He I line emission (2.06 um). Observations with the Chandra X-ray Observatory have revealed an X-ray point source at the core of the nebula that is heavily absorbed at energies <1 keV and has properties consistent with WR stars and colliding-wind binaries. We suggest that NaSt1 is a WR binary embedded in an equatorial outflow that formed as the result of non-conservative mass transfer. NaSt1 thus appears to be a rare and important example of a stripped-envelope WR forming through binary interaction, caught in the brief Roche-Lobe overflow phase.
Cosmic rays (CR), constrained by scattering on magnetic irregularities, are
believed to propagate diffusively. But a well-known defect of diffusive
approximation, whereby not all the particles propagate at realistic speeds,
causes attempts to justify an alternative approach based on the "telegraph"
equation. However, its derivations often lack rigor and transparency leading to
incorrect results.
The classic Chapman-Enskog method is applied to the pitch-angle averaged
spatial CR transport. We show that the convective term arises only from the
magnetic focusing effect and no telegraph (second order time derivative) term
emerges in any order of the proper asymptotic expansion with systematically
eliminated short time scales. However, this term may formally be converted from
the fourth order hyper-diffusive term in the Chapman-Enskog expansion. But,
within the method's validity range, it may only be important for a short
relaxation period associated with either strong pitch-angle anisotropy or
spatial inhomogeneity of the initial CR distribution. However, in neither of
these two cases a treatment based on merely the isotropic CR component is
sufficient, regardless of the correction type (telegraph or hyper-diffusive).
Moreover, in a long-time asymptotic regime these corrections are generally
insignificant, as the angular anisotropy decays rapidly.
The main accretion phase of star formation is investigated in clouds with different metallicities in the range of 0 \le Z \le Z_\odot, resolving the protostellar radius. Starting from a near-equilibrium prestellar cloud, we calculate the cloud evolution up to \sim100 yr after the first protostar formation. The star formation process considerably differs between clouds with lower (Z \le 10^-4 Z_\odot) and higher (Z > 10^-4 Z_\odot) metallicities. Fragmentation frequently occurs and many protostars appear without forming a stable circumstellar disc in lower-metallicity clouds. In these clouds, although protostars mutually interact and some are ejected from the cloud centre, many remain as a small stellar cluster. In contrast, higher-metallicity clouds produce a single protostar surrounded by a nearly stable rotation-supported disc. In these clouds, although fragmentation occasionally occurs in the disc, the fragments migrate inwards and finally fall onto the central protostar. The difference in cloud evolution is due to different thermal evolutions and mass accretion rates. The thermal evolution of the cloud determines the emergence and lifetime of the first core. The first core develops prior to the protostar formation in higher-metallicity clouds, whereas no (obvious) first core appears in lower-metallicity clouds. The first core evolves into a circumstellar disc with a spiral pattern, which effectively transfers the angular momentum outwards and suppresses frequent fragmentation. In lower-metallicity clouds, the higher mass accretion rate increases the disc surface density within a very short time, rendering the disc unstable to self-gravity and inducing vigorous fragmentation.
We obtain a sample of 87 radio-loud QSOs in the redshift range 3.6<z<4.4 by cross-correlating sources in the FIRST radio survey S{1.4GHz} > 1 mJy with star-like objects having r <20.2 in SDSS Data Release 7. Of these 87 QSOs, 80 are spectroscopically classified in previous work (mainly SDSS), and form the training set for a search for additional such sources. We apply our selection to 2,916 FIRST-DR7 pairs and find 15 likely candidates. Seven of these are confirmed as high-redshift quasars, bringing the total to 87. The candidates were selected using a neural-network, which yields 97% completeness (fraction of actual high-z QSOs selected as such) and an efficiency (fraction of candidates which are high-z QSOs) in the range of 47 to 60%. We use this sample to estimate the binned optical luminosity function of radio-loud QSOs at $z\sim 4$, and also the LF of the total QSO population and its comoving density. Our results suggest that the radio-loud fraction (RLF) at high z is similar to that at low-z and that other authors may be underestimating the fraction at high-z. Finally, we determine the slope of the optical luminosity function and obtain results consistent with previous studies of radio-loud QSOs and of the whole population of QSOs. The evolution of the luminosity function with redshift was for many years interpreted as a flattening of the bright end slope, but has recently been re-interpreted as strong evolution of the break luminosity for high-z QSOs, and our results, for the radio-loud population, are consistent with this.
Imaging a star's companion at multiple epochs over a short orbital arc provides only four of the six coordinates required for a unique orbital solution. Probability distributions of possible solutions are commonly generated by Monte Carlo (MCMC) analysis, but these are biased by priors and may not probe the full parameter space. We suggest alternative methods to characterise possible orbits, which compliment the MCMC technique. Firstly the allowed ranges of orbital elements are prior-independent, and we provide means to calculate these ranges without numerical analyses. Hence several interesting constraints (including whether a companion even can be bound, its minimum possible semi-major axis and its minimum eccentricity) may be quickly computed using our relations as soon as orbital motion is detected. We also suggest an alternative to posterior probability distributions as a means to present possible orbital elements, namely contour plots of elements as functions of line of sight coordinates. These plots are prior-independent, readily show degeneracies between elements and allow readers to extract orbital solutions themselves. This approach is particularly useful when there are other constraints on the geometry, for example if a companion's orbit is assumed to be aligned with a disc. As examples we apply our methods to several imaged sub-stellar companions including Fomalhaut b, and for the latter object we show how different origin hypotheses affect its possible orbital solutions. We also examine visual companions of A- and G-type main sequence stars in the Washington Double Star Catalogue, and show that $\gtrsim50$ per cent must be unbound.
We present detailed chemical evolution models for the Milky Way and M31 in presence of radial gas flows. These models follow in detail the evolution of several chemical elements (H, He, CNO, $\alpha$ elements, Fe-peak elements) in space and time. The contribution of supernovae of different type to chemical enrichment is taken into account. We find that an inside-out formation of the disks coupled with radial gas inflows of variable speed can reproduce very well the observed abundance gradients in both galaxies. We also discuss the effects of other parameters, such as a threshold in the gas density for star formation and efficiency of star formation varying with galactic radius. Moreover, for the first time we compute the galactic habitable zone in our Galaxy and M31 in presence of radial gas flows. The main effect is to enhance the number of stars hosting a habitable planet with respect to the models without radial flow, in the region of maximum probability for this occurrence. In the Milky Way the maximum number of stars hosting habitable planets is at 8 kpc from the Galactic center, and the model with radial gas flows predicts a number of planets which is 38% larger than that predicted by the classical model.
Aims: To investigate the relationship between surges and magnetic reconnection during the emergence of small-scale active regions. In particular, to examine how the large-scale geometry of the magnetic field, shaped by different phases of reconnection, guides the flowing of surges. Methods: We present three flux emergence models. The first model, and the simplest, consists of a region emerging into a horizontal ambient field that is initially parallel to the top of the emerging region. The second model is the same as the first but with an extra smaller emerging region which perturbs the main region. This is added to create a more complex magnetic topology and to test how this complicates the development of surges compared to the first model. The last model has a non-uniform ambient magnetic field to model the effects of emergence near a sunspot field and impose asymmetry on the system through the ambient magnetic field. At each stage, we trace the magnetic topology to identify the locations of reconnection. This allows for field lines to be plotted from different topological regions, highlighting how their geometry affects the development of surges. Results: In the first model, we identify distinct phases of reconnection. Each phase is associated with a particular geometry for the magnetic field and this determines the paths of the surges. The second model follows a similar pattern to the first but with a more complex magnetic topology and extra eruptions. The third model highlights how an asymmetric ambient field can result in preferred locations for reconnection, subsequently guiding the direction of surges.
As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.
We report trace element analyses from mineral phases in chondrules from carbonaceous chondrites (Vigarano, Renazzo and Acfer 187), carried out by laser ablation inductively coupled plasma mass spectrometry. Results are similar in all three meteorites. Mesostasis has Rare Earth Element (REE) concentrations of 10-20 x CI. Low-Ca pyroxene has light REE (LREE) concentrations near 0.1 x CI and heavy REE (HREE) near 1 x CI respectively. Olivine has HREE concentrations at 0.1-1 x CI and LREE around 10-2 x CI. The coarsest olivine crystals tend to have the most fractionated REE patterns, indicative of equilibrium partitioning. Low-Ca pyroxene in the most pyroxene-rich chondrules tends to have the lowest REE concentrations. Type I chondrules seem to have undergone a significant degree of batch crystallization (as opposed to fractional crystallization), which requires cooling rates slower than 1-100 K/h. This would fill the gap between igneous CAIs and type II chondrules. The anticorrelation between REE abundances and pyroxene mode may be understood as due to dilution by addition of silica to the chondrule melt, as in the gas-melt interaction scenario of Libourel et al. (2006). The rapid cooling rate (of the order of 1000 K/h) which seems recorded by low-Ca pyroxene, contrasted with the more diverse record of olivine, may point to a nonlinear cooling history or suggest that formation of pyroxene-rich chondrule margins was an event distinct from the crystallization of the interior.
We report trace element analyses by laser ablation inductively coupled plasma mass spectrometry of metal grains from 9 different CR chondrites, distinguishing grains from chondrule interior ("interior grains"), chondrule surficial shells ("margin grains") and the matrix ("isolated grains"). Save for a few anomalous grains, Ni-normalized trace element patterns are similar for all three petrographical settings, with largely unfractionated refractory siderophile elements and depleted volatile Au, Cu, Ag, S. All types of grains are interpreted to derive from a common precursor approximated by the least melted, fine-grained objects in CR chondrites. This also excludes recondensation of metal vapor as the origin of the bulk of margin grains. The metal precursors presumably formed by incomplete condensation, with evidence for high-temperature isolation of refractory platinum-group-element (PGE)-rich condensates before mixing with lower temperature PGE-depleted condensates. The rounded shape of the Ni-rich, interior grains indicates melting and equilibration with silicates upon slow cooling (1-100 K/h), largely by oxidation/evaporation of Fe. We propose that Ni-poorer, amoeboid margin grains, often included in the pyroxene-rich periphery common to type I chondrules, result from less intense processing of a rim accreted onto the chondrule subsequent to the melting event recorded by the interior. This means either that there were two separate heating events, which formed olivine/interior grains and pyroxene/margin grains, respectively, between which dust was accreted around the chondrule, or there was a single high-temperature event, of which the chondrule margin record a late "quenching phase", in which case dust accreted onto chondrules while they were molten. In the latter case, high dust concentrations in the chondrule-forming region are indicated.
PolyChord is a novel nested sampling algorithm tailored for high dimensional parameter spaces. In addition, it can fully exploit a hierarchy of parameter speeds such as is found in CosmoMC and CAMB. It utilises slice sampling at each iteration to sample within the hard likelihood constraint of nested sampling. It can identify and evolve separate modes of a posterior semi-independently and is parallelised using openMPI. PolyChord is available for download at: this http URL
We present recent improvements of the modeling of the disruption of strength dominated bodies using the Smooth Particle Hydrodynamics (SPH) technique. The improvements include an updated strength model and a friction model, which are successfully tested by a comparison with laboratory experiments. In the modeling of catastrophic disruptions of asteroids, a comparison between old and new strength models shows no significant deviation in the case of targets which are initially non-porous, fully intact and have a homogeneous structure (such as the targets used in the study by Benz&Asphaug (1999). However, for many cases (e.g. initially partly or fully damaged targets, rubble-pile structures, etc.) we find that it is crucial that friction is taken into account and the material has a pressure dependent shear strength. Our investigations of the catastrophic disruption threshold $Q^*_{D}$ as a function of target properties and target sizes up to a few 100 km show that a fully damaged target modeled without friction has a $Q^*_{D}$ which is significantly (5-10 times) smaller than in the case where friction is included. When the effect of the energy dissipation due to compaction (pore crushing) is taken into account as well, the targets become even stronger ($Q^*_{D}$ is increased by a factor of 2-3). On the other hand, cohesion is found to have an negligible effect at large scales and is only important at scales $\lesssim$ 1km. Our results show the relative effects of strength, friction and porosity on the outcome of collisions among small ($\lesssim$ 1000 km) bodies. These results will be used in a future study to improve existing scaling laws for the outcome of collisions (e.g. Leinhardt&Stewart, 2012).
Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogs that include new deep near-infrared observations. Nonetheless we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1-b=0.76+-0.05(stat)}+-0.06(syst)}, which does not resolve the tension with the measurements from the primary cosmic microwave background.
A thermal relic, often referred to as a weakly interacting massive particle (WIMP),is a particle produced during the early evolution of the Universe whose relic abundance (e.g., at present) depends only on its mass and its thermally averaged annihilation cross section (annihilation rate factor) sigma*v_ann. Late time WIMP annihilation has the potential to affect the cosmic microwave background (CMB) power spectrum. Current observational constraints on the absence of such effects provide bounds on the mass and the annihilation cross section of relic particles that may, but need not be dark matter candidates. For a WIMP that is a dark matter candidate, the CMB constraint sets an upper bound to the annihilation cross section, leading to a lower bound to their mass that depends on whether or not the WIMP is its own antiparticle. For a self-conjugate WIMP, m_min = 50f GeV, where f is an electromagnetic energy efficiency factor. For a non self-conjugate WIMP, the minimum mass is a factor of two larger. For a WIMP that is a subdominant component of the dark matter density there is no bound on its mass and the upper bound to its annihilation cross section imposed by the CMB transforms into a lower bound to its annihilation cross section. These results are outlined and quantified here using the latest CMB constraints for a stable, symmetric (equal number of particles and antiparticles), WIMP whose annihilation is s-wave dominated, and for particles that are, or are not, their own antiparticle.
A viscous instability in shearing laminar axisymmetric hydrodynamic flows around a gravitating center is described. In the linearized hydrodynamic equations written in the Boussinesq approximation with microscopic molecular transport coefficients, the instability arises when the viscous dissipation is taken into account in the energy equation. Using the local WKB approximation, we derive a third-order algebraic dispersion equation with two modes representing the modified Rayleigh modes R+ and R-, and the third X-mode. We show that in thin accretion flows the viscosity destabilizes one of the Rayleigh modes in a wide range of wavenumbers, while the X-mode always remains stable. In Keplerian flows, the instability increment is found to be a few Keplerian rotational periods at wavelengths with $kr\sim 10-50$. This instability may cause turbulence in astrophysical accretion discs even in the absence of magnetic field.
Using images from the Hubble Space Telescope Wide-Field Camera 3, we measure the rate of diffusion of stars through the core of the globular cluster 47 Tucanae using a sample of young white dwarfs identified in these observations. This is the first direct measurement of diffusion due to gravitational relaxation. We find that the diffusion rate $\kappa\approx 10-13$ arcsecond$^2$ Myr$^{-1}$ is consistent with theoretical estimates of the relaxation time in the core of 47 Tucanae of about 70 Myr.
The ARGO-YBJ experiment is a full coverage air shower detector operated at the Yangbajing International Cosmic Ray Observatory. The detector has been in stable data taking in its full configuration since November 2007 to February 2013. The high altitude and the high segmentation and spacetime resolution offer the possibility to explore the cosmic ray energy spectrum in a very wide range, from a few TeV up to the PeV region. The high segmentation allows a detailed measurement of the lateral distribution, which can be used in order to discriminate showers produced by light and heavy elements. In this work we present the measurement of the cosmic ray light component spectrum in the energy range 3-3000 TeV. The analysis has been carried out by using a two-dimensional unfolding method based on the Bayes' theorem.
Despite ablation and drag processes associated with atmospheric entry of meteoroids were a subject of intensive study over the last century, little attention was devoted to interpret the observed fireball terminal height. This is a key parameter because it not only depends on the initial mass, but also on the bulk physical properties of the meteoroids and hence of their ability to ablate in the atmosphere. In this work we have developed a new approach that is tested using the fireball terminal heights observed by the Meteorite Observation and Recovery Project operated in Canada between 1970-1985 (hereafter referred as MORP). We then compare them to the calculation made. Our results clearly show that the new methodology is able to forecast the degree of deepening of meteoroids in the Earth's atmosphere. Then, this approach has important applications in predicting the impact hazard from cm- to meter-sized bodies that are represented, in part, in the MORP bolide list.
Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.
Observational data from the Cassini spacecraft are used to obtain a chemical model of ocean water on Enceladus. The model indicates that Enceladus' ocean is a Na-Cl-CO3 solution with an alkaline pH of ~11-12. The dominance of aqueous NaCl is a feature that Enceladus' ocean shares with terrestrial seawater, but the ubiquity of dissolved Na2CO3 suggests that soda lakes are more analogous to the Enceladus ocean. The high pH implies that the hydroxide ion should be relatively abundant, while divalent metals should be present at low concentrations owing to buffering by clays and carbonates on the ocean floor. The high pH is interpreted to be a key consequence of serpentinization of chondritic rock, as predicted by prior geochemical reaction path models; although degassing of CO2 from the ocean may also play a role depending on the efficiency of mixing processes in the ocean. Serpentinization leads to the generation of H2, a geochemical fuel that can support both abiotic and biological synthesis of organic molecules such as those that have been detected in Enceladus' plume. Serpentinization and H2 generation should have occurred on Enceladus, like on the parent bodies of aqueously altered meteorites; but it is unknown whether these critical processes are still taking place, or if Enceladus' rocky core has been completely altered by past hydrothermal activity. The high pH also suggests that the delivery of oxidants from the surface to the ocean has not been significant, and the rocky core did not experience partial melting and igneous differentiation. On the other hand, the pH is compatible with life as we know it; life on Earth may have begun under similar conditions, and serpentinites on Earth support microbial communities that are centered on H2 that is provided by water-rock reactions.
In the present paper, we compare the predictions of two well known mechanisms considered able to solve the cusp/core problem (a. supernova feedback; b. baryonic clumps-DM interaction) by comparing their theoretical predictions to recent observations of the inner slopes of galaxies with masses ranging from dSphs to normal spirals. We compare the $\alpha$-$V_{\rm rot}$ and the $\alpha$-$M_{\ast}$ relationships, predicted by the two models with high resolution data coming from \cite{Adams2014}, \cite{Simon2005}, LITTLE THINGS \citep{Oh2014}, THINGS dwarves \citep{Oh2011a,Oh2011b}, THINGS spirals \citep{Oh2014}, Sculptor, Fornax and the Milky Way. The comparison of the theoretical predictions with the complete set of data shows that the two models perform similarly, while when we restrict the analysis to a smaller subsample of higher quality, we show that the method presented in this paper (baryonic clumps-DM interaction) performs better than the one based on supernova feedback. We also show that, contrarily to the first model prediction, dSphs of small mass could have cored profiles. This means that observations of cored inner profiles in dSphs having a stellar mass $<10^6 M_{\odot}$ not necessarily imply problems for the $\Lambda$CDM model.
Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than 0.5-0.7Msun. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.
We present results from a fifteen-month campaign of high-cadence (~ 3 days) mid-infrared Spitzer and optical (B and V ) monitoring of the Seyfert 1 galaxy NGC 6418, with the objective of determining the characteristic size of the dusty torus in this active galactic nucleus (AGN). We find that the 3.6 $\mu$m and 4.5 $\mu$m flux variations lag behind those of the optical continuum by $37.2^{+2.4}_{-2.2}$ days and $47.1^{+3.1}_{-3.1}$ days, respectively. We report a cross-correlation time lag between the 4.5 $\mu$m and 3.6 $\mu$m flux of $13.9^{+0.5}_{-0.1}$ days. The lags indicate that the dust emitting at 3.6 $\mu$m and 4.5 $\mu$m is located at a distance of approximately 1 light-month (~ 0.03 pc) from the source of the AGN UV-optical continuum. The reverberation radii are consistent with the inferred lower limit to the sublimation radius for pure graphite grains at 1800 K, but smaller by a factor of ~ 2 than the corresponding lower limit for silicate grains; this is similar to what has been found for near-infrared (K-band) lags in other AGN. The 3.6 and 4.5 $\mu$m reverberation radii fall above the K-band $\tau \propto L^{0.5}$ size-luminosity relationship by factors $\lesssim 2.7$ and $\lesssim 3.4$, respectively, while the 4.5 $\mu$m reverberation radius is only 27% larger than the 3.6 $\mu$m radius. This is broadly consistent with clumpy torus models, in which individual optically thick clouds emit strongly over a broad wavelength range.
We update the forecasts for the measurement of the tensor-to-scalar ratio $r$ for various ground-based experiments (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloons (EBEX and Spider) and satellites (CMBPol, COrE and LiteBIRD), taking into account the recent Planck data on polarized dust and using a component separation method. The forecasts do not change significantly with respect to previous estimates when at least three frequencies are available, provided foregrounds can be accurately described by few parameters. We argue that a theoretically motivated goal for future experiments is $r\sim2\times10^{-3}$, and that this is achievable if the noise is reduced to $\sim1\,\mu$K-arcmin and lensing is reduced to $10\%$ in power. We study the constraints experiments will be able to put on the frequency and $\ell$-dependence of the tensor signal as a check of its primordial origin. Futuristic ground-based experiments can have good constraints on these parameters, even for $r\sim2\times10^{-3}$. For the same value of $r$, satellites will marginally be able to detect the presence of the recombination bump, the most distinctive feature of the primordial signal.
Chemically peculiar stars of the upper part of the main sequence show periodical variability in line intensities and continua, modulated by the stellar rotation, which is attributed to the existence of chemical spots on the surface of these stars. The flux variability is caused by the changing redistribution rate of the radiative flux predominantly from the short-wavelength part of the spectra to the long-wavelength part, which is a result of abundance anomalies. We study the nature of the multi-spectral variability of one of the brightest chemically peculiar stars, $\theta$ Aur. We predict the flux variability of $\theta$ Aur from the emerging intensities calculated for individual surface elements of the star taking into account horizontal variation of chemical composition derived from Doppler abundance maps. The simulated optical variability in the Str\"omgren photometric system and the ultraviolet flux variability agree well with observations. The IUE flux distribution is reproduced in great detail by our models. The resonance lines of magnesium and possibly also some lines of silicon are relatively weak in the ultraviolet domain, which indicates non-negligible vertical abundance gradients in the atmosphere. We also derive a new period of the star, $P=3.618\,664(10)$ d, from all available photometric and magnetic measurements and show that the observed rotational period is constant over decades. The ultraviolet and visual variability of $\theta$ Aur is mostly caused by silicon bound-free absorption and chromium and iron line absorption. These elements redistribute the flux mainly from the far-ultraviolet region to the near-ultraviolet and optical regions in the surface abundance spots. The light variability is modulated by the stellar rotation. The ultraviolet domain is key for understanding the properties of chemically peculiar stars. (abridged)
We report observations of four sub-damped Lyman-alpha (sub-DLA) quasar absorbers at z<0.5 obtained with the Hubble Space Telescope Cosmic Origins Spectrograph. We measure the available neutrals or ions of C, N, O, Si, P, S, Ar, Mn, Fe, and/or Ni. Our data have doubled the sub-DLA metallicity samples at z<0.5 and improved constraints on sub-DLA chemical evolution. All four of our sub-DLAs are consistent with near-solar or super-solar metallicities and relatively modest ionization corrections; observations of more lines and detailed modeling will help to verify this. Combining our data with measurements from the literature, we confirm previous suggestions that the N(HI)-weighted mean metallicity of sub-DLAs exceeds that of DLAs at all redshifts studied, even after making ionization corrections for sub-DLAs. The absorber toward PHL 1598 shows significant dust depletion. The absorbers toward PHL 1226 and PKS 0439-433 show the S/P ratio consistent with solar, i.e., they lack a profound odd-even effect. The absorber toward Q0439-433 shows super-solar Mn/Fe. For several sub-DLAs at z<0.5, [N/S] is below the level expected for secondary N production, suggesting a delay in the release of the secondary N or a tertiary N production mechanism. We constrain the electron density using Si II* and C II* absorption. We also report different metallicity vs. Delta V_90 relations for sub-DLAs and DLAs. For two sub-DLAs with detections of emission lines from the underlying galaxies, our measurements of the absorption-line metallicities are consistent with the emission-line metallicities, suggesting that metallicity gradients are not significant in these galaxies.
If dark matter is unstable and the mass is within GeV-TeV regime, its decays produce high-energy photons that give contribution to the extragalactic gamma-ray background (EGRB). We constrain dark matter decay by analyzing the 50-month EGRB data measured with Fermi satellite, for different decay channels motivated with several supersymmetric scenarios featuring R-parity violation. We adopt the latest astrophysical models for various source classes such as active galactic nuclei and star-forming galaxies, and take associated uncertainties properly into account. The lower limits for the lifetime are very stringent for a wide range of dark matter mass, excluding the lifetime shorter than 10^28 s for mass between a few hundred GeV and ~1TeV, e.g., for b\bar{b} decay channel. Furthermore, most dark matter models that explain the anomalous positron excess are also excluded. These constraints are robust, being little dependent on astrophysical uncertainties, unlike other probes such as Galactic positrons or anti-protons.
As part of a program to identify the physical conditions in the jets of gamma-ray-flaring blazars detected by Fermi, including the role of shocks in the production of high-energy flaring, we obtained 4 years of 3-frequency, centimeter-band total flux density and linear polarization monitoring observations of the radio-bright blazar S5 0716+714 with the University of Michigan 26-m paraboloid. Light curves constructed from these data exhibit a series of rapid, high-amplitude, centimeter-band total flux density outbursts, and changes in the linear polarization consistent with the passage of shocks during the gamma-ray flaring. The observed spectral evolution of the radio-band flares, in combination with radiative transfer simulations incorporating propagating shocks, was used to constrain the shock and jet flow conditions in the parsec-scale regions of the jet. Eight forward-moving, transverse shocks with unusually-strong shock compression factors, a very fast Lorentz factor of the shocks of 77, a bulk Lorentz factor of the flow of 20, a viewing angle of 12 degrees, and an intrinsic opening angle of the radio jet of 5.2 degrees were identified.
We discuss new constraints on the epoch of cosmic reionization and test the assumption that most of the ionizing photons responsible arose from high redshift star-forming galaxies. Good progress has been made in charting the end of reionization through spectroscopic studies of z~6-8 QSOs, gamma-ray bursts and galaxies expected to host Lyman-alpha emission. However, the most stringent constraints on its duration have come from the integrated optical depth, tau, of Thomson scattering to the cosmic microwave background. Using the latest data on the abundance and luminosity distribution of distant galaxies from Hubble Space Telescope imaging, we simultaneously match the reduced value tau=0.066 +/- 0.012 recently reported by the Planck collaboration and the evolving neutrality of the intergalactic medium with a reionization history within 6 <~ z <~ 10, thereby reducing the requirement for a significant population of very high redshift (z>>10) galaxies. Our analysis strengthens the conclusion that star-forming galaxies dominated the reionization process and has important implications for upcoming 21cm experiments and searches for early galaxies with James Webb Space Telescope.
We re-evaluate the prompt atmospheric neutrino flux, using the measured charm cross sections at RHIC and the Large Hadron Collider to constrain perturbative QCD parameters such as the factorization and renormalization scales, as well as modern parton distribution functions and recent estimates of the cosmic-ray spectra. We find that our result for the prompt neutrino flux is lower than previous perturbative QCD estimates and, consequently, alters the signal-to-background statistics of the recent IceCube measurements at high energies.
It is known that in single scalar field inflationary models the standard curvature perturbation \zeta, which is supposedly conserved at superhorizon scales, diverges during reheating at times d\Phi/dt=0, i.e. when the time derivative of the background inflaton field vanishes. This happens because the comoving gauge \phi=0, where \phi\ denotes the inflaton perturbation, breaks down when d\Phi/dt=0. The issue is usually bypassed by averaging out the inflaton oscillations but strictly speaking the evolution of \zeta\ is ill posed mathematically. We solve this problem by introducing a family of smooth gauges that still eliminates the inflaton fluctuation \phi\ in the Hamiltonian formalism and gives a well behaved curvature perturbation \zeta, which is now rigorously conserved at superhorizon scales. In the linearized theory, this conserved variable can be used to unambiguously propagate the inflationary perturbations from the end of inflation to subsequent epochs. We discuss the implications of our results for the inflationary predictions.
We consider indirect detection signals of atomic dark matter, with a massive dark photon which mixes kinetically with hypercharge. In significant regions of parameter space, dark matter remains at least partially ionized today, and dark atom formation can occur efficiently in dense regions, such as the centers of galactic halos. The formation of dark atoms is accompanied by emission of a dark photon, which can subsequently decay into Standard Model particles. We discuss the expected signal strength and compare it to that of annihilating dark matter. As a case study, we explore the possibility that dark atom formation can account for the observed 511 keV line and outline the relevant parameter space.
Void induced coronene C23H12++ was suggested to be a possible carrier of the astronomically observed polycyclic aromatic hydrocarbon (PAH), which shows unique molecular structure with carbon two pentagons connected with five hexagons. Well observed astronomical infrared spectrum from 3-15 micron could be almost reproduced based on density functional theory. However, there remain several discrepancies with observed spectra, especially on 11-15 micron band weaker intensity. Observed 11.2 micron intensity is comparable to 7.6-7.8 micron one. Methyl-modified molecule C24H14++ revealed that calculated peak height of 11.4 micron show fairly large intensity up to 70-90% compared with that of 7.6-7.8 micron band. Also, nitrogen atom was substituted to peripheral C-H site of void coronene to be C22H11N1++. Pentagon site substituted case show 60% peak height. This molecule also reproduced well 12-15 micron peak position and relative intensity. Vibration mode analysis demonstrated that 11.3 micron mode comes from C-H out of plane bending. Heavy nitrogen plays as like an anchor role for molecule vibration.
In a model of 3-brane with matter embedded in 5D space-time conditions for gravitons emitted from the brane to the bulk to return back to the brane are found. For a given 5-momentum of graviton falling back to the brane the interval between emission and detection times is calculated. A method to calculate contribution to the energy-momentum tensor from multiple graviton bouncings is developed. Explicit expressions for the energy-momentum tensors of gravitons which have made one, two and three bounces are obtained and their magnitudes are numerically calculated. These expressions are used to solve the evolution equation for dark radiation. A relation connecting reheating temperature and the scale of extra dimension is obtained. For the reheating temperature $T_R\sim 10^6 GeV$ we estimate the scale of extra dimension to be of order $10^{-9}\div 10^{-10} GeV$.
We propose a large class of nonsingular cosmologies of arbitrary spatial curvature whose cosmic history is determined by a primeval dynamical $\Lambda (t)$-term. For all values of the curvature, the models evolve between two extreme de Sitter phases driven by the relic time-varying vacuum energy density. The transition from inflation to the radiation phase is universal and points to a natural solution of the graceful exit problem regardless of the values of the curvature parameter. The flat case recovers the scenario recently discussed in the literature (Perico et al., Phys. Rev. D88, 063531, 2013). The early de Sitter phase is characterized by an arbitrary energy scale $H_I$ associated to the primeval vacuum energy density. If $H_I$ is fixed to be nearly the Planck scale, the ratio between the relic and the present observed vacuum energy density is $\rho_{vI}/\rho_{v0} \simeq 10^{123}$.
Observational data have put unprecedentedly tight constraints on the behavior of the universe. Majority of the observational data still support the concordance $\Lambda$CDM model as a preferred explanation for the late time acceleration, but it may also allows a deviation from the corresponding cosmological constant $\Lambda$. As we know, the cosmological constant corresponds to the pressure $P_{\Lambda}=P_0$ (in which $P_0$ is a constant), ignoring the pressure from radiation, then the total pressure $P=P_{\Lambda}=P_0$. In this paper, we propose two parametric models for the total pressure $P(z)=P_a+P_b z$ and $P(z)=P_c+\frac{P_d}{1+z}$ to study the universe evolution at low redshift. We mimic our phenomenological models in the scenario of two type of scalar fields---quintessence and phantom. We constrain our model parameters with Supernova type Ia dataset and BAO dataset, and we show that data fitting results of the two models both mildly support $\omega_{de}<-1$ which implies a phantom dark energy scenario at the present time.
Fermi Normal Coordinates (FNC) are a useful frame for isolating the locally observable, physical effects of a long-wavelength spacetime perturbation. Their cosmological application, however, is hampered by the fact that they are only valid on scales much smaller than the horizon. We introduce a generalization that we call Conformal Fermi Coordinates (CFC). CFC preserve all the advantages of FNC, but in addition are valid outside the horizon. They allow us to calculate the coupling of long- and short-wavelength modes on all scales larger than the sound horizon of the cosmological fluid, starting from the epoch of inflation until today, by removing the complications of the second order Einstein equations to a large extent, and eliminating all gauge ambiguities. As an application, we present a calculation of the effect of long-wavelength tensor modes on small scale density fluctuations. We recover previous results, but clarify the physical content of the individual contributions in terms of locally measurable effects and "projection" terms.
We have studied the dynamics of an equal-mass magnetized neutron-star binary within a resistive magnetohydrodynamic (RMHD) approach in which the highly conducting stellar interior is matched to an electrovacuum exterior. Because our analysis is aimed at assessing the modifications introduced by resistive effects on the dynamics of the binary after the merger and through to collapse, we have carried out a close comparison with an equivalent simulation performed within the traditional ideal-MHD (IMHD) approximation. We have found that there are many similarities between the two evolutions, but also one important difference: the survival time of the hypermassive neutron star increases in a RMHD simulation. This difference is due to a less efficient magnetic-braking mechanism in the resistive regime, in which matter can move across magnetic-field lines, thus reducing the outward transport of angular momentum. Interestingly, a longer-lived magnetized hypermassive neutron star brings support to the recent modelling of short gamma-ray bursts in terms of such objects. Both the RMHD and the IMHD simulations carried here have been performed at higher resolutions and with a different grid structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L. Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732, L6 (2011)], but confirm the formation of a magnetic-jet structure in the low-density funnel produced by the black-hole--torus system. In both regimes the magnetic field is predominantly toroidal in the highly conducting torus and predominantly poloidal in the nearly evacuated funnel. Reconnection processes or neutrino annihilation occurring in the funnel, none of which we model here, could potentially increase the internal energy in the funnel and launch a relativistic outflow.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
Reverberation mapping probes the structure of the broad emission-line region (BLR) in active galactic nuclei (AGN). The radius of the BLR along with the virial velocity of the BLR gas can be used to measure the mass of the supermassive black hole. The main systematic uncertainty affecting reverberation mapping is the unknown structure of the BLR. We develop a new method for analysing reverberation mapping data based on regularized linear inversion (RLI) that includes statistical modelling of the AGN continuum light curves. This method enables fast, flexible, and robust calculation of velocity-resolved response maps to probe BLR structure. Contrary to other methods, RLI allows for negative response in the BLR, such as when some areas of the BLR respond in inverse proportion to a change in ionizing continuum luminosity. We present time delays, integrated response functions, and velocity-delay maps for the H{\beta} broad emission line in five nearby AGN, as well as H{\alpha} and H{\gamma} broad emission lines in Arp 151, using data from the Lick AGN Monitoring Project 2008. We find indications of prompt response in three of the objects (Arp 151, NGC 5548 and SBS 1116+583A) with additional prompt response in the red wing of H{\beta}. In SBS 1116+583A we find evidence for a multimodal broad prompt response followed by a second narrow response at 10 days. All time delays are consistent with results from cross correlation and the velocity-delay maps are consistent with results from maximum entropy and dynamical modelling methods. There is no clear indication of negative response in any of the objects. We conclude that regularized linear inversion with statistical light curve modelling provides a fast, complementary method for velocity-resolved reverberation mapping that makes very few assumptions regarding the shape of the response function and requires minimal user input.
We present a spectroscopic study of 287 Planetary Nebulas (PNs) in a total area of ~0.4 deg^2 around the BCG M87 in Virgo A. With these data we can distinguish the stellar halo from the co-spatial intracluster light (ICL). PNs were identified from their narrow and symmetric redshifted lambda 5007\4959 Angstrom [OIII] emission lines, and the absence of significant continuum. We implement a robust technique to measure the halo velocity dispersion from the projected phase-space to identify PNs associated with the M87 halo and ICL. The velocity distribution of the spectroscopically confirmed PNs is bimodal, containing a narrow component centred on the systemic velocity of the BCG and an off-centred broader component, that we identify as halo and ICL, respectively. Halo and ICPN have different spatial distributions: the halo PNs follow the galaxy's light, whereas the ICPNs are characterised by a shallower power-law profile. The composite PN number density profile shows the superposition of different PN populations associated with the M87 halo and the ICL, characterised by different PN alpha-parameters, the ICL contributing ~3 times more PNs per unit light. Down to m_5007=28.8, the M87 halo PN luminosity function (PNLF) has a steeper slope towards faint magnitudes than the IC PNLF, and both are steeper than the standard PNLF for the M31 bulge. Moreover, the IC PNLF has a dip at ~1-1.5 mag fainter than the bright cutoff, reminiscent of the PNLFs of systems with extended star formation history. The M87 halo and the Virgo ICL are dynamically distinct components with different density profiles and velocity distribution. The different alpha values and PNLF shapes of the halo and ICL indicate distinct parent stellar populations, consistent with the existence of a gradient towards bluer colours at large radii. These results reflect the hierarchical build-up of the Virgo cluster.
We report the discovery of the resolved disk around HD 131835 and present the analysis and modeling of its thermal emission. HD 131835 is a ~15 Myr A2 star in the Scorpius-Centaurus OB association at a distance of 122.7 +16.2 -12.8 parsec. The extended disk has been detected to ~1.5" (200 AU) at 11.7 {\mu}m and 18.3 {\mu}m with T-ReCS on Gemini South. The disk is inclined at an angle of ~75{\deg} with the position angle of ~61{\deg}. The flux of HD 131835 system is 49.3+-7.6 mJy and 84+-45 mJy at 11.7 {\mu}m and 18.3 {\mu}m respectively. A model with three grain populations gives a satisfactory fit to both the spectral energy distribution and the images simultaneously. This best-fit model is composed of a hot continuous power-law disk and two rings. We characterized the grain temperature profile and found that the grains in all three populations are emitting at temperatures higher than blackbodies. In particular, the grains in the continuous disk are unusually warm; even when considering small graphite particles as the composition.
We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. For the first time we demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^6M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.
\We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive.
Solitary stars that wander too close to their galactic centres can become tidally disrupted, if the tidal forces due to the supermassive black hole (SMBH) residing there overcome the self-gravity of the star. If the star is only partially disrupted, so that a fraction survives as a self-bound object, this remaining core will experience a net gain in specific orbital energy, which translates into a velocity "kick" of up to $\sim 10^3$ km/s. In this paper, we present the result of smoothed particle hydrodynamics (SPH) simulations of such partial disruptions, and analyse the velocity kick imparted on the surviving core. We compare $\gamma$ = 5/3 and $\gamma$ = 4/3 polytropes disrupted in both a Newtonian potential, and a generalized potential that reproduces most relativistic effects around a Schwarzschild black hole either exactly or to excellent precision. For the Newtonian case, we confirm the results of previous studies that the kick velocity of the surviving core is virtually independent of the ratio of the black hole to stellar mass, and is a function of the impact parameter $\beta$ alone, reaching at most the escape velocity of the original star. For a given $\beta$, relativistic effects become increasingly important for larger black hole masses. In particular, we find that the kick velocity increases with the black hole mass, making larger kicks more common than in the Newtonian case, as low-$\beta$ encounters are statistically more likely than high-$\beta$ encounters. The analysis of the tidal tensor for the generalized potential shows that our results are robust lower limits on the true relativistic kick velocities, and are generally in very good agreement with the exact results.
Future direct observations of extrasolar Earth-sized planets in the habitable zone could be hampered by a worrisome source of noise, starlight-reflecting exozodiacal dust. Mid-infrared surveys are currently underway to constrain the amount of exozodiacal dust in the habitable zones around nearby stars. However, at visible wavelengths another source of dust, invisible to these surveys, may dominate over exozodiacal dust. For systems observed near edge-on, a cloud of dust with face-on optical depth 10^-7 beyond ~5 AU can mimic the surface brightness of a cloud of exozodiacal dust with equal optical depth if the dust grains are sufficiently forward-scattering. We posit that dust migrating inward from cold debris belts via Poynting-Robertson drag could produce this "pseudo-zodiacal" effect, potentially making it ~50% as common as exozodiacal clouds. We place constraints on the disk radii and scattering phase function required to produce the effect.
With a substantial nuclear molecular gas reservoir and broad, high-velocity CO molecular line wings previously interpreted as an outflow, NGC 1266 is a rare SB$0$ galaxy. Previous analyses of interferometry, spectrally resolved low-$J$ CO emission lines, and unresolved high-$J$ emission lines have established basic properties of the molecular gas and the likely presence of an AGN. Here, new spectrally resolved CO $J = 5 - 4$ to $J = 8 - 7$ lines from {\it Herschel Space Observatory} HIFI observations are combined with ground-based observations and high-$J$ {\it Herschel} SPIRE observations to decompose the nuclear and putative outflow velocity components and to model the molecular gas to quantify its properties. Details of the modeling and results are described, with comparisons to previous results and exploration of the implications for the gas excitation mechanisms. Among the findings, like for other galaxies, the nuclear and putative outflow molecular gas are well represented by components that are cool ($T_{nuclear} = 6^{+10}_{-2}$ K and $T_{outflow} \sim 30$ K), comprising bulk of the mass (Log $M_{nuclear}/M_{\odot} = 8.3^{+0.5}_{-0.4}$ and Log $M_{outflow}/M_{\odot} = 7.6^{+0.3}_{-0.3}$), and the minority of the luminosity (Log $L_{nuclear}/L_{\odot} = 5.44^{+0.22}_{-0.18}$ and Log $L_{outflow}/L_{\odot} \sim 6.5$) and warm ($T_{nuclear} = 74^{+130}_{-26}$ K and $T_{outflow} > 100$ K), comprising a minority of the mass (Log $M_{nuclear}/M_{\odot} = 7.3^{+0.5}_{-0.5}$ and Log $M_{outflow}/M_{\odot} \sim 6.3$) but the majority of the luminosity (Log $L_{nuclear}/L_{\odot} = 6.90^{+0.16}_{-0.16}$ and Log $L_{outflow}/L_{\odot} \sim 7.2$). The outflow has an anomalously high $L_\mathrm{CO}/L_\mathrm{FIR}$ of $1.7 \times 10^{-3}$ and is almost certainly shock excited.
Context. Gravitational collapse theory and numerical simulations suggest that
the velocity field within large-scale galaxy filaments is dominated by motions
along the filaments.
Aims. Our aim is to check whether observational data reveal any preferred
orientation of galaxy pairs with respect to the underlying filaments as a
result of the expectedly anisotropic velocity field.
Methods. We use galaxy pairs and galaxy filaments identified from the Sloan
Digital Sky Survey data. For filament extraction, we use the Bisous model that
is based the marked point process technique. During the filament detection, we
use the centre point of each pair instead of the positions of galaxies to avoid
a built-in influence of pair orientation on the filament construction. For
pairs lying within filaments (3012 cases), we calculate the angle between the
line connecting galaxies of each pair and their host filament. To avoid
redshift-space distortions, the angle is measured in the plain of the sky.
Results. The alignment analysis shows that the orientation of galaxy pairs
correlates strongly with their host filaments. The alignment signal is stronger
for loose pairs, with at least 25% excess of aligned pairs compared to a random
distribution.
Conclusions. We conclude that the velocity field within large scale filaments
is indeed dominated by motions along the filaments (gravitational collapses,
streaming motions), bringing along the alignment of pair member trajectories
with galaxy filaments.
The accretion of satellites onto central galaxies along vast cosmic filaments is an apparent outcome of the anisotropic collapse of structure in our Universe. Numerical work (based on gravitational dynamics of N-body simulations) indicates that satellites are beamed towards hosts along preferred directions imprinted by the velocity shear field. Here we use the Sloan Digital Sky Survey to observationally test this claim. We construct 3D filaments and sheets and examine the relative position of satellites galaxies. A statistically significant alignment between satellite galaxy position and filament axis is confirmed. We find a similar (but stronger) signal by examining satellites and filaments similarly identified in the Millennium simulation, semi-analytical galaxy catalogue. We also examine the dependence of the alignment strength on galaxy properties such as colour, magnitude and (relative) satellite magnitude, finding that the alignment is strongest for the reddest and brightest central and satellite galaxies. Our results confirm the theoretical picture and the role of the cosmic web in satellite accretion. Furthermore our results suggest that filaments identified on larger scales can be reflected in the positions of satellite galaxies that are quite close to their hosts.
We have used the Robert C. Byrd Green Bank Telescope to perform the most sensitive search to date for neutral atomic hydrogen (HI) in the circumstellar envelope (CSE) of the carbon star IRC+10216. Our observations have uncovered a low surface brightness HI shell of diameter ~1300" (~0.8 pc), centered on IRC+10216. The HI shell has an angular extent comparable to the far ultraviolet-emitting astrosphere of IRC+10216 previously detected with the GALEX satellite, and its kinematics are consistent with circumstellar matter that has been decelerated by the local interstellar medium. The shell appears to completely surround the star, but the highest HI column densities are measured along the leading edge of the shell, near the location of a previously identified bow shock. We estimate a total mass of atomic hydrogen associated with IRC+10216 CSE of M_HI~3x10e-3 M_sun. This is only a small fraction of the expected total mass of the CSE (<1%) and is consistent with the bulk of the stellar wind originating in molecular rather than atomic form, as expected for a cool star with an effective temperature T_eff<~2200 K. HI mapping of a 2 deg x 2 deg region surrounding IRC+10216 has also allowed us to characterize the line-of-sight interstellar emission in the region and has uncovered a link between diffuse FUV emission southwest of IRC+10216 and the Local Leo Cold Cloud.
The linear growth rate is commonly defined through a simple deterministic relation between the velocity divergence and the matter overdensity in the linear regime. Here we introduce a formalism that extends this to a nonlinear, stochastic relation between $\theta = \nabla \cdot v({\bf x},t)/aH$ and $\delta$. This provides a new phenomenological approach that examines the conditional mean $\langle \theta|\delta\rangle$, together with fluctuations of $\theta$ around this mean. We measure these stochastic components of the velocity power spectra using simulations and find they are non-negative and increase with decreasing scale from $\sim$10% at $k<0.2 h $Mpc$^{-1}$ to 25% at $k\sim0.45h$Mpc$^{-1}$. We find that both the stochastic relation and nonlinearity are more pronounced for halos of mass $M \le 5 \times 10^{12}M_\odot h^{-1}$ compared to the dark matter at $z=0$ and $1$. Nonlinear growth effects manifest themselves in the formalism as a rotation of the mean $\langle \theta|\delta\rangle$ away from the linear theory prediction $-f_{\tiny \rm LT}\delta$, where $f_{\tiny \rm LT}$ is the linear theory growth rate. This rotation increases with increasing $k$ and can be well-described by 2nd order Lagrangian perturbation theory (2LPT) for $k < 0.1 h$Mpc$^{-1}$. The stochasticity in the $\theta$ - $\delta$ relation is not so simply described by 2LPT, and we discuss its impact on measurements of $f_{\tiny \rm LT}$ from clustering statistics in redshift space. In the presence of either nonlinearity or a stochastic relation, the correspondence between the $\mu^2$ and $\mu^4$ coefficients, and $f_{\rm LT}$ demands a more complex treatment. Given that the relationship between $\delta$ and $\theta$ is stochastic and nonlinear, this will have implications for the interpretation and precision of the linear growth rate extracted using models which assume a linear, deterministic expression.
The process by which the Sun affects the terrestrial environment on short timescales is predominately driven by the amount of magnetic reconnection between the solar wind and Earth's magnetosphere. Reconnection occurs most efficiently when the solar wind magnetic field has a southward component. The most severe impacts are during the arrival of a coronal mass ejection (CME) when the magnetosphere is both compressed and magnetically connected to the heliospheric environment, leading to disruptions to, for example, power grids and satellite navigation. Unfortunately, forecasting magnetic vectors within coronal mass ejections remains elusive. Here we report how, by combining a statistically robust helicity rule for a CME's solar origin with a simplified flux rope topology the magnetic vectors within the Earth-directed segment of a CME can be predicted. In order to test the validity of this proof-of-concept architecture for estimating the magnetic vectors within CMEs, a total of eight CME events (between 2010 and 2014) have been investigated. The angular rotation in the predicted magnetic field closely follows the broad rotational structure seen within the in situ data. This time varying field estimate is implemented into a process to quantitatively predict a time-varying Kp index. It is expected that future statistical work to better quantify the uncertainties may help the heuristic approach of the early forecasting systems used by forecasters.
We report follow-up spectroscopy of 29 cataclysmic variables from the Sloan Digital Sky Survey (SDSS), 22 of which were discovered by SDSS and seven other previously known systems that were recovered in SDSS. The periods for 16 of these objects were included in the tabulation by Gaensicke et al. (2009). While most of the systems have periods less than 2 hours, only one has a period in the 80-86 minute 'spike' found by Gaensicke et al. (2009), and 11 have periods longer than 3 hours, indicating that the present sample is skewed toward longer-period, higher-luminosity objects. Seven of the objects have spectra resembling dwarf novae, but have apparently never been observed in outburst, suggesting that many cataclysmics with relatively low variability amplitude remain to be discovered. Some of the objects are notable. SDSS J07568+0858 and SDSS J08129+1911 were previously known to have deep eclipses; in addition to spectroscopy, we use archival data from the CRTTS to refine their periods. We give a parallax-based distance of 195 (+54, -39) pc for LV Cnc (SDSS J09197+0857), which at Porb = 81 m has the shortest orbital period in our sample. SDSS J08091+3814 shows both the spectroscopic phase offset and phase-dependent absorption found in SW Sextantis stars. The average spectra of SDSS J08055+0720 and SDSS J16191+1351 show contributions from K-type secondaries, and SDSS J080440+0239 shows a contribution from an early M star. We use these to constrain the distances. SDSS J09459+2922 has characteristics typical of a magnetic system. SDSS11324+6249 may be a novalike variable, and if so, its orbital period (99 min) is unusually short for that subclass.
We investigate the star formation activity in a young star forming cluster embedded at the edge of the RCW 41 HII region. As a complementary goal, we aim at demonstrating the gain provided by Wide-Field Adaptive Optics instruments to study young clusters. We used deep, JHKs images from the newly commissioned Gemini-GeMS/GSAOI instrument, complemented with Spitzer IRAC observations, in order to study the photometric properties of the young stellar cluster. GeMS is an AO instrument, delivering almost diffraction limited images over a field of 2' across. The exquisite angular resolution allows us to reach a limiting magnitude of J = 22 for 98% completeness. The combination of the IRAC photometry with our JHKs catalog is used to build color-color diagrams, and select Young Stellar Objects (YSOs) candidates. We detect the presence of 80 Young Stellar Object (YSO) candidates. Those YSOs are used to infer the cluster age, which is found to be in the range 1 to 5 Myr. We find that 1/3 of the YSOs are in a range between 3 to 5 Myr, while 2/3 of the YSO are < 3 Myr. When looking at the spatial distribution of these two populations, we evidence a potential age gradient across the field, suggesting sequential star formation. We construct the IMF, and show that we can sample the mass distribution well into the brown dwarf regime (down to 0.01 Msun). The logarithmic mass function rises to peak at 0.3 Msun, before turning over and declining into the brown dwarf regime. The total cluster mass derived is estimated to be 78 +/- 18 Msun, while the ratio of brown dwarfs to star derived is 18 p/- 5 %. When comparing with other young clusters, we find that the IMF shape of the young cluster embedded within RCW 41 is consistent with those of Trapezium, IC 348 or Chamaeleon I, except for the IMF peak, which happens to be at higher mass. This characteristic is also seen in clusters like NGC 6611 or even Taurus.
We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.
We present the calculation of coherent radio pulses emitted by extensive air showers induced by ultra-high energy cosmic rays accounting for reflection on the Earth's surface. Results have been obtained with a simulation program that calculates the contributions from shower particles after reflection at a surface plane. The properties of the radiation are discussed in detail emphasizing the effects of reflection. The shape of the frequency spectrum is shown to be closely related to the angle of the observer with respect to shower axis, becoming hardest in the Cherenkov direction. The intensity of the flux at a fixed observation angle is shown to scale with the square of the primary particle energy to very good accuracy indicating the coherent aspect of the emission. The simulation methods of this paper provide the foundations for energy reconstruction of experiments looking at the Earth from balloons and satellites. They can also be used in dedicated studies of existing and future experimental proposals.
Currently, the analysis of transmission spectra is the most successful technique to probe the chemical composition of exoplanet atmospheres. But the accuracy of these measurements is constrained by observational limitations and the diversity of possible atmospheric compositions. Here we show the UV-VIS-IR transmission spectrum of Jupiter, as if it were a transiting exoplanet, obtained by observing one of its satellites, Ganymede, while passing through Jupiter's shadow i.e., during a solar eclipse from Ganymede. The spectrum shows strong extinction due to the presence of clouds (aerosols) and haze in the atmosphere, and strong absorption features from CH4. More interestingly, the comparison with radiative transfer models reveals a spectral signature, which we attribute here to a Jupiter stratospheric layer of crystalline H2O ice. The atomic transitions of Na are also present. These results are relevant for the modeling and interpretation of giant transiting exoplanets. They also open a new technique to explore the atmospheric composition of the upper layers of Jupiter's atmosphere.
Since the discovery of a neutrino flux in excess of the atmospheric background by the IceCube Collaboration, searches for the astrophysical sources have been ongoing. Due to the steeply falling background towards higher energies, the PeV events detected in three years of IceCube data are the most likely ones to be of extraterrestrial origin. Even excluding the PeV events detected so far, the neutrino flux is well above the atmospheric background, so it is likely that a number of sub-PeV events originate from the same astrophysical sources that produce the PeV events. We study the high-energy properties of AGN that are positionally coincident with the neutrino events from three years of IceCube data and show the results for event number 4. IC 4 is a event with a low angular error (7.1$^\circ$) and a large deposited energy of 165 TeV. We use multiwavelength data, including Fermi/LAT and X-ray data, to construct broadband spectra and present parametrizations of the broadband spectral energy distributions with logarithmic parabolas. Assuming the X-ray to {\gamma}-ray emission in blazars originates in the photoproduction of pions by accelerated protons, their predicted neutrino luminosity can be estimated. The measurements of the diffuse extragalactic background by Fermi/LAT gives us an estimate of the flux contributions from faint unresolved blazars. Their contribution increases the number of expected events by a factor of $\sim$2. We conclude that the detection of the IceCube neutrinos IC4, IC14, and IC20 can be explained by the integral emission of blazars, even though no individual source yields a sufficient energy output.
Long-slit medium-resolution spectra of the Galactic globular clusters (GCs) NGC6229 and NGC6779, obtained with the CARELEC spectrograph at the 1.93-m telescope of the Haute-Provence observatory, have been used to determine the age, helium abundance (Y), and metallicity [Fe/H] as well as the first estimate of the abundances of C, N, O, Mg, Ca, Ti, and Cr for these objects. We solved this task by comparing the observed spectra and the integrated synthetic spectra, calculated with the use of the stellar atmosphere models with the parameters preset for the stars from these clusters. The model mass estimates, $T_{\rm eff}$, and log~g were derived by comparing the observed "color-magnitude" diagrams and the theoretical isochrones. The summing-up of the synthetic blanketed stellar spectra was conducted according to the Chabrier mass function. To test the accuracy of the results, we estimated the chemical abundances, [Fe/H], $\log~t$, and $Y$ for the NGC5904 and NGC6254 clusters, which, according to the literature, are considered to be the closest analogues of the two GCs of our study. Using the medium-resolution spectra from the library of Schiavon et al., we obtained for these two clusters a satisfactory agreement with the reported estimates for all the parameters within the errors. We derived the following cluster parameters. NGC6229: [Fe/H]=-1.65 dex, t=12.6 Gyr, Y=0.26, [\alpha/Fe]=0.28 dex; NGC6779: [Fe/H]=-1.9 dex, t=12.6 Gyr, Y=0.23, [\alpha/Fe]=0.08 dex; NGC5904: [Fe/H]=-1.6 dex, t=12.6 Gyr, Y=0.30, [\alpha/Fe]=0.35 dex; NGC6254: [Fe/H]=-1.52 dex, t=11.2 Gyr, Y=0.30 [\alpha/Fe]=0.025 dex. The value [\alpha/Fe] denotes the average of the Ca and Mg abundances.
The well-known correlation between the radio luminosity ($L_R$) and the X-ray luminosity ($L_X$) $L_R / L_X \simeq 10^{-5}$ holds for a variety of objects like active galactic nuclei, galactic black holes, solar flares and cool stars. Here we extend the relation to gamma-ray bursts (GRBs), and find the GRBs also lay on the same $L_R-L_X$ relation, with a slightly different slope as $L_R \propto L_X^{1.1}$. This relation implies the explosions in different scales may have a common underlying origin.
The behavior of strong gravitational lens model software in the analysis of lens models is not necessarily consistent among the various software available, suggesting that the use of several models may enhance the understanding of the system being studied. Among the publicly available codes, the model input files are heterogeneous, making the creation of multiple models tedious. An enhanced method of creating model files and a method to easily create multiple models, may increase the number of comparison studies. HydraLens simplifies the creation of model files for four strong gravitational lens model software packages, including Lenstool, Gravlens/Lensmodel, glafic and PixeLens, using a custom designed GUI for each of the four codes that simplifies the entry of the model for each of these codes, obviating the need for user manuals to set the values of the many flags and in each data field. HydraLens is designed in a modular fashion, which simplifies the addition of other strong gravitational lens codes in the future. HydraLens can also translate a model generated for any of these four software packages into any of the other three. Models created using HydraLens may require slight modifications, since some information may be lost in the translation process. However the computer generated model greatly simplifies the process of developing multiple lens models. HydraLens may enhance the number of direct software comparison studies, and also assist in the education of young investigators in gravitational lens modeling. Future development of HydraLens will further enhance its capabilities.
ASTRI SST-2M is the end-to-end prototype telescope of the Italian National Institute of Astro- physics, INAF, designed to investigate the 10-100 TeV band in the framework of the Cherenkov Telescope Array, CTA. The ASTRI SST-2M telescope has been installed in Italy in September 2014, at the INAF ob- serving station located at Serra La Nave on Mount Etna. The telescope is foreseen to be completed and fully operative in spring 2015 including auxiliary instrumentation needed to support both operations and data anal- ysis. In this contribution we present the current status of a sub-set of the auxiliary instruments that are being used at the Serra La Nave site, namely an All Sky Camera, an Electric Field Meter and a Raman Lidar devoted, together with further instrumentation, to the monitoring of the atmospheric and environmental conditions. The data analysis techniques under development for these instruments could be applied at the CTA sites, where similar auxiliary instrumentation will be installed.
Blazars show rapid and violent variabilities, which timescale are often less than a day. We studied intraday variations by applying a shot analysis technique to Kepler monitoring of blazar W2R 1926+42 in Quarter 14. We obtained a mean profile calculated from 195 rapid variations. The mean profile shows three components; one is a sharp structure distributed within $\pm$0.1 day of the peak, and two slow-varying components. This spiky-peak component reflects features of rapid variations directly. The profile of peak component shows an exponential rise and decay of which timescales are different, 0.0416 and 0.0588 day respectively. This component is too sharp to represent a standard function which is often used to express blazar variations. This asymmetric profile at the peak is difficult to be explained by a simple variation of the Doppler factor by changing a geometry of the emitting region. This result indicates that intraday variations arise from a production of high-energy accelerated particles in the jet.
In this paper we operate under the assumption that no tensors from inflation will be measured in the future by the dedicated experiments and argue that, while for single-field slow-roll models of inflation the running of the spectral index will be hard to be detected, in multi-field models the running can be large due to its strong correlation with non-Gaussianity. A detection of the running might therefore be related to the presence of more than one active scalar degree of freedom during inflation.
We present a novel method to implement time-delayed propagation of radiation fields in cosmological radiative transfer simulations. Time-delayed propagation of radiation fields requires construction of retarded-time fields by tracking the location and lifetime of radiation sources along the corresponding light-cones. Cosmological radiative transfer simulations have, until now, ignored this "light-cone effect" or implemented ray-tracing methods that are computationally demanding. We show that radiative transfer calculation of the time-delayed fields can be easily achieved in numerical simulations when periodic boundary conditions are used, by calculating the time-discretized retarded-time Green's function using the Fast Fourier Transform (FFT) method and convolving it with the source distribution. We also present a direct application of this method to the long-range radiation field of Lyman-Werner band photons, which is important in the high-redshift astrophysics with first stars.
We report the discovery of two new low-mass, thermally bloated, hot white dwarfs among the Kepler sample of eclipsing binaries. These are KIC 9164561 and KIC 10727668 with orbital periods of 1.2670 and 2.3058 days, respectively. The current primary in both systems is an A star of about 2 Msun. This brings the number of similar binaries among the Kepler sample to six, and the two new systems have the shortest orbital periods among them. The white dwarf in KIC 9164561 has the largest thermal bloating, compared to its cold degenerate radius, of about a factor of 14. We utilize RV measurements of the A star in KIC 9164561 to determine the white dwarf mass rather accurately: 0.197 +/- 0.005 Msun. The mass of the white dwarf in KIC 10727668 is based on the Doppler boosting signal in the Kepler photometry, and is less accurately determined to be 0.266 +/- 0.035 Msun. Based on the inferred radii and effective temperatures of these two white dwarfs we are able to make an independent theoretical estimate of their masses to within ~0.01 Msun based on evolutionary models of their cooling history after they lose their hydrogen-rich envelopes. We also present evidence that there is a third body in the KIC 9164561 system with an orbital period of 8-14 years.
The HR 8799 system, with its four giant planets and two debris belts, has an architecture closely mirroring that of our Solar system where the inner, warm asteroid belt and outer, cool Edgeworth-Kuiper belt bracket the giant planets. As such, it is a valuable laboratory for examining exoplanetary dynamics and debris disk-exoplanet interactions. Whilst the outer debris belt of HR 8799 has been well resolved by previous observations, the spatial extent of the inner disk remains unknown. This leaves a significant question mark over both the location of the planetesimals responsible for producing the belt's visible dust and the physical properties of those grains. We have performed the most extensive simulations to date of the inner, unresolved debris belt around HR 8799, using UNSW Australia's Katana supercomputing facility to follow the dynamical evolution of a model inner disk comprising 300,298 particles for a period of 60 million years. These simulations have enabled the characterisation of the extent and structure of the inner disk in detail, and will in future allow us to provide a first estimate of the small-body impact rate and water delivery prospects for possible (as-yet undetected) terrestrial planet(s) in the inner system.
Small hydrocarbons, such as C2H, C3H and C3H2 are more abundant in photo-dissociation regions (PDRs) than expected based on gas-phase chemical models. To explore the hydrocarbon chemistry further, we observed a key intermediate species, the hydrocarbon ion l-C3H+, in the Horsehead PDR with the Plateau de Bure Interferometer at high-angular resolution (6''). We compare with previous observations of C2H and c-C3H2 at similar angular resolution and new gas-phase chemical model predictions to constrain the dominant formation mechanisms of small hydrocarbons in low-UV flux PDRs. We find that, at the peak of the HCO emission (PDR position), the measured l-C3H+, C2H and c-C3H2 abundances are consistent with current gas-phase model predictions. However, in the first PDR layers, at the 7.7 mum PAH band emission peak, which are more exposed to the radiation field and where the density is lower, the C2H and c-C3H2 abundances are underestimated by an order of magnitude. At this position, the l-C3H+ abundance is also underpredicted by the model but only by a factor of a few. In addition, contrary to the model predictions, l-C3H+ peaks further out in the PDR than the other hydrocarbons, C2H and c-C3H2. This cannot be explained by an excitation effect. Current gas-phase photochemical models thus cannot explain the observed abundances of hydrocarbons, in particular in the first PDR layers. Our observations are consistent with a top-down hydrocarbon chemistry, in which large polyatomic molecules or small carbonaceous grains are photo-destroyed into smaller hydrocarbon molecules/precursors.
We investigate the spatial distribution of charged particles accelerated by non-relativistic oblique fast collisionless shocks using three-dimensional test-particle simulations. We find that the density of low-energy particles exhibit a localised enhancement at the shock, resembling the "spike" measured at interplanetary shocks. In contrast to previous results based on numerical solutions to the focused transport equation, we find a shock spike for any magnetic obliquity, from quasi-perpendicular to parallel. We compare the pitch-angle distribution with respect to the local magnetic field and the momentum distribution far downstream and very near the shock within the spike; our findings are compatible with predictions from the scatter-free shock drift acceleration (SDA) limit in these regions. The enhancement of low-energy particles measured by Voyager 1 at solar termination shock is comparable with our profiles. Our simulations allow for predictions of supra-thermal protons at interplanetary shocks within ten solar radii to be tested by Solar Probe Mission. They also have implications for the interpretation of ions accelerated at supernova remnant shocks.
We investigate the effect of the axion cooling on the nucleosynthesis in a massive star with $16M_{\odot}$ by standard stellar evolution calculation. We find that the axion cooling suppresses the nuclear reactions in carbon, oxygen and silicon burning phases because of the extraction of the energy. As a result, larger amounts of the already synthesized neon and magnesium remain without being consumed to produce further heavier elements. Even in the case with the axion-photon coupling constant $g_{a\gamma}= 10^{-11}$ GeV$^{-1}$, which is six times smaller than the current upper limit, the amount of neon and magnesium that remain just before the core-collapse supernova explosion is considerably larger than the standard value. This implies that we could give a more stringent constraint on $g_{a\gamma}$ from the nucleosynthesis of heavy elements in massive stars.
Some asteroids eject dust, producing transient, comet-like comae and tails; these are the active asteroids. The causes of activity in this newly-identified population are many and varied. They include impact ejection and disruption, rotational instabilities, electrostatic repulsion, radiation pressure sweeping, dehydration stresses and thermal fracture, in addition to the sublimation of asteroidal ice. These processes were either unsuspected or thought to lie beyond the realm of observation before the discovery of asteroid activity. Scientific interest in the active asteroids lies in their promise to open new avenues into the direct study of asteroid destruction, the production of interplanetary debris, the abundance of asteroid ice and the origin of terrestrial planet volatiles.
During a search for gamma-ray emission from NGC 3628 (Arp 317), two new unidentified gamma-ray sources, Fermi J1049.7+0435 and J1103.2+1145 have been discovered \cite{ATel}. The detections are made in data from the Large Area Telescope (LAT), on board the Fermi Gamma-ray Space Telescope, in the 100\,MeV to 300\,GeV band during the period between 2008 August 5 and 2012 October 27. Neither is coincident with any source listed in the 2FGL catalogue \cite{Nolan2012}. Fermi J1049.7+0435 is at Galactic coordinates $(l,b) = (245.34^\circ, 53.27^\circ)$, $(\alpha_{J2000}, \delta_{J2000}) = (162.43^\circ, 4.60^\circ)$. Fermi J1103.2+1145 is at Galactic coordinates $(l,b) = (238.85^\circ, 60.33^\circ)$, $(\alpha_{J2000},\delta_{J2000})= (165.81^\circ, 11.75^\circ)$. Possible radio counterparts are found for both sources, which show flat radio spectra similar to other Fermi LAT detected AGN, and their identifications are discussed. These identification have been supoorted by snap-shot observations with the Australia Telescope Compact Array at several epochs in 2013 and 2014,
We propose a model to explain the ultra-bright GeV gamma-ray flares observed
from the blazar \c454. The model is based on the concept of a relativistic jet
interacting with compact gas condensations produced when a star (red giant)
crosses the jet close to the central black hole. The study includes an
analytical treatment of the evolution of the envelop lost by the star within
the jet, and calculations of the related high-energy radiation.
The model readily explains the day-long, variable on timescales of hours, GeV
gamma-ray flare from \c454, observed during November 2010 on top of a
weeks-long plateau. In the proposed scenario, the plateau state is caused by a
strong wind generated by the heating of the star atmosphere by nonthermal
particles accelerated at the jet-star interaction region. The flare itself
could be produced by a few clouds of matter lost by the red giant after the
initial impact of the jet. In the framework of the proposed scenario, the
observations constrain the key model parameters of the source, including the
mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total
jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor
of the gamma-ray emitting clouds, $\delta\simeq 20$. Whereas we do not specify
the particle acceleration mechanisms, the potential gamma-ray production
processes are discussed and compared in the context of the proposed model. We
argue that synchrotron radiation of protons has certain advantages compared to
other radiation channels of directly accelerated electrons.
Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent.
The lightest Kaluza-Klein particle (LKP), which appears in the theory of universal extra dimensions, is one of the good candidates for cold dark matter. The gamma-ray spectrum from annihilation of LKP dark matter shows a characteristic peak structure around the LKP mass. We investigate the detectability of this peak structure by considering energy resolution of near-future detectors, and calculate the expected count spectrum of the gamma-ray signal. In order to judge whether the count spectrum contains the LKP signal, the {\chi} squared test is employed. If the signal is not detected, we set some constraints on the boost factor that is an uncertain factor dependent on the substructure of the LKP distribution in the galactic halo. Detecting such peak structure would be conclusive evidence that dark matter is made of LKP.
We are carrying out a series of photometric monitoring to measure the rotation periods of members in the young $\beta$ Pictoris Association, as part of the RACE-OC project (Rotation and ACtivity Evolution in Open Clusters). In this paper, we present the results for HD 155555C which is believed to be physically associated to the spectroscopic binary V824 Ara (HD155555) and thus constituting a triple system. We collected B, V, and R-band photometric data timeseries and discovered from periodogram analysis the rotation period P = 4.43d. Combined with stellar radius and projected rotational velocity, we find this star almost equator-on with an inclination $i$ $\simeq$ 90$^{\circ}$. The rotational properties of HD155555C fit well into the period distribution of other $\beta$ Pic members, giving further support to the suggested membership to the association and to its physical association to V824 Ara. A comparison with Pre-Main-Sequence isochrones from various models allows us to estimate an age of 20$\pm$15 Myr for this triple system.
Since August 2009, MAXI experiment on the ISS has been performing all-sky X-ray monitoring. With MAXI, we detected flaring activities of some blazers, including Mrk 421, Mrk 501, and 3C 273. Recently, new X-ray flaring activities were detected from two blazers, MAXI J1930+093 = 2FGL J1931.1+0938 (Atel#5943) and 2MAXI J0243-582 = BZB J0244-5819 (Atel#6012). The MAXI monitoring also covers black hole binaries, including Cyg X-1 and Cyg X-3 which emit GeV gamma-rays. Their gamma-ray emission was found to coincide with their X-ray state transitions. We present light curves and outstanding events of these sources.
Angular momentum transport and particle acceleration during the magnetorotational instability (MRI) in a collisionless accretion disk are investigated using three-dimensional particle-in-cell (PIC) simulation. We show that the kinetic MRI can provide not only high energy particle acceleration but also enhancement of angular momentum transport. We find that the plasma pressure anisotropy inside the channel flow with $p_{\|} > p_{\perp}$ induced by active magnetic reconnection suppresses the onset of subsequent reconnection, which in turn leads to high magnetic field saturation and enhancement of Maxwell stress tensor of angular momentum transport. Meanwhile, during the quiescent stage of reconnection the plasma isotropization progresses in the channel flow, and the anisotropic plasma with $p_{\perp} > p_{\|}$ due to the dynamo action of MRI outside the channel flow contributes to rapid reconnection and strong particle acceleration. This efficient particle acceleration and enhanced angular momentum transport in a collisionless accretion disk may explain the origin of high energy particles observed around massive black holes.
We study the density profile and shape of the Galactic halo using deep multicolour images from the MENeaCS and CCCP projects, over 33 fields selected to avoid overlap with the Galactic plane. Using multicolour selection and PSF homogenization techniques we obtain catalogues of F stars (near-main sequence turnoff stars) out to Galactocentric distances up to 60kpc. Grouping nearby lines of sight, we construct the stellar density profiles through the halo in eight different directions by means of photometric parallaxes. Smooth halo models are then fitted to these profiles. We find clear evidence for a steepening of the density profile power law index around R=20 kpc, from -2.50 +- 0.04 to -4.85 +- 0.04, and for a flattening of the halo towards the poles with best-fit axis ratio 0.63 +- 0.02. Furthermore, we cannot rule out a mild triaxiality (w>=0.8). We recover the signatures of well-known substructure and streams that intersect our lines of sight. These results are consistent with those derived from wider but shallower surveys, and augur well for upcoming, wide-field surveys of comparable depth to our pencil beam surveys.
Elastic properties of the solid regions of neutron star crusts and white dwarfs play an important role in theories of stellar oscillations. Matter in compact stars is presumably polycrystalline and, since the elastic properties of single crystals of such matter are very anisotropic, it is necessary to relate elastic properties of the polycrystal to those of a single crystal. We calculate the effective shear modulus of polycrystalline matter with randomly oriented crystallites using a self-consistent theory that has been very successful in applications to terrestrial materials and show that previous calculations overestimate the shear modulus by approximately 28%.
The SEDs of 18 GeV FSRQs are collected and compiled from literature, in which both the jet emission and the accretion disk radiation can be observed, in order to investigate the correlations among their jet power (P_jet), accretion disk luminosity (L_disk), and luminosity of broad line region (BLR, L_BLR). On the basis of the SED fits with the jet radiation and accretion disk radiation models, we calculate P_jet and L_disk. No correlation between P_jet with either L_disk or L_BLR is found. With a sub-sample of L_BLR for 13 GeV FSRQs, it is observed that L_BLR is strongly correlated with their L_disk. We also study the BLR covering factors of the GeV FSRQs in our sample, averagely which are smaller than that of the large samples of radio-loud and radio-quiet quasars. P_jet of some GeV FSRQs is higher than L_disk, but P_jet of all the GeV FSRQs is lower than the accretion power of black hole (BH), which is estimated by \dot{M}c^2=L_disk/0.1, indicating that the total accretion power of BH is sufficient to drive the jets in these sources; however the uncorrelation between L_disk and P_jet of the GeV FSRQs may suggest that their jets are launched by the BZ process via extracting the rotational energy of BH. Using the L_BLR-L_disk relation of the GeV FSRQs, we estimate L_disk of a BL Lac sample with their L_BLR. A comparison of L_BLR and Eddington ratio (L_disk/L_Edd) among BL Lacs, very radio-loud NLS1 galaxies, and FSRQs is also presented. It is found that along with the BL Lac-NLS1-FSRQ sequence L_BLR and L_disk/L_Edd increase, which may correspond to the change of the accretion disk structure and the transformation of the dominant mechanism for jet launching.
We employ different shapes of apodizing windows in the local correlation tracking (LCT) routine to retrieve horizontal velocities using numerical simulations of convection. LCT was applied on a time sequence of temperature maps generated by the Nirvana code with four different apodizing windows, namely--Gaussian, Lorentzian, trapezoidal and triangular, with varying widths. In terms of correlations (between the LCT-retrieved and simulated flow field), the triangular and the trapezoidal perform the best and worst, respectively. On segregating the intrinsic velocities in the simulations on the basis of their magnitudes, we find that for all windows, a significantly higher correlation is obtained for the intermediate and high-velocity bins and only modest or weak values in the low-velocity bins. The differences between the LCT-retrieved and simulated flow fields were determined spatially which show large residuals at or close to the boundary of granules. The extent to which the horizontal flow vectors retrieved by LCT compare with the simulated values, depends entirely on the width of the central peak of the apodizing window for a given $\sigma$. Even though LCT suffers from a lack of spatial content as seen in simulations, its simplicity and speed could serve as a viable first-order tool to probe horizontal flows--one that is ideal for large data sets.
We investigate the impact on the classical dynamics of dark matter particles and dark energy of a non-minimal coupling in the dark sector, assuming that the mass of the dark matter particles is coupled to a dark energy scalar field. We show that standard results can only be recovered if the space-time variation of the dark energy scalar field is sufficiently smooth on the characteristic length scale of the dark matter particles, and we determine the associated constraint dependent on both the mass and radius of the dark matter particles and the coupling to the dark energy scalar field. We further show, using field theory numerical simulations, that a violation of such constraint results in a microscopic feedback effect strongly affecting the dynamics of dark matter particles, with a potential impact on structure formation and on the space-time evolution of the dark energy equation of state.
In this paper, we constrain the tilt of the power spectrum of relic gravitational waves by combining the data from BICEP2/Keck array and Planck (BKP) and the Laser Interferometer Gravitational-Waves Observatory (LIGO). From the data of BKP B-modes, the constraint on the tensor tilt is $n_t=0.66^{+1.83}_{-1.44}$ at the $68%$ confidence level. By further adding the LIGO upper limit on the energy density of gravitational waves, the constraint becomes $n_t=-0.76^{+1.37}_{-0.52}$ at the $68%$ confidence level. We conclude that there is no evidence for a blue-tilted power spectrum of relic gravitational waves and either sign of the index of tensor power spectrum is compatible with the data.
The Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics
Observatory (SDO) observes the Sun at the Fe I 6173 {\AA} line and returns
full-disk maps of line-of-sight (LOS) observables including the magnetic flux
density, velocities, Fe I line width, line depth, and continuum intensity.
These data are estimated through an algorithm (the MDI-like algorithm,
hereafter), which combines observables obtained at six wavelength positions
within the Fe I 6173 {\AA} line. To properly interpret such data it is
important to understand any effects of the instrument and the pipeline that
generates these data products. We tested the accuracy of the line width, line
depth, and continuum intensity returned by the MDI-like algorithm using various
one-dimensional (1D) atmosphere models. It was found that HMI estimates of
these quantities are highly dependent on the shape of the line, therefore on
the LOS angle and the magnetic flux density associated with the model, and less
to line shifts with respect to the central positions of the instrument
transmission profiles. In general, the relative difference between synthesized
values and HMI estimates increases toward the limb and with the increase of the
field; the MDI-like algorithm seems to fail in regions with fields larger than
approximately 2000 G.
Instrumental effects were investigated by the analysis of HMI data obtained
at daily intervals for a span of three years at disk center in the quiet Sun
and hourly intervals for a span of 200 hours. The analysis revealed
periodicities induced by the variation of the orbital velocity of the
observatory with respect to the Sun, and long-term trends attributed to
instrument adjustments, re-calibrations and instrumental degradation.
In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4<z<6.5$ and absolute magnitude interval $-22.5\lesssim M_{1450} \lesssim -18.5$ typical of local Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at $z\sim 4-6$ than that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at $z=4-6.5$ photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed in the Lyman-$\alpha$ forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization.
We present the first interpretation of the new isotropic diffuse gamma-ray background (IGRB), measured by the Fermi Large Area Telescope (LAT), based on a statistical analysis. We demonstrate that the gamma-ray emission from unresolved active galactic nuclei and star forming galaxies is consistent with the Fermi-LAT IGRB data within the uncertainties both on the choice of the Galactic diffuse emission model and on the gamma-ray emission mechanism of these sources. Furthermore, adding to the extragalactic sources the contribution from a smooth Galactic halo of annihilating weakly interacting dark matter (DM) particles, we are able to set stringent limits on the DM annihilation cross section. Finally, we demonstrate that the addition of DM can significantly improve the fit to IGRB data.
The First G-APD Cherenkov Telescope (FACT) was built on the Canary Island of La Palma in October 2011 as a proof of principle for silicon based photosensors in Cherenkov Astronomy. The scientific goal of the project is to study the variability of active galatic nuclei (AGN) at TeV energies. Observing a small sample of TeV blazars whenever possible, an unbiased data sample is collected. This allows to study the variability of the selected objects on timescales from hours to years. Results from the first three years of monitoring will be presented. To provide quick flare alerts to the community and trigger multi-wavelength observations, a quick look analysis has been installed on-site providing results publicly online within the same night. In summer 2014, several flare alerts were issued. Results of the quick look analysis are summarized.
We have obtained high-resolution spectra of 32 giants in the open cluster NGC 7789 using the Wisconsin-Indiana-Yale-NOAO Hydra spectrograph. We explore differences in atmospheric parameters and elemental abundances caused by the use of the linelist developed for the Gaia-ESO Survey (GES) compared to one based on Arcturus used in our previous work. [Fe/H] values decrease when using the GES linelist instead of the Arcturus-based linelist; these differences are probably driven by systematically lower (~ -0.1 dex) GES surface gravities. Using the GES linelist we determine abundances for 10 elements - Fe, Mg, Si, Ca, Ti, Na, Ni, Zr, Ba, and La. We find the cluster's average metallicity [Fe/H] = 0.03 +/- 0.07 dex, in good agreement with literature values, and a lower [Mg/Fe] abundance than has been reported before for this cluster (0.11 +/- 0.05 dex). We also find the neutron-capture element barium to be highly enhanced - [Ba/Fe] = +0.48 +/- 0.08 - and disparate from cluster measurements of neutron-capture elements La and Zr (-0.08 +/- 0.05 and 0.08 +/- 0.08, respectively). This is in accordance with recent discoveries of supersolar Ba enhancement in young clusters along with more modest enhancement of other neutron-capture elements formed in similar environments.
Soft lags from the emission of the lower kilohertz quasi-periodic oscillations (kHz QPOs) of neutron star low mass X-ray binaries have been reported from 4U1608-522 and 4U1636-536. Those lags hold prospects for constraining the origin of the QPO emission, including the location at which the oscillation takes place, a stepping stone before we can use the kHz QPOs to probe strong field General Relativity. In this paper, we investigate the spectral-timing properties of both the lower and upper kHz QPOs from the neutron star binary 4U1728-34, in which the duty cycles of both QPOs are comparable, using the entire Rossi X-ray Timing Explorer archive on this source. We show that the lag-energy spectra of the two QPOs are systematically different: while the lower kHz QPO shows soft lags, the upper kHz QPO shows either a flat lag-energy spectrum or hard variations lagging softer variations. This suggests two different QPO-generation mechanisms. We also computed the first covariance spectra for both kHz QPOs and performed a spectral deconvolution. The QPO spectra are consistent with Comptonized blackbody emission, similar to the one found from the decomposition of the time-averaged continuum, but with a higher seed-photon temperature, which could suggest that a more compact inner region of the Comptonization layer (boundary/spreading layer, corona) is responsible for the QPO emission. Considering our results together with other recent findings, we suggest that the lower kHz QPO signal is generated by coherent oscillations of the compact boundary layer region itself. The upper kHz QPO signal is linked to less-coherent accretion-rate variations produced in the inner accretion disk. However, the disk emission being extended, the modulation is only observed when the accretion variations reach the boundary layer.
Infrared spectroscopic observations have established the presence of hydrocarbon ices on Pluto and other TNOs, but the abundances of such molecules cannot be deduced without accurate optical constants (n, k) and reference spectra. In this paper we present our recent measurements of near- and mid-infrared optical constants for ethane (C$_2$H$_6$) and ethylene (C$_2$H$_4$) in multiple ice phases and at multiple temperatures. As in our recent work on acetylene (C$_2$H$_2$), we also report new measurements of the index of refraction of each ice at 670 nm. Comparisons are made to earlier work where possible, and electronic versions of our new results are made available.
Accurate rate coefficients for molecular vibrational transitions due to collisions with H$_2$, critical for interpreting infrared astronomical observations, are lacking for most molecules. Quantum calculations are the primary source of such data, but reliable values that consider all internal degrees of freedom of the collision complex have only been reported for H$_2$-H$_2$ due to the difficulty of the computations. Here we present essentially exact full-dimensional dynamics computations for rovibrational quenching of CO due to H$_2$ impact. Using a high-level six-dimensional potential surface, time-independent scattering calculations, within a full angular-momentum-coupling formulation, were performed for the deexcitation of vibrationally excited CO. Agreement with experimentally-determined results confirms the accuracy of the potential and scattering computations, representing the largest of such calculations performed to date. This investigation advances computational quantum dynamics studies representing initial steps toward obtaining CO-H$_2$ rovibrational quenching data needed for astrophysical modeling.
A full energy and flavor-dependent analysis of the three-year high-energy IceCube neutrino events is presented. By means of multi-dimensional fits, we derive the current preferred values of the high-energy neutrino flavor ratios, the normalization and spectral index of the astrophysical fluxes, and the expected atmospheric background events, including a prompt component. A crucial assumption resides on the choice of the energy interval used for the analyses, which significantly biases the results. When restricting ourselves to the $\sim$30 TeV - 3 PeV energy range, which contains all the observed IceCube events, we find that the inclusion of the spectral information improves the fit to the canonical flavor composition at Earth, ($1:1:1$)$_\oplus$, with respect to a single-energy bin analysis. Increasing both the minimum and the maximum deposited energies has dramatic effects on the reconstructed flavor ratios as well as on the spectral index. Imposing a higher threshold of 60 TeV yields a slightly harder spectrum by allowing a larger muon neutrino component, since above this energy most atmospheric track-like events are effectively removed. Extending the high-energy cutoff to fully cover the Glashow resonance region leads to a softer spectrum and a preference for tau neutrino dominance, as none of the expected electron (anti)-neutrino induced showers have been observed so far. The lack of showers at energies above 2 PeV may point to a broken power-law neutrino spectrum. Future data may confirm or falsify whether or not the recently discovered high-energy neutrino fluxes and the longstanding detected cosmic rays have a common origin.
We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma
rays) on Mrk421 between January 2009 and June 2009, which included VLBA,
F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other
instruments and collaborations. Mrk421 was found in its typical (non-flaring)
activity state, with a VHE flux of about half that of the Crab Nebula, yet the
light curves show significant variability at all wavelengths, the highest
variability being in the X-rays. We determined the power spectral densities
(PSD) at most wavelengths and found that all PSDs can be described by
power-laws without a break, and with indices consistent with pink/red-noise
behavior. We observed a harder-when-brighter behavior in the X-ray spectra and
measured a positive correlation between VHE and X-ray fluxes with zero time
lag. Such characteristics have been reported many times during flaring
activity, but here they are reported for the first time in the non-flaring
state. We also observed an overall anti-correlation between optical/UV and
X-rays extending over the duration of the campaign.
The harder-when-brighter behavior in the X-ray spectra and the measured
positive X-ray/VHE correlation during the 2009 multi-wavelength campaign
suggests that the physical processes dominating the emission during non-flaring
states have similarities with those occurring during flaring activity. In
particular, this observation supports leptonic scenarios as being responsible
for the emission of Mrk421 during non-flaring activity. Such a temporally
extended X-ray/VHE correlation is not driven by any single flaring event, and
hence is difficult to explain within the standard hadronic scenarios. The
highest variability is observed in the X-ray band, which, within the one-zone
synchrotron self-Compton scenario, indicates that the electron energy
distribution is most variable at the highest energies.
Galaxy clusters are the most massive systems in the known universe. They host relativistic cosmic ray populations and are thought to be gravitationally bound by large amounts of Dark Matter, which under the right conditions could yield to a detectable $\gamma$-ray flux. Prior to the launch of the Fermi satellite, predictions were optimistic that Galaxy clusters would be established as $\gamma$-ray bright objects by observations through its prime instrument, the Large Area Telescope (LAT). Yet, despite numerous efforts, even a single cluster detection is still pending.
We have studied the dynamics of an equal-mass magnetized neutron-star binary within a resistive magnetohydrodynamic (RMHD) approach in which the highly conducting stellar interior is matched to an electrovacuum exterior. Because our analysis is aimed at assessing the modifications introduced by resistive effects on the dynamics of the binary after the merger and through to collapse, we have carried out a close comparison with an equivalent simulation performed within the traditional ideal-MHD (IMHD) approximation. We have found that there are many similarities between the two evolutions, but also one important difference: the survival time of the hypermassive neutron star increases in a RMHD simulation. This difference is due to a less efficient magnetic-braking mechanism in the resistive regime, in which matter can move across magnetic-field lines, thus reducing the outward transport of angular momentum. Interestingly, a longer-lived magnetized hypermassive neutron star brings support to the recent modelling of short gamma-ray bursts in terms of such objects. Both the RMHD and the IMHD simulations carried here have been performed at higher resolutions and with a different grid structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L. Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732, L6 (2011)], but confirm the formation of a magnetic-jet structure in the low-density funnel produced by the black-hole--torus system. In both regimes the magnetic field is predominantly toroidal in the highly conducting torus and predominantly poloidal in the nearly evacuated funnel. Reconnection processes or neutrino annihilation occurring in the funnel, none of which we model here, could potentially increase the internal energy in the funnel and launch a relativistic outflow.
We consider searches for dark matter annihilation to monoenergetic neutrinos in the core of the Sun. We find that liquid scintillation neutrino detectors have enhanced sensitivity to this class of dark matter models, due to the energy and angular resolution possible for electron neutrinos and antineutrinos that scatter via charged-current interactions. In particular we find that KamLAND, utilizing existing data, could provide the best sensitivity to such models for $m_X \lesssim 14$ GeV.
We suggest the universe is Finslerian in the stage of inflation. The Finslerian background spacetime breaks rotational symmetry and induces parity violation. The primordial power spectrum is given for quantum fluctuation of the inflation field. It depends not only on the magnitude of wavenumber but also on the preferred direction. We derive the gravitational field equations in the perturbed Finslerian background spacetime, and obtain a conserved quantity outside the Hubble horizon. The angular correlation coefficients are presented in our anisotropic inflation model. The parity violation feature of Finslerian background spacetime requires that the anisotropic effect only appears in angular correlation coefficients if $l'=l+1$. The numerical results of the angular correlation coefficients are given to describe the anisotropic effect.
We expand the Einstein-Hilbert action with a positive cosmological constant up to the fourth order in terms of gravity fluctuations, and then use the in-in formalism to calculate the four-point correlation function for gravitational waves, including both contact and exchange diagrams, generated during a period of exactly de Sitter expansion. In addition, we also present the general properties of the $n$-point function of graviton in terms of both circularly and linearly polarized states.
Loop quantum cosmology tries to capture the main ideas of loop quantum gravity and to apply them to the Universe as a whole. Two main approaches within this framework have been considered to date for the study of cosmological perturbations: the dressed metric approach and the deformed algebra approach. They both have advantages and drawbacks. In this article, we accurately compare their predictions. In particular, we compute the associated primordial tensor power spectra. We show -- numerically and analytically -- that the large scale behavior is similar for both approaches and compatible with the usual prediction of general relativity. The small scale behavior is, the other way round, drastically different. Most importantly, we show that in a range of wavenumbers explicitly calculated, both approaches do agree on predictions that, in addition, differ from standard general relativity and do not depend on unknown parameters. These features of the power spectrum at intermediate scales might constitute a universal loop quantum cosmology prediction that can hopefully lead to observational tests and constraints. We also present a complete analytical study of the background evolution for the bouncing universe that can be used for other purposes.
We point out that the inflaton inevitably couples to all non-conformally coupled matters gravitationally through an oscillation in the Hubble parameter or the cosmic scale factor. It leads to particle production during the inflaton oscillation regime, which is most efficient just after inflation. Moreover, the analysis is extended to the model with non-minimal inflaton couplings to gravity, in which the Hubble parameter oscillates more violently. We apply our results to the graviton production by the inflaton: gravitons are also produced just after inflation, but the non-minimal coupling does not induce inflaton decay into the graviton pair.
In this work, I present exact cosmological solutions from Wesson's Induced Matter Model application to a general 5D metric in f(R,T) theory of gravity. The non-conservation of the energy-momentum tensor, predicted by f(R,T) theory, allows the derivation of a relation that describes the time evolution of the extra coordinate, revealing its compactification. It is showed that such a compactification could induce the effects of an accelerated expansion in the observable universe.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models, and inference with MCMC. By generating $\approx10^7$ model streams, we show that the geometry of the Pal 5 debris yields powerful constraints on the solar position and motion, the Milky Way and Pal 5 itself. All 10 model parameters were allowed to vary over large ranges without additional prior information. Using only SDSS data and a few radial velocities from the literature, we find that the distance of the Sun from the Galactic Center is $8.30\pm0.25$ kpc, and the transverse velocity is $253\pm16$ km/s. Both estimates are in excellent agreement with independent measurements of these quantities. Assuming a standard disk and bulge model, we determine the Galactic mass within Pal 5's apogalactic radius of 19 kpc to be $(2.1\pm0.4)\times10^{11}$ M$_\odot$. Moreover, we find the potential of the dark halo with a flattening of $q_z = 0.95^{+0.16}_{-0.12}$ to be essentially spherical within the radial range that is effectively probed by Pal 5. We also determine Pal 5's mass, distance and proper motion independently from other methods, which enables us to perform vital cross-checks. We conclude that with more observational data and by using additional prior information, the precision of this method can be significantly increased.
We combine HST/WFC3 imaging and G141 grism observations from the CANDELS and 3D-HST surveys to produce a catalog of grism spectroscopic redshifts for galaxies in the CANDELS/GOODS-South field. The WFC3/G141 grism spectra cover a wavelength range of 1.1<lambda<1.7 microns with a resolving power of R~130 for point sources, thus providing rest-frame optical spectra for galaxies out to z~3.5. The catalog is selected in the H-band (F160W) and includes both galaxies with and without previously published spectroscopic redshifts. Grism spectra are extracted for all H-band detected galaxies with H<24 and a CANDELS photometric redshift z_phot > 0.6. The resulting spectra are visually inspected to identify emission lines and redshifts are determined using cross-correlation with empirical spectral templates. To establish the accuracy of our redshifts, we compare our results against high-quality spectroscopic redshifts from the literature. Using a sample of 411 control galaxies, this analysis yields a precision of sigma_NMAD=0.0028 for the grism-derived redshifts, which is consistent with the accuracy reported by the 3D-HST team. Our final catalog covers an area of 153 square arcmin and contains 1019 redshifts for galaxies in GOODS-S. Roughly 60% (608/1019) of these redshifts are for galaxies with no previously published spectroscopic redshift. These new redshifts span a range of 0.677 < z < 3.456 and have a median redshift of z=1.282. The catalog contains a total of 234 new redshifts for galaxies at z>1.5. In addition, we present 20 galaxy pair candidates identified for the first time using the grism redshifts in our catalog, including four new galaxy pairs at z~2, nearly doubling the number of such pairs previously identified.
The Fe Kalpha emission line is the most ubiquitous feature in the X-ray spectra of active galactic nuclei (AGN), but the origin of its narrow core remains uncertain. Here, we investigate the connection between the sizes of the Fe Kalpha core emission regions and the measured sizes of the dusty tori in 13 local Type 1 AGN. The observed Fe Kalpha emission radii (R_fe) are determined from spectrally resolved line widths in X-ray grating spectra, and the dust sublimation radii (R_dust) are measured either from optical/near-infrared reverberation time lags or from resolved near-infrared interferometric data. This direct comparison shows that the dust sublimation radius forms an outer envelope to the bulk of the Fe Kalpha emission. R_fe matches R_dust well in the AGN with the best constrained line widths currently. In a significant fraction of objects without a clear narrow line core, R_fe is similar to, or smaller than the radius of the optical broad line region. These facts place important constraints on the torus geometries for our sample. Extended tori in which the solid angle of fluorescing gas peaks at well beyond the dust sublimation radius can be ruled out. We also test for luminosity scalings of R_fe, finding that Eddington ratio is not a prime driver in determining the line location in our sample. Large uncertainties on the line core widths preclude more detailed investigations at present, a limitation which Astro-H will help to overcome.
We describe a new mechanism that leads to the destabilisation of non-axisymmetric waves in astrophysical discs with an imposed radial temperature gradient. This might apply, for example, to the outer parts of protoplanetary discs. We use linear density wave theory to show that non-axisymmetric perturbations generally do not conserve their angular momentum in the presence of a forced temperature gradient. This implies an exchange of angular momentum between linear perturbations and the background disc. In particular, when the disturbance is a low-frequency trailing wave and the disc temperature decreases outwards, this interaction is unstable and leads to the growth of the wave. We demonstrate this phenomenon through numerical hydrodynamic simulations of locally isothermal discs in 2D using the FARGO code and in 3D with the ZEUS-MP and PLUTO codes. We consider radially structured discs with a self-gravitating region which remains stable in the absence of a temperature gradient. However, when a temperature gradient is imposed we observe exponential growth of a one-armed spiral mode (azimuthal wavenumber $m=1$) with co-rotation radius outside the bulk of the spiral arm, resulting in a nearly-stationary one-armed spiral pattern. The development of this one-armed spiral does not require the movement of the central star, as found in previous studies. Because destabilisation by a forced temperature gradient does not explicitly require disc self-gravity, we suggest this mechanism may also affect low-frequency one-armed oscillations in non-self-gravitating discs.
We investigate the dynamical evolution of a star cluster in an external tidal field by using N-body simulations, with focus on the effects of the presence or absence of neutron star natal velocity kicks.We show that, even if neutron stars typically represent less than 2% of the total bound mass of a star cluster, their primordial kinematic properties may affect the lifetime of the system by up to almost a factor of four. We interpret this result in the light of two known modes of star cluster dissolution, dominated by either early stellar evolution mass loss or two-body relaxation. The competition between these effects shapes the mass loss profile of star clusters, which may either dissolve abruptly ("jumping"), in the pre-core-collapse phase, or gradually ("skiing"), after having reached core collapse.
Cataclysmic Variables (CVs) are close binary star systems with one component an accreting white dwarf (WD) and the other a larger cooler star that fills its Roche Lobe. One consequence of the WDs accreting material, is the possibility that they are growing in mass and will eventually reach the Chandrasekhar Limit. This evolution could result in a Supernova Ia (SN Ia) explosion and is designated the Single Degenerate Progenitor (SD) scenario. One problem with the single degenerate scenario is that it is generally assumed that the accreting material mixes with WD core material at some time during the accretion phase of evolution and, since the typical WD has a carbon-oxygen (CO) core, the mixing results in large amounts of carbon and oxygen being brought up into the accreted layers. The presence of enriched carbon causes enhanced nuclear fusion and a Classical Nova (CN)explosion. Thus, the WD in a Classical Nova system is decreasing in mass and cannot be a SN Ia progenitor. In new calculations reported here, the consequences to the WD of no mixing of accreted material with core material have been investigated and the material involved in the explosion has only a Solar composition. I find that once sufficient material has been accreted, nuclear burning continues until a thermonuclear runaway (TNR) occurs and the WD either ejects a small amount of material or its radius grows to about $10^{12}$ cm and the calculations are stopped. In all cases where mass ejection occurs, the mass of the ejecta is far less than the mass of the accreted material. Therefore, all the WDs are growing in mass. (Abridged)
The physical nature of thermal composite supernova remnants (SNRs) remains controversial. We have revisited the archival XMM-Newton and Chandra data of the thermal composite SNR Kesteven 41 (Kes 41 or G337.8-0.1) and performed a millimeter observation toward this source in the $^{12}$CO, $^{13}$CO, and C$^{18}$O lines. The X-ray emission, mainly concentrated toward the southwestern part of the SNR, is characterized by distinct S and Ar He-like lines in the spectra. The X-ray spectra can be fitted with an absorbed nonequilibrium ionization collisional plasma model at a temperature of 1.3-2.6 keV and an ionization timescale of 0.1-1.2$\times$10$^{12}$ cm$^{-3}$ s. The metal species S and Ar are overabundant, with 1.2-2.7 and 1.3-3.8 solar abundances, respectively, which strongly indicate the presence of a substantial ejecta component in the X-ray-emitting plasma of this SNR. Kes 41 is found to be associated with a giant molecular cloud (MC) at a systemic local standard of rest velocity of -50 km s$^{-1}$ and confined in a cavity delineated by a northern molecular shell, a western concave MC that features a discernible shell, and an HI cloud seen toward the southeast of the SNR. The birth of the SNR in a preexisting molecular cavity implies a mass of $\gtrsim$18 M$_{\odot}$ for the progenitor if it was not in a binary system. Thermal conduction and cloudlet evaporation seem to be feasible mechanisms to interpret the X-ray thermal composite morphology, and the scenario of gas reheating by the shock reflected from the cavity wall is quantitatively consistent with the observations. An updated list of thermal composite SNRs is also presented in this paper.
We propose a novel scenario for the formation of Globular Clusters (GCs) based on the merger of two or more atomic cooling halos at high-redshift (z>6). The model naturally fulfills several key observational constraints on GCs that have emerged in the last decade. Specifically, absolute and relative ages, widespread presence of multiple stellar populations, spatial distribution around host galaxies, and correlations between galactocentric radius and metallicity. In our framework, the oldest globular clusters form the first generation stars as an intense burst in the center of a minihalo that grows above the threshold for hydrogen cooling (halo mass M_h~1e8 Msun) and undergoes a major merger within the cooling timescale (~150 Myr). Subsequent minor mergers and sustained gas infall bring new supply of pristine gas at the halo center, diluting AGB ejecta, and triggering additional bursts of star formation which form multiple generation of stars in the majority of the clusters. The DM halo around the GC is then stripped during assembly of the host galaxy halo. Our modeling is based on the merging history of dark-matter halos, and thus has no free adjustable parameters within the concordance LCDM cosmology. As a first application, based on a high-resolution cosmological simulation, we make quantitative predictions for the age distribution of the old GC population (Age=13.0+/-0.2 Gyr). We suggest that a similar merging mechanism is responsible for forming the sequence of younger and progressively metal richer clusters, through subhalo-subhalo merging in the later stages of massive halo assembly.
The Lick-index spectrophotometric system is investigated in its inherent statistical and operational properties to ease a more appropriate use for astrophysical studies. Non-Gaussian effects in the index standardization procedure suggest that a minimum S/N ratio has to be reached by spectral data, such as S/N >= 5 px^{-1} for a spectral resolution R~2000. In addition, index (re-)definition in terms of narrow-band "color" should be preferred over the classical pseudo-equivalent width scheme. The overlapping wavelength range among different indices is also an issue, as it may lead the latter ones to correlate, beyond any strictly physical relationship. The nested configuration of the Fe5335, Fe5270 indices, and the so-called "Mg complex" (including Mg1, Mg2 and Mgb) is analysed, in this context, by assessing the implied bias when joining entangled features into "global" diagnostic meta-indices, like the perused [MgFe] metallicity tracer. The perturbing effect of [OIII](5007) and [NI](5199) forbidden gas emission on Fe5015 and Mgb absorption features is considered, and an updated correction scheme is proposed when using [OIII](5007) as a proxy to appraise Hbeta residual emission. When applied to present-day elliptical galaxy population, the revised Hbeta scale leads, on average, to 20-30% younger age estimates. Finally, the misleading role of the christening element in Lick-based chemical analyses is illustrated for the striking case of Fe4531. In fact, while Iron is nominally the main contributor to the observed feature in high-resolution spectra, we have shown that the Fe4531 index actually maximizes its responsiveness to Titanium abundance.
Gamma Ray Bursts (GRBs) are a powerful probe of the high redshift Universe. We present a tool to estimate the detection rate of high-z GRBs by a generic detector with defined energy band and sensitivity. We base this on a population model that reproduces the observed properties of GRBs detected by Swift, Fermi and CGRO in the hard X-ray and gamma-ray bands. We provide the expected cumulative distributions of the flux and fluence of simulated GRBs in different energy bands. We show that scintillator detectors, operating at relatively high energies (e.g. tens of keV to the MeV), can detect only the most luminous GRBs at high redshifts due to the link between the peak spectral energy and the luminosity (Ep-Liso) of GRBs. We show that the best strategy for catching the largest number of high-z bursts is to go softer (e.g. in the soft X-ray band) but with a very high sensitivity. For instance, an imaging soft X-ray detector operating in the 0.2-5 keV energy band reaching a sensitivity, corresponding to a fluence of ~10^-8 erg cm^-2, is expected to detect ~40 GRBs yr^-1 sr^-1 at z>5 (~3 GRBs yr^-1 sr^-1 at z>10). Once high-z GRBs are detected the principal issue is to secure their redshift. To this aim we estimate their NIR afterglow flux at relatively early times and evaluate the effectiveness of following them up and construct usable samples of events with any forthcoming GRB mission dedicated to explore the high-z Universe.
Hadronic emission from supernova remnant (SNR)--molecular cloud (MC) association systems has been widely regarded as a probe of the shock accelerated cosmic protons. We here report on the detection of a $\gamma$-ray emission source, with a significance of $24\sigma$ in 0.2--300 GeV, projectively on the northwest of SNR Kesteven 41, using 5.6 yr Fermi-LAT observation data. The $3\sigma$ error circle, 0.09 degree in radius, covers the 1720MHz OH maser and is essentially consistent with the location of the $V_{LSR}~-50 km/s $ MC with which the SNR interacts. The source emission has a power-law spectrum with a photon index $2.38\pm0.03$ and a 0.2--300 GeV luminosity $~1.6*10^{36} erg /s$ at a distance 12 kpc. There is no radio pulsar in the $3\sigma$ circle responsible for the high luminosity. While the inverse Compton scattering scenario would lead to a difficulty in the electron energy budget, the source emission can be naturally explained with the hadronic interaction between the relativistic protons accelerated by the shock of SNR~Kesteven~41 and the adjacent northwestern MC.
The apparent correlation between the specific star formation rate (sSFR) and total stellar mass (M_star) of galaxies is a fundamental relationship indicating how they formed their stellar populations. To attempt to understand this relation, we hypothesize that the relation and its evolution is regulated by the increase in the stellar and gas mass surface density in galaxies with redshift, which is itself governed by the angular momentum of the accreted gas, the amount of available gas, and by self-regulation of star formation. With our model, we can reproduce the specific SFR-M_star relations at z~1-2 by assuming gas fractions and gas mass surface densities similar to those observed for z=1-2 galaxies. We further argue that it is the increasing angular momentum with cosmic time that causes a decrease in the surface density of accreted gas. The gas mass surface densities in galaxies are controlled by the centrifugal support (i.e., angular momentum), and the sSFR is predicted to increase as, sSFR(z)=(1+z)^3/t_H0, as observed (where t_H0 is the Hubble time and no free parameters are necessary). At z>~2, we argue that star formation is self-regulated by high pressures generated by the intense star formation itself. The star formation intensity must be high enough to either balance the hydrostatic pressure (a rather extreme assumption) or to generate high turbulent pressure in the molecular medium which maintains galaxies near the line of instability (i.e. Toomre Q~1). The most important factor is the increase in stellar and gas mass surface density with redshift, which allows distant galaxies to maintain high levels of sSFR. Without a strong feedback from massive stars, such galaxies would likely reach very high sSFR levels, have high star formation efficiencies, and because strong feedback drives outflows, ultimately have an excess of stellar baryons (abridged).
Solar energetic particles (SEPs) affect the solar-terrestrial space environment and become a very important aspect in space weather research. In this work, we numerically investigate the transport processes of SEPs in three-dimensional interplanetary magnetic field, with emphasis on the longitudinal distribution of SEPs in the heliosphere. We find that there exists an east-west longitudinal asymmetry in the SEP intensities, i.e., with the same longitude separations between the solar source centers and the magnetic footpoint of the observer, the fluxes of SEP events originating from solar sources located at eastern side of the nominal magnetic footpoint of observer are systematically larger than those of the SEP events originating from sources located at western side. In combination with the empirical results in previous works, we discuss the formation mechanism of this phenomenon, and propose that the longitudinally asymmetric distribution of SEPs results from the east-west azimuthal asymmetry in the topology of the Parker interplanetary magnetic field as well as the effects of perpendicular diffusion on the transport of SEPs in the heliosphere. Our results would be valuable in understanding the Sun-Earth relations and useful in space weather forecasting.
We report new laboratory studies of the radiation-induced destruction of glycine-containing ices for a range of temperatures and compositions that allow extrapolation to Martian conditions. In-situ infrared spectroscopy was used to study glycine decay rates as a function of temperature (from 15 to 280 K) and initial glycine concentrations in six mixtures whose compositions ranged from dry glycine to H$_2$O + glycine (300:1). Results are presented in several systems of units, with cautions concerning their use. The half-life of glycine under the surface of Mars is estimated as an extrapolation of this data set to Martian conditions, and trends in decay rates are described as are applications to Mars's near-surface chemistry.
Recent application of the Bayesian algorithm BORG to the Sloan Digital Sky Survey (SDSS) main sample galaxies resulted in the physical inference of the formation history of the observed large-scale structure from its origin to the present epoch. In this work, we use these inferences as inputs for a detailed probabilistic cosmic web-type analysis. To do so, we generate a large set of data-constrained realizations of the large-scale structure using a fast, fully non-linear gravitational model. We then perform a dynamic classification of the cosmic web into four distinct components (voids, sheets, filaments and clusters) on the basis of the tidal field. Our inference framework automatically and self-consistently propagates typical observational uncertainties to web-type classification. As a result, this study produces highly detailed and accurate cosmographic classification of large-scale structure elements in the SDSS volume. By also providing the history of these structure maps, the approach allows an analysis of the origin and growth of the early traces of the cosmic web present in the initial density field and of the evolution of global quantities such as the volume and mass filling fractions of different structures. For the problem of web-type classification, the results described in this work constitute the first connection between theory and observations at non-linear scales including a physical model of structure formation and the demonstrated capability of uncertainty quantification. A connection between cosmology and information theory using real data also naturally emerges from our probabilistic approach. Our results constitute quantitative chrono-cosmography of the complex web-like patterns underlying the observed galaxy distribution.
We report the fifth confirmed Luminous Blue Variable/S Doradus variable in M31. In 2006, J004526.62+415006.3 had the spectrum of hot Fe II emission line star with strong P Cygni profiles in the Balmer lines. In 2010, its absorption line spectrum resembled an early A-type supergiant with H and Fe II emission lines with strong P Cygni profiles, and in 2013 the spectrum had fully transitioned to an F-type supergiant due to the formation of the optically thick, cool wind which characterizes LBVs at maximum light. The photometric record supports the LBV/S Dor nature of the variability. Its bolometric luminosity ~ -9.65 mag places it on the HR Diagram near the known LBVs, AE And, Var C in M33 and S Dor.
We show that a self-consistent and coupled treatment of the weak decoupling, big bang nucleosynthesis, and photon decoupling epochs can be used to provide new insights and constraints on neutrino sector physics from high-precision measurements of light element abundances and cosmic microwave background observables. Implications of beyond-standard-model physics in cosmology, especially within the neutrino sector, are assessed by comparing predictions against five observables: the baryon energy density, helium abundance, deuterium abundance, effective number of neutrinos, and sum of the light neutrino mass eigenstates. We give examples for constraints on dark radiation, neutrino rest mass, lepton numbers, and scenarios for light and heavy sterile neutrinos.
We present new kinetics data on the radiolytic destruction of amino acids measured in situ with infrared spectroscopy. Samples were irradiated at 15, 100, and 140 K with 0.8-MeV protons, and amino-acid decay was followed at each temperature with and without H$_2$O present. Observed radiation products included CO$_2$ and amines, consistent with amino-acid decarboxylation. The half-lives of glycine, alanine, and phenylalanine were estimated for various extraterrestrial environments. Infrared spectral changes demonstrated the conversion from the non-zwitterion structure NH$_2$-CH$_2$(R)-COOH at 15 K to the zwitterion structure $^+$NH$_3$-CH$_2$(R)-COO$^-$ at 140 K for each amino acid studied.
The redshifted 21-cm background is expected to be a powerful probe of the early Universe, carrying both cosmological and astrophysical information from a wide range of redshifts. In particular, the power spectrum of fluctuations in the 21-cm brightness temperature is anisotropic due to the line-of-sight velocity gradient, which in principle allows for a simple extraction of this information in the limit of linear fluctuations. However, recent numerical studies suggest that the 21-cm signal is actually rather complex, and its analysis likely depends on detailed model fitting. We present the first realistic simulation of the anisotropic 21-cm power spectrum over a wide period of early cosmic history. We show that on observable scales, the anisotropy is large and thus measurable at most redshifts, and its form tracks the evolution of 21-cm fluctuations as they are produced early on by Lyman-a radiation from stars, then switch to X-ray radiation from early heating sources, and finally to ionizing radiation from stars. In particular, we predict a redshift window during cosmic heating (at z ~ 15), when the anisotropy is small, during which the shape of the 21-cm power spectrum on large scales is determined directly by the average radial distribution of the flux from X-ray sources. This makes possible a model-independent reconstruction of the X-ray spectrum of the earliest sources of cosmic heating.
We identify three isotopic tracers that can be used to constrain the $^{13}C$-pocket and show the correlated isotopic ratios of Sr and Ba in single mainstream presolar SiC grains. These newly measured data can be explained by postprocess AGB model calculations with large $^{13}C$-pockets with a range of relatively low $^{13}C$ concentrations, which may suggest that multiple mixing processes contributed to the $^{13}C$-pocket formation in parent AGB stars.
We present the design and measurements on a 90GHz prototype of a millimeter-wave channelizing spectrometer realized in rectangular waveguide for astronomical instrumentation. The device was fabricated using conventional high-precision metal machining, and the spectrometer can be tiled into a 2D array to fill the focal plane of a telescope. Measurements of the fabricated five-channel device matched well with electromagnetic simulations using HFSS and a cascaded S-matrix approach. This motivated the design of a 54-channel R=200 spectrometer that fills the single-moded passband of rectangular waveguide in the 130-175 GHz and 190-250 GHz atmospheric windows for millimeter-wave spectroscopic mapping and multi-object spectroscopy.
We use a Cartesian grid to simulate the flow of gas in a barred Galactic potential and investigate the effects of varying the sound speed in the gas and the resolution of the grid. For all sound speeds and resolutions, streamlines closely follow closed orbits at large and small radii. At intermediate radii shocks arise and the streamlines shift between two families of closed orbits. The point at which the shocks appear and the streamlines shift between orbit families depends strongly on sound speed and resolution. For sufficiently large values of these two parameters, the transfer happens at the cusped orbit as hypothesised by Binney et al. over two decades ago. For sufficiently high resolutions the flow downstream of the shocks becomes unsteady. If this unsteadiness is physical, as appears to be the case, it provides a promising explanation for the asymmetry in the observed distribution of CO.
We build on previous efforts to model CCD sensors, during illumination and collection of conversions. We use a finite summation of simple, electrostatic field models. The upgraded functionality of our framework provides specific predictions for perturbations in pixel boundary enclosures (e.g., at the backside window) and the bookkeeping capability to stack those perturbations so that they may be utilized as Greens functions -- portable calculation results that may be generically applied to a range of precision astronomy related problems that naturally including astrometric, photometric and shape transfer issues. We approach the topic of using ancillary pixel data, derived from the Greens function and the registered image, to analyze sky data and constrain object parameters of astronomical targets.
The Crab pulsar is the only astronomical pulsed source detected above 100 GeV. The emission mechanism of very high energy gamma-ray pulsation is not yet fully understood, although several theoretical models have been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 4.6 years of Fermi-LAT data were also analyzed. The known two pulses per period were detected with a significance of 8.0 sigma and 12.6 sigma. In addition, significant bridge emission was found between the two pulses with 6.2 sigma. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models.
We investigate properties of Galactic microlensing events in which a stellar object is lensed by a neutron star. For an all-sky photometric microlensing survey, we determine the number of lensing events caused by $\sim10^{5}$ potentially-observable radio pulsars to be $\sim0.2\ \rm{yr^{-1}}$ for $10^{10}$ background stellar sources. We expect a few detectable events per year for the same number of background sources from an astrometric microlensing survey. We show that such a study could lead to precise measurements of radio pulsar masses. For instance, if a pulsar distance could be constrained through radio observations, then its mass would be determined with a precision of $\sim10\%$. We also investigate the time-scale distributions for neutron star events, finding that they are much shorter than had been previously thought. For photometric events towards the Galactic centre that last $\sim15$ days, around $7\%$ will have a neutron star lens. This fraction drops rapidly for longer time-scales. Away from the bulge region we find that neutron stars will contribute $\sim40\%$ of the events that last less than $\sim10$ days. These results are in contrast to earlier work which found that the maximum fraction of neutron star events would occur on time-scales of hundreds of days.
In Kang (2015) we calculated the acceleration of cosmic-ray electrons and the ensuing radio synchrotron emission at weak spherical shocks that are expected to form in the outskirts of galaxy clusters.There we demonstrated that, at decelerating spherical shocks, the volume integrated spectra of both electrons and radiation deviate significantly from the test-particle power-laws predicted for constant planar shocks, because the shock compression ratio and the flux of injected electrons decrease in time. In this study, we consider spherical blast waves propagating into a constant density core surrounded by an isothermal halo with a decreasing density profile in order to explore how the deceleration rate of the shock speed affects the radio emission from accelerated electrons. The surface brightness profile and the volume-integrated radio spectrum of the model shocks are calculated by assuming a ribbon-like shock surface on a spherical shell and the associated downstream region of relativistic electrons. If the postshock magnetic field strength is about 7 microgauss, at the shock age of ~50 Myr, the volume-integrated radio spectrum steepens gradually with the spectral index from alpha_{inj} to alpha_{inj}+0.5 over 0.1-10 GHz, where alpha_{inj} is the injection index at the shock positionexpected from the diffusive shock acceleration theory. Such gradual steepening could explain the curved radio spectrum of the radio relic in cluster A2266, which was interpreted as a broken power-law by Trasatti et al. (2014), if the relic shock is young enough so that the break frequency falls in around 1 GHz.
We describe Rabacus, a Python package for calculating the transfer of hydrogen ionizing radiation in simplified geometries relevant to astronomy and cosmology. We present example solutions for three specific cases: 1) a semi-infinite slab gas distribution in a homogeneous isotropic background, 2) a spherically symmetric gas distribution with a point source at the center, and 3) a spherically symmetric gas distribution in a homogeneous isotropic background. All problems can accommodate arbitrary spectra and density profiles as input. The solutions include a treatment of both hydrogen and helium, a self-consistent calculation of equilibrium temperatures, and the transfer of recombination radiation. The core routines are written in Fortran 90 and then wrapped in Python leading to execution speeds thousands of times faster than equivalent routines written in pure Python. In addition, all variables have associated units for ease of analysis. The software is part of the Python Package Index and the source code is available on Bitbucket at https://bitbucket.org/galtay/rabacus. In addition, installation instructions and a detailed users guide are available at this http URL
Blazars are thought to possess a relativistic jet that is pointing toward the
direction of the Earth and the elect of relativistic beaming enhances its
apparent brightness. They radiate in all wavebands from the radio to the
gamma-ray bands via the synchrotron and the inverse Compton scattering process.
Numerous observations are performed but the mechanism of variability, creation
and composition of jets are still controversial.
We performed multi-wavelength monitoring with optical polarization for 3C
66A, Mrk 421, CTA 102 and PMN J0948+0022 to investigate the mechanisms of
variability and research the emission region in the relativistic jets.
Consequently, an emergence of new emission component in flaring state is
suggested in each object. The most significant aspect of these results is its
wide range of sizes of emission regions from $10^{14}-10^{16}$ cm, which
implies the model with a number of independent emission regions with variety
sizes and randomly orientation. The "shock-in-jet" scenario can explain high PD
and direction of PA in each objects. It might reflect the common mechanism of
flares in the relativistic jets.
Recently, a spatially extended excess of gamma rays collected by the Fermi-LAT from the inner region of the Milky Way has been detected by different groups and with increasingly sophisticated techniques. Yet, any final conclusion about the morphology and spectral properties of such an extended diffuse emission are subject to a number of potentially critical uncertainties, related to the high density of cosmic rays, gas, magnetic fields and abundance of point sources. We will present a thorough study of the systematic uncertainties related to the modelling of diffuse background and to the propagation of cosmic rays in the inner part of our Galaxy. We will test a large set of models for the Galactic diffuse emission, generated by varying the propagation parameters within extreme conditions. By using those models in the analysis of Fermi-LAT data as Galactic foreground, we will show that the gamma-ray excess survives and we will quantify the uncertainties affecting the excess morphology and energy spectrum.
Astronomical photometry is the science of measuring the flux of a celestial object. Since its introduction, the CCD has been the principle method of measuring flux to calculate the apparent magnitude of an object. Each CCD image taken must go through a process of cleaning and calibration prior to its use. As the number of research telescopes increases the overall computing resources required for image processing also increases. Existing processing techniques are primarily sequential in nature, requiring increasingly powerful servers, faster disks and faster networks to process data. Existing High Performance Computing solutions involving high capacity data centres are complex in design and expensive to maintain, while providing resources primarily to high profile science projects. This research describes three distributed pipeline architectures, a virtualised cloud based IRAF, the Astronomical Compute Node (ACN), a private cloud based pipeline, and NIMBUS, a globally distributed system. The ACN pipeline processed data at a rate of 4 Terabytes per day demonstrating data compression and upload to a central cloud storage service at a rate faster than data generation. The primary contribution of this research is NIMBUS, which is rapidly scalable, resilient to failure and capable of processing CCD image data at a rate of hundreds of Terabytes per day. This pipeline is implemented using a decentralised web queue to control the compression of data, uploading of data to distributed web servers, and creating web messages to identify the location of the data. Using distributed web queue messages, images are downloaded by computing resources distributed around the globe. Rigorous experimental evidence is presented verifying the horizontal scalability of the system which has demonstrated a processing rate of 192 Terabytes per day with clear indications that higher processing rates are possible.
Fermi-LAT reveals two huge gamma-ray bubbles existing in the Galactic Center, called 'Fermi Bubbles'. The existence of two microwave bubbles at the same region are also reported by the observation by WMAP, dubbed 'WMAP haze'. In order to explain these components, It has been argued that the gamma-rays arise from Inverse-Compton scattering of relativistic electrons accelerated by plasma turbulence, and the microwaves are radiated by synchrotron radiation. But no previous research reproduces both the Fermi Bubbles and WMAP haze under typical magnetic fields in the galaxy. We assume that shocks present in the bubbles and the efficiency of the acceleration by plasma turbulence, 'stochastic acceleration', changes with the distance from the shock front. The distance from the shock front increases with time, accordingly the efficiency of the acceleration changes with time. We also consider the time development of the electrons escape from the turbulence by diffusive loss. Our model succeed to reproduce both the observed characteristics of the Fermi Bubbles and WMAP haze under typical magnetic fields.
We performed simultaneous observations in 3 bands (UBV) of the flickering variability of the recurrent novae RS Oph and T CrB at quiescence. Using new and published data, we compare the colours of the flickering in cataclysmic variables and symbiotic recurrent novae. We find a difference between the colours of the flickering source in these two types of accreting white dwarfs. The detected difference is highly significant with $ p-value \approx 2 \times 10^{-6}$ for the distributions of $(U-B)_0$ colour and $p \approx 3 \times 10^{-5}$ on (U-B) versus (B-V) diagram. The possible physical reasons are briefly discussed. The data are available upon request from the authors.
In this letter we adopt a new method of estimating Hubble constant, which combines observation of supernovae and model dependent result from the latest report of Planck project. The process of our estimation is neither pure bayesian nor pure frequent, so the result can also acts as null hypothesis test of standard cosmological model and consistency check of the associated data-sets. We mainly report the Hubble constant is estimated as $h = 0.700\pm 0.118$, according to Union$2.1$ data-set and matter density estimated from CMB observation. We find no trend of deviation from standard cosmological model, and estimate the sensitivity of our method for future observation.
We review the current understanding of the diffuse gamma-ray background (DGRB). The DGRB is what remains of the total measured gamma-ray emission after the subtraction of the resolved sources and of the diffuse Galactic foregrounds. It is interpreted as the cumulative emission of sources that are not bright enough to be detected individually. Yet, its exact composition remains unveiled. Well-established astrophysical source populations (e.g. blazars, misaligned AGNs, star-forming galaxies and millisecond pulsars) all represent guaranteed contributors to the DGRB. More exotic scenarios, such as dark matter annihilation or decay, may contribute as well. In this review, we describe how these components have been modeled in the literature and how the DGRB can be used to provide valuable information on each of them. We summarize the observational information currently available on the DGRB, paying particular attention to the most recent measurement of its intensity energy spectrum by the Fermi LAT Collaboration. We also discuss the novel analyses of the auto-correlation angular power spectrum of the DGRB and of its cross-correlation with tracers of the large-scale structure of the Universe. New data sets already (or soon) available are expected to provide further insight on the nature of this emission. By summarizing where we stand on the current knowledge of the DGRB, this review is intended both as a useful reference for those interested in the topic and as a means to trigger new ideas for further research.
We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by $\sim60\,000$ secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy lensing and the stellar mass function. We interpret our measurements within the halo occupation distribution (HOD) framework, separating the contributions from central and satellite galaxies. We find that the SHMR for the central galaxies peaks at $M_{\rm h, peak} = 1.9^{+0.2}_{-0.1}\times10^{12} M_{\odot}$ with an amplitude of $0.025$, which decreases to $\sim0.001$ for massive halos ($M_{\rm h} > 10^{14} M_{\odot}$). Compared to central galaxies only, the total SHMR (including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size halos), a result consistent with cluster analyses from the literature based on fully independent methods. After properly accounting for differences in modelling, we have compared our results with a large number of results from the literature up to $z=1$: we find good general agreement, independently of the method used, within the typical stellar-mass systematic errors at low to intermediate mass (${M}_{\star} < 10^{11} M_{\odot}$) and the statistical errors above. We have also compared our SHMR results to semi-analytic simulations and found that the SHMR is tilted compared to our measurements in such a way that they over- (under-) predict star formation efficiency in central (satellite) galaxies.
The jet-counterjet system of the closest radio-loud active galaxy Centaurus A (Cen A) can be studied with Very Long Baseline Interferometry (VLBI) on unprecedented small linear scales of ~0.018pc. These high-resolution observations provide essential information on jet emission and propagation within the inner parsec of an AGN jet. We present the results of a kinematic study performed within the framework of the Southern-hemisphere AGN monitoring program TANAMI. Over 3.5years, the evolution of the central-parsec jet structure of Cen A was monitored with VLBI. These observations reveal complex jet dynamics which are well explained by a spine-sheath structure supported by the downstream acceleration occurring where the jet becomes optically thin. Both moving and stationary jet features are tracked. A persistent local minimum in surface brightness suggests the presence of an obstacle interrupting the jet flow, which can be explained by the interaction of the jet with a star at a distance of ~0.4pc from the central black hole. We briefly discuss possible implications of such an interaction regarding the expected neutrino and high-energy emission and the effect on a putative planet.
There is a maximum for the gravity of a black hole in the vertical direction in the accretion disc. Outflows may probably be driven from the disc if the radiation flux of the disc is greater than a critical value corresponding to the maximal vertical gravity. We find that outflows are driven by the radiation force from the disc if the accretion rate is greater than the Eddington rate. The radiation of the disc is therefore limited by such outflows. The disc luminosity, L=L_Edd\propto ln mdot, at large-mdot cases. The Eddington ratio of the disc is ~3 for mdot~100, which is significantly lower than that of a conventional slim disc without outflows. This implies that the emission from some ultra-luminous X-ray sources with highly super Eddington luminosity should be Doppler beamed, or intermediate mass black holes are in these sources instead of stellar mass black holes. The spectra of the discs with outflows are saturated in the high frequency end provided mdot>2. We suggest that the saturated emission can be observed to estimate the masses of the black holes accreting at high rates, such as the narrow-line Seyfert galaxies, with the model calculations. The rate of the mass accreted by the black hole is always around the Eddington rate even if the mass accretion rate at the outer radius is very high, because most of the gas is removed into the outflows. This implies that the luminous quasars at high redshifts z>6 should have grown up through persistent accretion at a rate close to the Eddington rate.
The spectrum of the recently discovered cataclysmic variable star (CV) ASAS-SN 13cl shows that a secondary star with spectral type K4 (+- 2 subclasses) contributes roughly half the optical light. The radial velocities of the secondary are modulated on an orbital period P_orb = 4.86 hr with a velocity semiamplitude K = 246 +- 9 km/s, and the light curve shows ellipsoidal variations and an apparent grazing eclipse. At this orbital period, the secondary stars in most CVs are substantially cooler, with spectral types near M3. ASN-13cl therefore joins the small group of CVs with anomalously warm secondary stars, which apparently form when the onset of mass transfer occurs after the secondary has undergone significant nuclear evolution.
Due to its extraordinarily high concentration of known relativistic particle accelerators such as pulsar wind nebula, supernova remnants, dense molecular cloud regions, and the supermassive black hole (Sgr A*); the center of the Milky Way galaxy has long been an ideal target for high energy (HE, 0.1-100 GeV) and very high energy ( VHE, 50 GeV-50 TeV) gamma-ray emission. Indeed, detections of Sgr A* and other nearby regions of gamma-ray emission have been reported by EGRET and Fermi-LAT in the HE band, as well as CANGAROO, Whipple, HESS, VERITAS, and MAGIC in the VHE band. Here we report on the results of extended observations of the region with VERITAS between 2010-2014. Due to the visibility of the source for VERITAS in the Northern Hemisphere, these observations provide the most sensitive probe of gamma-ray emission above 2 TeV in one of the most complicated and interesting regions of our home galaxy.
The number of extragalactic sources detected at very hight energy (VHE, E$>$100GeV) has dramatically increased during the past years to reach more than fifty. The High Energy Stereoscopic System (H.E.S.S.) had observed the sky for more than 10 years now and discovered about twenty objects. With the advent of the fifth 28 meters telescope, the H.E.S.S. energy range extends down to ~30 GeV. When H.E.S.S. data are combined with the data of the Fermi Large area Telescope, the covered energy range is of several decades allowing an unprecedented description of the spectrum of extragalactic objects. In this talk, a review of the extragalactic sources studied with H.E.S.S. will be given together with first H.E.S.S. phase II results on extragalactic sources.
The blazar PKS~2155-304 was the target of a multiwavelength campaign from June to October 2013 which widely improves our knowledge of its spectral energy distribution. This campaign involved the NuSTAR satellite (3-79 keV), the Fermi Large Area Telescope (LAT, 100~MeV-300~GeV) and the High Energy Stereoscopic System (H.E.S.S.) array phase II (with an energy threshold of few tens of GeV). While the observations with NuSTAR extend the X-ray spectrum to higher energies than before, H.E.S.S. phase II, together with the use of the LAT PASS 8, enhance the coverage of the $\gamma$-ray regime with an unprecedented precision. In this work, preliminary results from the multi-wavelength analysis are presented.
We present a new method to construct fully self-consistent equilibrium models of multi-component disc galaxies similar to the Milky Way. We define distribution functions for the stellar disc and dark halo that depend on phase space position only through action coordinates. We then use an iterative approach to find the corresponding gravitational potential. We study the adiabatic response of the initially spherical dark halo to the introduction of the baryonic component and find that the halo flattens in its inner regions with final minor-major axis ratios $q$ = 0.75 - 0.95. The extent of the flattening depends on the velocity structure of the halo particles with radially biased models exhibiting a stronger response. In this latter case, which is according to cosmological simulations the most likely one, the new density structure resembles a "dark disc" superimposed on a spherical halo. We discuss the implications of these results for our recent estimate of the local dark matter density. The velocity distribution of the dark-matter particles near the Sun is very non-Gaussian. All three principal velocity dispersions are boosted as the halo contracts, and at low velocities a plateau develops in the distribution of $v_z$. For models similar to a state-of-the-art Galaxy model we find velocity dispersions around 155 km s$^{-1}$ for $v_z$ and the tangential velocity, $v_\varphi$, and 140 - 175 km s$^{-1}$ for the in-plane radial velocity, $v_R$, depending on the anisotropy of the model.
Clusters are excellent test benches for verification and improvement of stellar evolution theory. The recent detection of solar-like oscillations in G-K giants in the open cluster NGC6819 with Kepler provides us with independent constraints on the masses and radii of stars on the red giant branch, as well as on the distance to clusters and their ages. We present, for NGC6819, evolutionary models by considering rotation-induced mixing ; and the theoretical low-l frequencies of our stellar models.
The high frequency peaked BL Lac object PG 1553+113 underwent a flaring event in 2012. The High Energy Stereoscopic System (H.E.S.S.) observed this source for two consecutive nights at very high energies (VHE, $E>$100~GeV). The data show an increase of a factor of three of the flux with respect to archival measurements with the same instrument and hints of intra-night variability. The data set has been used to put constraints on possible Lorentz invariance violation (LIV), manifesting itself as an energy dependence of the velocity of light in vacuum, and to set limits on the energy scale at which Quantum Gravity effects causing LIV may arise. With a new method to combine H.E.S.S. and Fermi large area telescope data, the previously poorly known redshift of PG 1555+113 has been determined to be close to the value derived from optical measurements.
We analyze the high-energy neutrino events observed by IceCube, aiming to probe the initial flavor of cosmic neutrinos. We study the track-to-shower ratio of the subset with energy above 60 TeV, where the signal is expected to dominate and show that different production mechanisms give rise to different predictions even accounting for the uncertainties due to neutrino oscillations. We include for the first time the passing muons observed by IceCube in the analysis. They corroborate the hypotheses that cosmic neutrinos have been seen and their flavor matches expectations.
As part of a national scientific network 'Pathways to Habitability' the formation of planets and the delivery of water onto these planets is a key question as water is essential for the development of life. In the first part of the paper we summarize the state of the art of planet formation - which is still under debate in the astronomical community - before we show our results on this topic. The outcome of our numerical simulations depends a lot on the choice of the initial distribution of planetesimals and planetary embryos after gas disappeared in the protoplanetary disk. We also take into account that some of these planetesimals of sizes in the order of the mass of the Moon already contained water; the quantity depends on the distance from the Sun - close-by bodies are dry, but starting from a distance of about 2 AU they can contain substantial amounts of water. We assume that the gas giants and terrestrial planets are already formed when we check the collisions of the small bodies containing water (in the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the crust and on the surface of Earth. In the second part we discuss in more detail how the formation of larger bodies after a collision may happen as the outcome depends on parameters like collision velocity, impact angle, and the materials involved. We present results obtained by SPH (Smooth Particle Hydrodynamics) simulations. We briefly describe this method and show different scenarios with respect to the formed bodies, possible fragmentation and the water content before and after the collision. In an appendix we discuss detection methods for extrasolar planets (close to 2000 such objects have been discovered so far).
The log-normal distribution represents the probability of finding randomly distributed particles in a micro canonical ensemble with high entropy. To a first approximation, a modified form of this distribution with a truncated termination may represent an isolated galactic disk, and this disk density distribution model was therefore run to give the best fit to the observational rotation curves for 37 representative galaxies. The resultant curves closely matched the observational data for a wide range of velocity profiles and galaxy types with rising, flat or descending curves in agreement with Verheijen's classification of 'R', 'F' and 'D' type curves, and the corresponding theoretical total disk masses could be fitted to a baryonic Tully Fisher relation (bTFR). Nine of the galaxies were matched to galaxies with previously published masses, suggesting a mean excess dynamic disk mass of dex0.61+/-0.26 over the baryonic masses. Although questionable with regard to other measurements of the shape of disk galaxy gravitational potentials, this model can accommodate a scenario in which the gravitational mass distribution, as measured via the rotation curve, is confined to a thin plane without requiring a dark-matter halo or the use of MOND.
In the solar neighbourhood, there are moving groups of stars with similar ages and others of stars with heterogeneous ages as the field stars. To explain these facts, we have constructed a simple model of three phases. Phase A: a giant interstellar cloud is uniformly accelerated (or decelerated) with respect to the field stars during a relatively short period of time (10 Myr) and the cloud's mass is uniformly increased; phase B: the acceleration (or deceleration) and mass accretion of the cloud cease. The star formation spreads throughout the cloud, giving origin to stellar groups of similar ages; and phase C: the cloud loses all its gaseous component at a constant rate and in parallel is uniformly decelerated (or accelerated) until reaching the initial velocity of phase A (case 1) or the velocity of the gas cloud remains constant (case 2). Both cases give equivalent results. The system equations for the star motions governed by a time-dependent gravitational potential of the giant cloud and referred to a coordinate system co-moving with the cloud have been solved analytically. We have assumed a homogeneous spheroidal cloud of fixed semi-major axis a=300 pc and of an initial density of 7 at cm^{-3}, with a density increment of 100 per cent and a cloud's velocity variation of 30 km s^{-1}, from the beginning to the end of Phase A. The result is that about 4 per cent of the field stars that are passing within the volume of the cloud at the beginning of phase A is captured. The Sun itself could have been captured by the same cloud that originated the moving groups of the solar neighbourhood.
GAMMA-400 is a new space mission which will be installed on board the Russian space platform Navigator. It is scheduled to be launched at the beginning of the next decade. GAMMA-400 is designed to study simultaneously gamma rays (up to 3 TeV) and cosmic rays (electrons and positrons from 1 GeV to 20 TeV, nuclei up to 10$^{15}$-10$^{16}$ eV). Being a dual-purpose mission, GAMMA-400 will be able to address some of the most impelling science topics, such as search for signatures of dark matter, cosmic-rays origin and propagation, and the nature of transients. GAMMA-400 will try to solve the unanswered questions on these topics by high-precision measurements of the Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission and the spectra of cosmic-ray electrons + positrons and nuclei, thanks to excellent energy and angular resolutions.
We describe the microphysics, phenomenology, and astrophysical implication of a $B$-field induced unpairing effect that may occur in magnetars, if the local $B$-field in the core of a magnetar exceeds a critical value $H_{c2}$. Using the Ginzburg-Landau (GL) theory of superconductivity, we derive the $H_{c2}$ field for proton condensate taking into the correction ($\le 30\%$) which arises from its coupling to the background neutron condensate. The density dependence of pairing of proton condensate implies that $H_{c2}$ is maximal at the crust-core interface and decreases towards the center of the star. As a consequence, magnetar cores with homogenous constant fields will be partially superconducting for "medium-field" magnetars ($10^{15}\le B\le 5 \times 10^{16}$ G) whereas "strong-field" magnetars ($B>5\times 10^{16}$ G) will be void of superconductivity. The emissivity of a magnetar's core changes in twofold manner: (i)~the $B$-field assisted direct Urca process is enhanced by orders of magnitude, because of unpairing effect in regions where $B\ge H_{c2}$; (ii)~the Cooper-pair breaking processes on protons vanish in these regions and the overall emissivity by the pair-breaking processes is reduced by a factor of few.
We study the temporal evolution of downflows observed at the lateral edges of penumbral filaments in a sunspot located very close to the disk center. Our analysis is based on a sequence of nearly diffraction-limited spectropolarimetric scans of the Fe I 6173 A line taken with the CRISP instrument at the Swedish 1-m Solar Telescope. We compute Dopplergrams from the observed intensity profiles using line bisectors and filter the resulting velocity maps for subsonic oscillations. Lateral downflows appear everywhere in the center-side penumbra as small, weak patches of redshifts next to or along the edges of blueshifted flow channels. These patches have an intermittent life and undergo mergings and fragmentations quite frequently. The lateral downflows move together with the hosting filaments and react to their shape variations, very much resembling the evolution of granular convection in the quiet Sun. There is a good relation between brightness and velocity of the flow structures in the center-side penumbra, with downflows being darker than upflows on average, which is again reminiscent of convection in the quiet Sun. These results point to the existence of overturning convection in sunspot penumbrae, with elongated cells in the radial direction where the flow is upward but very inclined, and weak lateral downward flows. In general, the circular polarization profiles emerging from the lateral downflows do not show sign reversals, although sometimes we detect three-lobed profiles which are suggestive of opposite magnetic polarities in the pixel.
Gamma-ray blazars are among the most extreme astrophysical sources, harboring phenomena far more energetic than those attainable by terrestrial accelerators. These galaxies are understood to be active galactic nuclei that are powered by accretion onto supermassive black holes and have relativistic jets pointed along the Earth line of sight. The emission displayed is variable at all wavelengths and timescales probed thus far, necessitating contemporaneous broadband observations to disentangle the details of the emission processes within the relativistic jets. The very high energy (VHE; $E\ge$100 GeV) photons emitted by these sources are detectable with ground-based imaging atmospheric Cherenkov telescopes such as VERITAS. As these photons propagate extragalactic distances, the interaction with the diffuse starlight that pervades the entire Universe results in a distance and energy dependent gamma-ray opacity, offering a unique method for probing photon densities on cosmological scales. These galaxies have also been postulated to be potential sources of ultra-high-energy cosmic rays, a theory which can be examined through deep gamma-ray observations of sources which probe moderate gamma-ray opacities. Within this work, I will highlight ongoing research regarding the broadband emission from VERITAS-observed VHE blazars, as well as the potential to use them for cosmological insight.
One of the most important problem of the blazar astrophysics is to understand the physical origin of the blazar sequence. In this study, we focus on the GeV gamma-ray variability of blazars and evolution perspective we search the relation between the redshift and the variability amplitude of blazars for each blazar subclass. We analyzed the Fermi-LAT data of the TeV blazars and the bright AGNs (flux $\geq$ 4$\times10^{-9}$ cm$^{-2}$s$^{-1}$) selected from the 2LAC (the 2nd LAT AGN catalog) data base. As a result, we found a hint of the correlation between the redshift and the variability amplitude in the FSRQs. Furthermore the BL Lacs which have relatively lower peak frequency of the synchrotron radiation and relatively lower redshift, have a tendency to have a smaller variability amplitude.
The H.E.S.S. Galactic plane scan has revealed a large population of Galactic very high energy (VHE; E > 100 GeV) emitters. The majority of the galactic sources are extended and can typically be associated with pulsar wind nebulae (35%) and supernova remnants (21%), while some of the sources remain unidentified (31%). A much smaller fraction of point-like sources (5 in total, corresponding to 4%) are identified as gamma-ray binaries. Active galactic nuclei located behind the Galactic plane are also a potential source class. An active galaxy could be identified in the VHE regime by a point like extension, a high variability amplitude (up to a factor of 100) and a typically soft spectrum (due to absorption by the extra-galactic background light). Here we report on VERITAS observations of HESS J1943+213, an unidentified point source discovered to emit above 470 GeV during the extended H.E.S.S. Galactic plane scan. This source is thought to be a distant BL Lac object behind the Galactic plane and, though it exhibits a steep spectrum it is a weak GeV source, only recently detected using 5 years of Fermi-LAT data. Deep VERITAS observations at high elevations result in the most significant VHE detection of this object so far, with an excess above 200 GeV of more than 18 standard deviations. We use variability and spectral analyses of VERITAS data on HESS J1943+213 in a multi-wavelength context to address the source classification.
The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.
We present deep ($>$2.4 Ms) observations of the Cassiopeia A supernova remnant with {\it NuSTAR}, which operates in the 3--79 keV bandpass and is the first instrument capable of spatially resolving the remnant above 15 keV. We find that the emission is not entirely dominated by the forward shock nor by a smooth "bright ring" at the reverse shock. Instead we find that the $>$15 keV emission is dominated by knots near the center of the remnant and dimmer filaments near the remnant's outer rim. These regions are fit with unbroken power-laws in the 15--50 keV bandpass, though the central knots have a steeper ($\Gamma \sim -3.35$) spectrum than the outer filaments ($\Gamma \sim -3.06$). We argue this difference implies that the central knots are located in the 3-D interior of the remnant rather than at the outer rim of the remnant and seen in the center due to projection effects. The morphology of $>$15 keV emission does not follow that of the radio emission nor that of the low energy ($<$12 keV) X-rays, leaving the origin of the $>$15 keV emission as an open mystery. Even at the forward shock front we find less steepening of the spectrum than expected from an exponentially cut off electron distribution with a single cutoff energy. Finally, we find that the GeV emission is not associated with the bright features in the {\it NuSTAR} band while the TeV emission may be, suggesting that both hadronic and leptonic emission mechanisms may be at work.
The fraction of positrons and electrons in cosmic rays recently observed on the International Space Station unveiled an unexpected excess of the positrons, undermining the current foundations of cosmic rays sources. We provide a quantum electrodynamics phenomenological model explaining the observed data. This model incorporates electroproduction, in which cosmic ray electrons decelerating in the interstellar medium emit photons that turn into electron-positron pairs. These findings not only advance our knowledge of cosmic ray physics, but also pave the way for computationally efficient formulations of quantum electrodynamics, critically needed in physics and chemistry.
Accretion and ejection are tightly connected and represent the fundamental mechanisms regulating star formation. However, the exact physical processes involved are not yet fully understood. We present high angular and spectral resolution observations of the Br Gamma emitting region in the Herbig Ae star HD163296 (MWC275) in order to probe the origin of this line and constrain the physical processes taking place at sub-AU scales in the circumstellar region. By means of VLTI-AMBER observations at high spectral resolution (R~12000), we studied interferometric visibilities, wavelength-differential phases, and closure phases across the Br Gamma line of HD163296. To constrain the physical origin of the Br Gamma line in Herbig Ae stars, all the interferometric observables were compared with the predictions of a line radiative transfer disc wind model. The measured visibilities clearly increase within the Br Gamma line, indicating that the Br Gamma emitting region is more compact than the continuum. By fitting a geometric Gaussian model to the continuum-corrected Br Gamma visibilities, we derived a compact radius of the Br Gamma emitting region of ~0.07+/-0.02AU (Gaussian half width at half maximum; or a ring-fit radius of ~0.08+/-0.02AU). To interpret the observations, we developed a magneto-centrifugally driven disc wind model. Our best disc wind model is able to reproduce, within the errors, all the interferometric observables and it predicts a launching region with an outer radius of ~0.04AU. However, the intensity distribution of the entire disc wind emitting region extends up to ~0.16AU. Our observations, along with a detailed modelling of the Br Gamma emitting region, suggest that most of the Br Gamma emission in HD163296 originates from a disc wind with a launching region that is over five times more compact than previous estimates of the continuum dust rim radius.
The gamma-ray binary LS I +61{\deg}303 shows a discontinuity of the periodicity in its GeV emission. In this paper, we show that during the epochs when the timing analysis fails to determine the orbital periodicity, the periodicity is in fact present in the two orbital phase intervals $\Phi = 0.0-0.5$ and $\Phi = 0.5-1.0$. That is, there are two periodic signals, one towards periastron (i.e., $\Phi = 0.0-0.5$) and another one towards apastron ($\Phi = 0.5-1.0$). The apastron peak shows the same orbital shift as the radio outburst and, in addition, reveals the same two periods $P_1$ and $P_2$ that are present in the radio data. The gamma-ray emission of the apastron peak normally just broadens the emission of the peak around periastron. Only when it appears at $\Phi = 0.8-1.0$, because of the orbital shift, it is detached enough from the first peak to become recognizable as a second orbital peak, which is the reason why the timing analysis fails. Two gamma-ray peaks along the orbit are predicted by the two-peak accretion model for an eccentric orbit that was proposed by several authors for LS I +61{\deg}303.
The Cygnus Loop is a nearby supernova remnant (SNR) observed across the electromagnetic spectrum. With the analysis of 6 years of Fermi/LAT data we find that, what previous studies had considered a single source, consists of an extended source plus a point-like source south-east of the SNR. The extended gamma-ray emission is well correlated with the thermal X-ray emission of the SNR, and the energy spectrum displays a pronounced maximum at $\sim0.6$\,GeV. However, in a region where the radio emission is strongly and distinctly polarized, the gamma-ray spectrum shows no sign of a break. Therefore, the spatially resolved gamma-ray emission permits the study of different interaction conditions of the SNR and the surrounding medium.
One of the most important unresolved issues in gamma-ray burst physics is the origin of the prompt gamma-ray spectrum. Its general non-thermal character and the softness in the X-ray band remain unexplained. We tackle these issues by performing Monte Carlo simulations of radiation-matter interactions in a scattering dominated photon-lepton plasma. The plasma -- initially in equilibrium -- is driven to non-equilibrium conditions by a sudden energy injection in the lepton population, mimicking the effect of a shock wave or the dissipation of magnetic energy. Equilibrium restoration occurs due to energy exchange between the photons and leptons. While the initial and final equilibrium spectra are thermal, the transitional photon spectra are characterized by non-thermal features such as power-law tails, high energy bumps, and multiple components. Such non-thermal features are observed at infinity if the dissipation occurs at small to moderate optical depths, and the spectrum is released before thermalization is complete. We model the synthetic spectra with a Band function and show that the resulting spectral parameters are similar to observations for a frequency range of 2-3 orders of magnitude around the peak. In addition, our model predicts correlations between the low-frequency photon index and the peak frequency as well as between the low- and high-frequency indices. We explore baryon and pair dominated fireballs and reach the conclusion that baryonic fireballs are a better model for explaining the observed features of gamma-ray burst spectra.
NASA's Swift satellite has completed ten years of amazing discoveries in time domain astronomy. Its primary mission is to chase gamma-ray bursts (GRBs), but due to its scheduling flexibility it has subsequently become a prime discovery machine for new types of behavior. The list of major discoveries in GRBs and other transients includes the long-lived X-ray afterglows and flares from GRBs, the first accurate localization of short GRBs, the discovery of GRBs at high redshift (z>8), supernova shock break-out from SN Ib, a jetted tidal disruption event, an ultra-long class of GRBs, high energy emission from flare stars, novae and supernovae with unusual characteristics, magnetars with glitches in their spin periods, and a short GRB with evidence of an accompanying kilonova. Swift has developed a dynamic synergism with ground based observatories. In a few years gravitational wave observatories will come on-line and provide exciting new transient sources for Swift to study.
We report an investigation of the extremely metal-poor and C-rich planetary nebula (PN) K648 in the globular cluster M15 using the UV to far-IR data obtained using the Subaru, HST, FUSE, Spitzer, and Herschel. We determined the nebular abundances of ten elements. The enhancement of F ([F/H]=+0.96) is comparable to that of the halo PN BoBn1. The central stellar abundances of seven elements are determined. The stellar C/O ratio is similar to the nebular C/O ratios from recombination line and from collisionally excited line (CEL) within error, and the stellar Ne/O ratio is also close to the nebular CEL Ne/O ratio. We found evidence of carbonaceous dust grains and molecules including Class B 6-9 um and 11.3 um polycyclic aromatic hydrocarbons and the broad 11 um feature. The profiles of these bands are similar to those of the C-rich halo PNe H4-1 and BoBn1. Based on the theoretical model, we determined the physical conditions of the gas and dust and their masses, i.e., 0.048 Msun and 4.95x10^{-7} Msun, respectively. The observed chemical abundances and gas mass are in good agreement with an asymptotic giant branch nucleosynthesis model prediction for stars with an initial 1.25 Msun plus a 2.0x10^{-3} Msun partial mixing zone (PMZ) and stars with an initial mass of 1.5 Msun without a PMZ. The core-mass of the central star is approximately 0.61-0.63 Msun. K648 is therefore likely to have evolved from a progenitor that experienced coalescence or tidal disruption during the early stages of evolution, and became a ~1.25-1.5 Msun blue straggler.
Within the "nuclear medium cooling" scenario of neutron stars all reliably known temperature - age data, including those of the central compact objects in the supernova remnants of Cassiopeia A and XMMU-J1732, can be comfortably explained by a set of cooling curves obtained by variation of the star mass within the range of typical observed masses. The recent measurements of the masses of the pulsars PSR J1616-2230, PSR J0348-0432 and J00737-3039B and the companion of J1756-2251 provide independent proof for the existence of neutron stars with masses in a broad range from 1.2 to 2 $M_\odot$. The values $M>2 M_{\odot}$ call for sufficiently stiff equations of state for neutron star matter. We investigate the response of the set of neutron star cooling curves to a stiffening of the nuclear equation of state so that maximum masses of about $2.4~M_\odot$ would be accessible and to a deconfinement phase transition from such stiff nuclear matter in the outer core to color superconducting quark matter in the inner core. Without readjustment of cooling inputs the mass range required to cover all cooling data for the stiff DD2 equation of state should include masses of $2.426~M_\odot$ for describing the fast cooling of CasA while the existence of a quark matter core accelerates the cooling so that CasA cooling data are described with a hybrid star of mass $1.674~M_\odot$.
The joint likelihood is a simple extension of the standard likelihood formalism that enables the estimation of common parameters across disjoint datasets. Joining the likelihood, rather than the data itself, means nuisance parameters can be dealt with independently. Application of this technique, particularly to Fermi-LAT dwarf spheroidal analyses, has already been met with great success. We present a description of the method's general implementation along with a toy Monte-Carlo study of its properties and limitations.
The rapidly declining population of bright quasars at z~3 appears to make an increasingly small contribution to the ionising background at the HI Lyman limit. It is then generally though that massive stars in (pre-)galactic systems may provide the additional ionising flux needed to complete HI reionisation by z>6. A galaxy dominated background, however, may require that the escape fraction of Lyman continuum radiation from high redshift galaxies is as high as 10%, a value somewhat at odds with (admittedly scarce) observational constraints. High escape fractions from dwarf galaxies have been advocated, or, alternatively, a so-far undetected (or barely detected) population of unobscured, high-redshift faint AGNs. Here we question the latter hypothesis, and show that such sources, to be consistent with the measured level of the unresolved X-ray background at z=0, can provide a fraction of the HII filling factor not larger than 13% by z=6. The fraction rises to <27% in the somewhat extreme case of a constant comoving redshift evolution of the AGN emissivity. This still calls for a mean escape fraction of ionising photons from high-z galaxies >10%.
Kinetic mixing between the metric and scalar degrees of freedom is an essential ingredient in contemporary scalar-tensor theories. This often makes hard to understand their physical content, especially when derivative mixing is present, as it is the case for Horndeski action. In this work we develop a method that allows to write a Ricci curvature-free scalar field equation and discuss some of the advantages of such rephrasing in the study of stability issues in the presence of matter, the existence of an Einstein frame and the generalization of the disformal screening mechanism. For quartic Horndeski theories, such procedure leaves, in general, a residual coupling to curvature, given by the Weyl tensor. This gives rise to a binary classification of scalar-tensor theories into stirred theories, for which the curvature can be substituted for, and shaken theories for which a residual coupling to curvature remains. Quite remarkably, we have found that generalized DBI Galileons belong to the first class. Finally, we discuss kinetic mixing in quintic theories for which non-linear mixing terms appears and in the recently proposed theories beyond Horndeski which display a novel form of kinetic mixing, in which the field equation is sourced by derivatives of the energy-momentum tensor.
We present LUXCalc, a new utility for calculating likelihoods and deriving WIMP-nucleon coupling limits from the recent results of the LUX direct search dark matter experiment. After a brief review of WIMP-nucleon scattering, we derive, for the first time, LUX limits on the spin-dependent WIMP-nucleon couplings over a broad range of WIMP masses, under standard assumptions on the relevant astrophysical parameters. We find that, under these and other common assumptions, LUX excludes the entire spin-dependent parameter space consistent with a dark matter interpretation of DAMA's anomalous signal, the first time a single experiment has been able to do so. We also revisit the case of spin-independent couplings, and demonstrate good agreement between our results and the published LUX results. Finally, we derive constraints on the parameters of an effective dark matter theory in which a spin-1 mediator interacts with a fermionic WIMP and Standard Model fermions via axial-vector couplings. A detailed appendix describes the use of LUXCalc with standard codes to place constraints on generic dark matter theories.
If forthcoming measurements of cosmic photon polarization restrict the primordial tensor-to-scalar ratio to $r < 0.01$, small field inflation will be a principal candidate for the origin of the universe. Here we show that small multifield inflation, without the hybrid mechanism, typically results in large squeezed nongaussianity. Small multifield potentials contain multiple flat field directions, often identified with the gauge invariant field directions in supersymmetric potentials. We find that unless these field directions have equal slopes, large nongaussianity arises. After identifying relevant differences between large and small two-field potentials, we demonstrate that the latter naturally fulfill the Byrnes-Choi-Hall large nongaussianity conditions. Computations of the primordial power spectrum, spectral index, and squeezed bispectrum, reveal that small two-field models which otherwise match observed primordial perturbations, produce excludably large nongaussianity if the inflatons' field directions have unequal slopes.
In a universe where, according to the standard cosmological models, some 97% of the total mass-energy is still "missing in action" it behooves us to spend at least a little effort critically assessing and exploring radical alternatives. Among possible, (dare we say plausible), nonstandard but superficially viable models, those spacetimes conformal to the standard Friedmann-Lemaitre-Robertson-Walker class of cosmological models play a very special role --- these models have the unique and important property of permitting large non-perturbative geometric deviations from Friedmann-Lemaitre-Robertson-Walker cosmology without unacceptably distorting the cosmic microwave background. Performing a "cosmographic" analysis, (that is, temporarily setting aside the Einstein equations, since the question of whether or not the Einstein equations are valid on galactic and cosmological scales is essentially the same question as whether or not dark matter/dark energy actually exist), and using both supernova data and information about galactic structure, one can nevertheless place some quite significant observational constraints on any possible conformal mode --- however there is still an extremely rich range of phenomenological possibilities for both cosmologists and astrophysicists to explore.
We present a numerical code for multi-component simulation of the galactic evolution. Our code includes the following parts: $N$-body is used to evolve dark matter, stellar dynamics and dust grains, gas dynamics is based on TVD-MUSCL scheme with the extra modules for thermal processes, star formation, magnetic fields, chemical kinetics and multi-species advection. We describe our code in brief, but we give more details for the magneto-gas dynamics. We present several tests for our code and show that our code have passed the tests with a reasonable accuracy. Our code is parallelized using the MPI library. We apply our code to study the large scale dynamics of galactic discs.
We present a nonlinear post-Friedmann framework for structure formation,
generalizing to cosmology the weak-field (post-Minkowskian) approximation,
unifying the treatment of small and large scales. We consider a universe filled
with a pressureless fluid and a cosmological constant $\Lambda$, the theory of
gravity is Einstein's general relativity and the background is the standard
flat $\Lambda$CDM cosmological model.
We expand the metric and the energy-momentum tensor in powers of $1/c$,
keeping the matter density and peculiar velocity as exact fundamental
variables. We assume the Poisson gauge, including scalar and tensor modes up to
$1/c^4$ order and vector modes up to $1/c^5$ terms. Through a redefinition of
the scalar potentials as a resummation of the metric contributions at different
orders, we obtain a complete set of nonlinear equations, providing a unified
framework to study structure formation from small to superhorizon scales, from
the nonlinear Newtonian to the linear relativistic regime. We explicitly show
the validity of our scheme in the two limits: at leading order we recover the
fully nonlinear equations of Newtonian cosmology; when linearized, our
equations become those for scalar and vector modes of first-order relativistic
perturbation theory in the Poisson gauge. Tensor modes are non-dynamical at the
$1/c^4$ order we consider: they are purely nonlinear and describe a distortion
of the spatial slices determined at this order by a constraint, quadratic in
the scalar and vector variables. The main results of our analysis are: at
leading order a purely Newtonian nonlinear energy current sources a
frame-dragging gravitomagnetic vector potential; in the leading-order Newtonian
regime and in the linear relativistic regime the two scalar metric potentials
are the same, while the nonlinearity of general relativity makes them
different.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
Large solar proton events (SPEs) affect the solar-terrestrial space environment and become a very important aspect in space weather research. In this work, we statistically investigate 78 solar proton events of 1996-2011 and find that there exists a longitudinally asymmetric distribution of flare sources of the solar proton events observed near 1 AU, namely, with the same longitude separation between magnetic field line footpoint of observer and flare sources, the number of the solar proton events originating from sources located at eastern side of the nominal magnetic footpoint of observer is much larger than that of the solar proton events originating from sources located at western side. A complete model calculation of solar energetic particle (SEP) propagation in the three-dimensional Parker interplanetary magnetic field is presented to give a numerical explanation for this longitudinally asymmetric distribution phenomenon. We find that the longitudinally asymmetric distribution of solar proton events results from the east-west azimuthal asymmetry in the topology of the Parker interplanetary magnetic field as well as the effects of perpendicular diffusion on the transport of SEPs in the heliosphere. Our results would be valuable in understanding the solar-terrestrial relations and useful in space weather forecasting.
We present a new algorithm designed to improve the signal to noise ratio (SNR) of point and extended source detections in direct imaging data. The novel part of our method is that it finds the linear combination of the science images that best match counterpart images with signal removed from suspected source regions. The algorithm, based on the Locally Optimized Combination of Images (LOCI) method, is called Matched LOCI or MLOCI. We show using data obtained with the Gemini Planet Imager (GPI) and Near-Infrared Coronagraphic Imager (NICI) that the new algorithm can improve the SNR of point source detections by 30-400% over past methods. We also find no increase in false detections rates. No prior knowledge of candidate companion locations is required to use MLOCI. While non-blind applications may yield linear combinations of science images which seem to increase the SNR of true sources by a factor > 2, they can also yield false detections at high rates. This is a potential pitfall when trying to confirm marginal detections or to re-detect point sources found in previous epochs. Our findings are relevant to any method where the coefficients of the linear combination are considered tunable, e.g. LOCI and Principal Component Analysis (PCA). Thus we recommend that false detection rates be analyzed when using these techniques.
The apparent anisotropies of the galaxy clustering in observable redshift space provide a unique opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales via the Alcock--Paczynski effect and redshift-space distortions. While the improved theoretical models have been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear scales, the applicability of these models is still limited in the presence of the non--perturbative smearing effect caused by the randomness of the relative velocities. Although the cosmological constraint from the anisotropic clustering will be certainly improved with a more elaborate theoretical model, we here consider an alternative way by using the statistical power of both the power spectrum and bispectrum at large scales. Based on the Fisher matrix analysis, we estimate the benefit of combining the power spectra and bispectra, finding that the constraints on the cosmic expansion and growth of structure will be improved by a factor of two. This compensates for the loss of constraining power using the power spectrum alone due to the randomness of the relative velocities.
We present Orion A giant molecular cloud core catalogs, which are based on 1.1 mm map with an angular resolution of 36 arcsec (sim 0.07 pc) and C18O (1-0) data with an angular resolution of 26.4 arcsec (sim 0.05 pc). We have cataloged 619 dust cores in the 1.1 mm map using the Clumpfind method. The ranges of the radius, mass, and density of these cores are estimated to be 0.01 - 0.20 pc, 0.6 - 1.2 times 10^2 Msun, and 0.3 times 10^4 - 9.2 times 10^6 cm^{-3}, respectively. We have identified 235 cores from the C18O data. The ranges of the radius, velocity width, LTE mass, and density are 0.13 -- 0.34 pc, 0.31 - 1.31 km s^{-1}, 1.0 - 61.8 Msun, and (0.8 - 17.5) times 10^3 cm^{-3}, respectively. From the comparison of the spatial distributions between the dust and C18O cores, four types of spatial relations were revealed: (1) the peak positions of the dust and C18O cores agree with each other (32.4% of the C18O cores), (2) two or more C18O cores are distributed around the peak position of one dust core (10.8% of the C18O cores), (3) 56.8% of the C18O cores are not associated with any dust cores, and (4) 69.3% of the dust cores are not associated with any C18O cores. The data sets and analysis are public.
Massive black hole binaries are naturally predicted in the context of the hierarchical model of structure formation. The binaries that manage to lose most of their angular momentum can coalesce to form a single remnant. In the last stages of this process, the holes undergo an extremely loud phase of gravitational wave emission, possibly detectable by current and future probes. The theoretical effort towards obtaining a coherent physical picture of the binary path down to coalescence is still underway. In this paper, for the first time, we take advantage of observational studies of active galactic nuclei evolution to constrain the efficiency of gas-driven binary decay. Under conservative assumptions we find that gas accretion toward the nuclear black holes can efficiently lead binaries of any mass forming at high redshift (> 2) to coalescence within the current time. The observed "downsizing" trend of the accreting black hole luminosity function further implies that the gas inflow is sufficient to drive light black holes down to coalescence, even if they bind in binaries at lower redshifts, down to z~0.5 for binaries of ~10 million solar masses, and z~0.2 for binaries of ~1 million solar masses. This has strong implications for the detection rates of coalescing black hole binaries of future space-based gravitational wave experiments.
The recent discovery of a diffuse neutrino flux up to PeV energies raises the question of which populations of astrophysical sources contribute to this diffuse signal. One extragalactic candidate source population to produce high-energy neutrinos are Blazars. We present results from a likelihood analysis searching for cumulative neutrino emission from Blazar populations selected with the 2nd Fermi-LAT AGN catalog (2LAC) using an IceCube data set that has been optimized for the detection of individual sources. In contrast to previous searches with IceCube, the investigated populations contain up to hundreds of sources, the biggest one being the entire Blazar sample measured by the Fermi-LAT. No significant neutrino signal was found from any of these populations. Some implications of this non-observation for the origin of the observed PeV diffuse signal will be discussed.
In a radiatively heated and cooled medium, the thermal instability is a plausible mechanism for forming clouds, while the radiation force provides a natural acceleration, especially when ions recombine and opacity increases. Here we extend Field's theory to self-consistently account for a radiation force resulting from bound-free and bound-bound transitions in the optically thin limit. We present physical arguments for clouds to be significantly accelerated by a radiation force due to lines during a nonlinear phase of the instability. To qualitatively illustrate our main points, we perform both one and two-dimensional (1-D/2-D) hydrodynamical simulations that allow us to study the nonlinear outcome of the evolution of thermally unstable gas subjected to this radiation force. Our 1-D simulations demonstrate that the thermal instability can produce long-lived clouds that reach a thermal equilibrium between radiative processes and thermal conduction, while the radiation force can indeed accelerate the clouds to supersonic velocities. However, our 2-D simulations reveal that a single cloud with a simple morphology cannot be maintained due to destructive processes, triggered by the Rayleigh-Taylor instability and followed by the Kelvin-Helmholtz instability. Nevertheless, the resulting cold gas structures are still significantly accelerated before they are ultimately dispersed.
We investigate the possibility to detect the integrated Sachs-Wolfe (ISW) effect by cross-correlating 21-cm surveys at high redshifts with galaxies, in a way similar to the usual CMB-galaxy cross-correlation. The high-redshift 21-cm signal is dominated by CMB photons that travel freely without interacting with the intervening matter, and hence its late-time ISW signature should correlate extremely well with that of the CMB at its peak frequencies. Using the 21-cm temperature brightness instead of the CMB would thus be a further check of the detection of the ISW effect, measured by different instruments at different frequencies and suffering from different systematics. We also study the ISW effect on the photons that are scattered by HI clouds. We show that a detection of the unscattered photons is achievable with planned radio arrays, while one using scattered photons will require advanced radio interferometers, either an extended version of the planned Square Kilometre Array or futuristic experiments such as a lunar radio array.
We present radio images within 30$''$ of Sgr A* based on recent VLA observations at 34 GHz with 7.8 microJy sensitivity and resolution $\sim88\times46$ milliarcseconds (mas). We report 44 partially resolved compact sources clustered in two regions in the E arm of ionized gas that orbits Sgr A*. These sources have size scales ranging between ~50 and 200 mas (400 to 1600 AUs), and a bow-shock appearance facing the direction of Sgr A*. Unlike the bow-shock sources previously identified in the near-IR but associated with massive stars, these 34 GHz sources do not appear to have near-IR counterparts at 3.8 $\mu$m. We interpret these sources as a candidate population of photoevaporative protoplanetary disks (proplyds) that are associated with newly formed low mass stars with mass loss rates ~10^{-7} - 10^{-6} solar mass per year and are located at the edge of a molecular cloud outlined by ionized gas. The disks are externally illuminated by strong Lyman continuum radiation from the ~100 OB and WR massive stars distributed within 10'' of Sgr A*. The presence of proplyds implies current in-situ star formation activity near Sgr A* and opens a window for the first time to study low mass star, planetary and brown dwarf formations near a supermassive black hole.
It is widely thought that core-collapse supernovae (CCSNe), the explosions of massive stars following the collapse of the stars' iron cores, is obtained due to energy deposition by neutrinos. So far, this scenario was not demonstrated from first principles. Kushnir and Katz (2014) have recently shown, by using one-dimensional simulations, that if the neutrinos failed to explode the star, a thermonuclear explosion of the outer shells is possible for some (tuned) initial profiles. However, the energy released was small and negligible amounts of ejected $^{56}$Ni were obtained, implying that these one-dimensional collapse induced thermonuclear explosions (CITE) are unlikely to represent typical CCSNe. Here I provide evidence supporting a scenario in which the majority of CCSNe are the result of CITE. I use two-dimensional simulations to show that collapse of stars that include slowly (few percent of breakup) rotating $\sim0.1-10\,M_{\odot}$ shells of mixed helium-oxygen, leads to an ignition of a thermonuclear detonation wave that unbinds the stars' outer layers. Simulations of massive stars with different properties show that CITE is a robust process, and results in explosions with kinetic energies in the range of $10^{49}-10^{52}\,\textrm{erg}$, and $^{56}$Ni yields of up to $\sim\,M_{\odot}$, which are correlated, in agreement with observations for the majority of CCSNe. Stronger explosions are predicted from higher mass progenitors that leave more massive remnants, in contrast to the neutrino mechanism. Neutron stars are produced in weak ($\lt10^{51}\,\textrm{erg}$) explosions, while strong ($\gt10^{51}\,\textrm{erg}$) explosions leave black hole remnants.
Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$\alpha$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302.
Graham et al. (2015) discovered a supermassive black hole binary (SMBHB) candidate and identified the detected 5.2 yr period of the optical variability as the orbital period of the binary. Hydrodynamical simulations predict multiple periodic components for the variability of SMBHBs, thus raising the possibility that the true period of the binary is different from 5.2 yr. We analyse the periodogram of PG1302 and find no compelling evidence for additional peaks. We derive upper limits on any additional periodic modulations in the available data, by modeling the light-curve as the sum of red noise and the known 5.2 yr periodic component, and injecting additional sinusoidal signals. We find that, with the current data, we would be able to detect with high significance (false alarm probability <1%) secondary periodic terms with periods in the range predicted by the simulations, if the amplitude of the variability was at least ~0.07 mag (compared to 0.14 mag for the main peak). A three-year follow-up monitoring campaign with weekly observations can increase the sensitivity for detecting secondary peaks ~50%, and would allow a more robust test of predictions from hydrodynamical simulations.
We discuss a new scale-discretised directional wavelet transform to analyse spin signals defined on the sphere, in particular the polarisation of the cosmic microwave background (CMB).
Observations of Gamma-Ray Bursts with the Fermi Large Area Telescope have prompted theoretical advances and posed big challenges in the understanding of such extreme sources, despite the fact that GRB emission above 100 MeV is a fairly rare event. The first Fermi/LAT GRB catalog, published a year ago, presented 28 detections out of ~300 bursts detected by the Fermi Gamma-Ray Burst Monitor (GBM) within the LAT field of view. Building on the results from that work and on recent development in the understanding of the systematic errors on GBM localizations, we developed a new detection algorithm which increased the number of detections by 40 %. Even more recently the development of the new event analysis for the LAT ("Pass 8") has increased the number of detections within the first 3 years of the mission to 45, up 50 % with respect to the published catalog. The second LAT GRB catalog, in preparation, will cover more than 6 years of the mission and will break the barrier of 100 detected GRBs, a more than 20-fold improvement with respect to observations before the Fermi era in the same energy range. We will review the main features of the new algorithm, as well as preliminary results from this investigation.
We describe how to apply the transport method to compute inflationary observables in a broad range of multiple-field models. The method is efficient and encompasses scenarios with curved field-space metrics, violations of slow-roll conditions and turns of the trajectory in field space. It can be used for an arbitrary mass spectrum, including massive modes and models with quasi-single-field dynamics. In this note we focus on practical issues. It is accompanied by a Mathematica code which can be used to explore suitable models, or as a basis for further development.
The far-infrared (FIR) lines are key tracers of the physical conditions of the interstellar medium (ISM) and are becoming workhorse diagnostics for galaxies throughout the universe. Our goal is to explain the differences and trends observed in the FIR line emission of dwarf galaxies compared to more metal-rich galaxies. We present Herschel PACS spectroscopic observations of the CII157um, OI63 and 145um, OIII88um, NII122 and 205um, and NIII57um fine-structure cooling lines in a sample of 48 low-metallicity star-forming galaxies of the guaranteed time key program Dwarf Galaxy Survey. We correlate PACS line ratios and line-to-LTIR ratios with LTIR, LTIR/LB, metallicity, and FIR color, and interpret the observed trends in terms of ISM conditions and phase filling factors with Cloudy radiative transfer models. We find that the FIR lines together account for up to 3 percent of LTIR and that star-forming regions dominate the overall emission in dwarf galaxies. Compared to metal-rich galaxies, the ratios of OIII/NII122 and NIII/NII122 are high, indicative of hard radiation fields. In the photodissociation region (PDR), the CII/OI63 ratio is slightly higher than in metal-rich galaxies, with a small increase with metallicity, and the OI145/OI63 ratio is generally lower than 0.1, demonstrating that optical depth effects should be small on the scales probed. The OIII/OI63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and is found ~4 times higher in the dwarfs than in metal-rich galaxies. The high CII/LTIR, OI/LTIR, and OIII/LTIR ratios, which decrease with increasing LTIR and LTIR/LB, are interpreted as a combination of moderate FUV fields and low PDR covering factor. Harboring compact phases of low filling factor and a large volume filling factor of diffuse gas, the ISM of low-metallicity dwarf galaxies has a more porous structure than that in metal-rich galaxies.
We use three dimensional magnetohydrodynamic (MHD) simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter $Q$, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acceleration mechanisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter $Q$ maps with observations proves to be a valuable tool to discern unambiguously which mechanism is taking place in the remnant of SN 1006, giving strong support to the quasi-parallel model.
The Fermi-LAT has observed new gamma-ray flares from the blazar B0218+357 during July 2014. While no significant change in the gamma-ray spectrum has been previously observed through the flaring phase in late-2012, during this recent high activity the source displayed an exceptionally hard spectrum. The latter led to the detection of very high energy (VHE, E > 100 GeV) gamma-rays from B0218+357 by the MAGIC telescopes, establishing this source as the most distant TeV emitter known to date. In addition to the detection of VHE emission, this blazar is of particular interest since it is known to be a double-image gravitationally lensed system with a lens delay of 11.46 +/- 0.16 days measured in gamma-rays. We present the Fermi-LAT study of the July 2014 flares and discuss them in the context of previous measurements.
The period-luminosity sequences and the multiple periods of luminous red giant stars are examined using the OGLE III catalogue of long-period variables in the Large Magellanic Cloud. It is shown that the period ratios in individual multimode stars are systematically different from the ratios of the periods at a given luminosity of different period-luminosity sequences. This leads to the conclusion that the masses of stars at the same luminosity on the different period-luminosity sequences are different. An evolutionary scenario is used to show that the masses of stars on adjacent sequences differ by about 16-26% at a given luminosity, with the shorter period sequence being more massive. The mass is also shown to vary across each sequence by a similar percentage, with the mass increasing to shorter periods. On one sequence, sequence B, the mass distribution is shown to be bimodal. It is shown that the small amplitude variables on sequences A', A and B pulsate in radial and nonradial modes of angular degree l=0, 1 and 2, with the l=1 mode being the most common. The stars on sequences C' and C are predominantly radial pulsators (l=0). Matching period ratios to pulsation models shows that the radial pulsation modes associated with sequences A', A, B, C' and C are the 4th, 3rd, 2nd and 1st overtones and the fundamental mode, respectively.
The observed properties of high redshift galaxies depend on the underlying foreground distribution of large scale structure, which distorts their intrinsic properties via gravitational lensing. We focus on the regime where the dominant contribution originates from a single lens and examine the statistics of gravitational lensing by a population of virialized and non-virialized structures using sub-mm galaxies at z ~ 2.6 and Lyman-break galaxies at redshifts z ~ 6-15 as the background sources. We quantify the effect of lensing on the luminosity function of the high redshift sources, focusing on the intermediate and small magnifications (mu < 3) which affect the majority of the background galaxies. We show that depending on the intrinsic properties of the background galaxies, gravitational lensing can significantly affect the observed luminosity function even when no obvious strong lenses are present. Finally, we find that in the case of the Lyman-break galaxies it is important to account for the surface brightness profiles of both the foreground and the background galaxies when computing the lensing statistics, which introduces a selection criterion for the background galaxies that can actually be observed. Not taking this criterion into account leads to an overestimation of the number densities of very bright galaxies by nearly two orders of magnitude.
This review summarises the invited presentation I gave on the Milky Way disc. The idea underneath was to touch those topics that can be considered hot nowadays in the Galactic disk research: the reality of the thick disk, the spiral structure of the Milky Way, and the properties of the outer Galactic disk. A lot of work has been done in recent years on these topics, but a coherent and clear picture is still missing. Detailed studies with high quality spectroscopic data seem to support a dual Galactic disk, with a clear separation into a thin and a thick component. Much confusion and very discrepant ideas still exist concerning the spiral structure of the Milky Way. Our location in the disk makes it impossible to observe it, and we can only infer it. This process of inference is still far from being mature, and depends a lot on the selected tracers, the adopted models and their limitations, which in many cases are neither properly accounted for, nor pondered enough. Finally, there are very different opinions on the size (scale length, truncation radius) of the Galactic disk, and on the interpretation of the observed outer disk stellar populations in terms either of external entities (Monoceros, Triangulus-Andromeda, Canis Major), or as manifestations of genuine disk properties (e.g., warp and flare).
The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.
Among the more than 1000 gamma-ray bursts observed by the Fermi Gamma-ray Space Telescope, a large fraction show narrow and hard spectra inconsistent with non-thermal emission, signifying optically thick emission from the photosphere. However, only a few of these bursts have spectra consistent with a pure Planck function. We will discuss the observational features of photospheric emission in these GRBs as well as in the ones showing multi-component spectra. We interpret the observations in light of models of subphotospheric dissipation, geometrical broadening and multi-zone emission, and show what we can learn about the dissipation mechanism and properties of GRB jets.
For about three years, it was known that precision spectrophotometry with FORS2 suffered from systematic errors that made quantitative observations of planetary transits impossible. We identified the Longitudinal Atmospheric Dispersion Compensator (LADC) as the most likely culprit, and therefore engaged in a project to exchange the LADC prisms with the uncoated ones from FORS1. This led to a significant improvement in the depth of FORS2 zero points, a reduction in the systematic noise, and should make FORS2 again competitive for transmission spectroscopy of exoplanets.
We introduce a classification scheme of the post-merger dynamics and gravitational-wave emission in binary neutron star mergers, after identifying a new mechanism by which a secondary peak in the gravitational-wave spectrum is produced. It is caused by a spiral deformation, the pattern of which rotates slower with respect to the double-core structure in center of the remnant. This secondary peak is typically well separated in frequency from the secondary peak produced by a nonlinear interaction between a quadrupole and a quasi-radial oscillation. The new mechanism allows for an explanation of low-frequency modulations seen in a number of physical characteristics of the remnant, such as the central lapse function, the maximum density and the separation between the two cores. We find empirical relations for both types of secondary peaks between their gravitational-wave frequency and the compactness of nonrotating individual neutron stars, that exist for fixed total binary masses. Our classification scheme may form the basis for the construction of detailed gravitational-wave templates of the post-merger phase. We find that the quasi-radial oscillation frequency of the remnant decreases with the total binary mass. For a given merger event our classification scheme may allow to determine the proximity of the measured total binary mass to the threshold mass for prompt black hole formation, which can, in turn, yield an estimate of the maximum neutron-star mass.
The massive and luminous star-forming region W49A is a well known Galactic candidate to probe the physical conditions and chemistry similar to those expected in external starburst galaxies. We aim to probe the physical and chemical structure of W49A on a spatial scale of ~0.8 pc based on the JCMT Spectral Legacy Survey, which covers the frequency range between 330 and 373 GHz. The wide 2x2 arcminutes field and the high spectral resolution of the HARP instrument on JCMT provides information on the spatial structure and kinematics of the cloud. For species where multiple transitions are available, we estimate excitation temperatures and column densities. We detected 255 transitions corresponding to 60 species in the 330-373 GHz range at the center position of W49A. Excitation conditions can be probed for 16 molecules. The chemical composition suggests the importance of shock-, PDR-, and hot core chemistry. Many molecular lines show a significant spatial extent across the maps including high density tracers (e.g. HCN, HNC, CS, HCO+) and tracers of UV-irradiation (e.g. CN and C2H). Large variations are seen between the sub-regions with mostly blue-shifted emission toward the Eastern tail, mostly red-shifted emission toward the Northern clump, and emission peaking around the expected source velocity toward the South-west clump. A comparison of column density ratios of characteristic species observed toward W49A to Galactic PDRs suggests that while the chemistry toward the W49A center is driven by a combination of UV-irradiation and shocks, UV-irradiation dominates for the Northern Clump, Eastern tail, and South-west clump regions. A comparison to a starburst galaxy and an AGN suggests similar C2H, CN, and H2CO abundances (with respect to the dense gas tracer 34CS) between the ~0.8 pc scale probed for W49A and the >1 kpc regions in external galaxies with global star-formation.
The detection of olivine on Vesta is interesting because it may provide critical insights into planetary differentiation early in our Solar System's history. Ground-based and Hubble Space Telescope (HST) observations of asteroid (4) Vesta have suggested the presence of olivine on the surface. These observations were reinforced by the discovery of olivine-rich HED meteorites from Vesta in recent years. However, analysis of data from NASA's Dawn spacecraft has shown that this olivine-bearing unit is actually impact melt in the ejecta of Oppia crater. The lack of widespread mantle olivine, exposed during the formation of the 19 km deep Rheasilvia basin on Vesta's South Pole, further complicated this picture. Ammannito et al., (2013a) reported the discovery of local scale olivine-rich units in the form of excavated material from the mantle using the Visible and InfraRed spectrometer (VIR) on Dawn. Here we explore alternative sources for the olivine in the northern hemisphere of Vesta by reanalyzing the data from the VIR instrument using laboratory spectral measurements of meteorites. We suggest that these olivine exposures could be explained by the delivery of olivine-rich exogenic material. Based on our spectral band parameters analysis, the lack of correlation between the location of these olivine-rich terrains and possible mantle-excavating events, and supported by observations of HED meteorites, we propose that a probable source for olivine seen in the northern hemisphere are remnants of impactors made of olivine-rich meteorites. Best match suggests these units are HED material mixed with either ordinary chondrites, or with some olivine-dominated meteorites such as R-chondrites.
Stars close to the Eddington luminosity can have large low-density inflated envelopes. We show that the rise times of shock breakout signals from supernovae can be extended significantly if supernova progenitors have an inflated stellar envelope. If the shock breakout occurs in such inflated envelopes, the shock breakout signals diffuse in them, and their rise time can be significantly extended. Then, the rise times of the shock breakout signals are dominated by the diffusion time in the inflated envelope rather than the light-crossing time of the progenitors. We show that our inflated Wolf-Rayet star models whose radii are of the order of the solar radius can have shock breakout signals which are longer than ~100 sec. The existence of inflated envelopes in Wolf-Rayet supernova progenitors may be related to the mysterious long shock breakout signal observed in Type Ib SN 2008D. Extended shock breakout signals may provide evidence for the existence of inflated stellar envelopes and can be used to constrain the physical properties of these enigmatic structures.
We investigate the August 9, 2011 solar flare of X-ray class X6.9, the "hottest" flare from 2000 to 2012, with a peak plasma temperature according to GOES data of 32.5 MK. Our goal is to determine the cause of such an anomalously high plasma temperature and to investigate the energy balance in the flare region with allowance made for the presence of a super-hot plasma (>30 MK). We analyze the RHESSI, GOES, AIA/SDO, and EVE/SDO data and discuss the spatial structure of the flare region and the results of our spectral analysis of its X-ray emission. Our analysis of the RHESSI X-ray spectra is performed in the one-temperature and two-temperature approximations by taking into account the emission of hot (20 MK) and super-hot (45 MK) plasmas. The hard X-ray spectrum in both models is fitted by power laws. The observed peculiarities of the flare are shown to be better explained in terms of the two-temperature model, in which the super-hot plasma is located at the flare loop tops (or in the magnetic cusp region). The formation of the super-hot plasma can be associated with its heating through primary energy release and with the suppression of thermal conduction. The anomalously high temperature (32.5 MK according to GOES) is most likely to be an artefact of the method for calculating the temperature based on two-channel GOES measurements in the one-temperature approximation applied to the emission of a multi-temperature flare plasma with a minor contribution from the low-temperature part of the differential emission measure.
Large spectroscopic surveys have enabled in the recent years the computation of three-dimensional interstellar extinction maps thanks to accurate stellar atmospheric parameters and line-of-sight distances. Such maps are complementary to 3D maps extracted from photometry, allowing a more thorough study of the dust properties. Our goal is to use the high-resolution spectroscopic survey Gaia-ESO in order to obtain with a good distance resolution the interstellar extinction and its dependency as a function of the environment and the Galactocentric position. We use the stellar atmospheric parameters of more than 5000 stars, obtained from the Gaia-ESO survey second internal data release, and combine them with optical (SDSS) and near-infrared (VISTA) photometry as well as different sets of theoretical stellar isochrones, in order to calculate line-of-sight extinction and distances. The extinction coefficients are then compared with the literature to discuss their dependancy on the stellar parameters and position in the Galaxy. Within the errors of our method, our work does not show that there is any dependence of the interstellar extinction coefficient on the atmospheric parameters of the stars. We do not find any evidence of the variation of E(J-H)/E(J-K) with the angle from the Galactic centre nor with Galactocentric distance. This suggests that we are dealing with a uniform extinction law in the SDSS ugriz bands and the near-IR JHKs bands. Therefore, extinction maps using mean colour-excesses and assuming a constant extinction coefficient can be used without introducing any systematic errors.
The Local Group compact elliptical galaxy M32 hosts one of the nearest super-massive black holes (SMBHs), which has manifested itself only in X-rays to date. Based on sensitive observations taken with the {\it Karl G. Jansky} Very Large Array (VLA), we detect for the first time a compact radio source coincident with the nucleus of M32, which exhibits a flux density of $\sim$$47.3 \pm 5.9$ $\mu$Jy at 6.6 GHz. We discuss several possibilities for the nature of this source, favoring an origin of the long-sought radio emission from the central SMBH, for which we also revisit the X-ray properties based on recently acquired {\sl Chandra} and {\sl XMM-Newton} data. Our VLA observations also discover radio emission from three previously know optical planetary nebulae in the inner region of M32.
An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated ON/OFF observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of $\sim 9$ h of ON/OFF observations. Upper limits on the velocity averaged cross section, $<\sigma v >$, for the annihilation of dark matter particles with masses in the range of $\sim 300$ GeV to $\sim 10$ TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of $<\sigma v >$ that are larger than $3\cdot 10^{-24}\:\mathrm{cm^3/s}$ are excluded for dark matter particles with masses between $\sim 1$ and $\sim 4$ TeV at $95%$ CL if the radius of the central dark matter density core does not exceed $500$ pc. This is the strongest constraint that is derived on $<\sigma v>$ for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.
The influence of systematic errors on the calculation of the statistical significance of a $\gamma$-ray signal with the frequently invoked Li and Ma method is investigated. A simple criterion is derived to decide whether the Li and Ma method can be applied in the presence of systematic errors. An alternative method is discussed for cases where systematic errors are too large for the application of the original Li and Ma method. This alternative method reduces to the Li and Ma method when systematic errors are negligible. Finally, it is shown that the consideration of systematic errors will be important in many analyses of data from the planned Cherenkov Telescope Array.
The Fermi Large Area Telescope (LAT) is a powerful pulsar detector, as
demonstrated by the over one hundred objects in its second catalog of pulsars.
Pass 8 is a new reconstruction and event selection strategy developed by the
Fermi-LAT collaboration. Due to the increased acceptance at low energy, Pass 8
improves the pulsation detection sensitivity. Ten new pulsars rise above the 5
sigma threshold and are presented in this work, as well as one previously seen
with the former Pass 7 reconstruction.
More than 60$\%$ of the known pulsars with spin-down power ($\dot{E}$)
greater than $10^{36}$ erg/s show pulsations in gamma-rays, as seen with the
Fermi Large Area Telescope. Many non-detections of these energetic pulsars are
thought to be a consequence of a high background level, or a large distance
leading to a flux below the sensitivity limit of the instrument. The gamma-ray
beams of the others probably miss the Earth. The new Pass 8 data now allows the
detection of gamma ray pulsations from three of these high spin-down pulsars,
PSRs J1828$-$1101, J1831$-$0952 and J1837$-$0604, as well as three others with
$\dot{E}$ $\ge 10^{35}$ erg/s. We report on their properties and we discuss the
reasons for their detection with Pass 8.
We examine the applicability of the stochastic electron acceleration to two high synchrotron peaked blazars, Mrk 421 and Mrk 501, assuming synchrotron self-Compton emission of gamma-rays. Our model considers an emitting region moving at relativistic speed, where non-thermal electrons are accelerated and attain a steady-state energy spectrum together with the photons they emit. The kinetic equations of the electrons and photons are solved numerically, given a stationary wave number spectrum of the magnetohydrodynamic (MHD) disturbances, which are responsible for the electron acceleration and escape. Our simple formulation appears to reproduce the two well-sampled, long-term averaged photon spectra. In order to fit the model to the emission component from the radio to the X-ray bands, we need both a steeper wave spectral index than the Kolmogorov spectrum and efficient particle escape. Although the model provides a natural explanation for the high-energy cutoff of the electron energy distribution, the derived physical parameters raise a problem with an energy budget if the MHD waves with the Alfv{\'e}n velocity are assumed to be the acceleration agent.
We present the radial velocity confirmation of the extrasolar planet Kepler-447b, initially detected as a candidate by the Kepler mission. In this work, we analyze its transit signal and the radial velocity data obtained with the Calar Alto Fiber-fed Echelle spectrograph (CAFE). By simultaneously modeling both datasets, we obtain the orbital and physical properties of the system. According to our results, Kepler-447b is a Jupiter-mass planet ($M_p=1.37^{+0.48}_{-0.46} M_{\rm Jup}$), with an estimated radius of $R_p=1.65^{+0.59}_{-0.56} R_{\rm Jup}$ (uncertainties provided in this work are $3\sigma$ unless specified). This translates into a sub-Jupiter density. The planet revolves every $\sim7.8$ days around a G8V star with detected activity in the Kepler light curve. Kepler-447b transits its host with a large impact parameter ($b=1.076^{+0.112}_{-0.086}$), being one of the few planetary grazing transits confirmed so far and the first in the Kepler large crop of exoplanets. We estimate that only around 20% of the projected planet disk occults the stellar disk. The relatively large uncertainties in the planet radius are due to the large impact parameter and short duration of the transit. Within the transit time interval, we find the presence of large (somehow modulated) outliers during the transit. We propose and analyze different scenarios that could explain these brighter data points, including instrumental effects, additional perturbing bodies, stellar pulsations, rotation of a non-spherical planet, and spot-crossing events. However, short-cadence photometric data (at the 1 minute level) is still needed to unveil the nature of this observational effect.
Scaling relations in the LCDM Cosmology predict that for a given mass the
clusters formed at larger redshift are hotter, denser and therefore more
luminous in X-rays than their local z~0 counterparts. This effect overturns the
decrease in the observable X-ray flux so that it does not decrease at z > 1,
but instead raises slowly, similar to the SZ signal. Provided that scaling
relations remain valid at larger redshifts, X-ray surveys will not miss massive
clusters at any redshift, no matter how far they are.
At the same time, the difference in scaling with mass and distance of the
observable SZ and X-ray signals from galaxy clusters at redshifts z<2 offers a
possibility to crudely estimate the redshift and the mass of a cluster. This
might be especially useful for preselection of massive high-redshift clusters
and planning of optical follow-up for overlapping surveys in X-ray (e.g., by
SRG/eRosita) and SZ (e.g. Planck, SPT, ACT and CoRe+).
The primordial internal structures of gas giant planets are unknown. Often giant planets are modeled under the assumption that they are adiabatic, convective, and homogeneously mixed, but this is not necessarily correct. In this work, we present the first self-consistent calculation of convective transport of both heat and material as the planets evolve. We examine how planetary evolution depends on the initial composition and its distribution, whether the internal structure changes with time, and if so, how it affects the evolution. We consider various primordial distributions, different compositions, and different mixing efficiencies and follow the distribution of heavy elements in a Jupiter-mass planet as it evolves. We show that a heavy-element core cannot be eroded by convection if there is a sharp compositional change at the core-envelope boundary. If the heavy elements are initially distributed within the planet according to some compositional gradient, mixing occurs in the outer regions resulting in a compositionally homogeneous outer envelope. Mixing of heavy materials that are injected in a convective gaseous envelope are found to mix efficiently. Our work demonstrates that the primordial internal structure of a giant planet plays a substantial role in determining its long-term evolution and that giant planets can have non-adiabatic interiors. These results emphasize the importance of coupling formation, evolution, and internal structure models of giant planets self-consistently.
Context. It is generally agreed that hydrogenation reactions dominate chemistry on grain surfaces in cold, dense molecular cores, saturating the molecules present in ice mantles. Aims. We present a study of the low temperature reactivity of solid phase isocyanic acid (HNCO) with hydrogen atoms, with the aim of elucidating its reaction network. Methods. Fourier transform infrared spectroscopy and mass spectrometry were employed to follow the evolution of pure HNCO ice during bombardment with H atoms. Both multilayer and monolayer regimes were investigated. Results. The hydrogenation of HNCO does not produce detectable amounts of formamide (NH2CHO) as the major product. Experiments using deuterium reveal that deuteration of solid HNCO occurs rapidly, probably via cyclic reaction paths regenerating HNCO. Chemical desorption during these reaction cycles leads to loss of HNCO from the surface. Conclusions. It is unlikely that significant quantities of NH2CHO form from HNCO. In dense regions, however, deuteration of HNCO will occur. HNCO and DNCO will be introduced into the gas phase, even at low temperatures, as a result of chemical desorption.
Among the many intriguing aspects of optically discovered tidal disruption events is that their temperatures are lower than expected and that the temperature does not evolve as rapidly with decreasing fallback rate as would be expected in standard disk theory. We show that this can be explained qualitatively using an idea proposed by Laor & Davis in the context of normal active galactic nuclei: that larger accretion rates imply stronger winds and thus that the accretion rate through the inner disk only depends weakly on the inflow rate at the outer edge of the disk. We also show that reasonable quantitative agreement with data requires that, as has been suggested in recent papers, the circularization radius of the tidal stream is approximately equal to the semimajor axis of the most bound orbit of the debris rather than twice the pericenter distance as would be expected without rapid angular momentum redistribution. If this explanation is correct, it suggests that the evolution of tidal disruption events may test both non-standard disk theory and the details of the interactions of the tidal stream.
GAMMA-400 is a new space mission, designed as a dual experiment, capable to study both high energy gamma rays (from $\sim$100 MeV to few TeV) and cosmic rays (electrons up to 20 TeV and nuclei up to $\sim$10$^{15}$ eV). The full simulation framework of GAMMA-400 is based on the Geant4 toolkit. The details of the gamma-ray reconstruction pipeline in the three main instruments (Tracker, Imaging Calorimeter, Homogeneous Calorimeter) will be outlined. The performance of GAMMA-400 (PSF, effective area and sensitivity) have been obtained using this framework. The most updated results on them will be shown.
We use magnetohydrodynamical simulations of converging flows to investigate the process of molecular cloud formation and evolution out of the magnetised ISM. Here, we investigate whether the observed subcritical HI clouds can become supercritical and hence allow the formation of stars within them. To do so, we vary the turbulent Mach number of the flows, as well as the initial magnetic field strength. We show that dense cores are able to build up under all conditions, but that star formation in these cores is either heavily delayed or completely suppressed if the initial field strength is B>3 microGauss. To probe the effect of magnetic diffusion, we introduce a tilting angle between the flows and the uniform background magnetic field, which mimics non--ideal MHD effects. Even with highly inclined flows, the formed cores are devoid of star formation, because no magnetically supercritical regions are build up. Hence we conclude, that the problem of how supercritical cloud cores are generated still persists.
Small asteroids intersecting Earth's orbit can deliver extraterrestrial rocks to the Earth, called meteorites. This process is accompanied by a luminous phenomena in the atmosphere, called bolides or fireballs. Observations of bolides provide pre-atmospheric orbits of meteorites, physical and chemical properties of small asteroids, and the flux (i.e. frequency of impacts) of bodies at the Earth in the centimeter to decameter size range. In this chapter we explain the processes occurring during the penetration of cosmic bodies through the atmosphere and review the methods of bolide observations. We compile available data on the fireballs associated with 22 instrumentally observed meteorite falls. Among them are the heterogeneous falls Almahata Sitta (2008 TC$_3$) and Bene\v{s}ov, which revolutionized our view on the structure and composition of small asteroids, the P\v{r}\'{\i}bram-Neuschwanstein orbital pair, carbonaceous chondrite meteorites with orbits on the asteroid-comet boundary, and the Chelyabinsk fall, which produced a damaging blast wave. While most meteoroids disrupt into fragments during atmospheric flight, the Carancas meteoroid remained nearly intact and caused a crater-forming explosion on the ground.
Thin slumped glass foils are considered good candidates for the realization of future X-ray telescopes with large effective area and high spatial resolution. However, the hot slumping process affects the glass strength, and this can be an issue during the launch of the satellite because of the high kinematical and static loads occurring during that phase. In the present work we have investigated the possible use of Gorilla glass (produced by Corning), a chemical tempered glass that, thanks to its strength characteristics, would be ideal. The un-tempered glass foils were curved by means of an innovative hot slumping technique and subsequently chemically tempered. In this paper we show that the chemical tempering process applied to Gorilla glass foils does not affect the surface micro-roughness of the mirrors. On the other end, the stress introduced by the tempering process causes a reduction in the amplitude of the longitudinal profile errors with a lateral size close to the mirror length. The effect of the overall shape changes in the final resolution performance of the glass mirrors was studied by simulating the glass foils integration with our innovative approach based on glass reinforcing ribs. The preliminary tests performed so far suggest that this approach has the potential to be applied to the X-ray telescopes of the next generation.
We report a limit on the ultra-high-energy neutrino flux based on a non-detection of radio pulses from neutrino-initiated particle cascades in the Moon, in observations with the Parkes radio telescope undertaken as part of the LUNASKA project. Due to the improved sensitivity of these observations, which had an effective duration of 127 hours and a frequency range of 1.2-1.5 GHz, this limit extends to lower neutrino energies than those from previous lunar radio experiments, with a detection threshold below 10^20 eV. The calculation of our limit allows for the possibility of lunar-origin pulses being misidentified as local radio interference, and includes the effect of small-scale lunar surface roughness. The targeting strategy of the observations also allows us to place a directional limit on the neutrino flux from the nearby radio galaxy Centaurus A.
We describe the calculation of the stochastically heated dust emission using the 3D ray-tracing dust radiative transfer code DART-Ray, which is designed to solve the dust radiative transfer problem for galaxies with arbitrary geometries. In order to reduce the time required to derive the non-equilibrium dust emission spectra from each volume element within a model, we implemented an adaptive SED library approach, which we tested for the case of axisymmetric galaxy geometries. To show the capabilities of the code, we applied DART-Ray to a high-resolution N-body+SPH galaxy simulation to predict the appearance of the simulated galaxy at a set of wavelengths from the UV to the sub-mm. We analyse the results to determine the effect of dust on the observed radial and vertical profiles of the stellar emission as well as on the attenuation and scattering of light from the constituent stellar populations. We also quantify the proportion of dust re-radiated stellar light powered by young and old stellar populations, both bolometrically and as a function of infrared wavelength.
The efficiency of convection in stars affects many aspects of their evolution and remains one of the key-open questions in stellar modelling. In particular, the size of the mixed core in core-He-burning low-mass stars is still uncertain and impacts the lifetime of this evolutionary phase and, e.g., the C/O profile in white dwarfs. One of the known observables related to the Horizontal Branch (HB) and Asymptotic Giant Branch (AGB) evolution is the AGB bump. Its luminosity depends on the position in mass of the helium-burning shell at its first ignition, that is affected by the extension of the central mixed region. In this preliminary work we show how various assumptions on near-core mixing and on the thermal stratification in the overshooting region affect the luminosity of the AGB bump, as well as the period spacing of gravity modes in core-He-burning models.
We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch, assumed here to occur at a temperature $T_{\rm QCD} \sim 170$ MeV; this leaves an axion population of such low density that its late-time cosmological impact is negligible. For larger axion masses, $m_a \gtrsim 0.15$ eV, where axions remain in equilibrium until after the QCD phase transition, we find that a Euclid-like survey combined with Planck CMB data can detect $m_a$ at very high significance. Our conclusions are robust against assumptions about prior knowledge of the neutrino mass. Given that the proposed IAXO solar axion search is sensitive to $m_a\lesssim 0.2$ eV, the axion mass range probed by cosmology is nicely complementary.
This paper attempts to model possible Light Time Effect of TU UMa using a new code applied on formerly available and newly determined maxima timings in order to confirm binarity and refine parameters of the orbit of RRab component in binary system. The binary hypothesis is further tested also using radial velocity measurements. A new approach for determination of maxima timings based on template fitting which is also usable on sparse or scattered data is described. This approach was successfully applied on measurements from different sources. For determination of orbital parameters of a double star TU UMa we developed a new code for analysis of LiTE involving also secular variation in pulsation period. Its usability was successfully tested on CL Aur - an eclipsing binary with mass-transfer in a triple system showing similar changes in O-C diagram. Since orbital motion would cause systematic shifts in mean radial velocities (dominated by pulsations) we computed and compared our model with center-of-mass velocities. They were determined using high-quality templates of radial velocity curves of RRab stars. Maxima timings adopted from the GEOS database (225) together with those newly determined from sky surveys and new measurements (149) were used for construction of O-C diagram spanning more than five proposed orbital cycles. This data set is five times larger than data sets used by previous authors. Modelling of the O-C dependence resulted in 23.36-year orbital period which translates to the lowest mass of the second component of about 0.33 M_sun. Secular changes in pulsation period of TU UMa over the whole O-C diagram was satisfactorily approximated by parabolic trend with the rate of -2.33 ms/yr. For the confirmation of binarity we used radial velocity measurements from nine independent sources.
Context. In binary star systems, the winds from the two components impact
each other, leading to strong shocks and regions of enhanced density and
temperature. Potentially habitable circumbinary planets must continually be
exposed to these interactions regions.
Aims. We study, for the first time, the interactions between winds from
low-mass stars in a binary system, to show the wind conditions seen by
potentially habitable circumbinary planets.
Methods. We use the advanced 3D numerical hydrodynamic code Nurgush to model
the wind interactions of two identical winds from two solar mass stars with
circular orbits and a binary separation of 0.5 AU. As input into this model, we
use a 1D hydrodynamic simulation of the solar wind, run using the Versatile
Advection Code. We derive the locations of stable and habitable orbits in this
system to explore what wind conditions potentially habitable planets will be
exposed to during their orbits.
Results. Our wind interaction simulations result in the formation of two
strong shock waves separated by a region of enhanced density and temperature.
The wind-wind interaction region has a spiral shape due to Coriolis forces
generated by the orbital motions of the two stars. The stable and habitable
zone in this system extends from approximately 1.4 AU to 2.4 AU. (TRUNCATED)
This note is the first in a series of notes in which we aim to infer a model for the response of a Pierre Auger water Cherenkov detector to an ideal single muon. The main goal of the series is to provide analytical support for muon counting techniques. In this note we derive the probability distribution of the muon tracklength as a function of the zenith angle of the muon.
We draw attention to a problem with the alpha-Omega dynamo when it is applied to the origin of the galactic magnetic field under the assumption of perfect flux freezing. The standard theory involves the expulsion of undesirable flux and, because of flux freezing, the mass anchored on this flux also must be expelled. The strong galactic gravitational field makes this impossible on energetic grounds. It is shown that if only short pieces of the undesirable field lines are expelled, then mass can flow down along these field lines without requiring much energy. This expulsion of only short lines of force can be accomplished by a spike instability associated with gigantic astrophysical superbubbles. The physics of this instability is discussed and the results enable an estimate to be made of the number of spikes in the galaxy. It appears that there are probably enough spikes to cut all the undesirable lines into pieces as short as a couple of kiloparsecs during a dynamo time of a billion years. These cut pieces then may be randomly rotated in a dynamo time by alpha-Omega diffusion and there is enough rotation to get rid of the undesirable flux without expelling the fields themselves. The spike process seems strong enough to allows the alpha-Omega dynamo to create the galactic field without any trouble from the boundary condition problem.
Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.
The {\em Solar Orbiter} is the next solar physics mission of the European Space Agency, ESA, in collaboration with NASA, with a launch planned in 2018. The spacecraft is designed to approach the Sun to within 0.28\,AU at perihelion of a highly eccentric orbit. The proximity with the Sun will also allow its observation at uniformly high resolution at EUV and visible wavelengths. Such observations are central for learning more about the magnetic coupling of the solar atmosphere. At a later phase in the mission the spacecraft will leave the ecliptic and study the enigmatic poles of the Sun from a heliographic latitude of up to 33$^\circ$.
We report the discovery of 6.7 GHz methanol maser periodic flares in four massive star forming regions and the updated light curve for the known periodic source G22.357+0.066. The observations were carried out with the Torun 32 m radio telescope between June 2009 and April 2014. Flux density variations with period of 120 to 245 d were detected for some or all spectral features. A variability pattern with a fast rise and relatively slow fall on time-scale of 30-60 d dominated. A reverse pattern was observed for some features of G22.357+0.066, while sinusoidal-like variations were detected in G25.411+0.105. A weak burst lasting ~520 d with the velocity drift of 0.24 km/s/yr occurred in G22.357+0.066. For three sources for which high resolution maps are available, we found that the features with periodic behaviour are separated by more than 500 au from those without any periodicity. This suggests that the maser flares are not triggered by large-scale homogeneous variations in either the background seed photon flux or the luminosity of the exciting source and a mechanism which is able to produce local changes in the pumping conditions is required.
A diffuse flux of astrophysical neutrinos above $100\,\mathrm{TeV}$ has been observed at the IceCube Neutrino Observatory. Here we extend this analysis to probe the astrophysical flux down to $35\,\mathrm{TeV}$ and analyze its flavor composition by classifying events as showers or tracks. Taking advantage of lower atmospheric backgrounds for shower-like events, we obtain a shower-biased sample containing 129 showers and 8 tracks collected in three years from 2010 to 2013. We demonstrate consistency with the $(f_e:f_{\mu}:f_\tau)_\oplus\approx(1:1:1)_\oplus$ flavor ratio at Earth commonly expected from the averaged oscillations of neutrinos produced by pion decay in distant astrophysical sources. Limits are placed on non-standard flavor compositions that cannot be produced by averaged neutrino oscillations but could arise in exotic physics scenarios. A maximally track-like composition of $(0:1:0)_\oplus$ is excluded at $3.3\sigma$, and a purely shower-like composition of $(1:0:0)_\oplus$ is excluded at $2.3\sigma$.
Cosmic rays are an important source of heating in the interstellar medium, in particular in dense molecular cloud cores shielded from the external ultraviolet radiation field. The limits placed on the cosmic-ray ionization rate from measurements of the gas temperature in dense clouds are unaffected by the uncertainties associated to the traditional methods based on the analysis of molecular abundances. However, high-resolution data are required to determine with sufficient accuracy the spatial temperature distribution within prestellar cores. In this contribution we illustrate in detail the case of the well-studied prestellar core L1544, showing that both its thermal structure and chemical composition are consistent with a cosmic ray ionization rate of $\sim 10^{-17}$~s$^{-1}$, significantly smaller than the value measured in the diffuse interstellar medium. We also briefly discuss possible applications of this method to molecular clouds of other galaxies
Because of its nearness to Earth, the centre of the Milky Way is the only galaxy nucleus in which we can study the characteristics, distribution, kinematics, and dynamics of the stars on milli-parsec scales. We have accurate and precise measurements of the Galactic centre's central black hole, Sagittarius A*, and can study its interaction with the surrounding nuclear star cluster in detail. This contribution aims at providing a concise overview of our current knowledge about the Milky Way's central black hole and nuclear star cluster, at highlighting the observational challenges and limitations, and at discussing some of the current key areas of investigation.
For the first time, we use the Gamma-ray Burst Monitor (GBM) on-board the Fermi satellite to search for sterile neutrino decay lines in the energy range 10-25 keV corresponding to sterile neutrino mass range 20-50 keV. This energy range has been out of reach of traditional X-ray satellites such as Chandra, Suzaku, XMM-Newton, and gamma-ray satellites such as INTEGRAL. Furthermore, the extremely wide field of view of the GBM opens a large fraction of the Milky Way dark matter halo to be probed. We start with 1601 days worth of GBM data, implement stringent data cuts, and perform two simple line search analyses on the reduced data: in the first, the line flux is limited without background modeling, and in the second, the background is modeled as a power-law. We find no significant excess lines in both our searches. We set new limits on sterile neutrino mixing angles, improving on previous limits by approximately an order of magnitude. Better understanding of detector and astrophysical backgrounds, as well as detector response, can further improve the limit.
We measure the cross-power spectrum between galaxy density from Canada-France-Hawaii-Telescope Lensing Survey (CFHTLenS) catalogues and gravitational lensing convergence from Planck data release 1 (2013) and 2 (2015). We investigate three main galaxy samples: $18.0<i_{\rm AB}<22.0$, $18.0<i_{\rm AB}<23.0$, $18.0<i_{\rm AB}<24.0$ in the redshift range $0.2<z<1.3$ in each of the four CFHTLenS wide fields. By comparing the measured cross-spectrum with model predictions, linear galaxy-dark matter biases of $b=0.82^{+0.24}_{-0.23}, 0.83^{+0.19}_{-0.18}, 0.82^{+0.16}_{-0.14}$ are inferred at significances of $3.5, 4.5, 5.6\sigma$ using the Planck 2015 release. These measurements are marginally consistent with biases derived from galaxy-galaxy auto-correlations: $b=1.15^{+0.02}_{-0.01}, 1.08^{+0.01}_{-0.01}$ and $0.96^{+0.01}_{-0.01}$ respectively. Using the 2013 Planck release, we obtain biases of $b=1.33^{+0.29}_{-0.28}, 1.19^{+0.23}_{-0.23}, 1.16^{+0.19}_{-0.18}$, showing significant differences between the releases.
Radio and mm-wavelength observations of Sagittarius A* (Sgr A*), the radio
source associated with the supermassive black hole at the center of our Galaxy,
show that it behaves as a partially self-absorbed synchrotron-emitting source.
The measured size of Sgr A* shows that the mm-wavelength emission comes from a
small region and consists of the inner accretion flow and a possible collimated
outflow. Existing observations of Sgr A* have revealed a time lag between light
curves at 43 GHz and 22 GHz, which is consistent with a rapidly expanding
plasma flow and supports the presence of a collimated outflow from the
environment of an accreting black hole.
Here we wish to measure simultaneous frequency-dependent time lags in the
light curves of Sgr A* across a broad frequency range to constrain direction
and speed of the radio-emitting plasma in the vicinity of the black hole. Light
curves of Sgr A* were taken in May 2012 using ALMA at 100 GHz using the VLA at
48, 39, 37, 27, 25.5, and 19 GHz. As a result of elevation limits and the
longitude difference between the stations, the usable overlap in the light
curves is approximately four hours. Although Sgr A* was in a relatively quiet
phase, the high sensitivity of ALMA and the VLA allowed us to detect and fit
maxima of an observed minor flare where flux density varied by ~10%.
The fitted times of flux density maxima at frequencies from 100 GHz to 19
GHz, as well as a cross-correlation analysis, reveal a simple
frequency-dependent time lag relation where maxima at higher frequencies lead
those at lower frequencies. Taking the observed size-frequency relation of Sgr
A* into account, these time lags suggest a moderately relativistic (lower
estimates: 0.5c for two-sided, 0.77c for one-sided) collimated outflow.
Large surveys have shown that red galaxies are preferentially aligned with their halos while blue galaxies have a more isotropic distribution. Since halos generally align with their filaments this introduces a bias in the measurement of the cosmic shear from weak lensing. It is therefore vitally important to understand why this difference arises. We explore the stability of different disc orientations within triaxial halos. We show that, in the absence of gas, the disc orientation is most stable when its spin is along the minor axis of the halo. Instead when gas cools onto a disc it is able to form in almost arbitrary orientation, including off the main planes of the halo (but avoiding an orientation perpendicular to the halo's intermediate axis). Substructure helps gasless galaxies reach alignment with the halo faster, but have less effect on galaxies when gas is cooling onto the disc. Our results provide a novel and natural interpretation for why red, gas poor galaxies are preferentially aligned with their halo, while blue, star-forming, galaxies have nearly random orientations, without requiring a connection between galaxies' current star formation rate and their merger history.
This is an overview of pulsar accretion modeling. The physics of pulsar accretion, i.e., the process of plasma flow onto the neutron star surface, can be constrained from the spectral properties of the X-ray source. We discuss a new implementation of the physical continuum model developed by Becker and Wolff (2007, ApJ 654, 435). The model incorporates Comptonized blackbody, bremsstrahlung, and cyclotron emission. We discuss preliminary results of applying the new tool to the test cases of Suzaku data of Cen X-3 and XTE J1946+274. Cen X-3 is a persistent accreting pulsar with an O-star companion observed during a bright period. XTE J1946+274 is a transient accreting pulsar with a Be companion observed during a dim period. Both sources show spectra that are well described with an empirical Fermi Dirac cutoff power law model. We extend the spectral analysis by making the first steps towards a physical description of Cen X-3 and XTE J1946+274.
The IceCube neutrino observatory is composed of more than five thousand light sensors, Digital Optical Modules (DOMs), installed on the surface and at depths between 1450 and 2450 m in clear ice at the South Pole. Each DOM incorporates a 10-inch diameter photomultiplier tube (PMT) intended to detect light emitted when high energy neutrinos interact with atoms in the ice. Depending on the energy of the neutrino and the distance from secondary particle tracks, PMTs can be hit by up to several thousand photons within a few hundred nanoseconds. The number of photons per PMT and their time distribution is used to reject background events and to determine the energy and direction of each neutrino. The detector energy scale was established from previous lab measurements of DOM optical sensitivity, then refined based on observed light yield from stopping muons and calibration of ice properties. A laboratory setup has now been developed to more precisely measure the DOM optical sensitivity as a function of angle and wavelength. DOMs are calibrated in water using a broad beam of light whose intensity is measured with a NIST calibrated photodiode. This study will refine the current knowledge of the IceCube response and lay a foundation for future precision upgrades to the detector.
We consider Horndeski cosmological models able to screen the vacuum energy coming from any field theory assuming that after this screening the space should be in a de Sitter vacuum with a particular value of the cosmological constant specified by the theory of gravity itself. The most general scalar-tensor cosmological models without higher than second order derivatives in the field equations that have a spatially flat de Sitter critical point for any kind of material content or vacuum energy are, therefore, presented. These models could allow us to understand the current accelerated expansion of the universe as the result of a dynamical evolution towards a de Sitter attractor.
A hypothesis about the existence of semiweak interaction of electronic neutrinos with nucleons mediated by exchange of massless pseudoscalar bosons is stated. Owing to approximately 10 collisions of a solar neutrino with nucleons of the Sun, the fluxes of left- and right-handed solar neutrinos at the Earth surface are approximately equal, and their spectrum is changed in comparison with the one at the production moment. Good agreement is demonstrating between the calculated and experimental characteristics of the processes with solar neutrinos: ${}^{37}{\rm Cl} \rightarrow {}^{37}{\rm Ar}$, ${}^{71}{\rm Ga} \rightarrow {}^{71}{\rm Ge}$, $\nu_{e} e^{-}\rightarrow \nu_{e} e^{-}$, and $\nu_{e}D \rightarrow e^{-}pp$.
We present a predictive model of the nonlinear phase of the Weibel instability induced by two symmetric, counter-streaming ion beams in the non-relativistic regime. This self-consistent model combines the quasilinear kinetic theory of Davidson et al. [Phys. Fluids 15, 317 (1972)] with a simple description of current filament coalescence. It allows us to follow the evolution of the ion parameters up to a stage close to complete isotropization, and is thus of prime interest to understand the dynamics of collisionless shock formation. Its predictions are supported by 2-D and 3-D particle-in-cell simulations of the ion Weibel instability. The derived approximate analytical solutions reveal the various dependencies of the ion relaxation to isotropy. In particular, it is found that the influence of the electron screening can affect the results of simulations using an unphysical electron mass.
Assuming a double-well bare potential for a self-interacting scalar field, with the Higgs vacuum expectation value, it is shown that non-perturbative quantum corrections naturally lead to non-interacting ultra-light particles of mass $\simeq10^{-23}$eV, if these non-perturbative effects occur at a time consistent with the Electroweak phase transition. This mechanism could be relevant in the context of Bose Einstein Condensate studies for the description of cold Dark Matter.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
A new quantum effect connected with the late time behavior of decaying states is described and its possible observational consequences are analyzed: It is shown that charged unstable particles as well as neutral unstable particles with non--zero magnetic moment which live sufficiently long may emit electromagnetic radiation. This mechanism is due to the nonclassical behavior of unstable particles at late times (at the post exponential time region). Analyzing the transition times region between exponential and non-exponential form of the survival amplitude it is found that the instantaneous energy of the unstable particle can take very large values, much larger than the energy of this state at times from the exponential time region. Based on the results obtained for the model considered, it is shown that this new purely quantum mechanical effect may be responsible for causing unstable particles produced by astrophysical sources and moving with relativistic velocities to emit electromagnetic--, $X$-- or $\gamma$--rays at some time intervals from the transition time regions.
Exoplanetary systems closest to the Sun, with the brightest host stars, provide the most favorable opportunities for characterization studies of the host star and their planet(s). The Transit Ephemeris Refinement and Monitoring Survey uses both new radial velocity measurements and photometry in order to greatly improve planetary orbit uncertainties and the fundamental properties of the star, in this case HD 130322. The only companion, HD 130322b, orbits in a relatively circular orbit, e = 0.029 every ~10.7 days. Radial velocity measurements from multiple sources, including 12 unpublished from the Keck I telescope, over the course of ~14 years have reduced the uncertainty in the transit midpoint to ~2 hours. The transit probability for the b-companion is 4.7%, where M_p sin i = 1.15 M_J and a = 0.0925 AU. In this paper, we compile photometric data from the T11 0.8m Automated Photoelectric Telescope at Fairborn Observatory taken over ~14 years, including the constrained transit window, which results in a dispositive null result for both full transit exclusion of HD 130322b to a depth of 0.017 mag and grazing transit exclusion to a depth of ~0.001 mag. Our analysis of the starspot activity via the photometric data reveals a highly accurate stellar rotation period: 26.53 +/-0.70 days. In addition, the brightness of the host with respect to the comparison stars is anti-correlated with the Ca II H and K indices, typical for a young solar-type star.
We perform multi-plane ray-tracing using the GLAMER gravitational lensing code within high-resolution light-cones extracted from the CoDECS simulations: a suite of cosmological runs featuring a coupling between Dark Energy and Cold Dark Matter. We show that the presence of the coupling is evident not only in the redshift evolution of the normalisation of the convergence power spectrum, but also in differences in non-linear structure formation with respect to {\Lambda}CDM. Using a tomographic approach under the assumption of a {\Lambda}CDM cosmology, we demonstrate that weak lensing measurements would result in a {\sigma}8 value that changes with the source redshift if the true underlying cosmology is a coupled Dark Energy one. This provides a generic null test for these types of models. We also find that different models of coupled Dark Energy can show either an enhanced or a suppressed correlation between convergence maps with differing source redshifts as compared to {\Lambda}CDM. This would provide a direct way to discriminate between different possible realisations of the coupled Dark Energy scenario. Finally, we discuss the impact of the coupling on several lensing observables for different source redshifts and angular scales with realistic source redshift distributions for current ground-based and future space-based lensing surveys.
We examine the X-ray spectra of VII Zw 403, a nearby low-metallicity blue compact dwarf (BCD) galaxy. The galaxy has been observed to contain an X-ray source, likely a high mass X-ray binary (HMXB), with a luminosity of 1.3-23x10^38 erg s^-1 in the 0.3-8 keV energy range. A new Suzaku observation shows a transition to a luminosity of 1.7x10^40 erg s^-1 [0.3-8 keV], higher by a factor of 7-130. The spectra from the high flux state are hard, best described by a disk plus Comptonization model, and exhibit curvature at energies above 5 keV. This is consistent with many high-quality ultraluminous X-ray source spectra which have been interpreted as stellar mass black holes (StMBH) accreting at super-Eddington rates. However, this lies in contrast to another HMXB in a low-metallicity BCD, I Zw 18, that exhibits a soft spectrum at high flux, similar to Galactic black hole binaries and has been interpreted as a possible intermediate mass black hole. Determining the spectral properties of HMXBs in BCDs has important implications for models of the Epoch of Reionization. It is thought that the main component of X-ray heating in the early universe was dominated by HMXBs within the first galaxies. Early galaxies were small, metal-deficient, star forming galaxies with large H I mass fractions --- properties shared by local BCDs we see today. Understanding the spectral evolution of HMXBs in early universe analogue galaxies, such as BCDs, is an important step in estimating their contribution to the heating of the intergalactic medium during the Epoch of Reionization. The strong contrast between the properties of the only two spectroscopically studied HMXBs within BCDs motivates further study on larger samples of HMXBs in low metallicity environments in order to properly estimate the X-ray heating in the early universe.
We measure the stellar mass-star formation rate relation in star-forming disk galaxies at z<0.085, using Galaxy~Zoo morphologies to examine different populations of spirals as classified by their kiloparsec-scale structure. We examine the number of spiral arms, their relative pitch angle, and the presence of a galactic bar in the disk, and show that both the slope and dispersion of the M-SFR relation is constant when varying all the above parameters. We also show that mergers (both major and minor), which represent the strongest conditions for increases in star formation at a constant mass, only boost the SFR above the main relation by ~0.3 dex; this is significantly smaller than the increase seen in merging systems at z>1. Of the galaxies lying significantly above the M-SFR relation in the local Universe, more than 50% are mergers. We interpret this as evidence that the spiral arms, which are imperfect reflections of the galaxy's current gravitational potential, are either fully independent of the various quenching mechanisms or are completely overwhelmed by the combination of outflows and feedback. The arrangement of the star formation can be changed, but the system as a whole regulates itself even in the presence of strong dynamical forcing.
Black holes are believed to be one of the key ingredients of galaxy formation models, but it has been notoriously challenging to simulate them due to the very complex physics and large dynamical range of spatial scales involved. Here we address significant shortcomings of a Bondi-Hoyle-like prescription commonly invoked to estimate black hole accretion in cosmological hydrodynamic simulations of galaxy formation. We describe and implement a novel super-Lagrangian refinement scheme to increase, adaptively and 'on the fly', the mass and spatial resolution in targeted regions around the accreting black holes at limited computational cost. While our refinement scheme is generically applicable and flexible, for the purpose of this paper we select the smallest resolvable scales to match black holes' instantaneous Bondi radii, thus effectively resolving Bondi-Hoyle-like accretion in full galaxy formation simulations. This permits us to not only estimate gas properties close to the Bondi radius much more accurately, but also allows us to improve black hole accretion and feedback implementations. We thus devise a more generic feedback model where accretion and feedback depend on the geometry of the local gas distribution and where mass, energy and momentum loading are followed simultaneously. We present a series of tests of our refinement and feedback methods and apply them to models of isolated disc galaxies. Our simulations demonstrate that resolving gas properties in the vicinity of black holes is necessary to follow black hole accretion and feedback with a higher level of realism and that doing so allows us to incorporate important physical processes so far neglected in cosmological simulations.
We study how runaway stellar collisions in high redshift, metal poor star clusters form very massive stars (VMSs) that can directly collapse to intermediate mass black holes (IMBHs). We follow the evolution of a pair of two neighboring high-redshift mini-halos which are expected to host central nuclear star clusters (NSCs) with very high resolution, cosmological hydrodynamical zoom-in simulations with the adaptive mesh refinement code RAMSES. One of the two mini-halos enriches the central NSC of the other mini-halo to a critical metallicity, sufficient for Pop. II star formation at redshift z~27. We then use the spatial configuration of the flattened, asymmetrical gas cloud forming in the RAMSES simulations at the core of the metal enriched halo to set the initial conditions for simulations of an initially non-spherical star cluster with the direct summation code NBODY6 which are compared to about 2000 NBODY6 simulations of spherical star clusters for a wide range of star cluster parameters. In this way we establish that the results of our modeling are robust to changes in the assumptions for stellar initial mass function, binary fraction, and degree of initial mass mass segregation. Most of our simulations result in the formation of VMSs. The final mass of the VMS that forms depends thereby strongly on the initial mass of the NSC as well as the initial central density. For the initial central densities suggested by our RAMSES simulations, VMSs with M_VMS>400 M_sun can form in clusters with stellar masses of ~10^4 M_sun and this can increase to well over 1000 M_sun for more massive and denser clusters. The high probability we find for forming a VMS in these halos at such an early cosmic time makes collisional runaway of Pop. II star clusters a promising channel for producing large numbers of high-redshift IMBHs that may act as the seeds of super-massive black holes.
The Initial Mass Function (IMF) of early-type galaxies (ETGs) has been found to feature systematic variations by both dynamical and spectroscopic studies. In particular, spectral line strengths, based on gravity-sensitive features, suggest an excess of low-mass stars in massive ETGs, i.e. a bottom-heavy IMF. The physical drivers of IMF variations are currently unknown. The abundance ratio of alpha elements, such as [Mg/Fe], has been suggested as a possible driver of the IMF changes, although dynamical constraints do not support this claim. In this letter, we take advantage of the large SDSS database. Our sample comprises 24,781 high-quality spectra, covering a large range in velocity dispersion (100<sigma0<320 km/s) and abundance ratio (-0.1<[Mg/Fe]<+0.4). The large volume of data allows us to stack the spectra at fixed values of sigma0 and [Mg/Fe]. Our analysis -- based on gravity-sensitive line strengths -- gives a strong correlation with central velocity dispersion and a negligible variation with [Mg/Fe] at fixed sigma0. This result is robust against individual elemental abundance variations, and seems not to raise any apparent inconsistency with the alternative method based on galaxy dynamics.
Models for supernovae (SNe) arising from thermonuclear explosions of white dwarfs (WDs) have been extensively studied over the last few decades, mostly focusing on the single degenerate (accretion of material of a WD) and double degenerate (WD-WD merger) scenarios. In recent years it was suggested that WD-WD direct collisions provide an additional channel for such explosions. Here we extend the studies of such explosions, and explore the role of Helium-shells in affecting the thermonuclear explosions. We study both the impact of low-mass helium ($\sim0.01$ M$_{\odot})$ shells, as well as high mass shells ($\ge0.1$ M$_{\odot}$). We find that detonation of the massive helium layers precede the detonation of the WD Carbon-Oxygen (CO) bulk during the collision and can change the explosive evolution and outcomes for the cases of high mass He-shells. In particular, the He-shell detonation propagates on the WD surface and inefficiently burns material prior to the CO detonation that later follows in the central parts of the WD. Such evolution leads to larger production of intermediate elements, producing large yields of $^{44}{\rm Ti}$ and $^{48}{\rm Cr}$ relative to the pure CO-CO WD collisions. Collisions of WDs with a low-mass He-shell do not give rise to helium detonation, but helium burning does precede the CO bulk detonation. Such collisions produce a high velocity, low-mass of ejected burned material enriched with intermediate elements, with smaller changes to the overall explosion outcomes. The various effects arising from the contribution of low/high mass He layers change the kinematics and the morphological structure of collision-induced SNe and may thereby provide unique observational signatures for such SNe, and play a role in the chemical enrichment of galaxies and the production of intermediate elements and positrons from their longer-term decay.
Self-gravitating bosonic fields can support stable and localised field configurations. For real fields, these solutions oscillate in time and are known as oscillatons. The density profile is static, and is soliton. Such solitons should be ubiquitous in models of axion dark matter, with the soliton characteristic mass and size depending on some inverse power of the axion mass. Stable configurations of non-relativistic axions are studied numerically using the Schr\"{o}dinger-Poisson system. This method, and the resulting soliton density profiles, are reviewed. Using a scaling symmetry and the uncertainty principle, the core size of the soliton can be related to the central density and axion mass, $m_a$, in a universal way. Solitons have a constant central density due to pressure-support, unlike the cuspy profile of cold dark matter (CDM). One consequence of this fact is that solitons composed of ultra-light axions (ULAs) may resolve the `cusp-core' problem of CDM. In DM halos, thermodynamics will lead to a CDM-like Navarro-Frenk-White profile at large radii, with a central soliton core at small radii. Using Monte-Carlo techniques to explore the possible density profiles of this form, a fit to stellar-kinematical data of dwarf spheroidal galaxies is performed. In order for ULAs to resolve the cusp-core problem (without recourse to baryon feedback or other astrophysical effects) the axion mass must satisfy $m_a<1.1\times 10^{-22}\text{ eV}$ at 95\% C.L. On the other hand, ULAs with $m_a\lesssim 1\times 10^{-22}\text{ eV}$ are in some tension with cosmological structure formation. An axion solution to the cusp-core problem thus makes novel predictions for future measurements of the epoch of reionisation. On the other hand, this can be seen as evidence that structure formation could soon impose a \emph{Catch 22} on axion/scalar field DM, similar to the case of warm DM.
Non-attractor models of inflation are characterized by the super-horizon evolution of curvature perturbations, introducing a violation of the non-Gaussian consistency relation between the bispectrum's squeezed limit and the power spectrum's spectral index. In this work we show that the bispectrum's squeezed limit of non-attractor models continues to respect a relation dictated by the evolution of the background. We show how to derive this relation using only symmetry arguments, without ever needing to solve the equations of motion for the perturbations.
We report on recent optical observations of the stellar and the nebular remnants of 22 southern post-novae. In this study, for each of our targets, we obtained and analysed long-slit spectra in the spectral range 3500-6600 A and in H$\alpha$+NII narrow-band images. The changes in the emission lines' equivalent widths with the time since the outburst agree with earlier published results of other authors. We estimated an average value $\alpha$=2.37 for the exponent of the power law fitted to the post-novae continua. Our observations clearly show the two-component structure of the V842 Cen expanding nebulae, owing to the different velocities of the ejected matter. We discovered an expanding shell around V382 Vel with an outer diameter of about 12 arcsec.
The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $\Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is degenerate with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40% underestimate to a 5% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass inside $0.2R_{200}$.
The Interface Region Imaging Spectrograph (IRIS) reveals small-scale rapid brightenings in the form of bright grains all over coronal holes and the quiet sun. These bright grains are seen with the IRIS 1330 \AA, 1400 \AA\ and 2796 \AA\ slit-jaw filters. We combine coordinated observations with IRIS and from the ground with the Swedish 1-m Solar Telescope (SST) which allows us to have chromospheric (Ca II 8542 \AA, Ca II H 3968 \AA, H\alpha, and Mg II k 2796 \AA), and transition region (C II 1334 \AA, Si IV 1402) spectral imaging, and single-wavelength Stokes maps in Fe I 6302 \AA at high spatial (0.33"), temporal and spectral resolution. We conclude that the IRIS slit-jaw grains are the counterpart of so-called acoustic grains, i.e., resulting from chromospheric acoustic waves in a non-magnetic environment. We compare slit-jaw images with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 \AA\ slit-jaw filter comes from both the Mg II k core and wings. The signal in the C II and Si IV lines is too weak to explain the presence of grains in the 1300 and 1400 \AA\ slit-jaw images and we conclude that the grain signal in these passbands comes mostly from the continuum. Even though weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C II spectra. For some grains, spectral signature can be found in IRIS Si IV. This suggests that upward propagating acoustic waves sometimes reach all the way up to the transition region.
In the unification scheme of active galactic nuclei (AGN), Seyfert 1s and Seyfert 2s are intrinsically same, but they are viewed at different angles. However, the Fe K\alpha emission line luminosity of Seyfert 1s was found in average to be about twice of that of Seyfert 2s at given X-ray continuum luminosity in the previous work (Ricci et al. 2014). We construct an accretion disc-corona model, in which a fraction of energy dissipated in the disc is extracted to heat the corona above the disc. The radiation transfer equation containing Compton scattering processes is an integro-differential equation, which is solved numerically for the corona with a parallel plane geometry. We find that the specific intensity of X-ray radiation from the corona changes little with the viewing angle \theta when \theta is small (nearly face-on), and it is sensitive to \theta if the viewing angle is large (\theta> 40 degrees). The radiation from the cold disc, mostly in infrared/optical/UV bands, is almost proportional to cos\theta when \theta< 40 degrees, while it decreases more rapidly than cos\theta when \theta> 40 degrees because of strong absorption in the corona in this case. For seyfert galaxies, the Fe K\alpha line may probably be emitted from the disc irradiated by the X-ray continuum emission. The observed equivalent width (EW) difference between Seyfert 1s and Seyfert 2s can be reproduced by our model calculations, provided Seyfert 1s are observed in nearly face-on direction and the average inclination angle of Seyfert 2s ~65 dgrees.
The morphological types of 5840 galaxies were classified by a visual inspection of color images using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) to produce a morphology catalog of a representative sample of local galaxies with $z<0.01$. The sample galaxies are almost complete for galaxies brighter than $r_{pet}=17.77$. Our classification system is basically the same as that of the Third Reference Catalog of Bright Galaxies with some simplification for giant galaxies. On the other hand, we distinguish fine features of dwarf elliptical-like galaxies to classify 5 subtypes: dwarf ellipticals (dE), blue-cored dwarf ellipticals (dE$_{bc}$), dwarf spheroidals (dSph), blue dwarf ellipticals (dE$_{blue}$), and dwarf lenticulars (dS0). In addition, we denote the presence of nucleation in dE, dSph, and dS0. Elliptical galaxies and lenticular galaxies contribute only $\sim1.5\%$ and $\sim4.9\%$ of the local galaxies, respectively, whereas spirals and irregulars contribute $\sim32.1\%$ and $\sim42.8\%$, respectively. The dE$_{blue}$ galaxies, which are recently found populations of galaxies, contribute a significant fraction to the dwarf galaxies. There seems to be structural difference between dSph and dE galaxies. The dSph galaxies are fainter and bluer with shallower surface brightness gradient than dE galaxies. They also have lower fraction of galaxies with small axis ratios ($b/a \lesssim 0.4$) than dE galaxies. The mean projected distance to the nearest neighbor galaxy is $\sim260$kpc. About $1\%$ of local galaxies have no neighbors with comparable luminosity within a projected distance of 2Mpc.
Study on the small-scale structures and material flows of solar quiescent filaments is very important for understanding the formation and equilibrium of solar filaments. Using the high resolution H{\alpha} data observed by the New Vacuum Solar Telescope (NVST), we present the structures of the barbs and the material flows along the threads across the spine in two quiescent filaments on 2013 September 29 and on 2012 November 2, respectively. During the evolution of the filament barb, several parallel tube-shaped structures formed and the width of the structures ranges from about 2.3 Mm to 3.3 Mm. The parallel tube-shaped structures merged together accompanied with the material flows from the spine to the barb. Moreover, the boundary between the barb and surrounding atmosphere is very neat. The counter-streaming flows were not found to appear alternately in the adjacent threads of the filament. However, the large-scale patchy counter-streaming flows are detected in the filament. The flows in one patch of the filament have the same direction and the flows in the adjacent patch have opposite directions. The patches of two opposite flows with a size of about ten arcseconds exhibited alternately along the spine of the filament. The velocity of these material flows ranges from 5.6 km/s to 15.0 km/s. The material flows along the threads of the filament did not change their direction for about two hours and fourteen minutes during the evolution of the filament. Our results confirm that the large-scale counter-streaming flows with the certain width along the threads of solar filaments exist and are well coaligned with the threads.
Very-high-energy $\gamma$-ray observations of the active galaxy IC 310 with the MAGIC telescopes have revealed fast variability with doubling time scales of less than 4.8min. This implies that the emission region in IC 310 is smaller than 20% of the gravitational radius of the central supermassive black hole with a mass of $3\times 10^8 M_\odot$, which poses serious questions on the emission mechanism and classification of this enigmatic object. We report on the first quasi-simultaneous multi-frequency VLBI observations of IC 310 conducted with the EVN. We find a blazar-like one-sided core-jet structure on parsec scales, constraining the inclination angle to be less than $\sim 20^\circ$ but very small angles are excluded to limit the de-projected length of the large-scale radio jet.
We use a volume-limited sample of quasars in the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog to identify quasar groups and address their statistical significance. This quasar sample has a uniform selection function on the sky and nearly a maximum possible contiguous volume that can be drawn from the DR7 catalog. Quasar groups are identified by using the Friend-of-Friend algorithm with a set of fixed comoving linking lengths. We find that the richness distribution of the richest 100 quasar groups or the size distribution of the largest 100 groups are statistically equivalent with those of randomly-distributed points with the same number density and sky coverage when groups are identified with the linking length of 70 h-1Mpc. It is shown that the large-scale structures like the huge Large Quasar Group (U1.27) reported by Clowes et al. (2013) can be found with high probability even if quasars have no physical clustering, and does not challenge the initially homogeneous cosmological models. Our results are statistically more reliable than those of Nadathur (2013), where the test was made only for the largest quasar group. It is shown that the linking length should be smaller than 50 h-1Mpc in order for the quasar groups identified in the DR7 catalog not to be dominated by associations of quasars grouped by chance. We present 20 richest quasar groups identified with the linking length of 70 h-1Mpc for further analyses.
Atmospheres with a high C/O ratio are expected to contain an important quantity of hydrocarbons, including heavy molecules (with more than 2 carbon atoms). To study correctly these C-rich atmospheres, a chemical scheme adapted to this composition is necessary. We have implemented a chemical scheme that can describe the kinetics of species with up to 6 carbon atoms. This chemical scheme has been developed with specialists of combustion and validated through experiments on a wide range of T and P. This chemical network is available on the online database KIDA. We have created a grid of 12 models to explore different thermal profiles and C/O ratios. For each of them, we have compared the chemical composition determined with a C0-C2 chemical scheme (species with up to 2 carbon atoms) and with the C0-C6 scheme. We found no difference in the results obtained with the two schemes when photolyses are not included in the model, whatever the temperature of the atmosphere. In contrast, when there is photochemistry, differences can appear in the upper atmosphere. These differences are found for all the tested PT profiles in the case that the C/O ratio is above 1. When the C/O ratio of the atmosphere is solar, differences are only found at temperatures lower than 1000K. The differences linked to the use of different chemical schemes do not have important influence on the synthetic spectra. However, we have confirmed that C2H2 and HCN as possible tracers of warm C-rich atmospheres. The use of this new chemical scheme is mandatory to model atmospheres with a high C/O ratio and, in particular, if one is interested in studying in details the photochemistry. If one is just interested in the synthetic spectra, the use of a smaller scheme may be sufficient.
Microflares are small dynamic signatures observed in X-ray and extreme-ultraviolet channels. Because of their impulsive emission enhancements and wide distribution, they are thought to be closely related to coronal heating. By using the high resolution 171 {\AA} images from the Atmospheric Imaging Assembly and the lines-of-sight magnetograms obtained by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, we trace 10794 microflares in a quiet region near the disk center with a field of view of 960 arcsec $\times$ 1068 arcsec during 24 hr. The microflares have an occurrence rate of 4.4 $\times$ 10$^{3}$ hr$^{-1}$ extrapolated over the whole Sun. Their average brightness, size, and lifetime are 1.7 I$_{0}$(of the quiet Sun), 9.6 Mm$^{2}$, and 3.6 min, respectively. There exists a mutual positive correlation between the microflares' brightness, area and lifetime. In general, the microflares distribute uniformly across the solar disk, but form network patterns locally, which are similar to and matched with the magnetic network structures. Typical cases show that the microflares prefer to occur in magnetic cancellation regions of network boundaries. We roughly calculate the upper limit of energy flux supplied by the microflares and find that the result is still a factor of $\sim$15 below the coronal heating requirement.
We compare the stellar motion around a spiral arm created in two different scenarios, transient/co-rotating spiral arms and density-wave-like spiral arms. We generate Gaia mock data from snapshots of the simulations following these two scenarios using our stellar population code, SNAPDRAGONS, which takes into account dust extinction and the expected Gaia errors. We compare the observed rotation velocity around a spiral arm similar in position to the Perseus arm, and find that there is a clear difference in the velocity features around the spiral arm between the co-rotating spiral arm and the density-wave-like spiral arm. Our result demonstrates that the volume and accuracy of the Gaia data are sufficient to clearly distinguish these two scenarios of the spiral arms.
The era of high-precision astrometry has dawned upon us. The potential of Gaia $\mu$as-level precision in positional measurements is about to be unleashed in the field of extrasolar planetary systems. The Gaia data hold the promise for much improved global characterization of planetary systems around stars of all types, ages, and chemical composition, particularly when synergistically combined with other indirect and direct planet detection and characterization programs.
The Atacama Large Millimeter/submillimeter Array offers an unprecedented view of our Sun at sub-/millimeter wavelengths. The high spatial, temporal, and spectral resolution facilitates the measurement of gas temperatures and magnetic fields in the solar chromosphere with high precision. The anticipated results will revolutionize our understanding of the solar atmosphere and may in particular result in major steps towards solving the coronal heating problem. Based on state-of-the-art 3D radiation magnetohydrodynamic simulations, we calculate the emergent continuum intensity (and thus brightness temperature maps) in the wavelength range accessed by ALMA and simulate instrumental effects for different array configurations. First results show that the local gas temperature can be closely mapped with ALMA and that much of the complex small-scale chromospheric pattern can be resolved.
We have investigated the thermally induced proton/deuteron exchange in mixed amorphous H$_2$O:D$_2$O ices by monitoring the change in intensity of characteristic vibrational bending modes of H$_2$O, HDO, and D$_2$O with time and as function of temperature. The experiments have been performed using an ultra-high vacuum setup equipped with an infrared spectrometer that is used to investigate the spectral evolution of homogeneously mixed ice upon co-deposition in thin films, for temperatures in the 90 to 140 K domain. With this non-energetic detection method we find a significantly lower activation energy for H/D exchange -- $3840 \pm 125$ K -- than previously reported. Very likely this is due to the amorphous nature of the interstellar ice analogues involved. This provides reactive timescales ($\tau<10^4$ years at $T$ $>70$ K) fast enough for the process to be important in interstellar environments. Consequently, an astronomical detection of D$_2$O will be even more challenging because of its potential to react with H$_2$O to form HDO. Furthermore, additional experiments, along with previous studies, show that proton/deuteron swapping also occurs in ice mixtures of water with other hydrogen bonded molecules, in particular on the OH and NH moieties. We conclude that H/D exchange in ices is a more general process that should be incorporated into ice models that are applied to protoplanetary disks or to simulate the warming up of cometary ices in their passage of the perihelion, to examine the extent of its influence on the final deuteron over hydrogen ratio.
We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]
We present the light curve data of a remarkable blazer 3C 454.3 (z=0.859) in optical, X-ray, and gamma-ray bands. Since January 2008, we have been monitoring this object using the 50 cm MITSuME, a optical telescope, and detected several flares including extraordinary and simultaneous flares in the $\gamma$-ray and optical bands in November 2010. Additionally, the Monitor of All-sky Image (MAXI) has been observing 3C 454.3 continuously since August 2009. Using these data and gamma-ray flux observed with Fermi-LAT, we discuss features and correlations of flux variations between the energy bands.
We investigated the formation of arc-like structures in the infalling envelope around protostars, motivated by the recent ALMA observations of the high-density molecular cloud core, MC27/L1527F. We performed self-gravitational hydrodynamical numerical simulations with an adaptive mesh refinement code. A filamentary cloud with a 0.1~pc width fragments into cloud cores because of perturbations due to weak turbulence. The cloud core undergoes gravitational collapse to form multiple protostars, and gravitational torque from the orbiting protostars produces arc structures extending up to a 1000~AU scale. As well as on a spatial extent, the velocity ranges of the arc structures, $\sim0.5\,\mathrm{km\,s}^{-1}$, are in agreement with the ALMA observations. We also found that circumstellar disks are often misaligned in triple system. The misalignment is caused by the tidal interaction between the protostars when they undergo close encounters because of a highly eccentric orbit of the tight binary pair.
The dissipation of kinetic and magnetic energy in the interstellar medium (ISM) can proceed through viscous, Ohmic or ambipolar diffusion (AD). It occurs at very small scales compared to the scales at which energy is presumed to be injected. This localized heating may impact the ISM evolution but also its chemistry, thus providing observable features. Here, we perform 3D spectral simulations of decaying magnetohydrodynamic turbulence including the effects of AD. We find that the AD heating power spectrum peaks at scales in the inertial range, due to a strong alignment of the magnetic and current vectors in the dissipative range. AD affects much greater scales than the AD scale predicted by dimensional analysis. We find that energy dissipation is highly concentrated on thin sheets. Its probability density function follows a lognormal law with a power-law tail which hints at intermittency, a property which we quantify by use of structure function exponents. Finally, we extract structures of high dissipation, defined as connected sets of points where the total dissipation is most intense and we measure the scaling exponents of their geometric and dynamical characteristics: the inclusion of AD favours small sizes in the dissipative range.
We report on recent Ks-band interferometric observations of the young pre-main-sequence star AB Aurigae obtained with the Palomar Fiber Nuller (PFN). Reaching a contrast of a few 10-4 inside a field of view extending from 35 to 275 mas (5-40AU at AB Aur's distance), the PFN is able to explore angular scales that are intermediate between those accessed by coronagraphic imaging and long baseline interferometry. This intermediate region is of special interest given that many young stellar objects are believed to harbor extended halos at such angular scales. Using destructive interference (nulling) between two sub-apertures of the Palomar 200 inch telescope and rotating the telescope pupil, we measured a resolved circumstellar excess at all probed azimuth angles. The astrophysical null measured over the full rotation is fairly constant, with a mean value of 1.52%, and a slight additional azimuthal modulation of +/-0.2%. The isotropic astrophysical null is indicative of circumstellar emission dominated by an azimuthally extended source, possibly a halo, or one or more rings of dust, accounting for several percent of the total Ks-band flux. The modest azimuthal variation may be explained by some skewness or anisotropy of the spatially extended source, e.g., an elliptical or spiral geometry, or clumping, but it could also be due to the presence of a point source located at a separation of ~120 mas (17AU) with ~6*10-3 of the stellar flux. We combine our results with previous Infrared Optical Telescope Array observations of AB Aur at H band, and demonstrate that a dust ring located at ~30 mas (4.3AU) represents the best-fitting model to explain both sets of visibilities.We are also able to test a few previously hypothesized models of the incoherent component evident at longer interferometric baselines.
Episodic jets are usually observed in the intermediate state of black hole transients during their X-ray outbursts. Here we report the discovery of a strong positive correlation between the peak radio power of the episodic jet $P_{\rm jet}$ and the corresponding peak X-ray luminosity $L_{\rm x}$ of the soft state (in Eddington units) in a complete sample of the outbursts of black hole transients observed during the RXTE era of which data are available, which follows the relation $\log P_{\rm jet}=(2.17\pm{0.32})+(1.63\pm0.24)\times \log {L_{\rm x}}$. The transient ultra-luminous X-ray source in M31 and HLX-1 in EXO 243-49 fall on the relation if they contain stellar mass black hole and either stellar mass black hole or intermediate mass black hole, respectively. Besides, a significant correlation between the peak power of the episodic jet and the rate-of-increase of the X-ray luminosity $\rm dL_{x}/dt$ during the rising phase of those outbursts is also found, following $\log P_{\rm jet}=(1.97\pm{0.42})+(0.69\pm0.15)\times \log {\rm d}L_{\rm x}/{\rm d}t$. In GX 339$-$4 and H 1743$-$322 in which data for two outbursts are available, measurements of the peak radio power of the episodic jet and the X-ray peak luminosity (and its rate-of-change) shows similar positive correlations between outbursts, which demonstrate the dominant role of accretion over black hole spin in generating episodic jet power. On the other hand, no significant difference is seen among the systems with different measured black hole spin in current sample. This implies that the power of the episodic jet is strongly affected by non-stationary accretion characterised primarily by the rate-of-change of the mass accretion rate.
Routinely operating since July 2012, the APACHE survey has celebrated its second birthday. While the main goal of the Project is the detection of transiting planets around a large sample of bright, nearby M dwarfs in the northern hemisphere, the APACHE large photometric database for hundreds of different fields represents a relevant resource to search for and provide a first characterization of new variable stars. We celebrate here the conclusion of the second year of observations by reporting the discovery of 14 new variables.
We report the results of the spectroscopic observations carried out at the SAO RAS 6-m telescope for the optical components of nine new extended radio sources found in the NVSS catalog. The measured redshifts of the host galaxies are in the range of z=0.1-0.4. The physical sizes of radio sources were calculated within the standard cosmological model. The two most extended objects, 0003+1512 and 0422+0351 reach the sizes of 2.1 Mpc and 4.0 Mpc, respectively. This is close to the maximum size of known radio sources.
Wide-Field MAXI (WF-MAXI: Wide-Field Monitor of All-sky X-ray Image) is a proposed mission to detect and localize X-ray transients including electro-magnetic counterparts of gravitational-wave events such as gamma-ray bursts and supernovae etc., which are expected to be directly detected for the first time in late 2010's by the next generation gravitational telescopes such as Advanced LIGO and KAGRA. The most distinguishing characteristics of WF-MAXI are a wide energy range from 0.7 keV to 1 MeV and a large field of view (~25 % of the entire sky), which are realized by two main instruments: (i) Soft X-ray Large Solid Angle Camera (SLC) which consists of four pairs of crisscross coded aperture cameras using CCDs as one-dimensional fast-readout detectors covering 0.7 - 12 keV and (ii) Hard X-ray Monitor (HXM) which is a multi-channel array of crystal scintillators coupled with avalanche photo-diodes covering 20 keV - 1 MeV.
We investigate the circumstellar disc fraction as determined from L-band excess observations of the young, massive Arches and Quintuplet clusters residing in the central molecular zone of the Milky Way. The Quintuplet cluster was searched for L-band excess sources for the first time. We find a total of 26 excess sources in the Quintuplet cluster and 21 in the Arches cluster, of which 13 are new detections. With the aid of proper motion membership samples, the disc fraction of the Quintuplet cluster was derived for the first time to be 4.0 +/- 0.7%. There is no evidence for a radially varying disc fraction in this cluster. In the case of the Arches cluster, a disc fraction of 9.2 +/- 1.2% approximately out to the cluster's predicted tidal radius, r < 1.5 pc, is observed. This excess fraction is consistent with our previously found disc fraction in the cluster in the radial range 0.3 < r < 0.8 pc. In both clusters, the host star mass range covers late A- to early B-type stars, 2 < M < 15 Msun, as derived from J-band photospheric magnitudes. We discuss the unexpected finding of dusty circumstellar discs in these UV intense environments in the context of primordial disc survival and formation scenarios of secondary discs. We consider the possibility that the L-band excess sources in the Arches and Quintuplet clusters could be the high-mass counterparts to T Tauri pre-transitional discs. As such a scenario requires a long pre-transitional disc lifetime in a UV intense environment, we suggest that mass transfer discs in binary systems are a likely formation mechanism for the B-star discs observed in these starburst clusters.
We present a backward approach for the interpretation of the evolution of the near-infrared and the far-infrared luminosity functions across the redshift range 0<z<3. In our method, late-type galaxies are treated by means of a parametric phenomenological method based on PEP/HerMES data up to z~4, whereas spheroids are described by means of a physically motivated backward model. The spectral evolution of spheroids is modelled by means of a single-mass model, associated to a present-day elliptical with K-band luminosity comparable to the break of the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr across the redshift range 1< z < 5. The key parameter is represented by the redshift z_0.5 at which half proto-spheroids are already formed. A statistical study indicates for this parameter values between z_0.5=1.5 and z_0.5=3. We assume as fiducial value z_0.5~2, and show that this assumption allows us to describe accourately the redshift distributions and the source counts. By assuming z_0.5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources observed at z>2 can be explained as progenitors of local spheroids caught during their formation. We also test the effects of mass downsizing by dividing the spheroids into three populations of different present-day stellar masses. The results obtained in this case confirm the validity of our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model which describes the evolution of galaxies at the break of the present-day early type K-band LF.
For over a century, light echoes have been observed around variable stars and transients. The discovery of centuries-old light echoes from supernovae in the Large Magellanic Cloud has allowed the spectroscopic characterization of these events using modern instrumentation, even in the complete absence of any visual record of those events. Here we review the pivotal role the Blanco 4m telescope played in these discoveries.
In 'modified' gravity the observed acceleration of the universe is explained by changing the gravitational force law or the number of degrees of freedom in the gravitational sector. Both possibilities can be tested by measurements of cosmological structure formation. In this paper we elaborate the details of such tests using the Galileon model as a case study. We pay attention to the possibility that each new degree of freedom may have stochastically independent initial conditions, generating different types of potential well in the early universe and breaking complete correlation between density and velocity power spectra. This 'stochastic bias' can confuse schemes to parametrize the predictions of modified gravity models, such as the use of the growth parameter f alone. Using data from the WiggleZ Dark Energy Survey we show that it will be possible to obtain constraints using information about the cosmological-scale force law embedded in the multipole power spectra of redshift-space distortions. As an example, we obtain an upper limit on the strength of the conformal coupling to matter in the cubic Galileon model, giving |1/M| < 200 / Mp. This allows the fifth-force to be stronger than gravity, but is consistent with zero coupling.
Obscured AGNs provide an opportunity to study the material surrounding the central engine. Geometric and physical constraints on the absorber can be deduced from the reprocessed AGN emission. In particular, the obscuring gas may reprocess the nuclear X-ray emission producing a narrow Fe K$\alpha$ line and a Compton reflection hump. In recent years, models of the X-ray reflection from an obscuring torus have been computed; however, although the reflecting gas may be dusty, the models do not yet take into account the effects of dust on the predicted spectrum. We study this problem by analyzing two sets of models, with and without the presence of dust, using the one dimensional photo-ionization code Cloudy. The calculations are performed for a range of column densities ($22 <{\rm log}[N_H(\rm cm^{-2})]< 24.5$ ) and hydrogen densities ( $6 <{\rm log}[n_H(\rm cm^{-3})]< 8$). The calculations show the presence of dust can enhance the Fe K$\alpha$ equivalent width (EW) in the reflected spectrum by factors up to $\approx$ 8 for Compton thick (CT) gas and a typical ISM grain size distribution. The enhancement in EW with respect to the reflection continuum is due to the reduction in the reflected continuum intensity caused by the anisotropic scattering behaviour of dust grains. This effect will be most relevant for reflection from distant, predominately neutral gas, and is a possible explanation for AGNs which show a strong Fe K$\alpha$ EW and a relatively weak reflection continuum. Our results show it is an important to take into account dust while modeling the X-ray reflection spectrum, and that inferring a CT column density from an observed Fe K$\alpha$ EW may not always be valid. Multi-dimensional models are needed to fully explore the magnitude of the effect.
The DANCe survey provides photometric and astrometric (position and proper
motion) measurements for approximately 2 millions unique sources in a region
encompassing $\approx$80deg$^{2}$ centered around the Pleiades cluster.
We aim at deriving a complete census of the Pleiades, and measure the mass
and luminosity function of the cluster. Using the probabilistic selection
method described in Sarro+2014, we identify high probability members in the
DANCe ($i\ge$14mag) and Tycho-2 ($V\lesssim$12mag) catalogues, and study the
properties of the cluster over the corresponding luminosity range. We find a
total of 2109 high probability members, of which 812 are new, making it the
most extensive and complete census of the cluster to date. The luminosity and
mass functions of the cluster are computed from the most massive members down
to $\approx$0.025M$_{\odot}$. The size, sensitivity and quality of the sample
result in the most precise luminosity and mass functions observed to date for a
cluster. Our census supersedes previous studies of the Pleiades cluster
populations, both in terms of sensitivity and accuracy.
Dark Matter (DM) is a fundamental ingredient of our Universe and of structure formation, and yet its nature is elusive to astrophysical probes. Information on the nature and physical properties of the WIMP (neutralino) DM (the leading candidate for a cosmologically relevant DM) can be obtained by studying the astrophysical signals of their annihilation/decay. Among the various e.m. signals, secondary electrons produced by neutralino annihilation generate synchrotron emission in the magnetized atmosphere of galaxy clusters and galaxies which could be observed as a diffuse radio emission (halo or haze) centered on the DM halo. A deep search for DM radio emission with SKA in local dwarf galaxies, galaxy regions with low star formation and galaxy clusters (with offset DM-baryonic distribution, like e.g. the Bullet cluster) can be very effective in constraining the neutralino mass, composition and annihilation cross-section. For the case of a dwarf galaxy, like e.g. Draco, the constraints on the DM annihilation cross-section obtainable with SKA1-MID will be at least a factor $\sim 10^3$ more stringent than the limits obtained by Fermi-LAT in the $\gamma$-rays. These limits scale with the value of the B field, and the SKA will have the capability to determine simultaneously both the magnetic field in the DM-dominated structures and the DM particle properties. The optimal frequency band for detecting the DM-induced radio emission is around $\sim 1$ GHz, with the SKA1-MID Band 1 and 4 important to probe the synchrotron spectral curvature at low-$\nu$ (sensitive to DM composition) and at high-$\nu$ (sensitive to DM mass).
We have used the broadband backend available at the ATCA to study the fast interstellar scintillation of quasar PKS 1257-326, resolving the core shift as a function of frequency on scales less than 10 microarcseconds. In this short paper we discuss the jet direction implied from the microarcsecond-scale core shift in PKS 1257-326.
Extragalactic jets are a common feature of radio-loud active galaxies. The nature of the observed jets in relation to the bulk flow is still unclear. In particular it is not clear whether the observations of parsec-scale jets using the very long baseline interferometric technique (VLBI) reveal wave-like structures that develop and propagate along the jet, or trace the jet flow itself. In this contribution I review the evidence collected during the last years showing that the ridge-lines of helical radio-jets do not correspond to observational artifacts. This conclusion was reached by studying a number of VLBI observations of the radio jet in the quasar S5~0836+710 at different frequencies and epochs. The ridge-line of the emission in the jet coincides at all frequencies within the errors. Moreover, small differences between the ridge-lines as observed at different epochs reveal wave-like motion transversal to the jet propagation axis. I also discuss similar results, albeit with different interpretations, obtained by other authors. The current challenge is to measure the propagation velocities of these waves and to try to characterise them in terms of simple perturbations or Kelvin-Helmholtz instability, which would help understanding the physical conditions of the flow where the waves develop. This problem can only be tackled by high-resolution observations such as those that can be achieved by the space radio-antenna Radioastron.
Microscopic fluctuations inherent to the fuzziness of spacetime at the Planck scale might accumulate in wavefronts propagating a cosmological distance and lead to noticeable blurring in an image of a pointlike source. Distant quasars viewed in the optical and ultraviolet with Hubble Space Telescope (HST} may show this weakly, and if real suggests a stronger effect should be seen for Gamma-Ray Bursts (GRBs) in X-rays and gamma-rays. Those telescopes, however, operate far from their diffraction limits. A description of how Planck-scale-induced blurring could be sensed at high energy, even with cosmic rays, while still agreeing with the HST results is discussed. It predicts dilated apparent source size and inflated uncertainties in positional centroids, effectively a threshold angular accuracy restricting knowledge of source location on the sky. These outcomes are found to be consistent with an analysis of the 10 highest-redshift GRB detections reported for the Fermi satellite. Confusion with photon cascade and scattering phenomena is also possible; prospects for a definitive multiwavelength measurement are considered.
In this paper we present the Adaptive Optics Lucky Imager (AOLI), a
state-of-the-art instrument which makes use of two well proved techniques for
extremely high spatial resolution with ground-based telescopes: Lucky Imaging
(LI) and Adaptive Optics (AO).
AOLI comprises an AO system, including a low order non-linear curvature
wavefront sensor together with a 241 actuators deformable mirror, a science
array of four 1024x1024 EMCCDs, allowing a 120x120 down to 36x36 arcseconds
field of view, a calibration subsystem and a powerful LI software. Thanks to
the revolutionary WFS, AOLI shall have the capability of using faint reference
stars ({\it I\/} $\sim$ 16.5-17.5), enabling it to be used over a much wider
part of the sky than with common Shack-Hartmann AO systems.
This instrument saw first light in September 2013 at William Herschel
Telescope. Although the instrument was not complete, these commissioning
demonstrated its feasibility, obtaining a FWHM for the best PSF of
0.151$\pm$0.005 arcsec and a plate scale of 55.0$\pm$0.3 mas/pixel. Those
observations served us to prove some characteristics of the interesting
multiple T Tauri system LkH$\alpha$ 262-263, finding it to be gravitationally
bounded. This interesting multiple system mixes the presence of proto-planetary
discs, one proved to be double, and the first-time optically resolved pair
LkH$\alpha$ 263AB (0.42 arcsec separation).
We present the results of the pulse phase- and luminosity-resolved spectroscopy of the transient X-ray pulsar V0332+53, performed for the first time in a wide luminosity range (1-40)x10^{37} erg/s during a giant outburst observed by the RXTE observatory in Dec 2004 - Feb 2005. We characterize the spectra quantitatively and built the detailed "three-dimensional" picture of spectral variations with pulse phase and throughout the outburst. We show that all spectral parameters are strongly variable with the pulse phase, and the pattern of this variability significantly changes with luminosity directly reflecting the associated changes in the structure of emission regions and their beam patterns. Obtained results are qualitatively discussed in terms of the recently developed reflection model for the formation of cyclotron lines in the spectra of X-ray pulsars.
Glitches are sudden changes in rotation frequency and spin-down rate, observed from pulsars of all ages. Standard glitches are characterized by a positive step in angular velocity ($\Delta\Omega$ $ > $ $0$) and a negative step in the spin-down rate ($\Delta \dot \Omega$ $ < $ $0$) of the pulsar. There are no glitch-associated changes in the electromagnetic signature of rotation-powered pulsars in all cases so far. For the first time, in the last glitch of PSR J1119-6127, there is clear evidence for changing emission properties coincident with the glitch. This glitch is also unusual in its signature. Further, the absolute value of the spin-down rate actually decreases in the long term. This is in contrast to usual glitch behaviour. In this paper we extend the vortex creep model in order to take into account these peculiarities. We propose that a starquake with crustal plate movement towards the rotational poles of the star induces inward vortex motion which causes the unusual glitch signature. The component of the magnetic field perpendicular to the rotation axis will decrease, giving rise to a permanent change in the pulsar external torque.
We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $\alpha = -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} \Psi^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $\mu>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.
We present observations and initial analysis from an HST/STIS program to obtain the first co-spatial, UV-optical spectra of ten Galactic planetary nebulae (PNe). Our primary objective was to measure the critical emission lines of carbon and nitrogen with unprecedented S/N and spatial resolution over UV-optical range, with the ultimate goal of quantifying the production of these elements in low- and intermediate-mass stars. Our sample was selected from PNe with a near-solar metallicity, but spanning a broad range in N/O. This study, the first of a series, concentrates on the observations and emission-line measurements obtained by integrating along the entire spatial extent of the slit. We derived ionic and total elemental abundances for the seven PNe with the strongest UV line detections (IC~2165, IC~3568, NGC~2440, NGC~3242, NGC~5315, NGC~5882, and NGC~7662). We compare these new results with other recent studies of the nebulae, and discuss the relative merits of deriving the total elemental abundances of C, N, and O using ionization correction factors (ICFs) versus summed abundances. For the seven PNe with the best UV line detections, we conclude that summed abundances from direct diagnostics of ions with measurable UV lines gives the most accurate values for the total elemental abundances of C and N. In some cases where significant discrepancies exist between our abundances and those from other studies, we show that the differences can often be attributed to their use of fluxes that are not co-spatial. Finally, we examined C/O and N/O versus O/H and He/H in well-observed Galactic, LMC, and SMC PNe, and found that highly accurate abundances are essential for properly inferring elemental yields from their progenitor stars.
We present wide-field, deep narrowband H$_2$, Br$\gamma$, H$\alpha$, [S II], [O III], and broadband I and K-band images of the Carina star formation region. The new images provide a large-scale overview of all the H$_2$ and Br$\gamma$ emission present in over a square degree centered on this signature star forming complex. By comparing these images with archival HST and Spitzer images we observe how intense UV radiation from O and B stars affects star formation in molecular clouds. We use the images to locate new candidate outflows and identify the principal shock waves and irradiated interfaces within dozens of distinct areas of star-forming activity. Shocked molecular gas in jets traces the parts of the flow that are most shielded from the intense UV radiation. Combining the H$_2$ and optical images gives a more complete view of the jets, which are sometimes only visible in H$_2$. The Carina region hosts several compact young clusters, and the gas within these clusters is affected by radiation from both the cluster stars and the massive stars nearby. The Carina Nebula is ideal for studying the physics of young H II regions and PDR's, as it contains multiple examples of walls and irradiated pillars at various stages of development. Some of the pillars have detached from their host molecular clouds to form proplyds. Fluorescent H$_2$ outlines the interfaces between the ionized and molecular gas, and after removing continuum, we detect spatial offsets between the Br$\gamma$ and H$_2$ emission along the irradiated interfaces. These spatial offsets can be used to test current models of PDRs once synthetic maps of these lines become available.
This work presents an extensive study of the previously discovered formation of bipolar flux concentrations in a two-layer model. We relate the formation process to the negative effective magnetic pressure instability (NEMPI), which is a possible mechanism to explain the origin of sunspots. In our simulations we use a Cartesian domain of isothermal stratified gas which is divided into two layers. In the lower layer, turbulence is forced with transverse non-helical random waves, whereas in the upper layer no flow is induced. An initially weak uniform horizontal magnetic field is imposed in the entire domain. In this study we vary the stratification by changing the gravitational acceleration, magnetic Reynolds number, the strength of the imposed magnetic field and the size of the domain to investigate their influence on the formation process. Bipolar magnetic structure formation takes place over a large range of parameters. The magnetic structures become more intensive for higher stratification. The large fluid Reynolds numbers allow for the generation of flux concentrations when the magnetic Prandtln umber is between 0.1 and 1. The magnetic field in bipolar regions increases with higher imposed field strength until the field becomes comparable to the equipartition field strength of the turbulence. A larger horizontal extent enables the flux concentrations to become stronger and more coherent. The size of the bipolar structures turns out to be independent of the domain size. Bipolar flux concentrations are correlated with strong large-scale downward and converging flows and can therefore be explained by NEMPI.
We show that a positive signal in a dark matter (DM) direct detection experiment can be used to place a lower bound on the DM capture rate in the Sun, independent of the DM halo. For a given particle physics model and DM mass we obtain a lower bound on the capture rate independent of the local DM density, velocity distribution, galactic escape velocity, as well as the scattering cross section. We illustrate this lower bound on the capture rate by assuming that upcoming direct detection experiments will soon obtain a significant signal. When comparing the lower bound on the capture rate with limits on the high-energy neutrino flux from the Sun from neutrino telescopes, we can place upper limits on the branching fraction of DM annihilation channels leading to neutrinos. With current data from IceCube and Super-Kamiokande non-trivial limits can be obtained for spin-dependent interactions and direct annihilations into neutrinos. In some cases also annihilations into $\tau\tau$ or $b\bar b$ start getting constrained. For spin-independent interactions current constraints are weak, but they may become interesting for data from future neutrino telescopes.
We perform new longterm (15-16 orbits) simulations of coalescing binary neutron stars in numerical relativity using an updated Einstein's equation solver, employing low-eccentricity initial data, and modeling the neutron stars by a piecewise polytropic equation of state. A convergence study shows that our new results converge more rapidly than the third order and using the determined convergence order, we construct an extrapolated waveform for which the estimated total phase error should be less than 1 radian. We then compare the extrapolated waveforms with those calculated by the latest effective-one-body (EOB) formalism in which the so-called tidal deformability, higher post-Newtonian corrections, and gravitational self-force effects are taken into account. We show that for a binary of compact neutron stars with their radius 11.1 km, the waveform by the EOB formalism agrees quite well with the numerical waveform so that the total phase error is smaller than 1 radian for the total phase of $\sim 200$ radian up to the merger. By contrast, for a binary of less compact neutron stars with their radius 13.6 km, the EOB and numerical waveforms disagree with each other in the last few wave cycles, resulting in the total phase error of $\sim 3$ radian.
The universe may contain several decoupled matter sectors which primarily couple through gravity to the Standard Model degrees of freedom. We focus here on the description of astrophysical environments that allow for comparable densities and spatial distributions of visible matter and decoupled dark matter. We discuss four Wolf-Rayet galaxies (NGC 1614, NGC 3367, NGC 4216 and NGC 5430) which should contain comparable amounts of decoupled dark and visible matter in the star forming regions. This could lead to the observation of Gamma Ray Burst events with physics modified by jets of dark matter radiation.
We investigate the Affleck-Dine baryogenesis after D-term inflation with a positive Hubble-induced mass term for a B-L flat direction. It stays at a large field value during D-term inflation, and just after inflation ends it starts to oscillate around the origin of the potential due to the positive Hubble-induced mass term. The phase direction is kicked by higher-dimensional Kahler potentials to generate the B-L asymmetry. The scenario predicts nonzero baryonic isocurvature perturbations, which would be detected by future observations of CMB fluctuations. We also provide a D-term inflation model which naturally explain the coincidence of the energy density of baryon and dark matter.
We investigate two classes of models of quintessential inflation, based upon canonical as well as noncanonical scalar fields. In particular, introducing potentials steeper than the standard exponential, we construct models that can give rise to a successful inflationary phase, with signatures consistent with Planck 2015 results. Additionally, using nonminimal coupling of the scalar field with massive neutrino matter, we obtain the standard thermal history of the Universe, with late-time cosmic acceleration as the last stage of evolution. In both cases, inflation and late-time acceleration are connected by a tracker solution.
The Sun and the undisturbed interstellar magnetic field and plasma velocity vectors (Bis,Vis) define a mirror symmetry plane of the flow at large heliospheric distances. We show that for the Bis direction defined by IBEX Ribbon center, the radial direction of Voyager 2 over the last decade, and the (thermal proton) plasma velocity measured by the spacecraft since 2010.5, are almost parallel to the (Bis,Vis)-plane, which coincides in practice with the Hydrogen Deflection Plane. These facts can be simply explained if approximate mirror symmetry is also maintained on the inner side of the heliopause. Such approximate symmetry is possible since the solar wind ram pressure is almost spherically symmetric and the plasma beta value in the inner heliosheath is high. In the proposed symmetry, the plasma flow speed measured by Voyager 2 in the inner heliosheath is expected to rotate more in the transverse than in the polar direction (explanation alternative to McComas & Schwadron (2014)), in evident agreement with available spacecraft data (our Fig.1).
Gravitational waves can act like gravitational lenses, affecting the observed positions, brightnesses, and redshifts of distant objects. Exact expressions for such effects are derived here, allowing for arbitrarily-moving sources and observers in the presence of plane-symmetric gravitational waves. The commonly-used predictions of linear perturbation theory are shown to be generically overshadowed---even for very weak gravitational waves---by nonlinear effects when considering observations of sufficiently distant sources; higher-order perturbative corrections involve secularly-growing terms which cannot necessarily be neglected. Even on more moderate scales where linear effects remain at least marginally dominant, nonlinear corrections are qualitatively different from their linear counterparts. There is a sense in which they can, for example, mimic the existence of a third type of gravitational wave polarization.
The effective potential for the Standard Model Higgs field allows two quasi-degenerate vacua; one is our vacuum at the electroweak scale, while the other is at a much higher scale. The latter minimum may be at a scale much smaller than the Planck scale, if the potential is lifted by new physics. This gives rise to a possibility of domain wall formation after inflation. If the high-scale minimum is a local minimum, domain walls are unstable and disappear through violent annihilation processes, producing a significant amount of gravitational waves. We estimate the amount of gravitational waves produced from unstable domain walls in the Higgs potential and discuss detectability with future experiments.
If two DGP branes carry U(1) gauge theories and overlap, particles of one
brane can interact with the photons from the other brane. This coupling
modifies in particular the Coulomb potentials between charges from the same
brane in the overlapping regions. The coupling also introduces Coulomb
interactions between charges from the different branes which can generate
exotic bound states.
The effective modification of the fine structure constant in the overlap
region generates a trough in signals at the redshift of the overlap region and
an increase at smaller or larger redshift, depending on the value of the
crosstalk parameter. This implies potentially observable perturbations in the
Lyman-alpha forest if our 3-brane overlapped with another 3-brane in a region
with redshift z<6. Crosstalk can also affect structure formation by enhancing
or suppressing radiative cooling.
We hypothesize that in a non-metallic crystalline structure under extreme pressures, atomic wavefunctions deform to adopt a reduced rotational symmetry consistent with minimizing interstitial space in the crystal. We exemplify with a simple numeric variational calculation that yields the energy cost of this deformation for Helium to 25%. Balancing this with the free energy gained by tighter packing we obtain the pressures required to effect such deformation. The consequent modification of the structure suggests a decrease in the resistance to tangential stress, and an associated decrease of the crystal's shear modulus. The atomic form factor is also modified. We also compare with neutron matter in the interior of compact stars.
We have studied channeling effects in a Cesium Iodide (CsI) crystal that is similar in composition to the ones being used in a search for Weakly Interacting Massive Particles (WIMPs) dark matter candidates, and measured its energy-dependent quenching factor, the relative scintillation yield for electron and nuclear recoils. The experimental results are reproduced with a GEANT4 simulation that includes a model of the scintillation efficiency as a function of electronic stopping power. We present the measured and simulated quenching factors and the estimated effects of channeling.
Interstellar is the first Hollywood movie to attempt depicting a black hole
as it would actually be seen by somebody nearby. For this we developed a code
called DNGR (Double Negative Gravitational Renderer) to solve the equations for
ray-bundle (light-beam) propagation through the curved spacetime of a spinning
(Kerr) black hole, and to render IMAX-quality, rapidly changing images. Our
ray-bundle techniques were crucial for achieving IMAX-quality smoothness
without flickering.
This paper has four purposes: (i) To describe DNGR for physicists and CGI
practitioners . (ii) To present the equations we use, when the camera is in
arbitrary motion at an arbitrary location near a Kerr black hole, for mapping
light sources to camera images via elliptical ray bundles. (iii) To describe
new insights, from DNGR, into gravitational lensing when the camera is near the
spinning black hole, rather than far away as in almost all prior studies. (iv)
To describe how the images of the black hole Gargantua and its accretion disk,
in the movie \emph{Interstellar}, were generated with DNGR. There are no new
astrophysical insights in this accretion-disk section of the paper, but disk
novices may find it pedagogically interesting, and movie buffs may find its
discussions of Interstellar interesting.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)