We study how optical galaxy morphology depends on mass and star formation rate (SFR) in the Illustris Simulation. To do so, we measure automated diagnostics of galaxy structure in 10808 simulated galaxies at z=0 with stellar masses 10^9.7 < M_*/M_sun < 10^12.3. We add observational realism to idealized synthetic images and measure non-parametric statistics in rest-frame u, g, i, and H band images from four directions. We find that the Illustris simulation creates a morphologically diverse population of galaxies, occupying roughly the observed bulge strength locus, and reproducing median morphology trends versus stellar mass, SFR, and compactness. Optical morphology correlates realistically with rotational structure, following qualitative classification schemes put forth by kinematic surveys. Relative type fractions as a function of environment agree roughly with data. These results imply that connections among mass, star formation, and galaxy structure arise naturally from models matching global star formation and halo occupation functions when simulated with accurate numerical methods. This raises a question of how to construct the best experiments on large galaxy surveys to better distinguish between models. We predict that at fixed halo mass near 10^12 M_sun, galaxies with relatively more disc-like morphologies have higher stellar mass than those with bulge-like morphologies, a possible consequence of the Illustris feedback model acting on massive galaxies. While Illustris galaxies at M_* ~ 10^11 M_sun have a reasonable size distribution, those at M_* ~ 10^10 M_sun have half-light radii larger than observed by roughly a factor of 2. Furthermore, at M_* ~ 10^10.5 through 10^11 M_sun, a relevant fraction of Illustris galaxies have distinct "ring-like" features, such that the bright pixels have an unusually wide spatial extent (M_20 > -1).
We present the discovery of two z > 6 quasars, selected as i band dropouts in the VST ATLAS survey. Our first quasar has redshift, z = 6.31 \pm 0.03, z band magnitude, zAB = 19.63 \pm 0.08 and rest frame 1450A absolute magnitude, M1450 = -27.8 \pm 0.2, making it one of the two most luminous quasars known at z > 6. The second quasar has z = 6.02 \pm 0.03, zAB = 19.54 \pm 0.08 and M1450 = -27.0 \pm 0.1. We also recover a z = 5.86 quasar discovered by Venemans et al. (2015, in prep.). To select our quasars we use a new 3D colour space, combining the ATLAS optical colours with mid-infra-red data from the WISE Space Telescope. We use iAB - zAB colour to exclude main sequence stars, galaxies and lower redshift quasars, W1 - W2 to exclude L dwarfs and zAB - W2 to exclude T dwarfs. A restrictive set of colour cuts returns only our three high redshift quasars and no contaminants, albeit with a sample completeness of \sim50%. We discuss how less restrictive cuts in our 3D colour space can be used to reject the majority of contaminants from samples of bright 5.7 < z < 6.4 quasars, replacing follow-up near-infra-red photometry, whilst retaining high completeness.
We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.
"Green Bean" Galaxies (GBs) are the most [O III]-luminous type-2 active galactic nuclei (AGN) at z~0.3. However, their infrared luminosities reveal AGN in very low activity states, indicating that their gas reservoirs must be ionized by photons from a recent high activity episode - we are observing quasar ionization echoes. We use integral field spectroscopy from the Gemini Multi-Object Spectrograph to analyse the 3D kinematics, ionization state, temperature and density of ionized gas in the GB SDSS J224024.1-092748. We model the emission line spectrum of each spaxel as a superposition of up to three Gaussian components and analyse the physical properties of each component individually. Two narrow components, tracing the velocity fields of the disc and an ionized gas cloud, are superimposed over the majority of the galaxy. Fast shocks produce hot ($T_e$ $\geq$ 20,000 K), dense ($n_e$ $\geq$ 100 cm$^{-3}$), turbulent ($\sigma$ $\geq$ 600 km s$^{-1}$), [O III]-bright regions with enhanced [N II]/H$\alpha$ and [S II]/H$\alpha$ ratios. The most prominent such spot is consistent with a radio jet shock-heating the interstellar medium. However, the AGN is still responsible for $\geq$ 82 per cent of the galaxy's total [O III] luminosity, strengthening the case for previous quasar activity. The ionized gas cloud has a strong kinematic link to the central AGN and is co-rotating with the main body of the galaxy, suggesting that it may be the remnant of a quasar-driven outflow. Our analysis of J224024.1-092748 indicates that GBs provide a unique fossil record of the transformation from the most luminous quasars to weak AGN.
The 5th edition of the Roma-BZCAT Multifrequency Catalogue of Blazars is available in a printed version and online at the ASDC website (this http URL); it is also in the NED database. It presents several relevant changes with respect to the past editions which are briefly described in this paper.
In this study we investigate the relationship between the star formation rate, SFR, and AGN luminosity, L(AGN), for ~2000 X-ray detected AGN. The AGN span over three orders of magnitude in X-ray luminosity (10^(42) < L(2-8keV) < 10^(45.5) erg/s) and are in the redshift range z = 0.2 - 2.5. Using infrared (IR) photometry (8 - 500um), including deblended Spitzer and Herschel images and taking into account photometric upper limits, we decompose the IR spectral energy distributions into AGN and star formation components. Using the IR luminosities due to star formation, we investigate the average SFRs as a function of redshift and AGN luminosity. In agreement with previous studies, we find a strong evolution of the average SFR with redshift, tracking the observed evolution of the overall star forming galaxy population. However, we find that the relationship between the average SFR and AGN luminosity is flat at all redshifts and across all the AGN luminosities investigated. By comparing to empirical models, we argue that the observed flat relationship is due to short timescale variations in AGN luminosity, driven by changes in the mass accretion rate, which wash out any underlying correlations between SFR and L(AGN). Furthermore, we show that the exact form of the predicted relationship between SFR and AGN luminosity (and it's normalisation) is highly sensitive to the assumed intrinsic Eddington ratio distribution.
We seek to improve the accuracy of joint galaxy photometric redshift estimation and spectral energy distribution (SED) fitting. By simulating different sources of uncorrected systematic errors, we demonstrate that if the uncertainties on the photometric redshifts are estimated correctly, so are those on the other SED fitting parameters, such as stellar mass, stellar age, and dust reddening. Furthermore, we find that if the redshift uncertainties are over(under)-estimated, the uncertainties in SED parameters tend to be over(under)-estimated by similar amounts. These results hold even in the presence of severe systematics and provide, for the first time, a mechanism to validate the uncertainties on these parameters via comparison with spectroscopic redshifts. We propose a new technique (annealing) to re-calibrate the joint uncertainties in the photo-z and SED fitting parameters without compromising the performance of the SED fitting + photo-z estimation. This procedure provides a consistent estimation of the multidimensional probability distribution function in SED fitting + z parameter space, including all correlations. While the performance of joint SED fitting and photo-z estimation might be hindered by template incompleteness, we demonstrate that the latter is "flagged" by a large fraction of outliers in redshift, and that significant improvements can be achieved by using flexible stellar populations synthesis models and more realistic star formation histories. In all cases, we find that the median stellar age is better recovered than the time elapsed from the onset of star formation [abridged].
AGB stars are responsible for producing a variety of elements, including carbon, nitrogen, and the heavy elements produced in the slow neutron-capture process ($s$-elements). There are many uncertainties involved in modelling the evolution and nucleosynthesis of AGB stars, and this is especially the case at low metallicity, where most of the stars with high enough masses to enter the AGB have evolved to become white dwarfs and can no longer be observed. The stellar population in the Galactic halo is of low mass ($\lesssim 0.85M_{\odot}$) and only a few observed stars have evolved beyond the first giant branch. However, we have evidence that low-metallicity AGB stars in binary systems have interacted with their low-mass secondary companions in the past. The aim of this work is to investigate AGB nucleosynthesis at low metallicity by studying the surface abundances of chemically peculiar very metal-poor stars of the halo observed in binary systems. To this end we select a sample of 15 carbon- and $s$-element-enhanced metal-poor (CEMP-$s$) halo stars that are found in binary systems with measured orbital periods. With our model of binary evolution and AGB nucleosynthesis, we determine the binary configuration that best reproduces, at the same time, the observed orbital period and surface abundances of each star of the sample. The observed periods provide tight constraints on our model of wind mass transfer in binary stars, while the comparison with the observed abundances tests our model of AGB nucleosynthesis.
We present broadband (3 -- 78 keV) NuSTAR X-ray imaging and spectroscopy of the Crab nebula and pulsar. We show that while the phase-averaged and spatially integrated nebula + pulsar spectrum is a power-law in this energy band, spatially resolved spectroscopy of the nebula finds a break at $\sim$9 keV in the spectral photon index of the torus structure with a steepening characterized by $\Delta\Gamma\sim0.25$. We also confirm a previously reported steepening in the pulsed spectrum, and quantify it with a broken power-law with break energy at $\sim$12 keV and $\Delta\Gamma\sim0.27$. We present spectral maps of the inner 100\as\ of the remnant and measure the size of the nebula as a function of energy in seven bands. These results find that the rate of shrinkage with energy of the torus size can be fitted by a power-law with an index of $\gamma = 0.094\pm 0.018$, consistent with the predictions of Kennel and Coroniti (1984). The change in size is more rapid in the NW direction, coinciding with the counter-jet where we find the index to be a factor of two larger. NuSTAR observed the Crab during the latter part of a $\gamma$-ray flare, but found no increase in flux in the 3 - 78 keV energy band.
The broad line region (BLR) of luminous active galactic nuclei (AGN) is a prominent observational signature of the accretion flow around supermassive black holes, which can be used to measure their masses (M_BH) over cosmic history. Due to the <100 {\mu}as angular size of the BLR, current direct constraints on BLR kinematics are limited to those provided by reverberation mapping studies, which are most efficiently carried out on low-luminosity L and low-redshift z AGN. We analyze the possibility to measure the BLR size and study its kinematic structure using spectroastrometry, whereby one measures the spatial position centroid of emission line photons as a function of velocity. We calculate the expected spectroastrometric signal of a rotation-dominated BLR for various assumptions about the ratio of random to rotational motions, and the radial distribution of the BLR gas. We show that for hyper-luminous quasars at z < 2.5, the size of the low-ionization BLR can already be constrained with existing telescopes and adaptive optics systems, thus providing a novel method to spatially resolve the kinematics of the accretion flow at 10^3 -- 10^4 gravitational radii, and measure M_BH at the high-L end of the AGN family. With a 30m-class telescope, BLR spectroastrometry should be routinely detectable for much fainter quasars out to z ~ 6, and for various emission lines. This will enable kinematic M_BH measurements as a function of luminosity and redshift, providing a compelling science case for next generation telescopes.
The mobility of atoms, molecules and radicals in icy grain mantles regulate ice restructuring, desorption, and chemistry in astrophysical environments. Interstellar ices are dominated by H2O, and diffusion on external and internal (pore) surfaces of H2O-rich ices is therefore a key process to constrain. This study aims to quantify the diffusion kinetics and barrier of the abundant ice constituent CO into H2O dominated ices at low temperatures (15-23 K), by measuring the mixing rate of initially layered H2O(:CO2)/CO ices. The mixed fraction of CO as a function of time is determined by monitoring the shape of the infrared CO stretching band. Mixing is observed at all investigated temperatures on minute time scales, and can be ascribed to CO diffusion in H2O ice pores. The diffusion coefficient and final mixed fraction depend on ice temperature, porosity, thickness and composition. The experiments are analyzed by applying Fick's diffusion equation under the assumption that mixing is due to CO diffusion into an immobile H2O ice. The extracted energy barrier for CO diffusion into amorphous H2O ice is ~160 K. This is effectively a surface diffusion barrier. The derived barrier is low compared to current surface diffusion barriers in use in astrochemical models. Its adoption may significantly change the expected timescales for different ice processes in interstellar environments.
The core rotation rates of massive stars have a substantial impact on the nature of core collapse supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited via envelope convection during a red supergiant phase or during vigorous late time burning phases, can have a significant impact on the rotation rate of the pre-SN core. In typical ($10 \, M_\odot \lesssim M \lesssim 20 \, M_\odot$) supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods $P_{\rm min,Fe} \gtrsim 50 \, {\rm s}$. Angular momentum (AM) conservation during the supernova would entail minimum NS rotation periods of $P_{\rm min,NS} \gtrsim 3 \, {\rm ms}$. In most cases, the combined effects of magnetic torques and IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods of $P_{\rm max,Fe} \lesssim 10^4 \, {\rm s}$ in typical core collapse supernova progenitors, and a corresponding spin period of $P_{\rm max, NS} \lesssim 400 \, {\rm ms}$ for newborn neutron stars. This is comparable to the typical birth spin periods of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying envelope.
We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum, and soft X-ray features observed in this galaxy, pointing to a tight disk-halo connection regulated by star formation. On the other hand, the morphology of the colder dust component detected on larger scale in the SPIRE 250, 350, and 500 um data matches the extraplanar H~I streams previously reported in NGC 4631 and suggests a tidal origin. The PACS 70/160 ratios are elevated in the central ~3.0 kpc region above the nucleus of this galaxy (the "superbubble"). A pixel-by-pixel analysis shows that dust in this region has a higher temperature and/or an emissivity with a steeper spectral index (beta > 2) than the dust in the disk, possibly the result of the harsher environment in the superbubble. Star formation in the disk seems energetically insufficient to lift the material out of the disk, unless it was more active in the past or the dust-to-gas ratio in the superbubble region is higher than the Galactic value. Some of the dust in the halo may also have been tidally stripped from nearby companions or lifted from the disk by galaxy interactions.
Several mechanisms have been proposed to explain the transient seis- mic emission, i.e., sunquakes, from some solar flares. Some theories associate high-energy electrons and/or white-light emission with sunquakes. High-energy charged particles and their subsequent heating of the photosphere and/or chro- mosphere could induce acoustic waves in the solar interior. We carried out a correlative study of solar flares with emission in hard-X rays (HXRs), enhanced continuum emission at 6173{\AA}, and transient seismic emission. We selected those flares observed by RHESSI (Reuven Ramaty High Energy Solar Spectroscopic Imager) with a considerable flux above 50 keV between January 1, 2010 and June 26, 2014. We then used data from the Helioseismic and Magnetic Imager onboard the Solar Dynamic Observatory (SDO/HMI) to search for excess visible continuum emission and new sunquakes not previously reported. We found a total of 18 sunquakes out of 75 investigated. All of the sunquakes were associated with a enhancement of the visible continuum during the flare time. Finally, we calculated a coefficient of correlation for a set of dichotomic variables related to these observations. We found a strong correlation between two of the standard helioseismic detection techniques, and between sunquakes and visible continuum enhancements. We discuss the phenomenological connectivity between these physical quantities and the observational difficulties of detecting seismic signals and excess continuum radiation.
Helical magnetic flux rope (MFR) is a fundamental structure of corona mass ejections (CMEs) and has been discovered recently to exist as a sigmoidal channel structure prior to its eruption in the extreme ultraviolet (EUV) high temperature passbands of the Atmospheric Imaging Assembly (AIA). However, when and where the MFR is built up are still elusive. In this paper, we investigate two MFRs (MFR1 and MFR2) in detail, whose eruptions produced two energetic solar flares and CMEs on 2014 April 18 and 2014 September 10, respectively. The AIA EUV images reveal that for a long time prior to their eruption, both MFR1 and MFR2 are under formation, which is probably through magnetic reconnection between two groups of sheared arcades driven by the shearing and converging flows in the photosphere near the polarity inversion line. At the footpoints of the MFR1, the \textit{Interface Region Imaging Spectrograph} Si IV, C II, and Mg II lines exhibit weak to moderate redshifts and a non-thermal broadening in the pre-flare phase. However, a relatively large blueshift and an extremely strong non-thermal broadening are found at the formation site of the MFR2. These spectral features consolidate the proposition that the reconnection plays an important role in the formation of MFRs. For the MFR1, the reconnection outflow may propagate along its legs, penetrating into the transition region and the chromosphere at the footpoints. For the MFR2, the reconnection probably takes place in the lower atmosphere and results in the strong blueshift and non-thermal broadening for the Mg II, C II, and Si IV lines.
We present an advance in the use of Cassini observations of stellar occultations by the rings of Saturn for stellar studies. Stewart et al. (2013) demonstrated the potential use of such observations for measuring stellar angular diameters. Here, we use these same observations, and tomographic imaging reconstruction techniques, to produce two dimensional images of complex stellar systems. We detail the determination of the basic observational reference frame. A technique for recovering model-independent brightness profiles for data from each occulting edge is discussed, along with the tomographic combination of these profiles to build an image of the source star. Finally we demonstrate the technique with recovered images of the {\alpha} Centauri binary system and the circumstellar environment of the evolved late-type giant star, Mira.
Dissipative dark matter, where dark matter particles interact with a massless (or very light) boson, is studied. Such dark matter can arise in simple hidden sector gauge models, including those featuring an unbroken $U(1)'$ gauge symmetry, leading to a dark photon. Previous work has shown that such models can not only explain the LSS and CMB, but potentially also dark matter phenomena on small scales, such as the inferred cored structure of dark matter halos. In this picture, dark matter halos of disk galaxies not only cool via dissipative interactions but are also heated via ordinary supernovae (facilitated by an assumed photon - dark photon kinetic mixing interaction). This interaction between the dark matter halo and ordinary baryons, a very special feature of these types of models, plays a critical role in governing the physical properties of the dark matter halo. Here, we further study the implications of this type of dissipative dark matter for disk galaxies. Building on earlier work, we develop a simple formalism which aims to describe the effects of dissipative dark matter in a fairly model independent way. This formalism is then applied to generic disk galaxies. We also consider specific examples, including NGC 1560 and a sample of dwarf galaxies from the LITTLE THINGS survey. We find that dissipative dark matter, as developed here, does a fairly good job accounting for the rotation curves of the galaxies considered. Not only does dissipative dark matter explain the linear rise of the rotational velocity of dwarf galaxies at small radii, but it can also explain the observed wiggles in rotation curves which are known to be correlated with corresponding features in the disk gas distribution.
The discovery of rings around extrasolar planets ("exorings") is one of the next breakthroughs in exoplanetary research. Previous studies have explored the feasibility of detecting exorings with present and future photometric sensitivities by seeking anomalous deviations in the residuals of a standard transit light curve fit, at the level of ~100 ppm for Kronian rings. In this work, we explore two much larger observational consequences of exorings: (1) the significant increase in transit depth that may lead to misclassification of ringed planetary candidates as false-positives and/or the underestimation of planetary density; and (2) the so-called "photo-ring" effect, a new asterodensity profiling effect, revealed by a comparison of the light curve derived stellar density to that measured with independent methods (e.g. asteroseismology). Whilst these methods do not provide an unambiguous discovery of exorings, we show that the large amplitude of these effects combined with their relatively simple analytic description, makes them highly suited to large scale surveys to identify candidate ringed planets worthy of more detailed investigation. Moreover, these methods lend themselves to ensemble analyses seeking to uncover evidence for a population of ringed planets. We describe the method in detail, develop the basic underlying formalism and test it in the parameter space of rings and transit configuration. We discuss the prospects of using the method for the first systematic search of exoplanetary rings in the Kepler database and provide basic computational code for implementing it.
The polarization signatures of the blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, sometimes large (> 180^o) polarization angle swings are observed. We suggest that such p henomena can be interpreted as arising from light-travel-time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change of its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic-field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.
Tidal disruption events (TDEs) occur when a star, passing too close to a massive black hole, is ripped apart by tidal forces. A less dramatic event occurs if the star orbits just outside the tidal radius, resulting in a mild stripping of mass. Thus, if a star orbits a central black hole on one of these bound eccentric orbits, weaker outbursts will occur recurring every orbital period. Thanks to five Swift observations, we observed a recent flare from the close by (92 Mpc) galaxy IC 3599, where a possible TDE was already observed in December 1990 during the Rosat All-Sky Survey. By light curve modeling and spectral fitting, we account for all these events as the non-disruptive tidal stripping of a single star into a 9.5 yr highly eccentric bound orbit. This is the first example of periodic partial tidal disruptions, possibly spoon-feeding the central black hole.
The model of current-carrying string loop oscillations is tested to explain the special set of frequencies related to the high-frequency quasiperiodic oscillations (HF QPOs) observed recently in the low-mass X-ray binary XTE J1701-407 containing a neutron star. The external geometry of the neutron star is approximated by the Kerr geometry, introducing errors not exceeding $10~\%$ for slowly rotating massive neutron stars. The frequencies of the radial and vertical string loop oscillations are then governed by the mass $M$ and dimensionless spin $a$ of the neutron star, and by the dimensionless parameter $\omega$ describing combined effects of the string loop tension and its angular momentum. It is explicitly demonstrated that the string-loop oscillation model can explain the observed kHz frequencies for the neutron star parameters restricted to the intervals ${0.2<a<0.4}$ and ${2.1<M/{\rm M}_{\odot}<2.5}$. However, the stringy parameter $\omega$ cannot be the same for all the three HF QPO observations in the XTE J1701-407 source; the limits on the acceptable values of $\omega$ are given in dependence on the spacetime parameters $M$ and $a$.
The CMB temperature-redshift relation, T_CMB(z)=T_0(1+z), is a key prediction of the standard cosmology, but is violated in many non standard models. Constraining possible deviations to this law is an effective way to test the LambdaCDM paradigm and to search for hints of new physics. We have determined T_CMB(z), with a precision up to 3%, for a subsample (104 clusters) of the Planck SZ cluster catalog, at redshift in the range 0.01-- 0.94, using measurements of the spectrum of the Sunyaev Zel'dovich effect obtained from Planck temperature maps at frequencies from 70 to 353 GHz. The method adopted to provide individual determinations of T_CMB(z) at cluster redshift relies on the use of SZ intensity change, Delta I_SZ(nu), at different frequencies, and on a Monte-Carlo Markov Chain approach. By applying this method to the sample of 104 clusters, we limit possible deviations of the form T_CMB(z)=T_0(1+z)^(1-beta) to be beta= 0.022 +/- 0.018, at 1 sigma uncertainty, consistent with the prediction of the standard model. Combining these measurements with previously published results we get beta=0.016+/-0.012.
The internal shocks scenario in relativistic jets is used to explain the variability of the blazar emission. Recent studies have shown that the magnetic field significantly alters the shell collision dynamics, producing a variety of spectral energy distributions and light-curves patterns. However, the role played by magnetization in such emission processes is still not entirely understood. In this work we numerically solve the magnetohydodynamic evolution of the magnetized shells collision, and determine the influence of the magnetization on the observed radiation. Our procedure consists in systematically varying the shell Lorentz factor, relative velocity, and viewing angle. The calculations needed to produce the whole broadband spectral energy distributions and light-curves are computationally expensive, and are achieved using a high-performance parallel code.
The physical properties of galactic winds are of paramount importance for our understanding of galaxy formation. Fortunately, they can be constrained using background quasars passing near star-forming galaxies (SFGs). From the 14 quasar$-$galaxy pairs in our VLT/SINFONI Mgii Program for Line Emitters (SIMPLE) sample, we reobserved the 10 brightest galaxies in H$_{\alpha}$ with the VLT/SINFONI with 0.7" seeing and the corresponding quasar with the VLT/UVES spectrograph. Applying geometrical arguments to these ten pairs, we find that four are likely probing galactic outflows, three are likely probing extended gaseous disks, and the remaining three are not classifiable because they are viewed face-on. In this paper we present a detailed comparison between the line-of-sight kinematics and the host galaxy emission kinematics for the pairs suitable for wind studies. We find that the kinematic profile shapes (asymmetries) can be well reproduced by a purely geometrical wind model with a constant wind speed, except for one pair (towards J2357$-$2736) that has the smallest impact parameter b = 6 kpc and requires an accelerated wind flow. Globally, the outflow speeds are $\sim$ 100 km/s and the mass ejection rates (or $\dot M _{\rm out}$) in the gas traced by the low-ionization species are similar to the star formation rate (SFR), meaning that the mass loading factor, $\eta$ = $\dot M _{\rm out}$/SFR, is $\sim$1.0. The outflow speeds are also smaller than the local escape velocity, which implies that the outflows do not escape the galaxy halo and are likely to fall back into the interstellar medium.
FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.
The physical origin of the X-ray emission in powerful quasar jets has been a long-standing mystery. Though these jets start out on the sub-pc scale as highly relativistic flows, we do not have any direct measurement of their speeds on the kpc scale, where the vast distances from the core necessitate in situ particle acceleration. If the jets remain highly relativistic on kpc scales, then the X-rays could be due to inverse-Compton upscattering of CMB photons. However, the IC/CMB explanation predicts a high level of gamma-ray emission, which should be detectible by the Fermi/LAT. We have searched for and ruled out this emission at a high level of significance for the well-known sources 3C 273 and PKS 0637-752, suggesting the X-rays are synchrotron, though of unknown origin. These recent results with Fermi also suggest that the kpc-scale jets in powerful quasars are significantly slower than have been presumed under the IC/CMB model. I will discuss the surprising implications of these findings for the energetics and radiative output of powerful quasars as well as their impact on their environment.
We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.
A multi-epoch, multi-instrument analysis of the Seyfert 1 galaxy HE 0436-4717 is conducted using optical to X-ray data from XMM-Newton and Swift (including the BAT). Fitting of the UV-to-X-ray spectral energy distribution shows little evidence of extinction and the X-ray spectral analysis does not confirm previous reports of deep absorption edges from OVIII. HE 0436-4717 is a "bare" Seyfert with negligible line-of-sight absorption making it ideal to study the central X-ray emitting region. Three scenarios were considered to describe the X-ray data: partial covering absorption, blurred reflection, and soft Comptonization. All three interpretations describe the 0.5-10.0 keV spectra well. Extrapolating the models to 100 keV results in poorer fits for the the partial covering model. When also considering the rapid variability during one of the XMM-Newton observations, the blurred reflection model appears to describe all the observations in the most self-consistent manner. If adopted, the blurred reflection model requires a very low iron abundance in HE 0436-4717. We consider the possibilities that this is an artifact of the fitting process, but it appears possible that it is intrinsic to the object.
We present the analysis of the integrated spectral energy distribution (SED) from the ultraviolet (UV) to the far-infrared and H$\alpha$ of a sample of 29 local systems and individual galaxies with infrared (IR) luminosities between 10^11 Lsun and 10^11.8 Lsun. We have combined new narrow-band H$\alpha$+[NII] and broad-band g, r optical imaging taken with the Nordic Optical Telescope (NOT), with archival GALEX, 2MASS, Spitzer, and Herschel data. The SEDs (photometry and integrated H$\alpha$ flux) have been fitted with a modified version of the MAGPHYS code using stellar population synthesis models for the UV-near-IR range and thermal emission models for the IR emission taking into account the energy balance between the absorbed and re-emitted radiation. From the SED fits we derive the star-formation histories (SFH) of these galaxies. For nearly half of them the star-formation rate appears to be approximately constant during the last few Gyrs. In the other half, the current star-formation rate seems to be enhanced by a factor of 3-20 with respect to that occured ~1 Gyr ago. Objects with constant SFH tend to be more massive than starbursts and they are compatible with the expected properties of a main-sequence (M-S) galaxy. Likewise, the derived SFHs show that all our objects were M-S galaxies ~1 Gyr ago with stellar masses between 10^10.1 and 10^11.5 Msun. We also derived from our fits the average extinction (A_v=0.6-3 mag) and the polycyclic aromatic hydrocarbons (PAH) luminosity to L(IR) ratio (0.03-0.16). We combined the A_v with the total IR and H$\alpha$ luminosities into a diagram which can be used to identify objects with rapidly changing (increasing or decreasing) SFR during the last 100 Myr.
A new room temperature line list for $^{15}$NH$_3$ is presented. This line list comprised of transition frequencies and Einstein coefficients has been generated using the `spectroscopic' potential energy surface NH3-Y2010 and an ab initio dipole moment surface. The $^{15}$NH$_3$ line list is based on the same computational procedure used for the line list for $^{14}$NH$_3$ BYTe reported recently and should be as accurate. Comparisons with experimental frequencies and intensities are presented. The synthetic spectra show excellent agreement with experimental spectra.
TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.
We study the global SF law - the relation between gas and SFRs in a sample of 181 local galaxies with L_IR spanning almost five orders of magnitude, which includes 115 normal galaxies and 66 (U)LIRGs. We derive their atomic, molecular gas and dense molecular gas masses using newly available HI, CO and HCN data from the literature, and SFRs are determined both from total IR and 1.4 GHz radio continuum (RC) luminosities. In order to derive the disk-averaged surface densities of gas and SFRs, we have used high-resolution RC observations to measure the radio sizes for all galaxies. We find that dense molecular gas (as traced by HCN) has the tightest correlation with that of SFRs, and is linear in (N=1.01 +/- 0.02) across the full galaxy sample. The correlation between densities of molecular gas (traced by CO) and SFRs is sensitive to the adopted value of the alpha_CO used to infer molecular gas masses from CO luminosities. For a fixed value of alpha_CO, a slope of 1.14+/-0.02 is found. If instead we adopt values of 4.6 and 0.8 for disk galaxies and (U)LIRGs, respectively, we find the two distinct relations. If applying a continuously varying alpha_CO to our sample, we recover a single relation with slope of 1.60+/-0.03. The SFRs is a steeper function of total gas than that of molecular gas, and is tighter among low-luminosity galaxies. We find no correlation between SFRs and atomic gas.
Differential Optical Transfer Function (dOTF) is an image-based, non-iterative wavefront sensing method that uses two star images with a single small change in the pupil. We describe two possible methods for introducing the required pupil modification to the JWST, one using a small (<lambda/4) displacement of a single segment's actuator and another that uses small misalignments of NIRCam's filter wheel. While both methods should work with NIRCam, the actuator method will allow both MIRI and NIRISS to be used for segment phasing, which is new functionality. Since the actuator method requires only small displacements, it should provide a fast and safe phasing alternative that reduces mission risk and can be performed frequently for alignment monitoring and maintenance. Since a single actuator modification can be seen by all three cameras, it should be possible to calibrate the non-common-path aberrations between them. Large segment discontinuities can be measured using dOTFs in two filter bands. Using two images of a star field, aberrations along multiple lines of sight through the telescope can be measured simultaneously. Also, since dOTF gives the pupil field amplitude as well as phase, it could provide a first approximation or constraint to the planned iterative phase retrieval algorithms.
Photoluminescence spectra show that silicon impurity is present in lattice of some nanodiamond grains (ND) of various chondrites as a silicon-vacancy (SiV) defect. The relative intensity of the SiV band in the diamond-rich separates depends on chemical composition of meteorites and on size of ND grains. The strongest signal is found for the size separates enriched in small grains; thus confirming our earlier conclusion that the SiV defects preferentially reside in the smallest (less than 2 nm) grains. The difference in relative intensities of the SiV luminescence in the diamond-rich separates of individual meteorites are due to variable conditions of thermal metamorphism of their parent bodies and/or uneven sampling of nanodiamonds populations. Annealing of separates in air eliminates surface sp2-carbon, consequently, the SiV luminescence is enhanced. Strong and well-defined luminescence and absorption of the SiV defect is a promising feature to locate cold (< 250 {\deg}C) nanodiamonds in space.
For more than three years now we have been conducting a spectroscopic survey of detached eclipsing binaries (DEBs) from the All-Sky Automated Survey (ASAS) database. Thousands of high-resolution spectra of over 300 systems were secured, and used for radial velocity measurements and spectral analysis. In our sample we found a zoo of multiple systems, such as spectroscopic triples and quadruples, visual binaries with eclipsing components, and circumbinary low-mass companions, including sub-stellar-mass candidates
We critically review the role of cosmological moduli in determining the post-inflationary history of the universe. Moduli are ubiquitous in string and M-theory constructions of beyond the Standard Model physics, where they parametrize the geometry of the compactification manifold. For those with masses determined by supersymmetry breaking this leads to their eventual decay slightly before Big Bang Nucleosynthesis (without spoiling its predictions). This results in a matter dominated phase shortly after inflation ends, which can influence baryon and dark matter genesis, as well as observations of the Cosmic Microwave Background and the growth of large-scale structure. Given progress within fundamental theory, and guidance from dark matter and collider experiments, non-thermal histories have emerged as a robust and theoretically well-motivated alternative to a strictly thermal one. We review this approach to the early universe and discuss both the theoretical challenges and the observational implications.
A key component of explaining the array of galaxies observed in the Universe is the feedback of active galactic nuclei, each powered by a massive black hole's accretion disc. For accretion to occur, angular momentum must be lost by that which is accreted. Electromagnetic radiation must offer some respite in this regard, the contribution for which is quantified in this paper using solely general relativity under the thin-disc regime. Herein, I calculate extremised situations where photons are entirely responsible for energy removal in the disc and then extend and relate this to the standard relativistic accretion disc outlined by Novikov & Thorne that includes the effect of viscosity. While there is potential for the contribution of angular-momentum removal from photons to be >~1% out to ~10^4 Schwarzschild radii, especially if the disc is irradiated and is liberated of angular momentum through scattering, it is more likely of order 10^2 Schwarzschild radii if thermal emission from the disc itself is stronger. Near the horizons of fast-spinning black holes, these modes of angular-momentum liberation become dominant.
Simulating a binary black hole coalescence by solving Einstein's equations is computationally expensive, requiring days to months of supercomputing time. In this paper, we construct an accurate and fast-to-evaluate surrogate model for numerical relativity (NR) waveforms from non-spinning binary black hole coalescences with mass ratios from $1$ to $10$ and durations corresponding to about $15$ orbits before merger. Our surrogate, which is built using reduced order modeling techniques, is distinct from traditional modeling efforts. We find that the full multi-mode surrogate model agrees with waveforms generated by NR to within the numerical error of the NR code. In particular, we show that our modeling strategy produces surrogates which can correctly predict NR waveforms that were {\em not} used for the surrogate's training. For all practical purposes, then, the surrogate waveform model is equivalent to the high-accuracy, large-scale simulation waveform but can be evaluated in a millisecond to a second depending on the number of output modes and the sampling rate. Our model includes all spherical-harmonic ${}_{-2}Y_{\ell m}$ waveform modes that can be resolved by the NR code up to $\ell=8$, including modes that are typically difficult to model with other approaches. We assess the model's uncertainty, which could be useful in parameter estimation studies seeking to incorporate model error. We anticipate NR surrogate models to be useful for rapid NR waveform generation in multiple-query applications like parameter estimation, template bank construction, and testing the fidelity of other waveform models.
Inflation is nowadays a well-established paradigm consistent with all the observations. The precise nature of the inflaton is however unknown and its role could be played by any candidate able to imitate a scalar condensate in the slow-roll regime. The discovery of a fundamental scalar in the LHC provides the less speculative candidate. Could the Higgs field itself be responsible for inflation? Do we really need to advocate new physics to explain the properties of the Universe at large scales? Which is the relation between the Standard Model parameters and the inflationary observables? What happens if our vacuum becomes unstable below the scale of inflation? We present an overview of Higgs inflation trying to provide answers to the previous questions with special emphasis on the vacuum stability issue.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
The physical mechanisms of the quasar ultraviolet (UV)-optical variability are not well understood despite the long history of observations. Recently, Dexter & Agol presented a model of quasar UV-optical variability, which assumes large local temperature fluctuations in the quasar accretion discs. This inhomogeneous accretion disc model is claimed to describe not only the single-band variability amplitude, but also microlensing size constraints and the quasar composite spectral shape. In this work, we examine the validity of the inhomogeneous accretion disc model in the light of quasar UV-optical spectral variability by using five-band multi-epoch light curves for nearly 9 000 quasars in the Sloan Digital Sky Survey (SDSS) Stripe 82 region. By comparing the values of the intrinsic scatter $\sigma_{\text{int}}$ of the two-band magnitude-magnitude plots for the SDSS quasar light curves and for the simulated light curves, we show that Dexter & Agol's inhomogeneous accretion disc model cannot explain the tight inter-band correlation often observed in the SDSS quasar light curves. This result leads us to conclude that the local temperature fluctuations in the accretion discs are not the main driver of the several years' UV-optical variability of quasars, and consequently, that the assumption that the quasar accretion discs have large localized temperature fluctuations is not preferred from the viewpoint of the UV-optical spectral variability.
Candidates for the modest galaxies that formed most of the stars in the early universe, at redshifts $z > 7$, have been found in large numbers with extremely deep restframe-UV imaging. But it has proved difficult for existing spectrographs to characterise them in the UV. The detailed properties of these galaxies could be measured from dust and cool gas emission at far-infrared wavelengths if the galaxies have become sufficiently enriched in dust and metals. So far, however, the most distant UV-selected galaxy detected in dust emission is only at $z = 3.25$, and recent results have cast doubt on whether dust and molecules can be found in typical galaxies at this early epoch. Here we report thermal dust emission from an archetypal early universe star-forming galaxy, A1689-zD1. We detect its stellar continuum in spectroscopy and determine its redshift to be $z = 7.5\pm0.2$ from a spectroscopic detection of the Ly{\alpha} break. A1689-zD1 is representative of the star-forming population during reionisation, with a total star-formation rate of about 12M$_\odot$ yr$^{-1}$. The galaxy is highly evolved: it has a large stellar mass, and is heavily enriched in dust, with a dust-to-gas ratio close to that of the Milky Way. Dusty, evolved galaxies are thus present among the fainter star-forming population at $z > 7$, in spite of the very short time since they first appeared.
We investigate the evolution of the H$\beta$+[OIII] and [OII] luminosity functions from $z$ ~ 0.8 to ~ 5 in multiple redshift slices using data from the High-$z$ Emission Line Survey (HiZELS). This is the first time that the H$\beta$+[OIII] and [OII] luminosity functions have been studied at these redshifts in a self-consistent analysis. This is also the largest sample of [OII] and H$\beta$+[OIII] emitters (3484 and 3301 emitters, respectively) in this redshift range, with large co-moving volumes ~ $1 \times 10^6$ Mpc$^{3}$ in two independent volumes (COSMOS and UDS), greatly reducing the effects of cosmic variance. The emitters were selected by a combination of photometric redshift and color-color selections, as well as spectroscopic follow-up, including recent spectroscopic observations using DEIMOS and MOSFIRE on the Keck Telescopes and FMOS on Subaru. We find a strong increase in $L_\star$ and a decrease in $\phi_\star$ with increasing redshift up to $z \sim 2$ and $z \sim 5$ for H$\beta$+[OIII] and [OII] emitters, respectively. For H$\beta$+[OIII], this evolution then flattens by $z$ ~ 3. We derive the [OII] star-formation history of the Universe since $z$ ~ 5 and find that the cosmic SFRD rises from $z$ ~ 5 to ~ 3 and then drops towards $z$ ~ 0. We also find that our star-formation history is able to reproduce the evolution of the stellar mass density up to $z$ ~ 5. When comparing the H$\beta$+[OIII] SFRDs to the [OII] and H$\alpha$ SFRD measurements in the literature, we find that there is a remarkable agreement, suggesting that the H$\beta$+[OIII] sample is dominated by star-forming galaxies at high-$z$ rather than AGNs.
The formation of high mass stars and clusters occurs in giant molecular clouds. Objects in evolved stages of massive star formation such as protostars, hot molecular cores, and ultracompact HII regions have been studied in more detail than earlier, colder objects. With this in mind, the APEX Telescope Large Area Survey of the whole inner Galactic plane at 870 micron (ATLASGAL) has been carried out to provide a global view of cold dust and star formation at submillimetre wavelengths. To derive kinematic distances to a large sample of ATLASGAL clumps we divided them into groups of sources, which are located close together, mostly within a radius of 2 pc, and have velocities in a similar range with a median velocity dispersion of ~ 1 km/s. Using NH3, N2H+ and CS velocities we calculate near and far kinematic distances to 296 groups of ATLASGAL sources in the first quadrant and 393 groups in the fourth quadrant. We analyse HI self-absorption and HI absorption to resolve the kinematic distance ambiguity. We obtain a scale height of ~ 28+/-2 pc and displacement below the Galactic midplane of ~ -7+/-1 pc. Within distances from 2 to 18 kpc ATLASGAL clumps have a broad range of gas masses with a median of 1050 solar masses and a wide distribution of radii with a median of 0.4 pc. Their distribution in galactocentric radii is correlated with spiral arms. Using a statistically significant ATLASGAL sample we derive a power-law exponent of -2.2+/-0.1 of the clump mass function. This is consistent with the slope derived for clusters and with that of the stellar initial mass function. Examining the power-law index for different galactocentric distances and various source samples shows that it is independent of environment and evolutionary phase. Fitting the mass-size relationship by a power law gives a slope of 1.76+/-0.01 for cold sources such as IRDCs and warm clumps associated with HII regions.
New data are reported from the operation of a 2-liter C$_3$F$_8$ bubble chamber in the 2100 meter deep SNOLAB underground laboratory, with a total exposure of 211.5 kg-days at four different recoil energy thresholds ranging from 3.2 keV to 8.1 keV. These data show that C3F8 provides excellent electron recoil and alpha rejection capabilities at very low thresholds, including the first observation of a dependence of acoustic signal on alpha energy. Twelve single nuclear recoil event candidates were observed during the run. The candidate events exhibit timing characteristics that are not consistent with the hypothesis of a uniform time distribution, and no evidence for a dark matter signal is claimed. These data provide the most sensitive direct detection constraints on WIMP-proton spin-dependent scattering to date, with significant sensitivity at low WIMP masses for spin-independent WIMP-nucleon scattering.
This work explores the mathematical properties of a distribution introduced by Basu & Jones (2004), and applies it to model the stellar initial mass function (IMF). The distribution arises simply from an initial lognormal distribution, requiring that each object in it subsequently undergoes exponential growth but with an exponential distribution of growth lifetimes. This leads to a modified lognormal with a power-law tail (MLP) distribution, which can in fact be applied to a wide range of fields where distributions are observed to have a lognormal-like body and a power-law tail. We derive important properties of the MLP distribution, like the cumulative distribution, the mean, variance, arbitrary raw moments, and a random number generator. These analytic properties of the distribution can be used to facilitate application to modeling the IMF. We demonstrate how the MLP function provides an excellent fit to the IMF compiled by Chabrier (2005) and how this fit can be used to quickly identify quantities like the mean, median, and mode, as well as number and mass fractions in different mass intervals.
We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.
We construct the rest-frame 2--10 keV intrinsic X-ray luminosity function of Active Galactic Nuclei (AGNs) from a combination of X-ray surveys from the all-sky Swift BAT survey to the Chandra Deep Field-South. We use ~3200 AGNs in our analysis, which covers six orders of magnitude in flux. The inclusion of the XMM and Chandra COSMOS data has allowed us to investigate the detailed behavior of the XLF and evolution. In deriving our XLF, we take into account realistic AGN spectrum templates, absorption corrections, and probability density distributions in photometric redshift. We present an analytical expression for the overall behavior of the XLF in terms of the luminosity-dependent density evolution, smoothed two power-law expressions in 11 redshift shells, three-segment power-law expression of the number density evolution in four luminosity classes, and binned XLF. We observe a sudden flattening of the low luminosity end slope of the XLF slope at z>~0.6. Detailed structures of the AGN downsizing have been also revealed, where the number density curves have two clear breaks at all luminosity classes above log LX>43. The two break structure is suggestive of two-phase AGN evolution, consisting of major merger triggering and secular processes.
The Large Synoptic Survey Telescope (LSST) will photometrically monitor ~1 billion stars for ten years. The resulting light curves can be used to detect transiting exoplanets. In particular, as demonstrated by Lund et al. (2015), LSST will probe stellar populations currently undersampled in most exoplanet transit surveys, including out to extragalactic distances. In this paper we test the efficiency of the box-fitting least-squares (BLS) algorithm for accurately recovering the periods of transiting exoplanets using simulated LSST data. We model planets with a range of radii orbiting a solar-mass star at a distance of 7 kpc, with orbital periods ranging from 0.5 to 20 d. We find that typical LSST observations will be able to reliably detect Hot Jupiters with periods shorter than ~3 d. At the same time, we find that the LSST deep drilling cadence is extremely powerful: the BLS algorithm successfully recovers at least 30% of sub-Saturn-size exoplanets with orbital periods as long as 20 d.
We discuss methods for performing weak lensing using radio observations to recover information about the intrinsic structural properties of the source galaxies. Radio surveys provide unique information that can benefit weak lensing studies, such as HI emission, which may be used to construct galaxy velocity maps, and polarized synchrotron radiation; both of which provide information about the unlensed galaxy and can be used to reduce galaxy shape noise and the contribution of intrinsic alignments. Using a proxy for the intrinsic position angle of an observed galaxy, we develop techniques for cleanly separating weak gravitational lensing signals from intrinsic alignment contamination in forthcoming radio surveys. Random errors on the intrinsic orientation estimates introduce biases into the shear and intrinsic alignment estimates. However, we show that these biases can be corrected for if the error distribution is accurately known. We demonstrate our methods using simulations, where we reconstruct the shear and intrinsic alignment auto and cross-power spectra in three overlapping redshift bins. We find that the intrinsic position angle information can be used to successfully reconstruct both the lensing and intrinsic alignment power spectra with negligible residual bias.
We report the discovery of two transiting extrasolar planets by the HATSouth survey. HATS-9b orbits an old (10.8 $\pm$ 1.5 Gyr) V=13.3 G dwarf star, with a period P = 1.9153 d. The host star has a mass of 1.03 M$_{\odot}$, radius of 1.503 R$_\odot$ and effective temperature 5366 $\pm$ 70 K. The planetary companion has a mass of 0.837 M$_J$, and radius of 1.065 R$_J$ yielding a mean density of 0.85 g cm$^{-3}$ . HATS-10b orbits a V=13.1 G dwarf star, with a period P = 3.3128 d. The host star has a mass of 1.1 M$_\odot$, radius of 1.11 R$_\odot$ and effective temperature 5880 $\pm$ 120 K. The planetary companion has a mass of 0.53 M$_J$, and radius of 0.97 R$_J$ yielding a mean density of 0.7 g cm$^{-3}$ . Both planets are compact in comparison with planets receiving similar irradiation from their host stars, and lie in the nominal coordinates of Field 7 of K2 but only HATS-9b falls on working silicon. Future characterisation of HATS-9b with the exquisite photometric precision of the Kepler telescope may provide measurements of its reflected light signature.
RCW120 is a Galactic HII region which has a beautiful infrared ring. Previous studies on RCW120 provided a wealth of information on the second generation star formation, but the origin of the exciting O star located inside the ring structure has not been focused so far. Our new CO observations performed with the NANTEN2, Mopra, and ASTE telescopes have revealed that two molecular clouds with a velocity separation 20km/s are both physically associated with RCW120. The cloud at -8km/s apparently traces the infrared ring, while the other cloud at -28km/s is mainly distributed just outside the opening of the infrared ring, interacting with the HII region as supported by high kinetic temperature of the molecular gas and by the complementary distribution with the ionized gas. A spherically expanding shell driven by the HII region is usually discussed as the origin of the observed ring structure in RCW120. In this model, the neutral material which surrounds the HII region is expected to have an expanding motion. Our observations, however, indicate no evidence of the expanding motion in the velocity space, being inconsistent with the expanding shell scenario. We here postulate an alternative that, by applying the model introduced by Habe & Ohta (1992), the exciting O star in RCW120 was formed by a collision between the present two clouds at a colliding velocity of ~30km/s. In the model, the observed infrared ring can be interpreted as the cavity created in the larger cloud by the collision, whose inner surface is illuminated by the strong UV radiation after the birth of the O star. We argue that the present cloud-cloud collision scenario explains the observed signatures of RCW120, i.e., its ring morphology, coexistence of the two clouds and their large velocity separation, and absence of the expanding motion.
We present an analysis of elemental abundances of ejecta of the recurrent nova RS Oph using published optical and near-infrared spectra during the 2006 outburst. We use the CLOUDY photoionization code to generate synthetic spectra by varying several parameters, the model generated spectra are then matched with the observed emission line spectra obtained at two epochs. We obtain the best fit model parameters through the $\chi^{2}$ minimization technique. Our model results fit well with observed optical and near-infrared spectra. The best-fit model parameters are compatible with a hot white dwarf source with T$_{BB}$ of 5.5 - 5.8 $\times$ 10$^{5}$ K and roughly constant a luminosity of 6 - 8 $\times$ 10$^{36}$ ergs s$^{-1}$. From the analysis we find the following abundances (by number) of elements with respect to solar: He/H = 1.8 $\pm$ 0.1, N/H = 12.0 $\pm$ 1.0, O/H = 1.0 $\pm$ 0.4, Ne/H = 1.5 $\pm$ 0.1, Si/H = 0.4 $\pm$ 0.1, Fe/H = 3.2 $\pm$ 0.2, Ar/H = 5.1 $\pm$ 0.1, and Al/H = 1.0 $\pm$ 0.1, all other elements were set at the solar abundance. This shows the ejecta are significantly enhanced, relative to solar, in helium, nitrogen, neon, iron and argon. Using the obtained parameter values, we estimate an ejected mass in the range of 3.4 - 4.9 $\times$ 10$^{-6}$ M$_{\odot}$ which is consistent with observational results.
The brightness temperature of the radio free-free emission at millimeter range is an effective tool for characterizing the vertical structure of the solar chromosphere. In this paper, we report on the first single-dish observation of a sunspot at 85 and 115 GHz with sufficient spatial resolution for resolving the sunspot umbra using the Nobeyama 45 m telescope. We used radio attenuation material, i.e. a solar filter, to prevent the saturation of the receivers. Considering the contamination from the plage by the side-lobes, we found that the brightness temperature of the umbra should be lower than that of the quiet region. This result is inconsistent with the preexisting atmospheric models. We also found that the brightness temperature distribution at millimeter range strongly corresponds to the ultraviolet (UV) continuum emission at 1700 {\AA}, especially at the quiet region.
Since its launch in 2008 June, the Fermi Gamma-ray Space Telescope has opened
a new era in high-energy astrophysics. The unprecedented sensitivity, angular
resolution and effective area of the Large Area Telescope on board Fermi,
together with the nearly continuous observation of the entire gamma-ray sky
assures a formidable opportunity to study in detail gamma-ray emitting AGN of
various types. In this context the Swift satellite, thanks to its broad band
coverage and scheduling flexibility, creates a perfect synergy with Fermi.
Swift and Fermi coordinated monitoring campaigns of radio-loud AGN allowed us
to investigate correlated variability at different frequencies and to build
time-resolved spectral energy distributions from optical to gamma-rays,
constraining the emission mechanisms at work in these objects. The rapid Swift
follow-up observations of gamma-ray flaring AGN detected by Fermi-LAT were also
fundamental in firmly associating the gamma-ray sources with their low-energy
counterparts. We present some interesting results obtained from Fermi-LAT and
Swift observations of gamma-ray flaring AGN in the first six years of Fermi
operation.
The $\Lambda$CDM framework offers a remarkably good description of our
universe with a very small number of free parameters, which can be determined
with high accuracy from currently available data. However, this does not mean
that the associated physical quantities, such as the curvature of the universe,
have been directly measured. Similarly, general relativity is assumed, but not
tested. Testing the relevance of general relativity for cosmology at the
background level includes a verification of the relation between its energy
contents and the curvature of space. Using an extended Newtonian formulation,
we propose an approach where this relation can be tested. Using the recent
measurements on cosmic microwave background, baryonic acoustic oscillations and
the supernova Hubble diagram, we show that the prediction of general relativity
is well verified in the framework of standard $\Lambda$CDM assumptions, i.e. an
energy content only composed of matter and dark energy, in the form of a
cosmological constant or equivalently a vacuum contribution.
However, the actual equation of state of dark fluids cannot be directly
obtained from cosmological observations. We found that relaxing the equation of
state of dark energy opens a large region of possibilities, revealing a new
type of degeneracy between the curvature and the total energy content of the
universe.
The Burst Observer and Optical Transient Exploring System (BOOTES) is a network of telescopes that allows the continuous monitoring of transient astrophysical sources. It was originally devoted to the study of the optical emission from gamma-ray bursts (GRBs) that occur in the Universe. In this paper we show the initial results obtained using the spectrograph COLORES (mounted on BOOTES-2), when observing compact objects of diverse nature.
VLBI measurements of the absolute proper motions of 23 radio stars have been collected from published data. These are stars with maser emission, or very young stars, or asymptotic-giant-branch stars. By comparing these measurements with the stellar proper motions from the optical catalogs of the Hipparcos Celestial Reference Frame (HCRF), we have found the components of the residual rotation vector of this frame relative to the inertial coordinate system: (\omega_x,\omega_y,\omega_z) = (-0.39,-0.51,-1.25)+/-(0.58,0.57,0.56) mas/yr. Based on all the available data, we have determined new values of the components of the residual rotation vector for the optical realization of the HCRF relative to the inertial coordinate system: (\omega_x,\omega_y,\omega_z) = (-0.15,+0.24,-0.53)+/-(0.11,0.10,0.13) mas/yr.
Observations of low-mass stars reveal a variety of magnetic field topologies ranging from large-scale, axial dipoles to more complex magnetic fields. At the same time, three-dimensional spherical simulations of convectively driven dynamos reproduce a similar diversity, which is commonly obtained either with Boussinesq models or with more realistic models based on the anelastic approximation, which take into account the variation of the density with depth throughout the convection zone. Nevertheless, a conclusion from different anelastic studies is that dipolar solutions seem more difficult to obtain as soon as substantial stratifications are considered. In this paper, we aim at clarifying this point by investigating in more detail the influence of the density stratification on dipolar dynamos. To that end, we rely on a systematic parameter study that allows us to clearly follow the evolution of the stability domain of the dipolar branch as the density stratification is increased. The impact of the density stratification both on the dynamo onset and the dipole collapse is discussed and compared to previous Boussinesq results. Furthermore, our study indicates that the loss of the dipolar branch does not ensue from a specific modification of the dynamo mechanisms related to the background stratification, but could instead result from a bias as our observations naturally favour a certain domain in the parameter space characterized by moderate values of the Ekman number, owing to current computational limitations. Moreover, we also show that the critical magnetic Reynolds number of the dipolar branch is scarcely modified by the increase of the density stratification, which provides an important insight into the global understanding of the impact of the density stratification on the stability domain of the dipolar dynamo branch.
(shortened for arXiv) We aim to progress towards more efficient exoplanet detection around active stars by optimizing the use of Doppler Imaging in radial velocity measurements. We propose a simple method to simultaneously extract a brightness map and a set of orbital parameters through a tomographic inversion technique derived from classical Doppler mapping. Based on the maximum entropy principle, the underlying idea is to determine the set of orbital parameters that minimizes the information content of the resulting Doppler map. We carry out a set of numerical simulations to perform a preliminary assessment of the robustness of our method, using an actual Doppler map of the very active star HR 1099 to produce a realistic synthetic data set for various sets of orbital parameters of a single planet in a circular orbit. Using a simulated time-series of 50 line profiles affected by a peak-to-peak activity jitter of 2.5 km/s, we are able in most cases to recover the radial velocity amplitude, orbital phase and orbital period of an artificial planet down to a radial velocity semi-amplitude of the order of the radial velocity scatter due to the photon noise alone (about 50 m/s in our case). One noticeable exception occurs when the planetary orbit is close to co-rotation, in which case significant biases are observed in the reconstructed radial velocity amplitude, while the orbital period and phase remain robustly recovered. The present method constitutes a very simple way to extract orbital parameters from heavily distorted line profiles of active stars, when more classical radial velocity detection methods generally fail. It is easily adaptable to most existing Doppler Imaging codes, paving the way towards a systematic search for close-in planets orbiting young, rapidly-rotating stars.
Transient short-period <100s oscillations have been found in the X-ray light
curves of three novae during their SSS phase and in one persistent SSS. We
pursue an observational approach to determine possible driving mechanisms and
relations to fundamental system parameters such as the white dwarf mass.
We performed a systematic search for short-period oscillations in all
available XMM-Newton and Chandra X-ray light curves of persistent SSS and novae
during their SSS phase. To study time evolution, we divided each light curve
into short time segments and computed power spectra. We then constructed
dynamic power spectra from which we identified transient periodic signals even
when only present for a short time. From all time segments of each system, we
computed fractions of time when periodic signals were detected.
In addition to the previously known systems with short-period oscillations,
RS Oph (35s), KT Eri (35s), V339 Del (54s), and Cal 83 (67s), we found one
additional system, LMC 2009a (33s), and also confirm the 35s period from
Chandra data of KT Eri. The amplitudes of oscillations are of order <15% of the
respective count rates and vary without any clear dependence on the X-ray count
rate. The fractions of the time when the respective periods were detected at
2-sigma significance (duty cycle) are 11.3%, 38.8%, 16.9%, 49.2%, and 18.7% for
LMC 2009a, RS Oph, KT Eri, V339 Del, and Cal 83, respectively. The respective
highest duty cycles found in a single observation are 38.1%, 74.5%, 61.4%,
67.8%, and 61.8%.
We show that in the anticenter region, between Galactic longitudes of $110^\circ<l<229^\circ$, there is an oscillating asymmetry in the main sequence star counts on either side of the Galactic plane using data from the Sloan Digital Sky Survey. This asymmetry oscillates from more stars in the north at distances of about 2 kpc from the Sun to more stars in the south at 4-6 kpc from the Sun to more stars in the north at distances of 8-10 kpc from the Sun. We also see evidence that there are more stars in the south at distances of 12-16 kpc from the Sun. The three more distant asymmetries form roughly concentric rings around the Galactic center, opening in the direction of the Milky Way's spiral arms. The northern ring, 9 kpc from the Sun, is easily identified with the previously discovered Monoceros Ring. Parts of the southern ring at 14 kpc from the Sun (which we call the TriAnd Ring) have previously been identified as related to the Monoceros Ring and others have been called the Triangulum Andromeda Overdensity. The two nearer oscillations are approximated by a toy model in which the disk plane is offset by of the order 100 pc up and then down at different radii. We also show that the disk is not azimuthally symmetric around the Galactic anticenter and that there could be a correspondence between our observed oscillations and the spiral structure of the Galaxy. Our observations suggest that the TriAnd and Monoceros Rings (which extend to at least 25 kpc from the Galactic center) are primarily the result of disk oscillations.
We review five often used quad lens models, each of which has analytical solutions and can produce four images at most. Each lens model has two parameters, including one that describes the intensity of non-dimensional mass density, and the other one that describes the deviation from the circular lens. In our recent work, we have found that the cusp and the fold summations are not equal to 0, when a point source infinitely approaches a cusp or a fold from inner side of the caustic. Based on the magnification invariant theory, which states that the sum of signed magnifications of the total images of a given source is a constant, we calculate the cusp summations for the five lens models. We find that the cusp summations are always larger than 0 for source on the major cusps, while can be larger or smaller than 0 for source on the minor cusps. We also find that if these lenses tend to the circular lens, the major and minor cusp summations will have infinite values, and with positive and negative signs respectively. The cusp summations do not change significantly if the sources are slightly deviated from the cusps. In addition, through the magnification invariants, we also derive the analytical signed cusp relations on the axes for three lens models. We find that both on the major and the minor axes the larger the lenses deviated from the circular lens, the larger the signed cusp relations. The major cusp relations are usually larger than the absolute minor cusp relations, but for some lens models with very large deviation from circular lens, the minor cusp relations can be larger than the major cusp relations.
We present an analysis of the relation between star formation rate (SFR) surface density (sigmasfr) and mass surface density of molecular gas (sigmahtwo), commonly referred to as the Kennicutt-Schmidt (K-S) relation, at its intrinsic spatial scale, i.e. the size of giant molecular clouds (10-150 pc), in the central, high-density regions of four nearby low-luminosity active galactic nuclei (AGN). We used interferometric IRAM CO(1-0) and CO(2-1), and SMA CO(3-2) emission line maps to derive sigmahtwo and HST-Halpha images to estimate sigmasfr. Each galaxy is characterized by a distinct molecular SF relation at spatial scales between 20 to 200 pc. The K-S relations can be sub-linear, but also super-linear, with slopes ranging from 0.5 to 1.3. Depletion times range from 1 and 2Gyr, compatible with results for nearby normal galaxies. These findings are valid independently of which transition, CO(1-0), CO(2-1), or CO(3-2), is used to derive sigmahtwo. Because of star-formation feedback, life-time of clouds, turbulent cascade, or magnetic fields, the K-S relation might be expected to degrade on small spatial scales (<100 pc). However, we find no clear evidence for this, even on scales as small as 20 pc, and this might be because of the higher density of GMCs in galaxy centers which have to resist higher shear forces. The proportionality between sigmahtwo and sigmasfr found between 10 and 100 Msun/pc2 is valid even at high densities, 10^3 Msun/pc2. However, by adopting a common CO-to-H2 conversion factor (alpha_CO), the central regions of the galaxies have higher sigmasfr for a given gas column than those expected from the models, with a behavior that lies between the mergers/high-redshift starburst systems and the more quiescent star-forming galaxies, assuming that the first ones require a lower value of alpha_CO.
We model the abundance of haloes in the $\sim(3 \ \text{Gpc}/h)^3$ volume of the MICE Grand Challenge simulation by fitting the universal mass function with an improved Jack-Knife error covariance estimator that matches theory predictions. We present unifying relations between different fitting models and new predictions for linear ($b_1$) and non-linear ($c_2$ and $c_3$) halo clustering bias. Different mass function fits show strong variations in their overall poor performance when including the low mass range ($M_h \lesssim 3 \ 10^{12} \ M_{\odot}/h$) in the analysis, which indicates noisy friends-of-friends halo detection given the MICE resolution ($m_p \simeq 3 \ 10^{10} \ M_{\odot}$/h). Together with fits from the literature we find an overall variance in the amplitudes of around $10%$ in the low mass and up to $50%$ in the high mass (galaxy cluster) range ($M_h > 10^{14} \ M_{\odot}/h$). These variations propagate into a $10%$ change in $b_1$ predictions and a $50%$ change in $c_2$ or $c_3$. Despite these strong variations we find tight universal relations between $b_1$ and $c_2$ or $c_3$ for $b_1\gtrsim 1.5$ for which we provide simple fits. Their dependence on the mass function fit increases moderately for smaller $b_1$. Excluding low mass haloes, different models fitted with reasonable goodness in this analysis, show percent level agreement in their $b_1$ predictions, but are systematically $5-10%$ lower than the bias directly measured with two-point halo-mass clustering. This result confirms previous findings on larger volumes (and larger masses). Inaccuracies in the bias predictions propagate into the prediction of bias ratios at two redshifts, which would lead to $5-10%$ errors in growth measurements. They also affect any HOD fitting or (cluster) mass calibration from clustering measurements.
This work is aimed at quantifying the uncertainties in the 3D reconstruction of the location of coronal mass ejections (CMEs) obtained with the polarization ratio technique. The method takes advantage of the different distributions along the line of sight (LOS) of total (tB) and polarized (pB) brightnesses to estimate the average location of the emitting plasma. To this end, we assumed two simple electron density distributions along the LOS (a constant density and Gaussian density profiles) for a plasma blob and synthesized the expected tB and pB for different distances $z$ of the blob from the plane of the sky (POS) and different projected altitudes $\rho$. Reconstructed locations of the blob along the LOS were thus compared with the real ones, allowing a precise determination of uncertainties in the method. Independently of the analytical density profile, when the blob is centered at a small distance from the POS (i.e. for limb CMEs) the distance from the POS starts to be significantly overestimated. Polarization ratio technique provides the LOS position of the center of mass of what we call folded density distribution, given by reflecting and summing in front of the POS the fraction of density profile located behind that plane. On the other hand, when the blob is far from the POS, but with very small projected altitudes (i.e. for halo CMEs, $\rho < 1.4$ R$_\odot$), the inferred distance from that plane is significantly underestimated. Better determination of the real blob position along the LOS is given for intermediate locations, and in particular when the blob is centered at an angle of $20^\circ$ from the POS. These result have important consequences not only for future 3D reconstruction of CMEs with polarization ratio technique, but also for the design of future coronagraphs aimed at providing a continuous monitoring of halo-CMEs for space weather prediction purposes.
The possibility that star-forming galaxies may leak ionizing photons is at the heart of many present-day studies that investigate the reionization of the Universe. We test this hypothesis on local blue compact dwarf galaxies of very high excitation. We assembled a sample of such galaxies by examining the spectra from Data Releases 7 and 10 of the Sloan Digital Sky Survey. We argue that reliable conclusions cannot be based on strong lines alone, and adopt a strategy that includes important weak lines such as [OI] and the high-excitation HeII and [ArIV] lines. Our analysis is based on purely observational diagrams and on a comparison of photoionization models with well-chosen emission-line ratio diagrams. We show that spectral energy distributions from current stellar population synthesis models cannot account for all the observational constraints, which led us to mimick several scenarios that could explain the data. These include the additional presence of hard X-rays or of shocks. We find that only ionization-bounded models (or models with an escape fraction of ionizing photons lower than 10%) are able to simultaneously explain all the observational constraints.
The journey of the Sun through space carries the solar system through a dynamic interstellar environment that is presently characterized by Mach 1 motion between the heliosphere and the surrounding interstellar medium (ISM). The interaction between the heliosphere and ISM is an evolving process due to the variable solar wind and to interstellar turbulence. Frisch et al. presented a meta-analysis of the historical data on the interstellar wind flowing through the heliosphere and concluded that temporal changes in the ecliptic longitude of the wind were statistically indicated by the data available in the refereed literature at the time of that writing. Lallement and Bertaux disagree with this result, and suggested, for instance, that a key instrumental response function of IBEX-Lo was incorrect and that the STEREO pickup ion data are unsuitable for diagnosing the flow of interstellar neutrals through the heliosphere. Here we show that temporal variations in the interstellar wind through the heliosphere are consistent with our knowledge of local ISM. The statistical analysis of the historical helium wind data is revisited, and a recent correction of a typographical error in the literature is incorporated into the new fits. With this correction, and including no newer IBEX results, these combined data still indicate that a change in the longitude of the interstellar neutral wind over the past forty years is statistically likely, but that a constant flow longitude is now also statistically possible. It is shown that the IBEX instrumental response function is known, and that the STEREO pickup ion data have been correctly utilized in this analysis.
We present a physical basis for algorithms to replace mixing-length theory
(MLT) in stellar evolutionary computations. The 321D procedure is based on
three-dimensional (3D) time-dependent solutions of the Navier-Stokes equations,
including the Kolmogorov cascade as a sub-grid model of dissipation (implicit
large eddy simulations; ILES). We use Reynolds-averaged Navier-Stokes (RANS)
averaging to make 3D simulation data concise, and use 3D simulations to give
RANS closure. We sketch a simple algorithm, which is non-local and
time-dependent, with both MLT and the Lorenz convective roll as particular
subsets of solutions. The damping length is determined from a balance between
the large-scale driving and damping at the Kolmogorov scale.
We find that (1) braking regions (boundary layers in which mixing occurs)
automatically appear {\it beyond} the edges of convection as defined by the
Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero
turbulent kinetic energy flux (unlike MLT), (3) the effects of composition
gradients on flow are important, and (4) convective boundaries in
neutrino-cooled stages differ in nature from those in photon-cooled stages. The
321D approach may be easily generalized, and allows connections with modern
research on turbulent flow of solar and terrestrial fluids and plasmas.
Calibration to astronomical systems is unnecessary, so the approach can be
predictive rather than merely descriptive. Implications for solar abundances,
helioseismology, asteroseismology, nucleosynthesis yields, supernova
progenitors and core collapse are indicated.
We use the Planck LFI 70GHz data to further probe point source detection technique in the sky maps of the cosmic microwave background (CMB) radiation. The method developed by Tegmark et al. for foreground reduced maps and the Kolmogorov parameter as the descriptor are adopted for the analysis of Planck satellite CMB temperature data. Most of the detected points coincide with point sources already revealed by other methods. However, we have also found 9 source candidates for which still no counterparts are known.
Characterizing the local space density of double degenerate binary systems is a complementary approach to broad sky surveys of double degenerates to determine the expected rates of white dwarf binary mergers, in particular those that may evolve into other observable phenomena such as extreme helium stars, Am CVn systems, and supernovae Ia. However, there have been few such systems detected in local space. We report here the discovery that WD 1242$-$105, a nearby bright WD, is a double-line spectroscopic binary consisting of two degenerate DA white dwarfs of similar mass and temperature, despite it previously having been spectroscopically characterized as a single degenerate. Follow-up photometry, spectroscopy, and trigonometric parallax have been obtained in an effort to determine the fundamental parameters of each component of this system. The binary has a mass ratio of 0.7 and a trigonometric parallax of 25.5 mas, placing it at a distance of 39 pc. The system's total mass is 0.95 M$_\odot$ and has an orbital period of 2.85 hours, making it the strongest known gravitational wave source ($\log h = -20.78$) in the mHz regime. Because of its orbital period and total mass, WD 1242$-$105 is predicted to merge via gravitational radiation on a timescale of 740 Myr, which will most likely not result in a catastrophic explosion.
A recently discovered filament of polarized starlight that traces a coherent magnetic field is shown to have several properties that are consistent with an origin in the outer heliosheath of the heliosphere: (1) The magnetic field that provides the best fit to the polarization position angles is directed within 6.7+-11 degrees of the observed upwind direction of the flow of interstellar neutral helium gas through the heliosphere. (2) The magnetic field is ordered; the component of the variation of the polarization position angles that can be attributed to magnetic turbulence is small. (3) The axis of the elongated filament can be approximated by a line that defines an angle of 80+/-14 degrees with the plane that is formed by the interstellar magnetic field vector and the vector of the inflowing neutral gas (the "BV" plane). We propose that this polarization feature arises from aligned interstellar dust grains in the outer heliosheath where the interstellar plasma and magnetic field are deflected around the heliosphere. The proposed outer heliosheath location of the polarizing grains requires confirmation by modeling grain-propagation through three-dimensional MHD heliosphere models that simultaneously calculate torques on asymmetric dust grains interacting with the heliosphere.
In an effort to better understand the details of the stellar structure and evolution of metal poor stars, the Gemini North telescope was used on two occasions to take speckle imaging data of a sample of known spectroscopic binary stars and other nearby stars in order to search for and resolve close companions. The observations were obtained using the Differential Speckle Survey Instrument, which takes data in two filters simultaneously. The results presented here are of 90 observations of 23 systems in which one or more companions was detected, and 6 stars where no companion was detected to the limit of the camera capabilities at Gemini. In the case of the binary and multiple stars, these results are then further analyzed to make first orbit determinations in five cases, and orbit refinements in four other cases. Mass information is derived, and since the systems span a range in metallicity, a study is presented that compares our results with the expected trend in total mass as derived from the most recent Yale isochrones as a function of metal abundance. These data suggest that metal-poor main-sequence stars are less massive at a given color than their solar-metallicity analogues in a manner consistent with that predicted from the theory.
We examine the rise and sudden demise of comet C/2010 X1 (Elenin) on its approach to perihelion. Discovered inbound at 4.2 AU, this long-period comet was predicted to become very bright when near perihelion, at 0.48 AU on 2011 September 10. Observations starting 2011 February (heliocentric distance $\sim$3.5 AU) indeed show the comet to brighten by about 11 magnitudes, with most of the increase occurring inside 1 AU from the Sun. The peak brightness reached $m_R$ = 6 on UT 2011 August 12.95$\pm$0.50, when at $\sim$0.83 AU from the Sun. We find that most of the surge in brightness in mid-August resulted from dust particle forward-scattering, not from a sudden increase in the activity. A much smaller ($\sim$3 magnitudes) brightening reached a maximum on UT 2011 August 30$\pm$1 (at 0.56 AU), and reflects the true break-up of the nucleus. This second peak was matched by a change in the morphology from centrally condensed to diffuse. The estimated cross-section of the nucleus when at 1 AU inbound was $\sim$1 km$^2$, corresponding to an equal-area circle of radius 0.6 km. No surviving fragments were found to a limiting red magnitude $r'$ = 24.4, corresponding to radii $\lesssim$40 m (red geometric albedo = 0.04 assumed). Our observations are consistent with disintegration of the nucleus into a power law size distribution of fragments with index $q$ = 3.3$\pm$0.2 combined with the action of radiation pressure. We speculate about physical processes that might cause nucleus disruption in a comet when still 0.7 AU from the Sun. Tidal stresses and devolatilization of the nucleus by sublimation are both negligible at this distance. However, the torque caused by mass loss, even at the very low rates measured in comet Elenin, is potentially large enough to be responsible by driving the nucleus to rotational instability.
The HDO/H2O ratio in interstellar gas is often used to draw conclusions on the origin of water in star-forming regions and on Earth. In cold cores and in the outer regions of protoplanetary disks, gas-phase water comes from photodesorption of water ice. We present fitting formulae for implementation in astrochemical models using photodesorption efficiencies for all water ice isotopologues obtained using classical molecular dynamics simulations. We investigate if the gas-phase HDO/H2O ratio reflects that present in the ice or whether fractionation can occur during photodesorption. Probabilities for the top four monolayers are presented for photodesorption of X (X=H,D) atoms, OX radicals, and X2O and HDO molecules following photodissociation of H2O, D2O, and HDO in H2O amorphous ice at temperatures from 10-100 K. Isotope effects are found for all products: (1) H atom photodesorption probabilities from H2O ice are larger than those for D atom photodesorption from D2O ice by a factor of 1.1; the ratio of H and D photodesorbed upon HDO photodissociation is a factor of 2. This process will enrich the ice in deuterium atoms over time; (2) the OD/OH photodesorption ratio upon D2O and H2O photodissociation is on average a factor of 2, but the ratio upon HDO photodissociation is almost constant at unity for all temperatures; (3) D atoms are more effective in kicking out neighbouring water molecules than H atoms. However, the ratio of the photodesorbed HDO and H2O molecules is equal to the HDO/H2O ratio in the ice, therefore, there is no isotope fractionation upon HDO and H2O photodesorption. Nevertheless, the enrichment of the ice in D atoms due to photodesorption can over time lead to an enhanced HDO/H2O ratio in the ice, and, when photodesorbed, also in the gas. The extent to which the ortho/para ratio of H2O can be modified by the photodesorption process is also discussed. (Abridged)
NGC 5194 (M51a) is a grand-design spiral galaxy and undergoing interactions with its companion. Here we focus on investigating main properties of its star-formation history (SFH) by constructing a simple evolution model, which assumes that the disc builds up gradually by cold gas infall and the gas infall rate can be parameterizedly described by a Gaussian form. By comparing model predictions with the observed data, we discuss the probable range for free parameter in the model and then know more about the main properties of the evolution and SFH of M51a. We find that the model predictions are very sensitive to the free parameter and the model adopting a constant infall-peak time $t_{\rm p}\,=\,7.0{\rm Gyr}$ can reproduce most of the observed constraints of M51a. Although our model does not assume the gas infall time-scale of the inner disc is shorter than that of the outer disc, our model predictions still show that the disc of M51a forms inside-out. We find that the mean stellar age of M51a is younger than that of the Milky Way, but older than that of the gas-rich disc galaxy UGC 8802. In this paper, we also introduce a 'toy' model to allow an additional cold gas infall occurred recently to imitate the influence of the interaction between M51a and its companion. Our results show that the current molecular gas surface density, the SFR and the UV-band surface brightness are important quantities to trace the effects of recent interaction on galactic SF process.
The Dark Matter Particle Explorer (DAMPE) is an upcoming scientific satellite mission for high energy gamma-ray, electron and cosmic rays detection. The silicon tracker (STK) is a sub detector of the DAMPE payload with an excellent position resolution (readout pitch of 242um), which measures the incident direction of particles, as well as charge. The STK consists 12 layers of Silicon Micro-strip Detector (SMD), equivalent to a total silicon area of 6.5m$^2$. The total readout channels of the STK are 73728, which leads to a huge amount of raw data to be dealt. In this paper, we focus on the on-board data compression algorithm and procedure in the STK, which was initially verified by cosmic-ray measurements.
Warped accretion disks have attracted intensive attention because of their critical role on shaping the spin of supermassive massive black holes (SMBHs) through the Bardeen-Petterson effect, a general relativistic effect that leads to final alignments or anti-alignments between black holes and warped accretion disks. We study such alignment processes by explicitly taking into account the finite sizes of accretion disks and the episodic lifetimes of AGNs that delineate the duration of gas fueling onto accretion disks. We employ an approximate global model to simulate the evolution of accretion disks, allowing to determine the gravitomagnetic torque that drives the alignments in a quite simple way. We then track down the evolutionary paths for mass and spin of black holes both in a single activity episode and over a series of episodes. Given with randomly and isotropically oriented gas fueling over episodes, we calculate the spin evolution with different episodic lifetimes and find that it is quite sensitive to the lifetimes. We therefore propose that spin distribution of SMBHs can place constraints on the episodic lifetimes of AGNs and vice versa. Applications of our results on the observed spin distributions of SMBHs and the observed episodic lifetimes of AGNs are discussed, although both the measurements at present are yet ambiguous to draw a firm conclusion. Our prescription can be easily incorporated into semi-analytic models for black hole growth and spin evolution.
The production of high-energy astrophysical neutrinos is tightly linked to the emission of hadronic gamma-rays. I will discuss the recent observation of TeV to PeV neutrinos by the IceCube Cherenkov telescope in the context of gamma-ray astronomy. The corresponding energy range of hadronic gamma-rays is not directly accessible by extragalactic gamma-ray astronomy due to interactions with cosmic radiation backgrounds. Nevertheless, the isotropic sub-TeV gamma-ray background observed by the Fermi Large Area Telescope (LAT) contains indirect information from secondary emission produced in electromagnetic cascades and constrains hadronic emission scenarios. On the other hand, observation of PeV gamma-rays would provide a smoking-gun signal for Galactic emission. In general, the cross-correlation of neutrino emission with (extended) Galactic and extragalactic gamma-ray sources will serve as the most sensitive probe for a future identification of neutrino sources.
The relation between long gamma-ray bursts (LGRBs) and low-luminosity GRBs (llgrbs) is a long standing puzzle -- on the one hand their high energy emission properties are fundamentally different, implying a different gamma-ray source, yet both are associated with similar supernovae of the same peculiar type (broad-line Ic), pointing at a similar progenitor and a similar explosion mechanism. Here we analyze the multi-wavelength data of the particularly well-observed SN 2006aj, associated with llgrb 060218, finding that its progenitor star is sheathed in an extended ($>100R_\odot$), low-mass ($\sim 0.01M_\odot$) envelope. This progenitor structure implies that the gamma-ray emission in this llgrb is generated by a mildly relativistic shock breakout. It also suggests a unified picture for llgrbs and LGRBs, where the key difference is the existence of an extended low-mass envelope in llgrbs and its absence in LGRBs. The same engine, which launches a relativistic jet, can drive the two explosions, but, while in LGRBs the ultra-relativistic jet emerges from the bare progenitor star and produces the observed gamma-rays, in llgrbs the extended envelope smothers the jet and prevents the generation of a large gamma-ray luminosity. Instead, the jet deposits all its energy in the envelope, driving a mildly relativistic shock that upon breakout produces a llgrb. In addition for giving a unified view of the two phenomena, this model provides a natural explanation to many observed properties of llgrbs. It also implies that llgrbs are a viable source of the observed extra-galactic diffuse neutrino flux and that they are promising sources for future gravitational wave detectors.
High redshift blazars are among the most powerful objects in the Universe. Although they represent a significant fraction of the extragalactic hard X-ray sky, they are not commonly detected in gamma-rays. High redshift (z>2) objects represent <10 per cent of the AGN population observed by Fermi so far, and gamma-ray flaring activity from these sources is even more uncommon. The characterization of the radio-to-gamma-ray properties of high redshift blazars represent a powerful tool for the study of both the energetics of such extreme objects and the Extragalactic Background Light. We present results of a multi-band campaign on TXS 0536+145, which is the highest redshift flaring gamma-ray blazar detected so far. At the peak of the flare the source reached an apparent isotropic gamma-ray luminosity of 6.6x10^49 erg/s, which is comparable with the luminosity observed from the most powerful blazars. The physical properties derived from the multi-wavelength observations are then compared with those shown by the high redshift population. In addition preliminary results from the high redshift flaring blazar PKS 2149-306 will be discussed.
Electrons/positrons produced in a pulsar magnetosphere emit synchrotron radiation, which is widely believed as the origin of the non-thermal X-ray emission detected from pulsars. Particles are produced by curvature photons emitted from accelerated particles in the magnetosphere. These curvature photons are detected as pulsed $\gamma$-ray emissions from pulsars with age $\lesssim10^6$ yr. Using $\gamma$-ray observations and analytical model, we impose severe constraints on the synchrotron radiation as a mechanism of the non-thermal X-ray emission. In most middle-aged pulsars ($\sim10^5-10^6$ yr) which photon-photon pair production is less efficient in their magnetosphere, we find that the synchrotron radiation model is difficult to explain the observed non-thermal X-ray emission.
Using the DR12 public release of APOGEE data, we show that thin and thick disk separate very well in the space defined by [$\alpha$/Fe], [Fe/H] and [C/N]. Thick disk giants have both higher [C/N] and higher [$\alpha$/Fe] than do thin disk stars with similar [Fe/H]. We deduce that the thick disk is composed of lower mass stars than the thin disk. Considering the fact that at a given metallicity there is a one-to-one relation between stellar mass and age, we are then able to infer the chronology of disk formation. Both the thick and the thin disks - defined by [$\alpha$/Fe] -- converge in their dependance on [C/N] and [C+N/Fe] at [Fe/H]$\approx$-0.7. We conclude that 1) the majority of thick disk stars formed earlier than did the thin disk stars 2) the formation histories of the thin and thick disks diverged early on, even when the [Fe/H] abundances are similar 3) that the star formation rate in the thin disk has been lower than in the thick disk, at all metallicities. Although these general conclusions remain robust, we also show that current stellar evolution models cannot reproduce the observed C/N ratios for thick disk stars. Unexpectedly, reduced or inhibited canonical extra-mixing is very common in field stars. While subject to abundance calibration zeropoint uncertainties, this implies a strong dependence of non canonical extra-mixing along the red giant branch on the initial composition of the star and in particular on the $\alpha$ elemental abundance.
Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.
Context. Various Be stars exhibit intensity variations of the violet and red
emission peaks in their HI lines observed in emission. This so-called $V/R$
phenomenon is usually explained by the precession of a one-armed spiral density
perturbation in the circumstellar disk. That global-disk oscillation scenario
was confirmed, both observationally and theoretically, in the previous series
of two papers analyzing the Be shell star {\zeta} Tauri. The vertically
averaged (2D) global-disk oscillation model used at the time was able to
reproduce the $V/R$ variations observed in H{\alpha}, as well as the spatially
resolved interferometric data from AMBER/VLTI. Unfortunately, that model failed
to reproduce the $V/R$ phase of Br15 and the amplitude of the polarization
variation, suggesting that the inner disk structure predicted by the model was
incorrect.
Aims. The first aim of the present paper is to quantify the temporal
variations of the shell-line characteristics of {\zeta} Tauri. The second aim
is to better understand the physics underlying the $V/R$ phenomenon by modeling
the shell-line variations together with the $V/R$ and polarimetric variations.
The third aim is to test a new 2.5D disk oscillation model, which solves the
set of equations that describe the 3D perturbed disk structure but keeps only
the equatorial (i.e., 2D) component of the solution. This approximation was
adopted to allow comparisons with the previous 2D model, and as a first step
toward a future 3D model.
Results. The new 2.5D formalism improves the agreement with the observed
$V/R$ variations of H{\alpha} and Br15, under the proviso that a large value of
the viscosity parameter, {\alpha} = 0.8, be adopted. Nonetheless, it remains
challenging for the models to reproduce consistently the amplitude and the
average level of the polarization data, whatever formalism is adopted.
Post equinox imaging of Uranus by HST, Keck, and Gemini telescopes has enabled new measurements of winds over previously sampled latitudes as well as measurements at high northern latitudes that have recently come into better view. These new observations also used techniques to greatly improve signal to noise ratios, making possible the detection and tracking of more subtle cloud features. The 250 m/s prograde jet peaking near 60 N was confirmed and more accurately characterized. Several long-lived cloud features have also been tracked. The winds pole-ward of 60 N are consistent with solid body rotation at a westward (prograde) rate of 4.3 deg/h with respect to Uranus' interior. When combined with 2007 and other recent measurements, it is clear that a small but well-resolved asymmetry exists in the zonal profile at middle latitudes, peaking at 35 deg, where southern winds are 20 m/s more westward than corresponding northern winds. High S/N Keck II imaging of the north polar region of Uranus reveals a transition from streaky bands below 60 N to a region from 60 deg to nearly the north pole, where widely distributed small bright spots, resembling cumulus cloud fields, with several isolated dark spots, are the dominant style of cloud features. This presents a stark contrast to 2003 detailed views of the south polar region of Uranus when no discrete cloud features could be detected in comparable Keck II near-IR images. The pressure levels of discrete clouds estimated from spatial modulations in H and Hcont images indicate that the polar cloud features are generally in the 1.3 to 2-3 bar range, as are equatorial and several mid-latitude features. Several of the brighter mid latitude features are found above the 1.2-bar level of methane condensation.
In this paper searches for flaring astrophysical neutrino sources and sources with periodic emission with the IceCube neutrino telescope are presented. In contrast to time integrated searches, where steady emission is assumed, the analyses presented here look for a time dependent signal of neutrinos using the information from the neutrino arrival times to enhance the discovery potential. A search was performed for correlations between neutrino arrival times and directions as well as neutrino emission following time dependent lightcurves, sporadic emission or periodicities of candidate sources. These include active galactic nuclei, soft $\gamma$-ray repeaters, supernova remnants hosting pulsars, micro-quasars and X-ray binaries. The work presented here updates and extends previously published results to a longer period that covers four years of data from 2008 April 5 to 2012 May 16 including the first year of operation of the completed 86-string detector. The analyses did not find any significant time dependent point sources of neutrinos and the results were used to set upper limits on the neutrino flux from source candidates.
On 2015 March 20, a total solar eclipse will occur in the North Atlantic,
with the Kingdom of Denmark's Faroe Islands and Norway's Svalbard archipelago
(formerly Spitzbergen) being the only options for land-based observing. The
region is known for wild, unpredictable, and often cloudy conditions, which
potentially pose a serious threat for people hoping to view the spectacle.
We report on a citizen-science, weather-monitoring project, based in the
Faroe Islands, which was conducted in March 2014 - one year prior to the
eclipse. The project aimed to promote awareness of the eclipse among the local
communities, with the data collected providing a quantitative overview of
typical weather conditions that may be expected in 2015. It also allows us to
validate the usefulness of short-term weather forecasts, which may be used to
increase the probability of observing the eclipse.
In order to investigate the characteristics and influence of the magnetic field in evolved stars, we performed a follow-up investigation of our previous submillimeter analysis of the proto-planetary nebula (PPN) OH 231.8+4.2 (Sabin et al. 2014), this time at 1.3mm with the CARMA facility in polarisation mode for the purpose of a multi-scale analysis. OH 231.8+4.2 was observed at ~2.5" resolution and we detected polarised emission above the 3-sigma threshold (with a mean polarisation fraction of 3.5 %). The polarisation map indicates an overall organised magnetic field within the nebula. The main finding in this paper is the presence of a structure mostly compatible with an ordered toroidal component that is aligned with the PPN's dark lane. We also present some alternative magnetic field configuration to explain the structure observed. These data complete our previous SMA submillimeter data for a better investigation and understanding of the magnetic field structure in OH 231.8+4.2.
We present a procedure for efficiently compressing astronomical radio data for high performance applications. Integrated, post-correlation data are first passed through a nearly lossless rounding step which compares the precision of the data to a generalized and calibration-independent form of the radiometer equation. This allows the precision of the data to be reduced in a way that has an insignificant impact on the data. The newly developed Bitshuffle lossless compression algorithm is subsequently applied. When the algorithm is used in conjunction with the HDF5 library and data format, data produced by the CHIME Pathfinder telescope is compressed to 28% of its original size and decompression throughputs in excess of 1 GB/s are obtained on a single core.
We describe a set of simulations of super-critical accretion onto a non-rotating supermassive BH. The accretion flow is radiation pressure dominated and takes the form of a geometrically thick disk with twin low-density funnels around the rotation axis. For accretion rates $\gtrsim 10 \dot M_{\rm Edd}$, there is sufficient gas in the funnel to make this region optically thick. Radiation from the disk first flows into the funnel, after which it accelerates the optically thick funnel gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted velocity $\approx 0.3c$. Much of the radiative luminosity is converted into kinetic energy by the time the escaping gas becomes optically thin. For an observer viewing down the axis, the isotropic equivalent luminosity of total energy is as much as $10^{48}\,\rm erg\,s^{-1}$ for a $10^7 M_\odot$ BH accreting at $10^3$ Eddington. Therefore, energetically, the simulated jets are consistent with observations of the most powerful tidal disruption events, e.g., Swift J1644. The jet velocity is, however, too low to match the Lorentz factor $\gamma > 2$ inferred in J1644. There is no such conflict in the case of other tidal disruption events. Since favorably oriented observers see isotropic equivalent luminosities that are highly super-Eddington, the simulated models can explain observations of ultra-luminous X-ray sources, at least in terms of luminosity and energetics, without requiring intermediate mass black holes. Finally, since the simulated jets are baryon-loaded and have mildly relativistic velocities, they match well the jets observed in SS433. The latter are, however, more collimated than the simulated jets. This suggests that, even if magnetic fields are not important for acceleration, they may perhaps still play a role in confining the jet.
Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a search for variability using well sampled JHKs-band light curves from the 2MASS survey calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 candidates. All 7 were young stellar objects within the {\rho} Ophiuchus star forming region, five with previous X-ray detections. A significant {\gamma}-ray detection (5{\sigma}) for the known blazar using 4.5 years of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified {\gamma}-ray blazars, with low contamination from other AGN.
Decades of observations show that CMEs can deflect from a purely radial trajectory yet no consensus exists as to the cause of these deflections. Many of theories attribute the CME deflection to magnetic forces. We developed ForeCAT (Kay et al. 2013, Kay et al. 2015), a model for CME deflections based solely on magnetic forces, neglecting any reconnection effects. Here we compare ForeCAT predictions to the observed deflection of the 2008 December 12 CME and find that ForeCAT can accurately reproduce the observations. Multiple observations show that this CME deflected nearly 30{\deg} in latitude (Byrne et al. 2010, Gui et al. 2011) and 4.4{\deg} in longitude (Gui et al. 2011). From the observations, we are able to constrain all of the ForeCAT input parameters (initial position, radial propagation speed, and expansion) except the CME mass and the drag coefficient that affects the CME motion. By minimizing the reduced chi-squared, $\chi^2_{\nu}$, between the ForeCAT results and the observations we determine an acceptable mass range between 4.5x10$^{14}$ and 1x10$^{15}$ g and the drag coefficient less than 1.4 with a best fit at 7.5x10$^{14}$ g and 0 for the mass and drag coefficient. ForeCAT is sensitive to the magnetic background and we are also able to constrain the rate at which the quiet sun magnetic field falls to be similar or to or fall slightly slower than the Potential Field Source Surface model.
The K2 mission has recently begun to discover new and diverse planetary systems. In December 2014 Campaign 1 data from the mission was released, providing high-precision photometry for ~22000 objects over an 80 day timespan. We searched these data with the aim of detecting further important new objects. Our search through two separate pipelines led to the independent discovery of EPIC201505350, a two-planet system of Neptune sized objects (4.2 and 7.2 $R_\oplus$), orbiting a K dwarf extremely close to the 3:2 mean motion resonance. The two planets each show transits, sometimes simultaneously due to their proximity to resonance and alignment of conjunctions. We obtain further ground based photometry of the larger planet with the NITES telescope, demonstrating the presence of large transit timing variations (TTVs) of over an hour. These TTVs allows us to confirm the planetary nature of the system, and place a limit on the mass of the outer planet of $386M_\oplus$.
Do We Inhabit The Best O All Possible Worlds? German mathematician Gottfried Leibniz thought so, writing in 1710 that our planet, warts and all, must be the most optimal one imaginable. Leibniz's idea was roundly scorned as unscientific wishful thinking, most notably by French author Voltaire in his magnum opus, Candide. Yet Leibniz might find sympathy from at least one group of scientists - the astronomers who have for decades treated Earth as a golden standard as they search for worlds beyond our own solar system. Because earthlings still know of just one living world - our own - it makes some sense to use Earth as a template in the search for life elsewhere, such as in the most Earth-like regions of Mars or Jupiter's watery moon Europa. Now, however, discoveries of potentially habitable planets orbiting stars other than our sun - exoplanets, that is - are challenging that geocentric approach.
The methane-nitrogen phase diagram of Prokhvatilov and Yantsevich (1983) indicates that at temperatures relevant to the surfaces of icy dwarf planets like Pluto, two phases contribute to the methane absorptions: nitrogen saturated with methane $\bf{\bar{N_{2}}}$:CH$_{4}$ and methane saturated with nitrogen $\bf{\bar{CH_{4}}}$:N$_{2}$. No optical constants are available so far for the latter component limiting construction of a proper model, in compliance with thermodynamic equilibrium considerations. New optical constants for solid solutions of methane diluted in nitrogen (N$_{2}$:CH$_{4}$) and nitrogen diluted in methane (CH$_{4}$:N$_{2}$) are presented at temperatures between 40 and 90 K, in the wavelength range 1.1-2.7 $\mu$m at different mixing ratios. These optical constants are derived from transmission measurements of crystals grown from the liquid phase in closed cells. A systematic study of the changes of methane and nitrogen solid mixtures spectral behavior with mixing ratio and temperature is presented.
Scalar-tensor Cosmologies can be dealt under the standard of the Hojman conservation theorem that allows to fix the form of the coupling $F(\phi)$, of the potential $V(\phi)$ and to find out exact solutions for related cosmological models. Specifically, the existence of a symmetry transformation vector for the equations of motion gives rise to a Hojman conserved quantity on the corresponding minisuperpace and exact solutions for the cosmic scale factor $a$ and the scalar field $\phi$ can be achieved. In particular, we take advantage of the fact that minimally coupled solutions, previously obtained in the Einstein frame, can be conformally transformed in non-minimally coupled solutions in the Jordan frame. Some physically relevant examples are worked out.
A technique to estimate mass erosion rate of surface soil during landing of the Apollo Lunar Module (LM) and total mass ejected due to the rocket plume interaction is proposed and tested. The erosion rate is proportional to the product of the second moment of the lofted particle size distribution N(D), and third moment of the normalized soil size distribution S(D), divided by the integral of S(D)D^2/v(D), where D is particle diameter and v(D) is the vertical component of particle velocity. The second moment of N(D) is estimated by optical extinction analysis of the Apollo cockpit video. Because of the similarity between mass erosion rate of soil as measured by optical extinction and rainfall rate as measured by radar reflectivity, traditional NWS radar/rainfall correlation methodology can be applied to the lunar soil case where various S(D) models are assumed corresponding to specific lunar sites.
We study the asymptotic scaling properties of domain wall networks with three different tensions in various cosmological epochs. We discuss the conditions under which a scale-invariant evolution of the network (which is well established for simpler walls) still applies, and also consider the limiting case where defects are locally planar and the curvature is concentrated in the junctions. We present detailed quantitative predictions for scaling densities in various contexts, which should be testable by means of future high-resolution numerical simulations.
In the light of the history of researches on electromagnetic wave spectrum, a sharp emission line of gravitational-wave background (GWB) would be an interesting observational target. Here we study an efficient method to detect a line GWB by correlating data of multiple ground-based detectors. We find that the width of frequency bin for coarse graining is a critical parameter, and the commonly-used value 0.25 Hz is far from optimal, decreasing the signal-to-noise ratio by up to a factor of seven. By reanalyzing the existing data with a smaller bin width, we might detect a precious line signal from the early universe.
The neutron-proton effective mass splitting in neutron-rich nucleonic matter reflects the space-time nonlocality of the isovector nuclear interaction. It affects the neutron/proton ratio during the earlier evolution of the Universe, cooling of protoneutron stars, structure of rare isotopes and dynamics of heavy-ion collisions. While there is still no consensus on whether the neutron-proton effective mass splitting is negative, zero or positive and how it depends on the density as well as the isospin-asymmetry of the medium, significant progress has been made in recent yeas in addressing these issues. We first recall the connections among the neutron-proton effective mass splitting, the momentum dependence of the isovector potential and the density dependence of the symmetry energy. We then make a few observations about the progress in calculating the neutron-proton effective mass splitting using various nuclear many-body theories and its effects on the isospin-dependence of in-medium nucleon-nucleon cross sections. Perhaps, our most reliable knowledge so far about the neutron-proton effective mass splitting at saturation density of nuclear matter comes from optical model analyses of huge sets of nucleon-nucleus scattering data accumulated over the last five decades. The momentum dependence of the symmetry potential from these analyses provide a useful boundary condition at saturation density for calibrating nuclear many-body calculations. Several observables in heavy-ion collisions have been identified as sensitive probes of the neutron-proton effective mass splitting in dense neutron-rich matter based on transport model simulations. We review these observables and comment on the latest experimental findings.
We review inflationary cosmology in modified gravity such as $R^2$ gravity with its extensions in order to generalize the Starobinsky inflation model. In particular, we explore inflation realized by three kinds of effects: modification of gravity, the quantum anomaly, and the $R^2$ term in loop quantum cosmology. It is explicitly demonstrated that in these inflationary models, the spectral index of scalar modes of the density perturbations and the tensor-to-scalar ratio can be consistent with the Planck results. Bounce cosmology in $F(R)$ gravity is also explained.
The main background above 3\,MeV for in-beam nuclear astrophysics studies with $\gamma$-ray detectors is caused by cosmic-ray induced secondaries. The two commonly used suppression methods, active and passive shielding, against this kind of background were formerly considered only as alternatives in nuclear astrophysics experiments. In this work the study of the effects of active shielding against cosmic-ray induced events at a medium deep location is performed. Background spectra were recorded with two actively shielded HPGe detectors. The experiment was located at 148\,m below the surface of the Earth in the Reiche Zeche mine in Freiberg, Germany. The results are compared to data with the same detectors at the Earth's surface, and at depths of 45\,m and 1400\,m, respectively.
In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the solar wind, however, with a hyperbolic flux-tube or "$X$-line" extending over integral length-scales.
Cluster formation is a fundamental aspect of the equation of state (EOS) of warm and dense nuclear matter such as can be found in supernovae (SN). Similar matter can be studied in heavy-ion collisions (HIC). We use the experimental data of Qin et al. 2012 to test calculations of cluster formation and the role of in-medium modifications of cluster properties in SN EOSs. For the comparison between theory and experiment we use chemical equilibrium constants as the main observables. This reduces some of the systematic uncertainties and allows deviations from ideal gas behavior to be identified clearly. In the analysis, we carefully account for the differences between matter in SN and HIC. We find that, at the lowest densities, the experiment and all theoretical models are consistent with the ideal gas behavior. At higher densities ideal behavior is clearly ruled out and interaction effects have to be considered. The contributions of continuum correlations are of relevance in the virial expansion and remain a difficult problem to solve at higher densities. We conclude that at the densities and temperatures discussed mean-field interactions of nucleons, inclusion of all relevant light clusters, and a suppression mechanism of clusters at high densities have to be incorporated in the SN EOS.
An efficient method for calculating inclusive conventional and prompt atmospheric leptons fluxes is presented. The coupled cascade equations are solved numerically by formulating them as matrix equation. The presented approach is very flexible and allows the use of different hadronic interaction models, realistic parametrizations of the primary cosmic-ray flux and the Earth's atmosphere, and a detailed treatment of particle interactions and decays. The power of the developed method is illustrated by calculating lepton flux predictions for a number of different scenarios.
We investigate the prospects of indirect and direct dark matter searches within the minimal supersymmetric standard model with nine parameters (MSSM-9). These nine parameters include three gaugino masses, Higgs, slepton and squark masses, all treated independently. We perform a Bayesian Monte Carlo scan of the parameter space taking into consideration all available particle physics constraints such as the Higgs mass of 126 GeV, upper limits on the scattering cross-section from direct-detection experiments, and assuming that the MSSM-9 provides all the dark matter abundance through thermal freeze-out mechanism. Within this framework we find two most probable regions for dark matter: 1-TeV higgsino-like and 3-TeV wino-like neutralinos. We discuss prospects for future indirect (in particular the Cherenkov Telescope Array, CTA) and direct detection experiments. We find that for slightly contracted dark matter profiles in our Galaxy, which can be caused by the effects of baryonic infall in the Galactic center, CTA will be able to probe a large fraction of the remaining allowed region in synergy with future direct detection experiments like XENON-1T.
We predict the new galvano-rotational effect which consists in the appearance of the electric current along the axis of the matter rotation. This effect is caused by the parity violating electroweak interaction between massless charged particles in the rotating matter. We start with the exact solution of the Dirac equation for a fermion involved in the electroweak interaction in the rotating frame. This equation includes the noninertial effects. Then, using the obtained solution, we derive the induced electric current which turns out to flow along the rotation axis. We study the possibility of the appearance of the galvano-rotational effect in dense matter of compact astrophysical objects. The particular example of neutron and hypothetical quark stars is discussed. It is shown that, using this effect, one can expect the generation of toroidal magnetic fields comparable with poloidal ones in old millisecond pulsars. We also briefly discuss the generation of the magnetic helicity in these stars.
Particle velocity distributions measured in the weakly collisional solar wind are frequently found to be non-Maxwellian, but how these non-Maxwellian distributions impact the physics of plasma turbulence in the solar wind remains unanswered. Using numerical solutions of the linear dispersion relation for a collisionless plasma with a bi-Maxwellian proton velocity distribution, we present a unified framework for the four proton temperature anisotropy instabilities, identifying the associated stable eigenmodes, highlighting the unstable region of wavevector space, and presenting the properties of the growing eigenfunctions. Based on physical intuition gained from this framework, we address how the proton temperature anisotropy impacts the nonlinear dynamics of the \Alfvenic fluctuations underlying the dominant cascade of energy from large to small scales and how the fluctuations driven by proton temperature anisotropy instabilities interact nonlinearly with each other and with the fluctuations of the large-scale cascade. We find that the nonlinear dynamics of the large-scale cascade is insensitive to the proton temperature anisotropy, and that the instability-driven fluctuations are unlikely to cause significant nonlinear evolution of either the instability-driven fluctuations or the turbulent fluctuations of the large-scale cascade.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 M_Sun), highly eccentric (e ~ 0.9) binary star system Eta Carinae. We demonstrate the methodology used to incorporate 3D interactive figures into a PDF journal publication and the benefits of using 3D visualization and 3D printing as tools to analyze data from multidimensional numerical simulations. Using a consumer-grade 3D printer (MakerBot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics (SPH) simulations of Eta Carinae's inner (r ~ 110 au) wind-wind collision interface at multiple orbital phases. The 3D prints and visualizations reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (phi ~ 1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. thin-shell, Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the fast (3000 km/s), adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unrecognized physical features highlight the important role 3D printing and interactive graphics can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.
Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.
We report the discovery of five Local Volume dwarf galaxies uncovered during a comprehensive archival search for optical counterparts to ultra-compact high velocity clouds (UCHVCs). The UCHVC population of HI clouds are thought to be candidate gas-rich, low mass halos at the edge of the Local Group and beyond, but no comprehensive search for stellar counterparts to these systems has been presented. Careful visual inspection of all publicly available optical and ultraviolet imaging at the position of the UCHVCs revealed six blue, diffuse counterparts with a morphology consistent with a faint dwarf galaxy beyond the Local Group. Optical spectroscopy of all six candidate dwarf counterparts show that five have an H$\alpha$-derived velocity consistent with the coincident HI cloud, confirming their association; the sixth diffuse counterpart is likely a background object. The size and luminosity of the UCHVC dwarfs is consistent with other known Local Volume dwarf irregular galaxies. The gas fraction ($M_{HI}/M_{star}$) of the five dwarfs are generally consistent with that of dwarf irregular galaxies in the Local Volume, although ALFALFA-Dw1 (associated with ALFALFA UCHVC HVC274.68+74.70$-$123) has a very high $M_{HI}/M_{star}$$\sim$40. Despite the heterogenous nature of our search, we demonstrate that the current dwarf companions to UCHVCs are at the edge of detectability due to their low surface brightness, and that deeper searches are likely to find more stellar systems. If more sensitive searches do not reveal further stellar counterparts to UCHVCs, then the dearth of such systems around the Local Group may be in conflict with $\Lambda$CDM simulations.
We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$).
We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.
We present a large-scale, volume-limited companion survey of 245 late-K to mid-M (K7-M6) dwarfs within 15 pc. Infrared adaptive optics (AO) data were analysed from the Very Large Telescope, Subaru Telescope, Canada-France-Hawaii Telescope, and MMT Observatory to detect close companions to the sample from $\sim$1 au to 100 au, while digitised wide-field archival plates were searched for wide companions from $\sim$100 au to 10,000 au. With sensitivity to the bottom of the main sequence over a separation range of 3 au to 10,000 au, multiple AO and wide-field epochs allow us to confirm candidates with common proper motions, minimize background contamination, and enable a measurement of comprehensive binary statistics. We detected 65 co-moving stellar companions and find a companion star fraction of $23.5 \pm 3.2$ per cent over the 3 au to 10,000 au separation range. The companion separation distribution is observed to rise to a higher frequency at smaller separations, peaking at closer separations than measured for more massive primaries. The mass ratio distribution across the $q = 0.2 - 1.0$ range is flat, similar to that of multiple systems with solar-type primaries. The characterisation of binary and multiple star frequency for low-mass field stars can provide crucial comparisons with star forming environments and hold implications for the frequency and evolutionary histories of their associated disks and planets.
We examine the radial distributions of stellar populations in the globular cluster (GC) M15, using HST/WFC3 photometry of red giants in the nitrogen-sensitive F343N-F555W color. Surprisingly, we find that giants with "primordial" composition (i.e., N abundances similar to those in field stars) are the most centrally concentrated within the WFC3 field. We then combine our WFC3 data with SDSS u, g photometry and find that the trend reverses for radii >1' (3 pc) where the ratio of primordial to N-enhanced giants increases outwards, as already found by Lardo et al. The ratio of primordial to enriched stars thus has a U-shaped dependency on radius with a minimum near the half-light radius. N-body simulations show that mass segregation might produce a trend resembling the observed one, but only if the N-enhanced giants are ~0.25 Mo less massive than the primordial giants, which requires extreme He enhancement (Y~0.40). However, such a large difference in Y is incompatible with the negligible optical color differences between primordial and enriched giants which suggest Delta Y < 0.03 and thus a difference in turn-off mass of Delta M < 0.04 Mo between the different populations. The radial trends in M15 are thus unlikely to be of dynamical origin and presumably reflect initial conditions, a result that challenges all current GC formation scenarios. We note that population gradients in the central regions of GCs remain poorly investigated and may show a more diverse behavior than hitherto thought.
Numerical simulations indicate that black holes carrying linear momentum and/or orbital momentum can power jets. The jets extract the kinetic energy stored in the black hole's motion. This could provide an important electromagnetic counterpart to gravitational wave searches. We develop the theory underlying these jets. In particular, we derive the analogues of the Penrose process and the Blandford-Znajek jet power prediction for boosted black holes. The jet power we find is $(v/2M)^2 \Phi^2/(4\pi)$, where $v$ is the hole's velocity, $M$ is its mass, and $\Phi$ is the magnetic flux. We show that energy extraction from boosted black holes is conceptually similar to energy extraction from spinning black holes. However, we highlight two key technical differences: in the boosted case, jet power is no longer defined with respect to a Killing vector, and the relevant notion of black hole mass is observer-dependent. We derive a new version of the membrane paradigm in which the membrane lives at infinity rather than the horizon and we show that this is useful for interpreting jets from boosted black holes. Our jet power prediction and the assumptions behind it can be tested with future numerical simulations.
We propose a decision criterion for segmenting the cosmic web into different structure types (voids, sheets, filaments and clusters) on the basis of their respective probabilities and the strength of data constraints. Our approach is inspired by an analysis of games of chance where the gambler only plays if a positive expected net gain can be achieved based on some degree of privileged information. The result is a general solution for classification problems in the face of uncertainty, including the option of not committing to a class for a candidate object. As an illustration, we produce high-resolution maps of web-type constituents in the nearby Universe as probed by the Sloan Digital Sky Survey main galaxy sample. Other possible applications include the selection and labeling of objects in catalogs derived from astronomical survey data.
We utilize $\Lambda$CDM halo occupation models of galaxy clustering to investigate the evolving stellar mass dependent clustering of galaxies in the PRIsm MUlti-object Survey (PRIMUS) and DEEP2 Redshift Survey over the past eight billion years of cosmic time, between $0.2<z<1.2$. These clustering measurements provide new constraints on the connections between dark matter halo properties and galaxy properties in the context of the evolving large-scale structure of the universe. Using both an analytic model and a set of mock galaxy catalogs, we find a strong correlation between central galaxy stellar mass and dark matter halo mass over the range $M_\mathrm{halo}\sim10^{11}$-$10^{13}~h^{-1}M_\odot$, approximately consistent with previous observations and theoretical predictions. However, the stellar-to-halo mass relation (SHMR) and the mass scale where star formation efficiency reaches a maximum appear to evolve more strongly than predicted by other models. We find that the fraction of satellite galaxies in haloes of a given mass decreases by $\approx5\%$ from $z\sim0.5$ to $z\sim0.9$, and we find that the $M_1/M_\mathrm{min}$ ratio, which quantifies the critical mass above which haloes host at least one satellite, decreases from $\approx20$ at $z\sim0$ to $\approx13$ at $z\sim0.9$. Considering the evolution of the subhalo mass function vis-\`{a}-vis satellite abundances, this trend has implications for relations between satellite galaxies and halo substructures and for intracluster mass, which we argue has grown due to stripped and disrupted satellites between $z\sim0.9$ and $z\sim0.5$.
Hydrodynamical simulations of galaxy formation such as the Illustris simulations have progressed to a state where they approximately reproduce the observed stellar mass function from high to low redshift. This in principle allows self-consistent models of reionization that exploit the accurate representation of the diffuse gas distribution together with the realistic growth of galaxies provided by these simulations, within a representative cosmological volume. In this work, we apply and compare two radiative transfer algorithms implemented in a GPU-accelerated code to the $106.5\,{\rm Mpc}$ wide volume of Illustris in postprocessing in order to investigate the reionization transition predicted by this model. We find that the first generation of galaxies formed by Illustris is just about able to reionize the universe by redshift $z\sim 7$, provided quite optimistic assumptions about the escape fraction and the resolution limitations are made. Our most optimistic model finds an optical depth of $\tau\simeq 0.065$, which is in very good agreement with recent Planck 2015 determinations. Furthermore, we show that moment-based approaches for radiative transfer with the M1 closure give broadly consistent results with our angular-resolved radiative transfer scheme as far as the global reionization history is concerned. We also confirm earlier findings that the reduced speed-of-light approximation introduces non-neglibible inaccuracies. In our favoured fiducial model, 20% of the hydrogen is reionized by redshift $z=9.20$, and this rapidly climbs to 80% by redshift $z=6.92$. It then takes until $z=6.24$ before 99% of the hydrogen is ionized. On average, reionization proceeds `inside-out' in our models, with a size distribution of reionized bubbles that progressively features regions of ever larger size while the abundance of small bubbles stays fairly constant.
In the cosmological paradigm, cold dark matter (DM) dominates the mass content of the Universe and is present at every scale. Candidates for DM include many extensions of the standard model, with a weakly interacting massive particle (WIMP) in the mass range from $\sim$10 GeV to greater than 10 TeV. The self-annihilation or decay of WIMPs in astrophysical regions of high DM density can produce secondary particles including very high energy (VHE) gamma rays with energy up to the DM particle mass. VERITAS, an array of atmospheric Cherenkov telescopes, sensitive to VHE gamma rays in the 85 GeV - 30 TeV energy range, has been utilized for DM searches. The possible astrophysical objects considered to be candidates for indirect DM detection are VERITAS dwarf spheroidal galaxies (dSphs) of the Local Group and the Galactic Center, among others. This presentation reports on our extensive observations of these targets and constraints of the dark matter physics from these objects, including the methodology and preliminary results of a combined DM search of five dSphs.
We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. (2014), comparing its high energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light-cylinder, is assumed to be negligible, while outside the light-cylinder it rescales with a finite conductivity ({\sigma}). In our approach we calculate the corresponding high energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young {\gamma}-ray pulsars, including features of the phase resolved spectra. Second we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The {\sigma} values that best describe each of the pulsars in this group show an increase with the spin-down rate $(\dot{E})$ and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.
We implement physically motivated recipes for partitioning cold gas into different phases (atomic, molecular, and ionized) in galaxies within semi-analytic models of galaxy formation based on cosmological merger trees. We then model the conversion of molecular gas into stars using empirical recipes motivated by recent observations. We explore the impact of these new recipes on the evolution of fundamental galaxy properties such as stellar mass, star formation rate (SFR), and gas and stellar phase metallicity. We present predictions for stellar mass functions, stellar mass vs. SFR relations, and cold gas phase and stellar mass-metallicity relations for our fiducial models, from redshift $z\sim 6$ to the present day. In addition we present predictions for the global SFR, mass assembly history, and cosmic enrichment history. We find that the predicted stellar properties of galaxies (stellar mass, SFR, metallicity) are remarkably insensitive to the details of the recipes used for partitioning gas into HI and H$_2$. We see significant sensitivity to the recipes for H$_2$ formation only in very low mass halos, which host galaxies that are not detectable with current observational facilities except very nearby. The properties of low-mass galaxies are also quite insensitive to the details of the recipe used for converting H$_2$ into stars, while the formation epoch of massive galaxies does depend on this significantly. (Abridged)
Accurate and repeated photometric follow-up observations of planetary-transit events are important to precisely characterize the physical properties of exoplanets. A good knowledge of the main characteristics of the exoplanets is fundamental to trace their origin and evolution. Multi-band photometric observations play an important role in this process. By using new photometric data, we computed precise estimates of the physical properties of two transiting planetary systems. We present new broad-band, multi-colour, photometric observations obtained using three small class telescopes and the telescope-defocussing technique. For each of the two targets, one transit event was simultaneously observed through four optical filters. One transit of WASP-48 b was monitored with two telescopes from the same observatory. The physical parameters of the systems were obtained by fitting the transit light curves with {\sc jktebop} and from published spectroscopic measurements. We have revised the physical parameters of the two planetary systems, finding a smaller radius for both HAT-P-23 b and WASP-48 b, $R_{b}=1.224 \pm 0.037 R_{Jup}$ and $R_{b}=1.396 \pm 0.051 \, R_{Jup}$, respectively, than those measured in the discovery papers ($R_{b}=1.368 \pm 0.090 R_{Jup}$ and $R_{b}=1.67 \pm 0.10 R_{Jup}$). The density of the two planets are higher than those previously published ($\rho_{b}$ ~1.1 and ~0.3 $\rho_{jup}$ for HAT-P-23 and WASP-48 respectively) hence the two Hot Jupiters are no longer located in a parameter space region of highly inflated planets. An analysis of the variation of the planet's measured radius as a function of optical wavelength reveals flat transmission spectra within the experimental uncertainties. We also confirm the presence of the eclipsing contact binary NSVS-3071474 in the same field of view of WASP-48, for which we refine the value of the period to be 0.459 d.
In the upcoming synoptic all--sky survey era of astronomy, thousands of new multiply imaged quasars are expected to be discovered and monitored regularly. Light curves from the images of gravitationally lensed quasars are further affected by superimposed variability due to microlensing. In order to disentangle the microlensing from the intrinsic variability of the light curves, the time delays between the multiple images have to be accurately measured. The resulting microlensing light curves can then be analyzed to reveal information about the background source, such as the size of the quasar accretion disc. In this paper we present the most extensive and coherent collection of simulated microlensing light curves; we have generated $>2.5$ billion light curves using the GERLUMPH high resolution microlensing magnification maps. Our simulations can be used to: train algorithms to measure lensed quasar time delays, plan future monitoring campaigns, and study light curve properties throughout parameter space. Our data are openly available to the community and are complemented by online eResearch tools, located at this http URL .
We present an empirical model for the halo evolution and global gas dynamics of Fornax, the brightest Milky Way (MW) dwarf spheroidal galaxy (dSph). Assuming a global star formation rate psi(t)=lambda_*[M_g(t)/M_sun]^alpha consistent with observations of star formation in nearby galaxies and using the data on Fornax's psi(t), we derive the evolution of the total mass M_g(t) for cold gas in Fornax's star-forming disk and the rate Delta F(t) of net gas flow to or from the disk. We identify the onset of the transition in Delta F(t) from a net inflow to a net outflow as the time t_sat at which the Fornax halo became an MW satellite and estimate the evolution of its total mass M_h(t) at t<t_sat using the median halo growth history in the LambdaCDM cosmology and its present mass within the half-light radius derived from observations. We examine three different cases of alpha=1, 1.5, and 2, and justify the corresponding lambda_* by comparing the gas mass fraction f_g(t)=M_g(t)/M_h(t) with results from simulations of gas accretion by halos in a reionized universe. We find that the Fornax halo grew to M_h(t_sat)=1.8x10^9M_sun at t_sat= 4.8 Gyr and that its subsequent global gas dynamics was dominated by ram-pressure stripping and tidal interaction with the MW. Gas loss over a few orbital periods eventually terminated its star formation. Our approach can be extended to other dSphs and provide input for detailed studies of their formation and evolution.
Recent observations of tidal debris around galaxies have revealed that the structural properties of the spheroidal components of tidally disturbed galaxies are similar to those found in non-interacting early-type galaxies(ETGs), likely due to minor merging events that do not strongly affect the bulge region or to major mergers that happened a long time ago. We show that independently of merger events, tidal features like shells or rings can also arise if the the dark matter is an ultra light scalar field of mass ~10$^{-22}$eV/c$^2$. In the scalar field dark matter (SFDM) model the small mass precludes halo formation below ~10$^8$ M$_{\odot}$ reducing the number of small galaxies today, it produces shallow density profiles due to the uncertainty principle in contrast to the steep profiles found in the standard cold dark matter (CDM) paradigm, in addition to the usual soliton solution there exists dark matter haloes in multistates, characterized by ripples in their density profiles, which are stable provided that most of the halo mass resides in the ground state. We use the hydrodynamics code ZEUS to track the gas evolution in a background potential given by a superposition of the ground and first excited state of the scalar field, we study this configuration when it is initially unstable (excited state more massive than ground state) but by a population inversion in the states it eventually becomes stable, this could happen when haloes decoupled from the expansion of the universe and collapse to reach a state of equilibrium. We found that tidal structures like rings are formed at a particular radii as a direct consequence of the wavelike structure of the dark matter halo(abridged)
We performed single point [C I] $^3$P$_1$-$^3$P$_0$ and CO J=4-3 observations toward three T Tauri stars, DM Tau, LkCa 15, and TW Hya, using the Atacama Large Millimeter/submillimeter Array (ALMA) Band 8 qualification model receiver installed on the Atacama Submillimeter Telescope Experiment (ASTE). Two protostars in the Taurus L1551 region, L1551 IRS 5 and HL Tau, were also observed. We successfully detected [C I] emission from the protoplanetary disk around DM Tau as well as the protostellar targets. The spectral profile of the [C I] emission from the protoplanetary disk is marginally single-peaked, suggesting that atomic carbon (C) extends toward the outermost disk. The detected [C I] emission is optically thin and the column densities of C are estimated to be <~10$^{16}$ cm$^{-2}$ and ~10$^{17}$ cm$^{-2}$ for the T Tauri star targets and the protostars, respectively. We found a clear difference in the total mass ratio of C to dust, $M$(C)/$M$(dust), between the T Tauri stars and protostellar targets; the $M$(C)/$M$(dust) ratio of the T Tauri stars is one order of magnitude smaller than that of the protostars. The decrease of the estimated $M$(C)/$M$(dust) ratios for the disk sources is consistent with a theoretical prediction that the atomic C can survive only in the near surface layer of the disk and C$^+$/C/CO transition occurs deeper into the disk midplane.
We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are treated as rotating spherical shells filled with stratified gas. The hydrodynamic equations are discretized by a robust and efficient high-order Spectral Difference Method (SDM) on unstructured meshes. The computational stencil of the spectral difference method is compact and advantageous for parallel processing. CHORUS demonstrates excellent parallel performance for all test cases reported in this paper, scaling up to 12,000 cores on the Yellowstone High-Performance Computing cluster at NCAR. The code is verified by defining two benchmark cases for global convection in Jupiter and the Sun. CHORUS results are compared with results from the ASH code and good agreement is found. The CHORUS code creates new opportunities for simulating such varied phenomena as multi-scale solar convection, core convection, and convection in rapidly-rotating, oblate stars.
We report on third epoch VLBI observations of the radio-bright supernova SN 2011dh located in the nearby galaxy (7.8 Mpc) M51. The observations took place at $t=453$ d after the explosion and at a frequency of 8.4 GHz. We obtained a fairly well resolved image of the shell of SN 2011dh, making it one of only six recent supernovae for which resolved images of the ejecta are available. By fitting a spherical shell model directly to the visibility measurements we determine the angular radius of SN 2011dh's radio emission to be $636 \pm 29$ $\mu$as . At a distance of 7.8 Mpc, this angular radius corresponds to a linear radius of $(7.4 \pm 0.3) \times 10^{16}$ cm and an average expansion velocity since the explosion of $18900^{+2800}_{-2400}$ kms$^{-1}$. We also calculated more precise radius measurements for the earlier VLBI observations and we show that all the measured values of the radius of the emission region, up to $t=453$ d, are still almost perfectly consistent with those derived from fitting synchrotron self-absorbed models to the radio spectral energy distribution. We find that SN 2011dh's radius evolves in a power-law fashion, with $R \sim t^{0.961\pm0.011}$, implying almost free expansion.
We adopt an isotropic representation of the Fourier-transformed two-point correlation tensor of the magnetic field for estimating magnetic energy and helicity spectra as well as current helicity spectra of individual active regions and the change of their spectral indices with the solar cycle. The departure of the spectral index of current helicity from 5/3 is analyzed, and it is found that it is lower than that of magnetic energy. There is no obvious relationship between the change of the normalized magnetic helicity and the integral scale of the magnetic field for individual active regions. The evolution of the spectral index reflects the development and distribution of various scales of magnetic structures in active regions. It is found that around solar maximum the magnetic energy and helicity spectra are steeper.
Observations of an unprecedented quality made by ALMA on the Red Rectangle of CO(3-2) and CO(6-5) emissions are analysed jointly with the aim of obtaining as simple as possible a description of the gas morphology and kinematics. Evidence is found for polar conical outflows and for a broad equatorial torus in rotation and expansion. Simple models of both are proposed. Comparing CO(6-5) and CO(3-2) emissions provides evidence for a strong temperature enhancement over the polar outflows. Continuum emission (dust) is seen to be enhanced in the equatorial region. Observed asymmetries are briefly discussed.
Type II radio bursts are thought to be a signature of coronal shocks. In this paper, we analyze a short-lived type II burst that started at 07:40 UT on 2011 February 28. By carefully checking white-light images, we find that the type II radio burst is not accompanied by a coronal mass ejection, only with a C2.4 class flare and narrow jet. However, in the extreme-ultraviolet (EUV) images provided by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), we find a wave-like structure that propagated at a speed of $\sim$ 600 km s$^{-1}$ during the burst. The relationship between the type II radio burst and the wave-like structure is in particular explored. For this purpose, we first derive the density distribution under the wave by the differential emission measure (DEM) method, which is used to restrict the empirical density model. We then use the restricted density model to invert the speed of the shock that produces the observed frequency drift rate in the dynamic spectrum. The inverted shock speed is similar to the speed of the wave-like structure. This implies that the wave-like structure is most likely a coronal shock that produces the type II radio burst. We also examine the evolution of the magnetic field in the flare-associated active region and find continuous flux emergence and cancellation taking place near the flare site. Based on these facts, we propose a new mechanism for the formation of the type II radio burst, i.e., the expansion of the strongly-inclined magnetic loops after reconnected with nearby emerging flux acts as a piston to generate the shock wave.
In order to study the origin of the architectures of low mass planetary systems, we perform numerical surveys of the evolution of pairs of coplanar planets in the mass range $(1-4)\ \rmn{M}_{\oplus}.$ These evolve for up to $2\times10^7 \rmn{yr}$ under a range of orbital migration torques and circularization rates assumed to arise through interaction with a protoplanetary disc. Near the inner disc boundary, significant variations of viscosity, interaction with density waves or with the stellar magnetic field could occur and halt migration, but allow ircularization to continue. This was modelled by modifying the migration and circularization rates. Runs terminated without an extended period of circularization in the absence of migration torques gave rise to either a collision, or a system close to a resonance. These were mostly first order with a few $\%$ terminating in second order resonances. Both planetary eccentricities were small $< 0.1$ and all resonant angles liberated. This type of survey produced only a limited range of period ratios and cannot reproduce Kepler observations. When circularization alone operates in the final stages, divergent migration occurs causing period ratios to increase. Depending on its strength the whole period ratio range between $1$ and $2$ can be obtained. A few systems close to second order commensurabilities also occur. In contrast to when arising through convergent migration, resonant trapping does not occur and resonant angles circulate. Thus the behaviour of the resonant angles may indicate the form of migration that led to near resonance.
We present cosmological bounds on the thermal axion mass in an extended cosmological scenario in which the primordial power spectrum of scalar perturbations differs from the usual power-law shape predicted by the simplest inflationary models. The power spectrum is instead modeled by means of a "piecewise cubic Hermite interpolating polynomial" (PCHIP). When using Cosmic Microwave Background measurements combined with other cosmological data sets, the thermal axion mass constraints are degraded only slightly. The addition of the measurements of $\sigma_8$ and $\Omega_m$ from the 2013 Planck cluster catalogue on galaxy number counts relaxes the bounds on the thermal axion mass, mildly favouring a $\sim 1$~eV axion mass, regardless of the model adopted for the primordial power spectrum.
In the last few decades many efforts have been made to understand the effect of spiral arms on the gas and stellar dynamics in the Milky Way disc. One of the fundamental parameters of the spiral structure is its angular velocity, or pattern speed $\Omega_p$, which determines the location of resonances in the disc and the spirals' radial extent. The most direct method for estimating the pattern speed relies on backward integration techniques, trying to locate the stellar birthplace of open clusters. Here we propose a new method based on the interaction between the spiral arms and the stars in the disc. Using a sample of around 500 open clusters from the {\it New Catalogue of Optically Visible Open Clusters and Candidates}, and a sample of 500 giant stars observed by APOGEE, we find $\Omega_p = 23.0\pm0.5$ km s$^{-1}$ kpc$^{-1}$, for a local standard of rest rotation $V_0=220$~km s$^{-1}$ and solar radius $R_0=8.0$~kpc. Exploring a range in $V_0$ and $R_0$ within the acceptable values, 200-240 km s$^{-1}$ and 7.5-8.5 kpc, respectively, results only in a small change in our estimate of $\Omega_p$, that is within the error. Our result is in close agreement with a number of studies which suggest values in the range 20-25 km s$^{-1}$ kpc$^{-1}$. An advantage of our method is that we do not need knowledge of the stellar age, unlike in the case of the birthplace method, which allows us to use data from large Galactic surveys. The precision of our method will be improved once larger samples of disk stars with spectroscopic information will become available thanks to future surveys such as 4MOST.
The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. We find VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select HeI lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass in excess of 25Msun for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating `single' O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.
W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". Our observational results do not support cloud-cloud collision scenario as claimed by previous observations. The blueshifted absorption of HCO$^{+}$ (3-2) is in contrast to the global collapse scenario in this region, which was proposed based on the previous detection of redshifted absorption features in HCO$^{+}$ (1-0) and CS (2-1) lines. The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. The outflow jet is in precession, which might account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating a cold layer is expanding in front of the background continuum source. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.
We construct equilibrium configurations of magnetized, two-fluid neutron stars using an iterative numerical method. We assume that the neutron star has two regions: the core, which is modelled as a two-component fluid consisting of type-II superconducting protons and superfluid neutrons, and the crust, a region composed of normal matter. Taking a new step towards more complete equilibrium models, we include the effect of entrainment, which implies that a magnetic force acts on neutrons, too. We consider purely poloidal field cases and present improvements to an earlier numerical scheme for solving equilibrium equations, by introducing new convergence criteria. We find that entrainment results in qualitative differences in the structure of field lines across the crust-core boundary and along the magnetic axis.
We report the results of a multi-waveband analysis of the masses and luminosities of $\sim$600 galaxy groups and clusters identified in the maxBCG catalogue. These data are intended to form the basis of future work on the formation of the "$m_{12}$ gap" in galaxy groups and clusters. We use SDSS spectroscopy and $g$, $r$ and $i$ band photometry to estimate galaxy group/cluster virial radii, masses and total luminosities. In order to establish the robustness of our results, we compare them with literature studies that utilize a variety of mass determinations techniques (dynamical, X-ray, weak lensing) and total luminosities estimated in the $B$, $r$, $i$, and $K$ wavebands. We also compare our results to predictions derived from the Millennium Simulation. We find that, once selection effects are properly accounted for, excellent agreement exists between our results and the literature with the exception of a single observational study. We also find that the Millennium Simulation does an excellent job of predicting the effects of our selection criteria. Our results show that, over the mass range $\sim10^{13}-10^{15}$ M$_{\odot}$, variations in the slope of the mass-luminosity scaling relation with mass detected in this and many other literature studies is in part the result of selection effects. We show that this can have serious ramifications on attempts to determine how the mass-to-light ratio of galaxy groups and cluster varies with mass.
We study an extensive sample of 87 GRBs for which there are well sampled and simultaneous optical and X-ray light-curves. We extract the cleanest possible signal of the afterglow component, and compare the temporal behaviors of the X-ray light-curve, observed by Swift XRT, and optical data, observed by UVOT and ground-based telescopes for each individual burst. Overall we find 62\% GRBs that are consistent with the standard afterglow model. When more advanced modeling is invoked, up to 91\% of the bursts in our sample may be consistent with the external shock model. A large fraction of these bursts are consistent with occurring in a constant interstellar density medium (ISM) (61\%) while only 39\% of them occur in a wind-like medium. Only 9 cases have afterglow light-curves that exactly match the standard fireball model prediction, having a single power law decay in both energy bands which are observed during their entire duration. In particular, for the bursts with chromatic behavior additional model assumptions must be made over limited segments of the light-curves in order for these bursts to fully agree with the external shock model. Interestingly, for 54\% of the X-ray and 40\% of the optical band observations the end of the shallow decay ($t^{\sim-0.5}$) period coincides with the jet break ($t^{\sim-p}$) time, causing an abrupt change in decay slope. The fraction of the burst that consistent with the external shock model is independent of the observational epochs in the rest frame of GRBs. Moreover, no cases can be explained by the cooling frequency crossing the X-ray or optical band.
Solar X-ray spectra from the RESIK crystal spectrometer on the {\em CORONAS-F} spacecraft (spectral range $3.3-6.1$~\AA) are analyzed for thirty-three flares using a method to derive abundances of Si, S, Ar, and K, emission lines of which feature prominently in the spectra. For each spectrum, the method first optimizes element abundances then derives the differential emission measure as a function of temperature based on a procedure given by Sylwester et al. and Withbroe. This contrasts with our previous analyses of RESIK spectra in which an isothermal assumption was used. The revised abundances (on a logarithmic scale with $A({\rm H}) = 12$) averaged for all the flares in the analysis are $A({\rm Si}) = 7.53 \pm 0.08$ (previously $7.89 \pm 0.13$), $A({\rm S}) = 6.91 \pm 0.07$ ($7.16 \pm 0.17$), $A({\rm Ar}) = 6.47 \pm 0.08$ ($6.45 \pm 0.07$), and $A({\rm K}) = 5.73 \pm 0.19$ ($5.86 \pm 0.20$), with little evidence for time variations of abundances within the evolution of each flare. Our previous estimates of the Ar and K flare abundances are thus confirmed by this analysis but those for Si and S are reduced. This suggests the flare abundances of Si and Ar are very close to the photospheric abundance or solar proxies, while S is significantly less than photospheric and the K abundance is much higher than photospheric. These estimates differ to some extent from those in which a single enhancement factor applies to elements with first ionization potential less than 10~eV.
We review our recent results on a unified normal mode approach to dynamic tides proposed in Ivanov, Papaloizou $\&$ Chernov (2013) and Chernov, Papaloizou $\&$ Ivanov (2013). Our formalism can be used whenever the tidal interactions are mainly determined by normal modes of a star with identifiable regular spectrum of low frequency modes. We provide in the text basic expressions for tidal energy and angular momentum transfer valid both for periodic and parabolic orbits, and different assumptions about efficiency of normal mode damping due to viscosity and/or non-linear effects and discuss applications to binary stars and close orbiting extrasolar planets.
We use a planetary albedo model to investigate variations in visible wavelength phase curves of exoplanets. The presence of clouds on these exoplanets significantly alters their planetary albedo spectra. We confirm that non-uniform cloud coverage on the dayside of tidally locked exoplanets will manifest as changes to the magnitude and shift of the phase curve. In this work, we first investigate a test case of our model using a Jupiter-like planet, at temperatures consistent to 2.0 AU insolation from a solar type star, to consider the effect of H2O clouds. We then extend our application of the model to the exoplanet Kepler-7b and consider the effect of varying cloud species, sedimentation efficiency, particle size, and cloud altitude. We show that, depending on the observational filter, the largest possible shift of the phase curve maximum will be 2-10 deg for a Jupiter-like planet, and up to 30 deg (0.08 in fractional orbital phase) for hot-Jupiter exoplanets at visible wavelengths as a function of dayside cloud distribution with a uniformly averaged thermal profile. Finally, we tailor our model for comparison with, and confirmation of, the recent optical phase-curve observations of Kepler-7b with the Kepler space telescope. The average planetary albedo can vary between 0.1-0.6 for the 1300 cloud scenarios that were compared to the observations. We observe that smaller particle size and increasing cloud altitude have a strong effect on increasing albedo. In particular, we show that a set of models where Kepler-7b has roughly half of its dayside covered in small-particle clouds high in the atmosphere, made of bright minerals like MgSiO3 and Mg2SiO4, provide the best fits to the observed offset and magnitude of the phase-curve, whereas Fe clouds are found to have too dark to fit the observations.
Analysis of solar radius measurements acquired by the Michelson Doppler Imager on the SOHO spacecraft supports previously reported evidence of solar internal r-mode oscillations in Mt Wilson radius data and in nuclear-decay data acquired at the Lomonosov Moscow State University. The frequencies of these oscillations are compatible with oscillations in a putative inner tachocline that separates a slowly rotating core from the radiative envelope.
The explosion of a supernovae (SN) represents the sudden injection of about 10^51 ergs of thermal and mechanical energy in a small region of space, causing the formation of powerful shock waves that propagate through the interstellar medium at speeds of several thousands of km/s. These waves sweep, compress and heat the interstellar material that they encounter, forming the supernova remnants. Their evolution over thousands of years change forever, irreversibly, not only the physical but also the chemical properties of a vast region of space that can span hundreds of parsecs. This contribution briefly analyzes the impact of these explosions, discussing the relevance of some phenomena usually associated with SNe and their remnants in the light of recent theoretical and observational results.
The aim of this PhD Thesis is to characterize a representative sample of Supergiant X-ray Binaries (SGXBs) formed by 4 sources: XTE J1855-026, a classical SGXB with long-term stable X-ray flux; AX J1841.0-0535 and AX J1845.0-0433, two supergiant fast X-ray transients (SFXTs) with the X-ray emission mostly dominated by flaring; and IGR J00370+6122, something in between these 2 sub-groups. The physical processes that produce these observable differences are still a matter of debate. In this PhD Thesis I performed a study of these 4 different systems to provide new data to constrain the models. This study consists of:(i) the determination of the orbital solution,(ii) a systematic study of the wind behavior along the orbit by the measure of Halpha variations,(iii) a model of stellar atmospheres of the donor star,(iv) establish whether there are X-ray flux variations modulated by the orbital period. The study of the wind shows that Halpha variations are dominated by intrinsic wind processes. The stellar atmospheres study shows that the supergiant stars that harbor these binaries have a higher projected rotational velocity, higher He abundance and higher N/C ratio than that of isolated supergiant stars. The results of this study show that the eccentricity of the binary does not have a simple correlation with the differences in the X-ray flux of the different sub-groups. Therefore, the idea of a more complex scenario is consolidated. The discovery of a new type of system, IGR J00370+6122, which properties do not fit in any of the established sub-groups, reinforces the idea of a continuum in the observed properties more than a strict classification of systems. Furthermore, I have developed a pipeline to reduce spectra of the FRODOSpec spectrograph at the Liverpool Telescope optimized for the reduction of the red spectra of the obscured supergiants that we find in these SGXBs.
As a solution to the well-known problem that the shock wave potentially responsible for the explosion of a supernova actually tends to stall, we propose a new energy source arising from our model for dark matter. Our earlier model proposed that dark matter should consist of cm-large white dwarf-like objects kept together by a skin separating two different sorts of vacua. These dark matter balls or pearls will collect in the middle of any star throughout its lifetime. At some stage during the development of a supernova the balls will begin to take in neutrons and then other surrounding material. By passing into a ball nucleons fall through a potential of order 10 MeV, causing a severe production of heat - of order 10 foe for a solar mass of material eaten by the balls. The temperature in the iron core will thereby be raised, splitting up the iron into smaller nuclei. This provides a mechanism for reviving the shock wave when it arrives and making the supernova explosion really occur. The onset of the heating due to the dark matter balls would at first stop the collapse of the supernova progenitor. This opens up the possibility of there being {\em two} collapses giving two neutrino outbursts, as apparently seen in the supernova SN1987A - one in Mont Blanc, and one 4 hours 43 minutes later in both IMB and Kamiokande.
We obtain conditions for the existence of an attractor in the system of equations describing a tachyon warm inflationary model with bulk viscosity taken into account. When these conditions are met the evolution approaches slow-roll regime. We present the primordial power spectrum for the tachyon field by considering a dissipation coefficient depending on the scalar field and temperature.
We derive an analytical expression for extracting the gravitational waveforms at null infinity using the Weyl scalar $\psi_4$ measured at a finite radius. Our expression is based on a series solution in orders of 1/r to the equations for gravitational perturbations about a spinning black hole. We compute this expression to order $1/r^2$ and include the spin parameter $a$ of the Kerr background. We test the accuracy of this extraction procedure by measuring the waveform for a merging black-hole binary at ten different extraction radii (in the range r/M=75-190) and for three different resolutions in the convergence regime. We find that the extraction formula provides a set of values for the radiated energy and momenta that at finite extraction radii converges towards the expected values with increasing resolution, which is not the case for the `raw' waveform at finite radius. We also examine the phase and amplitude errors in the waveform as a function of observer location and again observe the benefits of using our extraction formula. The leading corrections to the phase are ${\cal O}(1/r)$ and to the amplitude are ${\cal O}(1/r^2)$. This method provides a simple and practical way of estimating the waveform at infinity, and may be especially useful for scenarios such as well separated binaries, where the radiation zone is far from the sources, that would otherwise require extended simulation grids in order to extrapolate the `raw' waveform to infinity. Thus this method saves important computational resources and provides an estimate of errors.
We derive constraints facing models of axion inflation based on decay constant alignment from a string-theoretic and quantum gravitational perspective. In particular, we investigate the prospects for alignment and `anti-alignment' of $C_4$ axion decay constants in type IIB string theory, deriving a strict no-go result in the latter case. We discuss the relationship of axion decay constants to the weak gravity conjecture and demonstrate agreement between our string-theoretic constraints and those coming from the `generalized' weak gravity conjecture. Finally, we consider a particular model of decay constant alignment in which the potential of $C_4$ axions in type IIB compactifications on a Calabi-Yau three-fold is dominated by contributions from $D7$-branes, pointing out that this model evades some of the challenges derived earlier in our paper but is highly constrained by other geometric considerations.
Voyager 1 data from beyond the heliopause provide the first direct measurements of the interstellar cosmic ray spectra below 1 GeVnuc. In this paper we combine these Voyager measurements of H and He nuclei from 3-600 MeVnuc with higher energy measurements at 1 AU from the BESS and PAMELA experiments up to 100 GeVnuc. Using a Weighted Leaky Box Model for propagation in the galaxy, we obtain an excellent fit to these new Voyager observations and the much higher energy spectra up to 100 GeVnuc by using source spectra which are P-2.28, with the exponent independent of rigidity from low to high rigidities; along with a rigidity dependence of the diffusion path length which is P-0.5 at rigidities 1.00 GV, and possibly changing to P1.0 at lower rigidities.
The magnetic fields of the Ice Giant Planets Uranus and Neptune (U/N) are unique in the solar system. Based on a substantial database measured on Earth for representative planetary fluids at representative dynamic pressures up to 200 GPa (2 Mbar) and a few 1000 K, the complex magnetic fields of U/N are (i) probably made primarily by degenerate metallic fluid H (MFH) at or near the crossover from the H-He envelopes to Ice cores at ~100 GPa (Mbar) pressures and normalized radii of ~90% of the radii of U/N; (ii) because those magnetic fields are made relatively close to the surfaces of U/N, non-dipolar fields can be expected; (iii) the Ice cores are most probably a heterogeneous fluid mixture of H, N, O, C, Fe/Ni and silicate-oxides and their mutual reaction products at high pressures and temperatures; (iv) the shapes of the magnetic fields are probably caused by weak coupling between rotational motions of U/N and convective motions of conducting fluids in dynamos that make those magnetic fields. Ironically, there is probably little nebular Ice in the Ice Giant Planets.
Broadband suppression of quantum noise below the Standard Quantum Limit (SQL) becomes a top-priority problem for the future generation of large-scale terrestrial detectors of gravitational waves, as the interferometers of the Advanced LIGO project, predesigned to be quantum-noise-limited in the almost entire detection band, are phased in. To this end, among various proposed methods of quantum noise suppression or signal amplification, the most elaborate approach implies a so-called *xylophone* configuration of two Michelson interferometers, each optimised for its own frequency band, with a combined broadband sensitivity well below the SQL. Albeit ingenious, it is a rather costly solution. We demonstrate that changing the optical scheme to a Sagnac interferometer with weak detuned signal recycling and frequency dependent input squeezing can do almost as good a job, as the xylophone for significantly lower spend. We also show that the Sagnac interferometer is more robust to optical loss in filter cavity, used for frequency dependent squeezed vacuum injection, than an analogous Michelson interferometer, thereby reducing building cost even more.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
The presence of hot gaseous coronae around present-day massive spiral galaxies is a fundamental prediction of galaxy formation models. However, our observational knowledge remains scarce, since to date only four gaseous coronae were detected around spirals with massive stellar bodies ($\gtrsim2\times10^{11} \ \rm{M_{\odot}}$). To explore the hot coronae around lower mass spiral galaxies, we utilized Chandra X-ray observations of a sample of eight normal spiral galaxies with stellar masses of $(0.7-2.0)\times10^{11} \ \rm{M_{\odot}}$. Although statistically significant diffuse X-ray emission is not detected beyond the optical radii ($\sim20$ kpc) of the galaxies, we derive $3\sigma$ limits on the characteristics of the coronae. These limits, complemented with previous detections of NGC 1961 and NGC 6753, are used to probe the Illustris Simulation. The observed $3\sigma$ upper limits on the X-ray luminosities and gas masses exceed or are at the upper end of the model predictions. For NGC 1961 and NGC 6753 the observed gas temperatures, metal abundances, and electron density profiles broadly agree with those predicted by Illustris. These results hint that the physics modules of Illustris are broadly consistent with the observed properties of hot coronae around spiral galaxies. However, a shortcoming of Illustris is that massive black holes, mostly residing in giant ellipticals, give rise to powerful radio-mode AGN feedback, which results in under luminous coronae for ellipticals.
We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.
We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.
In this paper, we present the first results from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich AMR calculations of high redshift galaxy formation performed on the Blue Waters supercomputer. These simulations contain hundreds of well-resolved galaxies at $z \sim 25-8$, and make several novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity function (M$_{1600} \leq -17$), but at lower luminosities is essentially flat rather than rising steeply, as has been inferred by Schechter function fits to high-z observations. This flattening of the luminosity function is due to two factors: (i) the strong dependence of the stellar fraction on halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower stellar fractions and thus lower luminosities at a given halo virial mass; and (ii) the fact that halos with virial masses below $\simeq 2 \times 10^8$ M$_\odot$ do not universally contain stars, with the fraction of halos containing stars dropping to zero at $\simeq 7 \times 10^6$ M$_\odot$. Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep space-based surveys, particularly if lensed regions are chosen for observation.
Molecular outflows driven by protostellar cluster members likely impact their surroundings and contribute to turbulence, affecting subsequent star formation. The very young Serpens South cluster consists of a particularly high density and fraction of protostars, yielding a relevant case study for protostellar outflows and their impact on the cluster environment. We combined CO $J=1-0$ observations of this region using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Institut de Radioastronomie Millim\'{e}trique (IRAM) 30 m single dish telescope. The combined map allows us to probe CO outflows within the central, most active region at size scales of 0.01 pc to 0.8 pc. We account for effects of line opacity and excitation temperature variations by incorporating $^{12}$CO and $^{13}$CO data for the $J=1-0$ and $J=3-2$ transitions (using Atacama Pathfinder Experiment and Caltech Submillimeter Observatory observations for the higher CO transitions), and we calculate mass, momentum, and energy of the molecular outflows in this region. The outflow mass loss rate, force, and luminosity, compared with diagnostics of turbulence and gravity, suggest that outflows drive a sufficient amount of energy to sustain turbulence, but not enough energy to substantially counter the gravitational potential energy and disrupt the clump. Further, we compare Serpens South with the slightly more evolved cluster NGC 1333, and we propose an empirical scenario for outflow-cluster interaction at different evolutionary stages.
This is the first of a series of papers where we compare the expected performance of two of the largest stage IV next-generation surveys in the optical and infrared (LSST and Euclid), with a particular focus on cluster surveys. In this first paper, we introduce the mock catalogues we have utilized in this work, an N-body simulation+semi-analytical cone with a posterior modification with PhotReal, a technique which modifies the original photometry to make it more realistic by using an empirical library of spectral templates. We have confirmed the reliability of the mock catalogue by comparing the obtained color-magnitude relation, the luminosity and mass function and the angular correlation function with those of real data. We also analyze the behavior of the expected photometric redshifts for each different survey, in terms of photometric redshift resolution, photometric redshift bias and fraction of outliers. In addition, we discuss the benefits of using the BPZ \emph{odds} photometric redshift quality parameter to select the best quality data of the sample. We find that very deep near infrared surveys such as Euclid will provide very good performance ($\Delta z/(1+z) \sim 0.025-0.053$) down to H$\sim$24 AB mag and up to redshift $\sim 3$ depending on the optical observations available from the ground whereas extremely deep optical surveys such as LSST will obtain an overall lower photometric redshift resolution ($\Delta z/(1+z) \sim 0.045$) down to $i\sim27.5$ AB mag, being substantially improved ($\Delta z/(1+z) \sim 0.035$) if we restrict the sample down to i$\sim$24 AB mag. We highlight the fact that those numbers can be improved substantially by selecting a subsample of galaxies with the best quality photometric redshifts. We finally discuss the impact that these surveys will have for the community in terms of photometric redshift legacy once the data is available. (Abridged)
The coolest dwarf stars targeted by the Kepler Mission constitute a relatively small but scientifically valuable subset of the Kepler target stars, and provide a high-fidelity and nearby sample of transiting planetary systems. Using archival Kepler data spanning the entire primary mission we perform a uniform analysis to extract, confirm and characterize the transit signals discovered by the Kepler pipeline toward M-type dwarf stars. We recover all but two of the signals reported in a recent listing from the Exoplanet Archive resulting in 165 planet candidates associated with a sample of 106 low-mass stars. We fitted the observed light curves to transit models using Markov Chain Monte Carlo and we have made the posterior samples publicly available to facilitate further studies. We fitted empirical transit times to individual transit signals with significantly non-linear ephemerides for accurate recovery of transit parameters and measuring precise transit timing variations. We also provide the physical parameters for the stellar sample, including new measurements of stellar rotation, allowing the conversion of transit parameters into planet radii and orbital parameters.
One of the most compelling tasks of modern cosmology is to constrain the expansion history of the Universe, since this measurement can give insights on the nature of dark energy and help to estimate cosmological parameters. In this letter are presented two new measurements of the Hubble parameter H(z) obtained with the cosmic chronometer method up to $z\sim2$. Taking advantage of near-infrared spectroscopy of the few very massive and passive galaxies observed at $z>1.4$ available in literature, the differential evolution of this population is estimated and calibrated with different stellar population synthesis models to constrain H(z), including in the final error budget all possible sources of systematic uncertainties (star formation history, stellar metallicity, model dependencies). This analysis is able to extend significantly the redshift range coverage with respect to present-day constraints, crossing for the first time the limit at $z\sim1.75$. The new H(z) data are used to estimate the gain in accuracy on cosmological parameters with respect to previous measurements in two cosmological models, finding a small but detectable improvement ($\sim$5 %) in particular on $\Omega_{M}$ and $w_{0}$. Finally, a simulation of a Euclid-like survey has been performed to forecast the expected improvement with future data. The provided constraints have been obtained just with the cosmic chronometers approach, without any additional data, and the results show the high potentiality of this method to constrain the expansion history of the Universe at these redshifts.
We study the angular momentum of galaxies in the Illustris cosmological simulation, which captures gravitational and gas dynamics within galaxies, as well as feedback from stars and black holes. We find that the angular momentum of the simulated galaxies matches observations well, and in particular two distinct relations are found for late-type versus early-type galaxies. The relation for late-type galaxies corresponds to the value expected from full conservation of the specific angular momentum generated by cosmological tidal torques. The relation for early-type galaxies corresponds to retention of only ~30% of that, but we find that those early-type galaxies with low angular momentum at z=0 nevertheless reside at high redshift on the late-type relation. To gain further insight, we explore the scaling relations in simulations where the galaxy formation physics is modified with respect to the fiducial model. We find that galactic winds with high mass-loading factors are essential for obtaining the high angular momentum relation typical for late-type galaxies, while AGN feedback largely operates in the opposite direction. Hence, feedback controls the angular momentum of galaxies, and appears to be instrumental for establishing the Hubble sequence.
We investigate signatures of population III (PopIII) stars in the metal-enriched environment of GRBs originating from population II-I (PopII/I) stars by using abundance ratios derived from numerical simulations that follow stellar evolution and chemical enrichment. We find that at $z>10$ more than $10%$ of PopII/I GRBs explode in a medium previously enriched by PopIII stars (we refer to them as GRBII$\rightarrow$III). Although the formation of GRBII$\rightarrow$III is more frequent than that of pristine PopIII GRBs (GRBIIIs), we find that the expected GRBII$\rightarrow$III observed rate is comparable to that of GRBIIIs, due to the usually larger luminosities of these latter. GRBII$\rightarrow$III events take place preferentially in small proto-galaxies with stellar masses $\rm M_\star \sim 10^{4.5} - 10^7\,\rm M_\odot$, star formation rates $\rm SFR \sim 10^{-3}-10^{-1}\,\rm M_\odot/yr$ and metallicities $Z \sim 10^{-4}-10^{-2}\,\rm Z_\odot$. On the other hand, galaxies with $Z < 10^{-2.8}\,\rm Z_\odot$ are dominated by metal enrichment from PopIII stars and should preferentially host GRBII$\rightarrow$III. Hence, measured GRB metal content below this limit could represent a strong evidence of enrichment by pristine stellar populations. We discuss how to discriminate PopIII metal enrichment on the basis of various abundance ratios observable in the spectra of GRBs' afterglows. By employing such analysis, we conclude that the currently known candidates at redshift $z\simeq 6$ -- i.e. GRB 050904 \cite[][]{2006Natur.440..184K} and GRB 130606A \cite[][]{2013arXiv1312.5631C} -- are likely not originated in environments pre-enriched by PopIII stars.
We present new measurements of the evolution of the X-ray luminosity functions (XLFs) of unabsorbed and absorbed Active Galactic Nuclei (AGNs) out to z~5. We construct samples containing 2957 sources detected at hard (2-7 keV) X-ray energies and 4351 sources detected at soft (0.5-2 keV) energies from a compilation of Chandra surveys supplemented by wide-area surveys from ASCA and ROSAT. We consider the hard and soft X-ray samples separately and find that the XLF based on either (initially neglecting absorption effects) is best described by a new flexible model parametrization where the break luminosity, normalization and faint-end slope all evolve with redshift. We then incorporate absorption effects, separately modeling the evolution of the XLFs of unabsorbed ($20<\log N_H<22$) and absorbed ($22<\log N_H<24$) AGNs, seeking a model that can reconcile both the hard- and soft-band samples. We find that the absorbed AGN XLF generally has a lower break luminosity, a higher normalization, and a steeper faint-end slope than the unabsorbed AGN XLF. Hence, absorbed AGNs tend to dominate at low luminosities, with the absorbed fraction falling rapidly as luminosity increases. The XLFs of both populations undergo strong luminosity evolution which shifts the transition in the absorbed fraction to higher luminosities at higher redshifts. However, differences in the evolution of the two XLFs lead to a comparatively complex evolution in the shape of the total XLF of AGNs. Our work indicates that the evolution of AGNs may be driven by a combination of changes in the distributions of black hole mass and/or Eddington ratio, as well as the life cycles of unabsorbed and absorbed growth phases.
Ram pressure stripping can remove hot and cold gas from galaxies in the intracluster medium (ICM), as shown by observations of X-ray and HI galaxy wakes in nearby clusters of galaxies. However, ram pressure stripping, including pre-processing in group environments, does not remove all the hot coronal gas from cluster galaxies. Recent high-resolution Chandra observations have shown that $\sim 1 - 4$ kpc extended, hot galactic coronae are ubiquitous in group and cluster galaxies. To better understand this result, we simulate ram pressure stripping of a cosmologically motivated population of galaxies in isolated group and cluster environments. The galaxies and the host group and cluster are composed of collisionless dark matter and hot gas initially in hydrostatic equilibrium with the galaxy and host potentials. We show that the rate at which gas is lost depends on the galactic and host halo mass. Using synthetic X-ray observations, we evaluate the detectability of stripped galactic coronae in real observations by stacking images on the known galaxy centers. We find that coronal emission should be detected within $\sim 10$ arcsec, or $\sim 5$ kpc up to $\sim 2.3$ Gyr in the lowest (0.1 - 1.2 keV) energy band. Thus the presence of observed coronae in cluster galaxies significantly smaller than the hot X-ray halos of field galaxies indicates that at least some gas removal occurs within cluster environments for recently accreted galaxies. Finally, we evaluate the possibility that existing and future X-ray cluster catalogs can be used in combination with optical galaxy positions to detect galactic coronal emission via stacking analysis. We briefly discuss the effects of additional physical processes on coronal survival, and will address them in detail in future papers in this series.
Standard cosmological models based on general relativity (GR) with dark energy predict that the Universe underwent a transition from decelerating to accelerating expansion at a moderate redshift $z_{acc} \sim 0.7$. Clearly, it is of great interest to directly measure this transition in a model-independent way, without the assumption that GR is the correct theory of gravity. We explore to what extent supernova (SN) luminosity distance measurements provide evidence for such a transition: we show that, contrary to intuition, the well-known "turnover" in the SN distance residuals $\Delta\mu$ relative to an empty (Milne) model does not give firm evidence for such a transition within the redshift range spanned by SN data. The observed turnover in that diagram is predominantly due to the negative curvature in the Milne model, {\em not} the deceleration predicted by $\Lambda$CDM and relatives. We show that there are several advantages in plotting distance residuals against a flat, non-accelerating model $(w = -1/3)$, and also remapping the $z-$axis to $u = \ln(1+z)$; we outline a number of useful and intuitive properties of this presentation. We conclude that there are significant complementarities between SNe and baryon acoustic oscillations (BAOs): SNe offer high precision at low redshifts and give good constraints on the net {\em amount} of acceleration since $z \sim 0.7$, but are weak at constraining $z_{acc}$; while radial BAO measurements are probably superior for placing direct constraints on $z_{acc}$.
MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a 6-year
SDSS-IV survey that will obtain resolved spectroscopy from 3600 $\AA$ to 10300
$\AA$ for a representative sample of over 10,000 nearby galaxies. In this
paper, we derive spatially resolved stellar population properties and radial
gradients by performing full spectral fitting of observed galaxy spectra from
P-MaNGA, a prototype of the MaNGA instrument. These data include spectra for
eighteen galaxies, covering a large range of morphological type. We derive age,
metallicity, dust and stellar mass maps, and their radial gradients, using high
spectral-resolution stellar population models, and assess the impact of varying
the stellar library input to the models. We introduce a method to determine
dust extinction which is able to give smooth stellar mass maps even in cases of
high and spatially non-uniform dust attenuation.
With the spectral fitting we produce detailed maps of stellar population
properties which allow us to identify galactic features among this diverse
sample such as spiral structure, smooth radial profiles with little azimuthal
structure in spheroidal galaxies, and spatially distinct galaxy sub-components.
In agreement with the literature, we find the gradients for galaxies identified
as early-type to be on average flat in age, and negative (- 0.15 dex / R$_e$ )
in metallicity, whereas the gradients for late-type galaxies are on average
negative in age (- 0.39 dex / R$_e$ ) and flat in metallicity. We demonstrate
how different levels of data quality change the precision with which radial
gradients can be measured. We show how this analysis, extended to the large
numbers of MaNGA galaxies, will have the potential to shed light on galaxy
structure and evolution.
This article details the computation of the two-point correlators of the convergence, $E$- and $B$-modes of the cosmic shear induced by the weak-lensing by large scale structure assuming that the background spacetime is spatially homogeneous and anisotropic. After detailing the perturbation equations and the general theory of weak-lensing in an anisotropic universe, it develops a weak shear approximation scheme in which one can compute analytically the evolution of the Jacobi matrix. It allows one to compute the angular power spectrum of the $E$- and $B$-modes. In the linear regime, the existence of $B$-modes is a direct tracer of a late time anisotropy and their angular power spectrum scales as the square of the shear. It is then demonstrated that there must also exist off-diagonal correlations between the $E$-modes, $B$-modes and convergence that are linear in the geometrical shear and allow one to reconstruct the eigendirections of expansion. These spectra can be measured in future large scale surveys, such as Euclid and SKA, and offer a new tool to test the isotropy of the expansion of the universe at low redshift.
We compute the angular power spectrum of the $B$-modes of the weak-lensing shear in a spatially anisotropic spacetime. We find that there must also exist off-diagonal correlations between the $E$-modes, $B$-modes, and convergence that allow one to reconstruct the eigendirections of expansion. Focusing on future surveys such as Euclid and SKA, we show that observations can constrain the geometrical shear in units of the Hubble rate at the percent level, or even better, offering a new and powerful method to probe our cosmological model.
We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).
The structure of Gamma-Ray Burst (GRB) jets impacts on their prompt and afterglow emission properties. Insights into the still unknown structure of GRBs can be achieved by studying how different structures impact on the luminosity function (LF): i) we show that low ($10^{46} < L_{\rm iso} < 10^{48}$ erg/s) and high (i.e. with $L_{\rm iso} > 10^{50}$ erg/s) luminosity GRBs can be described by a unique LF; ii) we find that a uniform jet (seen on- and off-axis) as well as a very steep structured jet (i.e. $\epsilon(\theta) \propto \theta^{-s}$ with $s > 4$) can reproduce the current LF data; iii) taking into account the emission from the whole jet (i.e. including contributions from mildly relativistic, off-axis jet elements) we find that $E_{\rm iso}(\theta_{\rm v})$ (we dub this quantity "apparent structure") can be very different from the intrinsic structure $\epsilon(\theta)$: in particular, a jet with a Gaussian intrinsic structure has an apparent structure which is more similar to a power law. This opens a new viewpoint on the quasi-universal structured jet hypothesis.
We present spectroscopy of 28 SNR candidates as well as one H II region in M81, and two SNR candidates in M82. Twenty six out of the M81 candidates turn out to be genuine SNRs, and two in M82 may be shocked condensations in the galactic outflow or SNRs. The distribution of [N II]/H{\alpha} ratios of M81 SNRs is bimodal. M81 SNRs are divided into two groups in the spectral line ratio diagrams: an [O III]-strong group and an [O III]-weak group. The latter have larger sizes, and may have faster shock velocity. [N II]/H{\alpha} ratios of the SNRs show a strong correlation with [S II]/H{\alpha} ratios. They show a clear radial gradient in [N II]/H{\alpha} and [S II]/H{\alpha} ratios: dLog ([N II]/H{\alpha})/dLog R = -0.018 {\pm} 0.008 dex/kpc and dLog ([S II]/H{\alpha})/dLog R = -0.016 {\pm} 0.008 dex/kpc where R is a deprojected galactocentric distance. We estimate the nitrogen and oxygen abundance of the SNRs from the comparison with shock-ionization models. We obtain a value for the nitrogen radial gradient, dLog(N/H)/dLogR = -0.023 {\pm} 0.009 dex/kpc, and little evidence for the gradient in oxygen. This nitrogen abundance shows a few times flatter gradient than those of the planetary nebulae and H II regions. We find that five SNRs are matched with X-ray sources. Their X-ray hardness colors are consistent with thermal SNRs.
Aims: In this work, we discuss a way to combine High Dispersion Spectroscopy
and High Contrast Imaging (HDS+HCI). For a planet located at a resolvable
angular distance from its host star, the starlight can be reduced up to several
orders of magnitude using adaptive optics and/or coronography. In addition, the
remaining starlight can be filtered out using high-dispersion spectroscopy,
utilizing the significantly different (or Doppler shifted) high-dispersion
spectra of the planet and star. In this way, HDS+HCI can in principle reach
contrast limits of ~1e-5 x 1e-5, although in practice this will be limited by
photon noise and/or sky-background.
Methods: We present simulations of HDS+HCI observations with the E-ELT, both
probing thermal emission from a planet at infrared wavelengths, and starlight
reflected off a planet atmosphere at optical wavelengths. For the infrared
simulations we use the baseline parameters of the E-ELT and METIS instrument,
with the latter combining extreme adaptive optics with an R=100,000 IFS. We
include realistic models of the adaptive optics performance and atmospheric
transmission and emission. For the optical simulation we also assume R=100,000
IFS with adaptive optics capabilities at the E-ELT.
Results: One night of HDS+HCI observations with the E-ELT at 4.8 um (d_lambda
= 0.07 um) can detect a planet orbiting alpha Cen A with a radius of R=1.5
R_earth and a twin-Earth thermal spectrum of T_eq=300 K at a signal-to-noise
(S/N) of 5. In the optical, with a Strehl ratio performance of 0.3, reflected
light from an Earth-size planet in the habitable zone of Proxima Centauri can
be detected at a S/N of 10 in the same time frame. Recently, first HDS+HCI
observations have shown the potential of this technique by determining the
spin-rotation of the young massive exoplanet beta Pictoris b. [abridged]
In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the AGN continuum, and measurements of the broad H-beta line widths in mean and root-mean square (rms) spectra. For the most highly variable AGNs we also measured broad H-beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H-beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H-beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593 the broad H-beta velocity shifted by ~250 km/s over a one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise.
We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $\alpha$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$\alpha$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$\alpha$ radiative transfer and escape. Ly$\alpha$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$\alpha$ escape. We find that Ly$\alpha$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$\alpha$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model.
Primordial Black Holes (PBHs) are of interest in many cosmological contexts. PBHs lighter than about 1012 kg are predicted to be directly detectable by their Hawking radiation. This radiation should produce both a diffuse extragalactic gamma-ray background from the cosmologically-averaged distribution of PBHs and gamma-ray burst signals from individual light black holes. The Fermi, Milagro, Veritas, HESS and HAWC observatories, in combination with new burst recognition methodologies, offer the greatest sensitivity for the detection of such black holes or placing limits on their existence.
ASTRO-H, the sixth Japanese X-ray observatory, which is scheduled to be launched by the end of Japanese fiscal year 2015 has a capability to observe the prompt emission from Gamma-ray Bursts (GRBs) utilizing BGO active shields for the soft gamma-ray detector (SGD). The effective area of the SGD shield detectors is very large and its data acquisition system is optimized for short transients such as short GRBs. Thus, we expect to perform more detailed time-resolved spectral analysis with a combination of ASTRO-H and Fermi LAT/GBM to investigate the gamma-ray emission mechanism of short GRBs. In addition, the environment of the GRB progenitor should be a remarkable objective from the point of view of the chemical evolution of high-z universe. If we can maneuver the spacecraft to the GRBs, we can perform a high-resolution spectroscopy of the X-ray afterglow of GRBs utilizing the onboard micro calorimeter and X-ray CCD camera.
The NASA/IPAC Extragalactic Database (NED) has deployed a new rule-based cross-matching algorithm called Match Expert (MatchEx), capable of cross-matching very large catalogs (VLCs) with >10 million objects. MatchEx goes beyond traditional position-based cross-matching algorithms by using other available data together with expert logic to determine which candidate match is the best. Furthermore, the local background density of sources is used to determine and minimize the false-positive match rate and to estimate match completeness. The logical outcome and statistical probability of each match decision is stored in the database, and may be used to tune the algorithm and adjust match parameter thresholds. For our first production run, we cross-matched the GALEX All Sky Survey Catalog (GASC), containing nearly 40 million NUV-detected sources, against a directory of 180 million objects in NED. Candidate matches were identified for each GASC source within a 7.5 arcsecond radius. These candidates were filtered on position-based matching probability, and on other criteria including object type and object name. We estimate a match completeness of 97.6% and a match accuracy of 99.75%. MatchEx is being used to cross-match over 2 billion catalog sources to NED, including the Spitzer Source List, the 2MASS Point-Source Catalog, AllWISE, and SDSS DR 10. It will also speed up routine cross-matching of sources as part of the NED literature pipeline.
The baseline energy-resolution performance for the current generation of large-mass, low-temperature calorimeters is $>2$ orders of magnitude worse than theoretical predictions. A detailed study of several calorimetric detectors suggests that a mismatch between the sensor and signal bandwidths is the primary reason for suppressed sensitivity. With this understanding, we propose a detector design in which a thin-film Au pad is directly deposited onto a massive absorber that is then thermally linked to a separately fabricated TES chip via an Au wirebond, providing large electron-phonon coupling (i.e. high signal bandwidth), ease of fabrication, and cosmogenic background suppression. Interestingly, this design strategy is fully compatible with the use of hygroscopic crystals (NaI) as absorbers. An 80-mm diameter Si light detector based upon these design principles, with potential use in both dark matter and neutrinoless double-beta decay, has an estimated baseline energy resolution of 0.35eV, 20x better than currently achievable. A 1.75 kg ZnMoO$_{4}$ large-mass calorimeter would have a 3.5eV baseline resolution, 1000x better than currently achieved with NTDs with an estimated position dependence $\frac{\Delta E}{E}$ of 6$x$10$^{-4}$, near or below the variations found in absorber thermalization in ZnMoO$_{4}$ and TeO$_{2}$. Such minimal position dependence is made possible by forcing the sensor bandwidth to be much smaller than the signal bandwidth. Further, intrinsic event timing resolution is estimated to be $\sim$170 $\mu$s for 3 MeV recoils in the phonon detector, satisfying the event-rate requirements of next-generation neutrinoless double-beta decay experiments. Quiescent bias power for both of these designs is found to be significantly larger than parasitic power loads achieved in the SPICA/SAFARI infrared bolometers.
We discuss two multiple star systems that host known exoplanets: HD 2638 and 30 Ari B. Adaptive optics imagery revealed an additional stellar companion to both stars. We collected multi-epoch images of the systems with Robo-AO and the PALM-3000 adaptive optics systems at Palomar Observatory and provide relative photometry and astrometry. The astrometry indicates that the companions share common proper motion with their respective primaries. Both of the new companions have projected separations less than 30 AU from the exoplanet host star. Using the projected separations to compute orbital periods of the new stellar companions, HD 2638 has a period of 130 yrs and 30 Ari B has a period of 80 years. Previous studies have shown that the true period is most likely within a factor of three of these estimated values. The additional component to the 30 Ari makes it the second confirmed quadruple system known to host an exoplanet. HD 2638 hosts a hot Jupiter and the discovery of a new companion strengthens the connection between hot Jupiters and binary stars. We place the systems on a color-magnitude diagram and derive masses for the companions which turn out to be roughly 0.5 solar mass stars.
A fraction of the mid-$J$ ($J$= 14--13 to $J$= 24--23) CO emission detected by the \textit{Herschel}/PACS observations of embedded young stellar objects (YSOs) has been attributed to the UV-heated outflow cavity walls. We have applied our newly developed self-consistent models of Photon-Dominated Region (PDR) and non-local thermal equilibrium line Radiative transfer In general Grid (RIG) to the \textit{Herschel} FIR observations of 27 low mass YSOs and one intermediate mass YSO, NGC7129-FIRS2. When the contribution of the hot component (traced by transitions of $J> 24$) is removed, the rotational temperature of the warm component is nearly constant with $\sim250$ K. This can be reproduced by the outflow cavity wall ($n \geq 10^6\, \mathrm{cm}^{-3}$, $\log G_{0}/n \geq-4.5$, $\mathrm{log} G_0\ge 3$, $T_{\rm gas} \ge 300 $K, and X(CO)$ \ge 10^{-5}$) heated by a UV radiation field with a black body temperature of 15,000 K or 10,000 K. However, a shock model combined with an internal PDR will be required to determine the quantitative contribution of a PDR relative to a shock to the mid-$J$ CO emission.
We present medium resolution near-infrared (NIR) spectra, covering 1.1 to 3.4 microns, of the normal Type Ia supernova (SN Ia) SN 2014J in M82 obtained with the FLITECAM instrument aboard SOFIA approximately 17-25 days after maximum B light. Our 2.8-3.4 micron spectra may be the first ~3 micron spectra of a SN Ia ever published. The spectra spanning the 1.5-2.7 micron range are characterized by a strong emission feature at ~1.77 microns with a full width at half maximum of ~11,000-13,000 km/s. We compare the observed FLITECAM spectra to the recent non-LTE delayed detonation models of Dessart et al. (2014) and find that the models agree with the spectra remarkably well in the 1.5-2.7 micron wavelength range. Based on this comparison we identify the ~1.77 micron emission peak as a blend of permitted lines of Co II. Other features seen in the 2.0 - 2.5 micron spectra are also identified as emission from permitted transitions of Co II. However, the models are not as successful at reproducing the spectra in the 1.1 - 1.4 micron range or between 2.8 and 3.4 microns. These observations demonstrate the promise of SOFIA by allowing access to wavelength regions inaccessible from the ground, and serve to draw attention to the usefulness of the regions between the standard ground-based NIR passbands for constraining SN models.
The $\Lambda$ cold dark matter ($\Lambda\textrm{CDM}$) model is currently known as the simplest cosmology model that best describes observations with minimal number of parameters. Here we introduce a cosmology model that is preferred over the conventional $\Lambda\textrm{CDM}$ one by constructing dark energy as the sum of the cosmological constant $\Lambda$ and the additional fluid that is designed to have an extremely short transient spike in energy density during the radiation-matter equality era and the early scaling behavior with radiation and matter densities. The density parameter of the additional fluid is defined as a Gaussian function plus a constant in logarithmic scale-factor space. Searching for the best-fit cosmological parameters in the presence of such a dark energy spike gives a far smaller chi-square value by about five times the number of additional parameters introduced and narrower constraints on matter density and Hubble constant compared with the best-fit $\Lambda\textrm{CDM}$ model. The significant improvement in reducing chi-square mainly comes from the better fitting of Planck temperature power spectrum around the third ($\ell \approx 800$) and sixth ($\ell \approx 1800$) acoustic peaks. The likelihood ratio test and the Akaike information criterion suggest that the model of dark energy spike is strongly favored by the current cosmological observations over the conventional $\Lambda\textrm{CDM}$ model. However, based on the Bayesian information criterion which penalizes models with more parameters, the strong evidence supporting the presence of dark energy spike disappears. Our result emphasizes that the alternative cosmological parameter estimation with even better fitting of the same observational data is allowed in the Einstein's gravity.
The spectrum of the supernova relic neutrino (SRN) background from past stellar collapses including black hole formation (failed supernovae) is calculated. The redshift dependence of the black hole formation rate is considered on the basis of the metallicity evolution of galaxies. Assuming the mass and metallicity ranges of failed supernova progenitors, their contribution to SRNs is quantitatively estimated for the first time. Using this model, the dependences of SRNs on the cosmic star formation rate density, shock revival time and equation of state are investigated. The shock revival time is introduced as a parameter that should depend on the still unknown explosion mechanism of core collapse supernovae. The dependence on equation of state is considered for failed supernovae, whose collapse dynamics and neutrino emission are certainly affected. It is found that the low-energy spectrum of SRNs is mainly determined by the cosmic star formation rate density. These low-energy events will be observed in the Super-Kamiokande experiment with gadolinium-loaded water.
Optical interferometry is a powerful tool to investigate the close environment of AGB stars. With a spatial resolution of a few milli-arcseconds, it is even possible to image directly the surface of angularly large objects. This is of special interest forMira stars and red supergiants for which the dust-wind is initiated from or very close to the photosphere by an interplay between pulsation and convection. Based on two-epoch interferometric observations of the Mira star X Hya, we present how the variation of the angular size with wavelength challenges pulsation models and how reconstructed images can reveal the evolution of the object shape and of its asymmetric structures.
We report MAXI and Swift observations of short-term spectral softenings of the galactic black-hole X-ray binary Swift J1753.5-0127 in the low/hard state. These softening events are characterized by a simultaneous increase of soft X-rays (2-4 keV) and a decrease of hard X-rays (15-50 keV) lasting for a few tens of days. The X-ray energy spectra during the softening periods can be reproduced with a model consisting of a multi-color disk blackbody and its Comptonized component. The fraction of the Comptonized component decreased from 0.30 to 0.15 when the spectrum became softer; meanwhile the inner disk temperature (Tin) increased from 0.2 to 0.45 keV. These results imply that the softening events are triggered by a short-term increase of the mass accretion rate. During the observed spectral softening events, the disk flux (F) and Tin did not obey the relation: F is proportional to Tin^4, suggesting that the inner disk radius does not reach the innermost stable circular orbit.
We present HI line profiles for various models of circumstellar shells around
red giants. In the calculations we take into account the effect of the
background at 21 cm, and show that in some circumstances it may have an
important effect on the shape and intensity of the observed line profiles. We
show that self-absorption should also be considered depending on the mass loss
rate and the temperature reached by circumstellar gas.
HI emission from circumstellar shells has been mostly reported from stars
with mass loss rates around 10$^{-7}$ solar masses per year. We discuss the
possible reasons for the non detection of many sources with larger mass loss
rates that are hallmarks of the end of the AGB phase. Although radiative
transfer effects may weaken the line emission, they cannot alone account for
this effect. Therefore, it seems likely that molecular hydrogen, rather than
atomic hydrogen, dominates the composition of matter expelled by stars at the
end of their evolution on the Asymptotic Giant Branch. However sensitive HI
observations can still yield important information on the kinematics and
physical properties of the circumstellar material at large distances from
central stars with heavy mass loss, despite the low abundance of atomic
hydrogen.
We constrain the behavior of the radio luminosity function (RLF) of two classes of active galactic nuclei (AGN) namely AGN of low radio power (LRP) and BL Lac objects. The extrapolation of the observed steep RLFs to low power predicts a space density of such objects that exceeds that of the sources that can harbor them and this requires a break to a shallower slope. For LRP AGN we obtain P_br,LRP > 10^20.5 W/Hz at 1.4 GHz to limit their density to be smaller than that of elliptical galaxies with black hole masses M_BH > 10^7.5 solar masses. By combining this value with the limit derived by the observations the break must occur at P_br,LRP~10^20.5-10^21.5 W/Hz. For BL Lacs we find P_br,BLLAC > 10^23.3 W/Hz otherwise they would outnumber the density of weak-lined and compact radio sources, while the observations indicate P_br,BLLAC < 10^24.5 W/Hz. In the framework of the AGN unified model a low luminosity break in the RLF of LRP AGN must correspond to a break in the RLF of BL Lacs. The ratio between P_br,LRP and P_br,BLLAC is ~10^3, as expected for a jet Doppler factor of ~10.
We report the discovery of a new 21cm HI absorption system using commissioning data from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5-1015.5MHz band of ASKAP we were able to conduct a blind search for the 21cm line in a continuous redshift range between $z =$ 0.4-1.0, which has, until now, remained largely unexplored. The absorption line, detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKSB1740$-$517, is confirmed by optical spectroscopy, using the Gemini South telescope, to be intrinsic to the early-type host galaxy. We detect a broad component at 0.2 per cent of the continuum, demonstrating ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. The [OIII] and [OI] emission lines in the Gemini spectrum are broad and have double-peaked structures, pointing to outflowing ionised gas. Archival data from the XMM-Newton satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the AGN. The absorption profile is complex, with four distinct components ranging in width from 5-300kms$^{-1}$ and fractional depths from 0.2-20 per cent. In addition to systemic HI gas, likely in a regular disc or ring structure, we find evidence for one or two blue shifted clouds and a broad outflow of neutral gas moving at a radial velocity of $v \sim 300$kms$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$yr) is driving surrounding neutral gas in an outflow of $\sim1\mathrm{M}_{\odot}$yr$^{-1}$.
Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.
OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar envelopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO+ and H13CO+ molecular ions and the first detection of SO+, N2H+, and (tentatively) H3O+ in this source. SO+ and H3O+ have not been detected before in CSEs around evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 30 m radio telescope and with Herschel/HIFI. Except for H3O+, all the molecular ions detected in this work display emission lines with broad profiles (FWHM 50-90 km/s), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The narrow profile (FWHM 14 km/s) and high critical densities (>1e6cm-3 ) of the H3O+ transitions observed are consistent with this ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. From rotational diagram analysis, we deduce excitation temperatures of Tex 10-20 K for all ions except for H3O+, which is most consistent with Tex 100 K. Although uncertain, the higher excitation temperature suspected for H3O+ is similar to that recently found for H2O and a few other molecules, which selectively trace a previously unidentified, warm nebular component.The column densities of the molecular ions reported here are in the range Ntot [1-8]x1e13 cm-2, leading to beam-averaged fractional abundances relative to H2 of X(HCO+) 1e-8, X(H13CO+) 2e-9, X(SO+) 4e-9, X(N2H+) 2e-9, and X(H3O+) 7e-9 cm-2. We have performed chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated by cosmic-ray and UV-photon ionization. (abridged).
We constrain the explosion and circumstellar properties at the 2012b event of SN 2009ip based on its late-phase bolometric light curve recently reported. The explosion energy and ejected mass at the 2012b event are estimated as 0.02 Msun and 2e49 erg, respectively. The circumstellar medium is assumed to have two components: an inner shell and an outer wind. The inner shell which is likely created at the 2012a event has 0.2 Msun. The outer wind is created by the wind mass loss before the 2012a mass ejection, and the progenitor is estimated to have had the mass-loss rate about 0.1 Msun/yr with the wind velocity 550 km/s before the 2012a event. The estimated explosion energy and ejected mass indicate that the 2012b event is not caused by a regular supernova.
The Polarized Gamma-ray Observer (PoGOLite) is a balloon-borne instrument that can measure polarization in the energy range 25--240 keV. The instrument adopts an array of well-type "phoswich" detectors in order to suppress backgrounds. Based on the anisotropy of Compton scattering angles resulting from polarized gamma-rays, the polarization of the observed source can be reconstructed. During July 12-26 of 2013, a successful near-circumpolar pathfinder flight was conducted from Esrange, Sweden, to Norilsk, Russia. During this two-week flight, several observations of the Crab were conducted. Here, we present the PoGOLite instrument and summarize the 2013 flight.
Power-law tails are often seen in probability density functions (PDFs) of molecular cloud column densities, and have been attributed to the effect of gravity. We show that extinction PDFs of a sample of five molecular clouds obtained at a few tenths of a parsec resolution, probing extinctions up to A$_{{\mathrm{V}}}$ $\sim$ 10 magnitudes, are very well described by lognormal functions provided that the field selection is tightly constrained to the cold, molecular zone and that noise and foreground contamination are appropriately accounted for. In general, field selections that incorporate warm, diffuse material in addition to the cold, molecular material will display apparent core+tail PDFs. The apparent tail, however, is best understood as the high extinction part of a lognormal PDF arising from the cold, molecular part of the cloud. We also describe the effects of noise and foreground/background contamination on the PDF structure, and show that these can, if not appropriately accounted for, induce spurious tails or amplify any that are truly present.
In 2013 April a new magnetar, SGR 1745-2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the Centre of the Milky Way, Sagittarius A*. SGR 1745-2900 has a surface dipolar magnetic field of ~ 2x10^{14} G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands, with a peak X-ray luminosity L_X ~ 5x10^{35} erg s^{-1}. Here we report on the long-term Chandra (25 observations) and XMM-Newton (8 observations) X-ray monitoring campaign of SGR 1745-2900, from the onset of the outburst in April 2013 until September 2014. This unprecedented dataset allows us to refine the timing properties of the source, as well as to study the outburst spectral evolution as a function of time and rotational phase. Our timing analysis confirms the increase in the spin period derivative by a factor of ~2 around June 2013, and reveals that a further increase occurred between 2013 Oct 30 and 2014 Feb 21. We find that the period derivative changed from 6.6x10^{-12} s s^{-1} to 3.3x10^{-11} s s^{-1} in 1.5 yr. On the other hand, this magnetar shows a slow flux decay compared to other magnetars and a rather inefficient surface cooling. In particular, starquake-induced crustal cooling models alone have difficulty in explaining the high luminosity of the source for the first ~200 days of its outburst, and additional heating of the star surface from currents flowing in a twisted magnetic bundle is probably playing an important role in the outburst evolution.
Circumstellar dust disks have been observed around many nearby stars.
However, many stars are part of binary or multiple stellar systems. A natural
question arises regarding the presence and properties of such disks in systems
with more than one star. To address this, we consider a sample of 449 systems
(spectral types A-M) observed with the Herschel Space Observatory as part of
the DEBRIS program. We have examined the stellar multiplicity of this sample by
gathering information from the literature and performing an adaptive optics
imaging survey at Lick Observatory. Five new companions were revealed with our
program. In total, we identify 188 (42%) binary or multiple star systems. The
multiplicity of the sample is examined with regards to the detection of
circumstellar disks for stars of spectral types AFGK.
In general, disks are less commonly detected around binaries than single
stars, though the disk frequency is comparable among A stars regardless of
multiplicity. However, this sample reveals the period distribution of
disk-bearing binaries is consistent with that of non-disk binaries and with
comparison field samples. We find that the properties of disks in binary
systems are not statistically different from those around single stars.
Although the frequency of disk-bearing FGK binaries may be lower than in single
star systems, the processes behind disk formation and the characteristics of
these disks are comparable among both populations.
The late-type spiral galaxy NGC 6946 is a prime example of molecular gas dynamics driven by `bars within bars'. Here we use data from the BIMA SONG and HERACLES surveys to analyse the structure and stability of its molecular disc. Our radial profiles exhibit a clear transition at distance R ~ 1 kpc from the galaxy centre. In particular, the surface density profile breaks at R ~ 0.8 kpc and is well fitted by a double exponential distribution with scale lengths R_1 ~ 200 pc and R_2 ~ 3 kpc, while the 1D velocity dispersion sigma decreases steeply in the central kpc and is approximately constant at larger radii. The fact that we derive and use the full radial profile of sigma rather than a constant value is perhaps the most novel feature of our stability analysis. We show that the profile of the Q stability parameter traced by CO emission is remarkably flat and well above unity, while the characteristic instability wavelength exhibits clear signatures of the nuclear starburst and inner bar within bar. We also show that CO-dark molecular gas, stars and other factors can play a significant role in the stability scenario of NGC 6946. Our results provide strong evidence that gravitational instability, radial inflow and disc heating have driven the formation of the inner structures and the dynamics of molecular gas in the central kpc.
The demographics of the field of Astronomy, and the gender balance in particular, is an important active area of investigation. A piece of information missing from the discussion is the gender breakdown of the applicant pool for faculty positions. For a sample of 35 tenure-track faculty positions at 25 research universities advertised over the last few years in astronomy and astrophysics, I find that the ratio of female applicants to the total number of applicants is ~0.2, with little dispersion and with no strong dependence on the total number of applicants. Some discussion is provided in the context of the fraction of women at the graduate student, postdoctoral researcher, and assistant professor levels, but strong conclusions are not possible given the limitations of the study. Current and future faculty search committees will likely be interested to compare their numbers to this distribution to decide whether or not they could be doing more to attract an applicant pool that is representative of the community.
The Search for Extra-terrestrial Intelligence (SETI) using radio telescopes is an area of research that is now more than 50 years old. Thus far, both targeted and wide-area surveys have yet to detect artificial signals from intelligent civilisations. In this paper, I argue that the incidence of co-existing intelligent and communicating civilisations is probably small in the Milky Way. While this makes successful SETI searches a very difficult pursuit indeed, the huge impact of even a single detection requires us to continue the search. A substantial increase in the overall performance of radio telescopes (and in particular future wide-field instruments such as the Square Kilometre Array, SKA), provide renewed optimism in the field. Evidence for this is already to be seen in the success of SETI researchers in acquiring observations on some of the world's most sensitive radio telescope facilities via open, peer-reviewed processes. The increasing interest in the dynamic radio sky, and our ability to detect new and rapid transient phenomena such as Fast Radio Bursts (FRB) is also greatly encouraging. While the nature of FRBs is not yet fully understood, I argue they are unlikely to be the signature of distant extra-terrestrial civilisations. As astronomers face a data avalanche on all sides, advances made in related areas such as advanced Big Data analytics, and cognitive computing are crucial to enable serendipitous discoveries to be made. In any case, as the era of the SKA fast approaches, the prospects of a SETI detection have never have been better.
CONTEXT: Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres. AIMS: We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra. METHODS: We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses. RESULTS: Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.
We present the Morphology-Density and Morphology-Radius relations (T-Sigma and T-R, respectively) obtained from the WINGS database of galaxies in nearby clusters. Aiming to achieve the best statistics, we exploit the whole sample of galaxies brighter than MV=-19.5 (5,504 objects), stacking up the 76 clusters of the WINGS survey altogether. Using this global cluster sample, we find that the T-Sigma relation holds only in the inner cluster regions (R<1/3xR200), while the T-R relation keeps almost unchanged over the whole range of local density. A couple of tests and two sets of numerical simulations support the robustness of these results against the effects of the limited cluster area coverage of the WINGS imaging. The above mentioned results hold for all cluster masses (X-ray luminosity and velocity dispersion) and all galaxy stellar masses (M). The strength of the T-Sigma relation (where present) increases with increasing M, while this effect is not found for the T-R relation. Noticeably, the absence/presence of subclustering determines the presence/absence of the T-Sigma relation outside the inner cluster regions, leading us to the general conclusion that the link between morphology and local density is preserved just in dynamically evolved regions. We hypothesize that some mechanism of morphological broadening/redistribution operates in the intermediate/outer regions of substructured (non relaxed) clusters, producing a strong weakening of the T-Sigma relation.
One of the main challenges left for the present and future Cosmic Microwave Background (CMB) experiments is the high precision measurement of the CMB polarization anisotropies at large angular scales. The reionization bump in the CMB polarization power spectra encodes unique informations about the reionization history of the Universe and the inflationary epoch. Such valuable information can be accessed only with an unprecedented accuracy and care on each step of the data analysis and its interpretation. In this paper we present a cross-spectra based approach for the analysis of the CMB data at large angular scales to constrain the reionization optical depth, the tensor to scalar ratio and the amplitude of the primordial scalar perturbations. Using cross-spectra has the advantage to eliminate spurious noise bias and to give a better handle of residual systematics with respect to the pixel-based approach used so far, allowing to efficiently combine the cosmological information encoded in cross-frequency or cross-dataset spectra. We present two solutions to deal with the non-Gaussianity of the Cl estimator distributions at large angular scales: the first relies on an analytical parametrization of the estimator distribution, while the second is based on modification of the Hamimache&Lewis likelihood approximation at large angular scales. The modified HL method (oHL) is extremely powerful as it allows to easily deal with multipole and mode correlations for a combined temperature and polarization analysis. We validate our methods on realistic simulations generated with publicly available specifications, showing that they give consistent results for the constraints of the relevant cosmological parameters in the case of a realistic experimental settings that account for anisotropic correlated noise and incomplete sky coverage.
We highlight recent advances in neutrino astrophysics, the open issues and the interplay with neutrino properties. We emphasize the important progress in our understanding of neutrino flavor conversion in media. We discuss the case of solar neutrinos, of core-collapse supernova neutrinos and of SN1987A, and of the recently discovered ultra-high energy neutrinos whose origin is to be determined.
Highly ionized, z=0 metal absorption lines detected in the X-ray spectra of background active galactic nuclei (AGNs) provide an effective method to probe the hot ($T\sim10^6$ K) gas and its metal content in and around the Milky Way. We present an all-sky survey of the $K_{\alpha}$ transition of the local O VII absorption lines obtained by Voigt-profile fitting archival XMM-Newton observations. A total of 43 AGNs were selected, among which 12 are BL Lac-type AGNs, and the rest are Seyfert 1 galaxies. At above the $3\sigma$ level the local O VII absorption lines were detected in 21 AGNs, among which 7 were newly discovered in this work. The sky covering fraction, defined as the ratio between the number of detections and the sample size, increases from at about 40% for all targets to 100% for the brightest targets, suggesting a uniform distribution of the O VII absorbers. We correlate the line equivalent width with the Galactic coordinates and do not find any strong correlations between these quantities. Some AGNs have warm absorbers that may complicate the analysis of the local X-ray absorber since the recession velocity can be compensated by the outflow velocity, especially for the nearby targets. We discuss the potential impact of the warm absorbers on our analysis. A comprehensive theoretical modelling of the X-ray absorbers will be presented in a later paper.
Above tens of GeV, gamma-ray observations with the Fermi Large Area Telescope (LAT) can be dominated by statistical uncertainties due to the low flux of sources and the limited acceptance. We are developing a new event class which can improve the acceptance: the "Calorimeter-only (CalOnly)" event class. The LAT has three detectors: the tracker, the calorimeter, and the anti-coincidence detector. While the conventional event classes require information from the tracker, the CalOnly event class is meant to be used when there is no usable tracker information. Although CalOnly events have poor angular resolution and a worse signal/background separation compared to those LAT events with usable tracker information, they can increase the instrument acceptance above few tens of GeV, where the performance of Fermi-LAT is limited by low photon statistics. In these proceedings we explain the concept and report some preliminary characteristics of this novel analysis.
Ap Lib is one of the rare Low Synchrotron Peaked blazars detected so far at TeV energies. This type of source is not properly modelled by standard one-zone leptonic Synchrotron-self-Compton (SSC) emission scenarios. The aim of this paper is to study the relevance of additional components which should naturally occur in a SSC scenario for a better understanding of the emission mechanisms, especially at very high energies (VHE). Methods. We use simultaneous data from a multi-wavelength campaign of Planck, Swift-UVOT and Swift-XRT telescopes carried out in February 2010, as well as quasi-simultaneous data of WISE, Fermi and H.E.S.S. taken in 2010. The multi-lambda emission of Ap Lib is modelled by a blob-in-jet SSC scenario including the contribution of the base of the VLBI extended jet, the radiative blob-jet interaction, the accretion disk and its associated external photon field. We show that signatures of a strong parsec-scale jet and of an accretion disk emission are present in the SED. We can link the observationnal VLBI jet features from MOJAVE to parameters expected for a VHE emitting blob accelerated near the jet base. The VHE emission appears to be dominated by the inverse-Compton effect of the blob relativistic electrons interacting with the jet synchrotron radiation. In such scenario Ap Lib appears as an intermediate source between BL Lac objects and Flat Spectrum Radio Quasars. Ap Lib could be a bright representative of a specific class of blazars, in which the parsec-scale jet luminosity is no more negligible compared to the blob and contributes to the high energy emission via inverse Compton processes.
This paper describes a new real-time versatile backend, the Pulsar Ooty Radio Telescope New Digital Efficient Receiver (PONDER), which has been designed to operate along with the legacy analog system of the Ooty Radio Telescope (ORT). PONDER makes use of the current state of the art computing hardware, a Graphical Processing Unit (GPU) and sufficiently large disk storage to support high time resolution real-time data of pulsar observations, obtained by coherent dedispersion over a bandpass of 16 MHz. Four different modes for pulsar observations are implemented in PONDER to provide standard reduced data products, such as time-stamped integrated profiles and dedispersed time series, allowing faster avenues to scientific results for a variety of pulsar studies. Additionally, PONDER also supports general modes of interplanetary scintillation (IPS) measurements and very long baseline interferometry data recording. The IPS mode yields a single polarisation correlated time series of solar wind scintillation over a bandwidth of about four times larger (16 MHz) than that of the legacy system as well as its fluctuation spectrum with high temporal and frequency resolutions. The key point is that all the above modes operate in real time. This paper presents the design aspects of PONDER and outlines the design methodology for future similar backends. It also explains the principal operations of PONDER, illustrates its capabilities for a variety of pulsar and IPS observations and demonstrates its usefulness for a variety of astrophysical studies using the high sensitivity of the ORT.
We propose a scenario for the quasi-periodic oscillations observed in magnetar flares wherein a tangled component of the stellar magnetic field introduces nearly isotropic stress that gives the fluid core of the star an effective shear modulus. In a simple, illustrative model of constant density, the tangled field eliminates the problematic Alfv\'en continuum that would exist in the stellar core for an organized field. For a tangled field energy density comparable to that inferred from the measured dipole fields of $\sim 10^{15}$ G in SGRs 1806-20 and 1900+14, torsional modes exist with fundamental frequencies of about 20 Hz, and mode spacings of $\sim 10$ Hz. For fixed stellar mass and radius, the model has only one free parameter, and can account for {\em every} observed QPO under 160 Hz to within 3 Hz for both SGRs 1806-20 and 1900+14. The combined effects of stratification and crust stresses generally decrease the frequencies of torsional oscillations by $<10$% for overtones and increase the lowest-frequency fundamentals by up to 50%, and so the star can be treated as having constant density to a generally good first approximation. We address the issue of mode excitation by sudden readjustment of the stellar magnetosphere. While the total energy in excited modes is well within the energy budget of giant flares, the surface amplitude is $< 10^{-3}$ of the stellar radius for global oscillations, and decreases strongly with mode frequency. The 626 Hz QPO reported for SGR 1806-20 is particularly problematic to excite beyond a surface amplitude of $10^{-6}$ of the stellar radius.
We investigate the validity of the Cosmological Principle by mapping the cosmological parameters $H_0$ and $q_0$ through the celestial sphere. In our analysis, performed in a low-redshift regime to follow a model-independent approach, we use two compilations of type Ia Supernovae (SNe Ia), namely the Union2.1 and the JLA datasets. Firstly, we show that the angular distributions for both SNe Ia datasets are statistically anisotropic at high confidence level ($p$-value $<$ 0.0001), in particular the JLA sample. Then we find that the cosmic expansion and acceleration are mainly of dipolar type, with maximal anisotropic expansion [acceleration] pointing towards $(l,b) \simeq (326^{\circ},12^{\circ})$ [$(l,b) \simeq (174^{\circ},27^{\circ})$], and $(l,b) \simeq (58^{\circ},-60^{\circ})$ [$(l,b) \simeq (225^{\circ},51^{\circ})$] for the Union2.1 and JLA data, respectively. Secondly, we use a geometrical method to test the hypothesis that the non-uniformly distributed SNe Ia events could introduce anisotropic imprints on the cosmological expansion and acceleration. For the JLA compilation, we found significant correlations between the celestial distribution of data points and the directional studies of $H_0$ and $q_0$, suggesting that these results can be attributed to the intrinsic anisotropy of the sample. In the case of the Union2.1 data, nonetheless, these correlations are less pronounced, and we verify that the dipole asymmetry found in the $H_0$ analyses coincides with the well-known bulk-flow motion of our local group. From these analyses, we conclude that the directional asymmetry on the cosmological parameters maps are mainly either of local origin or due to celestial incompleteness of current SNe Ia samples.
The energy of a free gas of neutrons and protons is well known to be approximately isospin parabolic with a negligibly small quartic term of only $0.45$ MeV at the saturation density of nuclear matter $\rho_0=0.16/\rm{fm}^3$. Using an isospin-dependent single-nucleon momentum distribution including a high (low) momentum tail (depletion) with its shape parameters constrained by recent high-energy electron scattering and medium-energy nuclear photodisintegration experiments as well as the state-of-the-art calculations of the deuteron wave function and the equation of state of pure neutron matter near the unitary limit within several modern microscopic many-body theories, we show for the first time that the kinetic energy of interacting nucleons in neutron-rich nucleonic matter has a significant quartic term of $7.18\pm2.52\,\rm{MeV}$. Such a large quartic term has significant ramifications in determining the equation of state of neutron-rich nucleonic matter using both terrestrial and astrophysical observables.
The statistical properties of the primordial density perturbations has been considered in the past decade as a powerful probe of the physical processes taking place in the early universe. Within the inflationary paradigm, the properties of the bispectrum are one of the keys that serves to discriminate among competing scenarios concerning the details of the origin of cosmological perturbations. However, all of the scenarios, based on the conventional approach to the so-called "quantum-to-classical transition" during inflation, lack the ability to point out the precise physical mechanism responsible for generating the inhomogeneity and anisotropy of our universe starting from and exactly homogeneous and isotropic vacuum state associated with the early inflationary regime. In past works, we have shown that the proposals involving a spontaneous dynamical reduction of the quantum state provide plausible explanations for the birth of said primordial inhomogeneities and anisotropies. In the present manuscript we show that, when considering within the context of such proposals, the characterization of the spectrum and bispectrum turn out to be quite different from those found in the traditional approach, and in particular, some of the statistical features, must be treated in a different way leading to some rather different conclusions.
Many theories of modified gravity, including the well studied Horndeski models, are characterized by a screening mechanism that ensures that standard gravity is recovered near astrophysical bodies. In a recently introduced class of gravitational theories that goes beyond Horndeski, it has been found that new derivative interactions lead to a partial breaking of the Vainshtein screening mechanism inside any gravitational source, although not outside. We study the impact of this new type of deviation from standard gravity on the density profile of a spherically symmetric matter distribution, in the nonrelativistic limit. For simplicity, we consider a polytropic equation of state and derive the modifications to the standard Lane-Emden equations. We also show the existence of a universal upper bound on the amplitude of this type of modified gravity, independently of the details of the equation of state.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
Poynting-Robertson drag has been considered an ineffective mechanism for delivering dust to regions interior to the cool Kuiper belt analogues seen around other Sun-like stars. This conclusion is however based on the very large contrast in dust optical depth between the parent belt and the interior regions that results from the dominance of collisions over drag in systems with detectable cool belts. Here, we show that the levels of habitable zone dust arising from detectable Kuiper belt analogues can be tens to a few hundreds of times greater than the optical depth in the Solar Zodiacal cloud. Dust enhancements of more than a few tens of `zodi' are expected to hinder future Earth-imaging missions, but relatively few undetectable Kuiper belts result in such levels, particularly around stars older than a few Gyr. Thus, current mid to far-IR photometric surveys have already identified most of the 20-25% of nearby stars where P-R drag from outer belts could seriously impact Earth-imaging. The LBTI should easily detect such warm dust around many nearby stars with outer belts, and will provide insight into currently unclear details of the competition between P-R drag and collisions. Given sufficient confidence in future models, the inevitability of P-R drag means that the non-detection of warm dust where detectable levels were expected could be used to infer additional dust removal process, the most likely being the presence of intervening planets.
I will review recent advances in the field of blazars, highlighting the contribution of Swift. Together with other operating satellites (most notably Fermi, but also AGILE, WISE, Planck) and ground based facilities such as Cherenkov telescopes, Swift was (and is) crucial for improving our understanding of blazars. The main advances in the blazar field made possible by Swift includes the opening of the time domain investigation, since there are several sources with hundreds of simultaneous optical, UV and X-ray data taken at different times; the possibility to measure the black hole mass in very powerful blazars, that show clear signs of accretion disk emission; the possibility to classify blazar candidates, through X-ray observations; the finding of the most powerful and distant blazars, emitting strongly in the hard X-ray band accessible to Swift/BAT. All these improvements had and have a great impact on our understanding on how relativistic jets are formed and emit, on their power, and on how the heavy black holes in these systems first formed and grew.
We analyse the spatially-resolved stellar populations of 9 local ($z<0.1$) Brightest Cluster Galaxies (BCGs) observed with VIMOS in IFU mode. Our sample is composed of 7 slow-rotating and 2 fast-rotating BCGs. We do not find a connection between stellar kinematics and stellar populations in this small sample. The BCGs have shallow metallicity gradients (median $\Delta$[Fe/H] $= -0.11\pm0.1$), high central metallicities (median $[$Fe/H]$_{[\alpha/Fe]=0} = 0.13\pm0.07$), and a wide range of central ages (from 5 to 15 Gyr). We propose that the reason for this is diverse evolutionary paths in BCGs. 67 per cent of the sample (6/9) show $\sim 7$ Gyr old central ages, which reflects an active accretion history, and 33 per cent of the sample (3/9) have central ages older than 11 Gyr, which suggest no star formation since $z=2$. The BCGs show similar central stellar populations and stellar population gradients to early-type galaxies of similar mass (M$_{dyn}> 10^{11.3}$M$_{\odot}$) from the ATLAS$^{3D}$ survey (median [Z/H] $= 0.04\pm0.07$, $\Delta$[Z/H] $= -0.19\pm0.1$). However, massive early-type galaxies from ATLAS$^{3D}$ have consistently old ages (median Age $=12.0\pm3.8$Gyr). We also analyse the close massive companion galaxies of two of the BCGs. These galaxies have similar stellar populations to their respective BCGs.
Recent efforts in cosmic ray (CR) confinement and transport theory are discussed. Three problems are addressed as being crucial for understanding the present day observations and their possible telltale signs of the CR origin. The first problem concerns CR behavior right after their release from a source, such as a supernova remnant (SNR). At this phase the CRs are confined near the source by self-emitted Alfven waves. The second is the problem of diffusive propagation of CRs through the turbulent ISM. This is a seemingly straightforward and long-resolved problem, but it remains controversial and reveals paradoxes. A resolution based on the Chapman-Enskog asymptotic CR transport analysis, that also includes magnetic focusing, is suggested. The third problem is about a puzzling sharp ($\sim10^{\circ}$) anisotropies in the CR arrival directions that might bear on important clues of their transport between the source and observer. The overarching goal is to improve our understanding of all aspects of the CR's source escape and ensuing propagation through the galaxy to the level at which their sources can be identified observationally.
Rotating neutron stars, or pulsars, are plausibly the source of power behind many astrophysical systems, such as gamma-ray bursts, supernovae, pulsar wind nebulae and supernova remnants. In the past several years, 3D numerical simulations made it possible to compute pulsar spindown luminosity from first principles and revealed that oblique pulsar winds are more powerful than aligned ones. However, what causes this enhanced power output of oblique pulsars is not understood. In this work, using time-dependent 3D magnetohydrodynamic (MHD) and force-free simulations, we show that, contrary to the standard paradigm, the open magnetic flux, which carries the energy away from the pulsar, is laterally non-uniform. We argue that this non-uniformity is the primary reason for the increased luminosity of oblique pulsars. To demonstrate this, we construct simple analytic descriptions of aligned and orthogonal pulsar winds and combine them to obtain an accurate 3D description of the pulsar wind for any obliquity. Our approach describes both the warped magnetospheric current sheet and the smooth variation of pulsar wind properties outside of it. We find that generically the magnetospheric current sheet separates plasmas that move at mildly relativistic velocities relative to each other. This suggests that the magnetospheric reconnection is a type of driven, rather than free, reconnection. The jump in magnetic field components across the current sheet decreases with increasing obliquity, which could be a mechanism that reduces dissipation in near-orthogonal pulsars. Our analytical description of the pulsar wind can be used for constructing models of pulsar gamma-ray emission, pulsar wind nebulae, and magnetar-powered core-collapse gamma-ray bursts and supernovae.
Aims: We investigate the effect of ram pressure stripping (RPS) on
simulations of merging pairs of gas-rich spiral galaxies. Our goal is to
provide an estimate of the combined effect of merging and RPS on stripping
efficiency and star formation rate.
Methods: We make use of the combined N-body/hydrodynamic code GADGET-2. In
our simulations, we vary mass ratios between 1:4 and 1:8 in a binary merger. We
sample different geometric configurations of the merging systems (edge-on and
face-on mergers, different impact parameters). Furthermore, we vary the
properties of the intracluster medium (ICM) in rough steps: The speed of the
merging system relative to the ICM between 500 and 1000 km/s, the ICM density
between $10^{-29}$ and $10^{-27}$ g/cm$^3$, and the ICM direction relative to
the mergers' orbital plane. Ram pressure is kept constant within a simulation
time period, as is the ICM temperature of $10^7$ K. Each simulation in the ICM
is compared to simulations of the merger in vacuum and the non-merging galaxies
with acting ram pressure.
Results: Averaged over the simulation time (1 Gyr) the merging pairs show a
negligible 5% enhancement in SFR, when compared to single galaxies under the
same environmental conditions. The SFRs peak at the time of the galaxies first
fly-through. There, our simulations show SFRs of up to 20 M$_{\odot}$/yr
(compared to 3 M$_{\odot}$/yr of the non-merging galaxies in vacuum). In the
most extreme case, this constitutes a short-term ($<50$ Myr) SFR increase of
50% over the non-merging galaxies experiencing ram pressure. The wake of
merging galaxies in the ICM typically has a third to half the star mass seen in
the non-merging galaxies and 5% to 10 % less gas mass. The joint effect of RPS
and merging, according to our simulations, is not significantly different from
pure ram pressure effects.
X-ray bright quasars might be used to trace dust in the circumgalactic and intergalactic medium through the phenomenon of X-ray scattering, which is observed around Galactic objects whose light passes through a sufficient column of interstellar gas and dust. Of particular interest is the abundance of grey dust larger than 0.1 um, which is difficult to detect at other wavelengths. To calculate X-ray scattering from large grains, one must abandon the traditional Rayleigh-Gans approximation. The Mie solution for the X-ray scattering optical depth of the Universe is ~1%. This presents a great difficulty for distinguishing dust scattered photons from the point source image of Chandra, which is currently unsurpassed in imaging resolution. The variable nature of AGN offers a solution to this problem, as scattered light takes a longer path and thus experiences a time delay with respect to non-scattered light. If an AGN dims significantly (> 3 dex) due to a major feedback event, the Chandra point source image will be suppressed relative to the scattering halo, and an X-ray echo or ghost halo may become visible. I estimate the total number of scattering echoes visible by Chandra over the entire sky: N_ech ~ 10^3 (nu_fb / yr^-1), where nu_fb is the characteristic frequency of feedback events capable of dimming an AGN quickly.
Before the launch of the Fermi satellite only two classes of Active Galactic Nuclei (AGN) were known to generate relativistic jets and thus to emit up to the gamma-ray energy range: blazars and radio galaxies, both hosted in giant elliptical galaxies. The first four years of observations by the Large Area Telescope (LAT) on board Fermi confirmed that these two populations represent the most numerous identified sources in the extragalactic gamma-ray sky, but the discovery of variable gamma-ray emission from 5 radio-loud Narrow-Line Seyfert 1 (NLSy1) galaxies revealed the presence of a possible emerging third class of AGN with relativistic jets. Considering that NLSy1 are thought to be hosted in spiral galaxies, this finding poses intriguing questions about the nature of these objects, the knowledge of the development of relativistic jets, and the evolution of radio-loud AGN. In this context, the study of the radio-loud NLSy1 from radio to gamma-rays has received increasing attention. Here we discuss the radio-to-gamma-rays properties of the gamma-ray emitting NLSy1, also in comparison with the blazar scenario.
We present the best sensitivity and angular resolution maps of the molecular disk and outflow of Mrk231, obtained with IRAM/PdBI, and an analysis of archival Chandra and NuSTAR data. We constrain the physical properties of both the molecular disk and outflow, the presence of a highly-ionized ultra-fast nuclear wind, and their connection. The CO(2-1) outflow has a size of ~1 kpc, and extends in all directions around the nucleus, being more prominent along the south-west to north-east direction, suggesting a wide-angle biconical geometry. Its maximum projected velocity is nearly constant out to ~1 kpc, thus implying that the density of the outflowing material must decrease from the nucleus outwards as ~ r^-2. This suggests that either a large part of the gas leaves the flow during its expansion, or that the bulk of the outflow has not yet reached ~1 kpc, implying a limit on its age of ~ 1 Myr. The mass and energy rates of the molecular outflow are dM/dt(OF)=[500-1000] Msun/yr and dE(kin,OF)/dt=[7-10] 10^43 erg/s, its total kinetic energy is E(kin,OF)>~E(disk). Remarkably, our analysis reveals a nuclear ultra-fast outflow (UFO) with velocity ~-20000 km/s, dM(UFO)/dt=0.3-1.6 MSun/yr, and momentum load dP(UFO)/dt/(Lbol/c)=0.25-1.2. We find dE(kin,UFO)/dt~dE(kin,OF)/dt as predicted for outflows undergoing an energy conserving expansion. This suggests that most of the UFO kinetic energy is transferred to mechanical energy of the kpc-scale outflow, strongly supporting that the energy released during accretion of matter onto super-massive black holes is the ultimate driver of giant massive outflows.The momentum flux dP(OF)/dt derived for the large scale outflows in Mrk 231 enables us to estimate a momentum boost dP(OF)/dP(UFO)~ 30-50. The ratios dE(kin, UFO)/dt/L(bol,AGN) = 0.8-4% and dE(kin,OF)/dt/L(bol,AGN) = 1-3% agree with the requirements of the most popular models of AGN feedback.
We perform a 3D multi-probe analysis of the rich galaxy cluster A1689 by combining improved weak-lensing data from new BVRi'z' Subaru/Suprime-Cam observations with strong-lensing, X-ray, and Sunyaev-Zel'dovich effect (SZE) data sets. We reconstruct the projected matter distribution from a joint weak-lensing analysis of 2D shear and azimuthally integrated magnification constraints, the combination of which allows us to break the mass-sheet degeneracy. The resulting mass distribution reveals elongation with axis ratio ~0.7 in projection. When assuming a spherical halo, our full weak-lensing analysis yields a projected concentration of $c_{200c}^{2D}=8.9\pm 1.1$ ($c_{vir}^{2D}\sim 11$), consistent with and improved from earlier weak-lensing work. We find excellent consistency between weak and strong lensing in the region of overlap. In a parametric triaxial framework, we constrain the intrinsic structure and geometry of the matter and gas distributions, by combining weak/strong lensing and X-ray/SZE data with minimal geometric assumptions. We show that the data favor a triaxial geometry with minor-major axis ratio 0.39+/-0.15 and major axis closely aligned with the line of sight (22+/-10 deg). We obtain $M_{200c}=(1.2\pm 0.2)\times 10^{15} M_{\odot}/h$ and $c_{200c}=8.4\pm 1.3$, which overlaps with the $>1\sigma$ tail of the predicted distribution. The shape of the gas is rounder than the underlying matter but quite elongated with minor-major axis ratio 0.60+/-0.14. The gas mass fraction within 0.9Mpc is 10^{+3}_{-2}%. The thermal gas pressure contributes to ~60% of the equilibrium pressure, indicating a significant level of non-thermal pressure support. When compared to Planck's hydrostatic mass estimate, our lensing measurements yield a spherical mass ratio of $M_{Planck}/M_{GL}=0.70\pm 0.15$ and $0.58\pm 0.10$ with and without corrections for lensing projection effects, respectively.
We present analytic calculations of the electromagnetic torques acting on a magnetic neutron star rotating in vacuum, including near-zone torques associated with the inertia of dipole and quadrupole magnetic fields. We incorporate these torques into the rotational dynamics of a rigid-body neutron star, and show that the effects of the inertial torque can be understood as a modification of the moment of inertia tensor of the star. We apply our rotational dynamics equation to the Crab pulsar, including intrinsic distortions of the star and various electromagnetic torques, to investigate the possibility that the counter-alignment of the magnetic inclination angle, as suggested by recent observations, could be explained by pulsar precession. We find that if the effective principal axis of the pulsar is nearly aligned with either the magnetic dipole axis or the rotation axis, then precession may account for the observed counter-alignment over decade timescales. Over the spindown timescale of the pulsar, the magnetic inclination angle always decreases.
We present direct evidence for the detection of the main energy release site in a non-eruptive solar flare, SOL2013-11-09T06:38UT. This GOES C2.7 event was characterised by two flaring ribbons and a compact, bright coronal source located between them, which is the focus of our study. We use imaging from SDO/AIA, and imaging spectroscopy from RHESSI to characterise the thermal and non-thermal emission from the coronal source, and EUV spectroscopy from the Hinode/EIS, which scanned the coronal source during the impulsive peak, to analyse Doppler shifts in Fe XII and Fe XXIV emission lines, and determine the source density. The coronal source exhibited an impulsive emission lightcurve in all AIA filters during the impulsive phase. RHESSI hard X-ray images indicate both thermal and non-thermal emission at the coronal source, and its plasma temperature derived from RHESSI imaging spectroscopy shows an impulsive rise, reaching a maximum at 12-13 MK about 10 seconds prior to the hard X-ray peak. High redshifts associated with this bright source indicate downflows of 40-250 km/s at a broad range of temperatures, interpreted as loop shrinkage and/or outflows along the magnetic field. Outflows from the coronal source towards each ribbon are also observed by AIA images at 171, 193, 211, 304 and 1600 A. The electron density of the source obtained from a Fe XIV line pair is $10^{11.50}$ which is collisionally thick to electrons with energy up to 45-65 keV, responsible for the source's non-thermal X-ray emission. We conclude that the bright coronal source is the location of the main release of magnetic energy in this flare, with a geometry consistent with component reconnection between crossing, current-carrying loops. We argue that the energy that can be released via reconnection, based on observational estimates, can plausibly account for the non-thermal energetics of the flare.
The EXTraS project (Exploring the X-ray Transient and variable Sky) will harvest the hitherto unexplored temporal domain information buried in the serendipitous data collected by the European Photon Imaging Camera (EPIC) instrument onboard the ESA XMM-Newton X-ray observatory since its launch. This will include a search for fast transients, as well as a search and characterization of variability (both periodic and aperiodic) in hundreds of thousands of sources spanning more than nine orders of magnitude in time scale and six orders of magnitude in flux. X-ray results will be complemented by multiwavelength characterization of new discoveries. Phenomenological classification of variable sources will also be performed. All our results will be made available to the community. A didactic program in selected High Schools in Italy, Germany and the UK will also be implemented. The EXTraS project (2014-2016), funded within the EU/FP7 framework, is carried out by a collaboration including INAF (Italy), IUSS (Italy), CNR/IMATI (Italy), University of Leicester (UK), MPE (Germany) and ECAP (Germany).
The radially averaged metallicity distribution of the ISM and the young stellar population of a sample of 20 disk galaxies is investigated by means of an analytical chemical evolution model which assumes constant ratios of galactic wind mass loss and accretion mass gain to star formation rate. Based on this model the observed metallicities and their gradients can be described surprisingly well by the radially averaged distribution of the ratio of stellar mass to ISM gas mass. The comparison between observed and model predicted metallicity is used to constrain the rate of mass loss through galactic wind and accretion gain in units of the star formation rate. Three groups of galaxies are found: galaxies with either mostly winds and only weak accretion, or mostly accretion and only weak winds, and galaxies where winds are roughly balanced by accretion. The three groups are distinct in the properties of their gas disks. Galaxies with approximately equal rates of mass-loss and accretion gain have low metallicity, atomic hydrogen dominated gas disks with a flat spatial profile. The other two groups have gas disks dominated by molecular hydrogen out to 0.5 to 0.7 isophotal radii and show a radial exponential decline, which is on average steeper for the galaxies with small accretion rates. The rates of accretion (<1.0 x SFR) and outflow (<2.4 x SFR) are relatively low. The latter depend on the calibration of the zero point of the metallicity determination from the use of HII region strong emission lines.
Primordial black holes (PBHs) are black holes which may have formed very early on during the radiation dominated era in the early universe. We present here a method by which the large scale perturbations in the density of primordial black holes may be used to place tight constraints on non-gaussianity if PBHs account for dark matter (DM). The presence of local-type non-gaussianity is known to have a significant effect on the abundance of primordial black holes, and modal coupling from the observed CMB scale modes can significantly alter the number density of PBHs that form within different regions of the universe, which appear as DM isocurvature modes. Using the recent \emph{Planck} constraints on isocurvature perturbations, we show that PBHs are excluded as DM candidates for even very small local-type non-gaussianity, $|f_{NL}|\approx0.001$ and remarkably the constraint on $g_{NL}$ is almost as strong. Even small non-gaussianity is excluded if DM is composed of PBHs. If local non-Gaussianity is ever detected on CMB scales, the constraints on the fraction of the universe collapsing into PBHs (which are massive enough to have not yet evaporated) will become much tighter.
What is the meaning of the Fermi Paradox -- are we alone or is starfaring rare? Can general relativity be united with quantum mechanics? The searches for answers to these questions could intersect. It is known that an accelerator capable of energizing particles to the Planck scale requires cosmic proportions. The energy required to run a Planck accelerator is also cosmic, of order 100 M_sun c^2 for a hadron collider, because the natural cross section for Planck physics is so tiny. If aliens are interested in fundamental physics, they could resort to cosmic engineering for their experiments. These colliders are detectable through the vast amount of "pollution" they produce, motivating a YeV SETI program. I investigate what kinds of radiation they would emit in a fireball scenario, and the feasibility of detecting YeV radiation at Earth, particularly YeV neutrinos. Although current limits on YeV neutrinos are weak, Kardashev 3 YeV neutrino sources appear to be at least 30--100 Mpc apart on average, if they are long-lived and emit isotropically. I consider the feasibility of much larger YeV neutrino detectors, including an acoustic detection experiment that spans all of Earth's oceans, and instrumenting the entire Kuiper Belt. Any detection of YeV neutrinos implies an extraordinary phenomenon at work, whether artificial and natural. Searches for YeV neutrinos from any source are naturally commensal, so a YeV neutrino SETI program has value beyond SETI itself, particularly in limiting topological defects. I note that the Universe is very faint in all kinds of nonthermal radiation, indicating that cosmic engineering is extremely rare.
We present preliminary results of our analysis of the {\it Fermi}-LAT data from the direction of NGC 4151. We find a new gamma-ray source with a statistical significance sigma > 5, shifted by 0.5degr from the position of NGC 4151. Apparently, the source was bright only during a 1.5-year period between December 2011 and June 2013 and it strongly contaminated the signal from NGC 4151. Therefore, we neglect this period in our analysis. We find two additional, persistent gamma-ray sources with high sigma, shifted from NGC 4151 by ~1.5degr and 5degr, whose presence has been recently confirmed in the Third Fermi Catalog. After subtracting the above sources, we still see a weak residual, with sigma ~< 3, at the position of NGC 4151. We derive an upper limit (UL) for the gamma-ray flux from NGC 4151 and we compare it with predictions of the ADAF model which can explain the X-ray observations of this object. We find that the Fermi UL strongly constrains non-thermal acceleration processes in hot flows as well as the values of some crucial parameters. Here we present the comparison with the hot flow models in which heating of electrons is dominated by Coulomb interactions with hot protons. In such a version of the model, the gamma-ray UL, combined with the X-ray data, constrains the energy content in the non-thermal component of proton distribution to at most a few per cent, rules out a weak (sub-equipartition) magnetic field and favors a rapid rotation of the supermassive black hole.
Current and future astronomical survey facilities provide a remarkably rich opportunity for transient astronomy, combining unprecedented fields of view with high sensitivity and the ability to access previously unexplored wavelength regimes. This is particularly true of LOFAR, a recently-commissioned, low-frequency radio interferometer, based in the Netherlands and with stations across Europe. The identification of and response to transients is one of LOFAR's key science goals. However, the large data volumes which LOFAR produces, combined with the scientific requirement for rapid response, make automation essential. To support this, we have developed the LOFAR Transients Pipeline, or TraP. The TraP ingests multi-frequency image data from LOFAR or other instruments and searches it for transients and variables, providing automatic alerts of significant detections and populating a lightcurve database for further analysis by astronomers. Here, we discuss the scientific goals of the TraP and how it has been designed to meet them. We describe its implementation, including both the algorithms adopted to maximize performance as well as the development methodology used to ensure it is robust and reliable, particularly in the presence of artefacts typical of radio astronomy imaging. Finally, we report on a series of tests of the pipeline carried out using simulated LOFAR observations with a known population of transients.
We searched for a correlation between the two anomalous properties of K giants: Li enhancement and IR excess from an unbiased survey of a large sample of RGB stars. A sample of 2000 low-mass K giants with accurate astrometry from the Hipparcos catalog was chosen for which Li abundances have been determined from low-resolution spectra. Far-infrared data were collected from the $WISE$ and $IRAS$ catalogs. To probe the correlation between the two anomalies, we supplemented 15 Li-rich K giants discovered from this sample with 25 known Li-rich K giants from other studies. Dust shell evolutionary models and spectral energy distributions were constructed using the code DUSTY to estimate different dust shell properties, such as dust evolutionary time scales, dust temperatures, and mass-loss rates. Among 2000 K giants, we found about two dozen K giants with detectable far-IR excess, and surprisingly, none of them are Li-rich. Similarly, the 15 new Li-rich K giants that were identified from the same sample show no evidence of IR excess. Of the total 40 Li-rich K giants, only 7 show IR excess. Important is that K giants with Li enhancement and/or IR excess begin to appear only at the bump on the RGB. Results show that K giants with IR excess are very rare, similar to K giants with Li enhancement. This may be due to the rapid differential evolution of dust shell and Li depletion compared to RGB evolutionary time scales. We also infer from the results that during the bump evolution, giants probably undergo some internal changes, which are perhaps the cause of mass-loss and Li-enhancement events. However, the available observational results do not ascertain that these properties are correlated. That a few Li-rich giants have IR excess seems to be pure coincidence.
We have performed integral field spectroscopy of the planetary nebulae Hen 3-1333 (PNG332.9-09.9) and Hen 2-113 (PNG321.0+03.9), which are unusual in exhibiting dual-dust chemistry and multipolar lobes but also ionized by late-type [WC 10] central stars. The spatially resolved velocity distributions of the H$\alpha$ emission line were used to determine their primary orientations. The integrated H$\alpha$ emission profiles indicate that Hen 3-1333 and Hen 2-113 expand with velocities of ~ 32 and 23 km/s, respectively. The Hubble Space Telescope observations suggest that these planetary nebulae have two pairs of tenuous lobes extending upwardly from their bright compact cores. From three-dimensional geometric models, the primary lobes of Hen 3-1333 and Hen 2-113 were found to have inclination angles of about -30$^{\circ}$ and 40$^{\circ}$ relative to the line of sight, and position angles of -15$^{\circ}$ and 65$^{\circ}$ measured east of north in the equatorial coordinate system, respectively.
Astronomical optical interferometers (OI) sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid atmospheric perturbations, the phases of the complex Fourier samples (visibilities) cannot be directly exploited , and instead linear relationships between the phases are used (phase closures and differential phases). Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic OI instruments are now paving the way to multiwavelength imaging. This paper presents the derivation of a spatio-spectral ("3D") image reconstruction algorithm called PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm is able to solve large scale problems. It relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also from differential phases, which help to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm.
Long term observations by Brook et al. reveal that the derivative of rotational frequency of PSR J0738-4042 changed abruptly in 2005. Originally, the spin-down rate was relatively stable, with the rotational frequency derivative of $-1.14 \times 10^{-14}~\rm s^{-2}$. After September 2005, the derivative began to rise up. About 1000 days later, it arrived at another relatively stable value of about $-0.98 \times 10^{-14}~\rm s^{-2}$, indicating that the pulsar is spinning-down relatively slowly. To explain the observed spin-down rate change, we resort to an asteroid disrupted by PSR J0738-4042. In our model, the orbital angular momentum of the asteroid is assumed to be parallel to that of the rotating pulsar, so that the pronounced reduction in the spin-down rate can be naturally explained as due to the transfer of the angular momentum from the disrupted material to the central pulsar. The derived magnetospheric radius is about $4.0 \times 10^{9}$ cm, which is smaller than the tidal disruption radius ($4.9 \times 10^{10}$ cm). Our model is self-consistent. It is shown that the variability of the spin-down rate of PSR J0738-4042 can be quantitatively accounted for by the accretion from the asteroid disrupted by the central pulsar.
A large variation in 14 C around AD 775 has been considered to be caused by one or more solar super-flares within one year. We critically review all known aurora reports from Europe as well as the Near, Middle, and Far East from AD 731 to 825 and find 39 likely true aurorae plus four more potential aurorae and 24 other reports about halos, meteors, thunderstorms etc., which were previously misinterpreted as aurorae or misdated; we assign probabilities for all events according to five aurora criteria. We find very likely true aurorae in AD 743, 745, 762, 765, 772, 773, 793, 796, 807, and 817. There were two aurorae in the early 770s observed near Amida (now Diyarbakir in Turkey near the Turkish-Syrian border), which were not only red, but also green-yellow - being at a relatively low geo-magnetic latidude, they indicate a relatively strong solar storm. However, it cannot be argued that those aurorae (geo-magnetical latitude 43 to 50 deg, considering five different reconstructions of the geo-magnetic pole) could be connected to one or more solar super-flares causing the 14 C increase around AD 775: There are several reports about low- to mid-latitude aurorae at 32 to 44 deg geo-magnetical latitude in China and Iraq; some of them were likely observed (quasi-)simultaneously in two of three areas (Europe, Byzantium/Arabia, East Asia), one lasted several nights, and some indicate a particulary strong geo-magnetic storm (red colour and dynamics), namely in AD 745, 762, 793, 807, and 817 - always without 14 C peaks. We use 39 likely true aurorae as well as historic reports about sunspots together with the radiocarbon content from tree rings to reconstruct solar activity: From about AD 733 to 823, we see at least nine Schwabe cycles ...
Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.
X-ray spectra in the range $1.5-8.5$~keV have been analyzed for 526 large flares detected with the Solar Assembly for X-rays (SAX) on the Mercury {\em MESSENGER} spacecraft between 2007 and 2013. For each flare, the temperature and emission measure of the emitting plasma were determined from the spectrum of the continuum. In addition, with the SAX energy resolution of 0.6 keV (FWHM) at 6~keV, the intensities of the clearly resolved Fe-line complex at 6.7~keV and the Ca-line complex at 3.9~keV were determined, along with those of unresolved line complexes from S, Si, and Ar at lower energies. Comparisons of these line intensities with theoretical spectra allow the abundances of these elements relative to hydrogen to be derived, with uncertainties due to instrument calibration and the unknown temperature distribution of the emitting plasma. While significant deviations are found for the abundances of Fe and Ca from flare to flare, the abundances averaged over all flares are found to be enhanced over photospheric values by factors of $1.66 \pm 0.34$ (Fe), $3.89~\pm~0.76$ (Ca), $1.23~\pm~0.45$ (S), $1.64~\pm~0.66$ (Si), and $2.48~\pm~0.90$ (Ar). These factors differ from previous reported values for Fe and Si at least. They suggest a more complex relation of abundance enhancement with the first ionization potential (FIP) of the element than previously considered, with the possibility that fractionation occurs in flares for elements with a FIP of less than $\sim$7~eV rather than $\sim10$~eV.
The dipolar model \cite{Gordon:2005ai} has attracted much interest because it may phenomenologically explain the CMB hemispherical power asymmetry found in the WMAP and Planck data. Since such a model explicitly breaks isotropy at large angular scales it is natural to wonder whether it can also explain other CMB directional anomalies. Focusing on the low $\ell$ alignments and assuming $\Lambda$CDM, we confirm that the quadrupole/octupole and the dipole/quadrupole/octupole alignments are anomalous with a significance up to $99.9\%$ C.L., for both WMAP and Planck data. Moreover, we show for the first time that such features are anomalous also in the dipolar model, roughly at the same level as in $\Lambda$CDM. We conclude that the dipolar model does not provide a better fit to the data than the $\Lambda$CDM.
A paper about variable stars with possible multiple occurrence in the VSX catalogue is presented. Our main criteria for identification of such duplicities were the angular distance among stars (below 1 arcmin) and close periods of objects. In our approach, we also considered double or half values of periods to reveal possible misclassification among stars with similar light curve shapes. The probability of false identification is expressed by the parameter R giving the relative difference between periods. We found 1487 pairs of stars in angular distance lower than 1 arcmin with period difference R lower than 0.1 %, which are high-probable candidates on duplicates. From this sample, 354 pairs have exactly the same periods (R = 0.0 %) and should be considered as definite duplicates. The main contribution of certain duplicates comes from the Catalina Sky Survey (73 pairs have two names with CSS acronym) and from the BEST projects (71 pairs). Distribution of identified duplicates on the sky is not homogeneous but contains surprising depression in Galactic plane.
The discovery of OI atoms and CII ions in the upper atmosphere of HD 209458b, made with the Hubble Space Telescope Imaging Spectrograph (STIS) using the G140L grating, showed that these heavy species fill an area comparable to the planet's Roche lobe. The derived ~10% transit absorption depths require super-thermal processes and/or supersolar abundances. From subsequent Cosmic Origins Spectrograph (COS) observations, CII absorption was reported with tentative velocity signatures, and absorption by SiIII ions was also claimed in disagreement with a negative STIS G140L detection. Here, we revisit the COS dataset showing a severe limitation in the published results from having contrasted the in-transit spectrum against a stellar spectrum averaged from separate observations, at planetary phases 0.27, 0.72, and 0.49. We find variable stellar SiIII and CII emissions that were significantly depressed not only during transit but also at phase 0.27 compared to phases 0.72 and 0.49. Their respective off-transit 7.5 and 3.1% flux variations are large compared to their reported 8.2+/-1.4% and 7.8+/-1.3% transit absorptions. Significant variations also appear in the stellar line shapes, questioning reported velocity signatures. We furthermore present archive STIS G140M transit data consistent with no SiIII absorption, with a negative result of 1.7+/-18.7 including ~15% variability. Silicon may still be present at lower ionization states, in parallel with the recent detection of extended magnesium, as MgI atoms. In this frame, the firm detection of OI and CII implying solar or supersolar abundances contradicts the recent inference of potential x20-125 subsolar metallicity for HD 209458b.
Hot subdwarfs are considered to be the compact helium cores of red giants, which lost almost their entire hydrogen envelope. What causes this enormous mass loss is still unclear. Binary interactions are invoked and a significant fraction of the hot subdwarf population is indeed found in close binaries. In a large project we search for the close binary sdBs with the most and the least massive companions. Significantly enhancing the known sample of close binary sdBs we performed the first comprehensive study of this population. Triggered by the discovery of two sdB binaries with close brown dwarf companions in the course of this project, we were able to show that the interaction of stars with substellar companions is an important channel to form sdB stars. Finally, we discovered a unique and very compact binary system consisting of an sdB and a massive white dwarf, which qualifies as progenitor candidate for a supernova type Ia. In addition to that, we could connect those explosions to the class of hypervelocity hot subdwarf stars, which we consider as the surviving companions of such events. Being the stripped cores of red giants, hot subdwarfs turned out to be important markers of peculiar events in stellar evolution ranging all the way from star-planet interactions to the progenitors of stellar explosions used to measure the expansion of our Universe.
Detection of the cosmological 21-cm signal coming from the epoch of reionization (EoR) is challenging in the presence of astrophysical foregrounds and direction (in)dependent systematic errors of the instrument. Even after removing the foregrounds, the residual Stokes I maps contain, in addition to the system noise, the polarized foreground leaked from Stokes Q, U due to systematic and random errors which can mimic the EoR signal. Here we discuss the systematic errors, especially the primary beam, of LOFAR and present realistic simulations of the leakages caused by them. We made a Stokes Q, U sky model of the Galactic diffuse emission based on the LOFAR observations of the 3C196 field, simulated the full-Stokes visibilities that would be produced by LOFAR in the presence of its nominal model beam, created RM-cubes and the cylindrically and spherically averaged 3D power spectra (PS), and compared them with the PS of a simulated EoR signal. From the spherical PS, we found that at 134-166 MHz, within the central 4 deg of the field the (Q, U ) to I leakage power is lower than the EoR signal at k < 0.3 /Mpc. The leakage was found to be localized around a Faraday depth of 0, and in the cylindrical PS, the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3%, both of which could be utilized in the removal of leakage. Moreover, we could define an 'EoR window' in terms of the polarization leakage in the cylindrical PS above the PSF-induced wedge and below $k_\parallel\sim 0.5$ /Mpc, and the window extended up to $k_\parallel\sim 1$ /Mpc at all $k_\perp$ when 70% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim$ 4 deg in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 /Mpc.
Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.
We present simulations of collapsing filaments studying the impact of turbulence and magnetic field morphologies on their evolution and star formation properties. We vary the mass per unit length of the filaments as well as the orientation of the magnetic field with respect to the major axis. We find that the filaments, which have no or a perpendicular magnetic field, typically reveal a smaller width than the universal width of 0.1 pc proposed by e.g. Arzoumanian et al. 2011. We show that this also holds in the presence of supersonic turbulence and that accretion driven turbulence is too weak to stabilize the filaments along their radial direction. On the other hand, we find that a magnetic field that is parallel to the major axis can stabilize the filament against radial collapse resulting in widths of 0.1 pc. Furthermore, depending on the filament mass and magnetic field configuration, gravitational collapse and fragmentation in filaments occurs either in an edge-on way, uniformly distributed across the entire length, or in a mixed way. In the presence of initially moderate density perturbations, a centralized collapse towards a common gravitational centre occurs. Our simulations can thus reproduce different modes of fragmentation observed recently in star forming filaments. Moreover, we find that turbulent motions influence the distance between individual fragments along the filament, which does not always match the results of a Jeans analysis.
Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamical friction. However, in simulations where the CND is fragmented in high density gaseous clumps (clumpy CND), the orbits of the SMBHs are erratically perturbed by the gravitational interaction with these clumps, delaying, in some cases, the orbital decay of the SMBHs. The densities of these gaseous clumps in our simulations and in recent studies of clumpy CNDs are significantly higher than the observed density of molecular clouds in isolated galaxies or ULIRGs, thus, we expect that SMBH orbits are perturbed less in real CNDs than in the simulated CNDs of this study and other recent studies. We also find that the migration timescale has a weak dependence on the star formation rate of the CND. Furthermore, the migration timescale of a SMBH pair in a star-forming clumpy CND is at most a factor three longer than the migration timescale of a pair of SMBHs in a CND modeled with more simple gas physics. Therefore, we estimate that the migration timescale of the SMBHs in a clumpy CND is on the order of $10^7$ yrs.
NASA's Fermi space telescope has provided us with a bountiful new population of gamma-ray sources following its discovery of 150 new gamma-ray pulsars. One common feature exhibited by all of these pulsars is the form of their spectral energy distribution, which can be described by a power law followed by a spectral break occurring between $\sim$1 and $\sim$8 GeV. The common wisdom is that the break is followed by an exponential cut-off driven by radiation/reaction-limited curvature emission. The discovery of pulsed gamma rays from the Crab pulsar, the only pulsar so far detected at very high energies (E$>$100GeV), contradicts this "cutoff" picture. Here we present a new stacked analysis with an average of 4.2 years of data on 115 pulsars published in the 2nd LAT catalog of pulsars. This analysis is sensitive to low-level $\sim$100 GeV emission which cannot be resolved in individual pulsars but can be detected from an ensemble.
A systematic dynamical system approach is applied to study the cosmology of anisotropic Bianchi I universes in which a vector field is assumed to operate on a disformal frame. This study yields a number of new fixed points, among which anisotropic scaling solutions. Within the simplifying assumption of (nearly) constant-slope potentials these are either not stable attractors, do not describe accelerating expansion or else they feature too large anisotropies to be compatible with observations. Nonetheless, some solutions do have an appeal for cosmological applications in that isotropy is retained due to rapid oscillations of the vector field.
We use classical and quantum Monte Carlo techniques to study the static structure function $S(q)$ of a one-component ion lattice and use it to calculate the thermal conductivity $\kappa$ of high-density solid matter expected in the neutron star crust. We also calculate the phonon spectrum using the dynamic-matrix method and use it to obtain $\kappa$ in the one-phonon approximation. We compare the results obtained with these methods and assess the validity of some commonly used approximations in the literature. We find that quantum effects became relevant for the calculation of $\kappa$ when the temperature $T\lesssim 0.3~\Omega_\mathrm{P}$, where $\Omega_\mathrm{P}$ is the ion plasma frequency. Dynamical information beyond the static structure becomes relevant when $T\lesssim 0.1~\Omega_\mathrm{P}$. We discuss the implications of these findings for calculations of $\kappa$ in multi-component systems and identify strategies for using Monte Carlo techniques in future work.
We present X-ray, ultraviolet, and optical observations of 1RXS J154439.4-112820, the most probable counterpart of the unassociated Fermi LAT source 3FGL J1544.6-1125. The optical data reveal rapid variability, which is a feature of accreting systems. The X-ray data exhibit large-amplitude flux variations in the form of fast switching (within ~10 s) between two distinct flux levels that differ by a factor of $\approx$10. The detailed optical and X-ray behavior is virtually identical to that seen in the accretion-disk-dominated states of the transitional millisecond pulsar binaries PSR J1023+0038 and XSS J12270-4859, which are also associated with $\gamma$-ray sources. Based on the available observational evidence, we conclude that 1RXS J154439.4-112820 and 3FGL J1544.6-1125 are the same object, with the X-rays arising from intermittent low-luminosity accretion onto a millisecond pulsar and the $\gamma$-rays originating from an accretion-driven outflow. 1RXS J154439.4-112820 is only the fourth $\gamma$-ray emitting low-mass X-ray binary system to be identified and is likely to sporadically undergo transformations to a non-accreting rotation-powered pulsar system.
The known members of the class of hyperluminous X-ray sources (HLXs) are few in number, yet they are of great interest as they are regarded as the likeliest intermediate-mass black hole (IMBH) candidates amongst the wider population of ultraluminous X-ray sources (ULXs). Here we report optical photometry and spectroscopy of a HLX candidate associated with the galaxy IC 4320, that reveal it is a background AGN. We discuss the implications of the exclusion of this object from the small number of well-studied HLXs, that appears to accentuate the difference in characteristics between the good IMBH candidate ESO 243-49 HLX-1 and the small handful of other HLXs.
We introduce and test several novel approaches for periodicity detection in unevenly-spaced sparse datasets. Specifically, we examine five different kinds of periodicity metrics, which are based on non-parametric measures of serial dependence of the phase-folded data. We test the metrics through simulations in which we assess their performance in various situations, including various periodic signal shapes, different numbers of data points and different signal to noise ratios. One of the periodicity metrics we introduce seems to perform significantly better than the classical ones in some settings of interest to astronomers. We suggest that this periodicity metric - the Hoeffding-test periodicity metric - should be used in addition to the traditional methods, to increase periodicity detection probability.
We discuss the potential for using neutron stars to determine bounds on the Higgs-Kretschmann coupling by looking at peculiar shifts in gamma-ray spectroscopic features. In particular, we reanalyse multiple lines observed in GRB781119 detected by two gamma-ray spectrometers, and derive an upper bound on the Higgs-Kretschmann coupling that is much more constraining than the one recently obtained from white dwarfs. This calls for targeted analyses of spectra of gamma-ray bursts from more recent observatories, dedicated searches for differential shifts on electron-positron and proton-antiproton annihilation spectra in proximity of compact sources, and signals of electron and proton cyclotron lines from the same neutron star.
We account for recent observations of the binary system at the center of the bipolar planetary nebula Henize 2-428 by the presence of one degenerate core with a low-mass main sequence companion, rather than by two degenerate objects. We argue that the variability of the He II 5412A spectral line can be accounted for by a time-varying broad absorption line from the central star on top of which there is a time-varying narrow emission line from the compact nebula. The two (almost) symmetric broad minima in the light curve are attributed to tidal distortion caused by a companion. We find problems in the recently proposed and competing explanation of two equal-mass degenerate objects that supposedly will eventually merge, possibly leading to a SN Ia. We conclude that Henize 2-428 cannot be claimed yet to support the double-degenerate scenario for Type Ia supernovae.
Potentially hazardous asteroid (185851) 2000 DP107 was the first binary near-Earth asteroid to be imaged. Radar observations in 2000 provided images at 75 m resolution that revealed the shape, orbit, and spin-up formation mechanism of the binary. The asteroid made a more favorable flyby of the Earth in 2008, yielding images at 30 m resolution. We used these data to obtain shape models for the two components and to improve the estimates of the mutual orbit, component masses, and spin periods. The primary has a sidereal spin period of 2.7745 +/- 0.0007 h and is roughly spheroidal with an equivalent diameter of 863 m +/- 5%. It has a mass of 4.656 +/- 0.56 x 10^11 kg and a density of 1381 +/- 268 kg/m^3. It exhibits an equatorial ridge similar to the (66391) 1999 KW4 primary, however the equatorial ridge in this case is not as regular and has a ~300 m diameter concavity on one side. The secondary has a sidereal spin period of 1.77 +/- 0.02 days commensurate with the orbital period. The secondary is slightly elongated and has overall dimensions of 377 x 314 x 268 m (6% uncertainties). Its mass is 0.178 +/- 0.021 x 10^11 kg and its density is 1047 +/- 230 kg/m^3. The mutual orbit has a semi-major axis of 2.659 +/- 0.08 km, an eccentricity of 0.019 +/- 0.01, and a period of 1.7556 +/- 0.0015 days. The normalized total angular momentum of this system exceeds the amount required for the expected spin-up formation mechanism. An increase of angular momentum from non-gravitational forces after binary formation is a possible explanation.
We present the analysis of photometric and spectroscopic data on eight candidates for post-nova systems. Five post-novae, V528 Aql, HS Sge, BS Sgr, GR Sgr and V999 Sgr, are successfully recovered. We furthermore identify likely candidates for the fields of V1301 Aql, V1151 Sgr and V3964 Sgr. The spectroscopic properties of the confirmed post-novae are briefly discussed. We find that two of the oldest post-novae in our sample, GR Sgr and V999 Sgr, contain an optically thick accretion disc, and thus can be suspected to have a high mass-transfer rate, contrary to what one would expect from most models. HS Sge and V528 Aql show evidence for a (comparatively) high system inclination, which makes them attractive targets for time-series observations. Finally, the presence of particularly strong He II emission together with a small eruption amplitude suggests that BS Sgr is a good candidate for an intermediate polar.
The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy primarily because of their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($\beta$) as well as Ly$\alpha$ equivalent widths (EW). The spectroscopically measured $\beta$ are, on average, redder (less negative) compared to the photometrically measured $\beta$. The positive correlation of $\beta$ with the SED-based measurement of dust extinction E(B-V) emphasizes the importance of $\beta$ as an alternative dust indicator at high redshifts. For proper comparison, we divide these SFGs into three sub-groups based on their rest-frame Ly$\alpha$ EW: SFG_N (EW<=0A), SFG_L (EW>0A), and LAEs (EW=>20A). The fraction of LAEs at these redshifts is ~10%, which is consistent with previous observations. We compared best-fit SED estimated stellar parameters of the SFG_N, SFG_L and LAE samples. For the luminosities probed here, we find statistically significant correlations for dust and star-formation rates (SFR), such that, SFG_L (and LAEs) are less dusty and low star-forming compared to SFG_N, but the differences are small compared to the large dispersion in these stellar parameters. We do not observe any significant difference in stellar mass or UV absolute magnitude. We also observe similar trends of decreasing dust and SFR with increasing Ly$\alpha$ EW. When we divide the LAEs based on their Spitzer/IRAC 3.6$\mu$m fluxes, we find that the fraction of IRAC-detected (m$_{3.6}$<~25 mag) LAEs is much higher than the fraction of IRAC-detected NB-selected LAEs at z~2-3. [abridged]
The analysis of physical measurements often copes with highly correlated noises and interruptions caused by outliers, saturation events or transmission losses. We assess the impact of missing data on the performance of linear regression analysis involving the fit of modeled or measured time series. We show that data gaps can significantly alter the precision of the regression parameter estimation in the presence of colored noise, due to the frequency leakage of the noise power. We present a regression method which cancels this effect and estimates the parameters of interest with a precision comparable to the complete data case, even if the noise power spectral density (PSD) is not known a priori. The method is based on an autoregressive (AR) fit of the noise, which allows us to build an approximate generalized least squares estimator approaching the minimal variance bound. The method, which can be applied to any similar data processing, is tested on simulated measurements of the MICROSCOPE space mission, whose goal is to test the Weak Equivalence Principle (WEP) with a precision of $10^{-15}$. In this particular context the signal of interest is the WEP violation signal expected to be found around a well defined frequency. We test our method with different gap patterns and noise of known PSD and find that the results agree with the mission requirements, decreasing the uncertainty by a factor 60 with respect to ordinary least squares methods. We show that it also provides a test of significance to assess the uncertainty of the measurement.
The ratio of total mass $M$ to surface radius $R$ of spherical perfect fluid ball has an upper bound, $M/R < B$. Buchdahl obtained $B = 4/9$ under the assumptions; non-increasing mass density in outward direction, and barotropic equation of states. Barraco and Hamity decreased the Buchdahl's bound to a lower value $B = 3/8$ $(< 4/9)$ by adding the dominant energy condition to Buchdahl's assumptions. In this paper, we further decrease the Barraco-Hamity's bound to $B \simeq 0.3636403$ $(< 3/8)$ by adding the subluminal (slower-than-light) condition of sound speed. In our analysis, we solve numerically Tolman-Oppenheimer-Volkoff equations, and the mass-to-radius ratio is maximized by variation of mass, radius and pressure inside the fluid ball as functions of mass density.
We consider a pure scalar-Gauss-Bonnet gravitational theory without the Ricci scalar. We demonstrate that such a theory, with a quadratic coupling function between the scalar field and the Gauss-Bonnet term, naturally supports inflationary -- de Sitter -- solutions. During inflation, the scalar field decays exponentially and its effective potential remains always bounded. The theory contains also solutions where these de Sitter phases possess a natural exit mechanism and are replaced by linearly expanding -- Milne -- phases.
It is shown that Starobinsky-like inflation can be realized in non-geometric flux compactifications of string theory, where the inflaton is an axion whose shift symmetry can protect UV-corrections to the scalar potential. For that purpose we evaluate the backreacted, uplifted F-term axion-monodromy potential for large field values. Limitations due to the use of an effective field theory description are pointed out.
The direction angle of twilight sky background polarization in the celestial sphere points far from solar vertical is found to depend on the ratio of single and multiple scattering contributions. The polarization direction behavior during the twilight is related with general properties of background components and can be used to check the accuracy of their separation method. The basic assumptions of this method are confirmed by analysis made in this paper. This helps to hold the temperature and scattering medium study in the upper mesosphere, the least accessible layer of the Earth's atmosphere. The mesosphere temperature data obtained during the observations from 2011 till 2014 are also presented.
A brief description of the elements of noncommutative spectral geometry as an approach to unification is presented. The physical implications of the doubling of the algebra are discussed. Some high energy phenomenological as well as various cosmological consequences are presented. A constraint in one of the three free parameters, namely the one related to the coupling constants at unification, is obtained, and the possible role of scalar fields is highlighted. A novel spectral action approach based upon zeta function regularisation, in order to address some of the issues of the traditional bosonic spectral action based on a cutoff function and a cutoff scale, is discussed.
Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing traditionally distinct scientific boundaries. Following the paper of Erlykin et al. (Astropart. Phys. 57--58 (2014) 26--29) we develop the interpretation of observed changes in long-wave (LW) radiation (Aplin and Lockwood, Env. Res. Letts. 8, 015026 (2013)), by taking account of cosmic ray ionisation yields and atmospheric radiative transfer. To demonstrate this, we show that the thermal structure of the whole atmosphere needs to be considered along with the vertical profile of ionisation. Allowing for ionisation by all components of a cosmic ray shower and not just by the muons, reveals that the effect we have detected is certainly not inconsistent with laboratory observations of the LW absorption cross section. The analysis presented here, although very different from that of Erlykin et al., does come to the same conclusion that the events detected were not caused by individual cosmic ray primaries -- not because it is impossible on energetic grounds, but because events of the required energy are too infrequent for the 12/hr rate at which they were seen by the AL experiment. The present paper numerically models the effect of three different scenario changes to the primary GCR spectrum which all reproduce the required magnitude of the observed effect. However, they cannot solely explain the observed delay in the peak effect which, if confirmed, would appear to open up a whole new and interesting area in the study of water oligomers and their effects on LW radiation. We argue that a technical artefact in the earlier experiment is highly unlikely and that our initial observations merit both follow-up experiments and more rigorous, self-consistent, three-dimensional radiative transfer modelling.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)