We report on NuSTAR and Swift observations of a soft state of the neutron star low-mass X-ray binary GS 1826-24, commonly known as the "clocked" burster. The transition to the soft state was recorded in 2014 June through an increase of the 2-20 keV source intensity measured by MAXI, simultaneous with a decrease of the 15-50 keV intensity measured by Swift/BAT. The episode lasted approximately two months, after which the source returned to its usual hard state. We analyze the broad-band spectrum measured by Swift/XRT and NuSTAR, and estimate the accretion rate during the soft episode to be about 13% of Eddington, within the range of previous observations. However, the best fit spectral model, adopting the double Comptonization used previously, exhibits significantly softer components. We detect seven type-I X-ray bursts, all significantly weaker (and with shorter rise and decay times) than observed previously. The burst profiles and recurrence times vary significantly, ruling out the regular bursts that are typical for this source. One burst exhibited photospheric radius expansion, and we estimate the source distance at about (5.7 / xi_b^1/2) kpc, where xi_b parameterizes the possible anisotropy of the burst emission. Interpreting the soft state as a transition from an optically thin inner flow to an optically thick flow passing through a boundary layer, as is commonly observed in similar systems, is contradicted by the lower optical depth measured for the double Comptonization model we find for this soft state. The effect of a change in disk geometry on the burst behavior remains unclear.
We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of LCDM numerical simulations.
Metal-poor stars in the Milky Way are local relics of the epoch of the first stars and the first galaxies. However, a low metallicity does not prove that a star formed in this ancient era, as metal-poor stars form over a range of redshift in different environments. Theoretical models of Milky Way formation have shown that at constant metallicity, the oldest stars are those closest to the center of the Galaxy on the most tightly-bound orbits. For that reason, the most metal-poor stars in the bulge of the Milky Way provide excellent tracers of the chemistry of the high-redshift universe. We report the dynamics and detailed chemical abundances of three stars in the bulge with [Fe/H] $\lesssim-2.7$, two of which are the most metal-poor stars in the bulge in the literature. We find that with the exception of scandium, all three stars follow the abundance trends identified previously for metal-poor halo stars. These three stars have the lowest [Sc II/Fe] abundances yet seen in $\alpha$-enhanced giant stars in the Galaxy. Moreover, all three stars are outliers in the otherwise tight [Sc II/Fe]-[Ti II/Fe] relation observed among metal-poor halo stars. Theoretical models predict that there is a 30% chance that at least one of these stars formed at $z\gtrsim15$, while there is a 70% chance that at least one formed at $10 \lesssim z \lesssim 15$. These observations imply that by $z\sim10$, the progenitor galaxies of the Milky Way had both reached [Fe/H] $\sim-3.0$ and established the abundance pattern observed in extremely metal-poor stars.
We present J-band spectroscopy of passive galaxies focusing on the Na I doublet at 1.14 {\mu}m. Like the Na I 0.82 {\mu}m doublet, this feature is strong in low-mass stars and hence may provide a useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, we find that Na I 1.14 {\mu}m increases steeply with increasing velocity dispersion, {\sigma}, and for the most massive galaxies (\sigma > 300 km/s) is much stronger than predicted from synthetic spectra with Milky-Way-like IMFs and solar abundances. Reproducing Na I 1.14 {\mu}m at high {\sigma} likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of both. Using the Na D line to break the degeneracy between IMF and abundance, we infer [Na/H] $\approx$ +0.5 and a steep IMF (single-slope-equivalent x $\approx$ 3.2, where x = 2.35 for Salpeter), for the high-\sigma galaxies. At lower mass ({\sigma} = 50-100 km/s), the line strengths are compatible with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where strong gravitational lensing masses favour MW-like IMFs. Like the high-{\sigma} sample on average, these galaxies have strong Na I 1.14 \mu m; taken in isolation their sodium indices imply bottom-heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum-fitting approach, applied to the high-\sigma sample, recovers an IMF less heavy than Salpeter, but under-predicts the Na I 1.14 \mu m line at the 5{\sigma} level. We conclude that current models struggle to reproduce this feature in the most massive galaxies without breaking other constraints, and caution against over-reliance on the sodium lines in spectroscopic IMF studies.
We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.
We present the largest number of Milky Way sized dark matter halos simulated at very high mass ($\sim$$10^4$ M$_\odot$/particle) and temporal resolution ($\sim$5 Myrs/snapshot) done to date, quadrupling what is currently available in the literature. This initial suite consists of the first 24 halos of the $Caterpillar$ $Project$ (www.caterpillarproject.org) whose project goal of 60 - 70 halos will be made public when complete. We resolve $\sim$20,000 gravitationally bound subhalos within the virial radius of each host halo. Over the ranges set by our spatial resolution our convergence is excellent and improvements were made upon current state-of-the-art halo finders to better identify substructure at such high resolutions (e.g., on average we recover $\sim$4 subhalos in each host halo above 10$^8$ M$_\odot$ which would have otherwise not been found using conventional methods). For our relaxed halos, the inner profiles are reasonably fit by Einasto profiles ($\alpha$ = 0.169 $\pm$ 0.023) though this depends on the relaxed nature and assembly history of a given halo. Averaging over all halos, the substructure mass fraction is $f_{m,subs} = 0.121 \pm 0.041$, and mass function slope is d$N$/d$M\propto M^{-1.88 \pm 0.10}$ though we find scatter in the normalizations for fixed halo mass due to more concentrated hosts having less subhalos at fixed subhalo mass. There are no biases stemming from Lagrangian volume selection as all Lagrangian volume types are included in our sample. Our detailed contamination study of 264 low resolution halos has resulted in obtaining very large and unprecedented, high-resolution regions around our host halos for our target resolution (sphere of radius $\sim$$1.4 \pm 0.4$ Mpc) allowing for accurate studies of low mass dwarf galaxies at large galactocentric radii and the very first stellar systems at high redshift ($z \geq$ 10).
We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.
We use a sample of 151 local non-blazar AGN selected from the INTEGRAL all-sky hard X-ray survey to investigate if the observed declining trend of the fraction of obscured (i.e. showing X-ray absorption) AGN with increasing luminosity is mostly an intrinsic or selection effect. Using a torus-obscuration model, we demonstrate that in addition to negative bias, due to absorption in the torus, in finding obscured AGN in hard X-ray flux limited surveys, there is also positive bias in finding unobscured AGN, due to Compton reflection in the torus. These biases can be even stronger taking into account plausible intrinsic collimation of hard X-ray emission along the axis of the obscuring torus. Given the AGN luminosity function, which steepens at high luminosities, these observational biases lead to a decreasing observed fraction of obscured AGN with increasing luminosity even if this fraction has no intrinsic luminosity dependence. We find that if the central hard X-ray source in AGN is isotropic, the intrinsic (i.e. corrected for biases) obscured AGN fraction still shows a declining trend with luminosity, although the intrinsic obscured fraction is significantly larger than the observed one: the actual fraction is larger than $\sim 85$% at $L\lesssim 10^{42.5}$ erg/s (17--60 keV), and decreases to $\lesssim 60$% at $L\gtrsim 10^{44}$ erg/s. In terms of the half-opening angle, $\theta$, of an obscuring torus, this implies that $\theta\lesssim 30$ deg in lower-luminosity AGN, and $\theta\gtrsim 45$ deg in higher-luminosity ones. If, however, the emission from the central SMBH is collimated as $dL/d\Omega\propto\cos\alpha$, the intrinsic dependence of the obscured AGN fraction is consistent with a luminosity-independent torus half-opening angle $\theta\sim 30$ deg.
This is the fourth paper in a series that reports on our investigation of the clustering properties of active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the clustering of broad-line, luminous AGN at 0.16<z<0.36. We fit the H-alpha line profile in the SDSS spectra for all X-ray and optically-selected broad-line AGN, determine the mass of the super-massive black hole (SMBH), M_BH, and infer the accretion rate relative to Eddington (L/L_EDD). Since M_BH and L/L_EDD are correlated, we create AGN subsamples in one parameter while maintaining the same distribution in the other parameter. In both the X-ray and optically-selected AGN samples we detect a weak clustering dependence with M_BH and no statistically significant dependence on L/L_EDD. We find a difference of up to 2.7sigma when comparing the objects that belong to the 30% least and 30% most massive M_BH subsamples, in that luminous broad-line AGN with more massive black holes reside in more massive parent dark matter halos at these redshifts. These results provide evidence that higher accretion rates in AGN do not necessarily require dense galaxy environments in which more galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We also present semi-analytic models which predict a positive M_DMH dependence on M_BH, which is most prominent at M_BH ~ 10^{8-9} M_SUN.
This paper is the third in a series presenting the results of direct numerical integrations of the Fokker-Planck equation for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm of Paper II included diffusion coefficients that described the effects of random ("classical") and correlated ("resonant") relaxation. In this paper, the diffusion coefficients of Paper II have been generalized to account for the effects of "anomalous relaxation," the qualitatively different way in which eccentric orbits evolve in the regime of rapid relativistic precession. Two functional forms for the anomalous diffusion coefficients are investigated, based on power-law or exponential modifications of the resonant diffusion coefficients. The parameters defining the modified coefficients are first constrained by comparing the results of Fokker-Planck integrations with previously-published N-body integrations. Steady-state solutions are then obtained via the Fokker-Planck equation for models with properties similar to those of the Milky Way nucleus. Inclusion of anomalous relaxation leads to the formation of less prominent cores than in the case of resonant relaxation alone, due to the lengthening of diffusion timescales for eccentric orbits. Steady-state capture rates of stars by the SBH are found to always be less, by as much as an order of magnitude, than capture rates in the presence of resonant relaxation alone.
The Stripe 82 Massive Galaxy Catalog (S82-MGC) is the largest-volume stellar mass-limited sample of galaxies beyond z~1 constructed to date. Spanning 139.4 deg2, the S82-MGC includes a mass-limited sample of 41,770 galaxies with log Mstar > 11.2 to z~0.7, sampling a volume of 0.3 Gpc3, roughly equivalent to the volume of the Sloan Digital Sky Survey-I/II (SDSS-I/II) z < 0.15 MAIN sample. The catalog is built on three pillars of survey data: the SDSS Stripe 82 Coadd photometry which reaches r-band magnitudes of 23.5 AB, YJHK photometry at depths of 20th magnitude (AB) from the UK Infrared Deep Sky Survey (UKIDSS) Large Area Survey, and over 70,000 spectroscopic galaxy redshifts from SDSS-I/II and the Baryon Oscillation Spectroscopic Survey (BOSS). We describe the catalog construction and verification, the production of 9-band matched aperture photometry, tests of existing and newly estimated photometric redshifts required to supplement spectroscopic redshifts for 55% of the log Mstar > 11.2 sample, and geometric masking. We provide near-IR based stellar mass estimates and compare these to previous estimates. All catalog products are made publicly available. The S82-MGC not only addresses previous statistical limitations in high-mass galaxy evolution studies but begins tackling inherent data challenges in the coming era of wide-field imaging surveys.
We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how serious errors of interpretation can only be revealed through better sampling. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that strong feedback plays a key role in flattening metal gradients in early star-forming galaxies.
Previously we showed that a substantially misaligned viscous accretion disk with pressure that orbits around one component of a binary system can undergo global damped Kozai-Lidov (KL) oscillations. These oscillations produce periodic exchanges of the disk eccentricity with inclination. The disk KL mechanism is quite robust and operates over a wide range of binary and disk parameters. However, the effects of self-gravity, which are expected to suppress the KL oscillations for sufficiently massive disks, were ignored. Here, we analyze the effects of disk self-gravity by means of hydrodynamic simulations and compare the results with the expectations of analytic theory. The disk mass required for suppression in the simulations is a few percent of the mass of the central star and this roughly agrees with an analytical estimate. The conditions for suppression of the KL oscillations in the simulations are close to requiring that the disk be gravitationally unstable. We discuss some implications of our results for the dynamics of protoplanetary disks and the related planet formation.
Very high energy (VHE $>$100 GeV) gamma rays coming from blazars can produce pairs when interacting with the Extragalactic Background Light (EBL) and the Cosmic Microwave Background, generating an electromagnetic cascade. Depending on the Intergalactic Magnetic Field (IGMF) intensity, this cascade may result in an extended isotropic emission of photons around the source (halo), or in a broadening of the emission beam. The detection of these effects might lead to important constrains both on the IGMF intensity and the EBL density, quantities of great relevance in cosmological models. Using a Monte Carlo program, we simulate electromagnetic cascades for different values of the IGMF intensities and coming from a source similar to 1ES0229+200, a blazar with hard intrinsic spectrum at redshift $z=0.14$, which is an ideal distance for potentially observing the effect. We study the possible response of a generic future Cherenkov telescope using a simplified model for the sensitivity, effective area and angular resolution. Finally, combining these instrument properties, we calculate the angular distribution of photons and develop a method to test the statistical feasibility of detecting the effect in the near future.
In imaging atmospheric Cherenkov telescope (IACT) arrays, the standard method of statistically inferring the existence of a source is based on the maximum likelihood method of Li&Ma (1983). We present a new statistical approach, also based on maximum likelihood theory, which takes into account a priori knowledge of the source light curve. This approach is especially useful for observations of rapidly decaying gamma-ray bursts (GRBs). We also discuss results established by using this technique to analyze VERITAS GRB observations.
It is anticipated that the forthcoming Cherenkov Telescope Array (CTA) will include a number of medium-sized telescopes that are constructed using a dual-mirror Schwarzschild-Couder configuration. These telescopes will sample a wide ($8^{\circ}$) field of view using a densely pixelated camera comprising over $10^{4}$ individual readout channels. A readout frequency congruent with the expected single-telescope trigger rates would result in substantial data rates. To ameliorate these data rates, a novel, hardware-level Distributed Intelligent Array Trigger (DIAT) is envisioned. A copy of the DIAT operates autonomously at each telescope and uses reduced metadata from a limited subset of nearby telescopes to veto events prior to camera readout. We present the results of Monte-Carlo simulations that evaluate the efficacy of a "Parallax width" discriminator that can be used by the DIAT to efficiently distinguish between genuine gamma-ray initiated events and unwanted background events that are initiated by hadronic cosmic rays.
The Large Binocular Telescope Interferometer (LBTI) is a strategic instrument of the LBT designed for high-sensitivity, high-contrast, and high-resolution infrared (1.5-13 $\mu$m) imaging of nearby planetary systems. To carry out a wide range of high-spatial resolution observations, it can combine the two AO-corrected 8.4-m apertures of the LBT in various ways including direct (non-interferometric) imaging, coronagraphy (APP and AGPM), Fizeau imaging, non-redundant aperture masking, and nulling interferometry. It also has broadband, narrowband, and spectrally dispersed capabilities. In this paper, we review the performance of these modes in terms of exoplanet science capabilities and describe recent instrumental milestones such as first-light Fizeau images (with the angular resolution of an equivalent 22.8-m telescope) and deep interferometric nulling observations.
This is the third in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Our sample comprises 40 clusters identified as being dynamically relaxed and hot (i.e., massive) in Papers I and II of this series. Here we consider the thermodynamics of the intracluster medium, in particular the profiles of density, temperature and related quantities, as well as integrated measurements of gas mass, average temperature, total luminosity and center-excluded luminosity. We fit power-law scaling relations of each of these quantities as a function of redshift and cluster mass, which can be measured precisely and with minimal bias for these relaxed clusters. For the thermodynamic profiles, we jointly model the density and temperature and their intrinsic scatter as a function of radius, thus also capturing the behavior of the gas pressure and entropy. For the integrated quantities, we also jointly fit a multidimensional intrinsic covariance, providing the first observational constraints on many of the off-diagonal terms of the covariance matrix. Our results reinforce the view that self similar theory provides a good description of relaxed clusters outside their centers (radii $r > 0.5\,r_{2500} \approx 0.15\,r_{500}$), but that some combination of heating and cooling processes breaks self similarity within those centers. As cosmological tests using clusters continue to improve, statistical models with the power and flexibility to reflect the internal astrophysics of cluster formation and evolution will be needed; the analysis presented here is a step in that direction.
We present a new image of the 5.5 GHz radio emission from the extended Chandra Deep Field South. Deep radio observations at 5.5 GHz were obtained in 2010 and presented in the first data release. A further 76 hours of integration has since been obtained, nearly doubling the integration time. This paper presents a new analysis of all the data. The new image reaches 8.6 microJy rms, an improvement of about 40% in sensitivity. We present a new catalogue of 5.5 GHz sources, identifying 212 source components, roughly 50% more than were detected in the first data release. Source counts derived from this sample are consistent with those reported in the literature for S_{5.5GHz} > 0.1 mJy but significantly lower than published values in the lowest flux density bins (S_{5.5GHz} < 0.1 mJy), where we have more detected sources and improved statistical reliability. The 5.5 GHz radio sources were matched to 1.4 GHz sources in the literature and we find a mean spectral index of -0.35 +- 0.10 for S_{5.5GHz} > 0.5 mJy, consistent with the flattening of the spectral index observed in 5 GHz sub-mJy samples. The median spectral index of the whole sample is \alpha_{med} = -0.58, indicating that these observations may be starting to probe the star forming population. However, even at the faintest levels (0.05 < S_{5.5GHz} < 0.1 mJy), 39% of the 5.5 GHz sources have flat or inverted radio spectra. Four flux density measurements from our data, across the full 4.5 to 6.5 GHz bandwidth, are combined with those from literature and we find 10% of sources (S_{5.5GHz} >~ 0.1 mJy) show significant curvature in their radio spectral energy distribution spanning 1.4 to 9 GHz.
The fate of surface water on Venus is one of the most important outstanding problems in comparative planetology. Here a new concept is proposed to explain water removal on a steam-covered proto Venus, referred to as impact-driven planetary desiccation. Since a steam atmosphere is photochemically unstable, water vapor dissociates into hydrogen and oxygen. Then, hydrogen escapes easily into space through hydrodynamic escape driven by strong extreme ultraviolet radiation from the young Sun. The focus is on the intense impact bombardment during the terminal stage of planetary accretion as generators of a significant amount of reducing agent. The fine-grained ejecta remove the residual oxygen, the counter part of escaped hydrogen, via the oxidation of iron-bearing rocks in a hot atmosphere. Thus, hypervelocity impacts cause net desiccation of the planetary surface. I constructed a stochastic cratering model using a Monte Carlo approach to investigate the cumulative mass of nonoxidized, ejected rocks due to the intense impact bombardment. It is shown that a thick steam atmosphere with a mass equivalent to that of the terrestrial oceans would be removed. The cumulative mass of rocky ejecta released into the atmosphere reaches 1 wt% of the host planet, which is 10000 times of the current mass of the Earths atmosphere. These results strongly suggest that chemical reactions between such large amounts of ejecta and planetary atmospheres are among the key factors required to understand atmospheric mass and its composition, not only in the Solar System but also in extrasolar systems.
The Cherenkov Telescope Array (CTA), is an international project for the next generation ground- based observatory for gamma-ray astronomy in the energy range from 20 GeV to 300 TeV. The sensitivity in the core energy range will be dominated by up to 40 Medium-Sized Telescopes (MSTs). The MSTs, of Davies-Cotton type with a 12 m diameter reflector are currently in the prototype phase. A full-size mechanical telescope structure has been assembled in Berlin. The telescope is partially equipped with different mirror prototypes, which are currently being tested and evaluated for performances characteristics. A report concentrating on the details of the tele- scope structure, the drive assemblies and the optics of the MST prototype will be given.
In this work we propose the Amor-type asteroid 2009 WN25 as the likely progenitor of the November i-Draconids (NID, IAU#392), a recently detected weak annual meteoroid stream. We first describe our recovery and follow-up effort to obtain timely ground based astrometry with large aperture telescopes, and ensure that 2009 WN25 would not become lost. We then discuss the possible parent-stream association, using its updated orbit to model the ejection of dust particles from the surface of the parent body and match the observed properties of the stream.
The camera of the Large Size Telescopes (LSTs) of the Cherenkov Telescope Array (CTA) consists of 265 photosensor modules, each of them containing 7 photomultiplier tubes (PMTs), a slow control board (SCB), a readout board, and a trigger logic. We have developed the SCB, which is installed between the 7 PMTs and the readout board. The main task for SCBs is the controlling of the high voltages for the PMTs and the monitoring of their anode currents. In addition, the SCB provides the functionality to create test pulses that can be injected at the input of the PMT preamplifier in order to emulate a PMT signal without the need of setting a high voltage, or even without the PMT itself. The test pulses have a very similar width as the PMT pulses (less than 3 ns FWHM) and their amplitude can be adjusted in a wide dynamic range. These features allow us not only to test the functionality of the camera modules but also to fully characterize these. We report on the design and the functions of the SCB together with the results of test measurements.
Cosmic-ray (CR) protons can accumulate for cosmological times in clusters of galaxies. Their hadronic interactions with protons of the intra-cluster medium (ICM) generate secondary electrons, gamma-rays and high-energy neutrinos. In light of the high-energy neutrino events recently discovered by the IceCube observatory, we estimate the contribution from galaxy clusters to the diffuse gamma-ray and neutrino backgrounds. For the first time, we consistently take into account the synchrotron emission generated by secondary electrons and require the clusters radio counts to be respected. For a choice of parameters respecting current constraints from radio to gamma-rays, and assuming a proton spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the IceCube flux, and much less to the total extragalactic gamma-ray background observed by Fermi. They account for less than 1% for spectral indexes <-2. The high-energy neutrino flux observed by IceCube can be reproduced without violating radio constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the high-energy IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with the particle energy. We stress that our results are valid for all kind of sources injecting CR protons into the ICM, and that, while IceCube can test the most optimistic scenarios for spectral indexes >=-2.2 by stacking few nearby massive objects, clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds in any realistic scenario.
The Cherenkov Telescope Array (CTA) observatory for very high-energy gamma rays will consist of about a hundred of imaging atmospheric Cherenkov telescopes (IACTs) of different size with a total reflective area of about 10,000 m$^2$. Here we present a novel technology for the production of IACT mirrors that has been developed in the Institute of Nuclear Physics PAS in Krakow, Poland. The mirrors are made by cold-slumping of the front reflecting aluminium-coated panel and the rear panel interspaced with aluminium spacers. Each panel is built of two glass panels laminated with a layer of a fibreglass tissue in between for reinforcement of the structure against mechanical damage. The mirror structure is open and does not require a perfect sealing needed in closed-type designs. It prohibits water to be trapped inside and enables a proper ventilation of the mirror. Full-size hexagonal prototype mirrors produced for the medium-sized CTA telescopes will be presented together with the results of recent comprehensive optical and durability tests. Their design will be compared to the earlier technology developed at INP PAS that used a rigid flat open support structure with a reflective layer made by cold-slumping of the coated glass panel to the cast-in-mould spherical epoxy resin layer.
The number of galaxies at a given flux as a function of the redshift, $z$, is derived when the $z$-distance relation is non-standard. In order to compare different models, the same formalism is also applied to the standard cosmology. The observed luminosity function for galaxies of the zCOSMOS catalog at different redshifts is modelled by a new luminosity function for galaxies, which is derived by the truncated beta probability density function. Three astronomical tests, which are the photometric maximum as a function of the redshift for a fixed flux, the mean value of the redshift for a fixed flux, and the luminosity function for galaxies as a function of the redshift, compare the theoretical values of the standard and non-standard model with the observed value. The tests are performed on the FORS Deep Field (FDF) catalog up to redshift $z=1.5$ and on the zCOSMOS catalog extending beyond $z=4$. These three tests show minimal differences between the standard and the non-standard models.
We used ultra-deep UV observations obtained with the Hubble Space Telescope to search for optical companions to binary millisecond pulsars (MSPs) in the globular cluster 47 Tucanae. We identified four new counterparts (to MSPs 47TucQ, 47TucS, 47TucT and 47TucY) and confirmed those already known (to MSPs 47TucU and 47TucW). In the color magnitude diagram, the detected companions are located in a region between the main sequence and the CO white dwarf cooling sequences, consistent with the cooling tracks of He white dwarfs of mass between 0.15 Msun and 0.20 Msun. For each identified companion, mass, cooling age, temperature and pulsar mass (as a function of the inclination angle) have been derived and discussed. For 47TucU we also found that the past accretion history likely proceeded in a sub-Eddington rate. The companion to the redback 47TucW is confirmed to be a non degenerate star, with properties particularly similar to those observed for black widow systems. Two stars have been identified within the 2-sigma astrometric uncertainty from the radio positions of 47TucH and 47TucI, but the available data prevent us from firmly assessing whether they are the true companions of these two MSPs.
Penumbral microjets are short-lived, fine-structured and bright jets that are generally observed in chromospheric imaging of the penumbra of sunspots. Here we investigate their potential transition region signature, by combining observations with the Swedish 1-m Solar Telescope (SST) in the Ca II H and Ca II 8542{\AA} lines with ultraviolet imaging and spectroscopy obtained with the Interface Region Imaging Spectrograph (IRIS), which includes the C II 1334/1335{\AA}, Si IV 1394/1403{\AA} and Mg II h & k 2803/2796{\AA} lines. We find a clear corresponding signal in the IRIS Mg II k, C II and Si IV slit-jaw images, typically offset spatially from the Ca II signature in the direction along the jets: from base to top, the penumbral microjets are predominantly visible in Ca II, Mg II k and C II/Si IV, suggesting progressive heating to transition region temperatures along the jet extent. Hence, these results support the suggestion from earlier studies that penumbral microjets may heat to transition region temperatures.
We compared optical spectroscopic and photometric data for 18 AGN galaxies over 2 to 3 epochs, with time intervals of typically 5 to 10 years. We used the Multi-Object Double Spectrograph (MODS) at the Large Binocular Telescope (LBT) and compared the spectra to data taken from the SDSS database and the literature. We find variations in the forbidden oxygen lines as well as in the hydrogen recombination lines of these sources. For 4 of the sources we find that, within the calibration uncertainties, the variations in continuum and line spectra of the sources are very small. We argue that it is mainly the difference in black hole mass between the samples that is responsible for the different degree of continuum variability. In addition we find that for an otherwise constant accretion rate the total line variability (dominated by the narrow line contributions) reverberates the continuum variability with a dependency $\Delta L_{line} \propto (\Delta L_{cont.})^{\frac{3}{2}}$. Since this dependency is prominently expressed in the narrow line emission it implies that the luminosity dominating part of the narrow line region must be very compact with a size of the order of at least 10 light years. A comparison to literature data shows that these findings describe the variability characteristics of a total of 61 broad and narrow line sources.
For future space infrared astronomical coronagraphy, we perform experimental studies on the application of aluminum mirrors to a coronagraph. Cooled reflective optics is required for broad-band mid-infrared observations in space, while high-precision optics is required for coronagraphy. For the coronagraph instrument originally proposed for the next-generation infrared astronomical satellite project SPICA (SCI: SPICA Coronagraph Instrument), we fabricated and evaluated the optics consisting of high-precision aluminum off-axis mirrors with diamond-turned surfaces, and conducted a coronagraphic demonstration experiment using the optics with a coronagraph mask. We first measured the wave front errors (WFEs) of the aluminum mirrors with a He-Ne Fizeau interferometer to confirm that the power spectral densities of the WFEs satisfy the SCI requirements. Then we integrated the mirrors into an optical system and evaluated the overall performance of the system. As a result, we estimate the total WFE of the optics to be 33 nm (rms), each mirror contributing 10-20 nm (rms) for the central 14 mm area of the optics, and obtain a contrast of 10^(-5.4) as a coronagraph in the visible light. At a wavelength of 5 um, the coronagraphic system is expected to achieve a contrast of ~10^(-7) based on our model calculation with the measured optical performance. Thus our experiment demonstrates that aluminum mirror optics is applicable to a highly WFE-sensitive instrument such as a coronagraph in space.
The purpose of this HST white paper is to demonstrate that it is possible to monitor Jupiter's polar haze with HST/STIS without breaking the ground screening limit for bright objects. This demonstration rests on a thorough simulation of STIS output from an existing image obtained with HST/WFPC2. It is shown that the STIS NUV-MAMA + F25CIII filter assembly provides a count rate per pixel ~11 times smaller than that obtained for one pixel of WFPC2 WF3 CCD + F218W corresponding filter. This ratio is sufficiently large to cope with the bright solar light scattered by Jupiter's atmosphere, which was a lesser concern for WFPC2 CCD safety. These STIS images would provide unprecedented spatial and temporal resolution observations of small-scale stratospheric aerosol structures, possibly associated with Jupiter's complex FUV aurora.
G-band bright points (GBPs) are thought to be the foot-points of magnetic flux tubes. The aim of this paper is to investigate the relation between the diffusion regimes of GBPs and the associated longitudinal magnetic field strengths. Two high resolution observations of different magnetized environments were acquired with the Hinode/Solar Optical Telescope. Each observation was recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images. GBPs are identified and tracked automatically, and then categorized into several groups by their longitudinal magnetic field strengths, which are extracted from the calibrated NFI magnetograms using a point-by-point method. The Lagrangian approach and the distribution of diffusion indices approach are adopted separately to explore the diffusion regime of GBPs for each group. It is found that the values of diffusion index and diffusion coefficient both decrease exponentially with the increasing longitudinal magnetic field strengths whichever approach is used. The empirical formulas deduced from the fitting equations are proposed to describe these relations. Stronger elements tend to diffuse more slowly than weak elements, independently of the magnetic flux of the surrounding medium. This may be because the magnetic energy of stronger elements is not negligible compared with the kinetic energy of the gas, and therefore the flows cannot perturb them so easily.Yang
Based on its energy-dependent morphology the initially unidentified very high energy (VHE; E>100 GeV) gamma-ray source HESS J1303-631 was recently associated with the pulsar PSR J1301-6305. Subsequent detection of X-ray and GeV counterparts also supports the identification of the H.E.S.S. source as evolved pulsar wind nebula (PWN). We report here on recent radio observations of the PSR J1301-6305 field of view (FOV) with ATCA dedicated to search for the radio counterpart of this evolved PWN. Observations at 5.5 GHz and 7.5 GHz do not reveal any extended emission associated with the pulsar. The analysis of the archival 1.384 GHz and 2.368 GHz data also does not show any significant emission. The 1.384 GHz data reveal a hint of an extended shell-like emission in the PSR J1301-6305 FOV which might be a supernova remnant. We discuss the implications of the non-detection at radio wavelengths on the nature and evolution of the PWN as well as the possibility of the SNR candidate being the birth place of PSR J1301-6305.
The Cherenkov Telescope Array (CTA) is a ground-based $\gamma$-ray observatory that will observe the full sky in the energy range from 20 GeV to 100 TeV from facilities in both hemispheres. It is proposed to consist of more than 100 telescopes, producing large amounts of data. Apart from the storage system, there are also requirements on the software framework to allow efficient data processing, i.e. robustness, execution speed and coding efficiency. This contribution will present a plain and simple pipeline framework design prototype for CTA that builds upon well-known tools, allowing the users to focus on physics problems without learning complicated software paradigms.
The Cherenkov Telescope Array (CTA) is a ground-based $\gamma$-ray observatory that will observe the full sky in the energy range from 20 GeV to 100 TeV from facilities in both hemispheres. It is proposed to consist of more than 100 telescopes and the large amount of data produced will exceed the volume of current VHE Imaging Atmospheric Cherenkov Telescopes by $\sim$two orders of magnitude. This volume of data represents a new challenge to the community, which is looking for new data formats to transfer and store the CTA data. One of the prototypes currently under study is the ROI (Regions Of Interest) file format for camera images. It can store only those pixels of a camera image that are close to the shower, thus removing the major part of the night sky background (NSB) while keeping all pixels that might belong to the shower. Simple on-the-fly compression is used to reduce the file size even further. Here, we explain the ROI prototype in detail and present preliminary results when applied to simulations.
PSR B1259-63/LS 2883 is a very high energy (VHE; E > 100 GeV) gamma-ray emitting binary consisting of a 48 ms pulsar orbiting around a Be star with a period of 3.4 years. The Be star features a circumstellar disk which is inclined with respect to the orbit in such a way that the pulsar crosses it twice every orbit. The circumstellar disk provides an additional field of target photons which may contribute to inverse Compton scattering and gamma-gamma absorption, leaving a characteristic imprint in the observed spectrum and light curve of the high energy emission. We study the signatures of Compton-supported, VHE gamma-ray induced pair cascades in the circumstellar disc of the Be star and their possible contribution to the GeV flux. We also study a possible impact of the gamma-gamma absorption in the disk on the observed TeV light curve. We show that the cumulative absorption of VHE gamma-rays in stellar and disk photon fields can explain the modulation of the flux at the periastron passage.
The first Gamma-Ray Burst (GRB) catalog presented by the Fermi-Large Area Telescope (LAT) collaboration includes 28 GRBs, detected above 100 MeV over the first three years since the launch of the Fermi mission. However, more than 100 GRBs are expected to be found over a period of six years of data collection thanks to a new detection algorithm and to the development of a new LAT event reconstruction, the so-called "Pass 8." Our aim is to provide revised prospects for GRB alerts in the CTA era in light of these new LAT discoveries. We focus initially on the possibility of GRB detection with the Large Size Telescopes (LSTs). Moreover, we investigate the contribution of the Middle Size Telescopes (MSTs), which are crucial for the search of larger areas on short post trigger timescales. The study of different spectral components in the prompt and afterglow phase, and the limits on the Extragalactic background light are highlighted. Different strategies to repoint part of - or the entire array - are studied in detail.
Standing sausage modes in flare loops are important for interpreting quasi-periodic pulsations (QPPs) in solar flare lightcurves. We propose an inversion scheme that consistently uses their periods $P$ and damping times $\tau$ to diagnose flare loop parameters. We derive a generic dispersion relation governing linear sausage waves in pressure-less straight tubes, for which the transverse density inhomogeneity takes place in a layer of arbitrary width $l$ and is of arbitrary form. We find that $P$ and $\tau$ depend on the combination of $[R/v_{\rm Ai}, L/R, l/R, \rho_{\rm i}/\rho_{\rm e}]$, where $R$ is the loop radius, $L$ is the looplength, $v_{\rm Ai}$ is the internal Alfv\'en speed, and $\rho_{\rm i}/\rho_{\rm e}$ is the density contrast. For all the density profiles examined, $P$ and $\tau$ experience saturation when $L/R \gg 1$, yielding an inversion curve in the $[R/v_{\rm Ai}, l/R, \rho_{\rm i}/\rho_{\rm e}]$ space with a specific density profile when $L/R$ is sufficiently large. When applied to a spatially unresolved QPP event, the scheme yields that $R/v_{\rm Ai}$ is the best constrained, whereas $l/R$ corresponds to the other extreme. For spatially resolved QPPs, while $L/R \gg 1$ cannot be assumed beforehand, an inversion curve remains possible due to additional geometrical constraints. When a spatially resolved QPP event involves another mode, as is the case for a recent event, the full set of $[v_{\rm Ai}, l, \rho_{\rm i}/\rho_{\rm e}]$ can be inferred. We conclude that the proposed scheme provides a useful tool for magneto-seismologically exploiting QPPs.
Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two- dimensional Cartesian or polar coordinates. In Cartesian coordinates, the calcu- lations of the hydrodynamics and self-gravitational forces are relatively straight- forward for attaining second order accuracy. However, in polar coordinates, a second order calculation of self-gravitational forces is required for matching the second order accuracy of hydrodynamical schemes. We present a direct algorithm for calculating self-gravitational forces with second order accuracy without artifi- cial boundary conditions. The Poisson integral in polar coordinates is expressed in a convolution form and the corresponding numerical complexity is nearly lin- ear using a fast Fourier transform. Examples with analytic solutions are used to verify that the truncated error of this algorithm is of second order. The kernel integral around the singularity is applied to modify the particle method. The use of a softening length is avoided and the accuracy of the particle method is significantly improved.
The basic formulation describing quadratic mode coupling in rotating Newtonian stars is presented, focusing on polar modes. Due to the Chandrasekhar-Friedman-Schutz mechanism, the f-mode (fundamental oscillation) is driven unstable by the emission of gravitational waves. If the star falls inside the so-called instability window, the mode's amplitude grows exponentially, until it is halted by non-linear effects. Quadratic perturbations form three-mode networks inside the star, which evolve as coupled oscillators, exchanging energy. Coupling of the unstable f-mode to other (stable) modes can lead to a parametric resonance and the subsequent saturation of its amplitude, thus suppressing the instability. The saturation point determines the amplitude of the gravitational-wave signal obtained from an individual source, as well as the evolutionary path of the latter inside the instability window.
We present wide field near-infrared photometry of 12 Galactic globular clusters, typically extending from the tip of the cluster red giant branch (RGB) to the main sequence turnoff. Using recent homogenous values of cluster distance, reddening and metallicity, the resulting photometry is directly compared to the predictions of several recent libraries of stellar evolutionary models. Of the sets of models investigated, Dartmouth and Victoria-Regina models best reproduce the observed RGB morphology, albeit with offsets in J-Ks color which vary in their significance in light of all sources of observational uncertainty. Therefore, we also present newly recalibrated relations between near-IR photometric indices describing the upper RGB versus cluster iron abundance as well as global metallicity. The influence of enhancements in alpha elements and helium are analyzed, finding that the former affect the morphology of the upper RGB in accord with model predictions. Meanwhile, the empirical relations we derive are in good agreement with previous results, and minor discrepancies can likely be attributed to differences in the assumed cluster distances and reddenings. In addition, we present measurements of the horizontal branch (HB) and RGB bump magnitudes, finding a non-negligible dependence of the near-IR HB magnitude on cluster metallicity. Lastly, we discuss the influence of assumed cluster distances, reddenings and metallicities on our results, finding that our empirical relations are generally insensitive to these factors to within their uncertainties.
The Gamma Cherenkov Telescope (GCT) is proposed to be part of the Small Size Telescope (SST) array of the Cherenkov Telescope Array (CTA). The GCT dual-mirror optical design allows the use of a compact camera of diameter roughly 0.4 m. The curved focal plane is equipped with 2048 pixels of ~0.2{\deg} angular size, resulting in a field of view of ~9{\deg}. The GCT camera is designed to record the flashes of Cherenkov light from electromagnetic cascades, which last only a few tens of nanoseconds. Modules based on custom ASICs provide the required fast electronics, facilitating sampling and digitisation as well as first level of triggering. The first GCT camera prototype is currently being commissioned in the UK. On-telescope tests are planned later this year. Here we give a detailed description of the camera prototype and present recent progress with testing and commissioning.
Damped transverse oscillations of magnetic loops are routinely observed in the solar corona. This phenomenon is interpreted as standing kink magnetohydrodynamic waves, which are damped by resonant absorption owing to plasma inhomogeneity across the magnetic field. The periods and damping times of these oscillations can be used to probe the physical conditions of the coronal medium. Some observations suggest that interaction between neighboring oscillating loops in an active region may be important and can modify the properties of the oscillations compared to those of an isolated loop. Here we theoretically investigate resonantly damped transverse oscillations of interacting non-uniform coronal loops. We provide a semi-analytic method, based on the T-matrix theory of scattering, to compute the frequencies and damping rates of collective oscillations of an arbitrary configuration of parallel cylindrical loops. The effect of resonant damping is included in the T-matrix scheme in the thin boundary approximation. Analytic and numerical results in the specific case of two interacting loops are given as an application.
With over 1800 planets discovered outside of the Solar System in the past two
decades, the field of exoplanetology has broadened our perspective on planetary
systems. Research priorities are now moving from planet detection to planet
characterization. In this context, transiting exoplanets are of special
interest due to the wealth of data made available by their orbital
configuration. Here, I introduce two methods to gain new insights into the
atmospheric and interior properties of exoplanets.
The first method aims to map an exoplanet's atmosphere based on the scanning
obtained while it is occulted by its host star. I introduce the basics of
eclipse mapping, its caveats, and a framework to mitigate their effects via
global analyses including transits, phase curves, and radial velocity
measurements. I use this method to create the first 2D map and the first cloud
map of an exoplanet for the hot-Jupiters HD189733b and Kepler-7b, respectively.
Ultimately temperature, composition, and circulation patterns could be
constrained in 3D, a significant asset for informing atmospheric models.
The second method, MassSpec, aims to determine transiting planet masses and
atmospheric properties solely from transmission spectra. Determination of an
exoplanet's mass is key to understanding its basic properties, including its
potential for supporting life. To date, mass constraints for exoplanets are
mainly based on radial velocity measurements, which are not suited for planets
with low masses, large semi-major axes, or those orbiting faint or active
stars. I demonstrate that a planet's mass has to be accounted for by
atmospheric retrieval methods to avoid biases and that JWST could determine the
mass and atmospheric properties of half a dozen Earth-sized planets in their
host's habitable zones over its lifetime, which could lead to the first
identification of a habitable exoplanet.
The Tunka observatory is located close to Lake Baikal in Siberia, Russia. Its main detector, Tunka-133, is an array of photomultipliers measuring Cherenkov light of air showers initiated by cosmic rays in the energy range of approximately $10^{16}-10^{18}\,$eV. In the last years, several extensions have been built at the Tunka site, e.g., a scintillator array named Tunka-Grande, a sophisticated air-Cherenkov-detector prototype named HiSCORE, and the radio extension Tunka-Rex. Tunka-Rex started operation in October 2012 and currently features 44 antennas distributed over an area of about $3\,$km$^2$, which measure the radio emission of the same air showers detected by Tunka-133 and Tunka-Grande. Tunka-Rex is a technological demonstrator that the radio technique can provide an economic extension of existing air-shower arrays. The main scientific goal is the cross-calibration with the air-Cherenkov measurements. By this cross-calibration, the precision for the reconstruction of the energy and mass of the primary cosmic-ray particles can be determined. Finally, Tunka-Rex can be used for cosmic-ray physics at energies close to $1\,$EeV, where the standard Tunka-133 analysis is limited by statistics. In contrast to the air-Cherenkov measurements, radio measurements are not limited to dark, clear nights and can provide an order of magnitude larger exposure.
The simple assumption of an instantaneous reionization of the Universe may bias estimates of cosmological parameters. In this paper a model-independent principal component method for the reionization history is applied to give constraints on the cosmological parameters from recent Planck 2015 data. We find that the Universe are not completely reionized at redshifts $z \ge 8.5$ at 95% CL. Both the reionization optical depth and the matter fluctuation amplitude are higher than but consistent with those obtained in the standard instantaneous reionization scheme. The high estimated value of the matter fluctuation amplitude strengthens the tension between Planck CMB observations and some astrophysical data, such as cluster counts and weak lensing. The tension can significantly be relieved if the neutrino masses are allowed to vary. Thanks to a high scalar spectral index, the low-scale spontaneously broken SUSY inflationary model can fit the data well, which is marginally disfavored at 95% CL in the Planck analysis.
The CRESST-II experiment uses cryogenic detectors to search for nuclear recoil events induced by the elastic scattering of dark matter particles in CaWO$_4$ crystals. Given the low energy threshold of our detectors in combination with light target nuclei, low mass dark matter particles can be probed with high sensitivity. In this letter we present the results from data of a single detector module corresponding to 52 kg live days. A blind analysis is carried out. With an energy threshold for nuclear recoils of 307 eV we substantially enhance the sensitivity for light dark matter. Thereby, we extend the reach of direct dark matter experiments to the sub-region and demonstrate that the energy threshold is the key parameter in the search for low mass dark matter particles.
We study the distribution of the photometric rotation period (Prot), which is a direct measurement of the surface rotation at active latitudes, for three subsamples of Sun-like stars: one from CoRoT data and two from Kepler data. We identify the main populations of these samples and interpret their main biases particularly for a comparison with the solar Prot. Prot and variability amplitude (A) measurements were obtained from public CoRoT and Kepler catalogs, which were combined with public data of physical parameters. Because these samples are subject to selection effects, we computed synthetic samples with simulated biases to compare with observations, particularly around the Sun's HR-diagram location. Theoretical grids and empirical relations were used to combine physical parameters with Prot and A. Biases were simulated by performing cutoffs on the physical and rotational parameters in the same way as in each observed sample. A crucial cutoff is related with the detectability of the rotational modulation, which strongly depends on A. The synthetic samples explain the observed Prot distributions of Sun-like stars as having two main populations: one of young objects (group I, with ages below ~1 Gyr) and another of main-sequence and evolved stars (group II, with ages above ~1 Gyr). The proportions of groups I and II in relation to the total number of stars range within 64-84% and 16-36%, respectively. Hence, young objects abound in the distributions, producing the effect of observing a high number of short periods around the Sun's HR-diagram location. Differences in the Prot distributions between the CoRoT and Kepler Sun-like samples may be associated with different Galactic populations. Overall, the synthetic distribution around the solar period agrees with observations, which suggests that the solar rotation is normal with respect to Sun-like stars within the current data accuracy.
A dense neutrino medium such as that inside a core-collapse supernova can experience collective flavor conversion or oscillations because of the neutral-current weak interaction among the neutrinos. This phenomenon has been studied in a restricted, stationary supernova model which possesses the (spatial) spherical symmetry about the center of the supernova and the (directional) axial symmetry around the radial direction. Recently it has been shown that these spatial and directional symmetries can be broken spontaneously by collective neutrino oscillations. In this paper we analyze the neutrino flavor instabilities in a time-dependent supernova model. Our results show that collective neutrino oscillations start at approximately the same radius in both the stationary and time-dependent supernova models unless there exist very rapid variations in local physical conditions on timescales of a few microseconds or shorter. Our results also suggest that collective neutrino oscillations can vary rapidly with time in the regimes where they do occur which need to be studied in time-dependent supernova models.
Double-lobe radio galaxies in the local Universe have traditionally been found to be hosted in elliptical or lenticular galaxies. We report the discovery of four spiral-host double-lobe radio galaxies (J0836+0532, J1159+5820, J1352+3126 and J1649+2635) that are discovered by cross-matching a large sample of 187005 spiral galaxies from SDSS DR7 to the full catalogues of FIRST and NVSS. J0836+0532 is reported for the first time. The host galaxies are forming stars at an average rate of 1.7 $-$ 10 M$_{\odot}$ yr$^{-1}$ and possess Super Massive Black Holes (SMBHs) with masses of a few times 10$^{8}$ M$_{\odot}$. Their radio morphologies are similar to FR-II radio galaxies with total projected linear sizes ranging from 86 kpc to 420 kpc, but their total 1.4 GHz radio luminosities are only in the range 10$^{24}$ $-$ 10$^{25}$ W Hz$^{-1}$. We propose that the formation of spiral-host double-lobe radio galaxies can be attributed to more than one factor, such as the occurrence of strong interactions, mergers, and the presence of unusually massive SMBHs, such that the spiral structures are not destroyed. Only one of our sources (J1649+2635) is found in a cluster environment, indicating that processes other than accretion through cooling flows e.g., galaxy-galaxy mergers or interactions could be plausible scenarios for triggering radio-loud AGN activity in spiral galaxies.
We present a catalog of near-infrared (NIR) spectra and associated measurements for 886 nearby M dwarfs. The spectra were obtained with the NASA-Infrared Telescope Facility SpeX Spectrograph during a two-year observing campaign; they have high signal-to-noise ratios (SNR $>100-150$), span 0.8-2.4 $\mu$m and have $R\sim2000$. Our catalog of measured values contains useful T$_{\mathrm{eff}}$ and composition-sensitive features, empirical stellar parameter measurements, and kinematic, photometric, and astrometric properties compiled from the literature. We focus on measures of M dwarf abundances ([Fe/H] and [M/H]), capitalizing on the precision of recently published empirical NIR spectroscopic calibrations. We explore systematic differences between different abundance calibrations, and to other similar M dwarf catalogs. We confirm that the M dwarf abundances we measure show the expected inverse dependence with kinematic, activity, and color-based age indicators. Finally, we provide updated [Fe/H] and [M/H] for 16 M dwarf planet hosts. This catalog represents the largest published compilation of NIR spectra and associated parameters for M dwarfs. It provides a rich and uniform resource for the nearby M dwarfs, and will be especially valuable for measuring Habitable Zone locations and comparative abundances of the M dwarf planet hosts that will be uncovered by upcoming exoplanet surveys.
We present the graviton propagator for an infinite derivative, $D$-dimensional, non-local action, up to quadratic order in curvature around a Minkowski background, and discuss the conditions required for this class of gravity theory to be ghost-free. We then study the gravitational entropy for de-Sitter and Anti-de Sitter backgrounds, before comparing with a recently derived result for a Schwarzschild blackhole, generalised to arbitrary $D$-dimensions, whereby the entropy is given simply by the area law. A novel approach of decomposing the entropy into its $(r, t)$ and spherical components is adopted in order to illustrate the differences more clearly. We conclude with a discussion of de-Sitter entropy in the framework of a non-singular bouncing cosmology.
Theories with elementary scalar degrees of freedom seem nowadays required for simple descriptions of the Standard Model and of the Early Universe. It is then natural to embed theories of inflation in supergravity, also in view of their possible ultraviolet completion in String Theory. After some general remarks on inflation in supergravity, we describe examples of minimal inflaton dynamics which are compatible with recent observations, including higher-curvature ones inspired by the Starobinsky model. We also discuss different scenarios for supersymmetry breaking during and after inflation, which include a revived role for non-linear realizations. In this spirit, we conclude with a discussion of the link, in four dimensions, between "brane supersymmetry breaking" and the super--Higgs effect in supergravity.
We develop a framework that facilitates the study of the causal structure of spacetimes with a causally preferred foliation. Such spacetimes may arise as solutions of Lorentz-violating theories, e.g. Horava gravity. Our framework allows us to rigorously define concepts such as black/white holes and to formalize the notion of a `universal horizon', that has been previously introduced in the simpler setting of static and spherically symmetric geometries. We also touch upon the issue of development and prove that universal horizons are Cauchy horizons when evolution depends on boundary data or asymptotic conditions. We establish a local characterisation of universal horizons in stationary configurations. Finally, under the additional assumption of axisymmetry, we examine under which conditions these horizons are cloaked by Killing horizons, which can act like usual event horizons for low-energy excitations.
We perform a theoretical and numerical study of anti-parallel 2D magnetic reconnection with asymmetries in the density and reconnecting magnetic field strength in addition to a bulk flow shear across the reconnection site in the plane of the reconnecting fields, which commonly occurs at planetary magnetospheres. We predict the speed at which an isolated X-line is convected by the flow, the reconnection rate, and the critical flow speed at which reconnection no longer takes place for arbitrary reconnecting magnetic field strengths, densities, and upstream flow speeds, and confirm the results with two-fluid numerical simulations. The predictions and simulation results counter the prevailing model of reconnection at Earth's dayside magnetopause which says reconnection occurs with a stationary X-line for sub-Alfvenic magnetosheath flow, reconnection occurs but the X-line convects for magnetosheath flows between the Alfven speed and double the Alfven speed, and reconnection does not occur for magnetosheath flows greater than double the Alfven speed. We find that X-line motion is governed by momentum conservation from the upstream flows, which are weighted differently in asymmetric systems, so the X-line convects for generic conditions including sub-Alfvenic upstream speeds. For the reconnection rate, while the cutoff condition for symmetric reconnection is that the difference in flows on the two sides of the reconnection site is twice the Alfven speed, we find asymmetries cause the cutoff speed for asymmetric reconnection to be higher than twice the asymmetric form of the Alfven speed. The results compare favorably with an observation of reconnection at Earth's polar cusps during a period of northward interplanetary magnetic field, where reconnection occurs despite the magnetosheath flow speed being more than twice the magnetosheath Alfven speed, the previously proposed suppression condition.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.
Coherent radiation at radio frequencies from high-energy showers fully contained in a dense radio-transparent medium - like ice, salt or regolith - has been extensively investigated as a promising technique to search for ultra-high energy (UHE) neutrinos. Additional emission in the form of transition radiation may occur when a neutrino-induced shower produced close to the Earth surface emerges from the ground into atmospheric air. We present the first detailed evaluation of transition radiation from high-energy showers crossing the boundary between two different media. We found that transition radiation is sizable over a wide solid angle and coherent up to $\sim$ 1 GHz. These properties encourage further work to evaluate the potential of a large-aperture UHE neutrino experiment based on detection of transition radiation.
The progenitors of SNe Ia are not yet fully understood. The two leading progenitor scenarios are the single-degenerate (SD) scenario and the double-degenerate scenario. In the SD scenario, the collision of the SN Ia ejecta with its companion star is expected to produce detectable ultraviolet (UV) emission in the first few days after the SN explosion within certain viewing angles. A strong UV flash has recently been detected in an SN 2002es-like peculiar SN Ia iPTF14atg by Cao et al., which is interpreted as evidence of an early-time UV signature due to SN ejecta interacting with its companion star, supporting the SD scenario. In this paper, we present the expected luminosity distributions of early-time UV emission arising from SN Ia ejecta-companion interaction by performing binary population synthesis calculations for different progenitor systems in the SD scenario. Our theoretical predictions will be helpful for future early-time observations of SNe Ia to constrain their possible progenitors. Assuming the observed strong UV pulse of iPTF14atg was indeed produced by the SN ejecta-companion interaction, our population synthesis model suggests that the progenitor system of iPTF14atg is most likely a red-giant donor binary system, and it is unlikely to have been a main-sequence or helium-star donor system.
Hubble Space Telescope observations of the site of the supernova (SN) 2008ax obtained in 2011 and 2013 reveal that the possible progenitor object detected in pre-explosion images was in fact multiple. Four point sources are resolved in the new, higher-resolution images. We identify one of the sources with the fading SN. The other three objects are consistent with single supergiant stars. We conclude that their light contaminated the previously identified progenitor candidate. After subtraction of these stars, the progenitor appears to be significantly fainter and bluer than previously measured. Post-explosion photometry at the SN location indicates that the progenitor object has disappeared. If single, the progenitor is compatible with a supergiant star of B to mid-A spectral type, while a Wolf-Rayet (WR) star would be too luminous in the ultraviolet to account for the observations. Moreover, our hydrodynamical modelling shows the pre-explosion mass was $4-5$ $M_\odot$ and the radius was $30-50$ $R_\odot$, which is incompatible with a WR progenitor. We present a possible interacting binary progenitor computed with our evolutionary models that reproduces all the observational evidence. A companion star as luminous as an O9-B0 main-sequence star may have remained after the explosion.
Intensity mapping of neutral hydrogen (HI) is a promising observational probe of cosmology and large-scale structure. We present wide field simulations of HI intensity maps based on N-body simulations, the halo model, and a phenomenological prescription for assigning HI mass to halos. The simulations span a redshift range of 0.35 < z < 0.9 in redshift bins of width $\Delta z \approx 0.05$ and cover a quarter of the sky at an angular resolution of about 7'. We use the simulated intensity maps to study the impact of non-linear effects on the angular clustering of HI. We apply and compare several estimators for the angular power spectrum and its covariance. We verify that they agree with analytic predictions on large scales and study the validity of approximations based on Gaussian random fields, particularly in the context of the covariance. We discuss how our results and the simulated maps can be useful for planning and interpreting future HI intensity mapping surveys.
Using spectroscopic observations and photometric light curves of 238 nearby M dwarfs from the MEarth exoplanet transit survey, we examine the relationships between magnetic activity (quantified by H-alpha emission), rotation period, and stellar age. Previous attempts to investigate the relationship between magnetic activity and rotation in these stars were hampered by the limited number of M dwarfs with measured rotation periods (and the fact that vsini measurements probe only rapid rotation). However, the photometric data from MEarth allows us to probe a wide range of rotation periods for hundreds of M dwarf stars (from shorter than than one to longer than 100 days). Over all M spectral types that we probe, we find that the presence of magnetic activity is tied to rotation, including for late-type, fully convective M dwarfs. We also find evidence that the fraction of late-type M dwarfs that are active may be higher at longer rotation periods compared to their early-type counterparts, with several active, late-type, slowly rotating stars present in our sample. Additionally, we find that all M dwarfs with rotation periods shorter than 26 days (early-type; M1-M4) and 86 days (late-type; M5-M8) are magnetically active. This potential mismatch suggests that the physical mechanisms that connect stellar rotation to chromospheric heating may be different in fully convective stars. A kinematic analysis suggests that the magnetically active, rapidly rotating stars are consistent with a kinematically young population, while slow-rotators are less active or inactive and appear to belong to an older, dynamically heated stellar population.
We present a sample of 20 runaway M dwarf candidates (RdMs) within 1 kpc of the Sun whose Galactocentric velocities exceed 400 km s$^{-1}$. The candidates were selected from the SDSS DR7 M Dwarf Catalog of West et al. (2011). Our RdMs have SDSS+USNO-B proper motions that are consistent with those recorded in the PPMXL, LSPM, and combined WISE+SDSS+2MASS catalogs. Sixteen RdMs are classified as dwarfs, while the remaining four RdMs are subdwarfs. We model the Galactic potential using a bulge-disk-halo profile (Kenyon et al. 2008; Brown et al. 2014). Our fastest RdM, with Galactocentric velocity 658.5 $\pm$ 236.9 km s$^{-1}$, is a possible hypervelocity candidate, as it is unbound in 77% of our simulations. About half of our RdMs have kinematics that are consistent with ejection from the Galactic center. Seven of our RdMs have kinematics consistent with an ejection scenario from M31 or M32 to within 2{\sigma}, although our distance-limited survey makes such a realization unlikely. No more than four of our RdMs may have originated from the Leo stream. We propose that to within measurement errors, most of our bound RdMs are likely disk runaways or halo objects, and may have been accelerated through a series of multi-body interactions within the Galactic disk or possibly supernovae explosions.
Several types of extra-galactic high-energy transients have been discovered, which include high-luminosity and low-luminosity long-duration gamma-ray bursts (GRBs), short-duration GRBs, supernova shock breakouts (SBOs), and tidal disruption events (TDEs) without or with an associated relativistic jet. In this paper, we apply a unified method to systematically study the redshift-dependent event rate densities and the global luminosity functions (ignoring redshift evolution) of these transients. We introduce some empirical formulae for the redshift-dependent event rate densities for different types of transients, and derive the local specific event rate density, which also represents its global luminosity function. Long GRBs have a large enough sample to reveal features in the global luminosity function, which is best characterized as a triple power law. All the other transients are consistent with having a single power law luminosity function. The total event rate density depends on the minimum luminosity, and we obtain the following values in units of ${\rm Gpc^{-3}~yr^{-1}}$: $0.8^{+0.1}_{-0.1}$ for high-luminosity long GRBs above $ 10^{50}~{\rm erg~s^{-1}}$, $164^{+98}_{-65}$ for low-luminosity long GRBs above $5\times 10^{46}~{\rm erg~s^{-1}}$, $1.3^{+0.4}_{-0.3}$, $1.2^{+0.4}_{-0.3}$, and $3.3^{+1.0}_{-0.8}$ above $ 10^{50}~ {\rm erg~s^{-1}}$ for short GRBs with three different merger delay models (Gaussian, log-normal, and power law), $1.9^{+2.4}_{-1.2}\times 10^4$ above $ 10^{44}~{\rm erg~s^{-1}}$ for SBOs, $ 4.8^{+3.2}_{-2.1}\times10^2 $ for normal TDEs above $10^{44}~ {\rm erg~s^{-1}}$, and $0.03^{+0.04}_{-0.02}$ above $ 10^{48}~ {\rm erg~s^{-1}}$ for TDE jets as discovered by Swift. Intriguingly, the global luminosity functions of different kinds of transients, which cover over 12 orders of magnitude, are consistent with a single power law with an index of -1.6.
We report the detection of a strong Milky Way-type 2175 \AA$ $ extinction bump at $z$ = 2.1166 in the quasar spectrum towards SDSS J121143.42+083349.7 from the Sloan Digital Sky Survey (SDSS) Data Release 10. We conduct follow up observations with the Echelle Spectrograph and Imager (ESI) onboard the Keck-II telescope and the Ultraviolet and Visual Echelle Spectrograph (UVES) on the VLT. This 2175 \AA$ $ absorber is remarkable in that we simultaneously detect neutral carbon (C I), neutral chlorine (Cl I), and carbon monoxide (CO). It also qualifies as a damped Lyman alpha system. The J1211+0833 absorber is found to be metal-rich and has a dust depletion pattern resembling that of the Milky Way disk clouds. We use the column densities of the C I fine structure states and the C II/C I ratio (under the assumption of ionization equilibrium) to derive the temperature and volume density in the absorbing gas. A Cloudy photoionization model is constructed, which utilizes additional atoms/ions to constrain the physical conditions. The inferred physical conditions are consistent with a canonical cold (T $\sim$ 100 K) neutral medium with a high density ($n$(H I) $\sim$ 100 cm$^{-3}$) and a slightly higher pressure than the local interstellar medium. Given the simultaneous presence of C I, CO, and the 2175 \AA$ $ bump, combined with the high metallicity, high dust depletion level and overall low ionization state of the gas, the absorber towards J1211+0833 supports the scenario that the presence of the bump requires an evolved stellar population.
Direct imaging of extrasolar planets with future space-based coronagraphic telescopes may provide a means of detecting companion moons at wavelengths where the moon outshines the planet. We propose a detection strategy based on the positional variation of the center of light with wavelength, "spectroastrometry." This new application of this technique could be used to detect an exomoon, to determine the exomoon's orbit and the mass of the host exoplanet, and to disentangle of the spectra of the planet and moon. We consider two model systems, for which we discuss the requirements for detection of exomoons around nearby stars. We simulate the characterization of an Earth-Moon analog system with spectroastrometry, showing that the orbit, the planet mass, and the spectra of both bodies can be recovered. To enable the detection and characterization of exomoons we recommend that coronagraphic telescopes should extend in wavelength coverage to 3 micron, and should be designed with spectroastrometric requirements in mind.
We extend our effort to calibrate stellar isochrones in the Johnson-Cousins ($BVI_C$) and the 2MASS ($JHK_s$) filter systems based on observations of well-studied open clusters. Using cool main-sequence (MS) stars in Praesepe, we define empirical corrections to the Lejeune et al. color-effective temperature ($T_{\rm eff}$) relations down to $T_{\rm eff} \sim 3600$ K, complementing our previous work based on the Hyades and the Pleiades. We apply empirically corrected isochrones to existing optical and near-infrared photometry of cool ($T_{\rm eff} \leq 5500$ K) and metal-rich ([Fe/H]$=+0.37$) MS stars in NGC 6791. The current methodology relies on an assumption that color-$T_{\rm eff}$ corrections are independent of metallicity, but we find that estimates of color-excess and distance from color-magnitude diagrams with different color indices converge on each other at the precisely known metallicity of the cluster. Along with a satisfactory agreement with eclipsing binary data in the cluster, we view the improved internal consistency as a validation of our calibrated isochrones at super-solar metallicities. For very cool stars ($T_{\rm eff} \leq 4800$ K), however, we find that $B - V$ colors of our models are systematically redder than the cluster photometry by $\sim0.02$ mag. We use color-$T_{\rm eff}$ transformations from the infrared flux method (IRFM) and alternative photometry to examine a potential color-scale error in the input cluster photometry. After excluding $B - V$ photometry of these cool MS stars, we derive $E(B - V)=0.105\pm0.014$, [M/H]$=+0.42\pm0.07$, $(m - M)_0 = 13.04\pm0.08$, and the age of $9.5\pm0.3$ Gyr for NGC 6791.
Characterization of transiting planets with transit timing variations (TTVs) requires understanding how to translate the observed TTVs into masses and orbital elements of the planets. This can be challenging in multi-planet transiting systems, but fortunately these systems tend to be nearly plane-parallel and low eccentricity. Here we present a novel derivation of analytic formulae for TTVs that are accurate to first order in the planet-star mass ratios and in the orbital eccentricities. These formulae are accurate in proximity to first order resonances, as well as away from resonance, and compare well with more computationally expensive N-body integrations in the low eccentricity, low mass-ratio regime when applied to simulated and to actual multi-transiting Kepler planet systems. We make code available for implementing these formulae.
We present visible spectroscopic and albedo data of the 2.3 Gyr old Themis family and the <10 Myr old Beagle sub-family. The slope and albedo variations between these two families indicate C-complex asteroids become redder and darker in response to space weathering. Our observations of Themis family members confirm previously observed trends where phyllosilicate absorption features are less common among small diameter objects. Similar trends in the albedos of large (> 15 km) and small (< 15 km) Themis members suggest these phyllosilicate feature and albedo trends result from regolith variations as a function of diameter. Observations of the Beagle asteroids show a small, but notable fraction of members with phyllosilicate features. The presence of phyllosilicates and the dynamical association of the main-belt comet 133P/Elst-Pizarro with the Beagle family imply the Beagle parent body was a heterogenous mixture of ice and aqueously altered minerals.
We present Keck/NIRC2 and OSIRIS near-infrared imaging and spectroscopy of 2M0441+2301 AabBab, a young (1--3 Myr) hierarchical quadruple system comprising a low-mass star, two brown dwarfs, and a planetary-mass companion in Taurus. All four components show spectroscopic signs of low surface gravity, and both 2M0441+2301 Aa and Ab possess Pa$\beta$ emission indicating they each harbor accretion subdisks. Astrometry spanning 2008--2014 reveals orbital motion in both the Aab (0.23" separation) and Bab (0.095" separation) pairs, although the implied orbital periods of $>$300 years means dynamical masses will not be possible in the near future. The faintest component (2M0441+2301 Bb) has an angular $H$-band shape, strong molecular absorption (VO, CO, H$_2$O, and FeH), and shallow alkali lines, confirming its young age, late spectral type (L1 $\pm$ 1), and low temperature ($\approx$1800~K). With individual masses of 200$^{+100}_{-50}$ Mjup, 35 $\pm$ 5 Mjup, 19 $\pm$ 3 Mjup, and 9.8 $\pm$ 1.8 Mjup, 2M0441+2301 AabBab is the lowest-mass quadruple system known. Its hierarchical orbital architecture and mass ratios imply that it formed from the collapse and fragmentation of a molecular cloud core, demonstrating that planetary-mass companions can originate from a stellar-like pathway analogous to higher-mass quadruple star systems. More generally, cloud fragmentation may be an important formation pathway for the massive exoplanets that are now regularly being imaged on wide orbits.
The paper presents the basic principles of formation of a database (DB) with information about objects and their physical characteristics from observations carried out at the Crimean Astrophysical Observatory (CrAO) and published in "Izvestiya Krymskoi Astrofizicheskoi Observatorii" and other publications. The emphasis is placed on DBs that are not present in the most complete global library catalogs and data tables - VizieR (supported by the Strasbourg ADC). Separately, we consider the formation of a digital archive of observational data obtained at CrAO - as the interactive DB related to the DB of objects and publications. Examples of all the above DB as elements integrated into the Crimean Astronomical Virtual Observatory are presented in the paper. The operation with CrAO database is illustrated using tools of the International Virtual Observatory - Aladin, VOPlot, VOSpec jointly with VizieR DB and Simbad.
Numerical simulations are a versatile tool providing insight into the complicated process of structure formation in cosmology. This process is mainly governed by gravity, which is the dominant force on large scales. To date, a century after the formulation of general relativity, numerical codes for structure formation still employ Newton's law of gravitation. This approximation relies on the two assumptions that gravitational fields are weak and that they are only sourced by non-relativistic matter. While the former appears well justified on cosmological scales, the latter imposes restrictions on the nature of the "dark" components of the Universe (dark matter and dark energy) which are, however, poorly understood. Here we present the first simulations of cosmic structure formation using equations consistently derived from general relativity. We study in detail the small relativistic effects for a standard {\Lambda}CDM cosmology which cannot be obtained within a purely Newtonian framework. Our particle-mesh N-body code computes all six degrees of freedom of the metric and consistently solves the geodesic equation for particles, taking into account the relativistic potentials and the frame-dragging force. This conceptually clean approach is very general and can be applied to various settings where the Newtonian approximation fails or becomes inaccurate, ranging from simulations of models with dynamical dark energy or warm/hot dark matter to core collapse supernova explosions.
Abundance of iron in the atmosphere of Arcturus has been determined from the profiles or regions of the profiles of the weak lines sensitive to iron abundance. The selected lines of Fe I and Fe II were synthesized with the MARCS theoretical models of the atmosphere. From the observed profiles of lines available with a high spectral resolution in the atlas by Hinkle and Wallace (2005), the values of the iron abundance $A = 6.95 \pm 0.03$ and the radial-tangential macroturbulent velocity $5.6 \pm 0.2$ km/s were obtained for Arcturus. The same physical quantities were found for the Sun as a star; they are $7.42 \pm 0.02$ and $3.4 \pm 0.3$ km/s, respectively. For Arcturus, the iron abundance relative to the solar one was determined with the differential method as [Fe/H] $=-0.48 \pm 0.02$.
The determination of the brown dwarf binary fraction may contribute to the
understanding of the substellar formation mechanisms. Unresolved brown dwarf
binaries may be revealed through their peculiar spectra or the discrepancy
between optical and near-infrared spectral type classification.
We obtained medium-resolution spectra of 22 brown dwarfs with these
characteristics using the X-Shooter spectrograph at the VLT.
We aimed to identify brown dwarf binary candidates, and to test if the
BT-Settl 2014 atmospheric models reproduce their observed spectra.
To find binaries spanning the L-T boundary, we used spectral indices and
compared the spectra of the selected candidates to single spectra and synthetic
binary spectra. We used synthetic binary spectra with components of same
spectral type to determine as well the sensitivity of the method to this class
of binaries.
We identified three candidates to be combination of L plus T brown dwarfs. We
are not able to identify binaries with components of similar spectral type. In
our sample, we measured minimum binary fraction of $9.1^{+9.9}_{-3.0}$.
From the best fit of the BT-Settl models 2014 to the observed spectra, we
derived the atmospheric parameters for the single objects. The BT-Settl models
were able to reproduce the majority of the SEDs from our objects, and the
variation of the equivalent width of the RbI (794.8 nm) and CsI (852.0 nm)
lines with the spectral type. Nonetheless, these models did not reproduce the
evolution of the equivalent widths of the NaI (818.3 nm and 819.5 nm) and KI
(1253 nm) lines with the spectral type.
We present extensive optical (UBVRI, g'r'i'z', and open CCD) and near-infrared (ZYJH) photometry for the very nearby Type IIP supernova (SN) 2013ej extending from +1 to +215 days after shock breakout. Substantial time series ultraviolet and optical spectroscopic observations obtained from +8 to +135 days after the explosion are also presented. We use early optical photometry to determine the time of shock breakout to be MJD 56496.9+/-0.3. We calibrate unfiltered CCD observations to pseudo-bolometric UBVRI light curve data of Type IIP supernovae (SNe), obtaining 4% precision with a B-V color-dependent correction. Considering well-observed SNe from the literature, we derive a UBVRIJHK bolometric calibration from the UBVRI measurement to better than 2% precision. We observe moderately strong Si~II {\lambda}6355 as early as 8 days after explosion. We model spectra in the vicinity of Fe~II {\lambda}5169 whenever observed to determine the photospheric velocity, v{\phi}, yielding 4500+/-500 km s^-1 50 days after explosion. We observe similar early ultraviolet spectra at 10--12 days after explosion for SNe~IIP, but see diverse behavior several days earlier. Using the expanding photosphere method, and combining with SN 2002ap, we estimate the host galaxy distance to be 9.0^-0.6_+0.4 Mpc, consistent with estimates from other methods. Photometric and spectroscopic analysis during the plateau phase, which we estimate to be 86+/-6 days long, yields an estimated explosion energy of 1.05+/-0.43 x 10^51 ergs, a total ejected mass of 14.3+/-4.5 MSol, and a pre-supernova radius of 200+/-82 RSol. Measurements from the radioactive tail phase yield a 56Ni mass of 0.019+/-0.001 MSol.
We report the first identification of a compact dense galaxy overdensity at $z\sim 8$ dubbed A2744z8OD. A2744z8OD consists of eight $Y$-dropout galaxies behind Abell 2744 that is originally pinpointed by Hubble Frontier Fields studies. However, no studies have, so far, derived basic physical quantities of structure formation or made comparisons with theoretical models. We obtain a homogeneous sample of dropout galaxies at $z\sim 8$ from eight field data of Hubble legacy images that are as deep as the A2744z8OD data. Using the sample, we find that a galaxy surface overdensity value of A2744z8OD is very high $\delta\simeq 130$ that is defined in a small circle of $6"$ ($\simeq 30$ physical kpc) radius. Because there is no such a large $\delta$ value reported for high-$z$ overdensities to date, A2744z8OD is a system clearly different from those found in previous high-$z$ overdensity studies. The total stellar mass of A2744z8OD is estimated to be $3.8\times 10^9 M_\odot$ that is as small as today's Large Magellanic Cloud. In the galaxy+structure formation models of Henriques et al. (2015), there exist a very similar overdensity, Modelz8OD, that is made of eight model dropout galaxies at $z\sim 8$ in a $6"$-radius circle. Eight out of seven galaxies in Modelz8OD form a filament within $\Delta z \sim 0.03$ that is elongated in the line of sight, which enhance $\delta$ of Modelz8OD. Modelz8OD is a progenitor of a today's cluster with $10^{14} M_\odot$, and more than a half of the seven Modelz8OD galaxies are merged into the brightest cluster galaxy (BCG) in the today's cluster. If Modelz8OD is a counterpart of A2744z8OD, the models suggest that A2744z8OD would be a forming cluster core of a today's $10^{14} M_\odot$ cluster that started formation earlier than the most of the other BCG progenitors at $z>12$.
The Extremely Large Telescopes currently under construction have a collecting area that is an order of magnitude larger than the present largest optical telescopes. For seeing-limited observations the performance will scale as the collecting area but, with the successful use of adaptive optics, for many applications it will scale as $D^4$ (where $D$ is the diameter of the primary mirror). Central to the success of the ELTs, therefore, is the successful use of multi-conjugate adaptive optics (MCAO) that applies a high degree correction over a field of view larger than the few arcseconds that limits classical adaptive optics systems. In this letter, we report on the analysis of crowded field images taken on the central region of the Galactic globular cluster NGC 1851 in $K_s$ band using GeMS at the Gemini South telescope, the only science-grade MCAO system in operation. We use this cluster as a benchmark to verify the ability to achieve precise near-infrared photometry by presenting the deepest $K_s$ photometry in crowded fields ever obtained from the ground. We construct a colour-magnitude diagram in combination with the F606W band from HST/ACS. As well as detecting the "knee" in the lower main sequence at $K_s\simeq20.5$, we also detect the double subgiant branch of NGC 1851, that demonstrates the high photometric accuracy of GeMS in crowded fields.
In gamma ray astronomy, the energy range from sub-100GeV to TeV is crucial due to where there is a gap between space experiments and ground-based ones. In addition, observations in this energy range are expected to provide more details about the high energy emission from GRBs,and thus to understand EBL. Based on the observation results and the related knowledge, scientists may be able to unveil the mysteries of galaxy formation and the evolution of early universe. One of the principal issues for next generation Imaging Atmospheric Cherenkov Telescopes (IACT) is to achieve larger field of view (FoV). In this work, we report a refractive water convex lens as light collector to test the feasibility of a new generation of IACT, and some preliminary test results on the optical properties (the focal length, spot size, transmittance, etc.) of a 0.9 m diameter water lens, the photodetectors and DAQ system of a prototype are presented and discussed.
We compare the absolute gain photometric calibration of the Planck/HFI and Herschel/SPIRE instruments on diffuse emission. The absolute calibration of each of HFI and SPIRE relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 $\mu$m and 857 GHz / 350 $\mu$m. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^2. We have selected these fields carefully to provide a high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. SPIRE maps are convolved to match the HFI beam and put on a common pixel grid. We measure the relative gain between the instruments using two methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. Adopting the current SPIRE point-source to extended-emission conversion, we find that the two calibrations are in very good agreement. The SPIRE / HFI relative gains are 1.047 ($\pm$ 0.0069) and 1.003 ($\pm$ 0.0080) at 545 and 857 GHz, respectively. These relative gains deviate from unity by much less than the current uncertainty of the absolute extended emission calibration, which is about 6.4% and 9.5% for HFI and SPIRE, respectively, but the deviations are comparable to the values 1.4% and 5.5% for HFI and SPIRE if the uncertainty from models of the common calibrator can be discounted. Of the 5.5% for SPIRE, 4% arises from the beam area, highlighting that as focus for refinement.
In this paper, we present observations and analysis of an interesting sigmoid formation, eruption and the associated flare that occurred on 2014 April 18 using multi-wavelength data sets. We discuss the possible role of the sigmoid eruption in triggering the flare, which consists of two different set of ribbons: parallel ribbons as well as a large-scale quasi-circular ribbon. Several observational evidence and nonlinear force-free field extrapolation results show the existence of a large-scale fan-spine type magnetic configuration with a sigmoid lying under a section of the fan dome. The event can be explained with the following two phases. During the pre-flare phase, we observed the formation and appearance of sigmoid via tether-cutting reconnection between the two sets of sheared fields under the fan dome. The second, main flare phase, features the eruption of the sigmoid, the subsequent flare with parallel ribbons, and a quasi-circular ribbon. We propose the following multi-stage successive reconnections scenario for the main flare. First, tether-cutting reconnection is responsible for the formation and the eruption of the sigmoid structure. Second, the reconnection occurred in the wake of the erupting sigmoid produces the parallel flare ribbons on the both sides of the circular polarity inversion line. Third, the null-type reconnection higher in the corona, possibly triggered by the erupting sigmoid, leads to the formation of a large quasi-circular ribbon. For the first time we suggest a mechanism for this type of flare consisting of a double set of ribbons triggered by an erupting sigmoid in a large scale fan-spine type magnetic configuration.
M82 X-1 is the brightest ultraluminous X-ray source in starburst galaxy M82 and is one of the best intermediate mass black hole candidates. Previous studies based on the Rossi X-ray Timing Explorer/Proportional Counter Array (RXTE/PCA) reported a regular X-ray flux modulation of M82 with a period of 62 days, and attributed this periodic modulation to M82 X-1. However, this modulation is not necessarily from M82 X-1 because RXTE/PCA has a very poor spatial resolution of ~1 degree. In this work, we analyzed 1000 days of monitoring data of M82 X-1 from the Swift/X-ray telescope (XRT), which has a much better spatial resolution than RXTE/PCA. The periodicity distribution map of M82 reveals that the 62-day periodicity is most likely not from M82 X-1, but from the summed contributions of several periodic X-ray sources 4 arcsec southeast of M82 X-1. However, Swift/XRT is not able to resolve those periodic sources and locate the precise origin of the periodicity of M82. Thus, more long-term observations with higher spatial resolution are required.
Two major features of the prestellar CMF are: 1) a broad peak below 1 Msun, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the CMF/IMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of filaments by interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in nearby molecular clouds observed with Herschel using a 1-D power spectrum analysis. The observed filament power spectra were fitted by a power-law function $(P_{true}(s) \propto s^{\alpha})$ after removing the effect of beam convolution at small scales. A Gaussian-like distribution of power-spectrum slopes was found centered at -1.6, close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). An empirical correlation, $P^{0.5}(s_0) \propto <N_{\rm H_2}>^{1.4 \pm 0.1} $, was also found between the amplitude of each filament power spectrum $P(s_0)$ and the mean column density along the filament $<N_{\rm H_2}>$. Finally, the dispersion of line-mass fluctuations along each filament $\sigma_{\rm M_{line}}$ was found to scale with the physical length $L$ of the filament, roughly as $\sigma_{M_{line}} \propto L^{0.7}$. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the gravitational fragmentation of filaments.
A single-mirror small-size (SST-1M) Davies-Cotton telescope with a dish diameter of 4 m has been built by a consortium of Polish and Swiss institutions as a prototype for one of the proposed small-size telescopes for the southern observatory of the Cherenkov Telescope Array (CTA). The design represents a very simple, reliable, and cheap solution. The mechanical structure prototype with its drive system is now being tested at the Institute of Nuclear Physics PAS in Krakow. Here we present the design of the prototype and results of the performance tests of the structure and the drive and control system.
We present moderate resolution CCD spectra and R photometry for seven KP2001 stars. We revised the spectral classification of the stars in the range 3900 - 8500 A. On the bases of light curves of the NSVS (Northern Sky Variability Survey) database we classify KP2001-18 as a semi regular and KP2001-176 as Mira type variables. For all observed objects NSVS phase - dependence light curve analysis and variability type classification was performed with the VStar Software. Using period luminosity relation we computed MK magnitudes and the distances to variables.
A recent claim by Lieu et al that beam splitter intensity subtraction (or homodyne with one vacuum port) followed by high resolution sampling can lead to detection of brightness of thermal light at the shot noise limit is reexamined here. We confirm the calculation of Zmuidzinas that the claim of Lieu et al was falsified by an incorrect assumption about the correlations in thermal noise.
We perform two-dimensional hydrodynamic simulations for the thermonuclear explosion of Chandrasekhar-mass white dwarfs with dark matter (DM) cores in Newtonian gravity. We include a 19-isotope nuclear reaction network and make use of the pure turbulent deflagration model as the explosion mechanism in our simulations. Our numerical results show that the general properties of the explosion depend quite sensitively on the mass of the DM core M$_{{\rm DM}}$: a larger M$_{{\rm DM}}$ generally leads to a weaker explosion and a lower mass of synthesized iron-peaked elements. In particular, the total mass of $^{56}$Ni produced can drop from about 0.3 to 0.03 $M_{\odot}$ as M$_{{\rm DM}}$ increases from 0.01 to 0.03 $M_{\odot}$. We have also constructed the bolometric light curves obtained from our simulations and found that our results match well with the observational data of sub-luminous Type-Ia supernovae.
Quasar catalog analyses have allowed to identify a huge large quasar group
(HLQG) with a mean redshift of 1.27 and characteristic comoving size of $\sim
500$ Mpc.
Because of some controversy about its actual existence we look for an
independent evidence of its gravitational field, in search for its effects on
the cosmic microwave background (CMB) photons propagating through it. We
analyze the CMB Planck temperature map of the region of sky corresponding to
the angular position of the HLQG, in search of an independent confirmation that
it is a real structure and that it actually marks a very extended mass
overdensity. We compute an inner and an outer temperature by averaging the CMB
map over respectively the region subtended by the HLQG on the sky, and over a
region surrounding it. It turns out that the inner region is warmer than the
outer one, with the measured temperature difference $\Delta T^{obs} \approx
43\mu K$. The temperature excess was then compared with the results of
Montecarlo simulations of random gaussian realizations of the CMB map
indicating, with a $2.5\sigma$ confidence level, that the measured $\Delta
T^{obs}$ for the HLQG cannot be attributed to primordial gaussian fluctuations.
This temperature anomaly was not detected before due to the irregular shape of
the inner region, and the need to compare the inner and outer temperature to
identify it correctly. Its angular extension in the longitudinal direction is
about three times that of the cold spot (CS), while its total angular area is
comparable, due to its elongated shape compared to the CS. Our results are
stable under the choice of different simulations methods and different
definitions of the inner and outer regions.
We present the detection of compact radio structures of fourteen radio-loud narrow line Seyfert 1 (NLS1) galaxies from Very Long Baseline Array observations at 5 GHz, which were performed in 2013. While 50\% of the sources of our sample show a compact core only, the remaining 50\% exhibit a core-jet structure. The measured brightness temperatures of the cores range from $10^{8.4}$ to $10^{11.4}$ K with a median value of $10^{10.1}$ K, indicating that the radio emission is from non-thermal jets, and that, likely, most sources are not strongly beamed, then implying a lower jet speed in these radio-loud NLS1 galaxies. In combination with archival data taken at multiple frequencies, we find that seven sources show flat or even inverted radio spectra, while steep spectra are revealed in the remaining seven objects. Although all these sources are very radio-loud with $R > 100$, their jet properties are diverse, in terms of their milli-arcsecond (mas) scale (pc scale) morphology and their overall radio spectral shape. The evidence for slow jet speeds (i.e., less relativistic jets), in combination with the low kinetic/radio power, may offer an explanation for the compact VLBA radio structure in most sources. The mildly relativistic jets in these high accretion rate systems are consistent with a scenario, where jets are accelerated from the hot corona above the disk by the magnetic field and the radiation force of the accretion disk. Alternatively, a low jet bulk velocity can be explained by low spin in the Blandford-Znajek mechanism.
We present the results of new Suzaku observations of the Coma Cluster, the X-ray brightest, nearby, merging system hosting a well studied, typical giant radio halo. It has been previously shown that, on the western side of the cluster, the radio brightness shows a much steeper gradient compared to other azimuths. XMM-Newton and Planck revealed a shock front along the southern half of the region associated with this steep radio gradient, suggesting that the radio emission is enhanced by particle acceleration associated with the shock passage. Suzaku demonstrates for the first time that this shock front extends northwards, tracing the entire length of the western edge of the Coma radio halo. The shock is detected both in the temperature and X-ray surface brightness distributions and has a Mach number of around $\mathcal{M}\sim1.5$. The locations of the surface brightness edges align well with the edge of the radio emission, while the obtained temperature profiles seem to suggest shocks located 125-185 kpc further out in radius. In addition, the shock strengths derived from the temperature and density jumps are in agreement when using extraction regions parallel to the radio halo edge, but become inconsistent with each other when derived from radial profiles centred on the Coma Cluster core. It is likely that, beyond mere projection effects, the geometry of the shock is more complex than a front with a single, uniform Mach number and an approximately spherically symmetric shape.
We present a photometric catalog of 8,735,004 proper motion selected low-mass stars (KML-spectral types) within the Sloan Digital Sky Survey (SDSS) footprint, from the combined SDSS Data Release 10 (DR10), Two-Micron All-Sky Survey (2MASS) Point Source Catalog (PSC), and Wide-field Infrared Survey Explorer (WISE) AllWISE catalog. Stars were selected using $r-i$, $i-z$, $r-z$, $z-J$, and $z-W1$ colors, and SDSS, WISE, and 2MASS astrometry was combined to compute proper motions. The resulting 3,518,150 stars were augmented with proper motions for 5,216,854 earlier type stars from the combined SDSS and United States Naval Observatory B1.0 catalog (USNO-B). We used SDSS+USNO-B proper motions to determine the best criteria for selecting a clean sample of stars. Only stars whose proper motions were greater than their $2$$\sigma$ uncertainty were included. Our Motion Verified Red Stars (MoVeRS) catalog is available through SDSS CasJobs and VizieR.
A few particles of presolar Al2O3 grains with sizes above 0.5 mum are believed to have been produced in the ejecta of core-collapse supernovae (SNe). In order to clarify the formation condition of such large Al2O3 grains, we investigate the condensation of Al2O3 grains for wide ranges of the gas density and cooling rate. We first show that the average radius and condensation efficiency of newly formed Al2O3 grains are successfully described by a non-dimensional quantity "Lambda_on" defined as the ratio of the timescale with which the supersaturation ratio increases to the collision timescale of reactant gas species at dust formation. Then, we find that the formation of submicron-sized Al2O3 grains requires at least ten times higher gas densities than those presented by one-dimensional SN models. This indicates that presolar Al2O3 grains identified as a SN origin might be formed in dense gas clumps, allowing us to propose that the measured sizes of presolar grains can be a powerful tool to constrain the physical conditions in which they formed. We also briefly discuss the survival of newly formed Al2O3 grains against the destruction in the shocked gas within the SN remnants.
VER J0521+211 (RGB J0521.8+2112) is one of the brightest and most powerful blazars detected in the TeV gamma-ray regime. It is located at a redshift of z=0.108 and since its discovery in 2009, VER J0521+211 has exhibited an average TeV flux exceeding 0.1 times that of the Crab Nebula, corresponding to an isotropic luminosity of $3\times10^{44}$ erg s$^{-1}$. We present data from a comprehensive multiwavelength campaign on this object extending between November 2012 and February 2014, including single-dish radio observations, optical photometry and polarimetry, UV, X-ray, GeV and TeV gamma-ray data (VERITAS, MAGIC). Significant flux variability was observed at all wavelengths, including a long-lasting high state at gamma-ray energies in Fall 2013. Nightly-resolved spectra at X-ray and TeV energies are be presented, and emission mechanisms explaining the observed flux and spectral variability are discussed.
Carbon monoxide (CO) is widely used as a tracer of molecular hydrogen (H2) in metal-rich galaxies, but is known to become ineffective in low metallicity dwarf galaxies. Atomic carbon has been suggested as a superior tracer of H2 in these metal-poor systems, but its suitability remains unproven. To help us to assess how well atomic carbon traces H2 at low metallicity, we have performed a series of numerical simulations of turbulent molecular clouds that cover a wide range of different metallicities. Our simulations demonstrate that in star-forming clouds, the conversion factor between [CI] emission and H2 mass, $X_{\rm CI}$, scales approximately as $X_{\rm CI} \propto Z^{-1}$. We recover a similar scaling for the CO-to-H2 conversion factor, $X_{\rm CO}$, but find that at this point in the evolution of the clouds, $X_{\rm CO}$ is consistently smaller than $X_{\rm CI}$, by a factor of a few or more. We have also examined how $X_{\rm CI}$ and $X_{\rm CO}$ evolve with time. We find that $X_{\rm CI}$ does not vary strongly with time, demonstrating that atomic carbon remains a good tracer of H2 in metal-poor systems even at times significantly before the onset of star formation. On the other hand, $X_{\rm CO}$ varies very strongly with time in metal-poor clouds, showing that CO does not trace H2 well in starless clouds at low metallicity.
Near-infrared observations of stellar orbits at the Galactic Center provide conclusive evidence for a massive black hole associated with the compact radio source Sgr A*. The astrometric reference frame for these observations is tied to a set of red giant stars, which are also detectable at radio wavelengths through SiO maser emission in their envelopes. We have improved the precision and long-term stability of this reference frame, in which Sgr A* is localized to within a factor 5 better than previously: ~0.17 mas in position (in 2009) and ~0.07 mas/yr in velocity. This improvement is the result of modeling and correcting optical distortion in the VLT/NACO imager to a sub-mas level and including new infrared and radio measurements, which now both span more than a decade in time. A further improvement will follow future observations and facilitate the detection of relativistic orbital effects.
The Cherenkov Telescope Array (CTA), the new generation very high-energy gamma-ray observatory, will improve the flux sensitivity of the current Cherenkov telescopes by an order of magnitude over a continuous range from about 10 GeV to above 100 TeV. With tens of telescopes distributed in the Northern and Southern hemispheres, the large effective area and field of view coupled with the fast pointing capability make CTA a crucial instrument for the detection and understanding of the physics of transient, short-timescale variability phenomena (e.g. Gamma-Ray Bursts, Active Galactic Nuclei, gamma-ray binaries, serendipitous sources). The key CTA system for the fast identification of flaring events is the Real-Time Analysis (RTA) pipeline, a science alert system that will automatically detect and generate science alerts with a maximum latency of 30 seconds with respect to the triggering event collection and ensure fast communication to/from the astrophysics community. According to the CTA design requirements, the RTA search for a true transient event should be performed on multiple time scales (from minutes to hours) with a sensitivity not worse than three times the nominal CTA sensitivity. Given the CTA requirement constraints on the RTA efficiency and the fast response ability demanded by the transient science, we perform a preliminary evaluation of the RTA sensitivity as a function of the CTA high-level technical performance (e.g. effective area, point spread function) and the observing time. This preliminary approach allows the exploration of the complex parameter space defined by the scientific and technological requirements, with the aim of defining the feasibility range of the input parameters and the minimum background rejection capability of the RTA pipeline.
The Cherenkov Telescope Array (CTA) is the next generation observatory for the study of very high-energy gamma rays from about 20 GeV up to 300 TeV. Thanks to the large effective area and field of view, the CTA observatory will be characterized by an unprecedented sensitivity to transient flaring gamma-ray phenomena compared to both current ground (e.g. MAGIC, VERITAS, H.E.S.S.) and space (e.g. Fermi) gamma-ray telescopes. In order to trigger the astrophysics community for follow-up observations, or being able to quickly respond to external science alerts, a fast analysis pipeline is crucial. This will be accomplished by means of a Real-Time Analysis (RTA) pipeline, a fast and automated science alert trigger system, becoming a key system of the CTA observatory. Among the CTA design key requirements to the RTA system, the most challenging is the generation of alerts within 30 seconds from the last acquired event, while obtaining a flux sensitivity not worse than the one of the final analysis by more than a factor of 3. A dedicated software and hardware architecture for the RTA pipeline must be designed and tested. We present comparison of OpenCL solutions using different kind of devices like CPUs, Graphical Processing Unit (GPU) and Field Programmable Array (FPGA) cards for the Real-Time data reduction of the Cherenkov Telescope Array (CTA) triggered data.
The Cherenkov Telescope Array (CTA) observatory will be one of the largest ground-based very high-energy gamma-ray observatories. The On-Site Analysis will be the first CTA scientific analysis of data acquired from the array of telescopes, in both northern and southern sites. The On-Site Analysis will have two pipelines: the Level-A pipeline (also known as Real-Time Analysis, RTA) and the level-B one. The RTA performs data quality monitoring and must be able to issue automated alerts on variable and transient astrophysical sources within 30 seconds from the last acquired Cherenkov event that contributes to the alert, with a sensitivity not worse than the one achieved by the final pipeline by more than a factor of 3. The Level-B Analysis has a better sensitivity (not be worse than the final one by a factor of 2) and the results should be available within 10 hours from the acquisition of the data: for this reason this analysis could be performed at the end of an observation or next morning. The latency (in particular for the RTA) and the sensitivity requirements are challenging because of the large data rate, a few GByte/s. The remote connection to the CTA candidate site with a rather limited network bandwidth makes the issue of the exported data size extremely critical and prevents any kind of processing in real-time of the data outside the site of the telescopes. For these reasons the analysis will be performed on-site with infrastructures co-located with the telescopes, with limited electrical power availability and with a reduced possibility of human intervention. This means, for example, that the on-site hardware infrastructure should have low-power consumption. A substantial effort towards the optimization of high-throughput computing service is envisioned to provide hardware and software solutions with high-throughput, low-power consumption at a low-cost.
For the proposed Cherenkov Telescope Array (CTA), a post-calibration point-source location accuracy of 3 seconds of arc is aimed for under favorable observing conditions and for gamma-ray energies exceeding 100 GeV. In this contribution, results of first studies on the location accuracy are presented. These studies are based on a toy Monte Carlo simulation of a typical CTA-South array layout, taking into account the expected trigger rates of the different CTA telescope types and the gamma-ray spectrum of the simulated source. With this simulation code it is possible to study the location accuracy as a function of arbitrary telescope mis-orientations and for typical observing patterns on the sky. Results are presented for various scenarios, including one for which all individual telescopes are randomly mis-oriented within their specified limits. The study provides solid lower limits for the expected source location accuracy of CTA, and can be easily extended to include various other important effects like atmospheric refraction or partial cloud coverage.
The Cherenkov Telescope Array is a world-wide project for a new generation of ground-based Cherenkov telescopes of the Imaging class with the aim of exploring the highest energy region of the electromagnetic spectrum. With two planned arrays, one for each hemisphere, it will guarantee a good sky coverage in the energy range from a few tens of GeV to hundreds of TeV, with improved angular resolution and a sensitivity in the TeV energy region better by one order of magnitude than the currently operating arrays. In order to cover this wide energy range, three different telescope types are envisaged, with different mirror sizes and focal plane features. In particular, for the highest energies a possible design is a dual-mirror Schwarzschild-Couder optical scheme, with a compact focal plane. A silicon photomultiplier (SiPM) based camera is being proposed as a solution to match the dimensions of the pixel (angular size of ~ 0.17 degrees). INFN is developing a camera demonstrator made by 9 Photo Sensor Modules (PSMs, 64 pixels each, with total coverage 1/4 of the focal plane) equipped with FBK (Fondazione Bruno Kessler, Italy) Near UltraViolet High Fill factor SiPMs and Front-End Electronics (FEE) based on a Target 7 ASIC, a 16 channels fast sampler (up to 2GS/s) with deep buffer, self-trigger and on-demand digitization capabilities specifically developed for this purpose. The pixel dimensions of $6\times6$ mm$^2$ lead to a very compact design with challenging problems of thermal dissipation. A modular structure, made by copper frames hosting one PSM and the corresponding FEE, has been conceived, with a water cooling system to keep the required working temperature. The actual design, the adopted technical solutions and the achieved results for this demonstrator are presented and discussed.
We present a sophisticated likelihood reconstruction algorithm for shower-image analysis of imaging Cherenkov telescopes. The reconstruction algorithm is based on the comparison of the camera pixel amplitudes with the predictions from a Monte Carlo based model. Shower parameters are determined by a maximisation of a likelihood function. Maximisation of the likelihood as a function of shower fit parameters is performed using a numerical non-linear optimisation technique. A related reconstruction technique has already been developed by the CAT and the H.E.S.S. experiments, and provides a more precise direction and energy reconstruction of the photon induced shower compared to the second moment of the camera image analysis. Examples are shown of the performance of the analysis on simulated gamma-ray data from the VERITAS array.
We present confirmation of the cluster MOO J1142+1527, a massive galaxy cluster discovered as part of the Massive and Distant Clusters of WISE Survey. The cluster is confirmed to lie at $z=1.19$, and using the Combined Array for Research in Millimeter-wave Astronomy we robustly detect the Sunyaev-Zel'dovich (SZ) decrement at 13.2$\sigma$. The SZ data imply a mass of $\mathrm{M}_{200m}=(1.1\pm0.2)\times10^{15}$ $\mathrm{M}_\odot$, making MOO J1142+1527 the most massive galaxy cluster known at $z>1.15$ and the second most massive cluster known at $z>1$. For a standard $\Lambda$CDM cosmology it is further expected to be one of the $\sim 5$ most massive clusters expected to exist at $z\ge1.19$ over the entire sky. Our ongoing Spitzer program targeting $\sim1750$ additional candidate clusters will identify comparably rich galaxy clusters over the full extragalactic sky.
Previously, [1207.0458] proposed an astrophysical explanation of narrow gamma-ray line-like feature(s) at ~100 GeV from Galactic Center region observed by Fermi/LAT [1205.4700]. The model of [1207.0458] is based on the inverse Compton scattering of external ultra-violet/X-ray radiation by a cold ultra-relativistic electron-positron pulsar wind. We show that the extra broad ~30 MeV component should arise from Comptonization of cosmic microwave background radiation. We estimate the main parameters of this component and show that it can be detectable with MeV telescopes such as CGRO/COMPTEL. The location of CGRO/COMPTEL unidentified source GRO J1823-12 close to excess of 105-120 GeV emission (Reg.1 of [1205.4700]) can be interpreted as an argument in favour of astrophysical model of the narrow feature(s) at ~100 GeV.
MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes (IACT) that observe Very High Energy (VHE) gamma ray sources. The PMTs in their cameras are designed to operate under moonlight, but they are limited to Moon phases below 93% (300 Moon hours per year), as they can get damaged if the amount of light they receive is too high. As a result, they cannot be used in the three to five nights around full Moon. We have selected commercial inexpensive UV-pass filters rejecting light above a wavelength of 420 nm, where the moonlight intensity is stronger. We mounted them on light-weight frames that can be easily installed on the telescope cameras. Test observations have been performed during the last nine months, from which a moonlight transmission of about 20% and a Cherenkov light transmission of about 45% are estimated. This allows the observation of sources down to an angular distance of 5 degrees to the Moon during Full Moon: essentially in the whole sky and all possible moonlight conditions. Therefore, the duty cycle of MAGIC can be extended by about 30%, including nights when VHE observations with IACTs are currently not feasible. Here we evaluate the preliminary performance, in terms of sensitivity and energy threshold, of the MAGIC telescopes equipped with the UV-pass filters under different moonlight intensities, as inferred from Crab Nebula observations and Monte Carlo simulations.
Gas inflow feeds galaxies with low metallicity gas from the cosmic web, sustaining star formation across the Hubble time. We make a connection between these inflows and metallicity inhomogeneities in star-forming galaxies, by using synthetic narrow-band images of the Halpha emission line from zoom-in AMR cosmological simulations of galaxies with stellar masses of $M \simeq 10^9 $Msun at redshifts z=2-7. In $\sim$50\% of the cases at redshifts lower than 4, the gas inflow gives rise to star-forming, Halpha-bright, off-centre clumps. Most of these clumps have gas metallicities, weighted by Halpha luminosity, lower than the metallicity in the surrounding interstellar medium by $\sim$0.3 dex, consistent with observations of chemical inhomogeneities at high and low redshifts. Due to metal mixing by shear and turbulence, these metallicity drops are dissolved in a few disc dynamical times. Therefore, they can be considered as evidence for rapid gas accretion coming from cosmological inflow of pristine gas.
We report for the first time a gamma-ray and multiwavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E>100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/- 0.08 year-period gamma-ray cycle, seen in 3.5 oscillation maxima observed, is corroborated by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in sim 10 years of data has a similar period, while the 15 GHz oscillation is less regular. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.
We propose an observational test for gravitationally recoiling supermassive black holes (BHs) in active galactic nuclei, based on a correlation between the velocities of BHs relative to their host galaxies, |\Delta v|, and their obscuring dust column densities, \Sigma_{dust} (both measured along the line of sight). Proxies for both quantities can be derived from spectral features of individual quasars. We use toy models for the distribution of recoil velocities, BH trajectories, and the geometry of obscuring dust tori in galactic centres, to simulate 2.5x10^5 random observations of recoiling quasars. BHs with recoil velocities comparable to the escape velocity from the galactic centre remain bound to the nucleus, and do not fully settle back to the centre of the torus due to dynamical friction in a typical quasar lifetime. We find that |\Delta v| and \Sigma_ {dust} for these BHs are positively correlated. For obscured (\Sigma_{dust}>0) and for partially obscured (0<\Sigma_{dust}<~2.3 g/m^2) quasars with |\Delta v|>=45 km/s, the sample correlation coefficient between log10(|\Delta v|) and \Sigma_{dust} is r_{45} = 0.28+/-0.02 and r_{45} = 0.13+/-0.02, respectively. Allowing for random +/-100 km/s errors in |\Delta v| unrelated to the recoil dilutes the correlation for the partially obscured quasars to r_{45} = 0.026+/-0.004 measured between |\Delta v| and \Sigma_{dust}. A random sample of >~3,500 obscured quasars with |\Delta v|>=45 km/s would allow rejection of the no-correlation hypothesis with 3 sigma significance 95% of the time. Finally, we find that the fraction of obscured quasars, F_{obs}(|\Delta v|), decreases with |\Delta v| from F_{obs}(<10 km/s)>~0.8 to F_{obs}(>10^3 km/s)<~0.4. This predicted trend can be compared to the observed fraction of type II quasars, and can further test combinations of recoil, trajectory, and dust torus models.
The axions produced during the QCD phase transition by vacuum realignment, string decay and domain wall decay thermalize as a result of their gravitational self-interactions when the photon temperature is approximately 500 eV. They then form a Bose-Einstein condensate (BEC). Because the axion BEC rethermalizes on time scales shorter than the age of the universe, it has properties that distinguish it from other forms of cold dark matter. The observational evidence for caustic rings of dark matter in galactic halos is explained if the dark matter is axions, at least in part, but not if the dark matter is entirely WIMPs or sterile neutrinos.
We report results of a spectropolarimetric and photometric monitoring of the
weak-line T Tauri stars (wTTSs) V819 Tau and V830 Tau within the MaTYSSE
programme, involving the ESPaDOnS spectropolarimeter at the
Canada-France-Hawaii Telescope. At ~3 Myr, both stars dissipated their discs
recently and are interesting objects for probing star and planet formation.
Profile distortions and Zeeman signatures are detected in the unpolarized and
circularly-polarized lines, whose rotational modulation we modelled using
tomographic imaging, yielding brightness and magnetic maps for both stars.
We find that the large-scale magnetic fields of V819 Tau and V830 Tau are
mostly poloidal and can be approximated at large radii by 350-400 G dipoles
tilted at ~30 degrees to the rotation axis. They are significantly weaker than
the field of GQ Lup, an accreting classical T Tauri star (cTTS) with similar
mass and age which can be used to compare the magnetic properties of wTTSs and
cTTSs. The reconstructed brightness maps of both stars include cool spots and
warm plages. Surface differential rotation is small, typically ~4.4x smaller
than on the Sun, in agreement with previous results on wTTSs.
Using our Doppler images to model the activity jitter and filter it out from
the radial velocity (RV) curves, we obtain RV residuals with dispersions of
0.033 and 0.104 km/s for V819 Tau and V830 Tau respectively. RV residuals
suggest that a hot Jupiter may be orbiting V830 Tau, though additional data are
needed to confirm this preliminary result. We find no evidence for close-in
giant planet around V819 Tau.
We explore the relationship between the nonlinear matter power spectrum and the various Lagrangian and Standard Perturbation Theories (LPT and SPT). We first look at it in the context of one dimensional (1-d) dynamics, where 1LPT is exact at the perturbative level and one can exactly resum the SPT series into the 1LPT power spectrum. Shell crossings lead to non-perturbative effects, and the PT ignorance can be quantified in terms of their ratio, which is also the transfer function squared in the absence of stochasticity. At the order of PT we work, this parametrization is equivalent to the results of effective field theory (EFT), and can thus be expanded in terms of the same parameters. We find that its radius of convergence is larger than the SPT loop expansion. The same EFT parametrization applies to all SPT loop terms and, if stochasticity can be ignored, to all N-point correlators. In 3-d, the LPT structure is considerably more complicated, and we find that LPT models with parametrization motivated by the EFT exhibit running with $k$ and that SPT is generally a better choice. Since these transfer function expansions contain free parameters that change with cosmological model their usefulness for broadband power is unclear. For this reason we test the predictions of these models on baryonic acoustic oscillations (BAO) and other primordial oscillations, including string monodromy models, for which we ran a series of simulations with and without oscillations. Most models are successful in predicting oscillations beyond their corresponding PT versions, confirming the basic validity of the model.
We present the results of high resolution co-temporal and co-spatial photospheric and chromospheric observations of sunspot penumbral intrusions. The data was taken with the Swedish Solar Telescope (SST) on the Canary Islands. Time series of Ca\,II H images show a series of transient jets extending roughly 3000 km above a penumbral intrusion into the umbra. For most of the time series jets were seen along the whole length of the intruding bright filament. Some of these jets develop a clear $\lambda$-shaped structure, with a small loop appearing at their footpoint and lasting for around a minute. In the framework of earlier studies, the observed transient $\lambda$ shape of these jets strongly suggests that they are caused by magnetic reconnection between a curved arcade-like or flux-rope like field in the lower part of the penumbral intrusion and the more vertical umbral magnetic field forming a cusp-shaped structure above the penumbral intrusion.
We present several arguments which favor the scenario of two coexisting families of compact stars: hadronic stars and quark stars. Besides the well known hyperon puzzle of the physics of compact stars, a similar puzzle exists also when considering delta resonances. We show that these particles appear at densities close to twice saturation density and must be therefore included in the calculations of the hadronic equation of state. Such an early appearance is strictly related to the value of the L parameter of the symmetry energy that has been found, in recent phenomenological studies, to lie in the range $40<L<62$ MeV. We discuss also the threshold for the formation of deltas and hyperons for hot and lepton rich hadronic matter. Similarly to the case of hyperons, also delta resonances cause a softening of the equation of state which makes it difficult to obtain massive hadronic stars. Quark stars, on the other hand, can reach masses up to $2.75 M_{\odot}$ as predicted by perturbative QCD calculations. We then discuss the observational constraints on the masses and the radii of compact stars. The tension between the precise measurements of high masses and the indications of the existence of very compact stellar objects (with radii of the order of $10$ km) is relieved when assuming that very massive compact stars are quark stars and very compact stars are hadronic stars. Finally, we discuss recent interesting measurements of the eccentricities of the orbits of millisecond pulsars in low mass X-ray binaries. The high values of the eccentricities found in some cases could be explained by assuming that the hadronic star, initially present in the binary system, converts to a quark star due to the increase of its central density.
We will follow the two-families scenario described in the accompanying paper, in which compact stars having a very small radius and masses not exceeding about 1.5$M_\odot$ are made of hadrons, while more massive compact stars are quark stars. In the present paper we discuss the dynamics of the transition of a hadronic star into a quark star. We will show that the transition takes place in two phases: a very rapid one, lasting a few milliseconds, during which the central region of the star converts into quark matter and the process of conversion is accelerated by the existence of strong hydrodynamical instabilities, and a second phase, lasting about ten seconds, during which the process of conversion proceeds till the surface of the star via production and diffusion of strangeness. We will show that these two steps play a crucial role in the phenomenological implications of the model. We will discuss the possible implications of this scenario both for long and for short Gamma Ray Bursts, using the proto-magnetar model as the reference frame of our discussion. We will show that the process of quark deconfinement can be connected to specific observed features of the GRBs. In the case of long GRBs we will discuss the possibility that quark deconfinement is at the origin of the second peak present in quite a large fraction of bursts. Also we will discuss the possibility that long GRBs can take place in binary systems without being associated with a SN explosion. Concerning short GRBs, quark deconfinement can play the crucial role in limiting their duration. Finally we will shortly revisit the possible relevance of quark deconfinement in some specific type of Supernova explosions, in particular in the case of very massive progenitors.
A number of Galactic sources emit GeV-TeV gamma rays that are produced through leptonic and/or hadronic mechanisms. Spectral analysis in this energy range is crucial in order to understand the emission mechanisms. The HAWC Gamma-Ray Observatory, with a large field of view and location at $19^\circ$ N latitude, is surveying the Galactic Plane from high Galactic longitudes down to near the Galactic Center. Data taken with partially-constructed HAWC array in 2013-2014 exhibit TeV gamma-ray emission along the Galactic Plane. A high-level analysis likelihood framework for HAWC, also presented at this meeting, has been developed concurrently with the Multi-Mission Maximum Likelihood (3ML) architecture to deconvolve the Galactic sources and to perform multi-instrument analysis. It has been tested on early HAWC data and the same method will be applied on HAWC data with the full array. I will present preliminary results on Galactic sources from TeV observations with HAWC and from joint analysis on Fermi and HAWC data in GeV-TeV energy range.
We consider the relation between exact solutions of cosmological models having minimally and non-minimally coupled scalar fields. This is done for a particular class of solvable models which, in the Einstein frame, have potentials depending on hyperbolic functions and in the Jordan frame, where the non-minimal coupling is conformal, possess a relatively simple dynamics. We show that a particular model in this class can be generalized to the cases of closed and open Friedmann universes and still exhibits a simple dynamics. Further we illustrate the conditions for the existences of bounces in some sub-classes of the set of integrable models we have considered.
In this work we combine information from relic abundance, direct detection, cosmic microwave background, positron fraction, gamma rays, and colliders to explore the existing constraints on effective couplings between Dark Matter and Standard Model constituents when no underlying model or correlation is assumed. Our results show that Dark Matter masses below 20 GeV are disfavoured at the $3 \sigma$ level by tension between the relic abundance requirement and upper constraints on the Dark Matter couplings. Furthermore, large couplings are typically only allowed in combinations which avoid effective couplings to the nuclei used in direct detection experiments.
Dark matter in the sub-GeV mass range is a theoretically motivated but largely unexplored paradigm. Such light masses are out of reach for conventional nuclear recoil direct detection experiments, but may be detected through the small ionization signals caused by dark matter-electron scattering. Semiconductors are well-studied and are particularly promising target materials because their ${\cal O}(1~\rm{eV})$ band gaps allow for ionization signals from dark matter as light as a few hundred keV. Current direct detection technologies are being adapted for dark matter-electron scattering. In this paper, we provide the theoretical calculations for dark matter-electron scattering rate in semiconductors, overcoming several complications that stem from the many-body nature of the problem. We use density functional theory to numerically calculate the rates for dark matter-electron scattering in silicon and germanium, and estimate the sensitivity for upcoming experiments such as DAMIC and SuperCDMS. We find that the reach for these upcoming experiments has the potential to be orders of magnitude beyond current direct detection constraints and that sub-GeV dark matter has a sizable modulation signal. We also give the first direct detection limits on sub-GeV dark matter from its scattering off electrons in a semiconductor target (silicon) based on published results from DAMIC. We make available publicly our code, QEdark, with which we calculate our results. Our results can be used by experimental collaborations to calculate their own sensitivities based on their specific setup. The searches we propose will probe vast new regions of unexplored dark matter model and parameter space.
We discuss with a rather critical eye the current situation of black hole (BH) solutions in $f(R)$ gravity and shed light about its geometrical and physical significance. We also argue about the meaning, existence or lack thereof of a Birkhoff's theorem in this kind of modified gravity. We focus then on the analysis and quest of $non-trivial$ (i.e. hairy) $asymptotically\,\,flat$ (AF) BH solutions in static and spherically symmetric (SSS) spacetimes in vacuum having the property that the Ricci scalar does $not$ vanish identically in the domain of outer communication. To do so, we provide and enforce the $regularity\,\,conditions$ at the horizon in order to prevent the presence of singular solutions there. Specifically, we consider several classes of $f(R)$ models like those proposed recently for explaining the accelerated expansion in the universe and which have been thoroughly tested in several physical scenarios. Finally, we report analytical and numerical evidence about the $absence$ of $geometric\,\,hair$ in AFSSSBH solutions in those $f(R)$ models. First, we submit the models to the available no-hair theorems, and in the cases where the theorems apply, the absence of hair is demonstrated analytically. In the cases where the theorems do not apply, we resort to a numerical analysis due to the complexity of the non-linear differential equations. Within that aim, a code to solve the equations numerically was built and tested using well know exact solutions. In a future investigation we plan to analyze the problem of hair in De Sitter and Anti-De Sitter backgrounds.
Light and spectroscopy are among the most important and frequently taught topics in introductory, college-level, general education astronomy courses. This is due to the fact that the vast majority of observational data studied by astronomers arrives at Earth in the form of light. While there are many processes by which matter can emit and absorb light, Astro 101 courses typically limit their instruction to the Bohr model of the atom and electron energy level transitions. In this paper, we report on the development of a new Lecture-Tutorial to help students learn about other processes that are responsible for the emission and absorption of light, namely molecular rotations, molecular vibrations, and the acceleration of charged particles by magnetic fields.
We propose a possibility that the inflaton potential is significantly modified after inflation due to heavy field dynamics. During inflation there may be a heavy scalar field stabilized at a value deviated from the low-energy minimum. As the heavy field moves to the low-energy minimum, the inflaton potential could be significantly modified. In extreme cases, the inflaton potential vanishes and the inflaton becomes almost massless at some time after inflation. Such transition of the inflaton potential has interesting implications for primordial density perturbations, reheating, creation of unwanted relics, dark radiation, and experimental search for light degrees of freedom. To be concrete, we consider a chaotic inflation in supergravity where the inflaton mass parameter is promoted to a modulus field, finding that the inflaton becomes stable after the transition and contributes to dark matter. Another example is the new inflation by the MSSM Higgs field which acquires a large expectation value just after inflation, but it returns to the origin after the transition and settles down at the electroweak vacuum. Interestingly, the smallness of the electroweak scale compared to the Planck scale is directly related to the flatness of the inflaton potential.
The mininal dark matter model, in which a real scalar singlet is added to the standard model, has been comprehensively revisited, by taking into account the constraints from perturbativity, electroweak vacuum stability in the early Universe, dark matter direct detection, and Higgs invisible decay at the LHC. We show that the {\it resonant mass region} is totally excluded and the {\it high mass region} is reduced to a narrow window $1.1$ ~TeV $\leq m_{s} \leq$ $ 2.55$~ TeV, which is slightly reduced to $1.1$~TeV $\leq m_{s} \leq$ $ 2.0$~ TeV if the perturbativity is further imposed. This {\it high mass range} can be fully detected by the Xenon 1T experiment.
Quantum gravity, as a fundamental theory of space-time, is expected to reveal how the universe may have started, perhaps during or before an inflationary epoch. It may then leave a potentially observable (but probably minuscule) trace in cosmic large-scale structures that seem to match well with predictions of inflation models. A systematic quest to derive such tiny effects using one approach, loop quantum gravity, has, however, led to unexpected obstacles. Such models remain incomplete, and it is not clear whether loop quantum gravity can be consistent as a full theory. But some surprising effects appear to be generic and would drastically alter our understanding of space-time at large density. These new high-curvature phenomena are a consequence of a widening gap between quantum gravity and ordinary quantum-field theory on a background.
In this paper we present the techniques for computing cosmological bounces in polynomial $f(R)$ theories, whose order is more than two, for spatially flat FLRW spacetime. In these cases the conformally connected Einstein frame shows up multiple scalar potentials predicting various possibilities of cosmological evolution in the physical Jordan frame where the $f(R)$ theory lives. We present a reasonable way in which one can associate the various possible potentials in the Einstein frame, for cubic $f(R)$ gravity, to the cosmological development in the Jordan frame. The issue concerning the energy conditions in $f(R)$ theories is presented. We also point out the very important relationships between the conformal transformations connecting the Jordan frame and the Einstein frame and the various instabilities of $f(R)$ theory. All the calculations are done for cubic $f(R)$ gravity but we hope the results are sufficiently general for higher order polynomial gravity.
We investigate perturbations of a class of spherically symmetric solutions in massive gravity and bi-gravity. The background equations of motion for the particular class of solutions we are interested in reduce to a set of the Einstein equations with a cosmological constant. Thus, the solutions in this class include all the spherically symmetric solutions in general relativity, such as the Friedmann-Lema\^{i}tre-Robertson-Walker solution and the Schwarzschild (-de Sitter) solution, though the one-parameter family of two parameters of the theory admits such a class of solutions. We find that the equations of motion for the perturbations of this class of solutions also reduce to the perturbed Einstein equations at first and second order. Therefore, the stability of the solutions coincides with that of the corresponding solutions in general relativity. In particular, these solutions do not suffer from non-linear instabilities which often appear in the other cosmological solutions in massive gravity and bi-gravity.
The equations of anomalous magnetohydrodynamics describe an Abelian plasma where conduction and chiral currents are simultaneously present and constrained by the second law of thermodynamics. At high frequencies the magnetic currents play the leading role and the spectrum is dominated by two-fluid effects. The system behaves instead as a single fluid in the low-frequency regime where the vortical currents induce potentially large hypermagnetic fields. After deriving the physical solutions of the generalized Appleton-Hartree equation, the corresponding dispersion relations are scrutinized and compared with the results valid for cold plasmas. Hypermagnetic knots and fluid vortices can be concurrently present at very low frequencies and suggest a qualitatively different dynamics of the hydromagnetic nonlinearities.
At present, neutron star radii from both observations and model predictions remain very uncertain. Whereas different models can predict a wide range of neutron star radii, it is not possible for most models to predict radii that are smaller than about 10 km, thus if such small radii are established in the future they will be very difficult to reconcile with model estimates. By invoking a new term in the equation of state that enhances the energy density, but leaves the pressure unchanged we simulate the current uncertainty in the neutron star radii. This new term can be possibly due to the exchange of the weakly interacting light U-boson with appropriate in-medium parameters, which does not compromise the success of the conventional nuclear models. The validity of this new scheme will be tested eventually by more precise measurements of neutron star radii.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
We compute robust lower limits on the spin temperature, $T_{\rm S}$, of the $z=8.4$ intergalactic medium (IGM), implied by the upper limits on the 21-cm power spectrum recently measured by PAPER-64. Unlike previous studies which used a single epoch of reionization (EoR) model, our approach samples a large parameter space of EoR models: the dominant uncertainty when estimating constraints on $T_{\rm S}$. Allowing $T_{\rm S}$ to be a free parameter and marginalising over EoR parameters in our Markov Chain Monte Carlo code 21CMMC, we infer $T_{\rm S}\ge3 {\rm K}$ for a mean IGM neutral fraction of $\bar{x}_{H{\scriptsize I}}\gtrsim0.1$. We further improve on these limits by folding-in additional EoR constraints based on: (i) the dark fraction in QSO spectra, which implies a strict upper limit of $\bar{x}_{H{\scriptsize I}}[z=5.9]\leq 0.06+0.05 \,(1\sigma)$; and (ii) the electron scattering optical depth, $\tau_{e}=0.066\pm0.016\,(1\sigma)$ measured by the Planck satellite. By restricting the allowed EoR models, these additional observations tighten the lower limits on the spin temperature to $T_{\rm S} \ge 6$ K. Thus, even such preliminary 21-cm observations begin to rule out extreme scenarios such as `cold reionization', implying at least some prior heating of the IGM. The analysis framework developed here can be applied to upcoming 21-cm observations, thereby providing unique insights into the sources which heated and subsequently reionized the very early Universe.
We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.
This is the fourth in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here, we use measurements of weak gravitational lensing from the Weighing the Giants project to calibrate Chandra X-ray measurements of total mass that rely on the assumption of hydrostatic equilibrium. This comparison of X-ray and lensing masses provides a measurement of the combined bias of X-ray hydrostatic masses due to both astrophysical and instrumental sources. Assuming a fixed cosmology, and within a characteristic radius (r_2500) determined from the X-ray data, we measure a lensing to X-ray mass ratio of 0.96 +/- 9% (stat) +/- 9% (sys). We find no significant trends of this ratio with mass, redshift or the morphological indicators used to select the sample. In accordance with predictions from hydro simulations for the most massive, relaxed clusters, our results disfavor strong, tens-of-percent departures from hydrostatic equilibrium at these radii. In addition, we find a mean concentration of the sample measured from lensing data of c_200 = $3.0_{-1.8}^{+4.4}$. Anticipated short-term improvements in lensing systematics, and a modest expansion of the relaxed lensing sample, can easily increase the measurement precision by 30--50%, leading to similar improvements in cosmological constraints that employ X-ray hydrostatic mass estimates, such as on Omega_m from the cluster gas mass fraction.
We present an end-to-end, two-phase model for the origin of globular clusters (GCs). In the model, populations of stellar clusters form in the high-pressure discs of high-redshift ($z>2$) galaxies (a rapid-disruption phase due to tidal perturbations from the dense interstellar medium), after which the galaxy mergers associated with hierarchical galaxy formation redistribute the surviving, massive clusters into the galaxy haloes, where they remain until the present day (a slow-disruption phase due to tidal evaporation). The high galaxy merger rates of $z>2$ galaxies allow these clusters to be `liberated' into the galaxy haloes before they are disrupted within the high-density discs. This physically-motivated toy model is the first to include the rapid-disruption phase, which is shown to be essential for simultaneously reproducing the wide variety of properties of observed GC systems, such as their universal characteristic mass-scale, the dependence of the specific frequency on metallicity and galaxy mass, the GC system mass-halo mass relation, the constant number of GCs per unit supermassive black hole mass, and the colour bimodality of GC systems. The model predicts that most of these observables were already in place at $z=1$-$2$, although under rare circumstances GCs may still form in present-day galaxies. In addition, the model provides important constraints on models for multiple stellar populations in GCs by putting limits on initial GC masses and the amount of pristine gas accretion. The paper is concluded with a discussion of these and several other predictions and implications, as well as the main open questions in the field.
We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.
We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed by the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White profile between $r\approx1.5~\rm{kpc}$ and $r\approx10~\rm{kpc}$, and a flatter slope at smaller radii. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.96$ within $r=500~\rm{pc}$) and there is no offset larger than $45~\rm{pc}$ between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.
Recent Spitzer/InfraRed Array Camera (IRAC) photometric observations have revealed that rest-frame optical emission lines contribute signficantly to the broadband fluxes of high-redshift galaxies. Specifically, in the narrow redshift range z~5.1-5.4 the [3.6]-[4.5] color is expected to be very red, due to contamination of the 4.5-micron band by the dominant Halpha line, while the 3.6-micron filter is free of nebular emission lines. We take advantage of new reductions of deep Spitzer/IRAC imaging over the GOODS-North+South fields (Labbe+2015) to obtain a clean measurement of the mean Halpha equivalent width from the [3.6]-[4.5] color in the redshift range z=5.1-5.4. The selected sources either have measured spectroscopic redshifts (13 sources) or lie very confidently in the redshift range z=5.1-5.4 based on the photometric redshift likelihood intervals (11 sources). Our z_{phot}=5.1-5.4 sample and z_{spec}=5.10-5.40 spectroscopic sample have a mean [3.6]-[4.5] color of 0.31+/-0.05 mag and 0.35+/-0.07 mag, implying a rest-frame equivalent width EW(Halpha+[NII]+[SII]) of 665+/-53 Angstroms and 707+/-74 Angstroms, respectively, for sources in these samples. These values are consistent albeit slightly higher than derived by Stark+2013 at z~4, suggesting an evolution to higher values of the Halpha+[NII]+[SII] EW at z>2. Using the 3.6micron band, which is free of emission line contamination, we perform robust SED fitting and find a median specific star formation rate of sSFR = 17_{-5}^{+2} Gyr^{-1}, 7_{-2}^{+1}x higher than at z~2. We find no strong correlation (<2sigma) between the Halpha+[NII]+[SII] EW and the stellar mass of sources. Before the advent of JWST, improvements in these results will come through an expansion of current spectroscopic samples and deeper Spitzer/IRAC measurements.
Investigating the molecular gas in the inner regions of protoplanetary disks provides insight into how the molecular disk environment changes during the transition from primordial to debris disk systems. We conduct a small survey of molecular hydrogen (H$_2$) fluorescent emission, using 14 well-studied Classical T Tauri stars at two distinct dust disk evolutionary stages, to explore how the structure of the inner molecular disk changes as the optically thick warm dust dissipates. We simulate the observed HI-Lyman $\alpha$-pumped H$_2$ disk fluorescence by creating a 2D radiative transfer model that describes the radial distributions of H$_{2}$ emission in the disk atmosphere and compare these to observations from the Hubble Space Telescope. We find the radial distributions that best describe the observed H$_2$ FUV emission arising in primordial disk targets (full dust disk) are demonstrably different than those of transition disks (little-to-no warm dust observed). For each best-fit model, we estimate inner and outer disk emission boundaries (r$_{in}$ and r$_{out}$), describing where the bulk of the observed H$_2$ emission arises in each disk, and we examine correlations between these and several observational disk evolution indicators, such as n$_{13-31}$, r$_{in,CO}$, and the mass accretion rate. We find strong, positive correlations between the H$_2$ radial distributions and the slope of the dust SED, implying the behavior of the molecular disk atmosphere changes as the inner dust clears in evolving protoplanetary disks. Overall, we find that H$_2$ inner radii are $\sim$4 times larger in transition systems, while the bulk of the H$_2$ emission originates inside the dust gap radius for all transitional sources.
Recent advances on gamma-ray observations from SuperNova Remnants and Molecular Clouds offer the possibility to study in detail the properties of the propagation of escaping Cosmic Rays (CR). However, a complete theory for CR transport outside the acceleration site has not been developed yet. Two physical processes are thought to be relevant to regulate the transport: the growth of waves caused by streaming instability, and possible wave damping mechanisms that reduce the growth of the turbulence. Only a few attempts have been made so far to incorporate these mechanisms in the theory of CR diffusion. In this work we present recent advances in this subject. In particular, we show results obtained by solving the coupled equations for the diffusion of CRs and the evolution of Alfven waves. We discuss the importance of streaming instabilities and wave damping in different ISM phases.
N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped off individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at this http URL
The T Tauri star PTFO 8-8695 exhibits periodic fading events that have been interpreted as the transits of a giant planet on a precessing orbit. Here we present three tests of the planet hypothesis. First, we sought evidence for the secular changes in light-curve morphology that are predicted to be a consequence of orbital precession. We observed 28 fading events spread over several years, and did not see the expected changes. Instead we found that the fading events are not strictly periodic. Second, we attempted to detect the planet's radiation, based on infrared observations spanning the predicted times of occultations. We ruled out a signal of the expected amplitude. Third, we attempted to detect the Rossiter-McLaughlin effect by performing high-resolution spectroscopy throughout a fading event. No effect was seen at the expected level, ruling out most (but not all) possible orientations for the hypothetical planetary orbit. Our spectroscopy also revealed strong, time-variable, high-velocity H{\alpha} and Ca H & K emission features. All these observations cast doubt on the planetary hypothesis, and suggest instead that the fading events represent starspots, eclipses by circumstellar dust, or occultations of an accretion hotspot.
We analysed 385 galactic spectra from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) that belong to the catalog of isolated pairs of galaxies by Karachentsev. The spectra corresponds to physical pairs of galaxies as defined by V $\leq$ 1200 Km/s and a pair separation $\leq$ 100 kpc. We search for the incidence of nuclear activity, both thermal (star-forming) and non-thermal -Active Galactic Nuclei (AGN). After a careful extraction of the nuclear spectra, we use diagnostic diagrams and find that the incidence of AGN activity is 48 \% in the paired galaxies with emission lines and 40\% for the total sample (as compared to $\sim$ 43 \% and 41\% respectively in a sample of isolated galaxies). These results remain after dissecting the effects of morphological type and galactic stellar mass (with only a small, non significant, enhancement of the AGN fraction in pairs of objects). These results suggest that weak interactions are not necessary or sufficient to trigger low-luminosity AGN. Since the fraction of AGN is predominant in early type spiral galaxies, we conclude that the role of a bulge, and a large gas reservoir are both essential for the triggering of nuclear activity. The most striking result is that type 1 galaxies are almost absent from the AGN sample. This result is in conflict with the Unified Model, and suggests that high accretion rates are essential to form the Broad Line Region in active galaxies.
The Cygnus region hosts one of the most remarkable star-forming regions in the Milky Way. Indeed, the total mass in molecular gas of the Cygnus X complex exceeds 10 times the total mass of all other nearby star-forming regions. Surveys at all wavelengths, from radio to gamma-rays, reveal that Cygnus contains such a wealth and variety of sources---supernova remnants (SNRs), pulsars, pulsar wind nebulae (PWNe), H II regions, Wolf-Rayet binaries, OB associations, microquasars, dense molecular clouds and superbubbles---as to practically be a galaxy in microcosm. The gamma-ray observations along reveal a wealth of intriguing sources at energies between 1 GeV and tens of TeV. However, a complete understanding of the physical phenomena producing this gamma-ray emission first requires us to disentangle overlapping sources and reconcile discordant pictures at different energies. This task is made more challenging by the limited angular resolution of instruments such as the Fermi Large Area Telescope, ARGO-YBJ, and HAWC and the limited sensitivity and field of view of current imaging atmospheric Cherenkov telescopes (IACTs). The Cherenkov Telescope Array (CTA), with its improved angular resolution, large field of view, and order of magnitude gain in sensitivity over current IACTs, has the potential to finally create a coherent and well-resolved picture of the Cygnus region between a few tens of GeV and a hundred TeV. We describe a proposed strategy to study the Cygnus region using CTA data, which combines a survey of the whole region at $65^{\circ} < l < 85^{\circ}$ and $-3.5^{\circ} < b < 3.5^{\circ} $ with deeper observations of two sub-regions that host rich groups of known gamma-ray sources.
We explore the phenomenological consequences of general late-time modifications of gravity in the quasi-static approximation, in the case where baryons and cold dark matter have distinct couplings to the gravitational sector. Assuming spectroscopic and photometric surveys with configuration parameters similar to those of the Euclid mission, we derive constraints on our effective description from three observables: the galaxy power spectrum in redshift space, tomographic weak-lensing shear power spectrum and the correlation spectrum between the integrated Sachs-Wolfe effect and the galaxy distribution. In particular, with $\Lambda$CDM as fiducial model and a specific choice for the time dependence of our effective functions, we perform a Fisher matrix analysis and find that the unmarginalized $68\%$ CL errors on the three parameters describing the modifications of gravity are of order $\sigma\sim10^{-3}$. We also consider two other fiducial models. A nonminimal coupling of CDM enhances the effects of modified gravity and reduces the above statistical errors accordingly. In all cases, we find that the parameters are highly degenerate, which prevents the inversion of the Fisher matrices. Although all three observational probes are complementary in breaking some of the degeneracies, the ISW-galaxy correlation stands out as a promising probe to constrain the modifications of gravity.
The energy dissipation mechanism within Gamma ray bursts' (GRBs) ultra-relativistic outflows that drive the prompt $\gamma$-ray emission remains uncertain. Two leading candidates are internal shocks and magnetic reconnection. While the emission from internal shocks has been extensively studied, that from reconnection still has few quantitative predictions. We study the prompt GRB emission from magnetic reconnection and compare its expected temporal and spectral properties to observations. The main difference from internal shocks is that for magnetic reconnection one expects relativistic bulk motions with a Lorentz factor of $\Gamma'\gtrsim$a few in the mean rest frame of the outflow - the comoving frame. We consider a thin spherical shell (or reconnection layer) expanding at a bulk Lorentz factor $\Gamma\gg 1$ in which the emitting material moves with $\Gamma'$ in the comoving frame along this layer in two anti-parallel directions (e.g. of the reconnecting field lines). The resulting relativistic beaming of the emitted radiation in the comoving frame affects the temporal and spectral properties of the GRB lightcurves. In particular, if the emission occurs in radii $R_0<R<R_0+\Delta R$ then (for a constant $\Gamma$) the observed pulse width is $\Delta T \sim T_0\max(1/\Gamma',\,\Delta R/R_0)$ where $T_0 = R_0/2c\Gamma^2$, i.e. it can be up to $\sim\Gamma'$ times shorter than for isotropic emission in the comoving frame. Such anisotropic emission from magnetic reconnection can potentially reproduce many of the observed prompt GRB properties: its variability, pulse asymmetry and the rapid decay phase at its end. It can also account for many of the observed correlations: luminosity-variability, peak luminosity$\,$--$\,$peak frequency, pulse width energy dependence/spectral lags, and both hard to soft evolution (for $\Gamma'\lesssim 2$) or intensity tracking (for $\Gamma'>2$).
We explore the systematic uncertainties in dark energy constraints using the latest observational data from Type Ia Supernovae (SNe Ia), galaxy clustering, and cosmic microwave background anisotropy (CMB) data. We use the Joint Lightcurve Analysis (JLA) set of 740 SNe Ia, galaxy clustering measurements of H(z)s and D_A(z)/s (where s is the sound horizon at the drag epoch) from the Sloan Digital Sky Survey (SDSS) at z=0.35 (SDSS DR7) and z=0.57 (BOSS DR11), and the distance priors that we have derived from the 2015 Planck data (we present the mean values and covariance matrices required for using these). We find that omitting the BOSS DR11 measurement of H(z)s at z=0.57 leads to more concordant cosmological constraints, indicative of possible systematic uncertainties that affect the measurement of the line-of-sight galaxy clustering. We also find that flux-averaging of SNe Ia at z>= 0.5 gives significantly tighter constraints on dark energy; this can be due to the reduction in the distance measurement bias from flux averaging SNe Ia. Taking into consideration the possible systematic uncertainties, current observational data continue to be consistent with a flat universe with a cosmological constant.
The process of gathering and associating data from multiple sensors or sub-detectors due to a common physical event (the process of event-building) is used in many fields, including high-energy physics and $\gamma$-ray astronomy. Fault tolerance in event-building is a challenging problem that increases in difficulty with higher data throughput rates and increasing numbers of sub-detectors. We draw on biological self-assembly models in the development of a novel event-building paradigm that treats each packet of data from an individual sensor or sub-detector as if it were a molecule in solution. Just as molecules are capable of forming chemical bonds, "bonds" can be defined between data packets using metadata-based discriminants. A database -- which plays the role of a beaker of solution -- continually selects pairs of assemblies at random to test for bonds, which allows single packets and small assemblies to aggregate into larger assemblies. During this process higher-quality associations supersede spurious ones. The database thereby becomes fluid, dynamic, and self-correcting rather than static. We will describe tests of the self-assembly paradigm using our first fluid database prototype and data from the VERITAS $\gamma$-ray telescope.
We present the discovery by the WASP-South survey of three planets transiting moderately bright stars (V ~ 11). WASP-120b is a massive (5.0MJup) planet in a 3.6-day orbit that we find likely to be eccentric (e = 0.059+0.025-0.018) around an F5 star. WASP-122b is a hot-Jupiter (1.37MJup, 1.79RJup) in a 1.7-day orbit about a G4 star. Our predicted transit depth variation cause by the atmosphere of WASP-122b suggests it is well suited to characterisation. WASP-123b is a hot-Jupiter (0.92MJup, 1.33RJup) in a 3.0-day orbit around an old (~ 7 Gyr) G5 star.
Mass is one of the most fundamental parameters characterizing the dynamics of a coronal mass ejection (CME). It has been found that CME apparent mass measured from the brightness enhancement in coronagraph images shows an increasing trend during its evolution in the corona. However, the physics behind it is not clear. Does the apparent mass gain come from the mass outflow from the dimming regions in the low corona, or from the pileup of the solar wind plasma around the CME when it propagates outwards from the Sun? We analyzed the mass evolution of six CME events. Their mass can increase by a factor of 1.6 to 3.2 from 4 to 15 Rs in the field of view (FOV) of the coronagraph on board the Solar Terrestrial Relations Observatory (STEREO). Over the distance about 7 to 15 Rs, where the coronagraph occulting effect can be negligible, the mass can increase by a factor of 1.3 to 1.7. We adopted the `snow-plough' model to calculate the mass contribution of the piled-up solar wind in the height range from about 7 to 15 Rs. For 2/3 of the events, the solar wind pileup is not sufficient to explain the measured mass increase. In the height range from about 7 to 15 Rs, the ratio of the modeled to the measured mass increase is roughly larger than 0.55. Although the ratios are believed to be overestimated, the result gives evidence that the solar wind pileup probably makes a non-negligible contribution to the mass increase. It is not clear yet whether the solar wind pileup is a major contributor to the final mass derived from coronagraph observations. However, our study suggests that the solar wind pileup plays increasingly important role in the mass increase as a CME moves further away from the Sun.
Future hard (10 -100 keV) X-ray telescopes (SIMBOL-X, Con-X, HEXIT-SAT, XEUS) will implement focusing optics with multilayer coatings: in view of the production of these optics we are exploring several deposition techniques for the reflective coatings. In order to evaluate the achievable optical performance X-Ray Reflectivity (XRR) measurements are performed, which are powerful tools for the in-depth characterization of multilayer properties (roughness, thickness and density distribution). An exact extraction of the stack parameters is however difficult because the XRR scans depend on them in a complex way. The PPM code, developed at ERSF in the past years, is able to derive the layer-by-layer properties of multilayer structures from semi-automatic XRR scan fittings by means of a global minimization procedure in the parameters space. In this work we will present the PPM modeling of some multilayer stacks (Pt/C and Ni/C) deposited by simple e-beam evaporation. Moreover, in order to verify the predictions of PPM, the obtained results are compared with TEM profiles taken on the same set of samples. As we will show, PPM results are in good agreement with the TEM findings. In addition, we show that the accurate fitting returns a physically correct evaluation of the variation of layers thickness through the stack, whereas the thickness trend derived from TEM profiles can be altered by the superposition of roughness profiles in the sample image.
Several celestial bodies in co-orbital configurations exist in the solar system. However, co-orbital exoplanets have not yet been discovered. This lack may result from a degeneracy between the signal induced by co-orbital planets and other orbital configurations. Here we determine a criterion for the detectability of quasi-circular co-orbital planets and develop a demodulation method to bring out their signature from the observational data. We show that the precision required to identify a pair of co-orbital planets depends only on the libration amplitude and on the planet's mass ratio. We apply our method to synthetic radial velocity data, and show that for tadpole orbits we are able to determine the inclination of the system to the line of sight. Our method is also valid for planets detected through the transit and astrometry techniques.
A proposed method for dealing with foreground emission in upcoming 21-cm observations from the epoch of reionization is to limit observations to an uncontaminated window in Fourier space. Foreground emission can be avoided in this way, since it is limited to a wedge-shaped region in $k_{\parallel}, k_{\perp}$ space. However, the power spectrum is anisotropic owing to redshift-space distortions from peculiar velocities. Consequently, the 21-cm power spectrum measured in the foreground avoidance window---which samples only a limited range of angles close to the line-of-sight direction---differs from the full spherically-averaged power spectrum which requires an average over \emph{all} angles. In this paper, we calculate the magnitude of this "wedge bias" for the first time. We find that the bias is strongest at high redshifts, where measurements using foreground avoidance will over-estimate the power spectrum by around 100 per cent, possibly obscuring the distinctive rise and fall signature that is anticipated for the spherically-averaged 21-cm power spectrum. In the later stages of reionization, the bias becomes negative, and smaller in magnitude ($\lesssim 20$ per cent). The effect shows only a weak dependence on spatial scale and reionization topology.
The disk fraction, the percentage of stars with disks in a young cluster, is widely used to investigate the lifetime of the protoplanetary disk, which can impose an important constraint on the planet formation mechanism. The relationship between the decay timescale of the disk fraction and the mass dissipation timescale of an individual disk, however, remains unclear. Here we investigate the effect of the disk mass function (DMF) on the evolution of the disk fraction. We show that the time variation in the disk fraction depends on the spread of the DMF and the detection threshold of the disk. In general, the disk fraction decreases more slowly than the disk mass if a typical initial DMF and a detection threshold are assumed. We find that, if the disk mass decreases exponentially, {the mass dissipation timescale of the disk} can be as short as $1\,{\rm Myr}$ even when the disk fraction decreases with the time constant of ${\sim}2.5\,{\rm Myr}$. The decay timescale of the disk fraction can be an useful parameter to investigate the disk lifetime, but the difference between the mass dissipation of an individual disk and the decrease in the disk fraction should be properly appreciated to estimate the timescale of the disk mass dissipation.
Far-IR 16-1000 $\mu$m spectra of Saturn's hydrogen-helium continuum measured by Cassini's Composite Infrared Spectrometer (CIRS) are inverted to construct a near-continuous record of upper tropospheric (70-700 mbar) temperatures and para-H$_2$ fraction as a function of latitude, pressure and time for a third of a Saturnian year (2004-2014, from northern winter to northern spring). The thermal field reveals evidence of reversing summertime asymmetries superimposed onto the belt/zone structure. The temperature structure that is almost symmetric about the equator by 2014, with seasonal lag times that increase with depth and are qualitatively consistent with radiative climate models. Localised heating of the tropospheric hazes (100-250 mbar) create a distinct perturbation to the temperature profile that shifts in magnitude and location, declining in the autumn hemisphere and growing in the spring. Changes in the para-H$_2$ ($f_p$) distribution are subtle, with a 0.02-0.03 rise over the spring hemisphere (200-500 mbar) perturbed by (i) low-$f_p$ air advected by both the springtime storm of 2010 and equatorial upwelling; and (ii) subsidence of high-$f_p$ air at northern high latitudes, responsible for a developing north-south asymmetry in $f_p$. Conversely, the shifting asymmetry in the para-H$_2$ disequilibrium primarily reflects the changing temperature structure (and the equilibrium distribution of $f_p$), rather than actual changes in $f_p$ induced by chemical conversion or transport. CIRS results interpolated to the same point in the seasonal cycle as re-analysed Voyager-1 observations show qualitative consistency, with the exception of the tropical tropopause near the equatorial zones and belts, where downward propagation of a cool temperature anomaly associated with Saturn's stratospheric oscillation could potentially perturb tropopause temperatures, para-H$_2$ and winds. [ABRIDGED]
By studying young magnetically active late-type stars, i.e. analogues to the young Sun, one can draw conclusions on the evolution of the solar dynamo. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric spots in three young late-type stars. High-resolution spectropolarimetry of the targets were obtained with the HARPSpol instrument mounted at the ESO 3.6 m telescope. The signal-to-noise ratio of the Stokes IV measurements were boosted by combining the signal from a large number of spectroscopic absorption lines through the least squares deconvolution technique. Surface brightness and magnetic field maps were calculated using the Zeeman-Doppler imaging technique. All the three targets show clear signs of both magnetic fields and cool spots. Only one of the targets, namely V1358 Ori, shows evidence of the dominance of non-axisymmetric modes. In two of the targets, the poloidal field is significantly stronger than the toroidal one, indicative of an $\alpha^2$-type of a dynamo, in which convective turbulence effects dominate over the weak differential rotation. In two of the cases there is a slight anti-correlation between the cool spots and the strength of the radial magnetic field. However, even in these cases the correlation is much weaker than in the case of sunspots. The weak correlation between the measured radial magnetic field and cool spots may indicate a more complex magnetic field structure in the spots or spot groups involving mixed magnetic polarities. Comparison with a previously published magnetic field map shows that on one of the stars, HD 29615, the underlying magnetic field has changed its polarity between 2009 and 2013.
The optics of a number of future X-ray telescopes will have very long focal lengths (10 - 20 m), and will consist of a number of nested/stacked thin, grazing-incidence mirrors. The optical quality characterization of a real mirror can be obtained via profile metrology, and the Point Spread Function of the mirror can be derived via one of the standard computation methods. However, in practical cases it can be difficult to access the optical surfaces of densely stacked mirror shells, after they have been assembled, using the widespread metrological tools. For this reason, the assessment of the imaging resolution of a system of mirrors is better obtained via a direct, full-illumination test in X-rays. If the focus cannot be reached, an intra-focus test can be performed, and the image can be compared with the simulation results based on the metrology, if available. However, until today no quantitative information was extracted from a full-illumination, intra-focal exposure. In this work we show that, if the detector is located at an optimal distance from the mirror, the intensity variations of the intra-focal, full-illumination image in single reflection can be used to reconstruct the profile of the mirror surface, without the need of a wavefront sensor. The Point Spread Function can be subsequently computed from the reconstructed mirror shape. We show the application of this method to an intra-focal (8 m distance from mirror) test performed at PANTER on an optical module prototype made of hot-slumped glass foils with a 20 m focal length, from which we could derive an expected imaging quality near 16 arcsec HEW.
We report the discovery of KELT-10b, the first transiting exoplanet discovered using the KELT-South telescope. KELT-10b is a highly inflated sub-Jupiter mass planet transiting a relatively bright $V = 10.7$ star (TYC 8378-64-1), with T$_{eff}$ = $5948\pm74$ K, $\log{g}$ = $4.319_{-0.030}^{+0.020}$ and [Fe/H] = $0.09_{-0.10}^{+0.11}$, an inferred mass M$_{*}$ = $1.112_{-0.061}^{+0.055}$ M$_{\odot}$ and radius R$_{*}$ = $1.209_{-0.035}^{+0.047}$ R$_{\odot}$. The planet has a radius R$_{P}$ = $1.399_{-0.049}^{+0.069}$ R$_{J}$ and mass M$_{P}$ = $0.679_{-0.038}^{+0.039}$ M$_{J}$. The planet has an eccentricity consistent with zero and a semi-major axis $a$ = $0.05250_{-0.00097}^{+0.00086}$ AU. The best fitting linear ephemeris is $T_{0}$ = 2457066.72045$\pm$0.00027 BJD$_{TDB}$ and P = 4.1662739$\pm$0.0000063 days. This planet joins a group of highly inflated transiting exoplanets with a radius much larger and a mass much less than those of Jupiter. The planet, which boasts deep transits of 1.4%, has a relatively high equilibrium temperature of T$_{eq}$ = $1377_{-23}^{+28}$ K, assuming zero albedo and perfect heat redistribution. KELT-10b receives an estimated insolation of $0.817_{-0.054}^{+0.068}$ $\times$ 10$^9$ erg s$^{-1}$ cm$^{-2}$, which places it far above the insolation threshold above which hot Jupiters exhibit increasing amounts of radius inflation. Evolutionary analysis of the host star suggests that KELT-10b is unlikely to survive beyond the current subgiant phase, due to a concomitant in-spiral of the planet over the next $\sim$1 Gyr. The planet transits a relatively bright star which is accessible to large telescopes and exhibits the third largest transit depth of all transiting exoplanets with V $<$ 11 in the southern hemisphere, making it a promising candidate for future atmospheric characterization studies.
The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement of the spatial distribution of density, however, is difficult and indirect seismology inversion methods have been suggested as an alternative. We applied Bayesian inference, model comparison, and model-averaging techniques to the inference of the cross-field density structuring in solar magnetic waveguides using information on periods and damping times for resonantly damped magnetohydrodynamic (MHD) transverse kink oscillations. Three commonly employed alternative profiles were used to model the variation of the mass density across the waveguide boundary. Parameter inference enabled us to obtain information on physical quantities such as the Alfv\'en travel time, the density contrast, and the transverse inhomogeneity length scale. The inference results from alternative density models were compared and their differences quantified. Then, the relative plausibility of the considered models was assessed by performing model comparison. Our results indicate that the evidence in favor of any of the three models is minimal, unless the oscillations are strongly damped. In such a circumstance, the application of model-averaging techniques enables the computation of an evidence-weighted inference that takes into account the plausibility of each model in the calculation of a combined inversion for the unknown physical parameters.
We present the development of a novel 11328 pixel silicon photomultiplier (SiPM) camera for use with a ground-based Cherenkov telescope with Schwarzschild-Couder optics as a possible medium-sized telescope for the Cherenkov Telescope Array (CTA). The finely pixelated camera samples air-shower images with more than twice the optical resolution of cameras that are used in current Cherenkov telescopes. Advantages of the higher resolution will be a better event reconstruction yielding improved background suppression and angular resolution of the reconstructed gamma-ray events, which is crucial in morphology studies of, for example, Galactic particle accelerators and the search for gamma-ray halos around extragalactic sources. Packing such a large number of pixels into an area of only half a square meter and having a fast readout directly attached to the back of the sensors is a challenging task. For the prototype camera development, SiPMs from Hamamatsu with through silicon via (TSV) technology are used. We give a status report of the camera design and highlight a number of technological advancements that made this development possible.
Variations in the photometric parameters of stellar systems as a function of their evolution and the stellar populations comprising them are investigated. A set of seven evolutionary models with an exponential decrease in the star-formation rate and 672 models with a secondary burst of star formation are considered. The occurrence of a secondary burst of star formation can shift the position of a stellar system on two-color diagrams to the right or left of the normal color sequence for galaxies and the extinction line. This makes it possible to estimate the composition of the stellar population of a galaxy with a nonmonotonic star formation history from its position on two-color diagrams. Surface photometry in both the optical (UBVRI) and near-IR (JHK) is used to study the stellar populations and star-formation histories in the structural components (nucleus, bulge, disk, spiral arms, bar, ring) of 26 galaxies of various morphological types (from S0 to Sd). Components (nucleus, bulge, bar) with color characteristics corresponding to stellar systems with secondary bursts of star formation are indicated in 10 of the 26 galaxies. The parameters of these secondary bursts are estimated. Five of the 10 galaxies with complex star-formation histories display clear structural perturbations. Appreciable differences in the photometric characteristics of relatively red early-type galaxies (S0-Sb) and relatively blue later-type galaxies (Sb-Sd) have been found. Galaxies of both early and late types are encountered among the Sb-galaxies. Lenticular galaxies do not display different photometric characteristics from early-type spiral galaxies.
Estimates of the mass distribution and dark-matter (DM) content of dwarf spheroidal galaxies (dSphs) are usually derived under the assumption that the effect of the tidal field of the host galaxy is negligible over the radial extent probed by kinematic data-sets. We assess the implications of this assumption in the specific case of the Fornax dSph by means of N-body simulations of a satellite orbiting around the Milky Way. We consider observationally-motivated orbits and we tailor the initial distributions of the satellite's stars and DM to match, at the end of the simulations, the observed structure and kinematics of Fornax. In all our simulations the present-day observable properties of Fornax are not significantly influenced by tidal effects. The DM component is altered by the interaction with the Galactic field (up to 20% of the DM mass within 1.6 kpc is lost), but the structure and kinematics of the stellar component are only mildly affected even in the more eccentric orbit (more than 99% of the stellar particles remain bound to the dwarf). In the simulations that successfully reproduce Fornax's observables, the dark-to-luminous mass ratio within 1.6 kpc is in the range 5-6, and up to 16-18 if measured within 3 kpc.
Minimization of charged particle background in X-ray telescopes is a well known issue. Charged particles (chiefly protons and electrons) naturally present in the cosmic environment constitute an important background source when they collide with the X-ray detector. Even worse, a serious degradation of spectroscopic performances of the X-ray detector was observed in Chandra and Newton-XMM, caused by soft protons with kinetic energies ranging between 100 keV and some MeV being collected by the grazing-incidence mirrors and funneled to the detector. For a focusing telescope like SIMBOL-X, the exposure of the soft X-ray detector to the proton flux can increase significantly the instrumental background, with a consequent loss of sensitivity. In the worst case, it can also seriously compromise the detector duration. A well-known countermeasure that can be adopted is the implementation of a properly-designed magnetic diverter, that should prevent high-energy particles from reaching the focal plane instruments of SIMBOL-X. Although Newton-XMM and Swift-XRT are equipped with magnetic diverters for electrons, the magnetic fields used are insufficient to effectively act on protons. In this paper, we simulate the behavior of a magnetic diverter for SIMBOL-X, consisting of commercially-available permanent magnets. The effects of SIMBOL-X optics is simulated through GEANT4 libraries, whereas the effect of the intense required magnetic fields is simulated along with specifically-written numerical codes in IDL.
If dark matter consists of hidden-sector photons which kinetically mix with regular photons, a tiny oscillating electric-field component is present wherever we have dark matter. In the surface of conducting materials this induces a small probability to emit single photons almost perpendicular to the surface, with the corresponding photon frequency matching the mass of the hidden photons. We report on a construction of an experimental setup with a large ~14 m2 spherical metallic mirror that will allow for searches of hidden-photon dark matter in the eV and sub-eV range by application of different electromagnetic radiation detectors. We discuss sensitivity and accessible regions in the dark matter parameter space.
The MAGIC telescopes are an array of two imaging atmospheric Cherenkov
telescopes (IACTs) studying the gamma ray sky at very high-energies (VHE; E>100
GeV). The observations are performed in stereoscopic mode, with both telescopes
pointing at the same position in the sky. The MAGIC field of view (FoV)
acceptance for hadrons and gamma rays has a complex shape, which depends on
several parameters such as the azimuth and zenith angle of the observations. In
the standard MAGIC analysis, the strategy adopted for estimating this
acceptance is not optimal in the case of complex FoVs.
In this contribution we present the results of systematic studies intended to
characterise the acceptance for the entire FoV. These studies open the
possibility to apply improved background estimation methods to the MAGIC data,
useful to investigate the morphology of extended or multiple sources.
We present a theory in 3 parts, of the long-term (or secular) evolution of stellar systems orbiting within the sphere of influence of massive black holes in galactic nuclei. Here we describe the secular collisionless dynamics of a (Keplerian) stellar system of mass $M$ orbiting a black hole of mass $M_\bullet \gg M$. The stellar distribution function (DF) $f$ obeys the collisionless Boltzmann equation (CBE) in 6-dim phase space. The small mass ratio, $\varepsilon = M/M_\bullet \ll 1$, implies a separation of time scales in the motions of stars: the fast Kepler orbital periods and the secular time scale which is longer by a factor $\varepsilon^{-1}$. We orbit-average the CBE over the fast Keplerian orbital phase using the Method of Multiple Scales. Then $f$ is expressed as the sum of a secular DF $F$ in a 5-dim (Gaussian Ring) space, and small fluctuations that remain of $O(\varepsilon)$ over secular times. $F$ obeys a secular CBE that includes stellar self-gravity, general relativistic corrections up to 1.5 post-Newtonian order, and external sources. Secular dynamics conserves the semi-major axis of every star. This additional integral of motion promotes extra regularity of the stellar orbits, and enables the construction of secular equilibrium DFs ($F_0$) through a Secular Jeans theorem. Secular equilibria allow for varied spatial geometries including figure rotation. A linearized secular CBE determines the linear response and stability of $F_0$. Spherical, non-rotating equilibria may support small-amplitude, long-lived, warp-like distortions. We also prove that an axisymmetric, zero-thickness, flat disc is secularly stable to all in-plane perturbations, when its DF $F_0$ is a monotonic function of the angular momentum at fixed energy.
We present a first-principles theory of Resonant Relaxation (RR) of stellar systems orbiting within the sphere of influence of massive black holes in galactic nuclei. We extend the rigorous kinetic theory of Gilbert (1968) to include the Keplerian field of a black hole of mass $M_\bullet$, and specialize to a (Keplerian) stellar system of mass $M \ll M_\bullet$. Using the results of the secular collisionless theory of Paper I, we orbit-average the kinetic equation through perturbative development in the small parameter $\varepsilon = M/M_\bullet$. This is supplemented with contributions from general relativistic corrections up to 1.5 post-Newtonian order and external gravitational sources. The result is a kinetic equation for a secular distribution function (DF) in 5-dim (Gaussian Ring) space, with explicit forms for the fluctuation and dissipation components of the collision integral. For general DFs, both apsidal and nodal precessions contribute to RR; so the traditional, physically-motivated distinction between scalar-RR and vector-RR disappears. Irreversible 2-particle correlations, that build up through secular gravitational interactions, are the driving agents of RR. The correlation function can be written in terms of the wake function, which is the linear response of the system to the perturbation offered by any chosen stellar orbit. The relationship includes direct interactions between pairs of stars as well as collective effects (gravitational polarization). We discuss the interplay of secular dynamics and RR in the evolution toward secular thermodynamic equilibria. In Paper III we apply our RR theory to axisymmetric discs, and provide explicit formulae for the loss cone rates at which mass, energy and angular momentum are fed to the black hole.
Accurate star formation histories (SFHs) of galaxies are fundamental for
understanding the build-up of their stellar content. However, the most accurate
SFHs - those obtained from colour-magnitude diagrams (CMDs) of resolved stars
reaching the oldest main sequence turnoffs (oMSTO) - are presently limited to a
few systems in the Local Group. It is therefore crucial to determine the
reliability and range of applicability of SFHs derived from integrated light
spectroscopy, as this affects our understanding of unresolved galaxies from low
to high redshift.
To evaluate the reliability of current full spectral fitting techniques in
deriving SFHs from integrated light spectroscopy by comparing SFHs from
integrated spectra to those obtained from deep CMDs of resolved stars.
We have obtained a high signal--to--noise (S/N $\sim$ 36.3 per \AA)
integrated spectrum of a field in the bar of the Large Magellanic Cloud (LMC)
using EFOSC2 at the 3.6 meter telescope at La Silla Observatory. For this same
field, resolved stellar data reaching the oMSTO are available. We have compared
the star formation rate (SFR) as a function of time and the age-metallicity
relation (AMR) obtained from the integrated spectrum using {\tt STECKMAP}, and
the CMD using the IAC-star/MinnIAC/IAC-pop set of routines. For the sake of
completeness we also use and discuss other synthesis codes ({\tt STARLIGHT} and
{\tt ULySS}) to derive the SFR and AMR from the integrated LMC spectrum.
We find very good agreement (average differences $\sim$ 4.1 $\%$) between the
SFR(t) and the AMR obtained using {\tt STECKMAP} on the integrated light
spectrum, and the CMD analysis. {\tt STECKMAP} minimizes the impact of the
age-metallicity degeneracy and has the advantage of preferring smooth solutions
to recover complex SFHs by means of a penalized $\chi^2$. [abridged]
The main goal of this work is to explore which elements carry the most information about the birth origin of stars and as such that are best suited for chemical tagging. We explored different techniques to minimize the effect of outlier value lines in the abundances by using Ni abundances derived for 1111 FGK type stars.We evaluated how the limited number of spectral lines can affect the final chemical abundance. Then we were able to make an efficient even footing comparison of the [X/Fe] scatter between the elements that have different number of observable spectral lines in the studied spectra. We found that the most efficient way of calculating the average abundance of elements when several spectral lines are available is to use a weighted mean (WM) where as a weight we considered the distance from the median abundance. This method can be effectively used without removing suspected outlier lines.We showed that when the same number of lines is used to determine chemical abundances, the [X/Fe] star-to-star scatter for iron group and $\alpha$-capture elements is almost the same. On top of this, but at a lower level the largest scatter was observed for Al and the smallest for Cr and Ni. We recommend caution when comparing [X/Fe] scatters among elements that have a different number of spectral lines available. A meaningful comparison is necessary to identify elements that show the largest intrinsic scatter and can be thus used for chemical tagging.
We investigate how the Milky Way tidal field can affect the spatial mixing of multiple stellar populations in the globular cluster NGC 6362. We use $N$-body simulations of multiple population clusters on the orbit of this cluster around the Milky Way. Models of the formation of multiple populations in globular clusters predict that the second population should initially be more centrally concentrated than the first. However, NGC 6362 is comprised of two chemically distinct stellar populations having the same radial distribution. We show that the high mass loss rate experienced on this cluster's orbit significantly accelerates the spatial mixing of the two populations expected from two body relaxation. We also find that for a range of initial second population concentrations, cluster masses, tidal filling factors and fraction of first population stars, a cluster with two populations should be mixed when it has lost 70-80 per cent of its initial mass. These results fully account for the complete spatial mixing of NGC 6362, since, based on its shallow present day mass function, independent studies estimate that the cluster has lost 85 per cent of its initial mass.
We used a sample of super-Earth-like planets detected by the Doppler spectroscopy and transit techniques to explore the dependence of orbital parameters of the planets on the metallicity of their host stars. We confirm the previous results that super-Earths orbiting around metal-rich stars are not observed to be as distant from their host stars as we observe their metal-poor counterparts to be. The orbits of these super-Earths with metal-rich hosts usually do not reach into the Habitable Zone (HZ), keeping them very hot and inhabitable. We found that most of the known planets in the HZ are orbiting their GK-type hosts which are metal-poor. The metal-poor nature of planets in the HZ suggests a high Mg abundance relative to Si and high Si abundance relative to Fe. These results lead us to speculate that HZ planets might be more frequent in the ancient Galaxy and had compositions different from that of our Earth.
The FlashCam group is currently preparing photomultiplier-tube based cameras proposed for the medium-sized telescopes (MST) of the Cherenkov Telescope Array (CTA). The cameras are designed around the FlashCam readout concept which is the first fully-digital readout system for Cherenkov cameras, based on commercial FADCs and FPGAs as key components for the front-end electronics modules and a high performance camera server as back-end. This contribution describes the progress of the full-scale FlashCam camera prototype currently under construction, as well as performance results also obtained with earlier demonstrator setups. Plans towards the production and implementation of FlashCams on site are also briefly presented.
Mounting evidence suggests that the TeV-PeV neutrino flux detected by the IceCube telescope has mainly an extragalactic origin. If such neutrinos are primarily produced by a single class of astrophysical sources via hadronuclear ($pp$) interactions, a similar flux of gamma-ray photons is expected. For the first time, we employ tomographic constraints to pinpoint the origin of the IceCube neutrino events by analyzing recent measurements of the cross correlation between the distribution of GeV gamma rays, detected by the Fermi satellite, and several galaxy catalogs in different redshift ranges. We find that the corresponding bounds on the neutrino luminosity density are up to one order of magnitude tighter than those obtained by using only the spectrum of the gamma-ray background, especially for sources with mild redshift evolution. In particular, our method excludes any hadronuclear source with a spectrum softer than $E^{-2.1}$ as a main component of the neutrino background, if its evolution is slower than $(1+z)^3$. Starburst galaxies, if able to accelerate and confine cosmic rays efficiently, satisfy both spectral and tomographic constraints.
The CDMS low ionization threshold experiment (CDMSlite) uses cryogenic germanium detectors operated at a relatively high bias voltage to amplify the phonon signal in the search for weakly interacting massive particles (WIMPs). Results are presented from the second CDMSlite run with an exposure of 70 kg days, which reached an energy threshold for electron recoils as low as 56 eV. A fiducialization cut reduces backgrounds below those previously reported by CDMSlite. New parameter space for the WIMP-nucleon spin-independent cross section is excluded for WIMP masses between 1.6 and 5.5 GeV/$c^2$.
We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARBEAR experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR non-detection translates into a 95% confidence level (C.L.) upper limit of 93 nano-Gauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at sub-degree scales.
The Cherenkov Telescope Array (CTA) is an international project for a next-generation ground-based gamma-ray observatory. CTA, conceived as an array of tens of imaging atmospheric Cherenkov telescopes, comprising small, medium and large-size telescopes, is aiming to improve on the sensitivity of current-generation experiments by an order of magnitude and provide energy coverage from 20 GeV to more than 300 TeV. The Schwarzschild-Couder (SC) medium-size candidate telescope model features a novel aplanatic two-mirror optical design capable of a wide field-of-view with significantly improved imaging resolution as compared to the traditional Davis-Cotton optics design. Achieving this imaging resolution imposes strict alignment requirements to be accomplished by a dedicated alignment system. In this contribution we present the status of the development of the SC optical alignment system, soon to be materialized in a full-scale prototype SC medium-size telescope at the Fred Lawrence Whipple Observatory in southern Arizona.
The Decadal IRAC Bootes Survey (DIBS) is a mid-IR variability survey of the ~9 sq. deg. of the NDWFS Bootes Field and extnds the time baseline of its predecessor, the Spitzer Deep, Wide-Field Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a million mostly extragalactic sources. In the mid-IR color-color plane, sources with quasar colors constitute the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% have the colors of galaxies. Adding the fifth epoch doubles the number of variable AGNs for the same false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR variability of ~1500 spectroscopically confirmed AGNs using single power-law structure functions, which we find to be steeper (index $\gamma=0.45$) than in the optical ($\gamma=0.3$), leading to much lower amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN structure function slope $\gamma$ seems to be uncorrelated with both the luminosity and rest-frame wavelength, while the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame wavelength.
The excess of positrons in cosmic rays above $\sim$10 GeV has been a puzzle since it was discovered. Possible interpretations of the excess have been suggested, including acceleration in a local supernova remnant or annihilation of dark matter particles. To discriminate between these scenarios, the positron fraction must be measured at higher energies. One technique to perform this measurement is using the Earth-Moon spectrometer: observing the deflection of positron and electron moon shadows by the Earth's magnetic field. The measurement has been attempted by previous imaging atmospheric Cherenkov telescopes without success. The Cherenkov Telescope Array (CTA) will have unprecedented sensitivity and background rejection that could make this measurement successful for the first time. In addition, the possibility of using silicon photomultipliers in some of the CTA telescopes could greatly increase the feasibility of making observations near the moon. Estimates of the capabilities of CTA to measure the positron fraction using simulated observations of the moon shadow will be presented.
We present an initial report on 160 simulations of a highly simplified model of the post-bounce supernova environment in three spatial dimensions (3D). We set different values of a parameter characterizing the impact of nuclear dissociation at the stalled shock in order to regulate the post-shock fluid velocity, thereby determining the relative importance of convection and the stationary accretion shock instability (SASI). While our convection-dominated runs comport with the paradigmatic notion of a `critical neutrino luminosity' for explosion at a given mass accretion rate (albeit with a nontrivial spread in explosion times just above threshold), the outcomes of our SASI-dominated runs are much more stochastic: a sharp threshold critical luminosity is `smeared out' into a rising probability of explosion over a $\sim 20\%$ range of luminosity. We also find that the SASI-dominated models are able to explode with 3 to 4 times less efficient neutrino heating, indicating that progenitor properties, and fluid and neutrino microphysics, conducive to the SASI would make the neutrino-driven explosion mechanism more robust.
Type I X-ray bursts are produced by thermonuclear runaways that develop on accreting neutron stars. Once one location ignites, the flame propagates across the surface of the star. Flame propagation is fundamental in order to understand burst properties like rise time and burst oscillations. Previous work quantified the effects of rotation on the front, showing that the flame propagates as a deflagration and that the front strongly resembles a hurricane. However the effect of magnetic fields was not investigated, despite the fact that magnetic fields strong enough to have an effect on the propagating flame are expected to be present on many bursters. In this paper we show how the coupling between fluid layers introduced by an initially vertical magnetic field plays a decisive role in determining the character of the fronts that are responsible for the Type I bursts. In particular, on a star spinning at 450 Hz (typical among the bursters) we test seed magnetic fields of $10^{7} - 10^{10}$ G and find that for the medium fields the magnetic stresses that develop during the burst can speed up the velocity of the burning front, bringing the simulated burst rise time close to the observed values. By contrast, in a magnetic slow rotator like IGR J17480--2446, spinning at 11 Hz, a seed field $\gtrsim 10^9$ G is required to allow localized ignition and the magnetic field plays an integral role in generating the burst oscillations observed during the bursts.
(abridged) Context: The envelopes of molecular gas around embedded low-mass
protostars show different chemistries, which can be used to trace their
formation history and physical conditions. The excitation of some molecular
species can also be used to trace these physical conditions, making it possible
to constrain e.g. sources of heating and excitation.
Aims: To study the range of influence of an intermediate-mass Herbig Be
protostar, and to find what chemical and physical impact feedback effects from
the environment may have on embedded protostars.
Methods: We follow up on an earlier line survey of the Class 0/I source R CrA
IRS7B in the 0.8 mm window with an unbiased line survey of the same source in
the 1.3 mm window using the APEX telescope. We also study the excitation of the
key species H2CO, CH3OH, and c-C3H2 in a complete sample of the 18 embedded
protostars in the Corona Australis star-forming region. Radiative transfer
models are used to establish abundances of the molecular species.
Results: We detect line emission from 20 molecular species (32 including
isotopologues) in the two surveys. The most complex species detected are CH3OH,
CH3CCH, CH3CHO, and CH3CN. Several complex organics are significantly
under-abundant in comparison with "hot corino" protostars. The H2CO
temperatures of the sources in the region decrease with the distance to the
Herbig Be star R CrA, whereas the c-C3H2 temperatures remain constant across
the star-forming region.
Conclusions: The high H2CO temperatures observed towards objects close to R
CrA suggest that this star has a sphere of influence of several 10000 AU in
which it increases the temperature of the molecular gas to 30-50 K through
irradiation. The chemistry in the IRS7B envelope differs significantly from
many other embedded protostars, which could be an effect of the external
irradiation from R CrA.
By extending our recent framework to describe the tidal deformations of a spinning compact object, we compute for the first time the tidal Love numbers of a spinning neutron star to linear order in the angular momentum. The spin of the object introduces couplings between electric and magnetic distortions and new classes of spin-induced ("rotational") tidal Love numbers emerge. We focus on stationary tidal fields, which induce axisymmetric perturbations. We present the perturbation equations for both electric-led and magnetic-led rotational Love numbers for generic multipoles and explicitly solve them for various tabulated equations of state and for a tidal field with an electric (even parity) and magnetic (odd parity) component with l=2,3,4. For a binary system close to the merger, various components of the tidal field become relevant. In this case we find that an octupolar magnetic tidal field can significantly modify the mass quadrupole moment of a neutron star. Preliminary estimates, assuming a spin parameter \chi~0.05, show modifications of the order of 10% relative to the static case. Furthermore, the rotational Love numbers as functions of the moment of inertia are much more sensitive to the equation of state than in the static case, where approximate universal relations at the percent level exist. For a neutron-star binary approaching the merger, we estimate that the approximate universality of the induced mass quadrupole moment deteriorates from 1% in the static case to roughly 5% when \chi~0.05. Our results suggest that spin-tidal couplings can introduce important corrections to the gravitational waveforms of spinning neutron-star binaries approaching the merger.
A frequentist asymptotic expansion method for error estimation is employed for a network of gravitational wave detectors to assess the capability of gravitational wave observations, with Adv. LIGO and Adv. Virgo, to distinguish between the post-Einsteinian (ppE) description of coalescing binary systems and that of GR. When such errors are smaller than the parameter value, there is possibility to detect these violations from GR. A parameter space with inclusion of dominant dephasing ppE parameters is used for a study of first- and second-order (co)variance expansions, focusing on the inspiral stage of a nonspinning binary system of zero eccentricity detectible through Adv. LIGO and Adv. Virgo. Our procedure is more reliable than frequentist studies based only on Fisher information estimates and complements Bayesian studies. Second-order asymptotics indicate the possibility of constraining deviations from GR in low-SNR ($\rho \sim 15-17$) regimes. The errors on $\beta$ also increase errors of other parameters such as the chirp mass $\mathcal{M}$ and symmetric mass ratio $\eta$. Application is done to existing alternative theories of gravity, which include modified dispersion relation of the waveform, non-spinning models of quadratic modified gravity, and dipole gravitational radiation (i.e., Brans-Dicke type) modifications.
It remains to be determined experimentally if massive neutrinos are Majorana or Dirac particles. In this connection, it has been recently suggested that the detection of cosmic neutrino background of left-handed neutrinos $\nu^{}_{\rm L}$ and right-handed antineutrinos $\overline{\nu}^{}_{\rm R}$ in future experiments of neutrino capture on beta-decaying nuclei (e.g., $\nu^{}_e + {^3{\rm H}} \to {^3}{\rm He} + e^-$ for the PTOLEMY experiment) is likely to distinguish between Majorana and Dirac neutrinos, since the capture rate is twice larger in the former case. In this paper, we investigate the possible impact of right-handed neutrinos on the capture rate, assuming that neutrinos are Dirac particles and both right-handed neutrinos $\nu^{}_{\rm R}$ and left-handed antineutrinos $\overline{\nu}^{}_{\rm L}$ can be efficiently produced in the early Universe. It turns out that the capture rate can be enhanced at most by $28\%$ due to the presence of relic $\nu^{}_{\rm R}$ and $\overline{\nu}^{}_{\rm L}$ with a total number density of $95~{\rm cm}^{-3}$, which should be compared to the number density $336~{\rm cm}^{-3}$ of cosmic neutrino background. The enhancement has actually been limited by the latest cosmological and astrophysical bounds on the effective number of neutrino generations $N^{}_{\rm eff} = 3.14^{+0.44}_{-0.43}$ at the $95\%$ confidence level. Moreover, two possible scenarios have been proposed for thermal production of right-handed neutrinos in the early Universe.
Spacetime with general linear vector distortion is introduced. Thus, the torsion and the nonmetricity of the affine connection are assumed to be proportional to a vector field (and not its derivatives). The resulting two-parameter family of non-Riemannian geometries generalises the conformal Weyl geometry and some other interesting special cases. Taking into account the leading order quadratic curvature correction to the Einstein-Hilbert action results uniquely in the one-parameter extension of the Starobinsky inflation known as the alpha-attractor.
We report the measurement of muons and muon-induced phosphorescence in DM-Ice17, a NaI(Tl) direct detection dark matter experiment at the South Pole. Muons are identified by the observed pulse shape and large energy deposition of their interaction in the crystals. The measured muon rate in DM-Ice17 is 2.93 $\pm$ 0.04 $\mu$/crystal/day with a modulation amplitude of 12.3 $\pm$ 1.7%, consistent with expectation. Following muon interactions, we observe long-lived phosphorescence in the NaI(Tl) crystals with a decay time of 5.5 $\pm$ 0.5 s. The prompt energy deposited by a muon is correlated to the amount of delayed phosphorescence, the brightest of which consist of tens of millions of photons. As they are distributed over tens of seconds, the rate and timing structure of photon arrivals do not mimic a scintillation signal above 2 keV$_\mathrm{ee}$. While the properties of phosphorescence vary between individual crystals, the annually-modulating signal observed by DAMA cannot be accounted for by phosphorescence with the characteristics observed in DM-Ice17.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
Ultraviolet and optical spectra of interstellar gas along the lines of sight to nearby stars have been interpreted by Redfield & Linsky (2008) and previous studies as a set of discrete warm, partially ionized clouds each with a different flow vector, temperature, and metal depletion. Recently, Gry & Jenkins (2014) have proposed a fundamentally different model consisting of a single cloud with nonrigid flows filling space out to 9 parsecs from the Sun that they propose better describes the local ISM. Here we test these fundamentally different morphological models against the spatially unbiased Malamut et al. (2014) spectroscopic data set, and find that the multiple cloud morphology model provides a better fit to both the new and old data sets. The detection of three or more velocity components along the lines of sight to many nearby stars, the presence of nearby scattering screens, the observed thin elongated structures of warm interstellar gas, and the likely presence of strong interstellar magnetic fields also support the multiple cloud model. The detection and identification of intercloud gas and the measurement of neutral hydrogen density in clouds beyond the Local Interstellar Cloud could provide future morphological tests.
In order to explain rapid light curve variability in the context of gamma-ray bursts (GRBs) and jets from active galactic nuclei (AGNs), several authors have proposed the existence of "blobs" or "minijets" that move with relativistic speed relative to the main flow of the jet. Here we consider the possibility that these minijets, instead of being isotropically distributed in the co-moving frame of the jet, form primarily perpendicular to the direction of the flow. This anisotropic collection of minijets yields two robust features. First, the main burst of emission is significantly delayed compared with the isotropic case. This delay allows for the peak of the afterglow emission to appear during the prompt emission, in contrast to the simplest isotropic model, where the afterglow peak appears at or after the end of the main burst. Second, the flux decline following the end of the main burst of emission will be steeper than the isotropic case. We find that these two features are realized in the case of GRBs: 1. The peak of most of the GeV light curves (ascribed to the external shock model) appear during the prompt emission phase. 2. Many X-ray light curves exhibit a period of steep decay, which is faster than that predicted by the standard isotropic case. This emphasizes the importance of anisotropy in the gamma-ray generation mechanism in GRBs, and possibly in other relativistic flows.
We present preliminary diameters and albedos for 7,959 asteroids detected in the first year of the NEOWISE Reactivation mission. 201 are near-Earth asteroids (NEAs). 7,758 are Main Belt or Mars-crossing asteroids. 17% of these objects have not been previously characterized using WISE or NEOWISE thermal measurements. Diameters are determined to an accuracy of ~20% or better. If good-quality H magnitudes are available, albedos can be determined to within ~40% or better.
Recent observations of the Lyman-alpha forest show large-scale spatial variations in the intergalactic Lyman-alpha opacity that grow rapidly with redshift at z>5, far in excess of expectations from empirically motivated models. Previous studies have attempted to explain this excess with spatial fluctuations in the ionizing background, but found that this required either extremely rare sources or problematically low values for the mean free path of ionizing photons. Here we report that much -- or potentially all -- of the observed excess likely arises from residual spatial variations in temperature that are an inevitable byproduct of a patchy and extended reionization process. The amplitude of opacity fluctuations generated in this way depends on the timing and duration of reionization. If the entire excess is due to temperature variations alone, the observed fluctuation amplitude favors a late-ending but extended reionization process that was roughly half complete by z~9 and that ended at z~6. In this scenario, the highest opacities occur in regions that reionized earliest, since they have had the most time to cool, while the lowest opacities occur in the warmer regions that reionized most recently. This correspondence potentially opens a new observational window into patchy reionization.
Merging galaxy systems provides observational evidence of the existence of dark matter and constraints on its properties. Therefore, statistical uniform samples of merging systems would be a powerful tool for several studies. In this work we presents a new methodology for merging systems identification and the results of its application to galaxy redshift surveys. We use as starting point a mock catalogue of galaxy systems, identified using traditional FoF algorithms, which experienced a major merger as indicated by its merger tree. Applying machine learning techniques in this training sample, and using several features computed from the observable properties of galaxy members, it is possible to select galaxy groups with a high probability of have been experienced a major merger. Next we apply clustering techniques on galaxy members in order to reconstruct the properties of the haloes involved in such merger. This methodology provides a highly reliable sample of merging systems with low contamination and precise recovered properties. We apply our techniques in samples of galaxy systems obtained from SDSS-DR7, WINGS and HeCS. Our results recover previously known merging systems and provide several new candidates. We present its measured properties and discuss future analysis on current and forthcoming samples.
Polarization measurements from relativistic outflows are a valuable tool to probe the geometry of the emission region and the microphysics of the particle distribution. Indeed, the polarization level depends on: (i) the local magnetic field orientation, (ii) the geometry of the emitting region with respect to the line of sight, and (iii) the electron pitch-angle distribution. Here we consider optically thin synchrotron emission and we extend the theory of circular polarization from a point source to an extended radially expanding relativistic jet. We present numerical estimates for both linear and circular polarization in such systems. We consider different configurations of the magnetic field, spherical and jetted outflows, isotropic and anisotropic pitch-angle distributions, and outline the difficulty in obtaining the reported high level of circular polarization observed in the afterglow of GRB 121024A. We conclude that the origin of the observed polarization cannot be intrinsic to an optically thin synchrotron process, even when the electron pitch-angle distribution is extremely anisotropic.
Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end.
We study transients produced by equatorial disk-like outflows from catastrophically mass-losing binary stars with an asymptotic velocity and energy deposition rate near the inner edge which are proportional to the binary escape velocity v_esc. As a test case, we present the first smoothed-particle radiation-hydrodynamics calculations of the mass loss from the outer Lagrange point with realistic equation of state and opacities. The resulting spiral stream becomes unbound for binary mass ratios 0.06 < q < 0.8. For synchronous binaries with non-degenerate components, the spiral-stream arms merge at a radius of ~10a, where a is the binary semi-major axis, and the accompanying shock thermalizes 10-20% of the kinetic power of the outflow. The mass-losing binary outflows produce luminosities proportional to the mass loss rate and v_esc, reaching up to ~10^6 L_Sun. The effective temperatures depend primarily on v_esc and span 500 < T_eff < 6000 K. Dust readily forms in the outflow, potentially in a catastrophic global cooling transition. The appearance of the transient is viewing angle-dependent due to vastly different optical depths parallel and perpendicular to the binary plane. The predicted peak luminosities, timescales, and effective temperatures of mass-losing binaries are compatible with those of many of the class of recently-discovered red transients such as V838 Mon and V1309 Sco. We predict a correlation between the peak luminosity and the outflow velocity, which is roughly obeyed by the known red transients. Outflows from mass-losing binaries can produce luminous (10^5 L_Sun) and cool (T_eff < 1500 K) transients lasting a year or longer, as has potentially been detected by Spitzer surveys of nearby galaxies.
We present a realistic modeling of the NVSS number count dipole across the sky. The modeling relies on mock catalogues generated within the context of $\Lambda$CDM cosmology in the linear regime of structure formation. After removal of the solar motion dipole from the observation, the mocks show that the remaining signal is mostly (70\%) due to contribution from large scale structure within $\sim 500$Mpc ($z\sim0.1$). The amplitude of this contribution depends on the bias factor of the NVSS galaxies. An effective bias factor $b(z<0.1)<2.0$ is ruled out at the $\sim 2.8\sigma$ significance by the comparison between the model and observed NVSS dipole. The mismatch is toned down to $\sim 2.3\sigma$ for $b(z)=3$.
Radio halos are synchrotron radio sources detected in some massive galaxy clusters. Their Mpc-size indicates that (re)acceleration processes are taking place in the host cluster. X-ray catalogues of galaxy clusters have been used in the past to search for radio halos and to understand their connection with cluster-cluster mergers and with the thermal component of the intra-cluster medium. More recently, the Sunyaev-Zel'dovich effect has been proven to be a better route to search for massive clusters in a wider redshift range. With the aim of discovering new radio halos and understanding their connection with cluster-cluster mergers, we have selected from the Planck Early source catalog the most massive clusters, and we have observed with the Giant Metrewave Radio Telescope at 323 MHz those objects for which deep observations were not available. We have discovered new peculiar radio emission in three of the observed clusters finding: (i) a radio halo in the cluster RXCJ0949.8+1708; (ii) extended emission in Abell 1443 that we classify as a radio halo plus a radio relic, with a bright filament embedded in the radio halo; (iii) low-power radio emission is found in CIZA J1938.3+5409 which is ten times below the radio - X-ray correlation, and represents the first direct detection of the radio emission in the "upper-limit" region of the radio - X-ray diagram. We discuss the properties of these new radio halos in the framework of theoretical models for the radio emission.
We conduct a comprehensive numerical study of the orbital dependence of harassment on early-type dwarfs consisting of 168 different orbits within a realistic, Virgo-like cluster, varying in eccentricity and pericentre distance. We find harassment is only effective at stripping stars or truncating their stellar disks for orbits that enter deep into the cluster core. Comparing to the orbital distribution in cosmological simulations, we find that the majority of the orbits (more than three quarters) result in no stellar mass loss. We also study the effects on the radial profiles of the globular cluster systems of early-type dwarfs. We find these are significantly altered only if harassment is very strong. This suggests that perhaps most early-type dwarfs in clusters such as Virgo have not suffered any tidal stripping of stars or globular clusters due to harassment, as these components are safely embedded deep within their dark matter halo. We demonstrate that this result is actually consistent with an earlier study of harassment of dwarf galaxies, despite the apparent contradiction. Those few dwarf models that do suffer stellar stripping are found out to the virial radius of the cluster at redshift=0, which mixes them in with less strongly harassed galaxies. However when placed on phase-space diagrams, strongly harassed galaxies are found offset to lower velocities compared to weakly harassed galaxies. This remains true in a cosmological simulation, even when halos have a wide range of masses and concentrations. Thus phase-space diagrams may be a useful tool for determining the relative likelihood that galaxies have been strongly or weakly harassed.
An equation of state of cold nuclear matter with an arbitrary isotopic composition is studied within a relativistic mean-field approach with hadron masses and coupling constants depending self-consistently on the scalar mean-field. All hadron masses decrease universally with the scalar field growth, whereas meson-nucleon coupling constants can vary differently. More specifically we focus on two modifications of the KVOR model studied previously. One extension of the model (KVORcut) demonstrates that the equation of state stiffens if the increase of the scalar-field magnitude with the density is bounded from above at some value for baryon densities above the saturation nuclear density. This can be realized if the nucleon vector-meson coupling constant changes rapidly as a function of the scalar field slightly above the desired value. The other version of the model (MKVOR) utilizes a smaller value of the nucleon effective mass at the nuclear saturation density and a saturation of the scalar field in the isospin asymmetric matter induced by a strong variation of the nucleon isovector-meson coupling constant as function of the scalar field. A possibility of hyperonization of the matter in neutron star interiors is incorporated. Our equations of state fulfill majority of known empirical constraints including the pressure-density constraint from heavy-ion collisions, direct Urca constraint, gravitational-baryon mass constraint and the constraint on the maximum mass of the neutron stars.
We explore a plausible mechanism that the hemispherical power asymmetry in the CMB is produced by the spatial variation of the primordial sound speed parameter. We suggest that in a generalized approach of the $\delta N$ formalism the local e-folding number may depend on some other primordial parameters besides the initial values of inflaton. Here the $\delta N$ formalism is extended by considering the effects of a spatially varying sound speed parameter caused by a super-Hubble perturbation of a light field. Using this generalized $\delta N$ formalism, we systematically calculate the asymmetric primordial spectrum in the model of multi-speed inflation by taking into account the constraints of primordial non-Gaussianities. We further discuss specific model constraints, and the corresponding asymmetry amplitudes are found to be scale-dependent, which can accommodate current observations of the power asymmetry at different length scales.
We have obtained ALMA Band 7 observations of the FU Ori outburst system at 0.6"x0.5" resolution to measure the link between the inner disk instability and the outer disk through sub-mm continuum and molecular line observations. Our observations detect continuum emission which can be well modeled by two unresolved sources located at the position of each binary component. The interferometric observations recover the entire flux reported in previous single-dish studies, ruling out the presence of a large envelope. Assuming that the dust is optically thin, we derive disk dust masses of $2\times 10^{-4}$M$_{\odot}$ and $8\times 10^{-5}$M$_{\odot}$, for the north and south components respectively. We place limits on the disks' radii of $r<$45 AU. We report the detection of molecular emission from $^{12}$CO(3-2), HCO$^{+}$(4-3) and from HCN(4-3). The $^{12}$CO appears widespread across the two binary components, and is slightly more extended than the continuum emission. The denser gas tracer HCO$^{+}$ peaks close to the position of the southern binary component, while HCN appears peaked at the position of the northern component. This suggests that the southern binary component is embedded in denser molecular material, consistent with previous studies that indicate a heavily reddened object. At this angular resolution any interaction between the two unresolved disk components cannot be disentangled. Higher resolution images are vital to understanding the process of star formation via rapid accretion FU Ori-type episodes.
We present new spectroscopic observations that are part of our continuing monitoring campaign of 88 quasars at z<0.7 whose broad H-beta lines are offset from their systemic redshifts by a few thousand km/s. These quasars have been considered candidates for hosting supermassive black hole binaries (SBHBs) by analogy with single-lined spectroscopic binary stars. We present the data and describe our improved analysis techniques, which include an extensive evaluation of uncertainties. We also present a variety of measurements from the spectra that are of general interest and will be useful in later stages of our analysis. Additionally, we take this opportunity to study the variability of the optical continuum and integrated flux of the broad H-beta line. We compare the variability properties of the SBHB candidates to those of a sample of typical quasars with similar redshifts and luminosities observed multiple times during the Sloan Digital Sky Survey. We find that the variability properties of the two samples are similar (variability amplitudes of 10-30% on time scales of approximately 1-7 years) and that their structure functions can be described by a common model with parameters characteristic of typical quasars. These results suggest that the broad-line regions of SBHB candidates have a similar extent as those of typical quasars. We discuss the implications of this result for the SBHB scenario and ensuing constraints on the orbital parameters.
We report the detection of steady radio emission from the known X-ray source X9 in the globular cluster 47 Tuc. With a double-peaked C IV emission line in its ultraviolet spectrum providing a clear signature of accretion, this source had been previously classified as a cataclysmic variable. In deep ATCA imaging from 2010 and 2013, we identified a steady radio source at both 5.5 and 9.0 GHz, with a radio spectral index (defined as $S_{\nu}\propto\nu^{\alpha}$) of $\alpha=-0.4\pm0.4$. Our measured flux density of $42\pm4$ microJy/beam at 5.5 GHz implies a radio luminosity ($\nu L_{\nu}$) of 5.8e27 erg/s, significantly higher than any previous radio detection of an accreting white dwarf. Transitional millisecond pulsars, which have the highest radio-to-X-ray flux ratios among accreting neutron stars (still a factor of a few below accreting black holes at the same X-ray luminosity), show distinctly different patterns of X-ray and radio variability than X9. When combined with archival X-ray measurements, our radio detection places 47 Tuc X9 very close to the radio/X-ray correlation for accreting black holes, and we explore the possibility that this source is instead a quiescent stellar-mass black hole X-ray binary. The nature of the donor star is uncertain; although the luminosity of the optical counterpart is consistent with a low-mass main sequence donor star, the mass transfer rate required to produce the high quiescent X-ray luminosity of 1e33 erg/s suggests the system may instead be ultracompact, with an orbital period of order 25 minutes. This is the fourth quiescent black hole candidate discovered to date in a Galactic globular cluster, and the only one with a confirmed accretion signature from its optical/ultraviolet spectrum.
Two types of correlations between the radio and X-ray luminosities ($L_R$ and $L_X$) of black hole sources has been found. For the traditional type of sources, the correlation can be described by a single power-law. For the other type of sources, while the correlation can still be described by power-law forms, it consists three branches according to the X-ray luminosity, with different power-law indexes. In this paper, we try to explain these correlations in the framework of the coupled accretion-jet model. We attribute the difference between these two types of sources to the difference in the value of viscous parameter $\alpha$. For the "single power-law" sources, their $\alpha$ is high; so their accretion is always in the mode of ADAF (advection-dominated accretion flow) for the whole range of X-ray luminosity. For those "hybrid power-law" sources, the value of $\alpha$ is small so their accretion modes change from ADAF to LHAF (luminous hot accretion flow) to two-phase accretion as the accretion rate increases. Because the radiative efficiency of the hot accretion flow on the mass accretion rate is different for these three accretion modes, they will lead to different power-law indexes in the $L_R$ -- $L_X$ correlation. The reason of the different $\alpha$ may be because of the different configuration of magnetic field in the accretion material coming from the companion stars. Constraints on the ratio of the mass lost rate into the jet and the mass accretion rate in the accretion flow have been obtained, which can be tested in future by radiative magnetohydrodynamic (MHD) numerical simulations on jet formation.
Context: Pre-merger interactions between galaxies can induce significant changes in the morphologies and kinematics of the stellar and ISM components. Large amounts of gas and stars are often found to be disturbed or displaced as tidal debris. This debris then evolves, sometimes forming stars and occasionally tidal dwarf galaxies. Here we present results from our HI study of Arp 65, an interacting pair hosting extended HI tidal debris. Aims: In an effort to understand the evolution of tidal debris produced by interacting pairs of galaxies, including in situ star and tidal dwarf galaxy formation, we are mapping HI in a sample of interacting galaxy pairs. The Arp 65 pair is one of them. Methods: Our resolved HI 21 cm line survey is being carried out using the Giant Metrewave Radio Telescope (GMRT). We used our HI survey data as well as available SDSS optical, Spitzer infra-red and GALEX UV data to study the evolution of the tidal debris and the correlation of HI with the star-forming regions within it. Results: In Arp 65 we see a high impact pre-merger interaction involving a pair of massive galaxies (NGC 90 and NGC 93) that have a stellar mass ratio of ~ 1:3. The interaction, which probably occurred ~ 1.0 -- 2.5 $\times$ 10$^8$ yr ago, appears to have displaced a large fraction of the HI in NGC 90 (including the highest column density HI) beyond its optical disk. We also find extended ongoing star formation in the outer disk of NGC 90. In the major star-forming regions, we find the HI column densities to be ~ 4.7 $\times$ 10$^{20}$ cm$^{-2}$ or lower. But no signature of star formation was found in the highest column density HI debris, SE of NGC 90. This indicates conditions within the highest column density HI debris remain hostile to star formation and it reaffirms that high HI column densities may be a necessary but not sufficient criterion for star formation.
In previous studies, it has been shown that the long term time average jet power, $\overline{Q}$, is correlated with the spectral index in the extreme ultraviolet (EUV), $\alpha_{EUV}$ (defined by $F_{\nu} \sim \nu^{-\alpha_{EUV}}$ computed between 700\AA\, and 1100\AA\,). Larger $\overline{Q}$ tends to decrease the EUV emission. This is a curious relationship because it connects a long term average over $\sim 10^{6}$ years with an instantaneous measurement of the EUV. The EUV appears to be emitted adjacent to the central supermassive black hole and the most straightforward explanation of the correlation is that the EUV emitting region interacts in real time with the jet launching mechanism. Alternatively stated, the $\overline{Q}$ - $\alpha_{EUV}$ correlation is a manifestation of a contemporaneous (real time) jet power, $Q(t)$, correlation with $\alpha_{EUV}$. In order to explore this possibility, this paper considers the time variability of the strong radio jet of the quasar 1442+101 that is not aberrated by strong Doppler enhancement. This high redshift (z = 3.55) quasar is uniquely suited for this endeavor as the EUV is redshifted into the optical observing window allowing for convenient monitoring. More importantly, it is bright enough to be seen through the Lyman forest and its radio flux is strong enough that it has been monitored frequently. Quasi-simultaneous monitoring (five epochs spanning $\sim 40$ years) show that increases in $Q(t)$ correspond to decreases in the EUV as expected.
Simultaneous observations at multiple frequencies have the potential to overcome the fundamental limitation imposed by the atmospheric propagation in mm-VLBI observations. The propagation effects place a severe limit in the sensitivity achievable in mm-VLBI, reducing the time over which the signals can be coherently combined, and preventing the use of phase referencing and astrometric measurements. We carried out simultaneous observations at 22, 43, 87 and 130\,GHz of a group of five AGNs, the weakest of which is $\sim$200\,mJy at 130\,GHz, with angular separations ranging from 3.6 to 11 degrees, with the KVN. We analysed this data using the Frequency Phase Transfer (FPT) and the Source Frequency Phase Referencing (SFPR) techniques, which use the observations at a lower frequency to correct those at a higher frequency. The results of the analysis provide an empirical demonstration of the increase in the coherence times at 130\,GHz from a few tens of seconds to about twenty minutes, with FPT, and up to many hours with SFPR. Moreover the astrometric analysis provides high precision measurements for a bona fide astrometric image registration of the sources at the four bands, including, for the first time, high precision astrometric measurements at 130 GHz. Finally we demonstrate a method for the generalised decomposition of the relative measurements into the position shifts of the individual sources for a bona fide astrometric registration of the maps at the four frequencies.
Superclusters of galaxies can be defined kinematically from local evaluations of the velocity shear tensor. The location where the smallest eigenvalue of the shear is positive and maximal defines the center of a basin of attraction. Velocity and density fields are reconstructed with Wiener Filter techniques. Local velocities due to the density field in a restricted region can be separated from external tidal flows, permitting the identification of boundaries separating inward flows toward a basin of attraction and outward flows. This methodology was used to define the Laniakea Supercluster that includes the Milky Way. Large adjacent structures include Perseus-Pisces, Coma, Hercules, and Shapley but current kinematic data are insufficient to capture their full domains. However there is a small region trapped between Laniakea, Perseus-Pisces, and Coma that is close enough to be reliably characterized and that satisfies the kinematic definition of a supercluster. Because of its shape, it is given the name the Arrowhead Supercluster. This entity does not contain any major clusters. A characteristic dimension is ~25 Mpc and the contained mass is only ~10^15 Msun.
We developed a dual-linear-polarization HEMT (High Electron Mobility Transistor) amplifier receiver system of the 45-GHz band (hereafter Z45), and installed it in the Nobeyama 45-m radio telescope. The receiver system is designed to conduct polarization observations by taking the cross correlation of two linearly-polarized components, from which we process full-Stokes spectroscopy. We aim to measure the magnetic field strength through the Zeeman effect of the emission line of CCS ($J_N=4_3-3_2$) toward pre-protostellar cores. A linear-polarization receiver system has a smaller contribution of instrumental polarization components to the Stokes $V$ spectra than that of the circular polarization system, so that it is easier to obtain the Stokes $V$ spectra. The receiver has an RF frequency of 42 $-$ 46 GHz and an intermediate frequency (IF) band of 4$-$8 GHz. The typical noise temperature is about 50 K, and the system noise temperature ranges from 100 K to 150K over the frequency of 42 $-$ 46 GHz. The receiver system is connected to two spectrometers, SAM45 and PolariS. SAM45 is a highly flexible FX-type digital spectrometer with a finest frequency resolution of 3.81 kHz. PolariS is a newly-developed digital spectrometer with a finest frequency resolution of 60 Hz, having a capability to process the full-Stokes spectroscopy. The Half Power Beam Width (HPBW) of the beam was measured to be 37$"$ at 43 GHz. The main beam efficiency of the Gaussian main beam was derived to be 0.72 at 43 GHz. The SiO maser observations show that the beam pattern is reasonably round at about 10 \% of the peak intensity and the side-lobe level was less than 3 \% of the peak intensity. Finally, we present some examples of astronomical observations using Z45.
We present the SAMI Pilot Survey, consisting of integral field spectroscopy
of 106 galaxies across three galaxy clusters, Abell 85, Abell 168 and Abell
2399. The galaxies were selected by absolute magnitude to have $M_r<-20.25$
mag. The survey, using the Sydney-AAO Multi-object Integral field spectrograph
(SAMI), comprises observations of galaxies of all morphological types with 75\%
of the sample being early-type galaxies (ETGs) and 25\% being late-type
galaxies (LTGs). Stellar velocity and velocity dispersion maps are derived for
all 106 galaxies in the sample.
The $\lambda_{R}$ parameter, a proxy for the specific stellar angular
momentum, is calculated for each galaxy in the sample. We find a trend between
$\lambda_{R}$ and galaxy concentration such that LTGs are less concentrated
higher angular momentum systems, with the fast-rotating ETGs (FRs) more
concentrated and lower in angular momentum. This suggests that some dynamical
processes are involved in transforming LTGs to FRs, though a significant
overlap between the $\lambda_{R}$ distributions of these classes of galaxies
implies that this is just one piece of a more complicated picture.
We measure the kinematic misalignment angle, $\Psi$, for the ETGs in the
sample, to probe the intrinsic shapes of the galaxies. We find the majority of
FRs (83\%) to be aligned, consistent with them being oblate spheroids (i.e.
disks). The slow rotating ETGs (SRs), on the other hand, are significantly more
likely to show kinematic misalignment (only 38\% are aligned). This confirms
previous results that SRs are likely to be mildly triaxial systems.
The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic Cap, using the 90Prime wide-field imager on the 2.3m Bok telescope at Kitt Peak. The survey observations started in 2010 and ended in 2013. The final survey area is about 5000 deg2 with a median 5-sigma point source limiting magnitude of about 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS photometric catalog with those of the SDSS and CFHTLS.
The electric current helicity density $\displaystyle \chi=\langle\epsilon_{ijk}b_i\frac{\partial b_k}{\partial x_j}\rangle$ contains six terms, where $b_i$ are components of the magnetic field. Due to the observational limitations, only four of the above six terms can be inferred from solar photospheric vector magnetograms. By comparing the results for simulation we distinguished the statistical difference of above six terms for isotropic and anisotropic cases. We estimated the relative degree of anisotropy for three typical active regions and found that it is of order 0.8 which means the assumption of local isotropy for the observable current helicity density terms is generally not satisfied for solar active regions. Upon studies of the statistical properties of the anisotropy of magnetic field of solar active regions with latitudes and with evolution in the solar cycle, we conclude that the consistency of that assumption of local homogeneity and isotropy requires further analysis in the light of our findings.
We present multi-epoch near-infrared photo-spectroscopic observations of Nova Cephei 2014 and Nova Scorpii 2015, discovered in outburst on 2014 March 8.79 UT and 2015 February 11.84 UT respectively. Nova Cep 2014 shows the conventional NIR characteristics of a Fe II class nova characterized by strong CI, HI and O I lines, whereas Nova Sco 2015 is shown to belong to the He/N class with strong He I, HI and OI emission lines. The highlight of the results consists in demonstrating that Nova Sco 2015 is a symbiotic system containing a giant secondary. Leaving aside the T CrB class of recurrent novae, all of which have giant donors, Nova Sco 2015 is shown to be only the third classical nova to be found with a giant secondary. The evidence for the symbiotic nature is three-fold; first is the presence of a strong decelerative shock accompanying the passage of the nova's ejecta through the giant's wind, second is the H$\alpha$ excess seen from the system and third is the spectral energy distribution of the secondary in quiescence typical of a cool late type giant. The evolution of the strength and shape of the emission line profiles shows that the ejecta velocity follows a power law decay with time ($t^{-1.13 \pm 0.17}$). A Case B recombination analysis of the H I Brackett lines shows that these lines are affected by optical depth effects for both the novae. Using this analysis we make estimates for both the novae of the emission measure $n_e^2L$, the electron density $n_e$ and the mass of the ejecta.
HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\, m^2sr$ for electrons and diffuse photons and $>2\, m^2sr$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $\gamma$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $\gamma$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.
By taking into account the contamination of foreground radiations, we employ the Fisher matrix to forecast the future sensitivity on the tilt of power spectrum of primordial tensor perturbations for several ground-based (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloon-borne (EBEX, Spider) and satellite (CMBPol, COrE, LiteBIRD) experiments of B-mode polarizations. For the fiducial model $n_t=0$, our results show that the satellite experiments give good sensitivity on the tensor tilt $n_t$ to the level $\sigma_{n_t}\lesssim0.1$ for $r\gtrsim2\times10^{-3}$, while the ground-based and balloon-borne experiments give worse sensitivity. By considering the BICEP2/Keck Array and Planck (BKP) constraint on the tensor-to-scalar ratio $r$, we see that it is impossible for these experiments to test the consistency relation $n_t=-r/8$ in the canonical single-field slow-roll inflation models.
The Cherenkov Telescope Array (CTA) is the next generation of ground-based very high energy gamma-ray instruments; the facility will be organized in two arrays, one for each hemisphere. The atmospheric calibration of the CTA telescopes is a critical task. The atmosphere affects the measured Cherenkov yield in several ways: the air-shower development itself, the variation of the Cherenkov angle with altitude, the loss of photons due to scattering and absorption of Cherenkov light out of the camera field-of-view and the scattering of photons into the camera. In this scenario, aerosols are the most variable atmospheric component in time and space and therefore need a continuous monitoring. Lidars are among the most used instruments in atmospheric physics to measure the aerosol attenuation profiles of light. The ARCADE Lidar system is a very compact and portable Raman Lidar system that has been built within the FIRB 2010 grant and is currently taking data in Lamar, Colorado. The ARCADE Lidar is proposed to operate at the CTA sites with the goal of making a first survey of the aerosol conditions of the selected site and to use it as a calibrated benchmark for the other Lidars that will be installed on site. It is proposed for CTA that the ARCADE Lidar will be first upgraded in Italy and then tested in parallel to a Lidar of the EARLINET network in L'Aquila. Upgrades include the addition of the water vapour Raman channel to the receiver and the use of new and better performing electronics. It is proposed that the upgraded system will travel to and characterize both CTA sites, starting from the first selected site in 2016.
We present new multi-band (UBVI) time-series data of helium burning variables in the Carina dwarf spheroidal galaxy. The current sample includes 92 RR Lyrae-six of them are new identifications-and 20 Anomalous Cepheids, one of which is new identification. The analysis of the Bailey diagram shows that the luminosity amplitude of the first overtone component in double-mode variables is located along the long-period tail of regular first overtone variables, while the fundamental component is located along the short-period tale of regular fundamental variables. This evidence further supports the transitional nature of these objects. Moreover, the distribution of Carina double-mode variables in the Petersen diagram (P_1/P_0 vs P_0) is similar to metal-poor globulars (M15, M68), to the dwarf spheroidal Draco and to the Galactic Halo. This suggests that the Carina old stellar population is metal-poor and affected by a small spread in metallicity. We use trigonometric parallaxes for five field RR Lyrae stars to provide an independent estimate of the Carina distance using the observed reddening free Period--Wesenheit [PW, (BV)] relation. Theory and observations indicate that this diagnostic is independent of metallicity. We found a true distance modulus of \mu=20.01\pm0.02 (standard error of the mean) \pm0.05 (standard deviation) mag. We also provided independent estimates of the Carina true distance modulus using four predicted PW relations (BV, BI, VI, BVI) and we found: \mu=(20.08\pm0.007\pm0.07) mag, \mu=(20.06\pm0.006\pm0.06) mag, \mu=(20.07\pm0.008\pm0.08) mag and \mu=(20.06\pm0.006\pm0.06) mag. Finally, we identified more than 100 new SX Phoenicis stars that together with those already known in the literature (340) make Carina a fundamental laboratory to constrain the evolutionary and pulsation properties of these transitional variables.
We report results on multiband observations from radio to gamma-rays of the two radio-loud narrow-line Seyfert 1 (NLSy1) galaxies PKS 2004-447 and J1548+3511. Both sources show a core-jet structure on parsec scale, while they are unresolved at the arcsecond scale. The high core dominance and the high variability brightness temperature make these NLSy1 galaxies good gamma-ray source candidates. Fermi-LAT detected gamma-ray emission only from PKS 2004-447, with a gamma-ray luminosity comparable to that observed in blazars. No gamma-ray emission is observed for J1548+3511. Both sources are variable in X-rays. J1548+3511 shows a hardening of the spectrum during high activity states, while PKS 2004-447 has no spectral variability. A spectral steepening likely related to the soft excess is hinted below 2 keV for J1548+3511, while the X-ray spectra of PKS 2004-447 collected by XMM-Newton in 2012 are described by a single power-law without significant soft excess. No additional absorption above the Galactic column density or the presence of an Fe line is detected in the X-ray spectra of both sources.
Both physical arguments and simulations of the global heliosphere indicate that the tailward heliopause is flattened considerably in the direction perpendicular to the large-scale interstellar magnetic field. Despite this fact, all existing global analytical models of the outer heliosheath's magnetic field assume a circular cross section of the heliotail. To eliminate this inconsistency, we introduce a mathematical procedure by which any analytically or numerically given magnetic field can be deformed in such a way that the cross sections along the heliotail axis attain freely prescribed, spatially dependent values for their total area and aspect ratio. The distorting transformation of this method honors both the solenoidality condition and the stationary induction equation with respect to an accompanying flow field, provided that both constraints were already satisfied for the original magnetic and flow fields prior to the transformation. In order to obtain realistic values for the above parameters, we present the first quantitative analysis of the heliotail's overall distortion as seen in state-of-the-art 3D hybrid MHD-kinetic simulations.
The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.
A multifrequency campaign on the BL Lac object PG 1553+113 was organized by the Whole Earth Blazar Telescope (WEBT) in 2013 April-August, involving 19 optical, two near-IR, and three radio telescopes. The aim was to study the source behaviour at low energies during and around the high-energy observations by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes in April-July. We also analyse the UV and X-ray data acquired by the Swift and XMM-Newton satellites in the same period. The WEBT and satellite observations allow us to detail the synchrotron emission bump in the source spectral energy distribution (SED). In the optical we found a general bluer-when-brighter trend. The X-ray spectrum remained stable during 2013, but a comparison with previous observations suggests that it becomes harder when the X-ray flux increases. The long XMM-Newton exposure reveals a curved X-ray spectrum. In the SED, the XMM-Newton data show a hard near-UV spectrum, while Swift data display a softer shape that is confirmed by previous HST-COS and IUE observations. Polynomial fits to the optical-X-ray SED show that the synchrotron peak likely lies in the 4-30 eV energy range, with a general shift towards higher frequencies for increasing X-ray brightness. However, the UV and X-ray spectra do not connect smoothly. Possible interpretations include: i) orientation effects, ii) additional absorption, iii) multiple emission components, and iv) a peculiar energy distribution of relativistic electrons. We discuss the first possibility in terms of an inhomogeneous helical jet model.
We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.
We present THERMAP, a mid-infrared (8-16 {\mu}m) spectro-imager for space missions to small bodies in the inner solar system, developed in the framework of the MarcoPolo-R asteroid sample return mission. THERMAP is very well suited to characterize the surface thermal environment of a NEO and to map its surface composition. The instrument has two channels, one for imaging and one for spectroscopy: it is both a thermal camera with full 2D imaging capabilities and a slit spectrometer. THERMAP takes advantage of the recent technological developments of uncooled microbolometers detectors, sensitive in the mid-infrared spectral range. THERMAP can acquire thermal images (8-18 {\mu}m) of the surface and perform absolute temperature measurements with a precision better than 3.5 K above 200 K. THERMAP can acquire mid-infrared spectra (8-16 {\mu}m) of the surface with a spectral resolution {\Delta}{\lambda} of 0.3 {\mu}m. For surface temperatures above 350 K, spectra have a signal-to-noise ratio >60 in the spectral range 9-13 {\mu}m where most emission features occur.
The International X-ray Observatory (IXO) is a very ambitious mission, aimed at the X-ray observation of the early Universe. This makes IXO extremely demanding in terms of effective area and angular resolution. In particular, the HEW requirement below 10 keV is 5 arcsec Half-Energy Width (HEW). At higher photon energies, the HEW is expected to increase, and the angular resolution to be correspondingly degraded, due to the increasing relevance of the X-ray scattering off the reflecting surfaces. Therefore, the HEW up to 40 keV is required to be better than 30 arcsec, even though the IXO goal is to achieve an angular resolution as close as possible to 5 arcsec also at this energy. To this end, the roughness of the reflecting surfaces has to not exceed a tolerance, expressed in terms of a surface roughness PSD (Power-Spectral-Density). In this work we provide such tolerances by simulating the HEW scattering term for IXO, assuming a specific configuration for the optical module and different hypotheses on the PSD of mirrors.
Future X-ray telescopes like SIMBOL-X will operate in a wide band of the X-ray spectrum (from 0.1 to 80 keV); these telescopes will extend the optical performances of the existing soft X-ray telescopes to the hard X-ray band, and in particular they will be characterized by a angular resolution (conveniently expressed in terms of HEW, Half-Energy- Width) less than 20 arcsec. However, it is well known that the microroughness of the reflecting surfaces of the optics causes the scattering of X-rays. As a consequence, the imaging quality can be severely degraded. Moreover, the X-ray scattering can be the dominant problem in hard X-rays because its relevance is an increasing function of the photon energy. In this work we consistently apply a numerical method and an analytical one to evaluate the X-ray scattering impact on the HEW of an X-ray optic, as a function of the photon energy: both methods can also include the effects of figure errors in determining the final HEW. A comparison of the results obtained with the two methods for the particular case of the SIMBOL-X X-ray telescope will be presented.
One of the most important parameters defining the angular resolution of an X-ray optical module is its Half-Energy Width (HEW) as a function of the photon energy. Future X-ray telescopes with imaging capabilities (SIMBOL-X, Constellation-X, NeXT, EDGE, XEUS,...) should be characterized by a very good angular resolution in soft (< 10 keV) and hard (> 10 keV) X-rays. As a consequence, an important point in the optics development for these telescopes is the simulation of the achievable HEW for a system of X-ray mirrors. This parameter depends on the single mirror profile and nesting accuracy, but also on the mirrors surface microroughness that causes X-ray Scattering (XRS). In particular, owing to its dependence on the photon energy, XRS can dominate the profile errors in hard X-rays: thus, its impact has to be accurately evaluated in every single case, in order to formulate surface finishing requirements for X-ray mirrors. In this work we provide with some simulations of the XRS term of the HEW for some future soft and hard X-ray telescopes.
Around 2% of all A stars have photospheres depleted in refractory elements. This is hypothesized to arise from a preferential accretion of gas rather than dust, but the specific processes and the origin of the material -- circum- or interstellar -- are not known. The same depletion is seen in 30% of young, disk-hosting Herbig Ae/Be stars. We investigate whether the chemical peculiarity originates in a circumstellar disk. Using a sample of systems for which both the stellar abundances and the protoplanetary disk structure are known, we find that stars hosting warm, flaring group I disks typically have Fe, Mg and Si depletions of 0.5 dex compared to the solar-like abundances of stars hosting cold, flat group II disks. The volatile, C and O, abundances in both sets are identical. Group I disks are generally transitional, having radial cavities depleted in millimetre-sized dust grains, while those of group II are usually not. Thus we propose that the depletion of heavy elements emerges as Jupiter-like planets block the accretion of part of the dust, while gas continues to flow towards the central star. We calculate gas to dust ratios for the accreted material and find values consistent with models of disk clearing by planets. Our results suggest that giant planets of ~0.1 to 10 M_Jup are hiding in at least 30% of Herbig Ae/Be disks.
We discuss the potential of Athena X-ray telescope, in particular of its X-ray Integral Field Unit (X-IFU), for detection of the signal from the light-weight decaying dark matter with mass in the keV range. We show that high energy resolution and large collection area of X-IFU will provide an improvement of sensitivity which will be sufficient for the full test of the neutrino Minimal extension of the Standard Model (nuMSM). Search for the narrow spectral line produced by the decay of the dark matter sterile neutrino in the spectra of dwarf spheroidal galaxies with X-IFU will explore the whole allowed range masses and mixing angles of the nuMSM lightest sterile neutrino and in this way either to find the dark matter signal or rule out the nuMSM model.
We present an analysis of X-ray, Ultraviolet and optical/near-IR photometric data of the transitional millisecond pulsar binary XSSJ12270-4859, obtained at different epochs after the transition to a rotation-powered radio pulsar state. The observations, while confirming the large-amplitude orbital modulation found in previous studies after the state change, also reveal an energy dependence of the amplitudes as well as variations on time scale of months. The amplitude variations are anti-correlated in the X-ray and the UV/optical bands. The average X-ray spectrum is described by a power law with \Gamma index of 1.07(8) without requiring an additional thermal component. The power law index \Gamma varies from 1.2 to 1.0 between superior and inferior conjunction of the neutron star. We interpret the observed X-ray behaviour in terms of synchrotron radiation emitted in an extended intrabinary shock, located between the pulsar and the donor star, which is eclipsed due to the companion orbital motion. The G5 type donor dominates the UV/optical and near-IR emission and is similarly found to be heated up to ? 6500K as in the disc state. The analysis of optical light curves gives a binary inclination 46 < i < 65deg and a mass ratio 0.11 < q <0.26. The donor mass is found to be 0.15 < M2 < 0.36Msun for a neutron star mass of 1.4Msun. The variations in the amplitude of the orbital modulation are interpreted in terms of small changes in the mass flow rate from the donor star. The spectral energy distribution from radio to gamma-rays is composed by multiple contributions that are different from those observed during the accretion-powered state.
Radio and X-ray emission from brown dwarfs suggest that an ionised gas and a
magnetic field with a sufficient flux density must be present. We perform a
reference study for late M-dwarfs, brown dwarfs and giant gas planet to
identify which ultra-cool objects are most susceptible to plasma and magnetic
processes. Only thermal ionisation is considered. We utilise the {\sc
Drift-Phoenix} model grid where the local atmospheric structure is determined
by the global parameters T$_{\rm eff}$, $\log(g)$ and [M/H].
Our results show that it is not unreasonable to expect H$_{\alpha}$ or radio
emission to origin from Brown Dwarf atmospheres as in particular the rarefied
upper parts of the atmospheres can be magnetically coupled despite having low
degrees of thermal gas ionisation. Such ultra-cool atmospheres could therefore
drive auroral emission without the need for a companion's wind or an outgassing
moon. The minimum threshold for the magnetic flux density required for
electrons and ions to be magnetised is well above typical values of the global
magnetic field of a brown dwarf and a giant gas planet. Na$^{+}$, K$^{+}$ and
Ca$^{+}$ are the dominating electron donors in low-density atmospheres (low
log(g), solar metallicity) independent of T$_{\rm eff}$. Mg$^{+}$ and Fe$^{+}$
dominate the thermal ionisation in the inner parts of M-dwarf atmospheres.
Molecules remain unimportant for thermal ionisation. Chemical processes (e.g.
cloud formation) affecting the most abundant electron donors, Mg and Fe, will
have a direct impact on the state of ionisation in ultra-cool atmospheres.
Zeta Ori A is a hot star claimed to host a weak magnetic field, but no clear magnetic detection was obtained so far. In addition, it was recently shown to be a binary system composed of a O9.5I supergiant and a B1IV star. We aim at verifying the presence of a magnetic field in zeta Ori A, identifying to which of the two binary components it belongs (or whether both stars are magnetic), and characterizing the field.Very high signal-to-noise spectropolarimetric data were obtained with Narval at the Bernard Lyot Telescope (TBL) in France. Archival HEROS, FEROS and UVES spectroscopic data were also used. The data were first disentangled to separate the two components. We then analyzed them with the Least-Squares Deconvolution (LSD) technique to extract the magnetic information. We confirm that zeta Ori A is magnetic. We find that the supergiant component zeta Ori Aa is the magnetic component: Zeeman signatures are observed and rotational modulation of the longitudinal magnetic field is clearly detected with a period of 6.829 d. This is the only magnetic O supergiant known as of today. With an oblique dipole field model of the Stokes V profiles, we show that the polar field strength is ~ 140 G. Because the magnetic field is weak and the stellar wind is strong, zeta Ori Aa does not host a centrifugally supported magnetosphere. It may host a dynamical magnetosphere. Its companion zeta Ori Ab does not show any magnetic signature, with an upper limit on the undetected field of $\sim$ 300 G.
During the last few years, the Geneva stellar evolution group has released new grids of stellar models, including the effect of rotation and with updated physical inputs (Ekstr\"om et al. 2012; Georgy et al. 2013a,b). To ease the comparison between the outputs of the stellar evolution computations and the observations, a dedicated tool was developed: the Syclist toolbox (Georgy et al. 2014). It allows to compute interpolated stellar models, isochrones, synthetic clusters, and to simulate the time-evolution of stellar populations.
Dwarf spheroidal galaxies (dSphs) are mostly investigated in the Local Group. DSphs are difficult targets for observations because of their small size and very low surface brightness. Here we measure spectroscopic and photometric parameters of three candidates for isolated dSphs, KKH65=BTS23, KK180, and KK227, outside the Local Group. The galaxies are found to be of low metallicity and low velocity dispersion. They are among the lowest surface brightness objects in the Local Universe. According to the measured radial velocities, metallicities, and structural and photometric parameters, KKH65 and KK227 are representatives of the ultra-diffuse quenched galaxies. KKH65 and KK227 belong to the outer parts of the groups NGC3414 and NGC5371, respectively. KK180 is located in the Virgo cluster infall region.
Motivated by the current constraints on the epoch of reionisation from recent cosmic microwave background observations, ionising background measurements of star-forming galaxies, and low redshifts line-of-sight probes, we propose a new data-motivated parameterisation of the history of the average ionisation fraction. This parameterisation describes a flexible redshift-asymmetric reionisation process in two regimes that is capable of fitting all the current constraints.
We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, in the smooth terrains of the Imhotep region. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.
For the current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs), with their large mirrors and their cameras with fine segmentation of photodetectors, the focusing capability is a relevant issue. The optical system of an IACT has a limited depth of field. Therefore, focusing the telescopes close to the shower maximum in the atmosphere has a significant impact on the data acquisition and analysis. As the distance of the shower maximum to the telescope depends (among others) on the zenith angle, an adjustable focus would be desirable. The fifth Cherenkov telescope of the H.E.S.S. II array is equipped with a focus system which allows to adjust the position of the camera along the optical axis, possibly during data taking. This impact has been studied on gamma-ray Monte Carlo simulations, and the results in terms of gamma-ray trigger rate, energy reconstruction and gamma-ray effective area will be shown.
V1387 Aql (the optical companion to the microquasar GRS1915+105) is a low mass giant. Such star consists of a degenerate, nearly isothermal helium core and a hydrogen rich envelope. Both components are separated by an hydrogen burning shell. The structure of such an object is relatively simple and easy to model. Making use of the observational values of the luminosity and of the radius of V1387 Aql, we determined the mass of this star as equal 0.28+-0.02 solar masses. This determination is relatively precise thanks to high sensitivity of the luminosity of such structure to the mass of the helium core and high sensitivity of its radius to the mass of the envelope. The estimate does not depend on the knowledge of the distance to the system (which is not precisely known). The main source of the uncertainty of my estimate is uncertainty of the effective temperature of V1387 Aql. When the effective temperature will be known more accurately, the mass of V1387 Aql could be determined even more precisely.
The observed rapid cooling of the Cassiopeia A neutron star (Cas A NS) can be interpreted as being triggered by the onset of neutron superfluidity in the core of the star, causing enhanced neutrino emission from neutron Cooper pair breaking and formation (PBF). This provides a unique possibility for probing the neutron condensate in the core. Using consistent neutron star core and crust equation of state and composition, I explore the sensitivity of this interpretation to the phase state of the triplet superfluid condensate. Modeling cooling within an expected range of neutron star masses and envelope compositions, I found that the fast cooling of the Cas A NS can not be explained by the PBF processes in the conventional one-component $^3$P$_2$ condensate with $m_{j}=0$. The best-fit solutions are obtained for the multicomponent superfluid phases listed in Table. The $O_1$ solution yields $M=1.52M_{Sun}$ (carbon envelope with $10^{-15}M_{Sun}$). The $O_2$ solution yield $M=1.47M_{Sun}$ (carbon envelope with $5\times 10^{-15}M_{Sun}$), and the $O_{\pm 3}$ solution $M=1.49M_{Sun}$ (iron envelope).
The spin of a neutron star at birth may be impacted by the asymmetric character of the explosion of its massive progenitor. During the first second after bounce, the spiral mode of the Standing Accretion Shock Instability (SASI) is able to redistribute angular momentum and spin-up a neutron star born from a non-rotating progenitor. Our aim is to assess the robustness of this process. We perform 2D numerical simulations of a simplified setup in cylindrical geometry to investigate the timescale over which the dynamics is dominated by a spiral or a sloshing mode. We observe that the spiral mode prevails only if the ratio of the initial shock radius to the neutron star radius exceeds a critical value. In that regime, both the degree of asymmetry and the average expansion of the shock induced by the spiral mode increase monotonously with this ratio, exceeding the values obtained when a sloshing mode is artificially imposed. With a timescale of the order of 2-3 SASI oscillations, the dynamics of SASI is able to take place fast enough to affect the spin of the neutron star before the explosion. The spin periods deduced from the simulations are compared favorably to analytical estimates evaluated from the measured saturation amplitude of the SASI wave. Incidentally, our study indicates that the analytical prediction of SASI amplitude is not satisfactory yet. Despite the simplicity of our setup, numerical simulations revealed unexpected stochastic variations, including a reversal of the direction of rotation of the accretion shock. Our results show that the spin up of neutron stars by SASI spiral modes is a viable mechanism even though it is not systematic.
We investigated the habitability of the Milky Way, making use of recent observational analysis on the prevalence of Earth-sized planets, in order to estimate where and when potentially habitable star systems may have formed over the course of the Galaxy's history. We were then able to estimate the age distribution of potential intelligent life in our Galaxy using our own evolution and the age of the Sun as a proxy. To do this we created a galactic chemical evolution model and applied the following habitability constraints to the Sun-like (G-type) stars formed in our model: an environment free from life-extinguishing supernovae, a high enough metallicity for Earth-sized planet formation and sufficient time for the evolution of complex life. We determined a galactic habitable zone as the region containing all the potentially habitable star systems in our model. Our galactic habitable zone contains stars formed between 11 and 3.8 billion years ago at radial distances of between 7 and 14 kiloparsecs. We found that most potentially habitable star systems are much older than the Sun and located farther from the galactic centre. By comparing the ages of these systems we estimated that 77% of potentially habitable star systems are on average 3.13 billion years older than the Sun. This suggests that any intelligent life in the Galaxy is likely to be incredibly more advanced than we are assuming that they have evolved under similar timescales than we have. Implications and limitations of our study are discussed.
Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars at $z\sim 2.3-3$ within a $11.5'\times13.5'$ area of the COSMOS field ($\sim 1200\,\mathrm{deg}^2$ projected area density or $\sim 2.4\,h^{-1}\,\mathrm{Mpc}$ mean transverse separation), we reconstruct a 3D tomographic map of the foreground Ly$\alpha$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $\sigma_{3d}\approx3.5\,h^{-1}\,\mathrm{Mpc}$ comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the tomographic map. We find an extended IGM overdensity with deep absorption troughs at $z=2.45$ associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be $M_{\rm dm} (z=2.45) = (9\pm4)\times 10^{13}\,h^{-1}\,\mathrm{M_\odot}$, and argue based on this mass and absorption strength that it will form at least one $z\sim0$ galaxy cluster with $M(z=0) = (3\pm 2) \times 10^{14}\,h^{-1}\mathrm{M_\odot}$, although its elongated nature suggests that it will likely collapse into two separate clusters. We also point out a compact overdensity of six MOSDEF galaxies at $z=2.30$ within a $r\sim 1\,h^{-1}\,\mathrm{Mpc}$ radius and $\Delta z\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$\alpha$ forest tomography on larger volumes to study galaxy properties as a function of environment, as well as revealing the large-scale IGM overdensities associated with protoclusters and other features of large-scale structure.
To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.
We present the discovery of the highest velocity CIV broad absorption line to date in the z=2.47 quasar SDSS J023011.28+005913.6, hereafter J0230. In comparing the public DR7 and DR9 spectra of J0230, we discovered an emerging broad absorption trough outflowing at~60,000 km/s. In pursuing follow up observations we discovered a second emergent broad absorption trough outflowing at ~40,000 km/s. We collected seven spectral epochs of J0230 that demonstrate emergent and rapidly (~10 days in the rest-frame) varying broad absorption. We investigate two possible scenarios that could cause these rapid changes: bulk motion and ionization variability. Given our multi-epoch data, a transverse motion scenario would likely be a flow-tube feature travelling across the emitting region at 8,000 < v (km/s) < 56,000. If ionization variability is the cause for the changes observed, the absorber either has n_e >=1540 cm^-3 and is at r_{eq} >= 1.37 kpc, or is at r < 1.37 kpc with no constraint on the density.
Aims. We study the emergence of a non-twisted flux tube from the solar interior into the solar atmosphere. We investigate whether the length of the buoyant part of the flux tube (i.e. {\lambda}) affects the emergence of the field and the dynamics of the evolving magnetic flux system. Methods. We perform three-dimensional (3D), time-dependent, resistive, compressible MHD simulations using the Lare3D code. Results. We find that there are considerable differences in the dynamics of the emergence of a magnetic flux tube when {\lambda} is varied. In the solar interior, for larger values of {\lambda}, the rising magnetic field emerges faster and expands more due to its lower magnetic tension. As a result, its field strength decreases and its emergence above the photosphere occurs later than in the smaller {\lambda} case. However, in both cases, the emerging field at the photosphere becomes unstable in two places, forming two magnetic bipoles that interact dynamically during the evolution of the system. Most of the dynamic phenomena occur at the current layer, which is formed at the interface between the interacting bipoles. We find the formation and ejection of plasmoids, the onset of successive jets from the interface, and the impulsive heating of the plasma in the solar atmosphere. We discuss the triggering mechanism of the jets and the atmospheric response to the emergence of magnetic flux in the two cases.
PSR B1259-63/LS 2883 is a gamma-ray binary system containing a radio pulsar in a highly elliptical ~3.4-year orbit around a Be star. In its 2010 periastron passage, multiwavelength emission from radio to TeV was observed, as well as an unexpected GeV flare measured by the Fermi Large Area Telescope (LAT). Here, we report the results of LAT monitoring of PSR B1259-63 during its most recent 2014 periastron passage. We compare the gamma-ray behavior in this periastron with the former in 2010 and find that PSR B1259-63 shows a recurrent GeV flare. The similarities and differences in the phenomenology of both periastron passages are discussed.
The shape of the probability distribution function (PDF) of molecular clouds is an important ingredient for modern theories of star formation and turbulence. Recently, several studies have pointed out observational difficulties with constraining the low column density (i.e. Av <1) PDF using dust tracers. In order to constrain the shape and properties of the low column density probability distribution function, we investigate the PDF of multiphase atomic gas in the Perseus molecular cloud using opacity-corrected GALFA-HI data and compare the PDF shape and properties to the total gas PDF and the N(H2) PDF. We find that the shape of the PDF in the atomic medium of Perseus is well described by a lognormal distribution, and not by a power-law or bimodal distribution. The peak of the atomic gas PDF in and around Perseus lies at the HI-H2 transition column density for this cloud, past which the N(H2) PDF takes on a powerlaw form. We find that the PDF of the atomic gas is narrow and at column densities larger than the HI-H2 transition the HI rapidly depletes, suggesting that the HI PDF may be used to find the HI-H2 transition column density. We also calculate the sonic Mach number of the atomic gas by using HI absorption line data, which yields a median value of Ms=4.0 for the CNM, while the HI emission PDF, which traces both the WNM and CNM, has a width more consistent with transonic turbulence.
The High Energy Stereoscopic System (H.E.S.S.) is an array of five Imaging Atmospheric Cherenkov Telescopes (IACTs) designed to detect cosmogenic gamma-rays with very high energies. Originally consisting of just four identical IACTs (CT1-4) with an effective mirror diameter of 12$\,$m each, it was expanded with a fifth IACT (CT5) with a mirror diameter of 28$\,$m in 2012. Being the largest IACT worldwide, CT5 allows to lower the energy threshold of H.E.S.S., making the array sensitive at energies where space-based detectors run out of statistics. Events can be analysed either monoscopically (i.e. using only information of CT5) or stereoscopically (requiring at least two triggered telescopes per event). To achieve a good performance, a sophisticated event reconstruction and analysis framework is indispensable. This is particularly important for H.E.S.S. since it is now the first IACT array that consists of different telescope types. An advanced reconstruction method is based on a semi-analytical model of electromagnetic particle showers in the atmosphere (model analysis). The properties of the primary particle are reconstructed by comparing the image recorded by each triggered telescope with the Cherenkov emission from the shower model using a log-likelihood maximisation. Due to its performance, this method has become one of the standard analysis techniques applied to CT1-4 data. Now it has been modified for use with the five-telescope array. We present the adapted model analysis and its performance in both monoscopic and stereoscopic analysis mode.
The High Energy Stereoscopic System (H.E.S.S.) phase I instrument was an array of four $100\,\mathrm{m}^2$ mirror area Imaging Atmospheric Cherenkov Telescopes (IACTs) that has very successfully mapped the sky at photon energies above $\sim 100\,$GeV. Recently, a $600\,\mathrm{m}^2$ telescope was added to the centre of the existing array, which can be operated either in standalone mode or jointly with the four smaller telescopes. The large telescope lowers the energy threshold for gamma-ray observations to several tens of GeV, making the array sensitive at energies where the Fermi-LAT instrument runs out of statistics. At the same time, the new telescope makes the H.E.S.S. phase II instrument. This is the first hybrid IACT array, as it operates telescopes of different size (and hence different trigger rates) and different field of view. In this contribution we present results of H.E.S.S. phase II observations of the Crab Nebula, compare them to earlier observations, and evaluate the performance of the new instrument with Monte Carlo simulations.
Up to ages of ~100 Myr, massive clusters are still swamped in large amounts of gas and dust, with considerable and uneven levels of extinction. At the same time, large grains (ices?) produced by type II supernovae profoundly alter the interstellar medium (ISM), thus resulting in extinction properties very different from those of the diffuse ISM. To obtain physically meaningful parameters of stars, from basic luminosities and effective temperatures to masses and ages, we must understand and measure the local extinction law. This problem affects all the massive young clusters discussed in his volume. We have developed a powerful method to unambiguously determine the extinction law in an uniform way across a cluster field, using multi-band photometry of red giant stars belonging to the red clump (RC). In the Large Magellanic Cloud, with about 20 RC stars per arcmin^2, we can easily derive a solid and self-consistent absolute extinction curve over the entire wavelength range of the photometry. Here, we present the extinction law of the Tarantula nebula (30 Dor) based on thousands of stars observed as part of the Hubble Tarantula Treasury Project.
We consider a synchrotron radiation from a charged particle moving in a bound orbit around a weakly magnetized Schwarzschild black hole (a static black hole immersed into a constant uniform magnetic field) in its equatorial plane, perpendicular to the magnetic field. In particular, we study the case when the Lorentz force acting on the charged particle is directed outward from the black hole. The particle is initially moving in a nongeodesic bound orbit which due to synchrotron radiation decays to a nongeodesic circular orbit. We study this transition and calculate the radiated power and energy loss of the particle.
It was recently shown that gravitons with a very small mass should have formed a Bose-Einstein condensate in the very early Universe, whose density and quantum potential can account for the dark matter and dark energy in the Universe respectively. Here we show that the condensation can also naturally explain the observed large scale homogeneity and isotropy of the Universe. Furthermore gravitons continue to fall into their ground state within the condensate at every epoch, accounting for the observed flatness of space at cosmological distances scales. Finally, we argue that the density perturbations due to quantum fluctuations within the condensate give rise to a scale invariant spectrum. This therefore provides a viable alternative to inflation, which is not associated with the well-known problems associated with the latter.
We present a self-gravitating Skyrmion, an analytic and globally regular solution of the Einstein- Skyrme system in presence of a cosmological constant with winding number w = 1. The static spacetime metric is the direct product R x S3 and the Skyrmion is the self-gravitating generalization of the static hedgehog solution of Manton and Ruback with unit topological charge. This solution can be promoted to a dynamical one in which the spacetime is a cosmology of the Bianchi type-IX with time-dependent scale and squashing coefficients. Remarkably, the Skyrme equations are still identically satisfied for all values of these parameters. Thus, the complete set of field equations for the Einstein-Skyrme-Lambda system in the topological sector reduces to a pair of coupled, autonomous, nonlinear differential equations for the scale factor and a squashing coefficient. These equations admit analytic bouncing cosmological solutions in which the universe contracts to a mini- mum non-vanishing size, and then expands. A non-trivial byproduct of this solution is that a minor modification of the construction gives rise to a family of stationary, regular configurations in General Relativity with negative cosmological constant supported by an SU(2) nonlinear sigma model. These solutions represent traversable wormholes with NUT parameter in which the only "exotic matter" required for their construction is a negative cosmological constant.
The description of the universe evolving in time according to general relativity is given in comparison with the quantum description of the same universe in terms of quasiclassical wave functions. The spacetime geometry is determined by the Robertson-Walker metric. It is shown that the main equation of the quantum geometrodynamics is reduced to the non-linear Hamilton-Jacobi equation. Its non-linearity is caused by a new source of the gravitational field, which has a purely quantum dynamical nature, and is additional to ordinary matter sources. In quasiclassical approximation, the non-linear equation of motion is linearized and reduces to the Friedmann equation with the additional quantum source of gravity (or anti-gravity) in the form of the stiff Zel'dovich matter. The semi-classical wave functions of the universe, in which different types of matter-energies dominate, are obtained. As examples, the cases of the domination of radiation, barotropic fluid, or new quantum matter-energy are discussed. The probability of the transition from the quantum state, where radiation dominates into the state, in which barotropic fluid in the form of dust is dominant, is calculated. In the era of matter-radiation equality, this probability has the same order of magnitude as the matter density contrast.
A phenomenological turbulence model for kinetic Alfv\'en waves in a magnetized collisionless plasma, able to reproduce the non-universal power-law spectra observed at the sub-ion scales in the solar wind and the terrestrial magnetosphere, is presented. Nonlocal interactions are retained, and critical balance, characteristic of a strong turbulence regime, establishes dynamically as the cascade proceeds. The process of temperature homogenization along distorted magnetic field lines, induced by Landau damping, affects the turbulence transfer time and results in a steepening of the sub-ion power-law spectrum of critically-balanced turbulence, whose exponent is in particular sensitive to the ratio between the Alfv\'en wave period and the nonlinear timescale.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)
We summarise some of the recent progress in understanding the formation and evolution of globular clusters (GCs) in the context of galaxy formation and evolution. It is discussed that an end-to-end model for GC formation and evolution should capture four different phases: (1) star and cluster formation in the high-pressure interstellar medium of high-redshift galaxies, (2) cluster disruption by tidal shocks in the gas-rich host galaxy disc, (3) cluster migration into the galaxy halo, and (4) the final evaporation-dominated evolution of GCs until the present day. Previous models have mainly focussed on phase 4. We present and discuss a simple model that includes each of these four steps - its key difference with respect to previous work is the simultaneous addition of the high-redshift formation and early evolution of young GCs, as well as their migration into galaxy haloes. The new model provides an excellent match to the observed GC mass spectrum and specific frequency, as well as the relations of GCs to the host dark matter halo mass and supermassive black hole mass. These results show (1) that the properties of present-day GCs are reproduced by assuming that they are the natural outcome of regular high-redshift star formation (i.e. they form according to same physical processes that govern massive cluster formation in the local Universe), and (2) that models only including GC evaporation strongly underestimate their integrated mass loss over a Hubble time.
K2 space observations recently found that three super-Earths transit the nearby M dwarf K2-3. The apparent brightness and the small physical radius of their host star rank these planets amongst the most favourable for follow-up characterisations. The outer planet orbits close to the inner edge of the habitable zone and might become one of the first exoplanets searched for biomarkers using transmission spectroscopy. We used the HARPS velocimeter to measure the mass of the planets. The mass of planet $b$ is $8.4\pm2.1$ M$_\oplus$, while our determination of those planets $c$ and $d$ are affected by the stellar activity. With a density of $4.32^{+2.0}_{-0.76}$ $\mathrm{g\;cm^{-3}}$, planet $b$ is probably mostly rocky, but it could contain up to 50% water.
Mergers of gas-rich galaxies are key events in the hierarchical built-up of cosmic structures, and can lead to the formation of massive black hole binaries. By means of high-resolution hydrodynamical simulations we consider the late stages of a gas-rich major merger, detailing the dynamics of two circumnuclear discs, and of the hosted massive black holes during their pairing phase. During the merger gas clumps with masses of a fraction of the black hole mass form because of fragmentation. Such high-density gas is very effective in forming stars, and the most massive clumps can substantially perturb the black hole orbits. After $\sim 10$ Myr from the start of the merger a gravitationally bound black hole binary forms at a separation of a few parsecs, and soon after, the separation falls below our resolution limit of $0.39$ pc. At the time of binary formation the original discs are almost completely disrupted because of SNa feedback, while on pc scales the residual gas settles in a circumbinary disc with mass $\sim 10^5 M_\odot$. We also test that binary dynamics is robust against the details of the SNa feedback employed in the simulations, while gas dynamics is not. We finally highlight the importance of the SNa time-scale on our results.
We present a comprehensive catalog and analysis of broad-band afterglow observations for 103 short-duration gamma-ray bursts (GRBs), comprised of all short GRBs from November 2004 to March 2015 with prompt follow-up observations in the X-ray, optical, near-infrared and/or radio bands. These afterglow observations have uncovered 71 X-ray detections, 30 optical/NIR detections, and 4 radio detections. Employing the standard afterglow synchrotron model, we perform joint probability analyses for a subset of 38 short GRBs with well-sampled light curves to infer the burst isotropic-equivalent energies and circumburst densities. For this subset, we find median isotropic-equivalent gamma-ray and kinetic energies of E_gamma,iso~2x10^51 erg, and E_K,iso~(1-3)x10^51 erg, respectively, depending on the values of the model input parameters. We further find that short GRBs occur in low-density environments, with a median density of n~(3-15)x10^-3 cm^-3, and that ~80-95% of bursts have densities of less than 1 cm^-3. We investigate trends between the circumburst densities and host galaxy properties, and find that events located at large projected offsets of >10 effective radii from their hosts exhibit particularly low densities of n<10^-4 cm^-3, consistent with an IGM-like environment. Using late-time afterglow data for 11 events, we find a median jet opening angle of theta_jet=16+/-10 deg. We also calculate a median beaming factor of f_b~0.04, leading to a beaming-corrected total energy release of E_true~1.6x10^50 erg. Furthermore, we calculate a beaming-corrected event rate of R_true=270 (+1580,-180) Gpc^-3 yr^-1, or ~8 (+47,-5) yr^-1 within a 200 Mpc volume, the Advanced LIGO/Virgo typical detection distance for NS-NS binaries.
Data selection and the determination of systematic uncertainties in the spectroscopic measurements of neutron star radii from thermonuclear X-ray bursts have been the subject of numerous recent studies. In one approach, the uncertainties and outliers were determined by a data-driven Bayesian mixture model, whereas in a second approach, data selection was performed by requiring that the observations follow theoretical expectations. We show here that, due to inherent limitations in the data, the theoretically expected trends are not discernible in the majority of X-ray bursts even if they are present. Therefore, the proposed theoretical selection criteria are not practical with the current data for distinguishing clean data sets from outliers. Furthermore, when the data limitations are not taken into account, the theoretically motivated approach selects a small subset of bursts with properties that are in fact inconsistent with the underlying assumptions of the method. We conclude that the data-driven selection methods do not suffer from the limitations of this theoretically motivated one.
Observations using the Fermi Large Area Telescope (Fermi-LAT) have found a significant gamma-ray excess surrounding the center of the Milky Way (GC). One possible interpretation of this excess invokes gamma-ray emission from an undiscovered population of either young or recycled pulsars densely clustered throughout the inner kiloparsec of the Milky Way. While these systems, by construction, have individual fluxes that lie below the point source sensitivity of the Fermi-LAT, they may already be observed in multiwavelength observations. Notably the Australia Telescope National Facility (ATNF) catalog of radio pulsars includes 270 sources observed in the inner 10 degrees around the GC. We calculate the gamma-ray emission observed from these 270 sources and obtain three key results: (1) point source searches in the GC region produce a plethora of highly significant gamma-ray "hotspots", compared to searches far from the Galactic plane, (2) there is no statistical correlation between the positions of these gamma-ray hotspots and the locations of ATNF pulsars, and (3) the spectrum of the most statistically significant gamma-ray hotspots is substantially softer than the spectrum of the GC gamma-ray excess. These results place strong constraints on models where young pulsars produce the majority of the gamma-ray excess, and disfavor some models where millisecond pulsars produce the gamma-ray excess.
We present a wide area (~ 8 x 8 kpc), sensitive map of CO (2-1) emission around the nearby starburst galaxy M82. Molecular gas extends far beyond the stellar disk, including emission associated with the well-known outflow as far as 3 kpc from M82's midplane. Kinematic signatures of the outflow are visible in both the CO and HI emission: both tracers show a minor axis velocity gradient and together they show double peaked profiles, consistent with a hot outflow bounded by a cone made of a mix of atomic and molecular gas. Combining our CO and HI data with observations of the dust continuum, we study the changing properties of the cold outflow as it leaves the disk. While H_2 dominates the ISM near the disk, the dominant phase of the cool medium changes as it leaves the galaxy and becomes mostly atomic after about a kpc. Several arguments suggest that regardless of phase, the mass in the cold outflow does not make it far from the disk; the mass flux through surfaces above the disk appears to decline with a projected scale length of ~ 1-2 kpc. The cool material must also end up distributed over a much wider angle than the hot outflow based on the nearly circular isophotes of dust and CO at low intensity and the declining rotation velocities as a function of height from the plane. The minor axis of M82 appears so striking at many wavelengths because the interface between the hot wind cavity and the cool gas produces Halpha, hot dust, PAH emission, and scattered UV light. We also show the level at which a face-on version of M82 would be detectable as an outflow based on unresolved spectroscopy. Finally, we consider multiple constraints on the CO-to-H$_2$ conversion factor, which must change across the galaxy but appears to be only a factor of ~ 2 lower than the Galactic value in the outflow.
The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.
Small-scale dark matter structure within the Milky Way is expected to affect pulsar timing. The change in gravitational potential induced by a dark matter halo passing near the line of sight to a pulsar would produce a varying delay in the light travel time of photons from the pulsar. Individual transits produce an effect that would either be too rare or too weak to be detected in 30-year pulsar observations. However, a population of dark matter subhalos would be expected to produce a detectable effect on the measured properties of pulsars if the subhalos constitute a significant fraction of the total halo mass. The effect is to increase the dispersion of measured period derivatives across the pulsar population. By statistical analysis of the ATNF pulsar catalogue, we place an upper limit on this dispersion of $\log \sigma_{\dot{P}} \leq -17.05$. We use this to place strong upper limits on the number density of ultracompact minihalos within the Milky Way. These limits are completely independent of the particle nature of dark matter.
Ultracompact Minihalos (UCMHs) have been proposed as a type of dark matter sub-structure seeded by large-amplitude primordial perturbations and topological defects. UCMHs are expected to survive to the present era, allowing constraints to be placed on their cosmic abundance using observations within our own Galaxy. Constraints on their number density can be linked to conditions in the early universe that impact structure formation, such as increased primordial power on small scales, generic weak non-Gaussianity, and the presence of cosmic strings. We use new constraints on the abundance of UCMHs from pulsar timing to place generalised limits on the parameters of each of these cosmological scenarios. At some scales, the limits are the strongest to date, exceeding those from dark matter annihilation. Our new limits have the added advantage of being independent of the particle nature of dark matter, as they are based only on gravitational effects.
We present Sunyaev-Zel'dovich (SZ) effect measurements from wide-field images towards the galaxy cluster RX J1347.5-1145 obtained from the Caltech Submillimeter Observatory with the Multiwavelength Submillimeter Inductance Camera (MUSIC) at 147, 213, 281, and 337 GHz and with Bolocam at 140 GHz. As part of our analysis, we have used higher frequency data from Herschel-SPIRE and previously published lower frequency radio data to subtract the signal from the brightest dusty star-forming galaxies behind RX J1347.5-1145 and from the AGN in RX J1347.5-1145's BCG. Using these five-band SZ effect images, combined with previously published X-ray spectroscopic measurements of the temperature of the intra-cluster medium (ICM) from Chandra, we constrain the ICM optical depth to be $\tau_e = 2.73^{+0.38}_{-0.39} \times 10^{-3}$ and the ICM line of sight peculiar velocity to be $v_{pec} = -1260^{+760}_{-530}$ km s$^{-1}$. The errors for both quantities are limited by measurement noise rather than calibration uncertainties or astrophysical contamination, and significant improvements are possible with deeper observations. Our best-fit velocity is in mild tension with the two previously published SZ effect measurements for RX J1347.5-1145, although some or all of the tension may be because those measurements sample the cluster differently.
When they first appear in the HR diagram, young stars rotate at a mere 10\% of their break-up velocity. They must have lost most of the angular momentum initially contained in the parental cloud, the so-called angular momentum problem. We investigate here a new mechanism by which large amounts of angular momentum might be shed from young stellar systems, thus yielding slowly rotating young stars. Assuming that planets promptly form in circumstellar disks and rapidly migrate close to the central star, we investigate how the tidal and magnetic interactions between the protostar, its close-in planet(s), and the inner circumstellar disk can efficiently remove angular momentum from the central object. We find that neither the tidal torque nor the variety of magnetic torques acting between the star and the embedded planet are able to counteract the spin up torques due to accretion and contraction. Indeed, the former are orders of magnitude weaker than the latter beyond the corotation radius and are thus unable to prevent the young star from spinning up. We conclude that star-planet interaction in the early phases of stellar evolution does not appear as a viable alternative to magnetic star-disk coupling to understand the origin of the low angular momentum content of young stars.
Broad Fe K emission lines have been widely observed in the X-ray spectra of black hole systems, and in neutron star systems as well. The intrinsically narrow Fe K fluorescent line is generally believed to be part of the reflection spectrum originating in an illuminated accretion disk, and broadened by strong relativistic effects. However, the nature of the lines in neutron star LMXBs has been under debate. We therefore obtained the longest, high-resolution X-ray spectrum of a neutron star LMXB to date with a 300 ks Chandra HETGS observation of Serpens X-1. The observation was taken under the "continuous clocking" mode and thus free of photon pile-up effects. We carry out a systematic analysis and find that the blurred reflection model fits the Fe line of Serpens X-1 significantly better than a broad Gaussian component does, implying that the relativistic reflection scenario is much preferred. Chandra HETGS also provides highest spectral resolution view of the Fe K region and we find no strong evidence for additional narrow lines.
IceCube, a cubic-kilometer sized neutrino detector at the Geographic South Pole, has recently discovered a diffuse all-flavor flux of astrophysical neutrinos. However, the corresponding astrophysical sources have not yet been identified in current IceCube analyses. We present a method to interpret the results of a recently published angular correlation analysis in IceCube searching for spatial clustering of muon neutrino events in terms of astrophysical models (given by an arbitrary source count distribution). We exemplarily show the resulting limits on the parameters of a class of source count distributions motivated by Fermi-LAT observations of resolved blazars.
What main mechanisms set the star formation rate (SFR) of galaxies? This PhD thesis is a quest into the influences of gas and Active Galactic Nuclei (AGNs) on the SFR, with particular focus on massive galaxies at z~2. First, a new code if presented; SImulator of GAlaxy Millimeter/submillimeter Emission (S\'IGAME) which can predict the atomic/molecular line emission from galaxies. By post-processing the outputs of cosmological simulations of galaxy formation with sub-grid physics recipes, S\'IGAME divides the Interstellar Medium (ISM) into different gas phases and derives density and temperature structure, with locally resolved radiation fields. This method is used to predict the strengths of CO rotational transitions as well as the [CII] emission line in normal star-forming galaxies at z~2. A CO ladder close to that of our own Galaxy is found, but with CO-H2 conversion factors about 3 times smaller. For a set of 7 simulated galaxies at z~2, the relation between [CII] luminosity and SFR displays a slope significantly steeper than that found for observed galaxies at z<0.5. A corresponding relation on kpc-scales is established for the first time theoretically. Finally, a separate study entails the analysis of CHANDRA CDF-S X-ray data with the aim of uncovering AGNs among massive galaxies at z~2. It is found that about every fifth massive galaxy, quenched or not, contains an X-ray luminous AGN. Interestingly, an even higher fraction of low-luminosity AGNs emerges in the X-ray undetected galaxies when performing a stacking analysis, and preferentially in the quenched ones, lending support to the importance of AGNs in impeding star formation during galaxy evolution.
In this paper, we study the homogeneity of the GRB distribution using a subsample of the Greiner GRB catalogue, which contains 314 objects with redshift $0<z<2.5$ (244 of them discovered by the Swift GRB Mission). We try to reconcile the dilemma between the new observations and the current theory of structure formation and growth. To test the results against the possible biases in redshift determination and the incompleteness of the Greiner sample, we also apply our analysis to the 244 GRBs discovered by Swift and the subsample presented by the Swift Gamma-Ray Burst Host Galaxy Legacy Survey (SHOALS). The real space two-point correlation function (2PCF) of GRBs, $\xi(r),$ is calculated using a Landy-Szalay estimator. We perform a standard least-$\chi^2$ fit to the measured 2PCFs of GRBs. We use the best-fit 2PCF to deduce a recently defined homogeneity scale. The homogeneity scale, $R_H$, is defined as the comoving radius of the sphere inside which the number of GRBs $N(<r)$ is proportional to $r^3$ within $1\%$, or equivalently above which the correlation dimension of the sample $D_2$ is within $1\%$ of $D_2=3$. For the Swift subsample of 244 GRBs, the correlation length and slope are $r_0= 387.51 \pm 132.75~h^{-1}$Mpc and $\gamma = 1.57\pm 0.65$ (at $1\sigma$ confidence level). The corresponding scale for a homogeneous distribution of GRBs is $r\geq 7,700~h^{-1}$Mpc. The results help to alleviate the tension between the new discovery of the excess clustering of GRBs and the cosmological principle of large-scale homogeneity. It implies that very massive structures in the relatively local Universe do not necessarily violate the cosmological principle and could conceivably be present.
We report on the identification of a new $\gamma$-ray-emitting narrow-line Seyfert 1 (NLS1) galaxy, SDSS J122222.55+041315.7, which increases the number of known objects of this remarkable but rare type of active galactic nuclei (AGN) to seven. Its optical spectrum, obtained in the Sloan Digital Sky Survey-Baryon Oscillation Spectroscopic Survey, reveals a broad H $\beta$ emission line with a width (FWHM) of 1734$\pm$104 km s$^{-1}$. This, along with strong optical Fe II multiplets [$R_{4570}=0.9$] and a weak [O III] $\lambda 5007$ emission line, makes the object a typical NLS1. On the other hand, the source exhibits a high radio brightness temperature, rapid infrared variability, and a flat X-ray spectrum extending up to $\sim$200 keV. It is associated with a luminous $\gamma$-ray source detected significantly with {\it Fermi}/LAT. Correlated variability with other wavebands has not yet been tested. The spectral energy distribution can be well modelled by a one-zone leptonic jet model. This new member is by far the most distant $\gamma$-ray-emitting NLS1, at a redshift of $z=0.966$.
Based on fundamental particle physics processes like the production and subsequent decay of pions in interactions of high-energy particles, close connections exist between the acceleration sites of high-energy cosmic rays and the emission of high-energy gamma rays and high-energy neutrinos. In most cases these connections provide both spatial and temporal correlations of the different emitted particles. The combination of the complementary information provided by these messengers allows to lift ambiguities in the interpretation of the data and enables novel and highly sensitive analyses. In this contribution the H.E.S.S. multi-messenger program is introduced and described. The current core of this newly installed program is the combination of high-energy neutrinos and high-energy gamma rays. The search for gamma-ray emission following gravitational wave triggers is also discussed. Furthermore, the existing program for following triggers in the electromagnetic regime was extended by the search for gamma-ray emission from Fast Radio Bursts (FRBs). An overview over current and planned analyses is given and recent results are presented.
Microquasars, Galactic binary systems showing extended and variable radio emission, are potential gamma-ray emitters. Indications of gamma-ray transient episodes have been reported in at least two systems, Cyg X-1 and Cyg X-3. The identification of additional gamma-ray emitting microquasars is key for a better understanding of these systems. Very-high-energy gamma-ray emission from microquasars has been predicted to happen during periods of transient outbursts potentially connected with the formation of a jet-like outflow. Contemporaneous observations using the H.E.S.S. telescope array and the RXTE satellite were obtained on three microquasars: GRS 1915+105, Circinus X-1 and V4641 Sgr with the aim of detecting a broadband flaring event in the very-high-energy gamma-ray and X-ray bands. We report here on the analysis of these data for each system, including a detailed X-ray analysis assessing the location of the sources in a hardness-intensity diagram during the observations. Finally we discuss the derived upper limits on their very-high-energy gamma-ray flux.
We present new ALMA continuum observations at 336 GHz of two transition disks, SR 21 and HD 135344B. In combination with previous ALMA observations from Cycle 0 at 689 GHz, we compare the visibility profiles at the two frequencies and calculate the spectral index ($\alpha_{\rm{mm}}$). The observations of SR 21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two wavelengths. Notable radial variations of the spectral index are also detected for SR 21 with values of $\alpha_{\rm{mm}}\sim3.8-4.2$ in the inner region ($r\lesssim35$ AU) and $\alpha_{\rm{mm}}\sim2.6-3.0$ outside. An axisymmetric ring ("ring model") or a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure ("vortex model"), is assumed for fitting the disk morphology. For SR 21, the ring model better fits the emission at 336 GHz, conversely the vortex model better fits the 689 GHz emission. For HD 135344B, neither a significant shift in the null of the visibilities nor radial variations of $\alpha_{\rm{mm}}$ are detected. Furthermore, for HD 135344B, the vortex model fits both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations of $\alpha_{\rm{mm}}$ remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR 21 have grown to millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear evidence of radial trapping in HD 135344B, although it cannot be excluded either.
In this paper we show that the phase mixing of continuum Alfv\'{e}n waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfv\'{e}n and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfv\'{e}n waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfv\'{e}n/slow frequency across the magnetic surfaces and is slower than the Alfv\'{e}n/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.
The total population of Wolf-Rayet (WR) stars in the Galaxy is predicted by models to be as many as $\sim$6000 stars, and yet the number of catalogued WR stars as a result of optical surveys was far lower than this ($\sim$200) at the turn of this century. When beginning our WR searches using infrared techniques it was not clear whether WR number predictions were too optimistic or whether there was more hidden behind interstellar and circumstellar extinction. During the last decade we pioneered a technique of exploiting the near- and mid-infrared continuum colours for individual point sources provided by large-format surveys of the Galaxy, including 2MASS and Spitzer/GLIMPSE, to pierce through the dust and reveal newly discovered WR stars throughout the Galactic Plane. The key item to the colour discrimination is via the characteristic infrared spectral index produced by the strong winds of the WR stars, combined with dust extinction, which place WR stars in a relatively depopulated area of infrared colour-colour diagrams. The use of the Spitzer/GLIMPSE 8$\mu$m and, more recently, WISE 22$\mu$m fluxes together with cross-referencing with X-ray measurements in selected Galactic regions have enabled improved candidate lists that increased our confirmation success rate, achieved via follow-up infrared and optical spectroscopy. To date a total of 102 new WR stars have been found with many more candidates still available for follow-up. This constitutes an addition of $\sim$16\% to the current inventory of 642 Galactic WR stars. In this talk we review our methods and provide some new results and a preliminary review of their stellar and interstellar medium environments. We provide a roadmap for the future of this search, including statistical modeling, and what we can add to star formation and high mass star evolution studies.
We consider renormalization of the Adhesion Model for cosmic structure formation. This is a simple model that shares many relevant features of recent approaches which add effective viscosity and noise terms to the fluid equations of Cold Dark Matter, offering itself as a pedagogical playground to study the removal of the cutoff dependence from loop integrals. We show in this context that if the viscosity and noise terms are treated as perturbative corrections to the standard eulerian perturbation theory, as is done for example in the Effective Field Theory of Large Scale Structure (EFToLSS) approach, they are necessarily non-local in time. To ensure Galilean Invariance higher order vertices related to the viscosity and the noise must be added. We explicitly show at one-loop that these terms act as counter terms for vertex diagrams, while the Ward Identities ensure that the non-local theory can be renormalized consistently. A local-in-time theory is renormalizable if the viscosity is included in the linear propagator. Compared to the EFToLSS approach, the local-in-time theory requires less free parameters for its renormalization.
The Cherenkov Telescope Array (CTA) is an international next-generation ground-based gamma-ray observatory. CTA will be implemented as southern and northern hemisphere arrays of tens of small, medium and large-sized imaging Cherenkov telescopes with the goal of improving the sensitivity over the current-generation experiments by an order of magnitude. CTA will provide energy coverage from ~20 GeV to more than 300 TeV. The Schwarzschild-Couder (SC) medium size (9.5m) telescopes will feature a novel aplanatic two-mirror optical design capable of accommodating a wide field-of-view with significantly improved angular resolution as compared to the traditional Davies-Cotton optical design. A full-scale prototype SC medium size telescope structure has been designed and will be constructed at the Fred Lawrence Whipple Observatory in southern Arizona during the fall of 2015. concentrate on the novel features of the design.
An extended observation campaign of the gamma-ray binary system PSR B1259$-$63 has been conducted with the H.E.S.S. (High Energy Stereoscopic System) II 5-telescope array during the system's periastron passage in 2014. We report on the outcome of this campaign, which consists of more than 85 h of data covering both pre- and post-periastron orbital phases. The lower energy threshold of the H.E.S.S. II array allows very-high-energy (VHE; $E \gtrsim 100$ GeV) gamma-ray emission from PSR B1259$-$63 to be studied for the first time down to 200 GeV. The new dataset partly overlaps with and extends in phase previous H.E.S.S. campaigns on this source in 2004, 2007 and 2011, allowing for a detailed long-term characterisation of the flux level at VHEs. In addition, the 2014 campaign reported here includes VHE observations during the exact periastron time, $t_{\rm per}$, as well as data taken simultaneously to the gamma-ray flare detected with the Fermi-LAT. Our results will be discussed in a multiwavelength context, thanks to the dense broad-band monitoring campaign conducted on the system during this last periastron passage.
The field of the globular cluster M12 (NGC 6218) was monitored between 1995 and 2009 in a search for variable stars. $BV$ light curves were obtained for 36 periodic or likely periodic variables. Thirty-four of these are new detections. Among the latter we identified 20 proper-motion members of the cluster: six detached or semi-detached eclipsing binaries, five contact binaries, five SX~Phe pulsators, and three yellow stragglers. Two of the eclipsing binaries are located in the turnoff region, one on the lower main sequence and the remaining three among the blue stragglers. Two contact systems are blue stragglers, and the remaining three reside in the turnoff region. In the blue straggler region a total of 103 objects were found, of which 42 are proper motion members of M12, and another four are field stars. Forty-five of the remaining objects are located within two core radii from the center of the cluster, and as such they are likely genuine blue stragglers. We also report the discoveries of a radial color gradient of M12, and the shortest period among contact systems in globular clusters in general
Gamma-ray burst GRB 140430A was detected by the Swift satellite and observed promptly with the imaging polarimeter RINGO3 mounted on the Liverpool Telescope, with observations beginning while the prompt $\gamma$-ray emission was still ongoing. In this paper, we present densely sampled (10-second temporal resolution) early optical light curves in 3 optical bands and limits to the degree of optical polarization. We compare optical, X-ray and gamma-ray properties and present an analysis of the optical emission during a period of high-energy flaring. The complex optical light curve cannot be explained merely with a combination of forward and reverse shock emission from a standard external shock, implying additional contribution of emission from internal shock dissipation. We estimate an upper limit for time averaged optical polarization during the prompt phase to be as low as P < 12% (1$\sigma$). This suggests that the optical flares and early afterglow emission in this GRB are not highly polarized. Alternatively, time averaging could mask the presence of otherwise polarized components of distinct origin at different polarization position angles.
Observational cosmology is entering an era in which high precision will be required in both measurement and data analysis. Accuracy, however, can only be achieved with a thorough understanding of potential sources of contamination from foreground effects. Our primary focus will be on non- Gaussian effects in foregrounds. This issue will be crucial for coming experiments to determine B-mode polarization. We propose a novel method for investigating a data set in terms of skewness and kurtosis in locally defined regions that collectively cover the entire sky. The method is demonstrated on two sky maps: (i) the SMICA map of Cosmic Microwave Background fluctuations provided by the Planck Collaboration and (ii) a version of the Haslam map at 408 MHz that describes synchrotron radiation. We find that skewness and kurtosis can be evaluated in combination to reveal local physical information. In the present case, we demonstrate that the local properties of both maps are predominantly Gaussian. This result was expected for the SMICA map; that it also applies for the Haslam map is surprising. The approach described here has a generality and flexibility that should make it useful in a variety of astrophysical and cosmological contexts.
Time variability of the photon flux is a known feature of active galactic nuclei (AGN) and in particular of blazars. The high frequency peaked BL Lac (HBL) object PKS 2155-304 is one of the brightest sources in the TeV band and has been monitored regularly with different instruments and in particular with the H.E.S.S. experiment above 200 GeV for more than 11 years. These data together with the observations of other instruments and monitoring programs like SMARTS (optical), Swift-XRT/RXTE/XMM-Newton (X-ray) and Fermi-LAT (100 MeV < E < 300 GeV) are used to characterize the variability of this object in the quiescent state over a wide energy range. Variability studies are made by looking at the lognormality of the light curves and at the fractional root mean square (rms) variability Fvar in several energy bands. Lognormality is found in every energy range and the evolution of Fvar with the energy shows a similar increase both in X-rays and in TeV bands.
We report on annual parallax and proper motion observations of H2O masers in S235AB-MIR, which is a massive young stellar object in the Perseus Arm. Using multi-epoch VLBI astrometry we measured a parallax of pi = 0.63 +- 0.03 mas, corresponding to a trigonometric distance of D = 1.56+-0.09 kpc, and source proper motion of ( u alpha cos d , u d) = (0.79 +- 0.12, -2.41 +- 0.14) mas/yr. Water masers trace a jet of diameter 15 au which exhibits a definite radial velocity gradient perpendicular to its axis. 3D maser kinematics were well modelled by a rotating cylinder with physical parameters: v_out = 45+-2 km/s, v_rot = 22+-3 km/s, i = 12+-2 degrees, which are the outflow velocity, tangential rotation velocity and line-of-sight inclination, respectively. One maser feature exhibited steady acceleration which may be related to the jet rotation. During our 15 month VLBI programme there were three `maser burst' events caught `in the act' which were caused by the overlapping of masers along the line of sight.
Tomography of the solar corona can provide cruicial constraints for models of the low corona, unique information on changes in coronal structure and rotation rates, and a valuable boundary condition for models of the heliospheric solar wind. This is the first of a series of three papers which aim to create a set of maps of the coronal density over an extended period (1996-present). The papers will describe the data processing and calibration (this paper), the tomography method (\paperii) and resulting atlas of coronal electron density at a height of 5\Rs\ between years 1996-2014 (\paperiii). This first paper presents a detailed description of data processing and calibration for the Large-Angle and Spectrometric Coronagraph (LASCO) C2 instrument onboard the Solar and Heliospheric Observatory (SOHO) and the COR2 instruments of the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) package aboard the Solar Terrestial Relations Observatory (STEREO) A \& B spacecraft. The methodology includes noise suppression, background subtraction, separation of large dynamic events, conversion of total brightness to K-coronal brightness and simple functions for crosscalibration between C2/LASCO and COR2/SECCHI. Comparison of the brightness of stars between LASCO C2 total and polarized brightness (\pB) observations provide in-flight calibration factors for the \pB\ observations, resulting in considerable improved agreement between C2 and COR2 A, and elimination of curious artifacts in the C2 \pB\ images. The crosscalibration between LASCO C2 and the STEREO coronagraphs allows, for the first time, the potential use of multi-spacecraft coronagraph data for tomography and for CME analysis.
The search for habitable exoplanets in the Universe is actively ongoing in the field of astronomy. The biggest future milestone is to determine whether life exists on such habitable exoplanets. In that context, oxygen in the atmosphere has been considered strong evidence for the presence of photosynthetic organisms. In this paper, we show that a previously unconsidered photochemical mechanism by titanium(IV) oxide (titania) can produce abiotic oxygen from liquid water under near ultraviolet (NUV) lights on the surface of exoplanets. Titania works as a photocatalyst to dissociate liquid water in this process. This mechanism offers a different source of a possibility of abiotic oxygen in atmospheres of exoplanets from previously considered photodissociation of water vapor in upper atmospheres by extreme ultraviolet (XUV) light. Our order-of-magnitude estimation shows that possible amounts of oxygen produced by this abiotic mechanism can be comparable with or even more than that in the atmosphere of the current Earth, depending on the amount of active surface area for this mechanism. We conclude that titania may act as a potential source of false signs of life on habitable exoplanets.
We report a development of a multi-color simultaneous camera for the 188cm telescope at Okayama Astrophysical Observatory in Japan. The instrument, named MuSCAT, has a capability of 3-color simultaneous imaging in optical wavelength where CCDs are sensitive. MuSCAT is equipped with three 1024x1024 pixel CCDs, which can be controlled independently. The three CCDs detect lights in $g'_2$ (400--550 nm), $r'_2$ (550--700 nm), and $z_{s,2}$ (820--920 nm) bands using Astrodon Photometrics Generation 2 Sloan filters. The field of view of MuSCAT is 6.1x6.1 arcmin$^2$ with the pixel scale of 0.358 arcsec per pixel. The principal purpose of MuSCAT is to perform high precision multi-color transit photometry. For the purpose, MuSCAT has a capability of self autoguiding which enables to fix positions of stellar images within ~1 pix. We demonstrate relative photometric precisions of 0.101%, 0.074%, and 0.076% in $g'_2$, $r'_2$, and $z_{s,2}$ bands, respectively, for GJ436 (magnitudes in $g'$=11.81, $r'$=10.08, and $z'$=8.66) with 30 s exposures. The achieved precisions meet our objective, and the instrument is ready for operation.
We present a set of four three-dimensional, general relativistic, radiation MHD simulations of black hole accretion at super-critical mass accretion rates, $\dot{M} > \dot{M}_{\rm Edd}$. We use these simulations to study how disk properties are modified when we vary the black hole mass, the black hole spin, or the mass accretion rate. In the case of a non-rotating black hole, we find that the total efficiency is of order $3\%\dot M c^2$, approximately a factor of two less than the efficiency of a standard thin accretion disk. The radiation flux in the funnel along the axis is highly super-Eddington, but only a small fraction of the energy released by accretion escapes in this region. The bulk of the $3\%\dot M c^2$ of energy emerges farther out in the disk, either in the form of photospheric emission or as a wind. In the case of a black hole with a spin parameter of 0.7, we find a larger efficiency of about $8\%\dot M c^2$. By comparing the relative importance of advective and diffusive radiation transport, we show that photon trapping is effective near the equatorial plane. However, near the disk surface, vertical transport of radiation by diffusion dominates. We compare the properties of our fiducial three-dimensional run with those of an equivalent two-dimensional axisymmetric model with a mean-field dynamo. The latter simulation runs nearly 100 times faster than the three-dimensional simulation, and gives very similar results for time-averaged properties of the accretion flow.
Observations show the presence, in the halo of our Galaxy, of stars moving at velocities so high to require an acceleration mechanism involving the presence of a massive central black hole. Thus, in the frame of a galaxy hosting a supermassive black hole ($10^8$ $M_{\odot}$) we investigated a mechanism for the production of high velocity stars, which was suggested by the results of N-body simulations of the close interaction between a massive, orbitally decayed, globular cluster and the super massive black hole. The high velocity acquired by some stars of the cluster comes from the transfer of gravitational binding energy into kinetic energy of the escaping star originally orbiting around the cluster. After the close interaction with the massive black hole, stars could reach a velocity sufficient to travel in the halo and even overcome the galactic gravitational well, while some of them are just stripped from the globular cluster and start orbiting on precessing loops around the galactic centre.
The so-called "early activity" of comet 67P/Churyumov-Gerasimenko has been
observed to originate mostly in parts of the concave region or "neck" between
its two lobes. Since activity is driven by the sublimation of volatiles, this
is a puzzling result because this area is less exposed to the Sun and is
therefore expected to be cooler on average (Sierks et al. 2015).
We used a thermophysical model that takes into account thermal inertia,
global self-heating, and shadowing, to compute surface temperatures of the
comet. We found that, for every rotation in the August--December 2014 period,
some parts of the neck region undergo the fastest temperature variations of the
comet's surface precisely because they are shadowed by their surrounding
terrains. Our work suggests that these fast temperature changes are correlated
to the early activity of the comet, and we put forward the hypothesis that
erosion related to thermal cracking is operating at a high rate on the neck
region due to these rapid temperature variations. This may explain why the neck
contains some ice --as opposed to most other parts of the surface (Capaccioni
et al. 2015)-- and why it is the main source of the comet's early activity
(Sierks et al. 2015).
In a broader context, these results indicate that thermal cracking can
operate faster on atmosphereless bodies with significant concavities than
implied by currently available estimates (Delbo' et al. 2014).
We present an application of the stereoscopic self-similar-expansion model (SSSEM) to Solar Terrestrial Relations Observatory (STEREO)/Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI) observations of the 03 April 2010 CME and its associated shock. The aim is to verify whether CME-driven shock parameters can be inferred from the analysis of j-maps. For this purpose we use the SSSEM to derive the CME and the shock kinematics. Arrival times and speeds, inferred assuming either propagation at constant speed or with uniform deceleration, show good agreement with Advanced Composition Explorer (ACE) measurements. The shock standoff distance $[\Delta]$, the density compression $[\frac{\rho_d}{\rho_u}]$ and the Mach number $[M]$ are calculated combining the results obtained for the CME and shock kinematics with models for the shock location. Their values are extrapolated to $\textrm{L}_1$ and compared to in-situ data. The in-situ standoff distance is obtained from ACE solar-wind measurements, and the Mach number and compression ratio are provided by the Harvard-Smithsonian Center for Astrophysics interplanetary shock database. They are $\frac{\rho_d}{\rho_u} =2.84$ and $M = 2.2$. The best fit to observations is obtained when the SSSEM half width $\lambda = 40 \deg$ and the CME and shock propagate with uniform deceleration. In this case we find $\Delta = 23 \textrm{R}_{\odot}$, $\frac{\rho_d}{\rho_u} =2.61$, and $M = 2.93$. The study shows that CME-driven shock parameters can be estimated from the analysis of time-elongation plots and can be used to predict their in-situ values.
A discharge of nitrogen gas, as created in a microwave-induced plasma, exhibits a very dense molecular emission line spectrum. Emission spectra of this kind could serve as wavelength calibrators for high-resolution astrophysical spectrographs in the near-infrared, where only very few calibration sources are currently available. The compilation of a spectral line list and the characterization of line intensities and line density belong to the initial steps when investigating the feasibility of potential wavelength calibration sources. Although the molecular nitrogen spectrum was extensively studied in the past, to our knowledge, no line list exists that covers a continuous range of several thousand wavenumbers in the near-infrared. We recorded three high-resolution ($\Delta \tilde{\nu} = 0.018$cm$^{-1}$) spectra of a nitrogen gas discharge operated at different microwave powers. The nitrogen gas is kept inside a sealed glass cell at a pressure of 2mbar. The emission lines in the spectra were fitted by a superposition of Gaussian profiles to determine their position, relative intensity, and width. The line parameters were corrected for an absolute wavelength scale, instrumental line broadening, and intensity modulation. Molecular and atomic transitions of nitrogen were identified with available line positions from the literature. We report line lists with more than 40000 emission lines in the spectral range $4500-11000$cm$^{-1}$ ($0.9-2.2$$\mu$m). The spectra exhibit emission lines over the complete spectral range under investigation with about $350-1300$ lines per 100cm$^{-1}$. Depending on the microwave power, a fraction of $35\% - 55\%$ of all lines are blended. The total dynamic range of the detected lines covers about four orders of magnitude.
The distribution of subsurface water ice on Mars is a key constraint on past climate, while the volumetric concentration of buried ice (pore-filling versus excess) provides information about the process that led to its deposition. We investigate the subsurface of Arcadia Planitia by measuring the depth of terraces in simple impact craters and mapping a widespread subsurface reflection in radar sounding data. Assuming that the contrast in material strengths responsible for the terracing is the same dielectric interface that causes the radar reflection, we can combine these data to estimate the dielectric constant of the overlying material. We compare these results to a three-component dielectric mixing model to constrain composition. Our results indicate a widespread, decameters-thick layer that is excess water ice ~10^4 km^3 in volume. The accumulation and long-term preservation of this ice is a challenge for current Martian climate models.
The intermittent dissipation of interstellar turbulence is an important energy source in the diffuse ISM. Though on average smaller than the heating rates due to cosmic rays and the photoelectric effect on dust grains, the turbulent cascade can channel large amounts of energy into a relatively small fraction of the gas that consequently undergoes significant heating and chemical enrichment. In particular, this mechanism has been proposed as a solution to the long-standing problem of the high abundance of CH+ along diffuse molecular sight lines, which steady-state, low temperature models under-produce by over an order of magnitude. While much work has been done on the structure and chemistry of these small-scale dissipation zones, comparatively little attention has been paid to relating these zones to the properties of the large-scale turbulence. In this paper, we attempt to bridge this gap by estimating the temperature and CH+ column density along diffuse molecular sight-lines by post-processing 3-dimensional MHD turbulence simulations. Assuming reasonable values for the cloud density (30 / cm^3), size (20 pc), and velocity dispersion (2.3 km / s), we find that our computed abundances compare well with CH+ column density observations, as well as with observations of emission lines from rotationally excited H2 molecules.
The solar chromosphere contains thin, highly dynamic strands of plasma known as spicules. Recently, it has been suggested that the smallest and fastest (Type II) spicules are identical to intermittent jets observed by the Interface Region Imaging Spectrograph. These jets appear to expand out along open magnetic field lines rooted in unipolar network regions of coronal holes. In this paper we revisit a thirty-year-old idea that spicules may be caused by upward forces associated with Alfven waves. These forces involve the conversion of transverse Alfven waves into compressive acoustic-like waves that steepen into shocks. The repeated buffeting due to upward shock propagation causes nonthermal expansion of the chromosphere and a transient levitation of the transition region. Some older models of wave-driven spicules assumed sinusoidal wave inputs, but the solar atmosphere is highly turbulent and stochastic. Thus, we model this process using the output of a time-dependent simulation of reduced magnetohydrodynamic turbulence. The resulting mode-converted compressive waves are strongly variable in time, with a higher transition region occurring when the amplitudes are large and a lower transition region when the amplitudes are small. In this picture, the transition region bobs up and down by several Mm on timescales less than a minute. These motions produce narrow, intermittent extensions of the chromosphere that have similar properties as the observed jets and Type II spicules.
A new class of ultra-long duration (>10,000 s) gamma-ray bursts has recently been suggested. They may originate in the explosion of stars with much larger radii than normal long gamma-ray bursts or in the tidal disruptions of a star. No clear supernova had yet been associated with an ultra-long gamma-ray burst. Here we report that a supernova (2011kl) was associated with the ultra-long duration burst 111209A, at z=0.677. This supernova is more than 3 times more luminous than type Ic supernovae associated with long gamma-ray bursts, and its spectrum is distinctly different. The continuum slope resembles those of super-luminous supernovae, but extends farther down into the rest-frame ultra-violet implying a low metal content. The light curve evolves much more rapidly than super-luminous supernovae. The combination of high luminosity and low metal-line opacity cannot be reconciled with typical type Ic supernovae, but can be reproduced by a model where extra energy is injected by a strongly magnetized neutron star (a magnetar), which has also been proposed as the explanation for super-luminous supernovae.
We explore the potential of using intensity mapping surveys (MeerKAT, SKA) and optical galaxy surveys (DES, LSST) to detect HI clustering and weak gravitational lensing of 21cm emission in auto- and cross-correlation. Our forecasts show that high precision measurements of the clustering and lensing signals can be made in the near future using the intensity mapping technique. Such studies can be used to test the intensity mapping method, and constrain parameters such as the HI density $\Omega_{\rm HI}$, the HI bias $b_{\rm HI}$ and the galaxy-HI correlation coefficient $r_{\rm HI-g}$.
Gamma-ray bursts (GRBs), the most luminous explosions in the universe, invoke relativistic jets beaming towards Earth with the highest velocities for bulk motion in the universe. Some of them are followed by softer, less energetic, X-ray flares, which also move with relativistic velocities towards Earth. Observations and theoretical modeling suggest that X-ray flares share a similar physical mechanism as GRB prompt emission itself. Here we show a clear observational evidence that the X-ray flare emission region is undergoing rapid acceleration as the photons are emitted. The observed X-ray flare light curves and photon index evolution can be interpreted within a simple toy model invoking synchrotron radiation in an accelerating emission region far from the GRB central engine. Such an acceleration process demands an additional energy dissipation source other than kinetic energy, which points towards a significant Poynting-flux in the emission region.
Acceleration of cosmic-ray electrons (CRe) in the intra-cluster-medium (ICM) is probed by radio observations that detect diffuse, Mpc-scale, synchrotron sources in a fraction of galaxy clusters. Giant radio halos are the most spectacular manifestations of non-thermal activity in the ICM and are currently explained assuming that turbulence driven during massive cluster-cluster mergers reaccelerates CRe at several GeV. This scenario implies a hierarchy of complex mechanisms in the ICM that drain energy from large-scales into electromagnetic fluctuations in the plasma and collisionless mechanisms of particle acceleration at much smaller scales. In this paper we focus on the physics of acceleration by compressible turbulence. The spectrum and damping mechanisms of the electromagnetic fluctuations, and the mean-free-path (mfp) of CRe are the most relevant ingredients that determine the efficiency of acceleration. These ingredients in the ICM are however poorly known and we show that calculations of turbulent acceleration are also sensitive to these uncertainties. On the other hand this fact implies that the non-thermal properties of galaxy clusters probe the complex microphysics and the weakly collisional nature of the ICM.
We generalize the model of solid inflation to an anisotropic cosmic solid. Barring fine tunings, the observed isotropy of the cosmological background and of the scalar two-point function isolate the icosahedral group as the only possible symmetry group of such a solid. In such a case, higher-point correlation functions---starting with the three-point one---are naturally maximally anisotropic, which makes the standard detection strategies highly inefficient and calls for a dedicated analysis of CMB data. The tensor two-point function can also be highly anisotropic, but only in the presence of sizable higher-derivative couplings.
We calculate collision strengths and their thermally-averaged Maxwellian values for electron excitation and de-excitation between the fifteen lowest levels of singly-ionised cobalt, Co+, which give rise to emission lines in the near- and mid-infrared. Transition probabilities are also calculated and relative line intensities predicted for conditions typical of supernova ejecta. The diagnostic potential of the 10.52, 15.46 and 14.74 micro-metre transition lines is briefly discussed.
We show that a self-interacting neutrino gas spontaneously acquires a non-stationary pulsating component in its flavor content, with a frequency that can exactly cancel the velocity-dependent refractive effects of dense matter. Spatial inhomogeneities can then lead to small-scale instabilities at large neutrino densities, where the system would have been stable in the homogeneous case. This would allow for large flavor conversions very close to a supernova core, with important consequences for the explosion dynamics and nucleosynthesis, as well as for neutrino observations.
We discuss properties of photons in extremely strong magnetic fields induced by the relativistic heavy-ion collisions. We investigate the vacuum birefringence, the real-photon decay, and the photon splitting which are all forbidden in the ordinary vacuum, but become possible in strong magnetic fields. These effects potentially give rise to anisotropies in photon and dilepton spectra.
We construct a theory in which the gravitational interaction is described only by torsion, but that generalizes the Teleparallel Theory still keeping the invariance of local Lorentz transformations. We show that our theory falls, to a certain limit of a real parameter, in the $f(R)$ Gravity or, to another limit of the same real parameter, in a modified $f(T)$ Gravity, interpolating between these two theories and still can fall on several other theories. We explicitly show the equivalence with $f(R)$ Gravity for cases of Friedmann-Lemaitre-Robertson-Walker flat metric for diagonal tetrads, and a metric with spherical symmetry for diagonal and non-diagonal tetrads.
We present methodology and mission results from orbit determination of Planet Labs nanosatellites and differential-drag control of their relative motion. Orbit determination (OD) is required on Planet Labs satellites to accurately predict the positioning of satellites during downlink passes and we present a scalable OD solution for large fleets of small satellites utilizing two-way ranging. In the second part of this paper, we present mission results from relative motion differential-drag control of a constellation of satellites deployed in the same orbit.
One proposal for dS/CFT is that the Hartle-Hawking (HH) wave function in the large volume limit is equal to the partition function of a Euclidean CFT deformed by various operators. All saddle points defining the semiclassical HH wave function in cosmology have a representation in which their interior geometry is part of a Euclidean AdS domain wall with complex matter fields. We compute the wave functions of scalar and tensor perturbations around homogeneous isotropic complex saddle points, turning on single scalar field matter only. We compare their predictions for the spectra of CMB perturbations with those of a different dS/CFT proposal based on the analytic continuation of inflationary universes to real asymptotically AdS domain walls. We find the predictions of both bulk calculations agree to first order in the slow roll parameters but they differ at higher order.
We determine the Hubble expansion and the general cosmic perturbations equations for a general system consisting of self-conserved matter and self-conserved dark energy (DE). While at the background level the two components are non-interacting, they do interact at the perturbations level. We show that the coupled system of matter and DE perturbations can be transformed into a single, third order, matter perturbation equation, which reduces to the (derivative of the) standard one in the case that the DE is just a cosmological constant. As a nontrivial application we analyze a class of dynamical models whose DE density $\rho_D$ consists of a constant term, $C_0$, and a series of powers of the Hubble rate. These models were previously analyzed from the point of view of dynamical vacuum models, but here we treat them as self-conserved DE models with a dynamical equation of state. We fit them to the wealth of expansion history and linear structure formation data and compare the obtained fit quality with that of the concordance $\Lambda$CDM model. Interestingly, we find that they can be phenomenologically advantageous, except for the generic models with $C_0=0$ and especially the pure linear model $\rho_D\sim H$ (advocated in several places in the literature), which appears strongly disfavored. The remaining models are promising dynamical DE candidates whose phenomenological performance can be highly competitive with the rigid $\Lambda$-term inherent to the $\Lambda$CDM.
The ekpyrotic phase (a slow contraction cosmic phase before the current expansion phase) manages to solve the main problems of the standard cosmology by means of a scalar field interpreted as an isotropic cosmic fluid in the Friedmann equation. Moreover, this phase generates a nearly scale-invariant spectrum of perturbations in agreement with the latest data. Then, the ekpyrotic mechanism is a serious possibility to the inflationary model. In this work, we point out that it is impossible to generate a black hole with spherical symmetry supported by an isotropic fluid in this scenario. Using the approach of deforming metrics to obtain solutions with an isotropic energy-momentum tensor, we show that the stiff fluid, dominant in the ekpyrotic phase, does not support these black holes.
Links to: arXiv, form interface, find, astro-ph, recent, 1509, contact, help (Access key information)