In this work we investigate the phenomenon of blue outliers in two complete samples of radio-emitting narrow-line Seyfert 1 galaxies (NLS1s), one radio-loud and one radio-quiet. By analysing their optical spectra and decomposing the [O III] lines, we investigate the different properties of the narrow line region (NLR) in the samples. This provides in turn useful information on the jet formation mechanism, and on its interaction with the NLR.
The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.
Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $\mu$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($\Delta I/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($\nu_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/\mu{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4\sigma$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,\mu$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.
A long-held vision has been to realize diffraction-limited optical aperture synthesis over kilometer baselines. This will enable imaging of stellar surfaces and their environments, and reveal interacting gas flows in binary systems. An opportunity is now opening up with the large telescope arrays primarily erected for measuring Cherenkov light in air induced by gamma rays. With suitable software, such telescopes could be electronically connected and also used for intensity interferometry. Second-order spatial coherence of light is obtained by cross correlating intensity fluctuations measured in different pairs of telescopes. With no optical links between them, the error budget is set by the electronic time resolution of a few nanoseconds. Corresponding light-travel distances are approximately one meter, making the method practically immune to atmospheric turbulence or optical imperfections, permitting both very long baselines and observing at short optical wavelengths. Previous theoretical modeling has shown that full images should be possible to retrieve from observations with such telescope arrays. This project aims at verifying diffraction-limited imaging experimentally with groups of detached and independent optical telescopes. In a large optics laboratory, artificial stars were observed by an array of small telescopes. Using high-speed photon-counting solid-state detectors, intensity fluctuations were cross-correlated over up to 180 baselines between pairs of telescopes, producing coherence maps across the interferometric Fourier-transform plane. These measurements were used to extract parameters about the simulated stars, and to reconstruct their two-dimensional images. As far as we are aware, these are the first diffraction-limited images obtained from an optical array only linked by electronic software, with no optical connections between the telescopes.
Radio observations at multiple frequencies have detected a significant isotropic emission component between 22 MHz and 10 GHz, commonly termed the ARCADE-2 Excess. The origin of this radio emission is unknown, as the intensity, spectrum and isotropy of the signal are difficult to model with either traditional astrophysical mechanisms or novel physics such as dark matter annihilation. We posit a new model capable of explaining the key components of the excess radio emission. Specifically, we show that the re-acceleration of non-thermal electrons via turbulent shocks in merging galaxy clusters are capable of explaining the intensity, spectrum, and isotropy of the ARCADE-2 data. We examine the parameter spaces of cluster re-acceleration, magnetic field, and merger rate, finding that the radio excess can be reproduced assuming reasonable assumptions for each. We additionally show that this model is compatible with existing observations of the Coma cluster. Finally, we point out that future observations will definitively confirm or rule-out the contribution of cluster mergers to the isotropic radio background.
The examination of two 2010 Chandra ACIS exposures of the Circinus galaxy resulted in the discovery of two pulsators: CXO J141430.1-651621 and CXOU J141332.9-651756. We also detected 26-ks pulsations in CG X-1, consistently with previous measures. For ~40 other sources, we obtained limits on periodic modulations. In CXO J141430.1-651621, which is ~2 arcmin outside the Circinus galaxy, we detected signals at 6120(1) s and 64.2(5) ks. In the longest observation, the source showed a flux of ~1.1e-13 erg/cm^2/s (absorbed, 0.5-10 keV) and the spectrum could be described by a power-law with photon index ~1.4. From archival observations, we found that the luminosity is variable by ~50 per cent on time-scales of weeks-years. The two periodicities pin down CXO J141430.1-651621 as a cataclysmic variable of the intermediate polar subtype. The period of CXOU J141332.9-651756 is 6378(3) s. It is located inside the Circinus galaxy, but the low absorption indicates a Galactic foreground object. The flux was ~5e-14 erg/cm^2/s in the Chandra observations and showed ~50 per cent variations on weekly/yearly scales; the spectrum is well fit by a power law ~0.9. These characteristics and the large modulation suggest that CXOU J141332.9-651756 is a magnetic cataclysmic variable, probably a polar. For CG X-1, we show that if the source is in the Circinus galaxy, its properties are consistent with a Wolf-Rayet plus black hole binary. We consider the implications of this for ultraluminous X-ray sources and the prospects of Advanced LIGO and Virgo. In particular, from the current sample of WR-BH systems we estimate an upper limit to the detection rate of stellar BH-BH mergers of ~16 events per yr.
The progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M Ch WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels. The ejecta structures of our simulated detonations in sub-M Ch ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (1.2 Msun) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.
We present the stellar population analysis of a sample of 12 dwarf elliptical galaxies, observed with the SAURON integral field unit, using the full-spectrum fitting method. We show that star formation histories (SFHs) resolved into two populations can be recovered even within a limited wavelength range, provided that high S/N data is used. We confirm that dEs have had complex SFHs, with star formation extending to (more) recent epochs: for the majority of our galaxies star formation activity was either still strong a few ($\lesssim$ 5) Gyr ago or they experienced a secondary burst of star formation roughly at that time. This latter possibility is in agreement with the proposed dE formation scenario where tidal harassment drives the gas remaining in their progenitors inwards and induces a star formation episode. For one of our field galaxies, ID0918, we find a correlation between its stellar population and kinematic properties, pointing to a possible merger origin of its kinematically-decoupled core. One of our cluster objects, VCC1431, appears to be composed exclusively of an old population ($\gtrsim$ 10-12 Gyr). Combining this with our earlier dynamical results, we conclude that the galaxy was either ram-pressure stripped early on in its evolution in a group environment and subsequently tidally heated, or that it evolved in situ in the cluster's central parts, compact enough to avoid tidal disruption. These are only two of the examples illustrating the SFH richness of these objects confirmed with our data.
Theoretical modeling of the redshift-space power spectrum of galaxies is crucially important to correctly extract cosmological information from redshift surveys. The task is complicated by the nonlinear biasing and redshift space distortion effects, which change with halo mass, and by the wide distribution of halo masses and their occupations by galaxies. One of the main modeling challenges is the existence of satellite galaxies that have both radial distribution and large virial velocities inside halos, a phenomenon known as the Finger-of-God effect. We present a model for the galaxy power spectrum of in which we decompose a given galaxy sample into central and satellite galaxies and relate different contributions to 1- and 2-halo terms in a halo model. Our primary goal is to ensure that any parameters that we introduce have physically meaningful values, and are not just fitting parameters. For the 2-halo terms we use the previously developed RSD modeling of halos in the context of distribution function and perturbation theory approach. This term needs to be multiplied by the effect of radial distances and velocities of satellites inside the halo. To this one needs to add the 1-halo terms, which are non-perturbative. We show that the real space 1-halo terms can be modeled as almost constant, with the finite extent of the satellites inside the halo inducing a small k^2R^2P(k) term, where R is related to the size of the halo. We adopt a similar model for FoG in redshift space, ensuring that FoG velocity dispersion is related to the halo mass. For FoG k^2 type expansions do not work and FoG resummation must be used instead. We test several damping functions to model the velocity dispersion FoG effect. Applying the formalism to mock galaxies modeled after the "CMASS" sample of the BOSS survey, we find that our predictions for the redshift-space power spectra are accurate up to k~0.4Mpc/h.
We used the Magellan adaptive optics system to image the 11 Myr substellar companion 1RXS 1609 B at the bluest wavelengths to date (z' and Ys). Comparison with synthetic spectra yields a higher temperature than previous studies of $T_\mathrm{eff}=2000\pm100\mathrm{K}$ and significant dust extinction of $A_V=4.5^{+0.5}_{-0.7}$ mag. Mass estimates based on the DUSTY tracks gives 0.012-0.015 Msun, making the companion likely a low-mass brown dwarf surrounded by a dusty disk. Our study suggests that 1RXS 1609 B is one of the 25% of Upper Scorpius low-mass members harboring disks, and it may have formed like a star and not a planet out at 320 AU.
The U.S. Naval Observatory Robotic Astrometric Telescope (URAT) is observing the northern sky since April 2012 for an astrometric survey. Multiple overlaps per year are performed in a single bandpass (680$-$750 nm) using the "redlens" 20 cm aperture astrograph and a mosaic of large CCDs. Besides the regular, deep survey to magnitude 18.5, short exposures with an objective grating are taken to access stars as bright as 3rd magnitude. A brief overview of the program, observing and reductions is given. Positions on the 8 to 20 mas level are obtained of 66,202 Hipparcos stars at current epochs. These are compared to the Hipparcos Catalog to investigate its accuracy. About 20\% of the observed Hipparcos stars are found to have inconsitent positions with the Hipparcos Catalog prediction on the 3 sigma level or over (about 75 mas or more discrepant position offsets). Some stars are now seen at an arcsec (or 25 sigma) off their Hipparcos Catalog predicted position.
We study cluster mass dark matter haloes, their progenitors and surroundings in an coupled Dark Matter-Dark Energy model and compare it to quintessence and $\Lambda$CDM models with adiabatic zoom simulations. When comparing cosmologies with different expansions histories, growth functions & power spectra, care must be taken to identify unambiguous signatures of alternative cosmologies. Shared cosmological parameters, such as $\sigma_8$, need not be the same for optimal fits to observational data. We choose to set our parameters to $\Lambda$CDM $z=0$ values. We find that in coupled models, where DM decays into DE, haloes appear remarkably similar to $\Lambda$CDM haloes despite DM experiencing an additional frictional force. Density profiles are not systematically different and the subhalo populations have similar mass, spin, and spatial distributions, although (sub)haloes are less concentrated on average in coupled cosmologies. However, given the scatter in related observables ($V_{\rm max},R_{V_{\rm max}}$), this difference is unlikely to distinguish between coupled and uncoupled DM. Observations of satellites of MW and M31 indicate a significant subpopulation reside in a plane. Coupled models do produce planar arrangements of satellites of higher statistical significance than $\Lambda$CDM models, however, in all models these planes are dynamically unstable. In general, the nonlinear dynamics within and near large haloes masks the effects of a coupled dark sector. The sole environmental signature we find is that small haloes residing in the outskirts are more deficient in baryons than their $\Lambda$CDM counterparts. The lack of a pronounced signal for a coupled dark sector strongly suggests that such a phenomena would be effectively hidden from view.
The main objective of this paper is to study the possibility of triggered star formation on the border of the HII region S233, which is formed by a B-star. Using high-resolution spectra we determine the spectral class of the ionizing star as B0.5 V and the radial velocity of the star to be -17.5(1.4) km/s. This value is consistent with the velocity of gas in a wide field across the S233 region, suggesting that the ionizing star was formed from a parent cloud belonging to the S233 region. By studying spatial-kinematic structure of the molecular cloud in the S233 region, we detected an isolated clump of gas producing CO emission red-shifted relative to the parent cloud. In the UKIDSS and WISE images, the clump of gas coincides with the infrared source containing a compact object and bright-rimmed structure. The bright-rimmed structure is perpendicular to the direction of the ionizing star. The compact source coincides in position with IRAS source 05351+3549. All these features indicate a possibility of triggering formation of a next-generation star in the S233 region. Within the framework of a theoretical one-dimensional model we conclude that the "collect-and-collapse" process is not likely to take place in the S233 region. The presence of the bright-rimmed structure and the compact infrared source suggest that the "collapse of the pre-existing clump" process is taking place.
Quasi-periodic oscillations (QPOs) are observed in the optical flux of some polars with typical periods of 1 to 3 s but none have been observed yet in X-rays where a significant part of the accreting energy is released. QPOs are expected and predicted from shock oscillations. Most of the polars have been observed by the XMM-Newton satellite. We made use of the homogeneous set of observations of the polars by XMM-Newton to search for the presence of QPOs in the (0.5-10 keV) energy range and to set significant upper limits for the brightest X-ray polars. We extracted high time-resolution X-ray light curves by taking advantage of the 0.07 sec resolution of the EPIC-PN camera. Among the 65 polars observed with XMM-Newton from 1998 to 2012, a sample of 24 sources was selected on the basis of their counting rate in the PN instrument to secure significant limits. We searched for QPOs using Fast Fourier Transform (FFT) methods and defined limits of detection using statistical tools. Among the sample surveyed, none shows QPOs at a significant level. Upper limits to the fractional flux in QPOs range from 7% to 71%. These negative results are compared to the detailed theoretical predictions of numerical simulations based on a 2D hydrodynamical code presented in Paper II. Cooling instabilities in the accretion column are expected to produce shock quasi-oscillations with a maximum amplitude reaching ~ 40% in the bremsstrahlung (0.5-10 keV) X-ray emission and ~ 20% in the optical cyclotron emission. The absence of X-ray QPOs imposes an upper limit of ~ (5-10) g.cm-2.s-1 on the specific accretion rate but this condition is found inconsistent with the value required to account for the amplitudes and frequencies of the observed optical QPOs. This contradiction outlines probable shortcomings with the shock instability model.
Magnetic cataclysmic variables are close binary systems containing a strongly magnetized white dwarf that accretes matter coming from an M-dwarf companion. High-energy radiation coming from those objects is emitted from the accretion column close to the white dwarf photosphere at the impact region. Its properties depend on the characteristics of the white dwarf and an accurate accretion column model allows the properties of the binary system to be inferred, such as the white dwarf mass, its magnetic field, and the accretion rate. We study the temporal and spectral behaviour of the accretion region and use the tools we developed to accurately connect the simulation results to the X-ray and optical astronomical observations. The radiation hydrodynamics code Hades was adapted to simulate this specific accretion phenomena. Classical approaches were used to model the radiative losses of the two main radiative processes: bremsstrahlung and cyclotron. The oscillation frequencies and amplitudes in the X-ray and optical domains are studied to compare those numerical results to observational ones. Different dimensional formulae were developed to complete the numerical evaluations. The complete characterization of the emitting region is described for the two main radiative regimes: when only the bremsstrahlung losses and when both cyclotron and bremsstrahlung losses are considered. The effect of the non-linear cooling in- stability regime on the accretion column behaviour is analysed. Variation in luminosity on short timescales (~ 1 s quasi-periodic oscillations) is an expected consequence of this specific dynamic. The importance of secondary shock instability on the quasi-periodic oscillation phenomenon is discussed. The stabilization effect of the cyclotron process is confirmed by our numerical simulations, as well as the power distribution in the various modes of oscillation.
Some short GRBs are followed by longer extended emission, lasting anywhere from ~10 to ~100 s. These short GRBs with extended emission (EE) can possess observational characteristics of both short and long GRBs (as represented by GRB 060614), and the traditional classification based on the observed duration places some of them in the long GRB class. While GRBs with EE pose a challenge to the compact binary merger scenario, they may therefore provide an important link between short and long duration events. To identify the population of GRBs with EE regardless of their initial classifications, we performed a systematic search of short GRBs with EE using all available data (up to February 2013) of both Swift/BAT and Fermi/GBM. The search identified 16 BAT and 14 GBM detected GRBs with EE, several of which are common events observed with both detectors. We investigated their spectral and temporal properties for both the spikes and the EE, and examined correlations among these parameters. Here we present the results of the systematic search as well as the properties of the identified events. Finally, their properties are also compared with short GRBs with EE observed with BATSE, identified through our previous search effort. We found several strong correlations among parameters, especially when all of the samples were combined. Based on our results, a possible progenitor scenario of two-component jet is discussed.
Star formation activities in the NGC 2024 FIR 4 region were studied by imaging centimeter continuum sources and water maser sources using several archival data sets from the Very Large Array. The continuum source VLA 9 is elongated in the northwest-southeast direction, consistent with the FIR 4 bipolar outflow axis, and has a flat spectrum in the 6.2-3.6 cm interval. The three water maser spots associated with FIR 4 are also distributed along the outflow axis. One of the spots is located close to VLA 9, and another one is close to an X-ray source. Examinations of the positions of compact objects in this region suggest that the FIR 4 cloud core contains a single low-mass protostar. VLA 9 is the best indicator of the protostellar position. VLA 9 may be a radio thermal jet driven by this protostar, and it is unlikely that FIR 4 contains a high-mass young stellar object (YSO). A methanol 6.7 GHz maser source is located close to VLA 9, at a distance of about 100 AU. The FIR 4 protostar must be responsible for the methanol maser action, which suggests that methanol class II masers are not necessarily excited by high-mass YSOs. Also discussed are properties of other centimeter continuum sources in the field of view and the water masers associated with FIR 6n. Some of the continuum sources are radio thermal jets, and some are magnetically active young stars.
High-energy neutrinos and photons are complementary messengers, probing violent astrophysical processes and structural evolution of the Universe. X-ray and neutrino observations jointly constrain conditions in active galactic nuclei (AGN) jets: their baryonic and leptonic contents, and particle production efficiency. Testing two standard neutrino production models for local source Cen A \citep{KT2008,BB2009}, we calculate the high-energy neutrino spectra of single AGN sources and derive the flux of high-energy neutrinos expected for the current epoch. Assuming that accretion determines both X-rays and particle creation, our parametric scaling relations predict neutrino yield in various AGN classes. We derive redshift-dependent number densities of each class, from {\it Chandra} and {\it Swift}/BAT X-ray luminosity functions \citep{SGB2008,ACS2009}. We integrate the neutrino spectrum expected from the cumulative history of AGN (correcting for cosmological and source effects, e.g. jet orientation and beaming). Both emission scenarios yield neutrino fluxes well above limits set by {\it IceCube} (by $\sim 4$--$10^6 \times$ at 1 PeV, depending on the assumed jet models for neutrino production). This implies that: (i) Cen A might not be a typical neutrino source as commonly assumed; (ii) both neutrino production models overestimate the efficiency; (iii) neutrino luminosity scales with accretion power differently among AGN classes and hence does not follow X-ray luminosity universally; (iv) some AGN are neutrino-quiet (e.g. below a power threshold for neutrino production); (v) neutrino and X-ray emission have different duty cycles (e.g. jets alternate between baryonic and leptonic flows); or (vi) some combination of the above.
We present an estimate of the mass of the supermassive black hole (SMBH) in the nearby type-1 Seyfert galaxy \object{NGC 1097} using Atacamma Large Millimeter/Submillimeter Array (ALMA) observations of dense gas kinematics. Dense molecular gas dynamics are traced with ${\rm HCN} (J=1-0)$ and ${\rm HCO^{+}} (J=1-0)$ emission lines. Assuming a host galaxy inclination of $46^{\circ}$, we derive a SMBH mass, $M_{\rm BH}=1.40^{+0.27}_{-0.32} \times 10^{8}M_{\odot}$, and an I-band mass to light ratio to be $5.14^{+0.03}_{-0.04}$, using ${\rm HCN} (J=1-0)$. The estimated parameters are consistent between the two emission lines. The measured SMBH mass is in good agreement with the SMBH mass and bulge velocity dispersion relationship. Our result showcases ALMA's potential for deriving accurate SMBH masses, especially for nearby late-type galaxies. Larger samples and accurate SMBH masses will further elucidate the relationship between the black hole (BH) and host galaxy properties and constrain the coevolutionary growth of galaxies and BHs.
We construct equilibrium configurations of uniformly rotating neutron stars for selected relativistic mean-field nuclear matter equations of state (EOS). We compute in particular the gravitational mass ($M$), equatorial ($R_{\rm eq}$) and polar ($R_{\rm pol}$) radii, eccentricity, angular momentum ($J$), moment of inertia ($I$) and quadrupole moment ($M_2$) of neutron stars stable against mass-shedding and secular axisymmetric instability. By constructing the constant frequency sequence $f=716$ Hz of the fastest observed pulsar, PSR J1748-2446ad, and constraining it to be within the stability region, we obtain a lower mass bound for the pulsar, $M_{\rm min}=[1.2$-$1.4] M_\odot$, for the EOS employed. Moreover we give a fitting formula relating the baryonic mass ($M_b$) and gravitational mass of non-rotating neutron stars, $M_b/M_\odot=M/M_\odot+(13/200)(M/M_\odot)^2$ [or $M/M_\odot=M_b/M_\odot-(1/20)(M_b/M_\odot)^2$], which is independent on the EOS. We also obtain a fitting formula, although not EOS independent, relating the gravitational mass and the angular momentum of neutron stars along the secular axisymmetric instability line for each EOS. We compute the maximum value of the dimensionless angular momentum, $a/M\equiv c J/(G M^2)$ (or "Kerr parameter"), $(a/M)_{\rm max}\approx 0.7$, found to be also independent on the EOS. We compare and contrast then the quadrupole moment of rotating neutron stars with the one predicted by the Kerr exterior solution for the same values of mass and angular momentum. Finally we show that, although the mass quadrupole moment of realistic neutron stars never reaches the Kerr value, the latter is closely approached from above at the maximum mass value, as physically expected from the no-hair theorem. In particular the stiffer the EOS is, the closer the Kerr solution is approached.
The Arecibo Ultra Deep Survey (AUDS) combines the unique sensitivity of the telescope with the wide field of the Arecibo L-band Feed Array (ALFA) to directly detect 21cm HI emission from galaxies at distances beyond the local Universe bounded by the lower frequency limit of ALFA (z=0.16). AUDS has collected 1110 hours of integration time in two fields with a combined area of 1.35 square degrees. In this paper we present data from 50% of the total survey, corresponding to a sensitivity level of 80 micro-Jy. We discuss the data reduction, the search for galaxies, parametrisation, optical identification and completeness. We detect 102 galaxies in the mass range of log M_HI/M_sun-2log h=5.6-10.3. We compute the HI mass function (HIMF) at the highest redshifts so far measured. A fit of a Schechter function results in alpha=-1.37+-0.03, Phi=(7.72+-1.4)*10^3 h^3/Mpc^3 and log M_HI/M_sun=9.75+-0.041+2log h. Using the measured HIMF, we find a cosmic HI density of Omega_HI=(2.33+-0.07)*10^-4/h for the sample z=0.065. We discuss further uncertainties arising from cosmic variance. Because of its depth, AUDS is the first survey that can determine parameters for the HI mass function in independent redshift bins from a single homogeneous data set. The results indicate little evolution of the co-moving mass function and Omega_HI within this redshift range. We calculate a weighted average for Omega_HI in the range $0<z<0.2$, combining the results from AUDS as well as results from other 21cm surveys and stacking, finding a best combined estimate of Omega_HI=(2.63+-0.10)*10-4/h.
This PhD thesis (defended in 2014) is focused on the estimation of the CMB polarised anisotropies power spectra on a masked sky and on forecasts of constraints set on the primordial universe physics thanks to these anisotropies. After an introduction on the light polarisation, the standard model of cosmology and the CMB properties, I show the results obtained on the use and efficiency of pseudospectrum methods to correct for the so-called E-to-B leakage. Afterwards, I present the forecasts obtained on the detection of the tensor-to-scalar ratio r and on the detection of chiral gravity, using the pure pseudospectrum method. The study of forecasts of a primordial magnetic field detection using the CMB polarised anisotropies is finally briefly tackled.
Stellar magnetism, explorable via polarimetry, is a crucial driver of activity, ionization, photodissociation, chemistry and winds in stellar environments. Thus it has an important impact on the atmospheres and magnetospheres of surrounding planets. Modeling of stellar magnetic fields and their winds is extremely challenging, both from the observational and the theoretical points of view, and only recent ground breaking advances in observational instrumentation - as were discussed during this Symposium - and a deeper theoretical understanding of magnetohydrodynamic processes in stars enable us to model stellar magnetic fields and winds and the resulting influence on surrounding planets in more and more detail. We have initiated a national and international research network (NFN): 'Pathways to Habitability - From Disks to Active Stars, Planets to Life', to address questions on the formation and habitability of environments in young, active stellar/planetary systems. In this contribution we discuss the work we are carrying out within this project and focus on how stellar magnetic fields, their winds and the relation to stellar rotation can be assessed observationally with relevant techniques such as Zeeman Doppler Imaging (ZDI), field extrapolation and wind simulations.
We present multi-color CCD photometry of the neglected contact binary XZ Leo. Completely covered VRI band light curves and four times of minimum light were obtained. Combining the photometric and previously published radial velocity data, a revised photometric analysis was carried out for the binary system by applying the Wilson-Devinney code. With a hot spot placed on the massive primary component near the neck region of the common envelope, the light curves were satisfactorily modeled. The photometric solution combined with the radial velocity solution reveals that XZ Leo is an A-type contact binary with a degree of contact of 24($\pm1)\%$. The absolute parameters of the components were determined as M_1 = 1.74($\pm$0.06)M_\odot, M_2 = 0.61($\pm$0.02)M_\odot, R_1 = 1.69($\pm$0.01)R_\odot, R_2 = 1.07($\pm0.01$)R_\odot, L_1 = 6.73($\pm0.08$) L_\odot, L_2 = 2.40($\pm$0.04)L_\odot. Based on all the available data, the long-term orbital period behavior of the system was investigated. It indicates that the binary system was undergoing continuous orbital period increase in the past three decades with a rate of dP/dt = + 6.12 \times {10^{-8}} days yr^{-1}, which suggests a probable mass transfer from the secondary to the primary component at a rate of dM/dt= 3.92\times 10^{-8} M_\odot yr^{-1}. The binary system is expected to evolve into the broken-contact stage in $1.56 \times 10^6$ years. This could be evidence supporting the Thermal Relaxation Oscillation theory.
We present constraints on neutrino masses, the primordial fluctuation
spectrum from inflation, and other parameters of the $\Lambda$CDM model, using
the one-dimensional Ly$\alpha$-forest power spectrum measured by
Palanque-Delabrouille et al. (2013) from SDSS-III/BOSS, complemented by Planck
2015 cosmic microwave background (CMB) data and other cosmological probes. This
paper improves on the previous analysis by Palanque-Delabrouille et al. (2015)
by using a more powerful set of calibrating hydrodynamical simulations that
reduces uncertainties associated with resolution and box size, by adopting a
more flexible set of nuisance parameters for describing the evolution of the
intergalactic medium, by including additional freedom to account for systematic
uncertainties, and by using Planck 2015 constraints in place of Planck 2013.
Fitting Ly$\alpha$ data alone leads to cosmological parameters in excellent
agreement with the values derived independently from CMB data, except for a
weak tension on the scalar index $n_s$. Combining BOSS Ly$\alpha$ with Planck
CMB constrains the sum of neutrino masses to $\sum m_\nu < 0.12$ eV (95\% C.L.)
including all identified systematic uncertainties, tighter than our previous
limit (0.15 eV) and more robust. Adding Ly$\alpha$ data to CMB data reduces the
uncertainties on the optical depth to reionization $\tau$, through the
correlation of $\tau$ with $\sigma_8$. Similarly, correlations between
cosmological parameters help in constraining the tensor-to-scalar ratio of
primordial fluctuations $r$. The tension on $n_s$ can be accommodated by
allowing for a running ${\mathrm d}n_s/{\mathrm d}\ln k$. Allowing running as a
free parameter in the fits does not change the limit on $\sum m_\nu$. We
discuss possible interpretations of these results in the context of slow-roll
inflation.
The aim of this paper is to test the basic model of negative AGN feedback,
according to which once the central black hole accretes at the Eddington limit
and reaches a certain critical mass, AGN driven outflows blow out gas,
suppressing star formation in the host galaxy and self-regulating black hole
growth.
We consider a sample of 224 quasars selected from the Sloan Digital Sky
Survey(SDSS) at z < 1 observed in the infrared band by the Herschel Space
Observatory in point source photometry mode. We evaluate the star formation
rate in relation to several outflow signatures traced by [O III]4959,5007 and
[O III]3726,3729 emission line in about half of the sample with high quality
spectra. Most of the quasars show asymmetric and broad wings in [O III], clear
outflow signatures. We separate the quasars in two groups: "weakly" and
"strongly" outflowing using three different criteria. When we compare the mean
star formation rate in 5 bins of redshift in the two groups, we find that the
SFRs are comparable or slightly larger in the strongly outflowing quasars, in
contrast with what it is predicted by the basic negative AGN feedback model.
Moreover, for quasars dominated in the infrared by starburst or by AGN
emission, we do not find any correlation between the star formation rate and
the velocity of the outflow, a trend previously reported in literature for pure
starburst galaxies.
We conclude that the basic AGN negative feedback scenario seems not to be in
accordance with our results. Although we use a large sample of quasars, we did
not find any evidence that the star formation rate is suppressed in presence of
AGN driven outflows on large scale. A possibility is that feedback is effective
over much longer timescales than those of sigle episode of quasar activity.
Lithium abundance derived in metal-poor main sequence stars is about three times lower than the value of primordial Li predicted by the standard Big Bang nucleosynthesis when the baryon density is taken from the CMB or the deuterium measurements. This disagreement is generally referred as the lithium problem. We here reconsider the stellar Li evolution from the pre-main sequence to the end of the main sequence phase by introducing the effects of convective overshooting and residual mass accretion. We show that $^7$Li could be significantly depleted by convective overshooting in the pre-main sequence phase and then partially restored in the stellar atmosphere by a tail of matter accretion which follows the Li depletion phase and that could be regulated by EUV photo-evaporation. By considering the conventional nuclear burning and microscopic diffusion along the main sequence we can reproduce the Spite plateau for stars with initial mass $m_0=0.62 - 0.80 M_{\odot}$, and the Li declining branch for lower mass dwarfs, e.g, $m_0=0.57 - 0.60 M_{\odot}$, for a wide range of metallicities (Z=0.00001 to Z=0.0005), starting from an initial Li abundance $A({\rm Li}) =2.72$. This environmental Li evolution model also offers the possibility to interpret the decrease of Li abundance in extremely metal-poor stars, the Li disparities in spectroscopic binaries and the low Li abundance in planet hosting stars.
We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of -5/3 and -3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ~-3. A -2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.
We report on the results of the first XMM-Newton satellite observation of the luminous and helium-rich O-type subdwarf BD +37{\deg} 1977 carried out in April 2014. X-ray emission is detected with a flux of about 4*10^(-14) erg/cm2/s (0.2-1.5 keV), corresponding to a f_X/f_bol ratio about 10^(-7); the source spectrum is very soft, and is well fit by the sum of two plasma components at different temperatures. Both characteristics are in agreement with what is observed in the main-sequence early-type stars, where the observed X-ray emission is due to turbulence and shocks in the stellar wind. A smaller but still significant stellar wind has been observed also in BD +37{\deg} 1977; therefore, we suggest that also in this case the detected X-ray flux has the same origin.
The installation and operation of a telescope in Antarctica represent particular challenges, in particular the requirement to operate at extremely cold temperatures, to cope with rapid temperature fluctuations and to prevent frosting. Heating of electronic subsystems is a necessity, but solutions must be found to avoid the turbulence induced by temperature fluctua- tions on the optical paths. ASTEP 400 is a 40 cm Newton telescope installed at the Concordia station, Dome C since 2010 for photometric observations of fields of stars and their exoplanets. While the telescope is designed to spread star light on several pixels to maximize photometric stability, we show that it is nonetheless sensitive to the extreme variations of the seeing at the ground level (between about 0.1 and 5 arcsec) and to temperature fluctuations between --30 degrees C and --80 degrees C. We analyze both day-time and night-time observations and obtain the magnitude of the seeing caused by the mirrors, dome and camera. The most important effect arises from the heating of the primary mirror which gives rise to a mirror seeing of 0.23 arcsec K--1 . We propose solutions to mitigate these effects.
The Kuiper belt is a population of icy bodies beyond the orbit of Neptune. A particularly puzzling and up-to-now unexplained feature of the Kuiper belt is the so-called `kernel', a concentration of orbits with semimajor axes a~44 AU, eccentricities e~0.05, and inclinations i<5 deg. Here we show that the Kuiper belt kernel can be explained if Neptune's otherwise smooth migration was interrupted by a discontinuous change of Neptune's semimajor axis when Neptune reached ~28 AU. Before the discontinuity happened, planetesimals located at ~40 AU were swept into Neptune's 2:1 resonance, and were carried with the migrating resonance outwards. The 2:1 resonance was at ~44 AU when Neptune reached ~28 AU. If Neptune's semimajor axis changed by fraction of AU at this point, perhaps because Neptune was scattered off of another planet, the 2:1 population would have been released at ~44 AU, and would remain there to this day. We show that the orbital distribution of bodies produced in this model provides a good match to the orbital properties of the kernel. If Neptune migration was conveniently slow after the jump, the sweeping 2:1 resonance would deplete the population of bodies at ~45-47 AU, thus contributing to the paucity of the low-inclination orbits in this region. Special provisions, probably related to inefficiencies in the accretional growth of sizable objects, are still needed to explain why only a few low-inclination bodies have been so far detected beyond ~47 AU.
In the letter by Stadnik and Flambaum [Phys. Rev. Lett. 113, 151301 (2014)] it is claimed that topological defects passing through pulsars could be responsible for the observed pulsar glitches. Here, we show that, independently of the detailed network dynamics and defect dimensionality, defect networks cannot be at the origin of the pulsar glitch phenomenon.
Optical (BVRI) and infrared (JHK) photometry of the gamma-ray nova Nova Sco 2012 (V1324 Sco) is presented and the lightcurve reconstructed and discussed. An interstellar reddening E(B-V)=1.23 is derived. Dust begun to form at an early date in the nova, only one magnitude down and 20 days past maximum optical brightness and caused an extinction of at least 6 magnitudes in V band, that cleared some months later. This unusual early dust formation compromises the application of the magnitude at maximum versus rate of decline (MMRD) relations in estimating the distance to the nova.
The First Geiger-mode Avalanche photodiode (G-APD) Cherenkov Telescope (FACT)
has been operating on the Canary island of La Palma since October 2011.
Operations were automated so that the system can be operated remotely. Manual
interaction is required only when the observation schedule is modified due to
weather conditions or in case of unexpected events such as a mechanical
failure. Automatic operations enabled high data taking efficiency, which
resulted in up to two terabytes of FITS files being recorded nightly and
transferred from La Palma to the FACT archive at ISDC in Switzerland. Since
long term storage of hundreds of terabytes of observations data is costly, data
compression is mandatory. This paper discusses the design choices that were
made to increase the compression ratio and speed of writing of the data with
respect to existing compression algorithms.
Following a more detailed motivation, the FACT compression algorithm along
with the associated I/O layer is discussed. Eventually, the performances of the
algorithm is compared to other approaches.
In this paper we study the soft X-ray (SXR) signatures of one particular prominence. The X-ray observations used here were made by the Hinode/XRT instrument using two different filters. Both of them have a pronounced peak of the response function around 10 A. One of them has a secondary smaller peak around 170 A, which leads to a contamination of SXR images. The observed darkening in both of these filters has a very large vertical extension. The position and shape of the darkening corresponds nicely with the prominence structure seen in SDO/AIA images. First we have investigated the possibility that the darkening is caused by X-ray absorption. But detailed calculations of the optical thickness in this spectral range show clearly that this effect is completely negligible. Therefore the alternative is the presence of an extended region with a large emissivity deficit which can be caused by the presence of cool prominence plasmas within otherwise hot corona. To reproduce the observed darkening one needs a very large extension along the line-of-sight of the region amounting to around 10$^5$ km. We interpret this region as the prominence spine, which is also consistent with SDO/AIA observations in EUV.
Large Galactic stellar surveys and new generations of stellar atmosphere models and spectral line formation computations need to be subjected to careful calibration and validation and to benchmark tests. We focus on cool stars and aim at establishing a sample of 34 Gaia FGK Benchmark Stars with a range of different metallicities. The goal was to determine the effective temperature and the surface gravity independently from spectroscopy and atmospheric models as far as possible. Fundamental determinations of Teff and logg were obtained in a systematic way from a compilation of angular diameter measurements and bolometric fluxes, and from a homogeneous mass determination based on stellar evolution models. The derived parameters were compared to recent spectroscopic and photometric determinations and to gravity estimates based on seismic data. Most of the adopted diameter measurements have formal uncertainties around 1%, which translate into uncertainties in effective temperature of 0.5%. The measurements of bolometric flux seem to be accurate to 5% or better, which contributes about 1% or less to the uncertainties in effective temperature. The comparisons of parameter determinations with the literature show in general good agreements with a few exceptions, most notably for the coolest stars and for metal-poor stars. The sample consists of 29 FGK-type stars and 5 M giants. Among the FGK stars, 21 have reliable parameters suitable for testing, validation, or calibration purposes. For four stars, future adjustments of the fundamental Teff are required, and for five stars the logg determination needs to be improved. Future extensions of the sample of Gaia FGK Benchmark Stars are required to fill gaps in parameter space, and we include a list of suggested candidates.
PSR J1846-0258 is an object which straddles the boundary between magnetars and rotation powered pulsars. Though behaving for many years as a rotation-powered pulsar, in 2006, it exhibited distinctly magnetar-like behavior - emitting several short hard X-ray bursts, and a flux increase. Here we report on 7 years of post-outburst timing observations of PSR J1846-0258 using the Rossi X-ray Timing Explorer and the Swift X-ray Telescope. We measure the braking index over the post-magnetar outburst period to be $n=2.19\pm0.03$. This represents a change of $\Delta n=-0.46\pm0.03$ or a 14.5$\;\sigma$ difference from the pre-outburst braking index of $n=2.65\pm0.01$, which itself was measured over a span of 6.5 yr (Livingstone et al. 2006). So large and long-lived a change to a pulsar braking index is unprecedented and poses a significant challenge to models of pulsar spin-down.
We examine an interacting dark matter--variable vacuum energy model for a spatially flat Friedmann-Roberston-Walker spacetime, focusing on the appearance of cosmological singularities such as \emph{big rip, big brake, big freeze}, and \emph{ big separation} along with abrupt events (\emph{infinite $\gamma$- singularity} and \emph{new w-singularity}) at late times. We introduce a phenomenological interaction which has a nonlinear dependence on the total energy density of the dark sector and its derivative, solve exactly the source equation for the model and find the energy density as function of the scale factor as well as the time dependence of the approximate scale factor in the neighborhood of the singularities. We describe the main characteristics of these singularities by exploring the type of interaction that makes them possible along with behavior of dark components near them. We apply the geometric Tipler and Kr\'olak method for determining the fate of time-like geodesic curves around the singularities. We also explore the strength of them by analyzing the leading term in some geometric invariants such as the square Riemann scalar and the Ricci scalar.
Experimental studies of the Quark-Gluon Plasma (QGP) focus on two, in practice distinct, regimes: one in which the baryonic chemical potential $\mu_B$ is essentially zero, the other in which it is of the same order of magnitude as the temperature. The cosmic QGP which dominates the early Universe after reheating is normally assumed to be of the first kind, but recently it has been suggested that it might well be of the second: this is the case in the theory of "Little Inflation." If that is so, then it becomes a pressing issue to fix the trajectory of the Universe, as it cools, through the quark matter phase diagram: in particular, one wishes to know where in that diagram the cosmic plasma hadronizes, so that the initial conditions of the hadronic epoch can be determined. Here we combine various tools from strongly coupled QGP theory (the latest lattice results, together with gauge-gravity duality) in order to determine that trajectory, assuming that Little Inflation did occur. This can be done with surprising precision.
Reliable low-latency gravitational wave parameter estimation is essential to target limited electromagnetic followup facilities toward astrophysically interesting and electromagnetically relevant sources of gravitational waves. In this study, we examine the tradeoff between speed and accuracy. Specifically, we estimate the astrophysical relevance of systematic errors in the posterior parameter distributions derived using a fast-but-approximate waveform model, SpinTaylorF2 (STF2), in parameter estimation with lalinference_mcmc. Though efficient, the STF2 approximation to compact binary inspiral employs approximate kinematics (e.g., a single spin) and an approximate waveform (e.g., frequency domain versus time domain). More broadly, using a large astrophysically-motivated population of generic compact binary merger signals, we report on the effectualness and limitations of this single-spin approximation as a method to infer parameters of generic compact binary sources. For most low-mass compact binary sources, we find that the STF2 approximation estimates compact binary parameters with biases comparable to systematic uncertainties in the waveform. We illustrate by example the effect these systematic errors have on posterior probabilities most relevant to low-latency electromagnetic followup: whether the secondary is has a mass consistent with a neutron star; whether the masses, spins, and orbit are consistent with that neutron star's tidal disruption; and whether the binary's angular momentum axis is oriented along the line of sight.
The multipole moments and the tidal Love numbers of neutron stars and quark stars satisfy certain relations which are almost insensitive to the star's internal structure. A natural question is whether the same relations hold for different compact objects and how they possibly approach the black-hole limit. Here we consider "gravastars," which are hypothetical compact objects sustained by their internal vacuum energy. Such solutions have been proposed as exotic alternatives to the black-hole paradigm because they can be as compact as black holes and exist in any mass range. By constructing slowly-rotating, thin-shell gravastars to quadratic order in the spin, we compute the moment of inertia $I$, the mass quadrupole moment $Q$, and the tidal Love number $\lambda$ in exact form. When suitably normalized, these quantities are nonanalytical functions of the compactness of the object. The $I$-$\lambda$-$Q$ relations of a gravastar are dramatically different from those of an ordinary compact star, but the black-hole limit is continuous, i.e. these quantities approach their Kerr counterparts when the compactness is maximum. Therefore, such relations can be used to discern a gravastar from an ordinary compact star, but not to break the degeneracy with the black-hole case. Based on these results, we conjecture that the full multipolar structure and the tidal deformability of a spinning, ultracompact gravastar are identical to those of a Kerr black hole. The approach to the black-hole limit is nonanalytical, thus violating the critical behavior recently found for strongly anisotropic neutron stars.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
Kaib & Cowan (2015) recently used terrestrial planet formation simulations to conclude that the moon-forming impactor (Theia) had only a ~5% or less chance of having the same oxygen isotope composition as Earth, while Mastrobuono-Battisti et al. (2015) used seemingly similar simulations and methods to arrive at a higher value of ~20% or more. Here we derive the results of both papers from a single set of simulations. Compared to Kaib & Cowan (2015), the analysis of Mastrobuono-Battisti et al. (2015) systematically yields more massive Theia analogs and imposes flatter isotopic gradients across the original protoplanetary disk. Both of these effects diminish isotopic differences between Earth and Theia analogs. While it is notoriously difficult to produce systems resembling our actual terrestrial planets, the analysis of Kaib & Cowan (2015) more often selects and analyzes Earth and Mars analogs at orbital locations near the real planets. Given this, we conclude that the greater isotopic differences between Earth and Theia found in Kaib & Cowan (2015) better reflect the predictions of terrestrial planet formation models. Finally, although simulation uncertainties and a terrestrial contribution to Moon formation enhance the fraction of Theia analogs consistent with the canonical giant impact hypothesis, this fraction still remains in the 5-8% range.
We present significant evidence of halo assembly bias for redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$\sigma$ in the weak lensing and 4.6$\sigma$ in the auto-correlation data, respectively.
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high-order Discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall back to first order can be avoided. Finally, DG requires no deep stencils at high order, and offers an improved compute to memory access ratio compared with FV schemes, which is favorable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.
Circinus X-1 exhibited a bright X-ray flare in late 2013. Follow-up observations with Chandra and XMM-Newton from 40 to 80 days after the flare reveal a bright X-ray light echo in the form of four well-defined rings with radii from 5 to 13 arcminutes, growing in radius with time. The large fluence of the flare and the large column density of interstellar dust towards Circinus X-1 make this the largest and brightest set of rings from an X-ray light echo observed to date. By deconvolving the radial intensity profile of the echo with the MAXI X-ray lightcurve of the flare we reconstruct the dust distribution towards Circinus X-1 into four distinct dust concentrations. By comparing the peak in scattering intensity with the peak intensity in CO maps of molecular clouds from the Mopra Southern Galactic Plane CO Survey we identify the two innermost rings with clouds at radial velocity ~ -74 km/s and ~ -81 km/s, respectively. We identify a prominent band of foreground photoelectric absorption with a lane of CO gas at ~ -32 km/s. From the association of the rings with individual CO clouds we determine the kinematic distance to Circinus X-1 to be $D_{Cir X-1} = 9.4^{+0.8}_{-1.0}$ kpc. This distance rules out earlier claims of a distance around 4 kpc, implies that Circinus X-1 is a frequent super-Eddington source, and places a lower limit of $\Gamma \gtrsim 22$ on the Lorentz factor and an upper limit of $\theta_{jet} \lesssim 3^{\circ}$ on the jet viewing angle.
The Antlia Cluster is a nearby, dynamically young structure, and its proximity provides a valuable opportunity for detailed study of galaxy and group accretion onto clusters. We present a deep HI mosaic completed as part of spectral line commissioning of the Karoo Array Telescope (KAT-7), and identify infrared counterparts from the WISE extended source catalog to study neutral atomic gas content and star formation within the cluster. We detect 37 cluster members out to a radius of ~0.9 Mpc with M_HI > 5x10^7 M_Sun. Of these, 35 are new HI detections, 27 do not have previous spectroscopic redshift measurements, and one is the Compton thick Seyfert II, NGC 3281, which we detect in HI absorption. The HI galaxies lie beyond the X-ray emitting region 200 kpc from the cluster center and have experienced ram pressure stripping out to at least 600 kpc. At larger radii, they are distributed asymmetrically suggesting accretion from surrounding filaments. Combining HI with optical redshifts, we perform a detailed dynamical analysis of the internal substructure, identify large infalling groups, and present the first compilation of the large scale distribution of HI, and star forming galaxies within the cluster. We find that elliptical galaxy NGC 3268 is at the center of the oldest substructure and argue that NGC 3258 and its companion population are more recent arrivals. Through the presence of HI and on-going star formation, we rank substructures with respect to their relative time since accretion onto Antlia.
We characterise the radial density, metallicity and flattening profile of the Milky Way's stellar halo, based on the large sample of 1757 spectroscopically confirmed giant stars from SDSS/SEGUE-2 after excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream). Compared to BHB stars or RR Lyrae, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modelling of those data, based on ellipsoidal stellar density models, $\nu_* (R,z)$, with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of $n=3.1\pm 0.5$ with an effective radius $\rm r_{eff} = 15\pm2~$kpc and a flattening of $q=0.7\pm 0.02$; or comparably well by an equally flattened broken power-law, with radial slopes of $\alpha_{in}=2.1\pm 0.3$ and $\alpha_{out}=3.8\pm 0.1$, with a break-radius of $r_{break}=18\pm1$~kpc; this is largely consistent with earlier work. We find a modest, but significant metallicity gradient within the `outer' stellar halo, $\rm [Fe/H]$ decreasing outward. If we allow for a variable flattening $q = f(r_{GC} )$, we find the distribution of halo giants to be considerably more flattened at small radii, $q({\rm 10~kpc})\sim 0.57$, compared to $q(>30{\rm kpc})\sim 0.8$. Remarkably, the data are then very well fit by a single power-law of index $\rm \sim 4.2\pm0.1$ of the variable $r_q\equiv\sqrt{R^2+(z/q(r))^2}$. In this simple and better fitting model, there is a break in flattening at $\sim 20$~kpc, instead of a break in the radial density function.
We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements, wherein the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be most sensitive to the cosmological HI signal. In our recent paper, many signatures from the simulation which predicted this feature were validated against Murchison Widefield Array (MWA) data but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through coherent averaging of 12 independent snapshots with identical instrument settings, and provide the first confirmation of the prediction with a signal-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA, and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.
We present the first measurements of the abundances of alpha-elements (Mg, Si, and S) extending out to beyond the virial radius of a cluster of galaxies. Our results, based on Suzaku Key Project observations of the Virgo Cluster, show that the chemical composition of the intra-cluster medium is constant on large scales, with a flat distribution of the Si/Fe, S/Fe, and Mg/Fe ratios as a function of radius and azimuth out to 1.4 Mpc (1.3 r200). Chemical enrichment of the intergalactic medium due solely to core collapse supernovae (SNcc) is excluded with very high significance; instead, the measured S/Fe and Mg/Fe ratios are consistent with the Solar value, with a sub-solar Si/Fe ratio. The uniform metal abundance ratios observed today are likely the result of an early phase of enrichment and mixing, with both SNcc and type Ia supernovae (SNIa) contributing to the metal budget during the period of peak star formation activity at redshifts of 2-3. We estimate the ratio between the number of SNIa and the total number of supernovae enriching the inter galactic medium to be between 15-20%, generally consistent with the metal abundance patterns in our own Galaxy and only marginally lower than the SNIa contribution estimated for the cluster cores.
We present results from a 244 ks $NuSTAR$ observation of 3C 273, obtained during a cross-calibration campaign with the $Chandra$, $INTEGRAL$, $Suzaku$, $Swift$, and $XMM-Newton$ observatories. We show that the spectrum, when fit with a power-law model using data from all observatories except $INTEGRAL$ over the 1-78 keV band, leaves significant residuals in the $NuSTAR$ data between 30-78 keV. The $NuSTAR$ 3-78 keV spectrum is well-described by an exponentially cutoff power-law ($\Gamma = 1.646 \pm 0.006$, E$_\mathrm{cutoff} = 202_{-34}^{+51}$ keV) with a weak reflection component from cold, dense material. There is also evidence for a weak ($EW = 23 \pm 11$ eV) neutral iron line. We interpret these features as arising from coronal emission plus reflection off an accretion disk or distant material. Beyond 80 keV $INTEGRAL$ data show clear excess flux relative to an extrapolation of the AGN model fit to $NuSTAR$. This high-energy power-law is consistent with the presence of a beamed jet, which begins to dominate over emission from the inner accretion flow at 30-40 keV. Modeling the jet as a power-law, we find the coronal component is fit by $\Gamma_\mathrm{AGN} = 1.638 \pm 0.045$, $E_\mathrm{cutfoff} = 47 \pm 15$ keV, and jet photon index by $\Gamma_\mathrm{jet} = 1.05 \pm 0.4$. We applied comptonizing coronal electron plasma models to place constraints on the plasma temperature and optical depth. Finally, we investigate the variability and find an inverse correlation between flux and $\Gamma$. We interpret this hardening of the spectrum with increasing flux to be due to the jet.
Bolton and Schlegel presented a promising deconvolution method to extract 1D spectra from a 2D optical fiber spectral CCD image. The method could eliminate the PSF difference between fibers, extract spectra to the photo noise level, as well as improve the resolution. But the method is limited by its huge computation requirement and thus cannot be implemented in actual data reduction. In this article, we develop a practical computation method to solve the computation problem. The new computation method can deconvolve a 2D fiber spectral image of any size with actual PSFs, which may vary with positions. Our method does not require large amounts of memory and can extract a 4k multi 4k noise-free CCD image with 250 fibers in 2 hr. To make our method more practical, we further consider the influence of noise, which is thought to be an intrinsic illposed problem in deconvolution algorithms. We modify our method with a Tikhonov regularization item to depress the method induced noise. Compared with the results of traditional extraction methods, our method has the least residual and influence by cross talk and noise, even for extreme situation. Our method can convergent in 2 to 4 iterations, and the computation times are about 3.5 hr for the extreme fiber distance and about 2 hr for nonextreme cases. Finally, we apply our method to real LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope, a.k.a. Guo Shou Jing Telescope) data. We find that the 1D spectra extracted by our method have both higher signal-to-noise ratio and resolution than the traditional methods.
We present multi-epoch spectroscopy of "heartbeat stars," eccentric binaries with dynamic tidal distortions and tidally induced pulsations originally discovered with the Kepler satellite. Optical spectra of six known heartbeat stars using the Wyoming Infrared Observatory 2.3 m telescope allow measurement of stellar effective temperatures and radial velocities from which we determine orbital parameters including the periods, eccentricities, approximate mass ratios, and component masses. These spectroscopic solutions confirm that the stars are members of eccentric binary systems with eccentricities e>0.34 and periods P=7-20 days, strengthening conclusions from prior works which utilized purely photometric methods. Heartbeat stars in this sample have A- or F-type primary components. Constraints on orbital inclinations indicate that four of the six systems have minimum mass ratios q=0.3-0.5, implying that most secondaries are probable M dwarfs or earlier. One system is an eclipsing, double-lined spectroscopic binary with roughly equal-mass mid-A components (q=0.95), while another shows double-lined behavior only near periastron, indicating that the F0V primary has a G1V secondary (q=0.65). This work constitutes the first measurements of the masses of secondaries in a statistical sample of heartbeat stars. The good agreement between our spectroscopic orbital elements and those derived using a photometric model support the idea that photometric data are sufficient to derive reliable orbital parameters for heartbeat stars.
The Wilkinson Microwave Anisotropy Probe (WMAP) team subtracted template-based foreground models to produce foreground-reduced maps, and masked point sources and uncertain sky regions directly; however, whether foreground residuals exist in the WMAP foreground-reduced maps is still an open question. Here, we use Pearson correlation coefficient (PCC) analysis with AS-$\gamma$ TeV cosmic ray (CR) data to probe possible foreground residuals in the WMAP nine-year data. The correlation results between the CR and foreground-contained maps (WMAP foreground-unreduced maps, WMAP template-based and MEM foreground models) suggest that: (1) CRs can trace foregrounds in the WMAP data; (2) at least some TeV CRs originate from the Milky Way; (3) foregrounds may be related to the existence of CR anisotropy (loss-cone and tail-in structures); (4) there exist differences among different types of foregrounds in the declination range of $< 15^{\circ}$. Then, we generate 10,000 mock CMB sky maps to describe the cosmic variance, which is used to measure the effect of the fluctuations of all possible CMB maps to the correlations between CR and CMB maps. Finally, we do correlation analysis between the CR and WMAP foreground-reduced maps, and find that: (1) there are significant anticorrelations; and (2) the WMAP foreground-reduced maps are credible. However, the significant anticorrelations may be accidental, and the higher signal-to-noise ratio (SNR) Planck SMICA map cannot reject the hypothesis of accidental correlations. We therefore can only conclude that the foreground residuals exist with $\sim$ 95\% probability.
Random errors in any data set are expected to follow the Gaussian distribution with zero mean. We propose an elegant method based on Kolmogorov-Smirnov statistic to test the above and apply it on the measurement of Hubble constant which determines the expansion rate of the Universe. The measurements were made using Hubble Space Telescope. Our analysis shows that the errors in the above measurement are non-Gaussian.
Neutrinos propagate in astrophysical and cosmological environments modifying their flavor in intriguing ways. The study of neutrino propagation in media is based on the mean-field, extended mean-field and Boltzmann equations. We summarise salient features of these evolution equations and the methods employed so far to derive them. We emphasize applications to situations of observational interest.
In this paper we study periodic variability of the magnetic field in the Ap star 33 Lib. We found that its most probable period equals 83.5 years. There exist also possible shorter periods: 11.036 days, 7.649 days and 4.690 days. Analysis of the magnetic behavior of 33 Lib allows us to conclude, that the star shows the second longest period the slow rotator gamma Equ, the latter star with the period of 97 years.
We determined rotational period of practically nonrotating Ap star gamma Equ and claim $P_{\rm rot}= 97$ years. This period is about 35000 times larger than the average rotational period among stars of the same spectral class. Paper discusses possible mechanism explaining the origin of this phenomenon.
During moonlit nights, observations with ground-based Cherenkov telescopes at very high energies (VHE, $E>100$ GeV) are constrained since the photomultiplier tubes (PMTs) in the telescope camera are extremely sensitive to the background moonlight. Observations with the VERITAS telescopes in the standard configuration are performed only with a moon illumination less than 35$\%$ of full moon. Since 2012, the VERITAS collaboration has implemented a new observing mode under bright moonlight, by either reducing the voltage applied to the PMTs (reduced-high-voltage configuration, RHV), or by utilizing UV-transparent filters. While these operating modes result in lower sensitivity and increased energy thresholds, the extension of the available observing time is useful for monitoring variable sources such as blazars and sources requiring spectral measurements at the highest energies. In this paper we report the detection of $\gamma$-ray flaring activity from the BL Lac object 1ES 1727+502 during RHV observations. This detection represents the first evidence of VHE variability from this blazar. The integral flux is $(1.1\pm0.2)\times10^{-11}\mathrm{cm^{-2}s^{-1}}$ above 250 GeV, which is about five times higher than the low-flux state. The detection triggered additional \veritas\ observations during standard dark-time. Multiwavelength observations with the FLWO 48" telescope, and the Swift and Fermi satellites are presented and used to produce the first spectral energy distribution (SED) of this object during $\gamma$-ray flaring activity. The SED is then fitted with a standard synchrotron-self-Compton model, placing constraints on the properties of the emitting region and of the acceleration mechanism at the origin of the relativistic particle population in the jet.
By using the current photometric rotational data on eight galactic open clusters, we show that the evolutionary stellar model (isochrone) ages of these clusters are tightly correlated with the period shifts applied to the (B-V)_0 - P_rot ridges that optimally align these ridges to the one defined by Praesepe and the Hyades. On the other hand, when the traditional Skumanich-type multiplicative transformation is used, the ridges become far less aligned due to the age-dependent slope change introduced by the period multiplication. Therefore, we employ our simple additive gyro-age calibration on various datasets of Galactic field stars to test its applicability. We show that, in the overall sense, the gyro-ages are systematically greater than the isochrone ages. The difference could exceed several giga years, depending on the stellar parameters. Although the age overlap between the open clusters used in the calibration and the field star samples is only partial, the systematic difference indicates the limitation of the currently available gyro-age methods and suggests that the rotation of field stars slows down with a considerably lower speed than we would expect from the simple extrapolation of the stellar rotation rates in open clusters.
We investigate interaction effects on the dynamics and morphology of the galaxy pairs AM\,2058-381 and AM\,1228-260. This work is based on $r'$ images and long-slit spectra obtained with the Gemini Multi-Object Spectrograph at the Gemini South Telescope. The luminosity ratio between the main (AM\,2058A) and secondary (AM\,2058B) components of the first pair is a factor of $\sim$ 5, while for the other pair, the main (AM\,1228A) component is 20 times more luminous than the secondary (AM\,1228B). The four galaxies have pseudo-bulges, with a S\'ersic index $n<2$. Their observed radial velocities profiles (RVPs) present several irregularities. The receding side of the RVP of AM\,2058A is displaced with respect to the velocity field model, while there is a strong evidence that AM\,2058B is a tumbling body, rotating along its major axis. The RVPs for AM\,1228A indicate a misalignment between the kinematic and photometric major axes. The RVP for AM\,1228B is quite perturbed, very likely due to the interaction with AM\,1228A. NFW halo parameters for AM\,2058A are similar to those of the Milky Way and M\,31. The halo mass of AM\,1228A is roughly 10\% that of AM\,2058A. The mass-to-light (M/L) of AM\,2058 agrees with the mean value derived for late-type spirals, while the low M/L for AM\,1228A may be due to the intense star formation ongoing in this galaxy.
The ability to characterize exoplanets by spectroscopy of their atmospheres
requires direct imaging techniques to isolate planet signal from the bright
stellar glare. One of the limitations with the direct detection of exoplanets,
either with ground- or space-based coronagraphs, is pointing errors and other
low-order wavefront aberrations. The coronagraphic detection sensitivity at the
diffraction limit therefore depends on how well low-order aberrations upstream
of the focal plane mask are corrected. To prevent starlight leakage at the
inner working angle of a phase mask coronagraph, we have introduced a
Lyot-based low-order wavefront sensor (LLOWFS), which senses aberrations using
the rejected starlight diffracted at the Lyot plane. In this paper, we present
the implementation, testing and results of LLOWFS on the Subaru Coronagraphic
Extreme Adaptive Optics system (SCExAO) at the Subaru Telescope.
We have controlled thirty-five Zernike modes of a H-band vector vortex
coronagraph in the laboratory and ten Zernike modes on sky with an integrator
control law. We demonstrated a closed-loop pointing residual of 0.02 mas in the
laboratory and 0.15 mas on sky for data sampled using the minimal 2-second
exposure time of the science camera. We have also integrated the LLOWFS in the
visible high-order control loop of SCExAO, which in closed-loop operation has
validated the correction of the non-common path pointing errors between the
infrared science channel and the visible wavefront sensing channel with
pointing residual of 0.23 mas on sky.
Helium is a pivotal element in understanding multiple main sequences and extended horizontal branches observed in some globular clusters. Here we present a spectroscopic study of helium in the nearby globular cluster M4. We have obtained spectra of the chromospheric He I 10830 A line in 16 red horizontal branch, red giant branch, and asymptotic giant branch stars. Clear He I absorption or emission is present in most of the stars. Effective temperature is the principal parameter that correlates with 10830 A line strength. Stars with T_eff < 4450 K do not exhibit the helium line. Red horizontal branch stars, which are the hottest stars in our sample, all have strong He I line absorption. A number of these stars show very broad 10830 A lines with shortward extensions indicating outflows as high as 80-100 km/s and the possibility of mass loss. We have also derived [Na/Fe] and [Al/Fe] abundances to see whether these standard tracers of "second generation" cluster stars are correlated with He I line strength. Unlike the case for our previous study of Omega Cen, no clear correlation is observed. This may be because the sample does not cover the full range of abundance variations found in M4, or simply because the physical conditions in the chromosphere, rather than the helium abundance, primarily determine the He I 10830 A line strength. A larger sample of high-quality He I spectra of both "first" and "second" generation red giants within a narrow range of T_eff and luminosity is needed to test for the subtle spectroscopic variations in He I expected in M4.
We use the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey to place constraints on $f(R)$ and {\it Generalized Dilaton} models of modified gravity. This is highly complimentary to other probes since the constraints mainly come from the non-linear scales: maximal deviations with respects to the General-Relativity + $\Lambda$CDM scenario occurs at $k\sim1 h \mbox{Mpc}^{-1}$. At these scales, it becomes necessary to account for known degeneracies with baryon feedback and massive neutrinos, hence we place constraints jointly on these three physical effects. To achieve this, we formulate these modified gravity theories within a common tomographic parameterization, we compute their impact on the clustering properties relative to a GR universe, and propagate the observed modifications into the weak lensing $\xi_{\pm}$ quantity. Confronted against the cosmic shear data, we reject the $f(R)$ $\{ |f_{R_0}|=10^{-4}, n=1\}$ model with more than 99.9% confidence interval (CI) when assuming a $\Lambda$CDM dark matter only model. In the presence of baryonic feedback processes and massive neutrinos with total mass up to 0.2eV, the model is disfavoured with at least 94% CI in all different combinations studied. Constraints on the $\{ |f_{R_0}|=10^{-4}, n=2\}$ model are weaker, but nevertheless disfavoured with at least 89% CI. We identify several specific combinations of neutrino mass, baryon feedback and $f(R)$ or Dilaton gravity models that are excluded by the current cosmic shear data. Notably, universes with three massless neutrinos and no baryon feedback are strongly disfavoured in all modified gravity scenarios studied. These results indicate that competitive constraints may be achieved with future cosmic shear data.
The so-called "global polytropic model" is applied to the numerical study of the exoplanet systems HD 10180, Kepler-32, Kepler-33, Kepler-102, and Kepler-186. We compare computed distances of planets from their host stars with corresponding observations and discuss some further orbit predictions made by the model.
We test the models of vacuum energy interacting with cold dark matter, and try to probe the possible deviation from the $\Lambda$CDM model using current observations. We focus on two specific models, $Q=3\beta H\rho_{\Lambda}$ and $Q=3\beta H\rho_c$. The data combinations come from the Planck 2013 data, the baryon acoustic oscillations measurements, the Type-Ia supernovae data, the Hubble constant measurement, the redshift space distortions data and the galaxy weak lensing data. For the $Q=3\beta H\rho_c$ model, we find that it can be tightly constrained by all the data combinations, while for the $Q=3\beta H\rho_{\Lambda}$ model there still exist significant degeneracies between parameters. The tightest constraints for the coupling constant are $\beta=-0.026^{+0.036}_{-0.053}$ (for $Q=3\beta H\rho_{\Lambda}$) and $\beta=-0.00045\pm0.00069$ (for $Q=3\beta H\rho_c$) at $1\sigma$ level. For all the fit results, we find that the null interaction $\beta=0$ is always consistent with data. Our work completes the discussion on the interacting dark energy model in the recent planck 2015 papers. Combining with the Planck 2015 results, it is believed that there is no evidence for the models beyond the standard $\Lambda$CDM model from the point of view of possible interaction.
In several recent observational studies on Type Ib/c supernovae (SNe Ib/c), the inferred ejecta masses have a peak value of 2.0 -- 4.0 $M_\odot$, in favor of the binary scenario for their progenitors rather than the Wolf-Rayet star scenario. To investigate the observational properties of relatively low-mass helium stars in binary systems as SN Ib/c progenitors, we constructed atmospheric models with the non-LTE radiative transfer code CMFGEN, using binary star evolution models. We find that these helium stars can be characterized by relatively narrow helium emission lines if the mass-loss rate during the final evolutionary phase is significantly enhanced as implied by many SN Ib/c observations. The optical brightness of helium star progenitors can be meaningfully enhanced with a strong wind for $M \gtrsim 4.4 M_\odot$, but hardly affected or slightly weakened for relatively low-mass of $\sim 3.0 M_\odot$, compared to the simple estimate using blackbody approximation. We further confirm the previous suggestion that the optical brightness would be generally higher for a less massive SN Ib/c progenitor. In good agreement with previous studies, our results indicate that the optical magnitudes and colors of the recently detected progenitor of the SN Ib iPTF13bvn can be well explained by a binary progenitor with a final helium star mass of about 3.0 -- 4.4 $M_\odot$.
Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.
In this work we report a numerical study of the cosmic magnetic field amplification due to collisionless plasma instabilities. The collisionless magnetohydrodynamic equations derived account for the pressure anisotropy that leads, in specific conditions, to the firehose and mirror instabilities. We study the time evolution of seed fields in turbulence under the influence of such instabilities. An approximate analytical time evolution of magnetic field is provided. The numerical simulations and the analytical predictions are compared. We found that i) amplification of magnetic field was efficient in firehose unstable turbulent regimes, but not in the mirror unstable models, ii) the growth rate of the magnetic energy density is much faster than the turbulent dynamo, iii) the efficient amplification occurs at small scales. The analytical prediction for the correlation between the growth timescales with pressure anisotropy ratio is confirmed by the numerical simulations. These results reinforce the idea that pressure anisotropies - driven naturally in a turbulent collisionless medium, e.g. the intergalactic medium -, could efficiently amplify the magnetic field in the early Universe (post-recombination era), previous to the collapse of the first large-scale gravitational structures. This mechanism, though fast for the small scale fields ($\sim$kpc scales), is however unable to provide relatively strong magnetic fields at large scales. Other mechanisms that were not accounted here (e.g., collisional turbulence once instabilities are quenched, velocity shear, or gravitationally induced inflows of gas into galaxies and clusters) could operate afterwards to build up large scale coherent field structures in the long time evolution.
In this work we employ the pure-pseudo formalism devised to minimise the effects of the leakage on the variance of power spectrum estimates and discuss the limits on the tensor-to-scalar ratio, $r$, that could be realistically set by current and forthcoming measurements of the $B$-mode angular power spectrum. We compare those with the results obtained using other approaches: na\"{i}ve mode-counting, minimum-variance quadratic estimators, and re-visit the question of optimizing the sky coverage of small-scale, suborbital experiments in order to maximize the statistical significance of the detection of $r$. We show that the optimized sky coverage is largely insensitive to the adopted approach at least for reasonably compact sky patches. We find, however, that the mode-counting overestimates the detection significance by a factor $\sim1.17$ as compared to the lossless maximum variance approach and by a factor $\sim1.25$ as compared to the lossy pure pseudo-spectrum estimator. In a second time, we consider more realistic experimental configurations. With a pure pseudospectrum reconstruction of $B$-modes and considering only statistical uncertainties, we find that a detection of $r\sim0.11$, $r\sim0.0051$ and $r\sim0.0026$ at 99$\%$ of confidence level is within the reach of current sub-orbital experiments, future arrays of ground-based telescopes and a satellite mission, respectively. This means that an array of telescopes could be sufficient to discriminate between large- and small-field models of inflation, even if the $E$-to-$B$ leakage is consistently included but accounted for in the analysis. However, a satellite mission will be required to distinguish between different small-field models depending on the number of e-folds.
Thick fully depleted CCDs, while enabling wide spectral response, also present challenges in understanding the systematic errors due to 3D charge transport. This 2014 Workshop on Precision Astronomy with Fully Depleted CCDs covered progress that has been made in the testing and modeling of these devices made since a workshop by the same name in 2013. Presentations covered the science drivers, CCD characterization, laboratory measurements of systematics, calibration, and different approaches to modeling the response and charge transport. The key issue is the impact of these CCD sensor features on dark energy science, including astrometry and photometry. Successful modeling of the spatial systematics can enable first order correction in the data processing pipeline.
We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the northwest of M 87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. However, the northwestern part of the cold front is observed to have a nonzero width. The broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intra-cluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are ~ 10 per cent brighter than the surrounding gas and are separated by ~ 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches ~ 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.
Matter bounces are bouncing scenarios wherein the universe contracts as in a matter dominated phase at early times. Such scenarios are known to lead to a scale invariant spectrum of tensor perturbations just as de Sitter inflation does. In this work, we examine if the tensor bi-spectrum can discriminate between the inflationary and the bouncing scenarios. Using the Maldacena formalism, we analytically evaluate the tensor bi-spectrum in a matter bounce for an arbitrary triangular configuration of the wavevectors. We show that, over scales of cosmological interest, the non-Gaussianity parameter $h_{_{\rm NL}}$ that characterizes the amplitude of the tensor bi-spectrum is quite small when compared to the corresponding values in de Sitter inflation. During inflation, the amplitude of the tensor perturbations freeze on super-Hubble scales, a behavior that results in the so-called consistency condition relating the tensor bi-spectrum and the power spectrum in the squeezed limit. In contrast, in the bouncing scenarios, the amplitude of the tensor perturbations grow strongly as one approaches the bounce, which suggests that the consistency condition will not be valid in such situations. We explicitly show that the consistency relation is indeed violated in the matter bounce. We discuss the implications of the results.
We present a series of numerical simulations that explore how the `X-factor', $X_{CO}$ -- the conversion factor between the observed integrated CO emission and the column density of molecular hydrogen -- varies with the environmental conditions in which a molecular cloud is placed. Our investigation is centred around two environmental conditions in particular: the cosmic ray ionisation rate (CRIR) and the strength of the interstellar radiation field (ISRF). Since both these properties of the interstellar medium have their origins in massive stars, we make the assumption in this paper that both the strength of the ISRF and the CRIR scale linearly with the local star formation rate (SFR). The cloud modelling in this study first involves running numerical simulations that capture the cloud dynamics, as well as the time-dependent chemistry, and ISM heating and cooling. These simulations are then post-processed with a line radiative transfer code to create synthetic 12CO (1-0) emission maps from which $X_{CO}$ can be calculated. We find that for 1e4 solar mass virialised clouds with mean density 100 cm$^{-3}$, $X_{CO}$ is only weakly dependent on the local SFR, varying by a factor of a few over two orders of magnitude in SFR. In contrast, we find that for similar clouds but with masses of 1e5 solar masses, the X-factor will vary by an order of magnitude over the same range in SFR, implying that extra-galactic star formation laws should be viewed with caution. However, for denser ($10^4$ cm$^{-3}$), super-virial clouds such as those found at the centre of the Milky Way, the X-factor is once again independent of the local SFR.
Understanding whether M-dwarf stars may host habitable planets with Earth-like atmospheres and biospheres is a major goal in exoplanet research. If such planets exist, the question remains as to whether they could be identified via spectral signatures of biomarkers. Such planets may be exposed to extreme intensities of cosmic rays that could perturb their atmospheric photochemistry. Here, we consider stellar activity of M-dwarfs ranging from quiet up to strong flaring conditions and investigate one particular effect upon biomarkers, namely, the ability of secondary electrons caused by stellar cosmic rays to break up atmospheric molecular nitrogen (N2), which leads to production of nitrogen oxides in the planetary atmosphere, hence affecting biomarkers such as ozone. We apply a stationary model, that is, without a time-dependence, hence we are calculating the limiting case where the atmospheric chemistry response time of the biomarkers is assumed to be slow and remains constant compared with rapid forcing by the impinging stellar flares. This point should be further explored in future work with time-dependent models. For the flaring case O3 is mainly destroyed via direct titration with nitrogen oxides and not via the familiar catalytic cycle photochemistry, which occurs on Earth. For scenarios with low O3, Rayleigh scattering by the main atmospheric gases became more important for shielding the planetary surface from ultra-violet radiation. A major result of this work is that the biomarker O3 survived all the stellar-activity scenarios considered except for the strong case, whereas the biomarker nitrous oxide could survive in the planetary atmosphere under all conditions of stellar activity considered here, which clearly has important implications for missions that aim to detect spectroscopic biomarkers.
We present equilibrium models of relativistic magnetised, infinite,
axisymmetric jets with rotation propagating through an homogeneous,
unmagnetised ambient medium at rest. The jet models are characterised by six
functions defining the radial profiles of density, pressure, and the toroidal
and axial components of velocity and magnetic field. Fixing the ambient
pressure and the jet rest-mass density and axial components of the flow
velocity and magnetic field, we analyze the influence of the toroidal magnetic
field and several rotation laws on the structure of the equilibrium models. Our
approach excludes by construction the analysis of the self-consistently
magnetically launched jet models or the force-free equilibrium solutions.
Several forbidden regions in the magnetic pitch angle/magnetization plane are
found where models of the class considered in our study could not be settled.
These forbidden regions are associated with the existence of maximum axial and
toroidal magnetic field components compatible with the prescribed equilibrium
condition at the jet surface, and/or an excess of centrifugal force producing
gaps with negative pressures in the jet. The present study can be easily
extended to jet models with different transversal profiles and magnetic field
configurations.
In the last part of the paper, we test the ability of our RMHD code to
maintain steady equilibrium models of axisymmetric RMHD jets in one and two
spatial dimensions. The one dimensional numerical simulations serve also as a
consistency proof of the fidelity of the analytical steady solutions discussed
in the first part of the paper. The present study allows us to build initial
equilibrium jet models with selected properties for dynamical (and emission)
simulations of magnetised relativistic jets with rotation.
We study the origin and escape of catastrophically outgassed volatiles (H$_2$O, CO$_2$) from exomoons with Earth-like densities and masses of $0.1M_{\oplus}$, $0.5M_{\oplus}$ and $1M_{\oplus}$ orbiting an extra-solar gas giant inside the habitable zone of a young active solar-like star. We apply a radiation absorption and hydrodynamic upper atmosphere model to the three studied exomoon cases. We model the escape of hydrogen and dragged dissociation products O and C during the activity saturation phase of the young host star. Because the soft X-ray and EUV radiation of the young host star may be up to $\sim$100 times higher compared to today's solar value during the first 100 Myr after the system's origin, an exomoon with a mass $ < 0.25M_{\oplus}$ located in the HZ may not be able to keep an atmosphere because of its low gravity. Depending on the spectral type and XUV activity evolution of the host star, exomoons with masses between $\sim0.25-0.5M_{\oplus}$ may evolve to Mars-like habitats. More massive bodies with masses $ > 0.5M_{\oplus}$, however, may evolve to habitats that are a mixture of Mars-like and Earth-analogue habitats, so that life may originate and evolve at the exomoon's surface.
We report the result of a systematic methanol observation toward IRAS 19312+1950. The properties of the SiO, H2O and OH masers of this object are consistent with those of mass-losing evolved stars, but some other properties are difficult to explain in the standard scheme of stellar evolution in its late stage. Interestingly, a tentative detection of radio methanol lines was suggested toward this object by a previous observation. To date, there are no confirmed detections of methanol emission towards evolved stars, so investigation of this possible detection is important to better understand the circumstellar physical/chemical environment of IRAS 19312+1950. In this study, we systematically observed multiple methanol lines of IRAS 19312+1950 in the lambda=3mm, 7mm, and 13mm bands, and detected 6 lines including 4 thermal lines and 2 class I maser lines. We derived basic physical parameters including kinetic temperature and relative abundances by fitting a radiative transfer model. According to the derived excitation temperature and line profiles, a spherically expanding outflow lying at the center of the nebulosity is excluded from the possibilities for methanol emission regions. The detection of class I methanol maser emission suggests that a shock region is involved in the system of IRAS 19312+1950. If the central star of IRAS 19312+1950 is an evolved star as suggested in the past, the class I maser detected in the present observation is the first case detected in an interaction region between an evolved star outflow and ambient molecular gas.
We study the Extended Chandra Deep Field South (E-CDFS) Very Large Array sample, which reaches a flux density limit at 1.4 GHz of 32.5 microJy at the field centre and redshift ~ 4, and covers ~ 0.3 deg^2. Number counts are presented for the whole sample while the evolutionary properties and luminosity functions are derived for active galactic nuclei (AGN). The faint radio sky contains two totally distinct AGN populations, characterised by very different evolutions, luminosity functions, and Eddington ratios: radio-quiet (RQ)/radiative-mode, and radio-loud/jet-mode AGN. The radio power of RQ AGN evolves ~ (1+z)^2.5, similarly to star-forming galaxies, while the number density of radio-loud ones has a peak at ~ 0.5 and then declines at higher redshifts. The number density of radio-selected RQ AGN is consistent with that of X-ray selected AGN, which shows that we are sampling the same population. The unbiased fraction of radiative-mode RL AGN, derived from our own and previously published data, is a strong function of radio power, decreasing from ~ 0.5 at P_1.4GHz ~ 10^24 W/Hz to ~ 0.04$ at P_1.4GHz ~ 10^22 W/Hz. Thanks to our enlarged sample, which now includes ~ 700 radio sources, we also confirm and strengthen our previous results on the source population of the faint radio sky: star-forming galaxies start to dominate the radio sky only below ~ 0.1 mJy, which is also where radio-quiet AGN overtake radio-loud ones.
The European Virtual Observatory Euro-VO has been coordinating European VO activities through a series of projects co-funded by the European Commission over the last 15 years. The bulk of VO work in Europe is ensured by the national VO initiatives and those of intergovernmental agencies. VO activities at the European level coordinate the work in support of the three "pillars" of the Virtual Observatory: support to the scientific community, take-up by the data providers, and technological activities. Several Euro-VO projects have also provided direct support to selected developments and prototyping. This paper explains the methodology used by Euro-VO over the years. It summarizes the activities which were performed and their evolutions at different stages of the development of the VO, explains the Euro-VO role with respect to the international and national levels of VO activities, details the lessons learnt for best practices for the coordination of the VO building blocks, and the liaison with other European initiatives, documenting the added-value of European coordination. Finally, the current status and next steps of Euro-VO are briefly addressed.
The evolution and escape of the martian atmosphere and the planet's water inventory can be separated into an early and late evolutionary epoch. The first epoch started from the planet's origin and lasted $\sim$500 Myr. Because of the high EUV flux of the young Sun and Mars' low gravity it was accompanied by hydrodynamic blow-off of hydrogen and strong thermal escape rates of dragged heavier species such as O and C atoms. After the main part of the protoatmosphere was lost, impact-related volatiles and mantle outgassing may have resulted in accumulation of a secondary CO$_2$ atmosphere of a few tens to a few hundred mbar around $\sim$4--4.3 Gyr ago. The evolution of the atmospheric surface pressure and water inventory of such a secondary atmosphere during the second epoch which lasted from the end of the Noachian until today was most likely determined by a complex interplay of various nonthermal atmospheric escape processes, impacts, carbonate precipitation, and serpentinization during the Hesperian and Amazonian epochs which led to the present day surface pressure.
The energy spectrum of the cosmic microwave background (CMB) provides a powerful tool for constraining standard and non-standard physics in the primordial Universe. Previous studies mainly highlight spectral distortions (mu-, y- and r-type) created by episodes of early energy release; however, several processes also lead to copious photon production, which requires a different treatment. Here, we carry out a first detailed study for the evolution of distortions caused by photon injection at different energies in the CMB bands. We provide detailed analytical and numerical calculations illustrating the rich phenomenology of the associated distortion signals. We show that photon injection at very high and very low frequencies creates distortions that are similar to those from pure energy release. In the mu-era (z>3x10^5), a positive or negative chemical potential can be formed, depending on the balance between added photon energy and number. At lower redshifts (z<3x10^5), partial information about the photon injection process (i.e., injection time and energy) can still be recovered, with the distortion being found in a partially comptonized state. We briefly discuss current and future constraints on scenarios with photon production. We also argue that more detailed calculations for different scenarios with photon injection may be required to assess in which regimes these can be distinguished from pure energy release processes.
By considering martian-like planetary embryos inside the habitable zone of solar-like stars we study the behavior of the hydrodynamic atmospheric escape of hydrogen for small values of the Jeans escape parameter $\beta < 3$, near the base of the thermosphere, that is defined as a ratio of the gravitational and thermal energy. Our study is based on a 1-D hydrodynamic upper atmosphere model that calculates the volume heating rate in a hydrogen dominated thermosphere due to the absorption of the stellar soft X-ray and extreme ultraviolet (XUV) flux. We find that when the $\beta$ value near the mesopause/homopause level exceeds a critical value of $\sim$2.5, there exists a steady hydrodynamic solution with a smooth transition from subsonic to supersonic flow. For a fixed XUV flux, the escape rate of the upper atmosphere is an increasing function of the temperature at the lower boundary. Our model results indicate a crucial enhancement of the atmospheric escape rate, when the Jeans escape parameter $\beta$ decreases to this critical value. When $\beta$ becomes $\leq$2.5, there is no stationary hydrodynamic transition from subsonic to supersonic flow. This is the case of a fast non-stationary atmospheric expansion that results in extreme thermal atmospheric escape rates.
A stationary distribution function that describes the entire processes of propagation of relativistic particles, including the transition between the ballistic and diffusion regimes, is obtained. The spacial component of the constructed function satisfies to the first two moments of the Boltzmann equation. The angular part of the distribution provides accurate values for the angular moments derived from the Boltzmann equation, and gives a correct expression in the limit of small-angle approximation. Using the derived function, we studied the gamma-ray images produced through the $pp$ interaction of relativistic particles with gas clouds in the proximity of the accelerator. In general, the morphology and the energy spectra of gamma-rays significantly deviate from the "standard" results corresponding to the propagation of relativistic particles strictly in the diffusion regime.
The equations of smoothed particle magnetohydrodynamics (SPMHD), even with the various corrections to instabilities so far proposed, have been observed to be unstable when a very steep density gradient is necessarily combined with a variable smoothing length formalism. Here we consider in more detail the modifications made to the SPMHD equations in LBP2015 that resolve this instability by replacing the smoothing length in the induction and anisotropic force equations with an average smoothing length term. We then explore the choice of average used and compare the effects on a test `cylinder-in-a-box' problem and the collapse of a magnetised molecular cloud core. We find that, aside from some benign numerical effects at low resolutions for the quadratic mean, the formalism is robust as to the choice of average but that in complicated models it is essential to apply the average to both equations; in particular, all four averages considered exhibit similar conservation properties. This improved formalism allows for arbitrarily small sink particles and field geometries to be explored, vastly expanding the range of astronomical problems that can be modeled using SPMHD.
Redshift-space distortions are generally considered in the plane parallel limit, where the angular separation between the two sources can be neglected. Given that galaxy catalogues now cover large fractions of the sky, it becomes necessary to consider them in a formalism which takes into account the wide angle separations. In this article we derive an operational formula for the matter correlators in the Newtonian limit to be used in actual data sets, both in configuration and in Fourier spaces without relying on a plane-parallel approximation. We then recover the plane-parallel limit not only in configuration space where the geometry is simpler, but also in Fourier space, and we exhibit the first corrections that should be included in large surveys as a perturbative expansion over the plane-parallel results. We finally compare our results to existing literature, and show explicitly how they are related.
The outer solar system may provide a potential habitat for extraterrestrial life. Remote sensing data from the Galileo spacecraft suggest that the jovian icy moons, Europa, Ganymede, and possibly Callisto, may harbor liquid water oceans underneath their icy crusts. Although compositional information required for the discussion of habitability is limited because of significantly restricted observation data, organic molecules are ubiquitous in the universe. Recently, in-situ spacecraft measurements and experiments suggest that amino acids can be formed abiotically on interstellar ices and comets. These amino acids could be continuously delivered by meteorite or comet impacts to icy moons. Here, we show that polymerization of organic monomers, in particular amino acids and nucleotides, could proceed spontaneously in the cold environment of icy moons, in particular the Jovian icy moon Europa as a typical example, based on thermodynamic calculations, though kinetics of formation are not addressed. Observed surface temperature on Europa is 120 and 80 K in the equatorial region and polar region, respectively. At such low temperatures, Gibbs energies of polymerization become negative, and the estimated thermal structure of the icy crust should contain a shallow region (i.e., at a depth of only a few kilometers) favorable for polymerization. Investigation of the possibility of organic monomer polymerization on icy moons could provide good constraints on the origin and early evolution of extraterrestrial life.
We analyze 10 UV time series for 5 stars which fulfill specific sampling and spectral criteria to constrain the origin of large-scale wind structure in O stars. We argue that excited state lines must arise close to the stellar surface and are an excellent diagnostic complement to resonance lines which, due to radiative transfer effects, rarely show variability at low velocity. Consequently, we splice dynamic spectra of the excited state line, N IV1718, at low velocity to those of 1393 component of the Si IV 1400 doublet at high velocity in order to examine the temporal evolution of wind line features. These spliced time series reveal that nearly all of the features observed in the time series originate at or very near the stellar surface. Further, we positively identify the observational signature of equatorial co-rotating interaction regions in two of the five stars and possibly two others. In addition, we see no evidence for features originating further out in the wind. We use our results to consolidate the fact that the features seen in dynamic spectra must be huge, in order to remain in the line of sight for days, persisting to very large velocity and that the photospheric footprint of the features must also be quite large, ~ 15 - 20% of the stellar diameter.
We demonstrate the use of machine learning algorithms in combination with segmentation techniques in order to distinguish coronal holes and filaments in SDO/AIA EUV images of the Sun. Based on two coronal hole detection techniques (intensity-based thresholding, SPoCA), we prepared data sets of manually labeled coronal hole and filament channel regions present on the Sun during the time range 2011 - 2013. By mapping the extracted regions from EUV observations onto HMI line-of-sight magnetograms we also include their magnetic characteristics. We computed shape measures from the segmented binary maps as well as first order and second order texture statistics from the segmented regions in the EUV images and magnetograms. These attributes were used for data mining investigations to identify the most performant rule to differentiate between coronal holes and filament channels. We applied several classifiers, namely Support Vector Machine, Linear Support Vector Machine, Decision Tree, and Random Forest and found that all classification rules achieve good results in general, with linear SVM providing the best performances (with a true skill statistic of ~0.90). Additional information from magnetic field data systematically improves the performance across all four classifiers for the SPoCA detection. Since the calculation is inexpensive in computing time, this approach is well suited for applications on real-time data. This study demonstrates how a machine learning approach may help improve upon an unsupervised feature extraction method.
Rapidly rotating giant stars are relatively rare and may represent important stages of stellar evolution, resulting from stellar coalescence of close binary systems or accretion of sub-stellar companions by their hosting stars. In the present letter we report 17 giant stars observed in the scope of the Kepler space mission exhibiting rapid rotation behavior. For the first time the abnormal rotational behavior for this puzzling family of stars is revealed by direct measurements of rotation, namely from photometric rotation period, exhibiting very short rotation period with values ranging from 13 to 55 days. This finding points for remarkable surface rotation rates, up to 18 times the Sun rotation. These giants are combined with 6 other recently listed in the literature for mid-IR diagnostic based on WISE information, from which a trend for an infrared excess is revealed for at least a half of the stars, but at a level far lower than the dust excess emission shown by planet-bearing main-sequence stars.
(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...
We show that inflationary models with broken rotational invariance generate testable off-diagonal signatures in the correlation between the $\mu$-type distortion and temperature fluctuations of the Cosmic Microwave Background. More precisely, scenarios with a quadrupolar bispectrum asymmetry, usually generated by fluctuations of primordial vector fields, produce a non-vanishing $\mu$-$T$ correlation when $|\ell_1-\ell_2|=2$. Since spectral distortions are sensitive to primordial fluctuations up to very small scales, a cosmic variance limited spectral distortion experiment can detect such effects with high signal-to-noise.
Asteroseismology is a powerful tool that can precisely characterize the mass, radius, and other properties of field stars. However, our inability to properly model the near-surface layers of stars creates a frequency-dependent frequency difference between the observed and the modeled frequencies, usually referred to as the "surface term". This surface term can add significant errors to the derived stellar properties unless removed properly. In this paper we simulate surface terms across a significant portion of the HR diagram, exploring four different masses ($M=0.8, 1.0, 1.2$, and $1.5$ M$_\odot$) at five metallicities ($[\rm{Fe/H}]=0.5, 0.0, -0.5 ,-1.0, and -1.5$) from main sequence to red giants for stars with $T_{\rm{eff}}<6500 K$ and explore how well the most common ways of fitting and removing the surface term actually perform. We find that the two-term model proposed by Ball & Gizon (2014) works much better than other models across a large portion of the HR diagram, including the red giants, leading us to recommend its use for future asteroseismic analyses.
With all the new discoveries being made about nearby young stars, the ability to find new nearby young stars is as important as ever, and membership identification codes will continue to perform a vital role in scientific research. In the process of creating a new moving group membership identification code - LocAting Constituent mEmbers In Nearby Groups (LACEwING) - we have discovered a few pointers relevant to astronomers trying to use codes like LACEwING to locate young stars.
Satellites of asteroids have been discovered in nearly every known small body
population, and a remarkable aspect of the known satellites is the diversity of
their properties. They tell a story of vast differences in formation and
evolution mechanisms that act as a function of size, distance from the Sun, and
the properties of their nebular environment at the beginning of Solar System
history and their dynamical environment over the next 4.5 Gyr. The mere
existence of these systems provides a laboratory to study numerous types of
physical processes acting on asteroids and their dynamics provide a valuable
probe of their physical properties otherwise possible only with spacecraft.
Advances in understanding the formation and evolution of binary systems have
been assisted by: 1) the growing catalog of known systems, increasing from 33
to nearly 250 between the Merline et al. (2002) Asteroids III chapter and now,
2) the detailed study and long-term monitoring of individual systems such as
1999 KW4 and 1996 FG3, 3) the discovery of new binary system morphologies and
triple systems, 4) and the discovery of unbound systems that appear to be
end-states of binary dynamical evolutionary paths.
Specifically for small bodies (diameter smaller than 10 km), these
observations and discoveries have motivated theoretical work finding that
thermal forces can efficiently drive the rotational disruption of small
asteroids. Long-term monitoring has allowed studies to constrain the system's
dynamical evolution by the combination of tides, thermal forces and rigid body
physics. The outliers and split pairs have pushed the theoretical work to
explore a wide range of evolutionary end-states.
High-precision spectroscopy of large stellar samples plays a crucial role for several topical issues in astrophysics. Examples include studying the chemical structure and evolution of the Milky Way galaxy, tracing the origin of chemical elements, and characterizing planetary host stars. Data are accumulating from instruments that obtain high-quality spectra of stars in the ultraviolet, optical and infrared wavelength regions on a routine basis. These instruments are located at ground-based 2- to 10-m class telescopes around the world, in addition to the spectrographs with unique capabilities available at the Hubble Space Telescope. The interpretation of these spectra requires high-quality transition data for numerous species, in particular neutral and singly ionized atoms, and di- or triatomic molecules. We rely heavily on the continuous efforts of laboratory astrophysics groups that produce and improve the relevant experimental and theoretical atomic and molecular data. The compilation of the best available data is facilitated by databases and electronic infrastructures such as the NIST Atomic Spectra Database, the VALD database, or the Virtual Atomic and Molecular Data Centre (VAMDC). We illustrate the current status of atomic data for optical stellar spectra with the example of the Gaia-ESO Public Spectroscopic Survey. Data sources for 35 chemical elements were reviewed in an effort to construct a line list for a homogeneous abundance analysis of up to 100000 stars.
We present the results of our extensive search using the Bayesian block method for long tails following short bursts from a magnetar, SGR J1550-5418, over all RXTE observations of the source. We identified four bursts with extended tails, most of which occurred during its 2009 burst active episode. The durations of tails range between ~13 s and over 3 ks, which are much longer than the typical duration of bursts. We performed detailed spectral and temporal analysis of the burst tails. We find that the spectra of three tails show a thermal nature with a trend of cooling throughout the tail. We compare the results of our investigations with the properties of four other extended tails detected from SGR 1900+14 and SGR 1806-20 and suggest a scenario for the origin of the tail in the framework of the magnetar model.
We use a Monte Carlo code to calculate the geodesic orbits of test particles around Kerr black holes, generating a distribution function of both bound and unbound populations of dark matter particles. From this distribution function, we calculate annihilation rates and observable gamma-ray spectra for a few simple dark matter models. The features of these spectra are sensitive to the black hole spin, observer inclination, and detailed properties of the dark matter annihilation cross section and density profile. Confirming earlier analytic work, we find that for rapidly spinning black holes, the collisional Penrose process can reach efficiencies exceeding $600\%$, leading to a high-energy tail in the annihilation spectrum. The high particle density and large proper volume of the region immediately surrounding the horizon ensures that the observed flux from these extreme events is non-negligible.
Circumbinary bodies are objects that orbit around a more massive binary system. Here we show that, contrarily to the classical two-body problem, circumbinary bodies in planar quasi-circular orbits can present stable non-synchronous rotation. Denoting $n_b$ and $n$ the orbital mean motion of the binary and of the circumbinary body, respectively, there is an entirely new family of spin-orbit resonances at the frequencies $n\pm k\nu/2$, where $\nu = n_b - n$, and $k$ is an integer. In addition, when the natural rotational libration frequency has the same magnitude as $\nu$, the individual resonances overlap and the rotation becomes chaotic. We apply these results to the small satellites in the Pluto-Charon system. We conclude that the rotation of Nix and Styx can be chaotic, and that the rotation of Hydra and Kerberos is stable but not necessarily synchronous.
Carbonaceous chondrite meteorites are known for having high water and organic material contents, including amino acids. Here we address the origin of amino acids in the warm interiors of their parent bodies (planetesimals) within a few million years of their formation, and connect this with the astrochemistry of their natal protostellar disks. We compute both the total amino acid abundance pattern as well as the relative frequencies of amino acids within the CM2 (e.g. Murchison) and CR2 chondrite subclasses based on Strecker reactions within these bodies. We match the relative frequencies to well within an order of magnitude among both CM2 and CR2 meteorites for parent body temperatures $<$ 200$^{\circ}$C. These temperatures agree with 3D models of young planetesimal interiors. We find theoretical abundances of approximately 7x10$^5$ parts-per-billion (ppb), which is in agreement with the average observed abundance in CR2 meteorites of 4$\pm$7x10$^5$, but an order of magnitude higher than the average observed abundance in CM2 meteorites of 2$\pm$2x10$^4$. We find that the production of hydroxy acids could be favoured over the production of amino acids within certain meteorite parent bodies (e.g. CI1, CM2) but not others (e.g. CR2). This could be due to the relatively lower NH$_3$ abundances within CI1 and CM2 meteorite parent bodies, which leads to less amino acid synthesis. We also find that the water content in planetesimals is likely to be the main cause of variance between carbonaceous chondrites of the same subclass. We propose that amino acid abundances are primarily dependent on the ammonia and water content of planetesimals that are formed in chemically distinct regions within their natal protostellar disks.
The measured masses of the Higgs boson and top quark indicate that the effective potential of standard model either develops an unstable electroweak vacuum or stands stable all the way up to the Planck scale. In the latter case in which the top quark mass is about $2\sigma$ below its present central value, the Higgs boson can be the inflaton with the help of a large non-minimal coupling to curvature in four dimensions. We propose a scenario in which the Higgs boson can be the inflaton in five-dimensional Gauss-Bonnet brane-world model to solve both the unitarity and stability problems which usually plague Higgs inflation. We find that in order Higgs inflation to happen successfully in Gauss-Bonnet regime, the extra dimension scale must appear roughly in the range between the TeV scale and the instability scale of standard model. At the tree level, our model can give rise to a naturally small non-minimal coupling $\xi\sim\mathcal{O}(1)$ for the Higgs quartic coupling $\lambda\sim\mathcal{O}(0.1)$ if the extra dimension scale lies at the TeV scale. At the loop level, the inflationary predictions at the tree-level are preserved. Our model can be confronted with future experiments and observations from both particle physics and cosmology.
Linearized gravitational waves in de Sitter space-time are analyzed in detail to obtain guidance for constructing the theory of gravitational radiation in presence of a positive cosmological constant in full, nonlinear general relativity. Specifically: i) In the exact theory, the intrinsic geometry of $\scri$ is often assumed to be conformally flat in order to reduce the asymptotic symmetry group from $\Diff$ to the de Sitter group. Our {results show explicitly} that this condition is physically unreasonable; ii) We obtain expressions of energy-momentum and angular momentum fluxes carried by gravitational waves in terms of fields defined at $\scrip$; iii) We argue that, although energy of linearized gravitational waves can be arbitrarily negative in general, gravitational waves emitted by physically reasonable sources carry positive energy; and, finally iv) We demonstrate that the flux formulas reduce to the familiar ones in Minkowski space-time in spite of the fact that the limit $\Lambda \to 0$ is discontinuous (since, in particular, $\scri$ changes its space-like character to null in the limit).
We study numerically the high energy head-on collision of nonspinning equal mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the full numerical evolutions and computation of waveforms at infinity. We make use of new initial data with notably reduced spurious radiation content, allowing for initial speeds nearing the speed of light, i.e. $v\sim0.99c$. We thus estimate the maximum radiated energy from head-on collisions to be $E_{\text{max}}/M_{ADM}=0.13\pm0.01$. This value differs from the second order perturbative $(0.164)$ and zero-frequency-limit $(0.17)$ analytic computations, but is close to those obtained by thermodynamic arguments $(0.134)$ and to previous numerical estimates $(0.14\pm0.03)$.
A quintessential Inflation (QI)scenario from Mimetic Dark Matter(MDM) is presented in this paper. This scenario, which is based on the MDM model presented by Chamseddine and Mukhanov \cite{chams},uses a potential that is defined on two time intervals, one during inflation, and the other after it. The resulting energy density of the universe is constant during inflation, followed by that of a matter/radiation dominated universe, and finally ends with a constant energy density corresponding to dark energy. The scale factor has an accelerating expansion nature during and after inflation. It will be shown how this is still a viable scenario, even if the scale factor after inflation is not that of a decelerating De Sitter universe.
We respond to the criticisms of a recent paper of Buchert et al. [arXiv:1505.07800]
Collisions of particles near a rotating black hole can lead to unbound energies E_{c.m.} in their centre of mass frame. There are indications that the Killing energy of debris at infinity can also be unbound for some scenarios of collisions near the extremal black hole horizon (so-called super-Penrose process). They include participation of a particle that (i) has generic (not fine-tuned) parameters and (ii) moves away from a black hole after collision. We show that for any finite particle's mass, such a particle cannot be obtained as a result of preceding collision. However, this can be done if one of initial infalling particles has the mass of the order $N^{-2}$ that generalizes previous observation made in literature for radial infall in the Kerr background.
We consider the Standard Model with a new particle which is charged under $SU(2)_{L}$ with the hypercharge being zero. Such a particle is known as one of the dark matter (DM) candidates. We examine the realization of the multiple point criticality principle (MPP) in this class of models. Namely, we investigate whether the one-loop effective Higgs potential $V_{\text{eff}}(\phi)$ and its derivative $dV_{\text{eff}}(\phi)/d \phi$ can become simultaneously zero at around the string/Planck scale, based on the one/two-loop renormalization group equations. As a result, we find that only the $SU(2)_L$ triplet extensions can realize the MPP. More concretely, in the case of the triplet Majorana fermion, the MPP is realized at the scale $\phi=3.6\times10^{16}$ GeV if the top mass $M_{t}$ is $172.2$ GeV. On the other hand, for the real triplet scalar, the MPP can be satisfied for $10^{16}\text{ GeV}\lesssim\phi\lesssim10^{17}$GeV and $172\text{ GeV}\gtrsim M_{t}\gtrsim171$ GeV, depending on the coupling between the Higgs and DM.
In the initial stage of relativistic heavy-ion collisions, strong magnetic fields appear due to the large velocity of the colliding charges. The evolution of these fields appears as a novel and intriguing feature in the fluid-dynamical description of heavy-ion collisions. In this work, we study analytically the one-dimensional, longitudinally boost-invariant motion of an ideal fluid in the presence of a transverse magnetic field. Interestingly, we find that, in the limit of ideal magnetohydrodynamics, i.e., for infinite conductivity, and irrespective of the strength of the initial magnetization, the decay of the fluid energy density $e$ with proper time $\tau$ is the same as for the time-honored "Bjorken flow" without magnetic field. Furthermore, when the magnetic field is assumed to decay $\sim \tau^{-a}$, where $a$ is an arbitrary number, two classes of analytic solutions can be found depending on whether $a$ is larger or smaller than one. In summary, the analytic solutions presented here highlight that the Bjorken flow is far more general than formerly thought. These solutions can serve both to gain insight on the dynamics of heavy-ion collisions in the presence of strong magnetic fields and as testbeds for numerical codes.
We investigate cusp singularities in f(R) gravity, especially for Starobinsky and Hu-Sawicki dark energy models. We illustrate that, by using double-null numerical simulations, a cusp singularity can be triggered by gravitational collapses. This singularity can be cured by adding a quadratic term, but this causes a Ricci scalar bump that can be observed by an observer outside the event horizon. Comparing with cosmological parameters, it seems that it would be difficult to see super-Planckian effects by astrophysical experiments. On the other hand, at once there exists a cusp singularity, it can be a mechanism to realize a horizon scale curvature singularity that can be interpreted by a firewall.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
The purpose of this mock data and science challenge is to prepare the data analysis and science interpretation for the second generation of gravitational-wave experiments Advanced LIGO-Virgo in the search for a stochastic gravitational-wave background signal of astrophysical origin. Here we present a series of signal and data challenges, with increasing complexity, whose aim is to test the ability of current data analysis pipelines at detecting an astrophysically produced gravitational-wave background, test parameter estimation methods and interpret the results. We introduce the production of these mock data sets that includes a realistic observing scenario data set where we account for different sensitivities of the advanced detectors as they are continuously upgraded toward their design sensitivity. After analysing these with the standard isotropic cross-correlation pipeline we find that we are able to recover the injected gravitational-wave background energy density to within $2\sigma$ for all of the data sets and present the results from the parameter estimation. The results from this mock data and science challenge show that advanced LIGO and Virgo will be ready and able to make a detection of an astrophysical gravitational-wave background within a few years of operations of the advanced detectors, given a high enough rate of compact binary coalescing events.
Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asymmetric heating on local outflow structure. We use the Athena code for hydrodynamics to produce 3D simulations of hot Jupiter mass loss that jointly model wind launching and stellar heating via photoionization. Our fiducial planet is an inflated, hot Jupiter with radius $R_p=2.14 R_{\rm Jup}$ and mass $M_p = 0.53 M_{\rm Jup}$. We irradiate the initially neutral, atomic hydrogen atmosphere with 13.6 eV photons and compute the outflow's ionization structure. There are clear asymmetries in the atmospheric outflow, including a neutral shadow on the planet's nightside. Given an incident ionizing UV flux comparable to that of the Sun, we find a steady-state mass loss rate of ~$2\times10^{10}$ g s$^{-1}$. The total mass loss rate and the outflow substructure along the substellar ray show good agreement with earlier 1D models, for two different fluxes. Our 3D data cube can be used to generate the outflow's extinction spectrum during transit. As a proof of concept, we find absorption of stellar Lyman-alpha at Doppler-shifted velocities of up to $\pm 50$ km s$^{-1}$. Our work provides a starting point for further 3D models that can be used to predict observable signatures of hot Jupiter mass loss.
We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\log\rho/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms.
Recent chemical abundance measurements of damped Ly-alpha absorbers (DLAs) revealed an intrinsic scatter in their metallicity of ~0.5 dex out to z~5. In order to explore the origin of this scatter, we build a semi-analytic model which traces the chemical evolution of the interstellar matter in small regions of the Universe with different mean density, from over- to underdense regions. We show that the different histories of structure formation in these regions, namely halo abundance, mass and stellar content, is reflected in the chemical properties of the protogalaxies, and in particular of DLAs. We calculate mean metallicity-redshift relations and show that the metallicity dispersion arising from this environmental effect amounts to ~0.25 dex and is an important contributor to the observed overall intrinsic scatter.
We present a dynamical analysis to infer the structural parameters and properties of the two nearby, compact, high velocity dispersion galaxies MRK1216 & NGC1277. Combining deep HST imaging, wide-field IFU stellar kinematics, and complementary long-slit spectroscopic data out to 3 R_e, we construct orbit-based models to constrain their black hole masses, dark matter content and stellar mass-to-light ratios. We obtain a black hole mass of log(Mbh/Msun) = 10.1(+0.1/-0.2) for NGC1277 and an upper limit of log(Mbh/Msun) = 10.0 for MRK1216, within 99.7 per cent confidence. The stellar mass-to-light ratios span a range of Upsilon_V = 6.5(+1.5/-1.5) in NGC1277 and Upsilon_H = 1.8(+0.5/-0.8) in MRK1216 and are in good agreement with SSP models of a single power-law Salpeter IMF. Even though our models do not place strong constraints on the dark halo parameters, they suggest that dark matter is a necessary ingredient in MRK1216, with a dark matter contribution of 22(+30/-20) per cent to the total mass budget within 1 R_e. NGC1277, on the other hand, can be reproduced without the need for a dark halo, and a maximal dark matter fraction of 13 per cent within the same radial extent. In addition, we investigate the orbital structures of both galaxies, which are rotationally supported and consistent with photometric multi-S\'ersic decompositions, indicating that these compact objects do not host classical, non-rotating bulges formed during recent (z <= 2) dissipative events or through violent relaxation. Finally, both MRK 1216 and NGC 1277 are anisotropic, with a global anisotropy parameter delta of 0.33 and 0.58, respectively. While MRK 1216 follows the trend of fast-rotating, oblate galaxies with a flattened velocity dispersion tensor in the meridional plane of the order of beta_z = delta, NGC 1277 is highly tangentially anisotropic and seems to belong kinematically to a distinct class of objects.
We report detections of two candidate distant submillimeter galaxies (SMGs), MM J154506.4$-$344318 and MM J154132.7$-$350320, which are discovered in the AzTEC/ASTE 1.1 mm survey toward the Lupus-I star-forming region. The two objects have 1.1 mm flux densities of 43.9 and 27.1 mJy, and have Herschel/SPIRE counterparts as well. The Submillimeter Array counterpart to the former SMG is identified at 890 $\mu$m and 1.3 mm. Photometric redshift estimates using all available data from the mid-infrared to the radio suggest that the redshifts of the two SMGs are $z_{\rm photo} \simeq$ 4-5 and 3, respectively. Near-infrared objects are found very close to the SMGs and they are consistent with low-$z$ ellipticals, suggesting that the high apparent luminosities can be attributed to gravitational magnification. The cumulative number counts at $S_{\rm 1.1mm} \ge 25$ mJy, combined with other two 1.1-mm brightest sources, are $0.70 ^{+0.56}_{-0.34}$ deg$^{-2}$, which is consistent with a model prediction that accounts for flux magnification due to strong gravitational lensing. Unexpectedly, a $z > 3$ SMG and a Galactic dense starless core (e.g., a first hydrostatic core) could be similar in the mid-infrared to millimeter spectral energy distributions and spatial structures at least at $\gtrsim 1"$. This indicates that it is necessary to distinguish the two possibilities by means of broad band photometry from the optical to centimeter and spectroscopy to determine the redshift, when a compact object is identified toward Galactic star-forming regions.
We report rotational periods for 16 young brown dwarfs in the nearby Upper Scorpius association, based on 72 days of high-cadence, high-precision photometry from the Kepler space telescope's K2 mission. The periods range from a few hours to two days (plus one outlier at 5 days), with a median just above one day, confirming that brown dwarfs, except at the very youngest ages, are fast rotators. Interestingly, four of the slowest rotators in our sample exhibit mid-infrared excess emission from disks; at least two also show signs of disk eclipses and accretion in the lightcurves. Comparing these new periods with those for two other young clusters and simple angular momentum evolution tracks, we find little or no rotational braking in brown dwarfs between 1-10 Myr, in contrast to low-mass stars. Our findings show that disk braking, while still at work, is inefficient in the substellar regime, thus provide an important constraint on the mass dependence of the braking mechanism.
We simulate structure formation in the matter-dominated universe using a quasi-newtonian, one-dimensional model. In addition to dark matter, luminous matter is introduced to examine the potential bias in the distributions. We use multifractal analysis techniques to identify structures, including clusters and voids. Both dark matter and luminous matter exhibit fractal geometry as the universe evolves. We present the results on the generalized dimensions computed at various scales for each matter distribution.
(abridged) We report on multi-wavelength measurements of Swift J1753.5-0127 in the hard state at L=2.7e36 erg/s (assuming d=3 kpc) in 2014. The radio emission is optically thick synchrotron, presumably from a compact jet. We take advantage of the low extinction and model the near-IR to UV emission with a multi-temperature disk model. Assuming a BH mass of M_BH=5 Msun and a system inclination of 40 deg, the fits imply an inner radius for the disk of Rin/Rg>212 d_3 (5Msun/M_BH). The outer radius is R_out/R_g=90,000 d_3 (5Msun/M_BH), which corresponds to 6.6e10 d_3 cm, consistent with the expected size of the disk. The 0.5-240 keV spectrum measured by Swift/XRT, Suzaku, and NuSTAR is relatively well characterized by a power-law with a photon index of Gamma=1.722+/-0.003, but a significant improvement is seen when a second continuum component is added. Reflection is a possibility, but no iron line is detected, implying a low iron abundance. We are able to fit the entire SED with a multi-temperature disk component, a Comptonization component, and a broken power-law, representing the emission from the compact jet. The broken power-law cannot significantly contribute to the soft X-ray emission, and this may be related to why Swift J1753.5-0127 is an outlier in the radio/X-ray correlation. The broken power-law might dominate above 20 keV, which would constrain the break frequency to be between 2.4e10 Hz and 3.6e12 Hz. Although the fits to the full SED do not include significant thermal emission in the X-ray band, previous observations have consistently seen such a component, and we find that there is evidence at the 3.1-sigma level for a disk-blackbody component with a temperature of 150(+30)(-20) eV and an inner radius of 5-14 R_g. If this component is real, it might imply the presence of an inner accretion disk in addition to the strongly truncated (R_in>212 R_g) disk.
We present a model for the interstellar medium of NGC4214 with the objective to probe the physical conditions in the two main star-forming regions and their connection with the star formation activity of the galaxy. We used the spectral synthesis code Cloudy to model an HII region and the associated photodissociation region (PDR) to reproduce the emission of mid- and far-infrared fine-structure cooling lines from the Spitzer and Herschel space telescopes for these two regions. Input parameters of the model, such as elemental abundances and star formation history, are guided by earlier studies of the galaxy, and we investigated the effect of the mode in which star formation takes place (bursty or continuous) on the line emission. Furthermore, we tested the effect of adding pressure support with magnetic fields and turbulence on the line predictions. We find that this model can satisfactorily predict (within a factor of ~2) all observed lines that originate from the ionized medium ([SIV] 10.5um, [NeIII] 15.6um, [SIII] 18.7um, [SIII] 33.5um, and [OIII] 88um), with the exception of [NeII] 12.8um and [NII] 122um, which may arise from a lower ionization medium. In the PDR, the [OI] 63um, [OI] 145um, and [CII] 157um lines are matched within a factor of ~5 and work better when weak pressure support is added to the thermal pressure or when the PDR clouds are placed farther away from the HII regions and have covering factors lower than unity. Our models of the HII region agree with different evolutionary stages found in previous studies, with a more evolved, diffuse central region, and a younger, more compact southern region. However, the local PDR conditions are averaged out on the 175 pc scales that we probe and do not reflect differences observed in the star formation properties of the two regions.
The mass function of hydrogen-rich atmosphere white dwarfs has been frequently found to reveal a distinctive high-mass excess near 1Msun. However, a significant excess of massive white dwarfs has not been detected in the mass function of the largest white dwarf catalogue to date from the Sloan Digital Sky Survey. Hence, whether a high-mass excess exists or not has remained an open question. In this work we build the mass function of the latest catalogue of data release 10 SDSS hydrogen-rich white dwarfs, including the cool and faint population (i.e. effective temperatures 6,000 <~ Teff <~ 12,000 K, equivalent to 12 mag <~ Mbol <~ 13 mag). We show that the high-mass excess is clearly present in our mass function, and that it disappears only if the hottest (brightest) white dwarfs (those with Teff >~ 12,000 K, Mbol <~ 12 mag) are considered. This naturally explains why previous SDSS mass functions failed at detecting a significant excess of high-mass white dwarfs. Thus, our results provide additional and robust observational evidence for the existence of a distinctive high-mass excess near 1Msun. We investigate possible origins of this feature and argue that the most plausible scenario that may lead to an observed excess of massive white dwarfs is the merger of the degenerate core of a giant star with a main sequence or a white dwarf companion during or shortly after a common envelope event.
We report results of a high-resolution numerical-relativity simulation for the merger of black hole-magnetized neutron star binaries on Japanese supercomputer "K". We focus on a binary that is subject to tidal disruption and subsequent formation of a massive accretion torus. We find the launch of thermally driven torus wind, subsequent formation of a funnel wall above the torus and a magnetosphere with collimated poloidal magnetic field, and high Blandford-Znajek luminosity. We show for the first time this picture in a self-consistent simulation. The turbulence-like motion induced by the non-axisymmetric magnetorotational instability as well as the Kelvin-Helmholtz instability inside the accretion torus works as an agent to drive the mass accretion and converts the accretion energy to thermal energy, which results in the generation of a strong wind. By an in-depth resolution study, we reveal that high resolution is essential to draw such a picture. We also discuss the implication for the r-process nucleosynthesis, the radioactively-powered transient emission, and short gamma-ray bursts.
GRB 080503 is a short gamma-ray burst (GRB) detected by \emph{Swift} and has been classified as a compact-star-merger-origin GRB. The soft extended emission and the simultaneous late re-brightening in both the X-ray and optical afterglow lightcurves raise interesting questions regarding its physical origin. We show that the broad-band data of GRB 080503 can be well explained within the framework of the double neutron star merger model, provided that the merger remnant is a rapidly-rotating massive neutron star with an extremely high magnetic field (i.e. a millisecond magnetar). We show that the late optical re-brightening is consistent with the emission from a magnetar-powered "merger-nova". This adds one more case to the growing sample of merger-novae associated with short GRBs. The soft extended emission and the late X-ray excess emission are well connected through a magnetar dipole spin-down luminosity evolution function, suggesting that direct magnetic dissipation is the mechanism to produce these X-rays. The X-ray emission initially leaks from a hole in the merger ejecta pierced by the short GRB jet. The hole subsequently closes after the magnetar spins down and the magnetic pressure drops below ram pressure. The X-ray photons are then trapped behind the mergernova ejecta until the ejecta becomes optically thin at a later time. This explains the essentially simultaneous re-brightening in both the optical and X-ray lightcurves. Within this model, future gravitational wave sources could be associated with a bright X-ray counterpart along with the mergernova, even if the short GRB jet beams away from Earth.
We present a novel, observationally-based framework for the formation epochs and sites of globular clusters (GCs) in a cosmological context. Measuring directly the mean ages of the metal-poor and metal-rich GC subpopulations in our own Galaxy, and in other galaxies, is observationally challenging. Here we apply an alternative approach utilising the property that the galaxy mass-metallicity relation is a strong function of redshift (or look-back age) but is relatively insensitive to galaxy mass for massive galaxies. Assuming that GCs follow galaxy mass-metallicity relations that evolve with redshift, one can estimate the mean formation epochs of the two GC subpopulations by knowing their mean metallicities and the growth in host galaxy mass with redshift. Recently, the SLUGGS survey has measured the spectroscopic metallicities for over 1000 GCs in a dozen massive early-type galaxies. Here we use these measurements, and our new metallicity matching method, to infer a mean age for metal-rich GCs of 11.5 Gyr (z = 2.9) and a range of 12.2 to 12.8 Gyr (4.8 < z < 5.9) for the metal- poor GCs, depending on whether they mostly formed in accreted satellites or in-situ within the main host galaxy. We compare our values to direct age measurements for Milky Way GCs and predictions from cosmological models. Our findings suggest that reionisation preceded most GC formation, and that it is unlikely to be the cause of GC bimodal metallicity distributions.
We explore UV and optical variability signatures for several hundred members of NGC 2264 (3 Myr). We performed simultaneous u- and r-band monitoring over two full weeks with CFHT/MegaCam. About 750 young stars are probed; 40% of them are accreting. Statistically distinct variability properties are observed for accreting and non-accreting cluster members. The accretors exhibit a significantly higher level of variability than the non-accretors, especially in the UV. The amount of u-band variability correlates statistically with UV excess in disk-bearing objects, which suggests that accretion and star-disk interaction are the main sources of variability. Cool magnetic spots, several hundred degrees colder than the photosphere and covering from 5 to 30% of the stellar surface, appear to be the leading factor of variability for the non-accreting stars. In contrast, accretion spots, a few thousand degrees hotter than the photosphere and covering a few percent of the stellar surface, best reproduce the variability of accreting objects. The color behavior is also found to be different between accreting and non-accreting stars. Typical variability amplitudes for accreting members rapidly increase from r to u, which indicates a much stronger contrast at short wavelengths; a lower color dependence in the amplitudes is instead measured for diskless stars. We find that u-band variability on hour timescales is typically about 10% of the peak-to-peak variability on day timescales, while longer term (years) variability is consistent with amplitudes measured over weeks. We conclude that for both accreting and non-accreting stars, the mid-term rotational modulation by spots is the leading timescale for a variability of up to several years. In turn, this suggests that the accretion process is essentially stable over years, although it exhibits low-level shorter term variations in single accretion events.
Planets with sizes between those of Earth and Neptune divide into two populations: purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger gas-enveloped planets possessing voluminous and optically thick atmospheres. We show that whether a planet forms rocky or gas-enveloped depends on the solid surface density of its parent disk. Assembly times for rocky cores are sensitive to disk solid surface density. Lower surface densities spawn smaller planetary embryos; to assemble a core of given mass, smaller embryos require more mergers between bodies farther apart and therefore exponentially longer formation times. Gas accretion simulations yield a rule of thumb that a rocky core must be at least 2$M_\oplus$ before it can acquire a volumetrically significant atmosphere from its parent nebula. In disks of low solid surface density, cores of such mass appear only after the gas disk has dissipated, and so remain purely rocky. Higher surface density disks breed massive cores more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We test model predictions against observations, using planet radius as an observational proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface density. Theory can explain the observation that metal-rich stars host predominantly gas-enveloped planets.
We study the resonant decay of the primordial Standard Model Higgs condensate after inflation into $SU(2)$ gauge bosons on the lattice. We find that the non-Abelian interactions between the gauge bosons quickly extend the momentum distribution towards high values, efficiently destroying the condensate after the onset of backreaction. For the inflationary scale $H = 10^8$ GeV, we find that 90% of the Higgs condensate has decayed after $n \sim 10$ oscillation cycles. This differs significantly from the Abelian case where, given the same couplings strengths, most of the condensate would persist after the resonance.
The Cosmology Data Management system (CosmoDM) is an automated and flexible data management system for the processing and calibration of data from optical photometric surveys. It is designed to run on supercomputers and to minimize disk I/O to enable scaling to very high throughput during periods of reprocessing. It serves as an early prototype for one element of the ground-based processing required by the Euclid mission and will also be employed in the preparation of ground based data needed in the eROSITA X-ray all sky survey mission. CosmoDM consists of two main pipelines. The first is the single-epoch or detrending pipeline, which is used to carry out the photometric and astrometric calibration of raw exposures. The second is the co- addition pipeline, which combines the data from individual exposures into deeper coadd images and science ready catalogs. A novel feature of CosmoDM is that it uses a modified stack of As- tromatic software which can read and write tile compressed images. Since 2011, CosmoDM has been used to process data from the DECam, the CFHT MegaCam and the Pan-STARRS cameras. In this paper we shall describe how processed Pan-STARRS data from CosmoDM has been used to optically confirm and measure photometric redshifts of Planck-based Sunyaev-Zeldovich effect selected cluster candidates.
Scintillation noise due to the Earth's turbulent atmosphere can be a dominant noise source in high-precision astronomical photometry when observing bright targets from the ground. Here we describe the phenomenon of scintillation from its physical origins to its effect on photometry. We show that Young's (1967) scintillation-noise approximation used by many astronomers tends to underestimate the median scintillation noise at several major observatories around the world. We show that using median atmospheric optical turbulence profiles, which are now available for most sites, provides a better estimate of the expected scintillation noise and that real-time turbulence profiles can be used to precisely characterise the scintillation noise component of contemporaneous photometric measurements. This will enable a better understanding and calibration of photometric noise sources and the effectiveness of scintillation correction techniques. We also provide new equations for calculating scintillation noise, including for extremely large telescopes where the scintillation noise will actually be lower than previously thought. These equations highlight the fact that scintillation noise and shot noise have the same dependence on exposure time and so if an observation is scintillation limited, it will be scintillation limited for all exposure times. The ratio of scintillation noise to shot noise is also only weakly dependent on telescope diameter and so a bigger telescope may not yield a reduction in fractional scintillation noise.
The density-based clustering algorithm DBSCAN has been applied to the Fermi Large Area Telescope (LAT) dataset of $ E_{\gamma} \geqslant 100$~GeV events with $\lvert b\rvert>10^{\circ}$, in order to search for new very high energy (VHE) $\gamma$-ray sources. The clustering analysis returned 49 clusters, of which 21 correspond to already known VHE-emitting active galactic nuclei (AGN) within the TeVCat catalogue and a further 11 were found to be significant in a full Fermi analysis. Of these, 2 are previously detected Fermi VHE AGN, and 9 represent new VHE sources consisting of 6 BL Lac objects, one blazar of unknown type and 2 unassociated sources. Comparing these, along with the VHE AGN RBS 0679 and RBS 0970 previously detected with Fermi-LAT, to the current populations of AGN detected with ground-based instruments and Fermi suggests that the VHE-emitting AGN discovered in this study are very similar to the TeVCat AGN and therefore further observations with ground-based imaging atmospheric Cherenkov telescopes are recommended.
WISE J104915.57-531906.1 A+B and WISE J085510.83-071442.5 were recently discovered as the third and fourth closest known systems to the Sun, respectively (2.0 and 2.3 pc). The former consists of a L8+T0.5 binary and the latter is a probable Y dwarf and is the coldest known brown dwarf (~250 K). We present a search for common proper motion companions to these brown dwarfs using multi-epoch mid-infrared images from the Spitzer Space Telescope. We have also obtained near-infrared adaptive optics images of WISE J104915.57-531906.1 A+B with the Very Large Telescope to search for companions at smaller separations than reached by Spitzer. No new companions are detected in either system. At projected separations of 25-420" (50-840 AU) for WISE J104915.57-531906.1 A+B and 4-420" (9-970 AU) for WISE J085510.83-071442.5, the Spitzer images are sensitive to companions with M_4.5<21.6 and 21.9, respectively, which correspond to masses of >=1 M_Jup for ages of >=1 Gyr and temperatures of >=150 K. The detection limit in the adaptive optics images of WISE J104915.57-531906.1 A+B is dH~10 at 3-15" (6-30 AU), or >=7 M_Jup for >=1 Gyr.
IRC+10216 is a circumstellar envelope around a carbon-rich evolved star which contains a large variety of molecules. According to interferometric observations, molecules are distributed either concentrated around the central star or as a hollow shell with a radius of 15". We present ALMA Cycle 0 band 6 observations of the J=14-13 rotational transition of CH3CN in IRC+10216, obtained with an angular resolution of 0.76x0.61. The bulk of the emission is distributed as a hollow shell located at just 2" from the star, with a void of emission in the central region up to a radius of 1". This spatial distribution is markedly different from those found to date in this source for other molecules. Our analysis indicate that methyl cyanide is not formed neither in the stellar photosphere nor far in the outer envelope, but at radial distances as short as 1-2", reaching a maximum abundance of 5e-8 with respect to H2 at about 6" from the star. Standard chemical models of IRC+10216 predict that CH3CN should form farther out in the envelope, at a radius of 15", as other species such as polyyne radicals and cyanopolyynes. We discuss possible mechanims able to bring the region of formation of CH3CN to shorter radii.
We describe the SPIDER flight cryostat, which is designed to cool six millimeter-wavelength telescopes during an Antarctic long-duration balloon flight. The cryostat, one of the largest to have flown on a stratospheric payload, uses liquid helium-4 to deliver cooling power to stages at 4.2 and 1.6 K. Stainless steel capillaries facilitate a high flow impedance connection between the main liquid helium tank and a smaller superfluid tank, allowing the latter to operate at 1.6 K as long as there is liquid in the 4.2 K main tank. Each telescope houses a closed cycle helium-3 adsorption refrigerator that further cools the focal planes down to 300 mK. Liquid helium vapor from the main tank is routed through heat exchangers that cool radiation shields, providing negative thermal feedback. The system performed successfully during a 17 day flight in the 2014-2015 Antarctic summer. The cryostat had a total hold time of 16.8 days, with 15.9 days occurring during flight.
We present our preliminary results on the asteroseismological investigations of the ZZ Ceti star GD 244. We used literature values of the effective temperature and surface gravity and utilized the White Dwarf Evolution Code of Bischoff-Kim, Montgomery and Winget (2008, ApJ, 675, 1512) to build our model grid for the seismological analysis. Five observed pulsational modes published up to now were used to find acceptable model solutions. We found that the best model fits have masses between 0.61 and 0.74 M_solar and constitute two groups with hydrogen layer masses of either ~10^{-5} or 10^{-6} M_solar. Based on a statistical analysis of a larger sample of possible model solutions, we assume that the mass of the star is below ~0.68 M_solar and the oxygen content in the centre is less than 60 percent.
Previous studies of GD 244 revealed seven pulsation frequencies (two doublets and three single periods) in the light variations of the star. The data obtained at McDonald Observatory between 2003 and 2006, and our additional measurements in 2006 and 2007 at Konkoly Observatory, allow the investigation of the long-term pulsational behaviour of GD 244. We found that the 307.1 s period component of one of the doublets show long-term, periodic amplitude modulation with a time scale of ~740 days. Possible explanations are that nonlinear resonant mode coupling is operating among the rotationally split frequency components, or two modes, unresolved in the yearly data are excited at ~307.1 s. This is the first time that such long-term periodic amplitude modulation is published on a ZZ Ceti star.
Using the Hinode/EIS 2$"$ spectroscopic observations, we study the intensity, velocity, and FWHM variations of the strongest Fe XII 195.12 \AA\ line along the jet to find the signature of Alfv\'en waves. We simulate numerically the impulsively generated Alfv\'en waves within the vertical Harris current-sheet, forming the jet plasma flows, and mimicking their observational signatures. Using the FLASH code and the atmospheric model with embedded weakly expanding magnetic field configuration within a vertical Harris current-sheet, we solve the two and half-dimensional (2.5-D) ideal magnetohydrodynamic (MHD) equations to study the evolution of Alfv\'en waves and vertical flows forming the plasma jet. At a height of $\sim 5~\mathrm{Mm}$ from the base of the jet, the red-shifted velocity component of Fe XII 195.12 \AA\ line attains its maximum ($5~\mathrm{km\,s}^{-1}$) which converts into a blue-shifted one between the altitude of $5-10~\mathrm{Mm}$. The spectral intensity continously increases up to $10~\mathrm{Mm}$, while FWHM still exhibits the low values with almost constant trend. This indicates that the reconnection point within the jet's magnetic field topology lies in the corona $5-10~\mathrm{Mm}$ from its footpoint anchored in the Sun's surface. Beyond this height, FWHM shows a growing trend. This may be the signature of Alfv\'en waves that impulsively evolve due to reconnection and propagate along the jet. From our numerical data, we evaluate space- and time- averaged Alfv\'en waves velocity amplitudes at different heights in the jet's current-sheet, which contribute to the non-thermal motions and spectral line broadening. The synthetic width of Fe XII $195.12~\mathrm{\AA}$ line exhibits similar trend of increment as in the observational data, possibly proving the existence of impulsively generated (by reconnection) Alfv\'en waves which propagate along the jet.
The first Fermi Large Area Telescope (LAT) catalogue of sources (1FHL) emitting at high energies (above 10 GeV) reports the details of 514 objects detected in the first three years of the Fermi mission. Of these, 71 were reported as unidentified in the 1FHL catalogue, although six are likely to be associated with a supernova remnant (SNR), a Pulsar Wind Nebula (PWN) or a combination of both, thereby leaving a list of 65 still unassociated objects. Herein, we report a preliminary analysis on this sample of objects concentrating on nine 1FHL sources, which were found to have a clear optical extragalactic classification. They are all blazar, eight BL Lac and one flat spectrum radio quasar, typically at redshift greater than 0.1.
We compare the performance of several dust models in reproducing the dust
spectral energy distribution (SED) per unit extinction in the diffuse
interstellar medium (ISM). We use our results to constrain the variability of
the optical properties of big grains in the diffuse ISM, as published by the
Planck collaboration.
We use two different techniques to compare the predictions of dust models to
data from the Planck HFI, IRAS and SDSS surveys. First, we fit the far-infrared
emission spectrum to recover the dust extinction and the intensity of the
interstellar radiation field (ISRF). Second, we infer the ISRF intensity from
the total power emitted by dust per unit extinction, and then predict the
emission spectrum. In both cases, we test the ability of the models to
reproduce dust emission and extinction at the same time.
We identify two issues. Not all models can reproduce the average dust
emission per unit extinction: there are differences of up to a factor $\sim2$
between models, and the best accord between model and observation is obtained
with the more emissive grains derived from recent laboratory data on silicates
and amorphous carbons. All models fail to reproduce the variations in the
emission per unit extinction if the only variable parameter is the ISRF
intensity: this confirms that the optical properties of dust are indeed
variable in the diffuse ISM.
Diffuse ISM observations are consistent with a scenario where both ISRF
intensity and dust optical properties vary. The ratio of the far-infrared
opacity to the $V$ band extinction cross-section presents variations of the
order of $\sim20\%$ ($40-50\%$ in extreme cases), while ISRF intensity varies
by $\sim30\%$ ($\sim60\%$ in extreme cases). This must be accounted for in
future modelling.
I describe a method for synthesizing photometric passbands for use with current and future X-ray instruments. The method permits the standardisation of X-ray passbands and thus X-ray photometry between different instruments and missions. The method is illustrated by synthesizing a passband in the XMM-Newton EPIC pn which is similar to the ROSAT PSPC 0.5-2 keV band.
The existence of planets orbiting a central binary star system immediately
raises questions regarding their formation and dynamical evolution. Recent
discoveries of circumbinary planets by the Kepler space telescope has shown
that some of these planets reside close to the dynamical stability limit where
it is very difficult to form planets in situ. For binary systems with nearly
circular orbits, such as Kepler-38, the observed proximity of planetary orbits
to the stability limit can be understood by an evolutionary process in which
planets form farther out in the disk and migrate inward to their observed
position. The Kepler-34 system has a high orbital eccentricity of 0.52. Here,
we analyse evolutionary scenarios for the planet observed around this system
using two-dimensional hydrodynamical simulations.
The highly eccentric binary opens a wide inner hole in the disk which is also
eccentric, and displays a slow prograde precession. As a result of the large,
eccentric inner gap, an embedded planet settles in a final equilibrium position
that lies beyond the observed location of Kepler-34 b, but has the correct
eccentricity. In this configuration the planetary orbit is aligned with the
disk in a state of apsidal corotation.To account for the closer orbit of
Kepler-34 b to the central binary, we considered a two-planet scenario and
examined the evolution of the system through joint inward migration and capture
into mean-motion resonances. When the inner planet orbits inside the gap of the
disk, planet-planet scattering ensues. While often one object is thrown into a
large, highly eccentric orbit, at times the system is left with a planet close
to the observed orbit, suggesting that Kepler 34 might have had two
circumbinary planets where one might have been scattered out of the system or
into an orbit where it did not transit the central binary during the operation
of Kepler.
We investigate the pc-scale kinematics and kpc-scale radio morphology of the quasar PG 1302-102, which may harbour a sub-pc separation supermassive binary black hole system at its centre as inferred from optical variability. High-resolution radio interferometric measurements obtained with the Very Long Baseline Array (VLBA) in the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) programme at 15 GHz at 20 epochs spanning 17 years were analysed to investigate the pc-scale radio structure. Archival observations with the Very Large Array (VLA) at 1.4 GHz and 5 GHz were obtained to study the kpc-scale morphology. We find that the pc-scale jet is inclined within ~2 deg to the line of sight and has an apparent opening angle of about 10 deg. The parameters derived from the pc-scale radio jet are broadly consistent with those obtained from the analysis of the optical light curve of PG 1302-102. We obtain at least 0.08 for the mass ratio of the two black holes in the system. We find some indication for a helical jet structure on kpc-scale, but the directions of the inner and the extended radio jets are significantly different, obstructing a straightforward connection of the pc- and kpc-scale jets within the binary scenario.
We present a new analysis of ${\it Spitzer}$ observations of comet 29P/Schwassmann-Wachmann 1 taken on UT 2003 November 21, 23, and 24, similar to a previous investigation of the observations (Stansberry et al. 2004), but using the most recent ${\it Spitzer}$ data pipeline products and intensive image processing techniques. Analysis of images from the IRAC 5.8 & 8.0 $\mu$m bands and the MIPS 24.0 & 70.0 $\mu$m bands resulted in photometry measurements of the nucleus after a suite of coma modeling and removal processes were implemented. SW1 was not identified in the 5.8 $\mu$m image from the previous work so its incorporation into this analysis is entirely new. Using the Near Earth Asteroid Thermal Model (Harris 1998) resulted in a nucleus radius measurement of $R$ = 30.2 $^{+3.7}_{-2.9}$ km and an infrared beaming parameter value of $\eta = 0.99$ $^{+0.26}_{-0.19}$. We also measured an infrared geometric albedo, $p_{5.8}$ = 0.5 $\pm$ 0.5. Extrapolating a 0.04 V-band albedo and using a normalized reflectivity gradient $S' = 14.94 \pm 1.09$ [% (1000 \r{A})$^{-1}$] (Duffard et al. 2014) we recover an infrared albedo of $p_{5.8}$ = 0.31 in the near infrared consistent with the value recovered from thermal modeling. The dust composition extracted from IRS spectra are very comet-like, containing mainly amorphous ferromagnesian silicates (but with a minority of crystalline silicates as well), water ice, and metal sulfides.
We investigate the location of an ultra-hard X-ray selected sample of AGN from the Swift Burst Alert Telescope (BAT) catalog with respect to the main sequence (MS) of star-forming galaxies using Herschel-based measurements of the star formation rate (SFR) and stellar mass (\mstar) from Sloan Digital Sky Survey (SDSS) photometry where the AGN contribution has been carefully removed. We construct the MS with galaxies from the Herschel Reference Survey and Herschel Stripe 82 Survey using the exact same methods to measure the SFR and \mstar{} as the Swift/BAT AGN. We find a large fraction of the Swift/BAT AGN lie below the MS indicating decreased specific SFR (sSFR) compared to non-AGN galaxies. The Swift/BAT AGN are then compared to a high-mass galaxy sample (COLD GASS), where we find a similarity between the AGN in COLD GASS and the Swift/BAT AGN. Both samples of AGN lie firmly between star-forming galaxies on the MS and quiescent galaxies far below the MS. However, we find no relationship between the X-ray luminosity and distance from the MS. While the morphological distribution of the BAT AGN is more similar to star-forming galaxies, the sSFR of each morphology is more similar to the COLD GASS AGN. The merger fraction in the BAT AGN is much higher than the COLD GASS AGN and star-forming galaxies and is related to distance from the MS. These results support a model in which bright AGN tend to be in high mass star-forming galaxies in the process of quenching which eventually starves the supermassive black hole itself.
Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J=5-4 and J=10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6+/-4.4)e10 cm-2 in TMC-1 and (3.9+/-1.8)e10 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10)e-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of 1e-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10)e-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.
We present new multi-band transit photometry of three small (R$_{p}$ < 6 R$_{\oplus}$), short-period (P < 6 days) Kepler planet candidates acquired with the Gran Telescopio Canarias. These observations supplement the results presented in Col\'on & Ford (2011) and Col\'on et al. (2012), where we used multicolor transit photometry of five Kepler planet candidates to search for wavelength-dependent transit depths and either validate planet candidates or identify eclipsing binary false positives within our sample. In those previous studies, we provided evidence that three targets were false positives and two targets were planets. Here, we present observations that provide evidence supporting a planetary nature for KOI 439.01 and KOI 732.01, and we find that KOI 531.01, a 6 R$_{\oplus}$ planet candidate around an M dwarf, is likely a false positive. We also present a discussion of the purported "sub-Jovian desert" in the orbital period-planet radius plane, which cannot be easily explained by observational bias. Both KOI 439.01 and KOI 732.01 are likely planets located within the so-called desert and should be investigated with further follow-up observations. As only ~30 of the ~3600 currently active Kepler planet candidates are located within the sub-Jovian desert, it will be interesting to see if these candidates also survive the vetting process and fill in the gap in the period-radius plane. Confirming planets in this regime will be important for understanding planetary migration and evolution processes, and we urge additional follow-up observations of these planet candidates to confirm their nature.
Extrasolar planets that pass in front of their host star (transit) cause a temporary decrease in the apparent brightness of the star once per orbit, providing a direct measure of the planet's size and orbital period. In some systems with multiple transiting planets, the times of the transits are measurably affected by the gravitational interactions between neighbouring planets. In favorable cases, the departures from Keplerian orbits implied by the observed transit times permit planetary masses to be measured, which is key to determining bulk densities. Characterizing rocky planets is particularly difficult, since they are generally smaller and less massive than gaseous planets. Thus, few exoplanets near Earth's size have had their masses measured. Here we report the sizes and masses of three planets orbiting Kepler-138, a star much fainter and cooler than the Sun. We measure the mass of the Mars-sized inner planet based on on the transit times of its neighbour and thereby provide the first density measurement for an exoplanet smaller than Earth. The middle and outer planets are both slightly larger than Earth. The middle planet's density is similar to that of Earth, while the outer planet is less than half as dense, implying that it contains a greater portion of low density components such as H2O and/or H2.
This booklet contains a collection of contributions to the meeting of the JEts and Disks at INAF (JEDI) group, which took place at the Capodimonte Observatory during 9-10 April 2015. Scope of the meeting was to bring together the JEDI researchers of the Italian Istituto Nazionale di Astrofisica (INAF) working in the field of circumstellar disks and jets in young stars, to discuss together the different agents affecting the structure and the evolution of disks, namely accretion, jets and winds. More information on the JEDI group and its activities can be found at \texttt{this http URL}.
(Abridged) The chemical content of the planetary nebula NGC 3918 is investigated through deep, high-resolution UVES at VLT spectrophotometric data. We identify and measure more than 750 emission lines, making ours one of the deepest spectra ever taken for a planetary nebula. Among these lines we detect very faint lines of several neutron-capture elements (Se, Kr, Rb, and Xe), which enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up in the progenitor star of NGC 3918. We find that Kr is strongly enriched in NGC 3918 and that Se is less enriched than Kr, in agreement with the results of previous papers and with predicted s-process nucleosynthesis. We also find that Xe is not as enriched by the s-process in NGC 3918 as is Kr and, therefore, that neutron exposure is typical of modestly sub-solar metallicity AGB stars. A clear correlation is found when representing [Kr/O] vs. log(C/O) for NGC 3918 and other objects with detection of multiple ions of Kr in optical data, confirming that carbon is brought to the surface of AGB stars along with s-processed material during third dredge-up episodes, as predicted by nucleosynthesis models. We also detect numerous refractory element lines (Ca, K, Cr, Mn, Fe, Co, Ni, and Cu). We compute physical conditions from a large number of diagnostics. Thanks to the high ionization of NGC 3918 we detect a large number of recombination lines of multiple ionization stages of C, N, O and Ne. The abundances obtained for these elements by using recently-determined state-of-the-art ICF schemes or simply adding ionic abundances are in very good agreement, demonstrating the quality of the recent ICF scheme for high ionization planetary nebulae.
Using RXTE/ASM archival data, we investigate the behaviour of the spectral hardness ratio as a function of X-ray luminosity in a sample of six transient X-ray pulsars (EXO 2030+375, GX 304-1, 4U 0115+63, V 0332+63, A 0535+26 and MXB 0656-072). In all sources we find that the spectral hardness ratio defined as $F_{5-12\mathrm{keV}}/ F_{1.33-3\mathrm{keV}}$ increases with the ASM flux (1.33--12 keV) at low luminosities and then saturates or even slightly decreases above some critical X-ray luminosity falling into the range $\sim(3-7)\times10^{37}$~erg~s$^{-1}$. Two-dimensional structure of accretion columns in the radiation-diffusion limit is calculated for two possible geometries (filled and hollow cylinder) for mass accretion rates $\dot M$ ranging from $10^{17}$ to 1.2$\times 10^{18}$~g s$^{-1}$. The observed spectral behaviour in the transient X-ray pulsars with increasing $\dot M$ can be reproduced by a Compton saturated sidewall emission from optically thick magnetized accretion columns with taking into account the emission reflected from the neutron star atmosphere. At $\dot M$ above some critical value $\dot M_{cr}\sim (6-8)\times 10^{17}$~g~s$^{-1}$, the hight of the column becomes such that the contribution of the reflected component to the total emission starts decreasing, which leads to the saturation and even slight decrease of the spectral hardness. Hollow-cylinder columns have a smaller height than the filled-cylinder ones, and the contribution of the reflected component in the total emission does not virtually change with $\dot M$ (and hence the hardness of the continuum monotonically increases) up to higher mass accretion rates than $\dot M_{cr}$ for the filled columns.
We present the results of lifetime measurements made using beam-foil techniques on levels of astrophysical interest in Pb II producing lines at 1203.6\AA{} ($6s6p^{2}$ $^{2}D_{3/2}$) and 1433.9\AA{} ($6s^{2}6d$ $^{2}D_{3/2}$). We also report the first detection of the Pb II $\lambda1203$ line in the interstellar medium (ISM) from an analysis of archival spectra acquired by the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). The oscillator strengths derived from our experimental lifetimes for Pb II $\lambda\lambda1203$, $1433$ are generally consistent with recent theoretical results, including our own relativistic calculations. Our analysis of high-resolution HST/STIS spectra helps to confirm the relative strengths of the Pb II $\lambda\lambda 1203$, $1433$ lines. However, the oscillator strength that we obtain for Pb II $\lambda1433$ ($0.321\pm0.034$) is significantly smaller than earlier theoretical values, which have been used to study the abundance of Pb in the ISM. Our revised oscillator strength for $\lambda1433$ yields an increase in the interstellar abundance of Pb of 0.43 dex over determinations based on the value given by Morton, indicating that the depletion of Pb onto interstellar dust grains is less severe than previously thought.
Multiple terrestrial and space-based telescopes have been proposed for detecting and tracking near-Earth objects (NEOs). Detailed simulations of the search performance of these systems have used complex computer codes that are not widely available, which hinders accurate cross- comparison of the proposals and obscures whether they have consistent assumptions. Moreover, some proposed instruments would survey infrared (IR) bands, whereas others would operate in the visible band, and differences among asteroid thermal and visible light models used in the simulations further complicate like-to-like comparisons. I use simple physical principles to estimate basic performance metrics for the ground-based Large Synoptic Survey Telescope and three space-based instruments - Sentinel, NEOCam, and a Cubesat constellation. The performance is measured against two different NEO distributions, the Bottke et al. distribution of general NEOs, and the Veres et al. distribution of earth impacting NEO. The results of the comparison show simplified relative performance metrics, including the expected number of NEOs visible in the search volumes and the initial discovery rates expected for each system. Although these simplified comparisons do not capture all of the details, they give considerable insight into the physical factors limiting performance. Multiple asteroid thermal models are considered, including FRM, NEATM, and a new generalized for of FRM (GFRM). I describe issues with how IR albedo and emissivity have been estimated in previous studies, which may render them inaccurate. A thermal model for tumbling asteroids is also developed and suggests that tumbling asteroids may be surprisingly difficult for IR telescopes to observe.
In this paper we consider the question of observational signatures of a false vacuum decay event in the early universe followed by a period of inflation; in particular, motivated by the string landscape, we consider decays in which the parent vacuum has a smaller number of large dimensions than the current vacuum, which leads to an anisotropic universe. We go beyond previous studies, and examine the effects on the CMB temperature and polarisation power spectra, due to both scalar and tensor modes, and consider not only late-time effects but also the full cosmological perturbation theory at early times. We find that whilst the scalar mode behaves as one would expect, and the effects of anisotropy at early times are sub-dominant to the late-time effects already studied, for the tensor modes in fact the the early-time effects grow with multipole and can become much larger than one would expect, even dominating over the late-time effects. Thus these effects should be included if one is looking for such a signal in the tensor modes.
We demonstrate that a nonthermal distribution of particles described by a kappa distribution can be accurately approximated by a weighted sum of Maxwell-Boltzmann distributions. We apply this method to modeling collision processes in kappa-distribution plasmas, with a particular focus on atomic processes important for solar physics. The relevant collision process rate coefficients are generated by summing appropriately weighted Maxwellian rate coefficients. This method reproduces the rate coefficients for a kappa distribution to an estimated accuracy of better than 5%. This is equal to or better than the accuracy of rate coefficients generated using "reverse engineering" methods, which attempt to extract the needed cross sections from the published Maxwellian rate coefficient data and then reconvolve the extracted cross sections with the desired kappa distribution. Our approach of summing Maxwellian rate coefficients is easy to implement using existing spectral analysis software. Moreover, the weights in the sum of the Maxwell-Boltzmann distribution rate coefficients can be found for any value of the parameter kappa, thereby enabling one to model plasmas with a time-varying kappa. Tabulated Maxwellian fitting parameters are given for specific values of kappa from 1.7 to 100. We also provide polynomial fits to these parameters over this entire range. Several applications of our technique are presented, including the plasma equilibrium charge state distribution (CSD), predicting line ratios, modeling the influence of electron impact multiple ionization on the equilibrium CSD of kappa-distribution plasmas, and calculating the time-varying CSD of plasmas during a solar flare.
We test the statistical isotropy and Gaussianity of the cosmic microwave
background (CMB) anisotropies using observations made by the Planck satellite.
Our results are based mainly on the full Planck mission for temperature, but
also include some polarization measurements.
In particular, we consider the CMB anisotropy maps derived from the
multi-frequency Planck data by several component-separation methods. For the
temperature anisotropies, we find excellent agreement between results based on
these sky maps over both a very large fraction of the sky and a broad range of
angular scales, establishing that potential foreground residuals do not affect
our studies.
Tests of skewness, kurtosis, multi-normality, N-point functions, and
Minkowski functionals indicate consistency with Gaussianity, while a power
deficit at large angular scales is manifested in several ways, for example low
map variance. The results of a peak statistics analysis are consistent with the
expectations of a Gaussian random field. The "Cold Spot" is detected with
several methods, including map kurtosis, peak statistics, and mean temperature
profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting
it with several independent tests, and address the subject of a posteriori
correction. Tests of directionality suggest the presence of angular clustering
from large to small scales, but at a significance that is dependent on the
details of the approach. We perform the first examination of polarization data,
finding the morphology of stacked peaks to be consistent with the expectations
of statistically isotropic simulations. Where they overlap, these results are
consistent with the Planck 2013 analysis based on the nominal mission data and
provide our most thorough view of the statistics of the CMB fluctuations to
date.
We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto \nu^\alpha$, $\alpha$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $\Omega_{r\nu} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $\nu$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $\alpha$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.
The Arecibo L-Band Feed Array Zone of Avoidance (ALFA ZOA) Deep Survey is the deepest and most sensitive blind Hi survey undertaken in the ZOA. ALFA ZOA Deep will cover about 300 square degrees of sky behind the Galactic plane in both the inner (30 deg < l < 75 deg; b < |2 deg|) and outer (175 deg < l < 207 deg; -2 deg < b < +1 deg) Galaxy, using the Arecibo Radio Telescope. First results from the survey have found 61 galaxies within a 15 square degree area centered on l = 192 deg and b = -2 deg. The survey reached its expected sensitivity of rms = 1 mJy at 9 km/s channel resolution, and is shown to be complete above integrated flux, F_HI = 0.5 Jy km/s. The positional accuracy of the survey is 28 arcsec and detections are found out to a recessional velocity of nearly 19,000 km/s. The survey confirms the extent of the Orion and Abell 539 clusters behind the plane of the Milky Way and discovers expansive voids, at 10,000 km/s and 18,000 km/s. 26 detections (43%) have a counterpart in the literature, but only two of these have known redshift. Counterparts are 20% less common beyond v_hel = 10,000 km/s and 33% less common at extinctions higher than AB = 3.5 mag. ALFA ZOA Deep survey is able to probe large scale structure beyond redshifts that even the most modern wide-angle surveys have been able to detect in the Zone of Avoidance at any wavelength.
High-precision radial velocity surveys are exploring the population of low-mass exoplanets orbiting bright stars, allowing to derive their orbital parameters accurately, their occurrence rate as well as the statistical distribution of their properties with the aim of constraining planetary models of formation and evolution. The SOPHIE spectrograph has been continuously improved these last years and thanks to an appropriate correction of systematic instrumental drift is now reaching the 2 m/s precision in radial velocity measurements on all time scales. As part of a dedicated radial velocity survey devoting to search for low mass planets around a sample of 190 bright solar-type stars in the northern hemisphere, we report the detection of a warm Neptune with a minimum mass of 16.1 +- 2.7 Mearth orbiting in 40 +- 0.24 days the solar analog HD164595. We also revised the parameters of the muti-planetary system around HD190360. We discuss this new detection in the context of the upcoming space mission CHEOPS devoted to a transit search of bright stars harboring known exoplanets.
The scalar models with exponential interaction, introduced in arXiv:1506.00987, include theories with $\langle \phi(x)\rangle\neq0$. Here, we first consider the theory obtained by normal ordering the exponential of the integrated potential $\int d^Dx\mu^D \exp(-\alpha\phi)$, rather than of $V(\phi)$ itself. This corresponds to fill-in the vacuum of the free scalar theory coupled to the external source with the scalar modes. Next, we show that such a regularization prescription, that we are able to implement in the path-integral formulation, also cures some classical potentials which may be unbounded below. We focus on $V(\phi)=m^4\big(e^{-\phi/m}-e^{\phi/m}\big)$, whose regularized partition function $$ W_R[J]={}_J\langle 0| :e^{-\int d^4xV(\phi)}:|0\rangle_J $$ leads to the exact result $$ \langle\phi(x)\rangle=2m \ , $$ in agreement with the experimental data. Another test is that, while the $(2N+1)$-point function is non-trivial, the full propagator is the free one, so that $m^2$ also corresponds to the pole of the propagator. Such an investigation suggests a natural way to get the lagrangian of the Standard Model, with a different Higgs lagrangian, that may be tested in future experiments at LHC.
The putative recent indication of an unidentified 3.55 keV X-ray line in certain astrophysical sources is taken as a motivation for an improved theoretical computation of the cosmological abundance of 7.1 keV sterile neutrinos. If the line is interpreted as resulting from the decay of Warm Dark Matter, the mass and mixing angle of the sterile neutrino are known. Our computation then permits for a determination of the lepton asymmetry that is needed for producing the correct abundance via the Shi-Fuller mechanism, as well as for an estimate of the non-equilibrium spectrum of the sterile neutrinos. The latter plays a role in structure formation simulations. Results are presented for different flavour structures of the neutrino Yukawa couplings and for different types of pre-existing lepton asymmetries, accounting properly for the charge neutrality of the plasma and incorporating approximately hadronic contributions.
We study a number of $U(1)_X$ models featuring a Dirac fermion dark matter particle. We perform a comprehensive analysis which includes the study of corrections to the muon magnetic moment, dilepton searches with LHC data, as well as direct and indirect dark matter detection constraints. We consider four different coupling structures, namely $U(1)_{B-L}, U(1)_{d-u}, U(1)_{universal}$, and $U(1)_{10+\bar{5}}$, all motivated by compelling extensions to the standard model. We outline the viable and excluded regions of parameter space using a large set of probes. Our key findings are that (i) the combination of direct detection and collider constraints rule out dark matter particle masses lighter than $\sim 1$ TeV, unless rather suppressed Z'-fermion couplings exist, and that (ii) for several of the models under consideration, collider constraints rule out Z' masses up to $sim 3$ TeV.
The observed value of the Higgs mass indicates an instability of the Higgs scalar at large energy scales, and hence also at large field values. In the context of early universe cosmology, this is often considered to lead to problems. Here we point out that we can use the instability of the Higgs field to generate an Ekpyrotic phase of contraction. In the context of string theory it is possible that at very high energy densities extra states become massless, leading to an S-brane which leads to the transition between a contracting phase in the past and the current expanding phase. Thus, the Higgs field can be used to generate a non-singular bouncing cosmology in which the anisotropy problem of usual bouncing scenarios is mitigated.
We investigate the scalar and the tensor spectral indices of the quadratic inflation model in Eddington-inspired Born-Infeld (EiBI) gravity. We find the EiBI corrections to the spectral indices are of second and first order in the slow-roll approximation for the scalar and the tensor perturbations respectively. This is very promising since the quadratic inflation model in general relativity provides a very nice fit for the spectral indices. Together with the suppression of the tensor-to-scalar ratio EiBI inflation is well along with the observational data.
Scalar fields are among the possible candidates for dark energy. This paper is devoted to the scalar fields from the inert doublet model, where instead of one as in the standard model, two SU(2) Higgs doublets are used. The component fields of one SU(2) doublet ($\phi_1$) act in an identical way to the standard model Higgs while the component fields of the second SU(2) doublet ($\phi_2$) are taken to be the dark energy candidate (which is done by assuming that the phase transition in the field has not yet occurred). It is found that one can arrange for late time acceleration (dark energy) by using an SU(2) Higgs doublet in the inert Higgs doublet model, whose vacuum expectation value is zero, in the quintessential regime.
We present exact spherical black hole solutions in de Rham, Gabadadze and Tolley (dRGT) massive gravity for a generic choice of the parameters in the theory, and also discuss the thermodynamical and phase structure of the black hole in both the grand canonical and canonical ensembles (for charged case). It turns out that the dGRT black hole solutions includes the known solutions to the Einstein field equations, such as, the monopole-de Sitter-Schwarzschild ones with the coefficients for the third and fourth terms in the potential and the graviton mass in massive gravity naturally generates the cosmological constant and the global monopole term. Furthermore, we compute the mass, temperature, and entropy of dGRT black hole solutions and also perform thermodynamical stability. It turns out that the presence of the graviton mass completely changes the black hole thermodynamics, and it can provide the Hawking-Page phase transition which is also true for the obtained charged black holes. Interestingly, the entropy of a black hole is unaffected and still obeys area law. In particular, our results, in the limit $m_g \rightarrow 0$, reduced exactly to \emph{vis-$\grave{a}$-vis} the general relativity results.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
This thesis focuses on the two main applications where free-space optical communication (FSOC) can bring the most significant impact: interplanetary communications and quantum communications. Consequently, the dissertation is structured in two sections. In the first one, a novel proposal is suggested regarding to using Cherenkov telescopes as ground-station receivers. A feasibility study addresses the posibility of using the technology developed for the gamma-ray telescopes that will make up the Cherenkov Telescope Array (CTA) in the implementation of a new kind of ground station. Among the main advantages that these telescopes provide are the much larger apertures needed to overcome the power limitation that ground-based gamma-ray astronomy and deep-space optical communication both have. Also, the large number of big telescopes that will be built for CTA will make it possible to reduce unitary costs by economy-scale production. The second section of the thesis is framed in the field of free-space Quantum Key Distribution (QKD), which has become a new paradigm in the discipline of information security. This technique offers a theoretically-secure way to communicate over an insecure channel since the presence of an eventual eavesdropper can be detected. The main challenge of Free-space QKD is the need to operate both under strong atmospheric turbulence and daylight background noise. To mitigate these effects, a trade-off is usually required when designing the receiver's optics, since a narrow field-of-view improves background noise rejection, but increases turbulence-related losses and a wide field-of-view produces the opposite effect. A correction system for atmospheric turbulence is proposed to overcome both limitations at the same time, and different strategies are analyzed and experimented to carry out the implementation and integration within the QKD system.
We present J' and K' imaging linear polarimetric adaptive optics observations of NGC 1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the torus from a magnetohydrodynamical (MHD) framework. In a 0.5" (30 pc) aperture at K', we find that polarisation arising from the passage of radiation from the inner edge of the torus through magnetically aligned dust grains in the clumps is the dominant polarisation mechanism, with an intrinsic polarisation of 7.0%$\pm$2.2%. This result yields a torus magnetic field strength in the range of 4$-$82 mG through paramagnetic alignment, and 139$^{+11}_{-20}$ mG through the Chandrasekhar-Fermi method. The measured position angle (P.A.) of polarisation at K$'$ is found to be similar to the P.A. of the obscuring dusty component at few parsec scales using infrared interferometric techniques. We show that the constant component of the magnetic field is responsible for the alignment of the dust grains, and aligned with the torus axis onto the plane of the sky. Adopting this magnetic field configuration and the physical conditions of the clumps in the MHD outflow wind model, we estimate a mass outflow rate $\le$0.17 M$_{\odot}$ yr$^{-1}$ at 0.4 pc from the central engine for those clumps showing near-infrared dichroism. The models used were able to create the torus in a timescale of $\geq$10$^{5}$ yr with a rotational velocity of $\leq$1228 km s$^{-1}$ at 0.4 pc. We conclude that the evolution, morphology and kinematics of the torus in NGC 1068 can be explained within a MHD framework.
CR7 is the brightest $z=6.6 \, {\rm Ly}\alpha$ emitter (LAE) known to date, and spectroscopic follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III stars. We examine this interpretation using cosmological hydrodynamical simulations. Several simulated galaxies show the same "Pop III wave" pattern observed in CR7. However, to reproduce the extreme CR7 ${\rm Ly}\alpha$/HeII1640 line luminosities ($L_{\rm \alpha/He II}$) a top-heavy IMF and a massive ($>10^{7}{\rm M}_{\odot}$) PopIII burst with age $<2$ Myr are required. Assuming that the observed properties of ${\rm Ly}\alpha$ and HeII emission are typical for Pop III, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs at $z=6.6$ with $L_{\alpha} >10^{43.3}{\rm erg}\,{\rm s}^{-1}$ should also host Pop III stars producing an observable $L_{\rm He II}>10^{42.7}{\rm erg}\,{\rm s}^{-1}$. As an alternate explanation, we explore the possibility that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). Our model predicts $L_{\alpha}$, $L_{\rm He II}$, and X-ray luminosities that are in agreement with the observations. In any case, the observed properties of CR7 indicate that this galaxy is most likely powered by sources formed from pristine gas. We propose that further X-ray observations can distinguish between the two above scenarios.
The metal budgets in some white dwarf (WD) atmospheres reveal that volatile-rich circumstellar bodies must both exist in extrasolar systems and survive the giant branch phases of stellar evolution. The resulting behaviour of these active asteroids or comets which orbit WDs is not well-understood, but may be be strongly influenced by sublimation due to stellar radiation. Here we develop a model, generally applicable to any extrasolar system with a main sequence or WD star, that traces sublimation-induced orbital element changes in approximately km-sized extrasolar minor planets and comets traveling within hundreds of au. We derive evolution equations on orbital timescales and for arbitrarily steep power-law sublimation dependencies on distance, and place our model in a Solar system context. We also demonstrate the importance of coupling sublimation and general relativity, and the orbital consequences of outgassing in arbitrary directions. We prove that nongravitational accelerations alone cannot result in orbit crossing with the WD disruption radius, but may shrink or expand the orbit by up to several au after a single pericentre passage, potentially affecting subsequent interactions with remnant debris and planets. Our analysis suggests that extant planets must exist in polluted WD systems.
We present NuSTAR observations of the powerful radio galaxy Cygnus A, focusing on the central absorbed active galactic nucleus (AGN). Cygnus A is embedded in a cool-core galaxy cluster, and hence we also examine archival XMM-Newton data to facilitate the decomposition of the spectrum into the AGN and intracluster medium (ICM) components. NuSTAR gives a source-dominated spectrum of the AGN out to >70keV. In gross terms, the NuSTAR spectrum of the AGN has the form of a power law (Gamma~1.6-1.7) absorbed by a neutral column density of N_H~1.6x10^23 cm^-2. However, we also detect curvature in the hard (>10keV) spectrum resulting from reflection by Compton-thick matter out of our line-of-sight to the X-ray source. Compton reflection, possibly from the outer accretion disk or obscuring torus, is required even permitting a high-energy cutoff in the continuum source; the limit on the cutoff energy is E_cut>111keV (90% confidence). Interestingly, the absorbed power-law plus reflection model leaves residuals suggesting the absorption/emission from a fast (15,000-26,000km/s), high column-density (N_W>3x10^23 cm^-2), highly ionized (xi~2,500 erg cm/s) wind. A second, even faster ionized wind component is also suggested by these data. We show that the ionized wind likely carries a significant mass and momentum flux, and may carry sufficient kinetic energy to exercise feedback on the host galaxy. If confirmed, the simultaneous presence of a strong wind and powerful jets in Cygnus A demonstrates that feedback from radio-jets and sub-relativistic winds are not mutually exclusive phases of AGN activity but can occur simultaneously.
It has been suggested that the wiggle instability (WI) of spiral shocks in a galactic disk is responsible for the formation of gaseous feathers observed in grand-design spiral galaxies. We perform both a linear stability analysis and numerical simulations to investigate the effect of magnetic fields on the WI. The disk is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. We control the strengths of magnetic fields and spiral-arm forcing using the dimensionless parameters $\beta$ and $\mathcal{F}$, respectively. By solving the perturbation equations as a boundary-eigenvalue problem, we obtain dispersion relations of the WI for various values of $\beta=1-\infty$ and $\mathcal{F}=5\%$ and $10\%$. We find that the WI arising from the accumulation of potential vorticity at disturbed shocks is suppressed, albeit not completely, by magnetic fields. The stabilizing effect of magnetic fields is not from the perturbed fields but from the unperturbed fields that reduce the density compression factor in the background shocks. When $\mathcal{F}=5\%$ and $\beta\lesssim 10$ or $\mathcal{F}=10\%$ and $\beta\sim5-10$, the most unstable mode has a wavelength of $\sim0.1-0.2$ times the arm-to-arm separation, which appears consistent with a mean spacing of observed feathers.
Supernovae are the most energetic among stellar feedback processes, and are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing high resolution, 3D hydrodynamical simulations at various SN rates, $S$, and ISM average densities, $n$. We find that the evolution of a SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By traveling faster and further in the cooling-inefficient hot phase, the spatial-temporal domain of a SNR is enlarged by $>10^{2.5}$ in a hot-dominated multiphase medium (HDMM) compared to the uniform case. We then examine the resultant ISM as we vary $n$ and $S$, finding that a steady state can only be achieved when the hot gas volume fraction \fvh $\lesssim 0.6\pm 0.1$. Above that, overlapping SNRs render connecting topology of the hot gas, and such a HDMM is subjected to thermal runaway with growing pressure and \fvh. Photoelectric heating (PEH) has a surprisingly strong impact on \fvh. For $n \gtrsim 3 cm^{-3}$, a reasonable PEH rate is able to suppress the ISM from undergoing thermal runaway. Overall, we determine that the critical SN rate for the onset of thermal runaway is roughly $S_{crit} = 200 (n/1cm^{-3})^k (E_{SN}/10^{51} erg)^{-1} kpc^{-3} Myr^{-1}$, where k=(1.2,2.7) for $n$ < 1 and >1 cm$^{-3}$, respectively. We present a fitting formula of the ISM pressure $P(n, S)$, which can be used as an effective equation of state in cosmological simulations. The observed velocities of OB stars imply that the core collapse SN are almost randomly located on scales $\lesssim$ 150 pc. Despite the 5 orders of magnitude span of $(n,S)$, the average Mach number shows very small variations: $M \approx 0.5\pm 0.2, 1.2\pm 0.3, 2.3\pm 0.9$ for the hot, warm and cold phases, respectively.
Recent broad-band 34 and 44 GHz radio continuum observations of the Galactic center have revealed 41 massive stars identified with near-IR counterparts, as well as 44 proplyd candidates within 30" of Sgr A*. Radio observations obtained in 2011 and 2014 have been used to derive proper motions of eight young stars near Sgr A*. The accuracy of proper motion estimates based on near-IR observations by Lu et al. and Paumard et al. have been investigated by using their proper motions to predict the 2014 epoch positions of near-IR stars and comparing the predicted positions with those of radio counterparts in the 2014 radio observations. Predicted positions from Lu et al. show an rms scatter of 6 mas relative to the radio positions, while those from Paumard et al. show rms residuals of 20 mas, which is mainly due to uncertainties in the IR-based proper motions. Under the assumption of homogeneous ionized winds, we also determine the mass-loss rates of 11 radio stars, finding rates that are on average $\sim$2 times smaller than those determined from model atmosphere calculations and near-IR data. Clumpiness of ionized winds would reduce the mass loss rate of WR and O stars by additional factors of 3 and 10, respectively. One important implication of this is a reduction in the expected mass accretion rate onto Sgr A* from stellar winds by nearly an order of magnitude to a value of few$\times10^{-7}$ \msol\ yr$^{-1}$. Finally, we present the positions of 318 compact 34.5 GHz radio sources within 30\arcs\ of Sgr A*. At least 45 of these have stellar counterparts in the near-IR $K_s$ (2.18 $\mu$m) and $L'$ (3.8$\mu$m) bands.
We present unresolved single stellar population synthesis models in the near-infrared (NIR) range. The extension to the NIR is important for the study of early-type galaxies, since these galaxies are predominantly old and therefore emit most of their light in this wavelength range. The models are based on a library of empirical stellar spectra, the NASA infrared telescope facility (IRTF) spectral library. Integrating these spectra along theoretical isochrones, while assuming an initial mass function (IMF), we have produced model spectra of single age-metallicity stellar populations at a resolution R~2000. These models can be used to fit observed spectral of globular clusters and galaxies, to derive their age distribution, chemical abundances and IMF. The models have been tested by comparing them to observed colours of elliptical galaxies and clusters in the Magellanic Clouds. Predicted absorption line indices have been compared to published indices of other elliptical galaxies. The comparisons show that our models are well suited for studying stellar populations in unresolved galaxies. They are particularly useful for studying the old and intermediate-age stellar populations in galaxies, relatively free from contamination of young stars and extinction by dust. These models will be indispensable for the study of the upcoming data from JWST and extremely large telescopes, such as the E-ELT.
We present a detailed study of the stars of the IRTF spectral library to understand its full extent and reliability for use with Stellar Population (SP) modeling. The library consist of 210 stars, with a total of 292 spectra, covering the wavelength range of 0.94 to 2.41 micron at a resolution R = 2000. For every star we infer the effective temperature (Teff), gravity (logg) and metallicity ([Z/Zsun]) using a full-spectrum fitting approach in a section of the K band (2.19 to 2.34 micron) and temperature-NIR colour relations. We test the flux calibration of these stars by calculating their integrated colours and comparing them with the Pickles library colour-temperature relations. We also investigate the NIR colours as a function of the calculated effective temperature and compared them in colour-colour diagrams with the Pickles library. This latter test shows a good broad-band flux calibration, important for the SP models. Finally, we measure the resolution R as a function of wavelength. We find that the resolution increases as a function of lambda from about 6 angstrom in J to 10 angstrom in the red part of the K-band. With these tests we establish that the IRTF library, the largest currently available general library of stars at intermediate resolution in the NIR, is an excellent candidate to be used in stellar population models. We present these models in the next paper of this series.
What intrinsic properties shape the light curves of Type II supernovae (SNe)? To address this question we derive observational measures that are robust (i.e., insensitive to detailed radiative transfer) and constrain the contribution from 56Ni, as well as a combination of the ejecta mass, progenitor radius, and explosion energy. By applying our methods to a sample of type II SNe from the literature we find that 56Ni contribution is often significant. It is typically the source of about 20% of the radiated energy during the photospheric phase, and in extreme cases it even dominates. We find that the 56Ni relative contribution is anti-correlated with the luminosity decline rate, while it is not significantly correlated with other properties of the light curve. When added to other clues, this in turn suggests that the flat plateaus often observed in type II SNe are not a generic feature of the cooling envelope emission, and that without 56Ni many of the SNe that are classified as II-P would have shown a decline rate that is steeper by up to 1 mag/100 d. Furthermore, contrary to previous suggestions, we find that SNe with faster decline have at least comparable, and possibly larger, ejecta mass compared to those with flat plateaus.
We present a new model of the $\beta$ Pictoris disk-and-planet system that simulates both the planetesimal collisions and the dynamics of the resulting dust grains, allowing us to model features and asymmetries in both thermal and scattered light images of the disk. Given the observed inclination and eccentricity of the $\beta$ Pictoris b planet, the model neatly ties together several features of the disk: the central hole in the submillimeter images, the two-disk "x"-pattern seen in scattered light, the "wing-tilt" asymmetry, and possibly even the clumpy gas seen by ALMA. We also find that most of the dust in the $\beta$ Pictoris system is likely produced outside the ring at 60-100 AU. Instead of a birth ring, this disk has a "stirring ring" at 60-100 AU where the high-velocity collisions produced by the secular wave launched by the planet are concentrated. The two-disk x-pattern arises because collisions occur more frequently at the peaks and troughs of the secular wave. The perturbations of the disk in this region create an azimuthally and vertically asymmetric spatial distribution of collisions, which could yield an azimuthal clump of gas without invoking resonances or an additional planet.
The UV environment of a host star affects the photochemistry in the atmosphere, and ultimately the surface UV environment for terrestrial planets and therefore the conditions for the origin and evolution of life. We model the surface UV radiation environment for Earth-sized planets orbiting FGKM stars at the 1AU equivalent distance for Earth through its geological evolution. We explore four different types of atmospheres corresponding to an early Earth atmosphere at 3.9 Gyr ago and three atmospheres covering the rise of oxygen to present day levels at 2.0 Gyr ago, 0.8 Gyr ago and modern Earth (Following Kaltenegger et al. 2007). In addition to calculating the UV flux on the surface of the planet, we model the biologically effective irradiance, using DNA damage as a proxy for biological damage. We find that a pre-biotic Earth (3.9 Gyr ago) orbiting an F0V star receives 6 times the biologically effective radiation as around the early Sun and 3520 times the modern Earth-Sun levels. A pre-biotic Earth orbiting GJ 581 (M3.5V) receives 300 times less biologically effective radiation, about 2 times modern Earth-Sun levels. The UV fluxes calculated here provide a grid of model UV environments during the evolution of an Earth-like planet orbiting a range of stars. These models can be used as inputs into photo-biological experiments and for pre-biotic chemistry and early life evolution experiments.
We model the atmospheres and spectra of Earth-like planets orbiting the
entire grid of M dwarfs for active and inactive stellar models with $T_{eff}$ =
2300K to $T_{eff}$ = 3800K and for six observed MUSCLES M dwarfs with UV
radiation data. We set the Earth-like planets at the 1AU equivalent distance
and show spectra from the VIS to IR (0.4$\mu$m - 20$\mu$m) to compare
detectability of features in different wavelength ranges with JWST and other
future ground- and spaced-based missions to characterize exo-Earths. We focus
on the effect of UV activity levels on detectable atmospheric features that
indicate habitability on Earth, namely: H$_2$O, O$_3$, CH$_4$, N$_2$O and
CH$_3$Cl.
To observe signatures of life - O$_2$/O$_3$ in combination with reducing
species like CH$_4$, we find that early and active M dwarfs are the best
targets of the M star grid for future telescopes. The O$_2$ spectral feature at
0.76$\mu$m is increasingly difficult to detect in reflected light of later M
dwarfs due to low stellar flux in that wavelength region. N$_2$O, another
biosignature detectable in the IR, builds up to observable concentrations in
our planetary models around M dwarfs with low UV flux. CH$_3$Cl could become
detectable, depending on the depth of the overlapping N$_2$O feature.
We present a spectral database of Earth-like planets around cool stars for
directly imaged planets as a framework for interpreting future lightcurves,
direct imaging, and secondary eclipse measurements of the atmospheres of
terrestrial planets in the HZ to design and assess future telescope
capabilities.
In magnetic cataclysmic variables (CVs), X-ray emission regions are located close to the white dwarf surface, which is expected to reflect a significant fraction of intrinsic X-rays above 10 keV, producing a Compton reflection hump. However, up to now, a secure detection of this effect in magnetic CVs has largely proved elusive because of the limited sensitivity of non-imaging X-ray detectors. Here we report our analysis of joint NuSTAR/XMM-Newton observations of three magnetic CVs, V709 Cas, NY Lup, and V1223 Sgr. The improved hard X-ray sensitivity of the imaging NuSTAR data has resulted in the first robust detection of Compton hump in all three objects, with amplitudes of ~1 or greater in NY Lup, and likely <1.0 in the other two. We also confirm earlier report of a strong spin modulation above 10 keV in V709 Cas, and report the first detection of small spin amplitudes in the others. We interpret this as due to different height of the X-ray emitting region among these objects. A height of ~0.2 white dwarf radii provides a plausible explanation for the low reflection amplitude of V709 Cas. Since emission regions above both poles are visible at certain spin phases, this can also explain the strong hard X-ray spin modulation. A shock height of ~0.05 white dwarf radii can explain our results on V1223 Sgr, while the shock height in NY Lup appears negligible.
In recent years, a lot of work was done that has revealed some very interesting properties of neutron stars. One can relate the first multipole moments of a neutron star, or quantities that can be derived from them, with relations that are independent of the equation of state. This is a very significant result that has great implications for the description of neutron stars and in particular for the description of the spacetime around them. Additionally, it was recently shown that there is a four parameter analytic spacetime, known as two-soliton, which can accurately capture the properties of the geometry around neutron stars. This allows for the possibility of describing in a unified formalism the astrophysically relevant properties of the spacetime around a neutron star independently of the particulars of the equation of state for the matter of the star. In the present work we investigate properties such as the location of the innermost stable circular orbit $R_{ISCO}$, the various frequencies of perturbed circular equatorial geodesics, the efficiency of an accretion disk, its temperature distribution, and other properties associated to the emitted radiation from the disc, in a way independent of any choice of an equation of state for the neutron star. Furthermore, we provide proof of principle that if one were to measure the right combinations of pairs of these properties, with the additional knowledge of the mass of the neutron star, one could determine the equation of state of the star.
We present four new planetary nebulae (PNe) discovered in the Small Magellanic Cloud (SMC) from deep UK Schmidt telescope (UKST) narrow band H-alpha and broad-band short-red "SR" continuum images and confirmed spectroscopically. All new PNe show strong [NII]/H-alpha ratios in their spectra. We describe and detail the process of PN candidate selection based on wide-field multi-wavelength imaging of the SMC and our subsequent spectroscopic confirmation and classification. We carefully reviewed archived information and available imagery for previous SMC PN detections and various other types of emission objects in the SMC as a training set to help identify new PN candidates. These 4 preliminary discoveries provide a 4% increase to the previously known SMC PN population of ~100. Once spectroscopic follow-up of all our newly identified SMC PN candidates is complete, we expect to increase the total number of known SMC PNe by up to 50%. This will permit a significant improvement to determination of the SMC PN luminosity function and enable further insights into the chemical evolution and kinematics of the SMC PN population.
The formation process of circumstellar disks is still controversial because of the interplay of complex physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we investigate the effect of the Hall current term on the formation of circumstellar disk using three-dimensional simulations. In our simulations, all non-ideal effects as well as the radiation transfer are considered. We show that the size of the disk is significantly affected by a simple difference in the inherent properties of the prestellar core, namely whether the rotation vector and the magnetic field are parallel or anti-parallel. In the former case, only a very small disk ($< 1$ AU) is formed. On the other hand, in the latter case, a massive and large ($>20$ AU) disk is formed in the early phase of protostar formation. We also show that the anti-rotating envelopes against the disk-rotation appear with a size of $\gtrsim 200$ AU. We predict that the anti-rotating envelope will be found in the future observations.
We present the results of a multi-filter photometric monitoring of the wide binary HIP10680/HIP10679. We found both component to be variable with amplitude up to $\Delta$V = 0.03 mag in the case of HIP10680 and $\Delta$V = 0.07 mag in the case of HIP10679. We could measure the rotation periods P = 0.2396d of the hotter F5V component HIP10680 and P = 0.777d of the cooler G5V component HIP10679. We found that the rotation axes of both components are aligned with an inclination $i$ = $\sim$ 10$^{\circ}$. Although the two components have a mass difference not larger than 15%, they exhibit a significant difference between their rotation periods. Such difference may arise either from different initial rotation periods or to different disc life times. For instance, the slower rotating component HIP 10679 hosts a well know debris disc.
We model the inner knot of the Crab Nebula as a synchrotron emission coming from the non-spherical MHD termination shock of relativistic pulsar wind. The post-shock flow is mildly relativistic; as a result the Doppler-beaming has a strong impact on the shock appearance. The model can reproduce the knot location, size, elongation, brightness distribution, luminosity and polarization provided the effective magnetization of the section of the pulsar wind producing the knot is low, $\sigma \leq 1$. In the striped wind model, this implies that the striped zone is rather wide, with the magnetic inclination angle of the Crab pulsar $\ge 45^\circ$; this agrees with the previous model-dependent estimate based on the gamma-ray emission of the pulsar. We conclude that the tiny knot is indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock and that this shock is a site of efficient particle acceleration. On the other hand, the deduced low magnetization of the knot plasma implies that this is an unlikely site for the Crab's gamma-ray flares, if they are related to the fast relativistic magnetic reconnection events.
The distribution of warm molecular gas (1000--3000 K), traced by the near-IR H$_2$ 2.12 $\mu$m line, has been imaged with a resolution $<0.5$ arcsec in the central 1 kpc of seven nearby Seyfert galaxies. We find that this gas is highly concentrated towards the central 100 pc and that its morphology is often symmetrical. Lanes of warm H$_2$ gas are observed only in three cases (NGC\,1068, NGC\,1386 and Circinus) for which the morphology is much wider and extended than the dust filaments. We conclude that there is no one-to-one correlation between dust and warm gas. This indicates that, if the dust filaments and lanes of warm gas are radial streaming motions of fueling material, they must represent \textit{two different phases of accretion}: the dust filaments represent a colder phase than the gas close to the nucleus (within $\sim$100 pc). We predict that the morphology of the nuclear dust at these scales should resemble that of the cold molecular gas (e.g. CO at 10-40 K), as we show for CenA and NGC\,1566 by ALMA observations, whereas the inner H$_2$ gas traces a much warmer phase of material identified with warmer (40-500 K) molecular gas such as CO(6-5) or HCN (as shown by ALMA for NGC\,1068 and NGC\,1097). We also find that X-ray heating is the most likely dominant excitation mechanism of the H$_{2}$ gas for most sources.
To date, the radial velocity (RV) method has been one of the most productive techniques for detecting and confirming extrasolar planetary candidates. Unfortunately, stellar activity can induce RV variations which can drown out or even mimic planetary signals - and it is notoriously difficult to model and thus mitigate the effects of these activity-induced nuisance signals. This is expected to be a major obstacle to using next-generation spectrographs to detect lower mass planets, planets with longer periods, and planets around more active stars. Enter Gaussian processes (GPs) which, we note, have a number of attractive features that make them very well suited to disentangling stellar activity signals from planetary signals. We present here a GP framework we developed to model RV time series jointly with ancillary activity indicators (e.g. bisector velocity spans, line widths, chromospheric activity indices), allowing the activity component of RV time series to be constrained and disentangled from e.g. planetary components. We discuss the mathematical details of our GP framework, and present results illustrating its encouraging performance on both synthetic and real RV datasets, including the publicly-available Alpha Centauri B dataset.
We argue that the fine tuning problems of scalar-driven inflation may be worse than is commonly believed. The reason is that reheating requires the inflaton to be coupled to other matter fields whose vacuum fluctuations alter the inflaton potential. The usual response has been that even more fine-tuning of the classical potential $V(\varphi)$ can repair any damage done in this way. We point out that the effective potential in de Sitter background actually depends in a complicated way upon the dimensionless combination of $\varphi/H$. We also show that the factors of $H$ which occur in de Sitter do not even correspond to local functionals of the metric for general geometries, nor are they Planck-suppressed.
Two classes of GRBs, short and long, have been determined without any doubts, and are usually prescribed to different physical scenarios. A third class, intermediate in $T_{90}$ durations, has been reported to be present in the datasets of BATSE, Swift, RHESSI and possibly BeppoSAX. The latest release of $>1500$ GRBs observed by Fermi gives an opportunity to further investigate the duration distribution. The aim of this paper is to investigate whether a third class is present in the $\log T_{90}$ distribution, or is it described by a bimodal distribution. A standard $\chi^2$ fitting of a mixture of Gaussians is applied to 25 histograms with different binnings. Different binnings give various values of the fitting parameters, as well as the shape of the fitted curve. Among five statistically significant fits none is trimodal. Locations of the Gaussian components are in agreement with previous works. However, a trimodal distribution, understood in the sense of having three separated peaks, is not found for any binning. It is concluded that the duration distribution in Fermi data is well described by a mixture of three log-normal distributions, but it is intrinsically bimodal, hence no third class is present in the $T_{90}$ data of Fermi. It is suggested that the log-normal fit may not be an adequate model.
For the very first time we present 3D simulations of planets embedded in stellar irradiated discs. It is well known that thermal effects could reverse the direction of planetary migration from inwards to outwards, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. When considering stellar irradiation in addition to viscous friction as a source of heating, the outer disc changes from a shadowed to a flared structure. Using a suited analytical formula it has been shown that in the flared part of the disc the migration is inwards; planets can migrate outwards only in shadowed regions of the disc, { because the radial gradient of entropy is stronger there}. In order to confirm this result numerically, we have computed the total torque acting on planets held on fixed orbits embedded in stellar irradiated 3D discs using the hydrodynamical code FARGOCA. We find qualitatively good agreement between the total torque obtained with numerical simulations and the one predicted by the analytical formula. For large masses (>20 Earth masses) we find quantitative agreement, and we obtain outwards migration regions for planets up to 60 Earth masses in the early stages of accretional discs. We find nevertheless that the agreement with the analytic formula is quite fortuitous because the formula underestimates the size of the horseshoe region and therefore overestimates the amount of saturation of the corotation torque; this error is compensated by imperfect estimates of other terms, most likely for the cooling rate.
The sensitivity of Imaging Air Cherenkov Telescopes (IACTs) worsens significantly at low energies because the gamma/hadron separation becomes much more complex. In this paper we study the impact of the single electromagnetic subcascade events on the efficiency of the gamma/hadron separation for a system of four IACTs using Monte Carlo simulations. The studies are done for: two different altitudes of the observatory, three different telescope sizes and two hadron interaction models (GHEISHA and FLUKA). More than 90% of the single electromagnetic proton-induced subcascade events are showers with primary energy below 200 GeV, regardless on the trigger threshold. The estimated efficiency of the gamma/hadron separation using the FLUKA model is similar to results obtained using the GHEISHA model. Nevertheless, for at least one triggered telescope only, a higher fraction of single electromagnetic subcascade events was obtained from the FLUKA model. Finally, the calculated quality factors are anti-correlated with the contributions of the false gamma-ray events in the proton initiated showers. Therefore, the occurrence of single electromagnetic subcascade events is one of the main reasons of the worsening of the primary gamma-ray selection efficiency at low energies.
The ANTARES telescope is well-suited for detecting astrophysical transient neutrino sources as it can observe a full hemisphere of the sky at all times with a high duty cycle. The background due to atmospheric particles can be drastically reduced, and the point-source sensitivity improved, by selecting a narrow time window around possible neutrino production periods. Blazars, being radio-loud active galactic nuclei with their jets pointing almost directly towards the observer, are particularly attractive potential neutrino point sources, since they are among the most likely sources of the very high-energy cosmic rays. Neutrinos and gamma rays may be produced in hadronic interactions with the surrounding medium. Moreover, blazars generally show high time variability in their light curves at different wavelengths and on various time scales. This paper presents a time-dependent analysis applied to a selection of flaring gamma-ray blazars observed by the FERMI/LAT experiment and by TeV Cherenkov telescopes using five years of ANTARES data taken from 2008 to 2012. The results are compatible with fluctuations of the background. Upper limits on the neutrino fluence have been produced and compared to the measured gamma-ray spectral energy distribution.
Subject of this paper are intrinsic ellipticity correlations between galaxies, their statistical properties, their observability with future surveys and their interference with weak gravitational lensing measurements. Using an angular momentum-based, quadratic intrinsic alignment model we derive correlation functions of the ellipticity components and project them to yield the four non-zero angular ellipticity spectra $C^\epsilon_E(\ell)$, $C^\epsilon_B(\ell)$, $C^\epsilon_C(\ell)$ and $C^\epsilon_S(\ell)$ in their generalisation to tomographic surveys. For a Euclid-like survey, these spectra would have amplitudes smaller than the weak lensing effect on nonlinear structures, but would constitute an important systematic. Computing estimation biases for cosmological parameters derived from an alignment-contaminated survey suggests biases of $+5\sigma_w$ for the dark energy equation of state parameter $w$, $-20\sigma_{\Omega_m}$ for the matter density $\Omega_m$ and $-12\sigma_{\sigma_8}$ for the spectrum normalisation $\sigma_8$. Intrinsic alignments yield a signal which is easily observable with a survey similar to Euclid: While not independent, significances for estimates of each of the four spectra reach values of tens of $\sigma$ if weak lensing and shape noise are considered as noise sources, which suggests relative uncertainties on alignment parameters at the percent level.
We report the results from our analysis of a large set of archival data acquired with the X-ray telescope (XRT) onboard Swift, covering the sky region surrounding objects from the first Fermi Large Area Telescope (LAT) catalogue of high-energy sources (1FHL), which still lack an association. Of the 23 regions analysed, ten did not show any evidence of X-ray emission, but 13 were characterised by the presence of one or more objects emitting in the 0.3-10 keV band. Only in a couple of cases is the X-ray counterpart located outside the Fermi positional uncertainty, while in all other cases the associations found are compatible with the high-energy error ellipses. All counterparts we found have been studied in detail by means of a multi-waveband approach to evaluate their nature or class; in most cases, we have been able to propose a likely or possible association except for one Fermi source whose nature remains doubtful at the moment. The majority of the likely associations are extragalactic in nature, most probably blazars of the BL Lac type.
We study the power spectrum of super-Hubble fluctuations of an inflaton-like scalar field, the "system", coupled to another scalar field, the "environment" during de Sitter inflation. We obtain the reduced density matrix for the inflaton fluctuations by integrating out the environmental degrees of freedom. These are considered to be massless and conformally coupled to gravity as a \emph{proxy} to describe degrees of freedom that remain sub-Hubble all throughout inflation. The time evolution of the density matrix is described by a quantum master equation, which describes the decay of the vacuum state, the production of particles and correlated pairs and quantum entanglement between super and sub-Hubble degrees of freedom. The quantum master equation provides a non-perturbative resummation of secular terms from self-energy (loop) corrections to the inflaton fluctuations. In the case studied here these are Sudakov-type double logarithms which result in the \emph{decay} of the power spectrum of inflaton fluctuations upon horizon crossing with a concomitant violation of scale invariance. The reduced density matrix and its quantum master equation furnish a powerful non-perturbative framework to study the effective field theory of long wavelength fluctuations by tracing short wavelength degrees of freedom.
Zernike polynomials are commonly used to represent the wavefront phase on circular optical apertures, since they form a complete and orthonormal basis on the unit disk. In [Diaz et all, 2014] we introduced a new Zernike basis for elliptic and annular optical apertures based on an appropriate diffeomorphism between the unit disk and the ellipse and the annulus. Here, we present a generalization of this Zernike basis for a variety of important optical apertures, paying special attention to polygons and the polygonal facets present in segmented mirror telescopes. On the contrary to ad hoc solutions, most of them based on the Gram-Smith orthonormalization method, here we consider a piece-wise diffeomorphism that transforms the unit disk into the polygon under consideration. We use this mapping to define a Zernike-like orthonormal system over the polygon. We also consider ensembles of polygonal facets that are essential in the design of segmented mirror telescopes. This generalization, based on in-plane warping of the basis functions, provides a unique solution, and what is more important, it guarantees a reasonable level of invariance of the mathematical properties and the physical meaning of the initial basis functions. Both, the general form and the explicit expressions for a typical example of telescope optical aperture are provided.
We use the Herschel SPIRE color-color diagram to study the spectral energy distribution (SED) and the redshift estimation of high-z galaxies. We compiled a sample of 57 galaxies with spectroscopically confirmed redshifts and SPIRE detections in all three bands at $z=2.5-6.4$, and compared their average SPIRE colors with SED templates from local and high-z libraries. We find that local SEDs are inconsistent with high-z observations. The local calibrations of the parameters need to be adjusted to describe the average colors of high-z galaxies. For high-z libraries, the templates with an evolution from z=0 to 3 can well describe the average colors of the observations at high redshift. Using these templates, we defined color cuts to divide the SPIRE color-color diagram into different regions with different mean redshifts. We tested this method and two other color cut methods using a large sample of 783 Herschel-selected galaxies, and find that although these methods can separate the sample into populations with different mean redshifts, the dispersion of redshifts in each population is considerably large. Additional information is needed for better sampling.
This paper presents CCD multicolour photometry for the old open cluster NGC 188. The observations were carried out as a part of the Beijing--Arizona--Taiwan--Connecticut Multicolour Sky Survey from 1995 February to 2008 March, using 15 intermediate-band filters covering 3000--10000 \AA. By fitting the Padova theoretical isochrones to our data, the fundamental parameters of this cluster are derived: an age of $t=7.5\pm 0.5$ Gyr, a distant modulus of $(m-M)_0=11.17\pm0.08$, and a reddening of $E(B-V)=0.036\pm0.010$. The radial surface density profile of NGC 188 is obtained by star count. By fitting the King model, the structural parameters of NGC 188 are derived: a core radius of $R_{c}=3.80'$, a tidal radius of $R_{t}=44.78'$, and a concentration parameter of $C_{0}=\log(R_{t}/R_{c})=1.07$. Fitting the mass function to a power-law function $\phi(m) \propto m^{\alpha}$, the slopes of mass functions for different spatial regions are derived. We find that NGC 188 presents a slope break in the mass function. The break mass is $m_{\rm break}=0.885~M_{\odot}$. In the mass range above $m_{\rm break}$, the slope of the overall region is $\alpha=-0.76$. The slope of the core region is $\alpha=1.09$, and the slopes of the external regions are $\alpha=-0.86$ and $\alpha=-2.15$, respectively. In the mass range below $m_{\rm break}$, these slopes are $\alpha=0.12$, $\alpha=4.91$, $\alpha=1.33$, and $\alpha=-1.09$, respectively. The mass segregation in NGC 188 is reflected in the obvious variation of the slopes in different spatial regions of this cluster.
Results from a six year systematic observational survey on candidate eclipsing binaries with a $\delta$ delta Sct component are presented. More than a hundred systems with component(s) of A-F spectral types were observed in the frame of this survey in order to be checked for possible pulsational behaviour. The 14% (13 cases) of the currently known such systems were discovered during this survey. Using all the available information from the literature, an updated list with all the currently known systems of this type is presented, while possible correlations between their pulsational and binarity properties are discussed.
Shocks around clusters of galaxies accelerate electrons which upscatter the Cosmic Microwave Background photons to higher-energies. We use an analytical model to calculate this inverse Compton (IC) emission, taking into account the effects of additional energy losses via synchrotron and Coulomb scattering. We find that the surface brightness of the optical IC emission increases with redshift and halo mass. The IC emission surface brightness, 32--34~mag~arcsec$^{-2}$, for massive clusters is potentially detectable by the newly developed Dragonfly Telephoto Array.
We discuss the inflationary paradigm, how it can be tested, and how various models of inflation fare in the light of data from Planck and BICEP2. We introduce inflation and reheating, and discuss temperature and polarisation anisotropies in the cosmic microwave background radiation due to quantum fluctuations during inflation. Fitting observations of the anisotropies with theoretical realisations obtained by varying various parameters of the curvature power spectrum and cosmological parameters enables one to obtain the allowed ranges of these parameters. We discuss how to relate these parameters to inflation models which allows one to rule in or out specific models of inflation.
Horizontal Branch stars belong to an advanced stage in the evolution of the oldest stellar galactic population, occurring either as field halo stars or grouped in globular clusters. The discovery of multiple populations in these clusters, that were previously believed to have single populations gave rise to the currently accepted theory that the hottest horizontal branch members (the blue hook stars, which had late helium-core flash ignition, followed by deep mixing) are the progeny of a helium-rich "second generation" of stars. It is not known why such a supposedly rare event (a late flash followed by mixing) is so common that the blue hook of {\omega} Cen contains \sim 30% of horizontal branch stars 10 , or why the blue hook luminosity range in this massive cluster cannot be reproduced by models. Here we report that the presence of helium core masses up to \sim 0.04 solar masses larger than the core mass resulting from evolution is required to solve the luminosity range problem. We model this by taking into account the dispersion in rotation rates achieved by the progenitors, whose premain sequence accretion disc suffered an early disruption in the dense environment of the cluster's central regions where second-generation stars form. Rotation may also account for frequent late-flash-mixing events in massive globular clusters.
We present studies for optimizing the next generation of ground-based imaging atmospheric Cherenkov telescopes (IACTs). Results focus on mid-sized telescopes (MSTs) for CTA, detecting very high energy gamma rays in the energy range from a few hundred GeV to a few tens of TeV. We describe a novel, flexible detector Monte Carlo package, FAST (FAst Simulation for imaging air cherenkov Telescopes), that we use to simulate different array and telescope designs. The simulation is somewhat simplified to allow for efficient exploration over a large telescope design parameter space. We investigate a wide range of telescope performance parameters including optical resolution, camera pixel size, and light collection area. In order to ensure a comparison of the arrays at their maximum sensitivity, we analyze the simulations with the most sensitive techniques used in the field, such as maximum likelihood template reconstruction and boosted decision trees for background rejection. Choosing telescope design parameters representative of the proposed Davies-Cotton (DC) and Schwarzchild-Couder (SC) MST designs, we compare the performance of the arrays by examining the gamma-ray angular resolution and differential point-source sensitivity. We further investigate the array performance under a wide range of conditions, determining the impact of the number of telescopes, telescope separation, night sky background, and geomagnetic field. We find a 30-40% improvement in the gamma-ray angular resolution at all energies when comparing arrays with an equal number of SC and DC telescopes, significantly enhancing point-source sensitivity in the MST energy range. We attribute the increase in point-source sensitivity to the improved optical point-spread function and smaller pixel size of the SC telescope design.
Our tightest upper limit on the sum of neutrino mass eigenvalues $M_\nu$ comes from cosmological observations that will improve substantially in the near future, enabling a detection. The combination of the Baryon Acoustic Oscillation feature measured from the Dark Energy Spectroscopic Instrument and a Stage-IV Cosmic Microwave Background experiment has been forecasted to achieve $\sigma(M_\nu) < 1/3$ of the lower limit on $M_\nu$ from atmospheric and solar neutrino oscillations \citep{2013arXiv1309.5383A,2012PhRvD..86a3012F}. Here we examine in detail the physical effects of neutrino mass on cosmological observables that make these constraints possible. We also consider how these constraints would be improved to ensure at least a $5\sigma$ detection.
PSR J1357$-$6429 is a young and energetic radio pulsar detected in X-rays and $\gamma$-rays. It powers a compact pulsar wind nebula with a jet visible in X-rays and a large scale plerion detected in X-ray and TeV ranges. Previous multiwavelength studies suggested that the pulsar has a significant proper motion of about 180 mas yr$^{-1}$ implying an extremely high transverse velocity of about 2000 km s$^{-1}$. In order to verify that, we performed radio-interferometric observations of PSR J1357$-$6429 with the the Australia Telescope Compact Array (ATCA) in the 2.1 GHz band. We detected the pulsar with a mean flux density of $212\pm5$ $\mu$Jy and obtained the most accurate pulsar position, RA = 13:57:02.525(14) and Dec = $-$64:29:29.89(15). Using the new and archival ATCA data, we did not find any proper motion and estimated its 90 per cent upper limit $\mu < 106$ mas yr$^{-1}$. The pulsar shows a highly polarised single pulse, as it was earlier observed at 1.4 GHz. Spectral analysis revealed a shallow spectral index $\alpha_{\nu}$ = $0.5 \pm 0.1$. Based on our new radio position of the pulsar, we disclaim its optical counterpart candidate reported before.
We present a direct image of the innermost companion to the red giant delta Andromedae using the Stellar Double Coronagraph at the Palomar Observatory. We use a Markov chain Monte Carlo based algorithm to simultaneously reduce the data and perform astrometry and photometry of the companion. We determine that the companion is most likely a main-sequence K-type star and is certainly not the previously hypothesized white dwarf.
A fundamental problem in plasma physics, space science, and astrophysics is the transport of energetic particles interacting with stochastic magnetic fields. In particular the motion of particles across a large scale magnetic field is difficult to describe analytically. However, progress has been achieved in the recent years due to the development of the unified non-linear transport theory which can be used to describe magnetic field line diffusion as well as perpendicular diffusion of energetic particles. The latter theory agrees very well with different independently performed test-particle simulations. However, the theory is still based on different approximations and assumptions. In the current article we extend the theory by taking into account the finite gyroradius of the particle motion and calculate corrections in different asymptotic limits. We consider different turbulence models as examples such as the slab model, noisy slab turbulence, and the two-dimensional model. Whereas there are no finite gyroradius corrections for slab turbulence, the perpendicular diffusion coefficient is reduced in the other two cases. The matter investigated in this article is also related to the parameter a^2 occurring in non-linear diffusion theories.
The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries which mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5-20$ eV and $n_{e}\approx10^{11}-5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.
For light curve generation, a pre-planned photometry survey is needed
nowadays, where all of the exposure coordinates have to be given and don't
change during the survey. This thesis shows it is not required and we can
data-mine these light curves from astronomical data that was never meant for
this purpose. With this approach, we can recycle all of the photometric surveys
in the world and generate light curves of observed objects for them.
This thesis is addressing mostly the catalog generation process, which is
needed for creating the light curves. In practice, it focuses on one of the
most important problems in astroinformatics which is clustering data volumes on
Big Data scale where most of the traditional techniques stagger. We consider a
wide variety of possible solutions from the view of performance, scalability,
distributability, etc. We defined criteria for time and memory complexity which
we evaluated for all of the tested solutions. Furthermore, we created quality
standards which we also take into account when evaluating the results.
We are using relational databases as a starting point of our implementation
and compare them with the newest technologies potentially usable for solving
our problem. These are noSQL Array databases or transferring the heavy
computations of clustering towards supercomputers by using parallelism.
We find an exact formula for the thermally averaged cross section times the relative velocity $\langle \sigma v_{\text{rel}} \rangle$ with relativistic Maxwell-Boltzmann statistics. The formula is valid in the effective field theory approach when the masses of the annihilation products can be neglected compared with the dark matter mass and cut-off scale. The expansion at $x=m/T\gg 1$ directly gives the nonrelativistic limit of $\langle \sigma v_{\text{rel}}\rangle$ which is usually used to compute the relic abundance for heavy particles that decouple when they are nonrelativistic. We compare this expansion with the one obtained by expanding the total cross section $\sigma(s)$ in powers of the nonrelativistic relative velocity $v_r$. We show the correct invariant procedure that gives the nonrelativistic average $\langle \sigma_{nr} v_r \rangle_{nr}$ coinciding with the large $x$ expansion of $\langle \sigma v_{\text{rel}}\rangle$ in the comoving frame. We explicitly formulate flux, cross section, thermal average, collision integral of the Boltzmann equation in an invariant way using the true relativistic relative $v_\text{rel}$, showing the uselessness of the M\o{}ller velocity and further elucidating the conceptual and numerical inconsistencies related with its use.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
Small- and intermediate-scale galaxy clustering can be used to establish the galaxy-halo connection to study galaxy formation and evolution and to tighten constraints on cosmological parameters. With the increasing precision of galaxy clustering measurements from ongoing and forthcoming large galaxy surveys, accurate models are required to interpret the data and extract relevant information. We introduce a method based on high-resolution N-body simulations to accurately and efficiently model the galaxy two-point correlation functions (2PCFs) in projected and redshift spaces. The basic idea is to tabulate all information of haloes in the simulations necessary for computing the galaxy 2PCFs within the framework of halo occupation distribution or conditional luminosity function. It is equivalent to populating galaxies to dark matter haloes and using the mock 2PCF measurements as the model predictions. Besides the accurate 2PCF calculations, the method is also fast and therefore enables an efficient exploration of the parameter space. As an example of the method, we decompose the redshift-space galaxy 2PCF into different components based on the type of galaxy pairs and show the redshift-space distortion effect in each component. The generalizations and limitations of the method are discussed.
COSMOGRAIL is a long-term photometric monitoring of gravitationally lensed QSOs aimed at implementing Refsdal's time-delay method to measure cosmological parameters, in particular H0. Given long and well sampled light curves of strongly lensed QSOs, time-delay measurements require numerical techniques whose quality must be assessed. To this end, and also in view of future monitoring programs or surveys such as the LSST, a blind signal processing competition named Time Delay Challenge 1 (TDC1) was held in 2014. The aim of the present paper, which is based on the simulated light curves from the TDC1, is double. First, we test the performance of the time-delay measurement techniques currently used in COSMOGRAIL. Second, we analyse the quantity and quality of the harvest of time delays obtained from the TDC1 simulations. To achieve these goals, we first discover time delays through a careful inspection of the light curves via a dedicated visual interface. Our measurement algorithms can then be applied to the data in an automated way. We show that our techniques have no significant biases, and yield adequate uncertainty estimates resulting in reduced chi2 values between 0.5 and 1.0. We provide estimates for the number and precision of time-delay measurements that can be expected from future time-delay monitoring campaigns as a function of the photometric signal-to-noise ratio and of the true time delay. We make our blind measurements on the TDC1 data publicly available
We investigate the ISM properties of 13 star-forming galaxies within the z~2 COSMOS cluster. We show that the cluster members have [NII]/Ha and [OIII]/Hb emission-line ratios similar to z~2 field galaxies, yet systematically different emission-line ratios (by ~0.17 dex) from the majority of local star-forming galaxies. We find no statistically significant difference in the [NII]/Ha and [OIII]/Hb line ratios or ISM pressures among the z~2 cluster galaxies and field galaxies at the same redshift. We show that our cluster galaxies have significantly larger ionization parameters (by up to an order of magnitude) than local star-forming galaxies. We hypothesize that these high ionization parameters may be associated with large specific star formation rates (i.e. a large star formation rate per unit stellar mass). If this hypothesis is correct, then this relationship would have important implications for the geometry and/or the mass of stars contained within individual star clusters as a function of redshift.
Large-scale structure surveys in the coming years will measure the redshift-space power spectrum to unprecedented accuracy, allowing for powerful new tests of the LambdaCDM picture as well as measurements of particle physics parameters such as the neutrino masses. We extend the Time-RG perturbative framework to redshift space, computing the power spectrum P_s(k,mu) in massive neutrino cosmologies with time-dependent dark energy equations of state w(z). Time-RG is uniquely capable of incorporating scale-dependent growth into the P_s(k,mu) computation, which is important for massive neutrinos as well as modified gravity models. Although changes to w(z) and the neutrino mass fraction both affect the late-time scale-dependence of the non-linear power spectrum, we find that the two effects depend differently on the line-of-sight angle mu. Finally, we use the HACC N-body code to quantify errors in the perturbative calculations. For a LambdaCDM model at redshift z=1, our procedure predicts the monopole~(quadrupole) to 1% accuracy up to a wave number 0.19h/Mpc (0.28h/Mpc), compared to 0.08h/Mpc (0.07h/Mpc) for the Kaiser approximation and 0.19h/Mpc (0.16h/Mpc) for the current state-of-the-art perturbation scheme. Our calculation agrees with the simulated redshift-space power spectrum even for neutrino masses $\sum$ m_nu ~ 1eV, several times the current bound, as well as rapidly-evolving dark energy equations of state, |dw/dz| ~ 1. Along with this article, we make our redshift-space Time-RG implementation publicly available as the code redTime.
Strong gravitational lensing provides some of the deepest views of the Universe, enabling studies of high-redshift galaxies only possible with next-generation facilities without the lensing phenomenon. To date, 21 cm radio emission from neutral hydrogen has only been detected directly out to z~0.2, limited by the sensitivity and instantaneous bandwidth of current radio telescopes. We discuss how current and future radio interferometers such as the Square Kilometre Array (SKA) will detect lensed HI emission in individual galaxies at high redshift. Our calculations rely on a semi-analytic galaxy simulation with realistic HI disks (by size, density profile and rotation), in a cosmological context, combined with general relativistic ray tracing. Wide-field, blind HI surveys with the SKA are predicted to be efficient at discovering lensed HI systems, increasingly so at z > 2. This will be enabled by the combination of the magnification boosts, the steepness of the HI luminosity function at the high-mass end, and the fact that the HI spectral line is relatively isolated in frequency. These surveys will simultaneously provide a new technique for foreground lens selection and yield the highest redshift HI emission detections. More near term (and existing) cm-wave facilities will push the high redshift HI envelope through targeted surveys of known lenses.
Scaling relations between black hole (BH) masses and their host galaxy properties have been studied extensively over the last two decades, and point towards co-evolution of central massive BHs and their hosts. However, these relations remain poorly constrained for BH masses below $\sim10^{6}$ M_sun. Here we present optical and X-ray observations of the dwarf galaxy RGG 118 taken with the Magellan Echellette Spectrograph on the 6.5m Clay Telescope and Chandra X-ray Observatory. Based on Sloan Digital Sky Survey spectroscopy, RGG 118 was identified as possessing narrow emission line ratios indicative of photoionization partly due to an active galactic nucleus. Our higher resolution spectroscopy clearly reveals broad H$\alpha$ emission in the spectrum of RGG 118. Using virial BH mass estimate techniques, we calculate a BH mass of $\sim50,000$ \msun. We detect a nuclear X-ray point source in RGG 118, suggesting a total accretion powered luminosity of $L=4\times10^{40}~{\rm erg~s^{-1}}$, and an Eddington fraction of $\sim1$ per cent. The BH in RGG 118 is the smallest ever reported in a galaxy nucleus and we find that it lies on the extrapolation of the $M_{\rm BH}-\sigma_{\ast}$ relation to the lowest masses yet.
Young tidal dwarf galaxies (TDGs) are observed in the tidal debris of gas-rich interacting galaxies. In contrast to what is generally assumed to be the case for isolated dwarf galaxies, TDGs are not embedded in their own dark matter (DM) sub-halo. Hence, they are more sensitive to stellar feedback and could be disrupted on a short time-scale. Detailed numerical and observational studies demonstrate that isolated DM-dominated dwarf galaxies can have lifetimes of more than 10 Gyr. For TDGs that evolve in a tidal field with compressing accelerations equal to the gravitational acceleration within a DM sub-halo typical of an isolated dwarf galaxy, a similar survival time is expected. The tidal acceleration profile depends on the virial mass of the host galaxy and the distance between the TDG and its host. We analytically compare the tidal compression to the gravitational acceleration due to either cuspy or cored DM sub-halos of various virial masses. For example, the tidal field at a distance of 100 kpc to a host halo of 10^13 Msol can be as stabilizing as a 10^9 Msol DM sub-halo. By linking the tidal field to the equivalent gravitational field of a DM sub-halo, we can use existing models of isolated dwarfs to estimate the survivability of TDGs. We show that part of the unexpectedly high dynamical masses inferred from observations of some TDGs can be explained by tidal compression and hence TDGs require to contain less unobservable matter to understand their rotation curves.
Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their cores. We present a study of 191 unobscured quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years which was subsequently quenched. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.
We use the Very Large Array (VLA) in the A configuration with the Pie Town (PT) Very Long Baseline Array (VLBA) antenna to spatially resolve the extended atmosphere of Betelgeuse over multiple epochs at 0.7, 1.3, 2.0, 3.5, and 6.1 cm. The extended atmosphere deviates from circular symmetry at all wavelengths while at some epochs we find possible evidence for small pockets of gas significantly cooler than the mean global temperature. We find no evidence for the recently reported e-MERLIN radio hotspots in any of our multi-epoch VLA/PT data, despite having sufficient spatial resolution and sensitivity at short wavelengths, and conclude that these radio hotspots are most likely interferometric artefacts. The mean gas temperature of the extended atmosphere has a typical value of 3000 K at 2 $R_{\star}$ and decreases to 1800 K at 6 $R_{\star}$, in broad agreement with the findings of the single epoch study from Lim et al. (1998). The overall temperature profile of the extended atmosphere between $2 R_{\star} \lesssim r \lesssim 6 R_{\star}$ can be described by a power law of the form $T_{\mathrm{gas}}(r) \propto r^{-0.6}$, with temporal variability of a few 100 K evident at some epochs. Finally, we present over 12 years of V band photometry, part of which overlaps our multi-epoch radio data. We find a correlation between the fractional flux density variability at V band with most radio wavelengths. This correlation is likely due to shock waves induced by stellar pulsations, which heat the inner atmosphere and ionize the more extended atmosphere through radiative means. Stellar pulsations may play an important role in exciting Betelgeuse's extended atmosphere.
For X-ray imaging spectroscopy, high spatial resolution over a large field of view is often as important as high energy resolution, but current X-ray detectors do not provide both in the same device. Thermal Kinetic Inductance Detectors (TKIDs) are being developed as they offer a feasible way to combine the energy resolution of transition edge sensors with pixel counts approaching CCDs and thus promise significant improvements for many X-ray spectroscopy applications. TKIDs are a variation of Microwave Kinetic Inductance Detectors (MKIDs) and share their multiplexibility: working MKID arrays with 2024 pixels have recently been demonstrated and much bigger arrays are under development. In this work, we present our first working TKID prototypes which are able to achieve an energy resolution of 75 eV at 5.9 keV, even though their general design still has to be optimized. We further describe TKID fabrication, characterization, multiplexing and working principle and demonstrate the necessity of a data fitting algorithm in order to extract photon energies. With further design optimizations we expect to be able to improve our TKID energy resolution to less than 10 eV at 5.9 keV.
Exoplanets orbiting close to their parent stars could lose some fraction of their atmospheres because of the extreme irradiation. Atmospheric mass loss primarily affects low-mass exoplanets, leading to suggest that hot rocky planets might have begun as Neptune-like, but subsequently lost all of their atmospheres; however, no confident measurements have hitherto been available. The signature of this loss could be observed in the ultraviolet spectrum, when the planet and its escaping atmosphere transit the star, giving rise to deeper and longer transit signatures than in the optical spectrum. Here we report that in the ultraviolet the Neptune-mass exoplanet GJ 436b (also known as Gliese 436b) has transit depths of 56.3 +/- 3.5% (1 sigma), far beyond the 0.69% optical transit depth. The ultraviolet transits repeatedly start ~2 h before, and end >3 h after the ~1 h optical transit, which is substantially different from one previous claim (based on an inaccurate ephemeris). We infer from this that the planet is surrounded and trailed by a large exospheric cloud composed mainly of hydrogen atoms. We estimate a mass-loss rate in the range of ~10^8-10^9 g/s, which today is far too small to deplete the atmosphere of a Neptune-like planet in the lifetime of the parent star, but would have been much greater in the past.
We examine a broad range of mass flow diagnostics in a large sample of Herbig Ae/Be stars (HAEBES) using high resolution optical spectra. The H-beta and He I 5876 angstrom lines show the highest incidence of P-Cygni (30%) and inverse P-Cygni (14%) morphologies, respectively. The Fe II 4924 angstrom line also shows a large incidence of P-Cygni profiles (11%). We find support for many of the conclusions reached in a study based on the analysis of the He I 10830 angstrom line in a large sample of HAEBES. Namely, HAEBES exhibit smaller fractions of both blue-shifted absorption (i.e. mass outflow) and red-shifted absorption (i.e. mass infall or accretion) than their lower mass cousins, the classical T Tauri stars (CTTSs). In particular, the optical data supports the conclusion that HAEBES displaying red-shifted absorption, in general, show maximum red-shifted absorption velocities that are smaller fractions of their stellar escape velocities than is found for CTTSs. This suggests that HAEBE accretion flows are originating deeper in the gravitational potentials of their stars than in CTTS systems. In addition, we find a lack of inner disk wind signatures in the blue-shifted absorption objects; only stellar wind signatures are clearly observed. These findings, along with the lack of detected magnetic fields around HAEBES, support the idea that large magnetospheres are not prevalent around HAEBES and that accretion flows are instead mediated by significantly smaller magnetospheres with relatively smaller truncation radii (e.g. 1-2 stellar radii). Red-shifted absorption is much more common around Herbig Ae stars than Be stars, suggesting that Herbig Be stars may accrete via a boundary layer rather than along magnetic field lines.
The positions and velocities of galaxies in the Local Group (LG) measure the gravitational field within it. This is mostly due to the Milky Way (MW) and Andromeda (M31). We constrain their masses using a sample of 32 galaxies with measured distances and radial velocities (RVs). To do this, we follow the trajectories of several thousand simulated particles on a pure Hubble flow from redshift 9. For each observed galaxy, we obtain a trajectory which today is at the same position. Its final velocity is the model prediction for the velocity of that galaxy. We carefully consider the impact of tides raised by objects outside the LG. We directly include Centaurus A and try to account for IC 342 and M81. With our analysis, the total LG mass is $4.33^{+0.37}_{-0.32} \times {10}^{12} M_\odot$, with $0.20^{+0.05}_{-0}$ of this being in the MW. However, no plausible set of initial conditions yields a good match to the RVs of our sample of LG galaxies. We introduce a parameter $\sigma_{extra}$ to quantify the typical disagreement between observed RVs and those predicted by the best-fitting model. We find that $\sigma_{extra} \approx 45^{+7}_{-5}$ km/s. This seems too high to explain as a result of interactions between LG dwarf galaxies. We suggest that the observations may be explained by a past close flyby of the MW and M31, which arises in some modified gravity theories due to a shorter orbital period. Gravitational slingshot encounters of material in the LG with either of these massive fast-moving galaxies could plausibly explain why some non-satellite LG galaxies are racing away from the LG even faster than a pure Hubble flow (e.g. DDO 99, 125 and 190). A modification to gravity might also explain why some galaxies have RVs substantially below our model predictions.
Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. As an urgent response to these challenges, in this paper we propose an adaptive storage scheme for simulation data, inspired by the block time step integration scheme found in a number of direct $N$-body integrators available nowadays. The proposed scheme, namely the block time step storage scheme, works by minimizing the data redundancy with assignments of data with individual output frequencies as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct $N$-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used.
Due to the chaotic nature of the Solar System, the question of its long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner planets, leading to close encounters and/or collisions can be initiated through a large increase in Mercury's eccentricity, with a currently assumed likelihood of ~1%. However, little is known at present about the robustness of this number. Here I report ensemble integrations of the full equations of motion of the eight planets and Pluto over 5 Gyr, including contributions from general relativity. The results show that different numerical algorithms lead to statistically different results for the evolution of Mercury's eccentricity (eM). For instance, starting at present initial conditions (eM ~= 0.21), Mercury's maximum eccentricity achieved over 5 Gyr is on average significantly higher in symplectic ensemble integrations using heliocentricthan Jacobi coordinates and stricter error control. In contrast, starting at a possible future configuration (eM ~= 0.53), Mercury's maximum eccentricity achieved over the subsequent 500 Myr is on average significantly lower using heliocentric than Jacobi coordinates. For example, the probability for eM to increase beyond 0.53 over 500 Myr is >90% (Jacobi) vs. only 40-55% (heliocentric). This poses a dilemma as the physical evolution of the real system - and its probabilistic behavior - cannot depend on the coordinate system or numerical algorithm chosen to describe it. Some tests of the numerical algorithms suggest that symplectic integrators using heliocentric coordinates underestimate the odds for destabilization of Mercury's orbit at high initial eM.
By examining the locations of central black holes in two elliptical galaxies, M32 and M87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than eight percent, $|\eta_N|<0.08$ from M32. The constraints from M87 are substantially weaker but could improve dramatically with better astrometry.
The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15,000 objects during the cryogenic mission lifetime. Observations provided low-resolution (R~60-127) spectra over ~5-38um and high-resolution (R~600) spectra over ~10-37um. The Cornell Atlas of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. Low-resolution spectra have been available in CASSIS since 2011, and we present here the addition of the high-resolution spectra. The high-resolution observations represent approximately one third of all staring observations performed with the IRS instrument. While low-resolution observations are adapted to faint objects and/or broad spectral features (e.g., dust continuum, molecular bands), high-resolution observations allow more accurate measurements of narrow features (e.g., ionic emission lines) as well as a better sampling of the spectral profile of various features. Given the narrow aperture of the two high-resolution modules, cosmic ray hits and spurious features usually plague the spectra. Our pipeline is designed to minimize these effects through various improvements. A super sampled point-spread function was created in order to enable the optimal extraction in addition to the full aperture extraction. The pipeline selects the best extraction method based on the spatial extent of the object. For unresolved sources, the optimal extraction provides a significant improvement in signal-to-noise ratio over a full aperture extraction. We have developed several techniques for optimal extraction, including a differential method that eliminates low-level rogue pixels (even when no dedicated background observation was performed). The updated CASSIS repository now includes all the spectra ever taken by the IRS, with the exception of mapping observations.
We present an update of the CLUMPY code for the calculation of the astrophysical J-factors (from dark matter annihilation/decay) for any Galactic or extragalactic dark matter halo including substructures: the concentration-mass relationship may now be drawn from a distribution, boost factors can include several levels of substructures, and triaxiality is a new option for dark matter haloes. This new version takes advantage of the cfitsio and HEALPix libraries to propose FITS output maps using the HEALPix pixelisation scheme. Skymaps for $\gamma$-ray and neutrino signals from generic annihilation/decay spectra are now direct outputs of CLUMPY. Smoothing by a user-defined instrumental Gaussian beam is also possible. In addition to these improvements, the main novelty is the implementation of a Jeans analysis module, to obtain dark matter density profiles from kinematic data in relaxed spherical systems (e.g., dwarf spheroidal galaxies). The code is also interfaced with the GreAT toolkit designed for Markov Chain Monte Carlo analyses, from which probability density functions and credible intervals can be obtained for velocity dispersions, dark matter profiles, and J- factors.
Our previous semi-analytic treatment of f- and p-mode multiple scattering from ensembles of thin flux tubes (Hanson and Cally, Astrophys. J. 781, 125; 791, 129, 2014) is extended by allowing both sausage and kink waves to freely escape at the top of the model using a radiative boundary condition there. As expected, this additional avenue of escape, supplementing downward loss into the deep solar interior, results in substantially greater absorption of incident f- and p-modes. However, less intuitively, it also yields mildly to substantially smaller phase shifts in waves emerging from the ensemble. This may have implications for the interpretation of seismic data for solar plage regions, and in particular their small measured phase shifts.
Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as "Einstein Rings" in high resolution images. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope $\gamma'$ of the total mass density profile $\rho_{tot}\propto r^{-\gamma'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $\gamma'$ be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. However, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.
In this letter we explore the environment of Pluto and Charon in the far infrared with the main aim to identify the signs of any possible dust ring, should it exist in the system. Our study is based on observations performed at 70 um with the PACS instrument onboard the Herschel Space Observatory at 9 epochs between March 14 and 19, 2012. The far-infrared images of the Pluto-Charon system are compared to those of the point spread function (PSF) reference quasar 3C454.3. The deviation between the observed Pluto-Charon and reference PSFs are less then 1 sigma indicating that clear evidence for an extended dust ring around the system was not found. Our method is capable of detecting a hypothetical ring with a total flux of ~3.3 mJy at a distance of ~153 000 km (~8.2 Pluto-Charon distances) from the system barycentre. We place upper limits on the total disk mass and on the column density in a reasonable disk configuration and analyse the hazard during the flyby of NASAs New Horizons in July 2015. This realistic model configuration predicts a column density of 8.7x10^(-10) gcm^(-2) along the path of the probe and an impactor mass of 8.7x10^(-5) g.
Tidal dissipation in planets and stars is one of the key physical mechanisms driving the evolution of star-planet and planet-moon systems. Several signatures of its action are observed in planetary systems thanks to their orbital architecture and the rotational state of their components. Tidal dissipation inside the fluid layers of celestial bodies are intrinsically linked to the dynamics and the physical properties of the latter. This complex dependence must be characterized. We compute the tidal kinetic energy dissipated by viscous friction and thermal diffusion in a rotating local fluid Cartesian section of a star/planet/moon submitted to a periodic tidal forcing. The properties of tidal gravito-inertial waves excited by the perturbation are derived analytically as explicit functions of the tidal frequency and local fluid parameters (i.e. the rotation, the buoyancy frequency characterizing the entropy stratification, viscous and thermal diffusivities) for periodic normal modes. The sensitivity of the resulting possibly highly resonant dissipation frequency-spectra to a control parameter of the system is either important or negligible depending on the position in the regime diagram relevant for planetary and stellar interiors. For corresponding asymptotic behaviors of tidal gravito-inertial waves dissipated by viscous friction and thermal diffusion, scaling laws for the frequencies, number, width, height and contrast with the non-resonant background of resonances are derived to quantify these variations. We characterize the strong impact of the internal physics and dynamics of fluid planetary layers and stars on the dissipation of tidal kinetic energy in their bulk. We point out the key control parameters that really play a role and demonstrate how it is now necessary to develop ab-initio modeling for tidal dissipation in celestial bodies.
[Abbreviated] We present the results of a near-infrared (NIR) spectroscopic follow-up survey of 182 M4-L7 low-mass stars and brown dwarfs (BDs) from the BANYAN All-Sky Survey (BASS) for candidate members of nearby, young moving groups (YMGs). We confirm signs of low-gravity for 42 new BD discoveries with estimated masses between 8-75 $M_{Jup}$ and identify previously unrecognized signs of low gravity for 24 known BDs. This allows us to refine the fraction of low-gravity dwarfs in the high-probability BASS sample to $\sim$82%. We use this unique sample of 66 young BDs, supplemented with 22 young BDs from the literature, to construct new empirical NIR absolute magnitude and color sequences for low-gravity BDs. We obtain a spectroscopic confirmation of low-gravity for 2MASS J14252798-3650229, which is a new $\sim$27 $M_{Jup}$, L4 $\gamma$ bona fide member of AB Doradus. We identify a total of 19 new low-gravity candidate members of YMGs with estimated masses below 13 $M_{Jup}$, seven of which have kinematically estimated distances within 40 pc. These objects will be valuable benchmarks for a detailed atmospheric characterization of planetary-mass objects with the next generation of instruments. We find 16 strong candidate members of the Tucana-Horologium association with estimated masses between 12.5-14 $M_{Jup}$, a regime where our study was particularly sensitive. This would indicate that for this association there is at least one isolated object in this mass range for every $17.5_{-5.0}^{+6.6}$ main-sequence stellar member, a number significantly higher than expected based on standard log-normal initial mass function, however in the absence of radial velocity and parallax measurements for all of them, it is likely that this over-density is caused by a number of young interlopers from other moving groups. We identify 12 new L0-L5 field BDs, seven of which display peculiar properties.
Context: Transverse loop oscillations and loop contractions are commonly
associated with solar flares, but the two types of motion have traditionally
been regarded as separate phenomena.
Aims: We present an observation of coronal loops contracting and oscillating
following onset of a flare. We aim to explain why both behaviours are seen
together and why only some of the loops oscillate.
Methods: A time sequence of SDO/AIA 171 \r{A} images is analysed to identify
positions of coronal loops following the onset of M6.4 flare
SOL2012-03-09T03:53. We focus on five loops in particular, all of which
contract during the flare, with three of them oscillating as well. A simple
model is then developed for contraction and oscillation of a coronal loop.
Results: We propose that coronal loop contractions and oscillations can occur
in a single response to removal of magnetic energy from the corona. Our model
reproduces the various types of loop motion observed and explains why the
highest loops oscillate during their contraction while no oscillation is
detected for the shortest contracting loops. The proposed framework suggests
that loop motions can be used as a diagnostic for the removal of coronal
magnetic energy by flares, while rapid decrease of coronal magnetic energy is a
newly-identified excitation mechanism for transverse loop oscillations.
We report the first high angular resolution imaging (3.4\arcsec $\times$ 3.0\arcsec) of deuterated formaldehyde (HDCO) toward Orion--KL, carried out with the Submillimeter Array (SMA). We find that the spatial distribution of the formaldehyde emission systematically differs from that of methanol: while methanol is found towards the inner part of the region, HDCO is found in colder gas that wraps around the methanol emission on four sides. The HDCO/H$_2$CO ratios are determined to be 0.003--0.009 within the region, up to an order of magnitude higher than the D/H measured for methanol. These findings strengthen the previously suggested hypothesis that there are differences in the chemical pathways leading to HDCO (via deuterated gas phase chemistry) and deuterated methanol (through conversion of formaldehyde into methanol on the surface of icy grain mantles).
The identification of unsubtracted foreground residuals in the cosmic microwave background maps on large scales is of crucial importance for the analysis of polarization signals. These residuals add a non-Gaussian contribution to the data. We propose the Kullback-Leibler (KL) divergence as an effective, non-parametric test on the one-point probability distribution function of the data. With motivation in information theory, the KL divergence takes into account the entire range of the distribution and is highly non-local. We demonstrate its use by analyzing the large scales of the Planck 2013 SMICA temperature fluctuation map and find it consistent with the expected distribution at a level of 6%. Comparing the results to those obtained using the more popular Kolmogorov-Smirnov test, we find the two methods to be in general agreement.
We present an analysis of the high-resolution Chandra observation of the multiple system, HD 155555 (an RS CVn type binary system, HD 155555 AB, and its spatially resolved low-mass companion HD 155555 C). This is an intriguing system which shows properties of both an active pre-main sequence star and a synchronised (main sequence) binary. We obtain the emission measure distribution, temperature structures, plasma densities, and abundances of this system and compare them with the coronal properties of other young/active stars. HD 155555 AB and HD 155555 C produce copious X-ray emission with log Lx of 30.54 and 29.30, respectively, in the 0.3-6.0 keV energy band. The light curves of individual stars show variability on timescales of few minutes to hours. We analyse the dispersed spectra and reconstruct the emission measure distribution using spectral line analysis. The resulting elemental abundances exhibit inverse first ionisation potential effect in both cases. An analysis of He-like triplets yields a range of coronal electron densities ~10^10-10^13 cm-3. Since HD 155555 AB is classified both as an RS CVn and a PMS star, we compare our results with those of other slightly older active main-sequence stars and T Tauri stars, which indicates that the coronal properties of HD 155555 AB closely resemble that of an older RS CVn binary rather than a younger PMS star. Our results also suggests that the properties of HD 155555 C is very similar to those of other active M dwarfs.
The design of spaceborne high-energy (E>100 MeV) gamma-ray detectors depends on two principal factors: (1) the basic physics of detecting and measuring the properties of the gamma rays; and (2) the constraints of operating such a detector in space for an extended period. Improvements in technology have enabled major advances in detector performance, as illustrated by two successful instruments, EGRET on the Compton Gamma Ray Observatory and LAT on the Fermi Gamma-ray Space Telescope.
Einstein Probe is a small mission dedicated to time-domain high-energy astrophysics. Its primary goals are to discover high-energy transients and to monitor variable objects in the $0.5-4~$keV X-rays, at higher sensitivity by one order of magnitude than those of the ones currently in orbit. Its wide-field imaging capability, featuring a large instantaneous field-of-view ($60^\circ \times60^\circ$, $\sim1.1$sr), is achieved by using established technology of micro-pore (MPO) lobster-eye optics, thereby offering unprecedentedly high sensitivity and large Grasp. To complement this powerful monitoring ability, it also carries a narrow-field, sensitive follow-up X-ray telescope based on the same MPO technology to perform follow-up observations of newly-discovered transients. Public transient alerts will be downlinked rapidly, so as to trigger multi-wavelength follow-up observations from the world-wide community. Over three of its 97-minute orbits almost the entire night sky will be sampled, with cadences ranging from 5 to 25 times per day. The scientific objectives of the mission are: to discover otherwise quiescent black holes over all astrophysical mass scales by detecting their rare X-ray transient flares, particularly tidal disruption of stars by massive black holes at galactic centers; to detect and precisely locate the electromagnetic sources of gravitational-wave transients; to carry out systematic surveys of X-ray transients and characterize the variability of X-ray sources. Einstein Probe has been selected as a candidate mission of priority (no further selection needed) in the Space Science Programme of the Chinese Academy of Sciences, aiming for launch around 2020.
This white paper is a summarising report of the Forum on monitoring the transient X-ray Universe in the multi-messenger era organized by the International Space Science Institute in Beijing (ISSI-BJ) on May 6-7, 2014. Time-domain astronomy will enter a golden era towards the end of this decade with the advent of major facilities across the electromagnetic spectrum and in the multi-messenger realms of gravitational wave and neutrino. In the soft X-ray regime, the novel micro-pore lobster-eye optics provides a promising technology to realise, for the first time, focusing X-ray optics for wide-angle monitors to achieve a good combination of sensitivity and wide field of view. In this context, Einstein Probe - a soft X-ray all-sky monitor - has been proposed and selected as a candidate mission of priority in the space science programme of the Chinese Academy of Sciences. This report summarises the most important science developments in this field towards 2020 and beyond and how to achieve them technologically, which were discussed at this brainstorming forum. It also introduces briefly the Einstein Probe mission, including its key science goals and mission definition, as well as some of the key technological issues.
Accurate positional measurements of planets and satellites are used to improve our knowledge of their dynamics and to infer the accuracy of planet and satellite ephemerides. In the framework of the FP7 ESPaCE project, we provide the positions of Mars, Phobos, and Deimos taken with the U.S. Naval Observatory 26-inch refractor during the 1971 opposition of the planet. These plates were measured with the digitizer of the Royal Observatory of Belgium and reduced through an optimal process that includes image, instrumental, and spherical corrections to provide the most accurate data. We compared the observed positions of the planet Mars and its satellites with the theoretical positions from INPOP10 and DE430 planetary ephemerides, and from NOE and MAR097 satellite ephemerides. The rms residuals in RA and Dec. of one position is less than 60 mas, or about 20 km at Mars. This accuracy is comparable to the most recent CCD observations. Moreover, it shows that astrometric data derived from photographic plates can compete with those of old spacecraft (Mariner 9, Viking 1 and 2).
We present solution of the problem of wavelength calibration for Hamilton Echelle spectrograph using hollow cathode lamp, which was operated at Lick Observatory Shane telescope before June 9, 2011. The spectrum of the lamp claimed to be thorium-argon, contains, in addition to the lines of thorium and argon, a number of the unrecognized lines identified by us with tantalum. Using atomic data for measured lines of tantalum and thorium, we estimated the temperature of the gas in the lamp as T=3120+/-60 K. From the atomic line database VALD3 we selected all lines of TaI and TaII which can be seen in the spectrum of the lamp and compiled a list for the use in the processing of spectral observations. We note a limitation of the accuracy of calibration due to the influence of the hyperfine line splitting.
We perform a Bayesian analysis in the basis of a recently developed two-parameter class of hybrid equations of state that allow for high-mass compact star twins. While recently a wide range of radii, from 9 - 15 km, has been inferred for different neutron stars using different techniques, we perform our analysis under the supposition that the radii are towards the large end ($13-15$ km). We use this radius constraint together with the undebated statistically independent constraint for high masses ($\sim 2~M_\odot$) as priors in selecting the most probable hybrid equations of state from a family with two free parameters: the baryon excluded volume in the hadronic phase and the 8-quark vector channel interaction in the quark matter phase.
Recent studies of active galactic nuclei (AGN) found a statistical inverse linear scaling between the X-ray normalized excess variance $\sigma_{\rm rms}^2$ (variability amplitude) and the black hole mass spanning over $M_{\rm BH}=10^6- 10^9\ M_{\odot}$. Being suggested to have a small scatter, this scaling relation may provide a novel method to estimate the black hole mass of AGN. However, a question arises as to whether this relation can be extended to the low-mass regime below $\sim10^6\ M_{\odot}$. If confirmed, it would provide an efficient tool to search for AGN with low-mass black holes using X-ray variability. This paper presents a study of the X-ray excess variances for a sample of AGN with black hole masses in the range of $10^5- 10^6\ M_{\odot}$ observed with {\it XMM-Newton} and {\it ROSAT}, including data both from the archives and from newly preformed observations. It is found that the relation is no longer a simple extrapolation of the linear scaling; instead, the relation starts to flatten at $\sim10^6\ M_{\odot}$ toward lower masses. Our result is consistent with the recent finding of \citet{L15}. Such a flattening of the $M_{\rm BH}-\sigma_{\rm rms}^2$ relation is actually expected from the shape of the power spectrum density of AGN, whose break frequency is inversely scaled with the mass of black holes.
We introduce a designer approach for extended Brans-Dicke gravity and use it to construct a semi-analytic model for the effective equation of state, the effective Newton's constant at the background and at the linear level and the gravitational slip. By doing so, we are able to explore the dependence of these four phenomenological parameters on more fundamental parameters of the theory.
Long gamma-ray bursts (LGRBs) are associated with the death of massive stars. Their host galaxies therefore represent a unique class of objects tracing star formation across the observable Universe. Indeed, recently accumulated evidence shows that GRB hosts do not differ substantially from general population of galaxies at high (z > 2) redshifts. However, it has been long recognised that the properties of z < 1.5 hosts, compared to general star-forming population, are unusual. To better understand the reasons for the supposed difference in LGRB hosts properties at z < 1.5, we obtained VLT/X- Shooter spectra of six hosts lying in the redshift range of 0.8 < z < 1.3. Some of these hosts have been observed before, yet we still lack well constrained information on their characteristics such as metallicity, dust extinction and star formation rate. We search for emission lines in the VLT/X-Shooter spectra of the hosts and measure their fluxes. We perform a detailed analysis, estimating host average extinction, star-formation rates, metallicities and electron densities where possible. Measured quantities of our hosts are compared to a larger sample of previously observed GRB hosts at z < 2. Star-formation rates and metallicities are measured for all the hosts analyzed in this paper and metallicities are well determined for 4 hosts. The mass-metallicity relation, the fundamental metallicity relation and SFRs derived from our hosts occupy similar parameter space as other host galaxies investigated so-far at the same redshift. We therefore conclude that GRB hosts in our sample support the found discrepancy between the properties of low-redshift GRB hosts and the general population of star- forming galaxies.
We present a full description of the N-probability density function of the galaxy number density fluctuations. This N-pdf is given in terms, on the one hand, of the cold dark matter correlations and, on the other hand, of the galaxy bias parameter. The method relies on the assumption commonly adopted that the dark matter density fluctuations follow a local non-linear transformation of the initial energy density perturbations. The N-pdf of the galaxy number density fluctuations allows for an optimal estimation of the bias parameter (e.g., via maximum-likelihood estimation, or Bayesian inference if there exists any \emph{a priori} information on the bias parameter). It also provides the proper framework to perform model selection between two competitive hypotheses (e.g., galaxy biasing versus a one-to-one relation between the galaxy number density and the dark matter perturbations). The bias estimation capabilities of the N-pdf are proved by SDSS-like simulations, showing that our estimator is unbiased. We apply our formalism to the 7th release of the SDSS main sample (for a volume-limited subset with absolute magnitudes $M_r \leq -20$). We obtain a maximum-likelihood bias estimate $\hat{b} = 1.500 \pm 0.036$, for galaxy number density fluctuations in cells of the size of $30h^{-1}$Mpc. Different model selection criteria show that galaxy biasing is clearly favoured.
Two classes of GRBs have been identified thus far and are prescribed to different physical scenarios $-$ NS-NS or NS-BH mergers, and collapse of massive stars, for short and long GRBs, respectively. A third, intermediate in durations class, was suggested to be present in previous catalogs, such as BATSE and Swift, based on statistical tests regarding a mixture of two or three normal distributions. However, this might possibly not be an adequate model. This paper investigates whether the distribution of $\log T_{90}$ from Fermi shows evidence for a third, intermediate, class of GRBs. Mixtures of standard Gaussians, skew-normal, sinh-arcsinh and alpha-skew-normal distributions are fitted using a maximum likelihood method. The preferred model is chosen based on the Akaike information criterion. It is found that mixtures of two skew-normal or two sinh-arcsinh distributions are more likely to describe the observed duration distribution than a mixture of three standard Gaussians. Based on statistical reasoning, existence of a third (intermediate) class of GRBs in Fermi data is rejected.
We cross-match galaxy cluster candidates selected via their Sunyaev-Zel'dovich effect (SZE) signatures in 129.1 deg$^2$ of the South Pole Telescope 2500d SPT-SZ survey with optically identified clusters selected from the Dark Energy Survey (DES) science verification data. We identify 25 clusters between $0.1\lesssim z\lesssim 0.8$ in the union of the SPT-SZ and redMaPPer (RM) samples. RM is an optical cluster finding algorithm that also returns a richness estimate for each cluster. We model the richness $\lambda$-mass relation with the following function $\langle\ln\lambda|M_{500}\rangle\propto B_\lambda\ln M_{500}+C_\lambda\ln E(z)$ and use SPT-SZ cluster masses and RM richnesses $\lambda$ to constrain the parameters. We find $B_\lambda= 1.14^{+0.21}_{-0.18}$ and $C_\lambda=0.73^{+0.77}_{-0.75}$. The associated scatter in mass at fixed richness is $\sigma_{\ln M|\lambda} = 0.18^{+0.08}_{-0.05}$ at a characteristic richness $\lambda=70$. We demonstrate that our model provides an adequate description of the matched sample, showing that the fraction of SPT-SZ selected clusters with RM counterparts is consistent with expectations and that the fraction of RM selected clusters with SPT-SZ counterparts is in mild tension with expectation. We model the optical-SZE cluster positional offset distribution with the sum of two Gaussians, showing that it is consistent with a dominant, centrally peaked population and a sub-dominant population characterized by larger offsets. We also cross-match the RM catalog with SPT-SZ candidates below the official catalog threshold significance $\xi=4.5$, using the RM catalog to provide optical confirmation and redshifts for additional low-$\xi$ SPT-SZ candidates.In this way, we identify 15 additional clusters with $\xi\in [4,4.5]$ over the redshift regime explored by RM in the overlapping region between DES science verification data and the SPT-SZ survey.
Radiation damage to space-based Charge-Coupled Device (CCD) detectors creates defects which result in an increasing Charge Transfer Inefficiency (CTI) that causes spurious image trailing. Most of the trailing can be corrected during post-processing, by modelling the charge trapping and moving electrons back to where they belong. However, such correction is not perfect -- and damage is continuing to accumulate in orbit. To aid future development, we quantify the limitations of current approaches, and determine where imperfect knowledge of model parameters most degrade measurements of photometry and morphology. As a concrete application, we simulate $1.5\times10^{9}$ "worst case" galaxy and $1.5\times10^{8}$ star images to test the performance of the Euclid visual instrument detectors. There are two separable challenges: If the model used to correct CTI is perfectly the same as that used to add CTI, $99.68$ % of spurious ellipticity is corrected in our setup. This is because readout noise is not subject to CTI, but gets over-corrected during correction. Second, if we assume the first issue to be solved, knowledge of the charge trap density within $\Delta\rho/\rho\!=\!(0.0272\pm0.0005)$ %, and the characteristic release time of the dominant species to be known within $\Delta\tau/\tau\!=\!(0.0400\pm0.0004)$ % will be required. This work presents the next level of definition of in-orbit CTI calibration procedures for Euclid.
Two classes of GRBs have been identified thus far without doubt and are prescribed to different physical scenarios $-$ NS-NS or NS-BH mergers, and collapse of massive stars, for short and long GRBs, respectively. The existence of two distinct populations was inferred through a bimodal distribution of the observed durations $T_{90}$, and the commonly applied $2\,{\rm s}$ limit between short and long GRBs was obtained by fitting a parabola between the two peaks in binned data from BATSE 1B. Herein, by means of a maximum likelihood (ML) method a mixture of two Gaussians is fitted to the datasets from BATSE, Swift, BeppoSAX, RHESSI, and Fermi in search for a local minimum that might serve as a new, more proper, limit for the two GRB classes. It is found that Swift and BeppoSAX distributions are unimodal, hence no local minimum is present, Fermi is consistent with the conventional limit, whereas RHESSI and BATSE give the limits significantly shorter ($1.28\pm 0.14\,{\rm s}$) and longer ($3.38\pm 0.27\,{\rm s}$) than $2\,{\rm s}$, respectively. These new values change the fractions of short and long GRBs in the samples examined, and imply that the observed $T_{90}$ durations are detector dependent, hence no universal limiting value may be applied to all satellites due to their different instrumental specifications. Because of this, and due to the strong overlap of the two-Gaussian components, the straightforward association of short GRBs to mergers and long one to collapsars is ambiguous.
Here we aim to understand the origin of 47 Tuc's La-rich star Lee 4710. We report abundances for O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Co, Ni, Zn, Y, Zr, Ba, La, Ce, Pr, Nd, and Eu, and present a detailed abundance analysis of two 47 Tuc stars with similar stellar parameters but different slow neutron-capture (s-)process enrichment. Star Lee 4710 has the highest known La abundance ratio in this cluster ([La/Fe] = 1.14), and star Lee 4626 is known to have normal s-process abundances (e.g., [Ba/Eu]$<0$). The nucleosynthetic pattern of elements with Z$\gtrsim$56 for star Lee 4710 agrees with the predicted yields of a $1.3M_{\odot}$ asymptotic giant branch (AGB) star. Therefore, Lee 4710 may have been enriched by mass transfer from a more massive AGB companion, which is compatible with its location far away from the center of this relatively metal-rich ([Fe/H]$\sim-0.7$) globular cluster. A further analysis comparing the abundance pattern of Lee 4710 with data available in the literature reveals that nine out of the $\sim200$ 47 Tuc stars previously studied show strong s-process enhancements that point towards later enrichment by more massive AGB stars.
Self-induced flavor conversions of supernova (SN) neutrinos have been characterized in the spherically symmetric "bulb model", reducing the neutrino evolution to a one dimensional problem along a radial direction. We lift this assumption, presenting a two-dimensional model where neutrinos are launched from a spherical neutrino-sphere with many zenithal angles and two azimuthal angles. We also assume that self-induced conversions are not suppressed by large matter effects. In this situation we find that self-interacting neutrinos spontaneously break axial and spherical symmetries. As a result the flavor content and the lepton number of the neutrino gas would acquire seizable direction-dependent variations, breaking the coherent behavior found in the spherically symmetric case. This finding would suggest that the previous results of the self-induced flavor evolution obtained in one-dimensional models should be critically re-examined.
We investigate the thermalization process of the Universe after inflation to determine the evolution of the effective temperature. The time scale of thermalization is found to be so long that it delays the evolution of the effective temperature, and the resulting maximal temperature of the Universe can be significantly lower than the one obtained in the literature. Our results clarify the finite density corrections to the effective potential of a scalar field and also processes of heavy particle production. In particular, we find that the maximum temperature of the Universe may be at most electroweak scale if the reheating temperature is as low as ${\cal O} (1)$ MeV, which implies that the electroweak symmetry may be marginally restored. In addition, it is noticeable that the dark matter may not be produced from thermal plasma in such a low reheating scenario, since the maximum temperature can be smaller than the conventional estimation by five orders of magnitude. We also give implications to the Peccei-Quinn mechanism and the Affleck-Dine baryogenesis.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)